Skip to main content

Is Hepatitis Delta infections important in Brazil?



The Hepatitis Delta Virus (HDV) can increase the incidence of fulminant hepatitis. For this infection occurs, the host must also be infected with Hepatitis B Virus. Previous studies demonstrated the endemicity and near exclusivity of this infection in the Amazon region, and as a consequence of the difficulty in accessing this area we used dried blood spots (DBS) in sample collection. The aims of this study were to investigate the presence of recombination, to analyze the epidemiology, ancestry and evolutionary pressures on HDV in Brazil.


Blood samples from 50 individuals were collected using dried-blood spots (DBS 903, Whatman), and sent via regular mail to Retrovirology Laboratory from Federal University of São Paulo, where the samples were processed. In the analysis the following software were used: PhyML, RDP, BEAST, jModelTest and CODEML.


Our results confirm the prevalence of HDV-3 in the Amazon region of Brazil, with the absence of inter-genotypic recombination. It was identified a positive selection in probable epitopes of HDV on B lymphocytes that might indicate that the virus is changing to escape the humoral response of the host. The analysis of the time of the most common ancestor demonstrated the exponential growth of this virus in late 1970s that lasted until 1995, after which it remained constant. It was also observed a probable founder effect in two cities, which demonstrate the need to focus on prevention methods against HBV/HDV infection.


We confirmed the prevalence of HDV-3 in the Amazon region of Brazil, without inter-genotypic recombination. The analysis of the time of the most common ancestor showed that this infection remain constant in the studied area. Taking into account the probable founder effect established in the cities of Rio Branco and Porto Velho, a focus on preventive methods is recommended against these infections.

Peer Review reports


Hepatitis delta is a disease that has been neglected by the health system [1], despite being responsible for the infection of 15 to 20 million people worldwide [2, 3]. In Brazil, 77 % of the infections with Hepatitis Delta Virus (HDV) occur in the North region.

The Northern region of Brazil has locations that are geographically isolated, which hinders the samples transportation until the processing laboratory. Thereat, it is used the dried blood spot (DBS) in sample collection at this region, for monitoring the HDV viral load (unpublished data). The use of DBS for sample collection and transportation is widely used in the HIV field, for monitoring the HIV viral load and the selection of drug resistance, on samples from remote sites or with limited financial resources. The results obtained using DBS sample collection are the same as using vacuum tubes [46].

HDV is the only subviral satellite agent known to infect humans. This virus requires the protection provided by the envelope of the Hepatitis B Virus (HBV) to enable its budding process from infected cells. Therefore, hepatitis delta occurs only in patients previously infected with HBV [7]. This co-infection increases the incidence of fulminant hepatitis is ten times [810], particularly when the characteristic Amazon genotypes (HDV-3 and HBV-F) are found combined [11, 12].

An important consideration regarding HDV is the availability of only one treatment option. This happens because this virus does not encode specific enzymes that could be used as target for the development of new antiviral drugs [3]. The efficacy of the treatment done with Interferon-Pegylated, is of only 43 %, and it produces very severe side effects [13].

Another situation that complicates the immune control of this infection is viral recombination. This event can allow the virus to escape the immune pressure imposed by the host's immune system, lead to an increase in viral fitness or even changes in the pathogenicity of the hybrid virus [14]. The recombination can occur when the individual is infected with a mixed viral population, as reported between genotypes 1 and 2 and between genotypes 1 and 4 in Taiwan [15, 16]. The probability of recombination occurring is greater when similar HDV genotypes cause the infection [1719].

Perhaps because HDV-3 is the most distantly related among all HDVs genotypes, recombination with HDV-3 has never been reported [11, 19]. Furthermore, most studies with HDV-3 analyze only the coding portion of the genome and perhaps analyzing the entire genome, that is a much more detailed approach, may lead to different results in the identification of inter-genotypic recombination.

The adaptation of the virus to the environment can be measured by non synonymous (dN) and synonymous substitutions (dS) present in their genomes. The ratio between dN/dS is used to measure the strength and mode of natural selection acting on protein coding genes [20]. Results > 1 indicate that a protein is under positive selection, whereas ratio  < 1 denotes that the protein is evolving slowly under negative or purifying selection, and dN/dS 1 implies that the protein is evolving neutrally [21].

Another form to better understand an infection is to analyze its diversification, ancestry and phylogeography. These informations can be used to shed light on the demographic, social and biological factors that were responsible for the emergence and spread of the disease present in the population [22, 23]. This approach may also help in the development of preventive strategies against that specific pathogen.

Therefore, according to the above information, this study aims to investigate the presence of recombination and to analyze the epidemiology, ancestry and evolutionary pressures on HDV in samples collected with DBS from the Northern region of Brazil.


Samples and processing

