Compilation of human mtDNA control region sequences.

This paper describes the organisation of a database for human mitochondrial control-region sequences. The data are divided into three ASCII files that contain aligned sequences from the hypervariable region I (HVRI), from the hypervariable region II (HVRII), and the available information about the individuals, from whom the sequences stem. The current collection comprises 4079 HVRI and 969 HVRII sequences. From 728 individuals sequences of both HVRI and HVRII are available. For easy access, the collection is made available to the scientific community via World Wide Web at URL http://www.zi.biologie.uni-muenchen.de/[symbol: see text]meyers/mtdna.html     

HvrBase: compilation of mtDNA control region sequences from primates.

HvrBase is a compilation of human and ape mtDNA control region sequences. Sequences and related information on individuals, such as from where the sequences were obtained, is stored in three ASCII files as described previously. Moreover, the collection is also available as Mac/PC database application with a graphical user interface. It can be accessed through the WWW at URL http://www.eva.mpg.de/hvrbase. The current collection comprises 5846 human sequences from hypervariable region I (HVRI) and 2302 human sequences from hypervariable region II (HVRII). From apes, 295 HVRI sequences and 13 HVRII sequences are available.     

Frequency and Pattern of Heteroplasmy in the Control Region of Human Mitochondrial DNA

In this work, we present the results of the screening of human mitochondrial DNA (mtDNA) heteroplasmy in the control region of mtDNA from 210 unrelated Spanish individuals. Both hypervariable regions of mtDNA were amplified and sequenced in order to identify and quantify point and length heteroplasmy. Of the 210 individuals analyzed, 30% were fully homoplasmic and the remaining presented point and/or length heteroplasmy. The prevalent form of heteroplasmy was length heteroplasmy in the poly(C) tract of the hypervariable region II (HVRII), followed by length heteroplasmy in the poly(C) tract of hypervariable region I (HVRI) and, finally, point heteroplasmy, which was found in 3.81% of the individuals analyzed. Moreover, no significant differences were found in the proportions of the different kinds of heteroplasmy in the population when blood and buccal cell samples were compared. The pattern of heteroplasmy in HVRI and HVRII presents important differences. Moreover, the mutational profile in heteroplasmy seems to be different from the mutational pattern detected in population. The results suggest that a considerable number of mutations and, particularly, transitions that appear in heteroplasmy are probably eliminated by drift and/or by selection acting at different mtDNA levels of organization. Taking as a whole the results reported in this work, it is mandatory to perform a broad-scale screening of heteroplasmy to better establish the heteroplasmy profile which would be important for medical, evolutionary, and forensic proposes.     

Intraspecific DNA sequence variation of the mitochondrial control region of white sturgeon (Acipenser transmontanus)

Intraspecific sequence variation in the D-loop region of mtDNA in white sturgeon (Acipenser transmontanus), a relict North American fish species, was examined in 27 individuals from populations of the Columbia and Fraser rivers. Thirty-three varied nucleotide positions were present in a 462-nucleotide D-loop sequence, amplified using the polymerase chain reaction. Bootstrapped neighbor-joining and maximum- parsimony trees of sequences from 19 haplotypes suggest that the two populations have recently diverged. This is consistent with the hypothesis that the Columbia River, a Pleistocene refugium habitat, was the source of founders for the Fraser River after the last glacial recession. On the basis of a divergence time of 10-12 thousand years ago, the estimated substitution rate of the white sturgeon D-loop region is 1.1-1.3 x 10(-7) nucleotides/site/year, which is comparable to rates for hypervariable sequences in the human D-loop region. Furthermore, the ratio of mean percent nucleotide differences in the D- loop (2.27%) to that in whole mtDNA (0.54%, as estimated from restriction-enzyme data) is 4.3, which is similar to the fourfold-to- fivefold-higher substitution rate estimated for the human D-loop. The high nucleotide substitution rate of the hypervariable region indicates that the vertebrate D-loop has potential as a genetic marker in molecular population studies.      

Hepatitis C Virus Hypervariable Region 1 Variants Presented on Hepatitis B Virus Capsid-Like Particles Induce Cross-Neutralizing Antibodies

