High divergence among Drosophila simulans mitochondrial haplogroups arose in midst of long term purifying selection

We characterize the type of selection acting within and among mitochondrial lineages in five closely related Drosophila species. We focus on D. simulans, where three genetically distinct mitochondrial haplogroups show high interhaplogroup divergence and low intrahaplogroup polymorphism. Using maximum likelihood models we find that the branches leading to these three distinct mitochondrial groups show a significantly reduced rate of nonsynonymous relative to synonymous substitution. This interhaplogroup rate is significantly reduced compared to the intrahaplogroup rate, and closely resembles the rate observed between distinct species. The data suggest that slightly deleterious mutations segregating within D. simulans haplogroups are removed by selection prior to their fixation among haplogroups. We explore several hypotheses to explain how lineages within a single species can be compatible with this model of slightly deleterious mutation. The most likely hypothesis is that D. simulans haplogroups have persisted in isolation, perhaps due to association with the bacterial symbiont Wolbachia and/or demographic history, introducing a bias against the fixation of slightly deleterious mutations.     

Difference between evolutionarily effective and germ line mutation rate due to stochastically varying haplogroup size

Within a Y-chromosome haplogroup defined by unique event mutations, variation in microsatellites can accumulate due to their rapid mutation. Estimates based on pedigrees for the Y-chromosome microsatellite mutation rate are 3 or more times greater than the same estimates from evolutionary considerations. We show by simulation that the haplogroups that survive the stochastic processes of drift and extinction accumulate microsatellite variation at a lower rate than predicted from corresponding pedigree estimates; in particular, under constant total population size, the accumulated variance is on average 3-4 times smaller.     

Evolutionary rate variation at multiple levels of biological organization in plant mitochondrial DNA

We examined patterns of mitochondrial polymorphism and divergence in the angiosperm genus Silene and found substantial variation in evolutionary rates among species and among lineages within species. Moreover, we found corresponding differences in the amount of polymorphism within species. We argue that, along with our earlier findings of rate variation among genes, these patterns of rate heterogeneity at multiple phylogenetic scales are most likely explained by differences in underlying mutation rates. In contrast, no rate variation was detected in nuclear or chloroplast loci. We conclude that mutation rate heterogeneity is a characteristic of plant mitochondrial sequence evolution at multiple biological scales and may be a crucial determinant of how much polymorphism is maintained within species. These dramatic patterns of variation raise intriguing questions about the mechanisms driving and maintaining mutation rate heterogeneity in plant mitochondrial genomes. Additionally, they should alter our interpretation of many common phylogenetic and population genetic analyses.     

Refining the Y chromosome phylogeny with southern African sequences

The recent availability of large-scale sequence data for the human Y chromosome has revolutionized analyses of and insights gained from this non-recombining, paternally inherited chromosome. However, the studies to date focus on Eurasian variation, and hence the diversity of early-diverging branches found in Africa has not been adequately documented. Here, we analyze over 900 kb of Y chromosome sequence obtained from 547 individuals from southern African Khoisan- and Bantu-speaking populations, identifying 232 new sequences from basal haplogroups A and B. We identify new clades in the phylogeny, an older age for the root, and substantially older ages for some individual haplogroups. Furthermore, while haplogroup B2a is traditionally associated with the spread of Bantu speakers, we find that it probably also existed in Khoisan groups before the arrival of Bantu speakers. Finally, there is pronounced variation in branch length between major haplogroups; in particular, haplogroups associated with Bantu speakers have significantly longer branches. Technical artifacts cannot explain this branch length variation, which instead likely reflects aspects of the demographic history of Bantu speakers, such as recent population expansion and an older average paternal age. The influence of demographic factors on branch length variation has broader implications both for the human Y phylogeny and for similar analyses of other species.     

Factors affecting germline mutations in a hypervariable microsatellite: a comparative analysis of six species of swallows (Aves: Hirundinidae)