A total of 50 samples from patients in medical care in different regional clinics at the Amazon region of Brazil was collected between 2009 and 2010. The samples were originated from the states of Acre (Cruzeiro do Sul and Rio Branco Cities), Amazonas (Manaus City) and Rondônia (Porto Velho City), all located in the North region of Brazil. The samples were collected and transported in dried-blood spots (DBS 903- Whatman). The mean age of the patients was 32 years (maximum of 69 and minimum of 18), and 52 % were male.

The samples were processed only after the CEP/UNIFESP (# 72971) was approved, and analyzed anonymously to ensure the protection of the participants. Informed consent was not required because the samples had been previously collected and stored. Patients co-infected with HIV and/or HCV were excluded.

The sample RNA was purified by the method developed by Boom [24], and the cDNA were synthesized using SuperScript III reverse transcriptase (Invitrogen, Carlsbad, CA). HDV full genome was sequenced using the shotgun technique, which amplified five overlapping fragments named A - E. The PCR methodology used was based on the study performed by Nakano [25]. The primers used amplified all HDV genotypes, and the sensitivity of amplification was up to 10 3 copies/mL. For the D fragment a rescue strategy was necessary, where a nested PCR was implemented, combining the primers 1267_F with F2 [26] in the 1st round and A3 with 138_R in the 2nd round [26, 27].

The processing of all samples included a negative control to exclude the possibility of cross contamination and a positive control to ensure the reliability of the reaction. Moreover, both strands of the PCR were sequenced and thereafter a consensus sequence was created, which was used in the analysis to guarantee the best quality of the results obtained.

Sequencing was performed using the ABI Prism 3130 Genetic Analyzer (Applied Biosystems, USA). The sequences were edited using the Sequencher v.4.2 Program (Gene Code, Ann Arbor, USA) and analyzed using Blastn ), to verify their similarity with the database sequences.

Sequence alignment and phylogenetic analyses