Hepatitis C virus (HCV) infection is still a serious global health burden. Despite improved therapeutic options, a preventative vaccine would be desirable especially in undeveloped countries. Traditionally, highly conserved epitopes are targets for antibody-based prophylactic vaccines. In HCV-infected patients, however, neutralizing antibodies are primarily directed against hypervariable region I (HVRI) in the envelope protein E2. HVRI is the most variable region of HCV, and this heterogeneity contributes to viral persistence and has thus far prevented the development of an effective HVRI-based vaccine. The primary goal of an antibody-based HCV vaccine should therefore be the induction of cross-reactive HVRI antibodies. In this study we approached this problem by presenting selected cross-reactive HVRI variants in a highly symmetric repeated array on capsid-like particles (CLPs). SplitCore CLPs, a novel particulate antigen presentation system derived from the HBV core protein, were used to deliberately manipulate the orientation of HVRI and therefore enable the presentation of conserved parts of HVRI. These HVRI-CLPs induced high titers of cross-reactive antibodies, including neutralizing antibodies. The combination of only four HVRI CLPs was sufficient to induce antibodies cross-reactive with 81 of 326 (24.8%) naturally occurring HVRI peptides. Most importantly, HVRI CLPs with AS03 as an adjuvant induced antibodies with a 10-fold increase in neutralizing capability. These antibodies were able to neutralize infectious HCVcc isolates and 4 of 19 (21%) patient-derived HCVpp isolates. Taken together, these results demonstrate that the induction of at least partially cross-neutralizing antibodies is possible. This approach might be useful for the development of a prophylactic HCV vaccine and should also be adaptable to other highly variable viruses.     

New Genetic Evidence on the Evolution of Chimpanzee Populations and Implications for Taxonomy

Primatologists widely recognize chimpanzees as belonging to a single species, Pan troglodytes, which they traditionally have further divided into 3 subspecies: west African P. t. verus, central African P. t. troglodytes, and east African P. t. schweinfurthii. Previously, we suggested that the phylogeographic history of chimpanzees may be different from that implied by the widely used taxonomy of the species. We based the suggestion on only a limited sample of haplotypes from the first hypervariable region (HVRI) of mitochondrial (mt)DNA from chimpanzees in Nigeria. We have now compiled a more geographically comprehensive genetic database for chimpanzees, including samples obtained near the Niger and Sanaga Rivers. Our database is composed of 254 HVRI haplotypes from chimpanzees of known geographic origin, including 79 unique HVRI haplotypes from chimpanzees living in Nigeria and Cameroon. The genetic data provide clear evidence that a major phylogeographic break between chimpanzee lineages occurs near the Sanaga River in central Cameroon and suggest the need for a reclassification of chimpanzees.     

Substitution rate variation among sites in mitochondrial hypervariable region I of humans and chimpanzees

Mitochondrial D-loop hypervariable region I (HVI) sequences are widely used in human molecular evolutionary studies, and therefore accurate assessment of rate heterogeneity among sites is essential. We used the maximum-likelihood method to estimate the gamma shape parameter alpha for variable substitution rates among sites for HVI from humans and chimpanzees to provide estimates for future studies. The complete data of 839 humans and 224 chimpanzees, as well as many subsets of these data, were analyzed to examine the effect of sequence sampling. The effects of the genealogical tree and the nucleotide substitution model were also examined. The transition/transversion rate ratio (kappa) is estimated to be about 25, although much larger and biased estimates were also obtained from small data sets at low divergences. Estimates of alpha were 0.28-0.39 for human data sets of different sizes and 0.20-0.39 for data sets including different chimpanzee subspecies. The combined data set of both species gave estimates of 0.42-0.45. While all those estimates suggest highly variable substitution rates among sites, smaller samples tend to give smaller estimates of alpha. Possible causes for this pattern were examined, such as biases in the estimation procedure and shifts in the rate distribution along certain lineages. Computer simulations suggest that the estimation procedure is quite reliable for large trees but can be biased for small samples at low divergences. Thus, an alpha of 0.4 appears suitable for both humans and chimpanzees. Estimates of alpha can be affected by the nucleotide sites included in the data, the overall tree length (the amount of sequence divergence), the number of rate classes used for the estimation, and to a lesser extent, the included sequences. The genealogical tree, the substitution model, and demographic processes such as population expansion do not have much effect.     

Analysis of sequence diversity in hypervariable regions of the external glycoprotein of human immunodeficiency virus type 1.

Nucleotide sequences in three hypervariable regions of the human immunodeficiency virus type 1 (HIV-1) env gene were obtained by sequencing provirus present in peripheral blood mononuclear cells of HIV-infected individuals. Single molecules of target sequences were isolated by limiting dilution and amplified in two stages by the polymerase chain reaction, using nested primers. The product was directly sequenced to avoid errors introduced by Taq polymerase during the amplification process. There was extensive variation between sequences from the same individual as well as between sequences from different individuals. Interpatient variability was markedly less in individuals infected from a common source. A high proportion of amino acid substitutions in the hypervariable regions altered the number and positions of potential N-linked glycosylation sites. Sequences in two hypervariable regions frequently contained short (3- to 15-bp) duplications or deletions, and by amplifying peripheral blood mononuclear cell DNA containing 10(2) or 10(3) proviral molecules and analyzing the product by high-resolution electrophoresis, the total number and abundance of distinct length variants within an individual could be estimated, providing a more comprehensive analysis of the variants present than would be obtained by sequencing alone. Sequences from many individuals showed frequent amino acid substitutions at certain key positions for neutralizing-antibody and cytotoxic T-cell recognition in the immunodominant loop. The rates of synonymous and nonsynonymous nucleotide substitution in the region of this and flanking regions indicate that strong positive selection for amino acid change is operating in the generation of antigenic diversity.     