Microsatellites mutate frequently by replication slippage. Empirical evidence shows that the probability of such slippage mutations may increase with the length of the repeat region as well as exposure to environmental mutagens, but the mutation rate can also differ between the male and female germline. It has been hypothesized that more intense sexual selection or sperm competition can also lead to elevated mutation rates, but the empirical evidence is inconclusive. Here, we analyzed the occurrence of germline slippage mutations in the hypervariable pentanucleotide microsatellite locus HrU10 across six species of swallow (Aves: Hirundinidae). These species exhibit marked differences in the length range of the microsatellite, as well as differences in the intensity of sperm competition. We found a strong effect of microsatellite length on the probability of mutation, but no residual effect of species or their level of sperm competition when the length effect was accounted for. Neither could we detect any difference in mutation rate between tree swallows (Tachycineta bicolor) breeding in Hamilton Harbour, Ontario, an industrial site with previous documentation of elevated mutation rates for minisatellite DNA, and a rural reference population. However, our cross-species analysis revealed two significant patterns of sex differences in HrU10 germline mutations: (1) mutations in longer alleles occurred typically in the male germline, those in shorter alleles in the female germline, and (2) male germline mutations were more often expansions than contractions, whereas no directional bias was evident in the female germline. These results indicate some fundamental differences in male and female gametogenesis affecting the probability of slippage mutations. Our study also reflects the value of a comparative, multi-species approach for locus-specific mutation analyses, through which a wider range of influential factors can be assessed than in single-species studies.     

Quantification of the paternal allele bias for new germline mutations in the retinoblastoma gene

New germline mutations in the human retinoblastoma gene are known to arise preferentially on paternally derived chromosomes, but the magnitude of that bias has not been measured. We evaluated 49 cases with a new germline mutation and found that in 40 cases (82%) the mutation arose on the paternally derived allele. We also evaluated 48 cases likely to have a somatic initial mutation; in this group the initial mutation arose on paternal or maternal chromosomes with approximately equal frequency. There was no statistically significant difference in the average age of fathers of children with new paternal germline mutations from the average age of fathers of children with new maternal germline mutations or somatic initial mutations. Combining the data with that from previous reports from other groups, the proportion of new germline mutations arising on a paternally derived allele is 85% (based on 72 cases; 95% confidence interval = 76–93%). This number can be useful in the genetic counseling of some families with retinoblastoma. Received: 18 December 1996 / Accepted: 30 April 1997     

人类线粒体DNA变异的检测方法和思路

基于线粒体DNA（mtDNA)的研究对于人群源流迁移、线粒体相关疾病病因的探讨和法医鉴定等具有重意义,就检测人线粒体突变的一些常用方法,如RFLP、SSO和控制区测序等作一小结和归纳,并重点介绍目前mtDNA突变的筛选方法和思路,另外,还总结了近年来对人mtDNA方面的研究结果,对世界人群中主要单倍型类群（haplogroup)特征变异位点和相应的酶切检测引物作了归纳。     

Differentiation of mitochondrial DNA and Y chromosomes in Russian populations

The genetic composition of the Russian population was investigated by analyzing both mitochondrial DNA (mtDNA) and Y-chromosome loci polymorphisms that allow for the different components of a population gene pool to be studied, depending on the mode of DNA marker inheritance. mtDNA sequence variation was examined by using hypervariable segment I (HVSI) sequencing and restriction analysis of the haplogroup-specific sites in 325 individuals representing 5 Russian populations from the European part of Russia. The Y-chromosome variation was investigated in 338 individuals from 8 Russian populations (including 5 populations analyzed for mtDNA variation) using 12 binary markers. For both uniparental systems most of the observed haplogroups fell into major West Eurasian haplogroups (97.9% and 99.7% for mtDNA and Y-chromosome haplogroups, respectively). Multidimensional scaling analysis based on pairwise F(ST) values between mtDNA HVSI sequences in Russians compared to other European populations revealed a considerable heterogeneity of Russian populations; populations from the southern and western parts of Russia are separated from eastern and northern populations. Meanwhile, the multidimensional scaling analysis based on Y-chromosome haplogroup F(ST) values demonstrates that the Russian gene pool is close to central-eastern European populations, with a much higher similarity to the Baltic and Finno-Ugric male pools from northern European Russia. This discrepancy in the depth of penetration of mtDNA and Y-chromosome lineages characteristic for the most southwestern Russian populations into the east and north of eastern Europe appears to indicate that Russian colonization of the northeastern territories might have been accomplished mainly by males rather than by females.     