The HDV references sequences were obtained from the GenBank database (, described in Table 1. All the sequences were aligned using ClustalX software version 2.1 [28]. After the first alignment, the sequences were manually edited using the SE Al program version 2.0 (

Table 1 Hepatitis Delta Virus sequences used as reference

Maximum likelihood trees and bootstrap values were obtained using PhyML software [29]. The HKY model and gamma distribution (Γ) were selected and, based on the likelihood ratio test (LRT), implemented in the jModeltest software [30, 31]. The resulting trees were created and edited using FigTree (

Recombination analyses

The RDP v.4 software was used to determine the presence of recombination in the sequences [32]. 3Seq performs an exact nonparametric test to detect recombination in triplets of sequences [33]. The maximum χ 2 (MaxChi) is a method implemented by Maynard Smith [34], and it uses variable/invariable sites to detect recombination in pairs of sequences. The maximum match χ 2 (Chimaera) is a modification of Smith’s method that uses variable sites to calculate the maximum χ 2 match statistics [35]. Geneconv detects recombination by evaluating conserved substitutions in fragments between two sequences [36]. Although evolutionary methods are not explicitly implemented in Geneconv, it is robust and has low levels of false positive [37].

Initially, default parameters were used; then, these were optimized to avoid detection of false positive recombination. In addition, window sizes of 20 to 50 as well as Bonferroni correction were utilized.

Phylogeographic analysis

A Bayesian Markov chain Monte Carlo (BMCMC) coalescent framework was used to estimate the ancestral genealogy, phylogeographic and time to the most common ancestor (tMRCA). After Markov chain runs enough so that the parameter space have been sufficiently visited, the mean and the 95 % HPD intervals of each parameter was calculated. The HKY model plus a gamma correction (Γ) was applied to all analyses; the evolutionary and demographic parameters were iteratively adjusted. We used the Bayesian discrete methods to estimate the node probability in regard to the geographical location of samples [38].

In order to explore the correlation between sample location and tree topology, the Bayesian discrete model [38] was used because it allows the reconstruction of ancestral location states while accounting for the uncertainties of maximum posteriori (MP) tree topology and mapping states at nodes.

The Bayesian skyline plot method (BSP) that provides an unbiased description of genetic diversity over time was used (expressed as the product of the effective sample size Ne and generation time τ). BSP contemplates a wide varied of demographic scenarios without formally assuming any determined model [39]. Since absolute time was used (years) to scale branch lengths and a specific generation time was not assumed, our estimates of Ne.τ will reflect only the relative genetic diversity of the infections over time.

The Markov Chain Monte Carlo (MCMC) processes were run for 3x108 generations with the initial 10 % of each run discarded as burn in. The convergence of chains was evaluated using the TRACER software, version 1.5 (available at:, runs were accepted when all parameters presented the effective sample size number (ESS) greater than 200. Two independent chains were run for each dataset and combined with LogCombiner software [39, 40].

The MP tree was inferred using the BSP method and the maximum clade credibility parameter were annotated using the TreeAnnotator software [39, 40]. All these analyses were performed with the BMCMC approaches implemented in the BEAST package version 1.8.0 [40].

Estimates of the nonsynonymous (dN) and synonymous (dS) substitution rates

The codon based maximum likelihood method was used to estimate the selective pressures on the PTV sequences. This approach estimates the likelihood of distinct models of codon evolution and computes the ratio (dN/dS = ω) of the number of non-synonymous (dN) and synonymous (dS) substitution rates between sites considering the phylogenetic relationships of the sequences [41, 42]. To estimate the overall mean ratio, the one ratio model (M0) was used. This model assumes a single ω for all sites in the alignment and is the simplest model. The hypothesis of positive evolution is based on the likelihood ratio test (LRT) comparing pairs of nested models, with the corresponding statistics approximated by a chi square (χ 2) distribution. To detect positive selection, we compared the null hypothesis of model M1, which allows for different proportions of conserved (ω = 0) and neutral codons (ω = 1), against the alternative hypothesis of model M2, which has an additional class of codons with ω ratio > 1, thus assuming positive selection. These models were implemented in the CODEML program of the PAML package, version 4.8, available at


Phylogenetic overview

Out of a total of 50 samples, 46 demonstrated the necessary quality for further analysis. The demographic information and the GenBank ID of each sample used in the analysis are detailed in Table 2 (present in the end of the manuscript). A maximum likelihood tree was initially constructed using the full length genomes of 46 samples and reference genomes of the main genotypes. Figure 1 shows that all sequences generated in the current study belong to genotype 3 and clustered within a single sequence with a bootstrap proportion of 97. Although the samples belong to genotype 3, they formed two groups: one composed of isolates from Yucpa Indians in Venezuela [43] and another composed of Brazilian isolates plus one isolate from the Yavari River in the Amazon basin close to the Peru Brazil border [19]. All sequences used to infer the tree mentioned above lacked a recombination signal.

Table 2 Demographic information and GenBank ID of the samples from our study
Fig. 1

Phylogenetic analysis of the sequenced HDV genomes. Each triangle includes the sequences of a same genotype. The Maximum likelihood (ML) tree was constructed using the HKY + G + I model as implemented in PhyML software; values at the nodes of the tree are bootstrap support values obtained using 500 replicates. The tree was constructed using the genomes of 42 isolates generated in this study and the following reference genomes: genotype 1) IR_1_AY633627_Iran, TW_1_TWD2577-66_A_Taiwan, SO_1c_U81988_Somalia, ET_1C_HDU81989_Ethiopia; genotype 2) FR_2_AX74129_France, JP_2_AB118846_M37_Japan, TW_2_X60193_Taiwan; genotype 3) VE_AB037947_VnzD8375_Venezuela, VE_AB037949_VnzD8624_Venezuela, VE_AB037948_VnzD8349_Venezuela, PE_3_L22063_Peru; genotype 4) JP_4_AB118842_Japan, TW_4_AF209859_Taiwan; genotype 5) TG_5_AM183326_Togo, SN_5_AM183328_Senegal, GW_5_AM18333_Guinea-Bissau; genotype 6) NG_6_AM183329_Nigeria, CM_6_AJ584847_Cameroon; genotype 7) CM_7_AM183333_Cameroon, CM_7_AJ584844_Cameroon; genotype 8) SN_8_AM183327_Senegal, CD_8_AJ584849_Congo

Phylogeographic analysis of HDV-3 in the Amazon region

42 samples were used for this analysis, four of the 46 original sequences were excluded because they had significant deletions along the genome. The phylogeography of HDV-3 isolates (Fig. 2), show that Brazilian samples clustered into four clades (named A, B, C, D), one sample from Peru isolated in 1986 is located at the base of these clades. The MP tree presents two large clusters, A and B, containing samples from different geographic locations (Rio Branco, Porto Velho, Cruzeiro do Sul and Manaus) and two small clusters, C and D, with sequences from Rio Branco and Porto Velho, respectively.

Fig. 2

Phylogeography analysis of the HDV-3. The branch colors indicate the location from which samples were collected. The numbers at nodes indicate the location probability. The estimates were obtained using the full-length genomes of HDV-3 samples isolated from distinct locations in the Brazilian Amazon basin (i.e., Manaus, Porto Velho, Cruzeiro do Sul, Rio Branco). HDV-3 references from Venezuela and Peru were also used. The scale under the tree indicates the number of years before present time

The topology also indicates two sets of nearly identical isolates: one from Rio Branco and another from Porto Velho (filled triangles in Fig. 2). The branching pattern of intermingled samples from Rio Branco and Porto Velho suggests a bidirectional flow of infections between these two locations.

The pattern of geographical dispersal of HDV-3 is depicted by the highest location probability on the nodes of the MP tree (Fig. 2; the location is indicated by colored branches and the location probability by numbers at the nodes). Although the multi-location clades A and B indicate that Rio Branco contributes significantly to the geographical spread of HDV-3, this may occur because samples from this location were overrepresented. When random sets of sequences containing equivalent numbers of samples from Rio Branco and Porto Velho were evaluated, the location probability to support Rio Branco at internal nodes decreased to values below 0.7 (data not shown).