Substitution rate variation among sites in hypervariable region 1 of human mitochondrial DNA

More than an order of magnitude difference in substitution rate exists among sites within hypervariable region 1 of the control region of human mitochondrial DNA. A two-rate Poisson mixture and a negative binomial distribution are used to describe the distribution of the inferred number of changes per nucleotide site in this region. When three data sets are pooled, however, the two-rate model cannot explain the data. The negative binomial distribution always fits, suggesting that substitution rates are approximately gamma distributed among sites. Simulations presented here provide support for the use of a biased, yet commonly employed, method of examining rate variation. The use of parsimony in the method to infer the number of changes at each site introduces systematic errors into the analysis. These errors preclude an unbiased quantification of variation in substitution rate but make the method conservative overall. The method can be used to distinguish sites with highly elevated rates, and 29 such sites are identified in hypervariable region 1. Variation does not appear to be clustered within this region. Simulations show that biases in rates of substitution among nucleotides and non-uniform base composition can mimic the effects of variation in rate among sites. However, these factors contribute little to the levels of rate variation observed in hypervariable region 1.     

Understanding differences between phylogenetic and pedigree-derived mtDNA mutation rate: a model using families from the Azores Islands (Portugal)

We analyzed the control region of the mitochondrial DNA (mtDNA)from maternally related individuals originating from the AzoresIslands (Portugal) in order to estimate the mutation rate ofmtDNA and to gain insights into the process by which a new mutationarises and segregates into heteroplasmy. Length and/or pointheteroplasmies were found at least in one individual of 72%of the studied families. Eleven new point substitutions werefound, all of them in heteroplasmy, from which five appear tobe somatic mutations and six can be considered germinal, evidencingthe high frequency of somatic mutations in mtDNA in healthyyoung individuals. Different values of the mutation rate accordingto different assumptions were estimated. When considering allthe germinal mutations, the value of the mutation rate obtainedis one of the highest reported so far in family studies. However,when corrected for gender (assuming that the mutations presentin men have the same evolutionary weight of somatic mutationsbecause they will inevitably be lost) and for the probabilityof intraindividual fixation, the value for the mutation rateobtained for HVRI and HVRII (0.2415 mutations/site/Myr) wasin the upper end of the values provided by phylogenetic estimations.These results indicate that the discrepancy, that has been reportedpreviously, between the human mtDNA mutation rates observedalong evolutionary timescales and the estimations obtained usingfamily pedigrees can be minimized when corrections for genderproportions in newborn individuals and for the probability ofintraindividual fixation are introduced. The analyses performedsupport the hypothesis that (1) in a constant, tight bottleneckgenetic drift alone can explain different patterns of heteroplasmysegregation and (2) in neutral conditions, the destiny of anew mutation is strictly related to the initial proportion ofthe new variant. Another important point arising from the dataobtained is that, even in the absence of a paternal contributionof mtDNA, recombination may occur between mtDNA molecules presentin an individual, which is only observable if it occurs betweenmtDNA types that differ at two or more positions.     

The impact of population expansion and mutation rate heterogeneity on DNA sequence polymorphism

In order to study the effect of mutation rate heterogeneity on patterns of DNA polymorphism, we simulated samples of DNA sequences with gamma- distributed nucleotide substitution rates in stationary and expanding populations. We find that recent population expansions and mutation rate heterogeneity have similar effects on several polymorphism indicators, like the shape and the mean of the observed pairwise difference distribution, or the number of segregating sites. The inferred size of population expansion thus appears overestimated if nucleotides have dissimilar substitution rates. Interestingly, population expansion and uneven mutation rates have contrasting effects on Tajima's D statistic when acting separately, and the consequence on the associated test of selective neutrality is investigated. The patterns of polymorphism of several human populations analyzed for the mitochondrial control region are examined, mainly showing the difficulty in quantifying the respective contribution of past demographic history and uneven mutation rates from a single sampled evolutionary process. However, substitution rates appear more heterogeneous in the second hypervariable segment of the control region than in the first segment.      