No recombination of mtDNA after heteroplasmy for 50 generations in the mouse maternal germline

Variants of mitochondrial DNA (mtDNA) are commonly used as markers to track human evolution because of the high sequence divergence and exclusive maternal inheritance. It is assumed that the inheritance is clonal, i.e. that mtDNA is transmitted between generations without germline recombination. In contrast to this assumption, a number of studies have reported the presence of recombinant mtDNA molecules in cell lines and animal tissues, including humans. If germline recombination of mtDNA is frequent, it would strongly impact phylogenetic and population studies by altering estimates of coalescent time and branch lengths in phylogenetic trees. Unfortunately, this whole area is controversial and the experimental approaches have been widely criticized as they often depend on polymerase chain reaction (PCR) amplification of mtDNA and/or involve studies of transformed cell lines. In this study, we used an in vivo mouse model that has had germline heteroplasmy for a defined set of mtDNA mutations for more than 50 generations. To assess recombination, we adapted and validated a method based on cloning of single mtDNA molecules in the λ phage, without prior PCR amplification, followed by subsequent mutation analysis. We screened 2922 mtDNA molecules and found no germline recombination after transmission of mtDNA under genetically and evolutionary relevant conditions in mammals.     

Germline mutations of tetranucleotide DNA repeats in families with normal children and reproductive pathology

We have previously reported a high rate of tetranucleotide DNA repeat mutations, including mutations of both germline and somatic origin, in spontaneous human abortuses. To analyze in more detail mutational microsatellite (MS) variability in meiosis and its possible association with disturbed embryonic development, we have conducted a comparative study of mutation rates of a complex of 15 autosomal tetranucleotide MSs in 55 families with healthy children and in 103 families that have had spontaneous abortuses with normal karyotypes. In the families with miscarriage, the gametic MS mutation rate was higher than in the families with normal reproductive function (4.36 x 10(-3) versus 2.32 x 10(-3) per locus per gamete per generation), but this difference was statistically nonsignificant (P = 0.25). No association of MS mutations with familiar miscarriage was found. Mutations at the MS loci studied were recorded almost 3 times as often in spermatogenesis as in oogenesis, which is likely to result from a greater number of DNA replication cycles in male germline cell precursors than in female ones. Mutations increasing and reducing the MS sequence length appeared at virtually the same rate. Changes in MS DNA sequence length per one repeated element, i.e., single-step mutations (93% of cases) exceeded all other events of allele length change. The highest number of mutations (81.2%) was found in longer alleles. This distribution of mutations by size, direction, and parental origin corresponds to the multistep mutation model of their emergence via mechanism of DNA strand slippage during replication.     

Mutations in the calcium-binding motifs of CDH23 and the 35delG mutation in GJB2 cause hearing loss in one family