The population dynamics of HDV-3

Figure 3 shows the reconstructed dynamics of the genetic diversity (Ne.τ, y axis) plotted against time in years (x axis). The plot shows the mean (thick solid line) and the 95 % highest density interval (thin solid lines). The shape of the plot demonstrates a period of exponential growth that began in the late 1970s and lasted until 1995. The estimated tMRCA values for the HDV-3 and for the Amazon clade were 1934 (1927 until 1940) and 1952 (1948 until 1955), respectively.

Fig. 3

Bayesian skyline plot (BSP) of the population dynamics of HDV-3. The plot represents the genetic diversity (y-axis) against chronological time expressed in years (x-axis). The thick line indicates the mean and the thin lines the 95 % confidence posterior interval of the Ne.τ. The estimates were obtained using the full-length genomes of 42 HDV-3 isolates

Selective regimen in the large HDV antigen (L-HDAg)

A systematic codon based analysis was performed on HDV isolates using the 215 codons from the L-HDAg region. In Table 3, the estimates of the codon based analysis are summarized. To detect positive selection models M7 and M8 were compared using likelihood ratio test assuming two degrees of freedom between models. Initially an intra-subtype analysis was performed using a dataset with 47 isolates of HDV-3, 46 from this study and one reference sequence from Peru (L22063), results of Model M8 indicated 24 codons in the L-HDAg under positive selection with mean ω values of 5.04.

Table 3 Codon-based analysis of dN/dS rates in the L-HDAg gene of Hepatitis Delta virus

When this same analysis was performed with other HDV genotypes (78 sequences) the Model M8 indicated only nine codons under positive selection, four of which were similar with the result demonstrated with the HDV-3 genotype (8 K, 9 K, 37 T and 120S).


Chronic HDV/HBV infection is perhaps the most intriguing and difficult to treat among the human viral hepatitis and is one of the most neglected diseases worldwide [44]. The study area included was based upon previous articles that demonstrated the endemicity and near exclusivity of HDV infection in the Amazon region of Brazil [12, 19, 45, 46], only two studies have identified HDV cases outside the locations presented in this study [47, 48].

In Brazil, 77 % of the HDV infections occur in the North region [48]. The isolation of this virus in this location can be a consequence that the majority of HDV cases are present in the indigenous population, that generally lives in isolated regions and interacts preferentially among individuals of their own tribes.

The DBS was used to collect the samples because: distant areas were included, this method is less invasive for the patient once the blood is collected via puncture of the finger using a lancet, it requires minimal collector training, allows the samples to be stored at room temperature with transportation in envelopes sent via regular mail (less expensive than using dried ice transportation) [49, 50].

Based upon the phylogenetic tree created from the sequences obtained in this study, it was possible to conclude that all samples belongs to genotype 3, being clustered (97 in the bootstrap) along with one HDV-3 sequence from Peru. These results confirm previous studies, that also demonstrated the high prevalence of this genotype in the studied area [12, 19, 45, 46].

It was not detected any inter-genotypic recombination. This may have happened because we analyzed areas of difficult access, which may have hindered the introduction of different HDV genotypes.

The BMCMC method was used in the phylogeographic analysis because it incorporates the uncertainty of the measurements by considering the intrinsic errors both in the tree reconstruction and in the coalescent method. These considerations guarantee more confiability in the results obtained. This analysis demonstrated the presence of four clades in HDV-3 samples from the Amazon region (A, B, C and D). Two of them (A and B) contained samples from all the locations studied (Rio Branco, Porto Velho, Cruzeiro do Sul and Manaus), whereas clusters C and D were significantly smaller and represented the cities of Rio Branco and Porto Velho, respectively. This analysis also indicates two clusters with nearly identical sequences, one from Rio Branco (composed of four samples) and another from Porto Velho (composed of five samples), which suggests a bidirectional flow of infection in these two cities.

Rio Branco is the capital city of the Acre State, populated in the beginning of the 20th century by natives who already lived in the area, along with the migration of people from Northeast of Brazil and immigrants from Turkey, Portugal and Lebanon attracted by the enrichment of the rubber era [51]. The cluster identified at this location indicates a probable founder effect, that may suggest that the dissemination of this infection was introduced by an indian (HDV-3 is characteristic of the indigenous population), and not brought with the city foundation.

Porto Velho, the capital city of Rondônia State, was founded as a result of the railroad “Madeira Mamoré”, widely used in the rubber and gold eras. Nowadays several indigenous tribes still live in this region [51]. Five samples from this location were probably indigenous (surnames characteristic from tribes living in this area), and they grouped together into a single cluster. This situation demonstrates that probably a single virus caused the spread of the HDV infection in this region, as well.