On the estimation of the rate of nucleotide substitution for the control region of human mitochondrial DNA

To apply molecular clock for studying human evolution, the pattern of nucleotide substitution for the control region of human mtDNA was analyzed in detail. It is well known that the rate of nucleotide substitution for the control region is much higher than that for any other part of mtDNA. In this study, the higher substitution rate was attributed to the higher rate of transition-type substitution between pyrimidines within the D-loop part, whereas the rates of other types of substitution were essentially the same over the entire mtDNA molecule. Even within the control region, the rate and pattern of nucleotide substitution were different between the D-loop part and the rest. The rate and pattern for the non-D-loop part were very similar to those for fourfold-degenerate sites in the protein-coding region. In contrast, the D-loop and non-D-loop parts showed similarities in the base composition, whereas the base composition of fourfold-degenerate sites slightly different from that of the both parts of the control region. It is concluded, therefore, that the nucleotide frequencies of the control region should be used to estimate the number of substitutions (d) between the control region sequences. However, a method to verify the accuracy of the estimation of d by means of the transition/transversion (s/v) ratio was theoretically studied. It was suggested that the s/v ratio becomes constant over a wide range of d values only when the estimation of d is unbiased. On the basis of this result, the estimates of d previously obtained between human sequences were evaluated.     

Application of the hypervariable region of the 16S rDNA sequence as an index for the rapid identification of species in the genus Alicyclobacillus

A comparison of the 16S rRNA gene (rDNA) sequences of seven type strains belonging to different Alicyclobacillus species (i.e., five validated species, one proposed species and one genomic species) suggested that the 5' end hypervariable region (259-273 bases in length) of 16S rDNA was specific for the respective type strains. Further phylogenetic analysis based on DNA sequences of the hypervariable region using 24 Alicyclobacillus strains revealed that the strains could be categorized into five species and the A. acidocaldarius-Alicyclobacillus genomic species 1 group. The hypervariable region was highly conserved among the five species: A. acidiphilus, A. acidoterrestris, A. cycloheptanicus, A. herbarius, and A. hesperidum. The strains in the A. acidocaldarius-Alicyclobacillus genomic species 1 group were subdivided into two clusters (Clusters I and II) based on DNA sequences in the hypervariable region. On the basis of phenotypic characteristics, chemotaxonomic and phylogenetic analyses, and DNA-DNA hybridization data, strains in Cluster I were grouped as Alicyclobacillus genomic species 1 and strains in Cluster II were re-identified as A. acidocaldarius, thereby demonstrating that the hypervariable regions were also highly conserved within these two species. These results suggest that as is the case with Bacillus, the hypervariable region is significantly species-specific in the genus Alicyclobacillus to distinguish Alicyclobacillus species by DNA sequence comparison of the hypervariable region.     

Relative rates of evolution in the coding and control regions of African mtDNAs

Reduced median networks of African haplogroup L mitochondrial DNA (mtDNA) sequences were analyzed to determine the pattern of substitutions in both the noncoding control and coding regions. In particular, we attempted to determine the causes of the previously reported (Howell et al. 2004) violation of the molecular clock during the evolution of these sequences. In the coding region, there was a significantly higher rate of substitution at synonymous sites than at nonsynonymous sites as well as in the tRNA and rRNA genes. This is further evidence for the operation of purifying selection during human mtDNA evolution. For most sites in the control region, the relative rate of substitution was similar to the rate of neutral evolution (assumed to be most closely approximated by the substitution rate at 4-fold degenerate sites). However, there are a number of mutational hot spots in the control region, approximately 3% of the total sites, that have a rate of substitution greater than the neutral rate, at some sites by more than an order of magnitude. It is possible either that these sites are evolving under conditions of positive selection or that the substitution rate at some sites in the control region is strongly dependent upon sequence context. Finally, we obtained preliminary evidence for "nonideal" evolution in the control region, including haplogroup-specific substitution patterns and a decoupling between relative rates of substitution in the control and coding regions.     

Design of bioactive peptides based on antibody hypervariable region structures. Development of conformationally constrained and dimeric peptides with enhanced affinity