We have ascertained a multi-generation family with apparent autosomal recessive non-syndromic childhood hearing loss (DFNB). Failure to demonstrate linkage in a genome-wide scan with 300 polymorphic markers has suggested genetic heterogeneity for the hearing loss in this family. This heterogeneity could be demonstrated by analysis of candidate loci and genes for DFNB. Patients in one branch of the family (branch C) are homozygous for the 35delG mutation in the GJB2 gene (DFNB1). Patients in two other branches (A and B) carry two new mutations in the cadherin 23 ( CDH23) gene (DFNB12). A homozygous CDH23 c.6442G-->A (D2148N) mutation is present in branch A. Patients in branch B are compound heterozygous for this mutation and the c.4021G-->A (D1341N) mutation. The substituted aspartic acid residues are highly conserved and are part of the calcium-binding sites of the extracellular cadherin (EC) domains. Molecular modeling of the mutated EC domains of CDH23 based on the structure of E-cadherin indicates that calcium-binding is impaired. In addition, other aspartic and glutamic acid residue substitutions in the highly conserved calcium-binding sites reported to cause DFNB12 are also likely to result in a decreased affinity for calcium. Since calcium provides rigidity to the elongated structure of cadherin molecules enabling homophilic lateral interaction, these mutations are likely to impair interactions of CDH23 molecules either with CDH23 or with other proteins. DFNB12 is the first human disorder that can be attributed to inherited missense mutations in the highly conserved residues of the extracellular calcium-binding domain of a cadherin.     

Germline Mutations of Tetranucleotide DNA Repeats in Families with Normal Children and Reproductive Pathology

We have previously reported a high rate of tetranucleotide DNA repeat mutations, including mutations of both germline and somatic origin, in spontaneous human abortions. To analyze in more detail mutational microsatellite (MS) variability in meiosis and its possible association with disturbed embryonic development, we have conducted a comparative study of mutation rates of a panel of 15 autosomal tetranucleotide MSs in 55 families with healthy children and in 103 families that have had spontaneous abortions with normal karyotypes. In the families with miscarriage, the gametic MS mutation rate was higher than in the families with normal reproductive function (4.36 × 10⁻³ versus 2.32 × 10⁻³ per locus per gamete per generation), but this difference was statistically nonsignificant (P = 0.25). No association of MS mutations with familiar miscarriage was found. Mutations at the MS loci studied were recorded almost 3 times as often in spermatogenesis as in oogenesis, which is likely to result from a greater number of DNA replication cycles in male germline cell precursors than in female ones. Mutations increasing and reducing the MS sequence length appeared at virtually the same rate. Changes in MS DNA sequence length per one repeated element, i.e., single-step mutations (93% of cases) exceeded all other events of allele length change. The highest number of mutations (81.2%) was found in longer alleles. This distribution of mutations by size, direction, and parental origin corresponds to the multistep mutation model of their emergence via mechanism of DNA strand slippage during replication.__________Translated from Genetika, Vol. 41, No. 7, 2005, pp. 943–953.Original Russian Text Copyright © 2005 by Nikitina, Lebedev, Sukhanova, Nazarenko.     

The pedigree rate of sequence divergence in the human mitochondrial genome: there is a difference between phylogenetic and pedigree rates

We have extended our previous analysis of the pedigree rate of control-region divergence in the human mitochondrial genome. One new germline mutation in the mitochondrial DNA (mtDNA) control region was detected among 185 transmission events (generations) from five Leber hereditary optic neuropathy (LHON) pedigrees. Pooling the LHON pedigree analyses yields a control-region divergence rate of 1.0 mutation/bp/10(6) years (Myr). When the results from eight published studies that used a similar approach were pooled with the LHON pedigree studies, totaling >2,600 transmission events, a pedigree divergence rate of 0.95 mutations/bp/Myr for the control region was obtained with a 99.5% confidence interval of 0.53-1.57. Taken together, the cumulative results support the original conclusion that the pedigree divergence rate for the control region is approximately 10-fold higher than that obtained with phylogenetic analyses. There is no evidence that any one factor explains this discrepancy, and the possible roles of mutational hotspots (rate heterogeneity), selection, and random genetic drift and the limitations of phylogenetic approaches to deal with high levels of homoplasy are discussed. In addition, we have extended our pedigree analysis of divergence in the mtDNA coding region. Finally, divergence of complete mtDNA sequences was analyzed in two tissues, white blood cells and skeletal muscle, from each of 17 individuals. In three of these individuals, there were four instances in which an mtDNA mutation was found in one tissue but not in the other. These results are discussed in terms of the occurrence of somatic mtDNA mutations.     

Germline fitness-based scoring of cancer mutations