These probable founder effects demonstrate the necessity of focusing on prevention strategies against HBV and HDV infection, mainly in the North of Brazil. This region represents respectively 14 and 77 % of the total cases of these two viral infections. The most effective prevention against HDV is the HBV vaccination, since the HDV depends of the HBV for its budding from the infected cell, and if HDV is not contained it can be disseminated to other regions of Brazil, causing a public health problem.

The analysis of the most common ancestor time (tMRCA) showed that the infection of HDV experienced an exponential growth period in the late 1970s that lasted until 1995 and has been maintained since then. In 1989, the HBV vaccination was implemented in the Amazon region and in 1993, the minimal age for vaccination was decreased, these events could explains the plateau in the exponential growth of HDV after 1995; but this same plateau does not exclude the need to focus on preventive campaigns against this infection, since its treatment is not effective, no remission has been observed and the vaccination against HBV is not effective in all vaccinated individuals.

It was also observed that the HDV-3 has more codons under positive selection in the L-HDAg than the other HDV genotypes, but four codons (8 K, 9 K, 37 T and 120S) were present in both analysis indicating that these modification are not genotype specific. Some modifications were located in probable epitopes of HDV for B lymphocytes (amino acids 3–18; 107–122 and 170–196) [52, 53]. This finding may indicate that the virus is evolving to escape the humoral response of the host.

The variation of epitopes at amino acids 174 to 195, often emerges after severe hepatitis attacks during the chronic phase of infections [15]. Changes in epitopes are common during the chronic phase of infections because of the intense selection forces by the immune defense against the virus [54]. Unfortunately, we did not have information to confirm that these individuals were in this phase of infection.


Our results confirm the prevalence of HDV-3 in the Amazon region of Brazil. This geographically isolated region may have hindered the introduction of different HDV genotypes, which may account for the absence of inter-genotypic recombination. The analysis of the most common ancestor time showed that the infection with HDV experienced an exponential growth period in the late 1970s that lasted until 1995, after which it has remained constant. Taking into account the probable founder effect established in the cities of Rio Branco and Porto Velho, a focus on preventive methods is recommended against this infection, particularly in the North of Brazil. A positive selection was also found in probable epitopes of HDV for B lymphocytes, which may indicate that the virus is evolving to escape the humoral response of the host.



Bayesian Markov chain Monte Carlo


Bayesian skyline plot method


Complementary DNA


Dried blood spots


Non synonymous substitutions


Synonymous substitutions


Effective sample size number


Hepatitis B Virus


Hepatitis Delta Virus


Highest posterior density


Large HDV antigen region protein


Likelihood ratio test


Maximum posteriori


Polymerase chain reaction