The variable regions of antibody molecules bind antigens with high affinity and specificity. This binding is imparted largely by the hypervariable portions of the variable region. Hypervariable regions typically fold into reverse turn or loop structures. Peptides derived from antibody hypervariable region sequences can bind antigens with similar specificity, albeit with markedly lower affinity. In this study, cyclic and dimeric peptide analogs of an anti-idiotypic/antireceptor antibody hypervariable region were developed. This antibody (87.92.6) binds to reovirus type 3 receptors on cells as well as to a neutralizing anti-reovirus type 3 monoclonal antibody (9B.G5). The cyclic peptides were utilized to probe the optimal conformation for binding to both the receptor and 9B.G5. By dimerizing or constraining the conformation of these peptides, higher affinity binding was produced. By utilizing several different cyclic peptides, the optimal conformation for binding was established. The conformationally optimized cyclic peptide possessed greater than 40-fold higher affinity for the receptor and the idiotype than the linear analog. This study suggests that conformationally constrained and dimeric peptides derived from antibody hypervariable loop sequences can bind antigens (including receptors) with reasonable affinity. hypervariable loop sequences can bind antigens (including     

A naturally occurring single basic amino acid substitution in the V3 region of the human immunodeficiency virus type 1 env protein alters the cellular host range and antigenic structure of the virus.

Human immunodeficiency virus type 1 circulates in vivo as a mixture of heterologous populations (quasispecies). We previously analyzed the quasispecies of the third hypervariable region (V3) in the viral envelope glycoprotein gp120 in an infected individual and found that the species with a basic amino acid substitution (lysine for aspartic acid) at a particular position evolved and became a distinct population within a short period, followed by progression to the typical immunodeficiency stage (S. Oka et al., AIDS Res. Hum. Retroviruses 10:271-277, 1994). In the present study, we examined the biological significance of this amino acid substitution by constructing recombinant viruses with specific point mutations and comparing their replication capabilities in different cell types. The results demonstrated that the single basic amino acid substitution was sufficient to render a virus fully capable of replicating in human T-cell lines under certain conditions. With an acidic amino acid at the position, the virus grew much less fast or did not grow at all in the T-cell lines. Viral neutralization assay and peptide enzyme-linked immunosorbent assays further showed that this amino acid substitution resulted in different recognition by several of the serum specimens from human immunodeficiency virus type 1-infected individuals and thus could alter the antigenic structure. An additional finding worthy of note was that at the terminal stage, the proviral sequences of peripheral blood mononuclear cells and the viral isolates from them were without exception of the late type with the basic amino acid substitution, whereas the early sequence without the substitution was retained as a major subset in the spleen. These results support the notion that basic amino acid substitutions in V3 are a strong predictor of virus tropism and may be relevant to disease progression.     

South from Alaska: a pilot aDNA study of genetic history on the Alaska Peninsula and the eastern Aleutians

The Aleutian Islands were colonized, perhaps several times, from the Alaskan mainland. Earlier work documented transitions in the relative frequencies of mtDNA haplogroups over time, but little is known about potential source populations for prehistoric Aleut migrants. As part of a pilot investigation, we sequenced the mtDNA first hypervariable region (HVRI) in samples from two archaeological sites on the Alaska Peninsula (the Hot Springs site near Port Moller, Alaska; and samples from a cluster of sites in the Brooks River area near Katmai National Park and Preserve) and one site from Prince William Sound (Mink Island). The sequences revealed not only the mtDNA haplogroups typically found in both ancient and modern Aleut populations (A2 and D2) but also haplogroups B2 and D1 in the Brooks River samples and haplogroup D3 in one Mink Islander. These preliminary results suggest greater mtDNA diversity in prehistoric populations than previously observed and facilitate reconstruction of migration scenarios from the peninsula into the Aleutian archipelago in the past.     

Comparison of the amino acid sequences of the variable domains of two homogeneous rabbit antibodies to type III pneumococcal polysaccharide.

The amino acid sequences of the V (variable) regions of the H (heavy) and L (light) chains derived from rabbit antibody K-25, specific for type III pneumococci, were determined; this is the second homogeneous rabbit antibody besides antibody BS-5 whose complete sequence of the V domain has been established (Jaton, 1974d). The V regions of L chains BS-5 and K-25 (both of allotype b4) differ from each other by 19 amino acid residues; 11 of these 19 substitutions are located within the three hypervariable sections of the V region. On the basis of seven amino acid differences within the N-terminal 28 positions, it is suggested that L chain K-25 belongs to a different subgroup of rabbit K chains and L chain BS-5. H chain K-25 (allotype a2) differs from another H chain of the same allotype by one amino acid substitution within the N-terminal 70 positions in addition to interchanges occurring in the first two hypervariable sections. H chain K-25 was compared with H chain BS-5 (allotype a1) and with the known V-region rabbit sequences. Allotype-related differences between a1, a2 and a3 chains appear to occur within the N-terminal 16 positions and possibly in scattered positions throughout the V-region. In the hypervariable positions, variability between the two antibodies is remarkably more pronounced within the third hypervariable section of both H and L chains than within the first two.