A key goal in cancer research is to find the genomic alterations that underlie malignant cells. Genomics has proved successful in identifying somatic variants at a large scale. However, it has become evident that a typical cancer exhibits a heterogenous mutation pattern across samples. Cases where the same alteration is observed repeatedly seem to be the exception rather than the norm. Thus, pinpointing the key alterations (driver mutations) from a background of variations with no direct causal link to cancer (passenger mutations) is difficult. Here we analyze somatic missense mutations from cancer samples and their healthy tissue counterparts (germline mutations) from the viewpoint of germline fitness. We calibrate a scoring system from protein domain alignments to score mutations and their target loci. We show first that this score predicts to a good degree the rate of polymorphism of the observed germline variation. The scoring is then applied to somatic mutations. We show that candidate cancer genes prone to copy number loss harbor mutations with germline fitness effects that are significantly more deleterious than expected by chance. This suggests that missense mutations play a driving role in tumor suppressor genes. Furthermore, these mutations fall preferably onto loci in sequence neighborhoods that are high scoring in terms of germline fitness. In contrast, for somatic mutations in candidate onco genes we do not observe a statistically significant effect. These results help to inform how to exploit germline fitness predictions in discovering new genes and mutations responsible for cancer.     

Polarity and temporality of high-resolution y-chromosome distributions in India identify both indigenous and exogenous expansions and reveal minor genetic influence of Central Asian pastoralists

Although considerable cultural impact on social hierarchy and language in South Asia is attributable to the arrival of nomadic Central Asian pastoralists, genetic data (mitochondrial and Y chromosomal) have yielded dramatically conflicting inferences on the genetic origins of tribes and castes of South Asia. We sought to resolve this conflict, using high-resolution data on 69 informative Y-chromosome binary markers and 10 microsatellite markers from a large set of geographically, socially, and linguistically representative ethnic groups of South Asia. We found that the influence of Central Asia on the pre-existing gene pool was minor. The ages of accumulated microsatellite variation in the majority of Indian haplogroups exceed 10,000-15,000 years, which attests to the antiquity of regional differentiation. Therefore, our data do not support models that invoke a pronounced recent genetic input from Central Asia to explain the observed genetic variation in South Asia. R1a1 and R2 haplogroups indicate demographic complexity that is inconsistent with a recent single history. Associated microsatellite analyses of the high-frequency R1a1 haplogroup chromosomes indicate independent recent histories of the Indus Valley and the peninsular Indian region. Our data are also more consistent with a peninsular origin of Dravidian speakers than a source with proximity to the Indus and with significant genetic input resulting from demic diffusion associated with agriculture. Our results underscore the importance of marker ascertainment for distinguishing phylogenetic terminal branches from basal nodes when attributing ancestral composition and temporality to either indigenous or exogenous sources. Our reappraisal indicates that pre-Holocene and Holocene-era--not Indo-European--expansions have shaped the distinctive South Asian Y-chromosome landscape.     

The K tree score: quantification of differences in the relative branch length and topology of phylogenetic trees

SUMMARY: We introduce a new phylogenetic comparison method that measures overall differences in the relative branch length and topology of two phylogenetic trees. To do this, the algorithm first scales one of the trees to have a global divergence as similar as possible to the other tree. Then, the branch length distance, which takes differences in topology and branch lengths into account, is applied to the two trees. We thus obtain the minimum branch length distance or K tree score. Two trees with very different relative branch lengths get a high K score whereas two trees that follow a similar among-lineage rate variation get a low score, regardless of the overall rates in both trees. There are several applications of the K tree score, two of which are explained here in more detail. First, this score allows the evaluation of the performance of phylogenetic algorithms, not only with respect to their topological accuracy, but also with respect to the reproduction of a given branch length variation. In a second example, we show how the K score allows the selection of orthologous genes by choosing those that better follow the overall shape of a given reference tree. AVAILABILITY: http://molevol.ibmb.csic.es/Ktreedist.html     