Time to the most common ancestor


  1. 1.

    Genne D, Rossi I. Hepatitis delta in Switzerland: a silent epidemic. Swiss Med Wkly. 2011;141:w13176.

    PubMed  Google Scholar 

  2. 2.

    Cross TJ, Rizzi P, Horner M, Jolly A, Hussain MJ, Smith HM, Vergani D, Harrison PM. The increasing prevalence of hepatitis delta virus (HDV) infection in South London. J Med Virol. 2008;80:277–82.

    Article  PubMed  Google Scholar 

  3. 3.

    Wedemeyer H, Manns MP. Epidemiology, pathogenesis and management of hepatitis D: update and challenges ahead. Nat Rev Gastroenterol Hepatol. 2010;7:31–40.

    Article  PubMed  Google Scholar 

  4. 4.

    de Truchis P, Le MP, Daou M, Madougou B, Nouhou Y, Moussa Saley S, Sani A, Adehossi E, Rouveix E, Saidou M, et al. High efficacy of first-line ART in a West African cohort, assessed by dried blood spot virological and pharmacological measurements. J Antimicrob Chemother. 2016.

  5. 5.

    Gueye SB, Diop-Ndiaye H, Diallo MM, Ly O, Sow-Ndoye A, Diagne-Gueye ND, Kebe-Fall K, Diop F, Gaye-Diallo A, Belec L, et al. Performance of Roche CAP/CTM HIV-1 qualitative test version 2.0 using dried blood spots for early infant diagnosis. J Virol Methods. 2016;229:12–5.

    CAS  Article  PubMed  Google Scholar 

  6. 6.

    Wei X, Smith AJ, Forrest DW, Cardenas GA, Beck DW, LaLota M, Metsch LR, Sionean C, Owen SM, Johnson JA. Incident infection and resistance mutation analysis of dried blood spots collected in a field study of HIV risk groups, 2007–2010. PLoS One. 2016;11:e0159266.

    Article  PubMed  PubMed Central  Google Scholar 

  7. 7.

    Rizzetto M, Canese MG, Arico S, Crivelli O, Trepo C, Bonino F, Verme G. Immunofluorescence detection of new antigen-antibody system (delta/anti-delta) associated to hepatitis B virus in liver and in serum of HBsAg carriers. Gut. 1977;18:997–1003.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  8. 8.

    Smedile A, Farci P, Verme G, Caredda F, Cargnel A, Caporaso N, Dentico P, Trepo C, Opolon P, Gimson A, et al. Influence of delta infection on severity of hepatitis B. Lancet. 1982;2:945–7.

    CAS  Article  PubMed  Google Scholar 

  9. 9.

    Lai MM. The molecular biology of hepatitis delta virus. Annu Rev Biochem. 1995;64:259–86.

    CAS  Article  PubMed  Google Scholar 

  10. 10.

    Purcell RH, Gerin JL. Hepatitis Delta virus. In: Fields BN, Knipe DM, Howley PM, editors. Fields virology. 3rd ed. Philadelphia: Lippincott-Raven; 1996. p. 2819–29.

    Google Scholar 

  11. 11.

    Casey JL, Niro GA, Engle RE, Vega A, Gomez H, McCarthy M, Watts DM, Hyams KC, Gerin JL. Hepatitis B virus (HBV)/hepatitis D virus (HDV) coinfection in outbreaks of acute hepatitis in the Peruvian Amazon basin: the roles of HDV genotype III and HBV genotype F. The Journal of infectious diseases. 1996;174:920–6.

    CAS  Article  PubMed  Google Scholar 

  12. 12.

    Quintero A, Uzcategui N, Loureiro CL, Villegas L, Illarramendi X, Guevara ME, Ludert JE, Blitz L, Liprandi F, Pujol FH. Hepatitis delta virus genotypes I and III circulate associated with hepatitis B virus genotype F in Venezuela. J Med Virol. 2001;64:356–9.

    CAS  Article  PubMed  Google Scholar 

  13. 13.

    Han Z, Nogusa S, Nicolas E, Balachandran S, Taylor J. Interferon impedes an early step of hepatitis delta virus infection. PLoS One. 2011;6:e22415.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  14. 14.

    Chao M. RNA recombination in hepatitis delta virus: implications regarding the abilities of mammalian RNA polymerases. Virus Res. 2007;127:208–15.

    CAS  Article  PubMed  Google Scholar 

  15. 15.

    Wu JC, Chiang TY, Shiue WK, Wang SY, Sheen IJ, Huang YH, Syu WJ. Recombination of hepatitis D virus RNA sequences and its implications. Mol Biol Evol. 1999;16:1622–32.

    CAS  Article  PubMed  Google Scholar 

  16. 16.

    Wang TC, Chao M. RNA recombination of hepatitis delta virus in natural mixed-genotype infection and transfected cultured cells. J Virol. 2005;79:2221–9.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  17. 17.

    McCahon D, Slade WR, Priston RA, Lake JR. An extended genetic recombination map for foot-and-mouth diseases virus. J Gen Virol. 1977;35:555–65.

    CAS  Article  PubMed  Google Scholar 

  18. 18.

    McCahon D, King AM, Roe DS, Slade WR, Newman JW, Cleary AM. Isolation and biochemical characterization of intertypic recombinants of foot-and-mouth disease virus. Virus Res. 1985;3:87–100.

    CAS  Article  PubMed  Google Scholar 

  19. 19.

    Casey JL, Brown TL, Colan EJ, Wignall FS, Gerin JL. A genotype of hepatitis D virus that occurs in northern South America. Proc Natl Acad Sci U S A. 1993;90:9016–20.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  20. 20.

    Jeffares DC, Tomiczek B, Sojo V, Dos Reis M. A beginners guide to estimating the non-synonymous to synonymous rate ratio of all protein-coding genes in a genome. Methods Mol Biol. 2015;1201:65–90.

    CAS  Article  PubMed  Google Scholar 

  21. 21.

    Kimura M. Preponderance of synonymous changes as evidence for the neutral theory of molecular evolution. Nature. 1977;267:275–6.

    CAS  Article  PubMed  Google Scholar 

  22. 22.

    Markov PV, Pepin J, Frost E, Deslandes S, Labbe AC, Pybus OG. Phylogeography and molecular epidemiology of hepatitis C virus genotype 2 in Africa. The Journal of general virology. 2009;90:2086–96.

    CAS  Article  PubMed  Google Scholar 

  23. 23.

    Lago BV, Mello FC, Kramvis A, Niel C, Gomes SA. Hepatitis B virus subgenotype A1: evolutionary relationships between Brazilian, African and Asian isolates. PLoS One. 2014;9:e105317.

    Article  PubMed  PubMed Central  Google Scholar 

  24. 24.

    Boom R, Sol CJ, Salimans MM, Jansen CL, Wertheim-van Dillen PM, van der Noordaa J. Rapid and simple method for purification of nucleic acids. J Clin Microbiol. 1990;28:495–503.

    CAS  PubMed  PubMed Central  Google Scholar 

  25. 25.

    Nakano T, Lu L, Hu X, Mizokami M, Orito E, Shapiro C, Hadler S, Robertson B. Characterization of hepatitis B virus genotypes among Yucpa Indians in Venezuela. The Journal of general virology. 2001;82:359–65.

    CAS  Article  PubMed  Google Scholar 

  26. 26.

    Behzadian F, Sabahi F, Karimi M, Sadeghizadeh M, Maghsoudi N, Forooshani RS, Shahinsaz L. Molecular phylogenetic analysis of Iranian HDV complete genome. Virus Genes. 2005;30:383–93.

    CAS  Article  PubMed  Google Scholar 

  27. 27.

    Ma SP, Sakugawa H, Makino Y, Tadano M, Kinjo F, Saito A. The complete genomic sequence of hepatitis delta virus genotype IIb prevalent in Okinawa, Japan. The Journal of general virology. 2003;84:461–4.

    CAS  Article  PubMed  Google Scholar 

  28. 28.

    Larkin MA, Blackshields G, Brown NP, Chenna R, McGettigan PA, McWilliam H, Valentin F, Wallace IM, Wilm A, Lopez R, et al. Clustal W and Clustal X version 2.0. Bioinformatics. 2007;23:2947–8.

    CAS  Article  PubMed  Google Scholar 

  29. 29.

    Guindon S, Gascuel O. A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol. 2003;52:696–704.

    Article  PubMed  Google Scholar 

  30. 30.

    Posada D, Crandall KA. MODELTEST: testing the model of DNA substitution. Bioinformatics. 1998;14:817–8.

    CAS  Article  PubMed  Google Scholar 

  31. 31.

    Posada D. jModelTest: phylogenetic model averaging. Mol Biol Evol. 2008;25:1253–6.

    CAS  Article  PubMed  Google Scholar 

  32. 32.

    Martin DP, Williamson C, Posada D. RDP2: recombination detection and analysis from sequence alignments. Bioinformatics. 2005;21:260–2.

    CAS  Article  PubMed  Google Scholar 

  33. 33.

    Boni MF, Posada D, Feldman MW. An exact nonparametric method for inferring mosaic structure in sequence triplets. Genetics. 2007;176:1035–47.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  34. 34.

    Smith JM. Analyzing the mosaic structure of genes. J Mol Evol. 1992;34:126–9.

    CAS  PubMed  Google Scholar 

  35. 35.

    Posada D, Crandall KA. Evaluation of methods for detecting recombination from DNA sequences: computer simulations. Proc Natl Acad Sci U S A. 2001;98:13757–62.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  36. 36.

    Sawyer S. Statistical tests for detecting gene conversion. Mol Biol Evol. 1989;6:526–38.

    CAS  PubMed  Google Scholar 

  37. 37.

    Bay RA, Bielawski JP: Recombination Detection Under Evolutionary Scenarios Relevant to Functional Divergence. J Mol Evol. 2012

  38. 38.

    Lemey P, Rambaut A, Drummond AJ, Suchard MA. Bayesian phylogeography finds its roots. PLoS Comput Biol. 2009;5:e1000520.

    Article  PubMed  PubMed Central  Google Scholar 

  39. 39.

    Drummond AJ, Rambaut A, Shapiro B, Pybus OG. Bayesian coalescent inference of past population dynamics from molecular sequences. Mol Biol Evol. 2005;22:1185–92.

    CAS  Article  PubMed  Google Scholar 

  40. 40.

    Drummond AJ, Suchard MA, Xie D, Rambaut A. Bayesian phylogenetics with BEAUti and the BEAST 1.7. Mol Biol Evol. 2012;29:1969–73.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  41. 41.

    Yang Z, Nielsen R. Codon-substitution models for detecting molecular adaptation at individual sites along specific lineages. Mol Biol Evol. 2002;19:908–17.

    CAS  Article  PubMed  Google Scholar 

  42. 42.

    Yang Z, Nielsen R, Goldman N, Pedersen AM. Codon-substitution models for heterogeneous selection pressure at amino acid sites. Genetics. 2000;155:431–49.

    CAS  PubMed  PubMed Central  Google Scholar 

  43. 43.

    Nakano T, Shapiro CN, Hadler SC, Casey JL, Mizokami M, Orito E, Robertson BH. Characterization of hepatitis D virus genotype III among Yucpa Indians in Venezuela. The Journal of general virology. 2001;82:2183–9.

    CAS  Article  PubMed  Google Scholar 

  44. 44.

    Braga WS, de Oliveira CM, de Araujo JR, Castilho Mda C, Rocha JM, Gimaque JB, Silva ML, Vasconcelos HL, Ramasawmy R, Parana R. Chronic HDV/HBV co-infection: predictors of disease stage---a case series of HDV-3 patients. J Hepatol. 2014;61:1205–11.

    Article  PubMed  Google Scholar 

  45. 45.

    Gomes-Gouvea MS, Soares MC, Bensabath G, de Carvalho-Mello IM, Brito EM, Souza OS, Queiroz AT, Carrilho FJ, Pinho JR. Hepatitis B virus and hepatitis delta virus genotypes in outbreaks of fulminant hepatitis (Labrea black fever) in the western Brazilian Amazon region. The Journal of general virology. 2009;90:2638–43.

    CAS  Article  PubMed  Google Scholar 

  46. 46.

    Gomes-Gouvea MS, Pereira Soares Mdo C, de Carvalho Mello IM G, Brito EM, Pereira Moia Lde J, Bensabath G, Nunes HM, Carrilho FJ, Pinho JR. Hepatitis D and B virus genotypes in chronically infected patients from the Eastern Amazon Basin. Acta Trop. 2008;106:149–55.

    CAS  Article  PubMed  Google Scholar 

  47. 47.

    Barros LM, Gomes-Gouvea MS, Pinho JR, Alvarado-Mora MV, Dos Santos A, Mendes-Correa MC, Caldas AJ, Sousa MT, Santos MD, Ferreira AS. Hepatitis Delta virus genotype 8 infection in Northeast Brazil: inheritance from African slaves? Virus Res. 2011;160:333–9.

    CAS  Article  PubMed  Google Scholar 

  48. 48.

    Mendes-Correa MC, Gomes-Gouvea MS, Alvarado-Mora MV, Da Silva MH, Lazari C, Cavalcanti NC, Alonso FK, Carpinelli CC, Uip DE, Pinho JR. Hepatitis delta in HIV/HBV co-infected patients in Brazil: is it important? Int J Infect Dis. 2011;15:e828–832.

    Article  PubMed  Google Scholar 

  49. 49.

    Garcia Boy R, Henseler J, Mattern R, Skopp G. Determination of morphine and 6-acetylmorphine in blood with use of dried blood spots. Ther Drug Monit. 2008;30:733–9.

    CAS  Article  PubMed  Google Scholar 

  50. 50.

    la Marca G, Malvagia S, Filippi L, Luceri F, Moneti G, Guerrini R. A new rapid micromethod for the assay of phenobarbital from dried blood spots by LC-tandem mass spectrometry. Epilepsia. 2009;50:2658–62.

    Article  PubMed  Google Scholar 

  51. 51.

    Instituto Brasileiro de Geografia e Estatística (IBGE): Estimativa populacional 2011. Prefeitura Municipal de Porto Velho: Instituto Brasileiro de Geografia e Estatística (IBGE); 2011.

  52. 52.

    Nisini R, Paroli M, Accapezzato D, Bonino F, Rosina F, Santantonio T, Sallusto F, Amoroso A, Houghton M, Barnaba V. Human CD4+ T-cell response to hepatitis delta virus: identification of multiple epitopes and characterization of T-helper cytokine profiles. J Virol. 1997;71:2241–51.

    CAS  PubMed  PubMed Central  Google Scholar 

  53. 53.

    Wang JG, Jansen RW, Brown EA, Lemon SM. Immunogenic domains of hepatitis delta virus antigen: peptide mapping of epitopes recognized by human and woodchuck antibodies. J Virol. 1990;64:1108–16.

    CAS  PubMed  PubMed Central  Google Scholar 

  54. 54.

    Wang SY, Wu JC, Chiang TY, Huang YH, Su CW, Sheen IJ. Positive selection of hepatitis delta antigen in chronic hepatitis D patients. J Virol. 2007;81:4438–44.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