Patterns of Y-chromosome diversity intersect with the Trans-New Guinea hypothesis

The island of New Guinea received part of the first human expansion out of Africa (>40,000 years ago), but its human genetic history remains poorly understood. In this study, we examined Y-chromosome diversity in 162 samples from the Bird's Head region of northwest New Guinea (NWNG) and compared the results with previously obtained data from other parts of the island. NWNG harbors a high level of cultural and linguistic diversity and is inhabited by non-Austronesian (i.e., Papuan)-speaking groups as well as harboring most of West New Guinea's (WNG) Austronesian-speaking groups. However, 97.5% of its Y-chromosomes belong to 5 haplogroups that originated in Melanesia; hence, the Y-chromosome diversity of NWNG (and, according to available data, of New Guinea as a whole) essentially reflects a local history. The remaining 2.5% belong to 2 haplogroups (O-M119 and O-M122) of East Asian origin, which were brought to New Guinea by Austronesian-speaking migrants around 3,500 years ago. Thus, the Austronesian expansion had only a small impact on shaping Y-chromosome diversity in NWNG, although the linguistic impact of this expansion to this region was much higher. In contrast, the expansion of Trans-New Guinea (TNG) speakers (non-Austronesian) starting about 6,000-10,000 years ago from the central highlands of what is now Papua New Guinea, presumably in combination with the expansion of agriculture, played a more important role in determining the Y-chromosome diversity of New Guinea. In particular, we identified 2 haplogroups (M-P34 and K-M254) as suggestive markers for the TNG expansion, whereas 2 other haplogroups (C-M38 and K-M9) most likely reflect the earlier local Y-chromosome diversity. We propose that sex-biased differences in the social structure and cultural heritage of the people involved in the Austronesian and the TNG expansions played an important role (among other factors) in shaping the New Guinean Y-chromosome landscape.     

Heterogeneity in viral populations increases the rate of deleterious mutation accumulation

RNA viruses have high mutation rates, with the majority of mutations being deleterious. We examine patterns of deleterious mutation accumulation over multiple rounds of viral replication, with a focus on how cellular coinfection and heterogeneity in viral output affect these patterns. Specifically, using agent-based intercellular simulations we find, in agreement with previous studies, that coinfection of cells by viruses relaxes the strength of purifying selection and thereby increases the rate of deleterious mutation accumulation. We further find that cellular heterogeneity in viral output exacerbates the rate of deleterious mutation accumulation, regardless of whether this heterogeneity in viral output is stochastic or is due to variation in the cellular multiplicity of infection. These results highlight the need to consider the unique life histories of viruses and their population structure to better understand observed patterns of viral evolution.     

Rate and pattern of mutation at microsatellite loci in maize

Microsatellites are important tools for plant breeding, genetics, and evolution, but few studies have analyzed their mutation pattern in plants. In this study, we estimated the mutation rate for 142 microsatellite loci in maize (Zea mays subsp. mays) in two different experiments of mutation accumulation. The mutation rate per generation was estimated to be 7.7 x 10(-4) for microsatellites with dinucleotide repeat motifs, with a 95% confidence interval from 5.2 x 10(-4) to 1.1 x 10(-3). For microsatellites with repeat motifs of more than 2 bp in length, no mutations were detected; so we could only estimate the upper 95% confidence limit of 5.1 x 10(-5) for the mutation rate. For dinucleotide repeat microsatellites, we also determined that the variance of change in the number of repeats (sigma(m)2) is 3.2. We sequenced 55 of the 73 observed mutations, and all mutations proved to be changes in the number of repeats in the microsatellite or in mononucleotide tracts flanking the microsatellite. There is a higher probability to mutate to an allele of larger size. There is heterogeneity in the mutation rate among dinucleotide microsatellites and a positive correlation between the number of repeats in the progenitor allele and the mutation rate. The microsatellite-based estimate of the effective population size of maize is more than an order of magnitude less than previously reported values based on nucleotide sequence variation.