Download references


We would like to thank Dr. Alexanda Reis and Centro de Genomas for their methodological assistance and all the patients who enabled this study.


We would like to thank the financial support obtained from São Paulo Research Foundation (FAPESP, grant#2012/13285-0), National Council for Scientific and Technological Development (CNPQ, grant#135382/2011-4).

Availability of data and material

The datasets generated during and/or analysed during the current study are available in the GenBank repository ( and the GenBank ID information is present in Table 2.

Authors’ contributions

RSD and SVK were responsible for the conception and designed of the study; MCF and SVK wrote the manuscript; MCF, NMP, LCS, RA, RGA carried out the molecular studies and participated in the sequence alignment; MCF, LCS, ESL participated in the sequence alignment. All authors read and approved the final manuscript.

Competing interests

The authors declare that they have no competing interests.

Consent for publication

Not applicable.

Ethics approval and consent to participate

The Ethics approval for this study was number # 72971 (CEP/UNIFESP), and the informed consent was not required because the samples had been previously collected and stored in the laboratory.

Author information



Corresponding author

Correspondence to Shirley Vasconcelos Komninakis.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Cicero, M.F., Pena, N.M., Santana, L.C. et al. Is Hepatitis Delta infections important in Brazil?. BMC Infect Dis 16, 525 (2016).

Download citation


  • Hepatitis Delta virus
  • Dried-blood spot
  • Founder effect
  • Phylogeography
  • Phylogenetic analysis
  • tMRCA
  • Brazil