Linkage disequilibrium and recombination in hominid mitochondrial DNA

Epistatic selection on a selfish Segregation Distorter supergene - drive, recombination, and genetic load

Meiotic drive supergenes are complexes of alleles at linked loci that together subvert Mendelian segregation resulting in preferential transmission. In males, the most common mechanism of drive involves the disruption of sperm bearing one of a pair of alternative alleles. While at least two loci are important for male drive-the driver and the target-linked modifiers can enhance drive, creating selection pressure to suppress recombination. In this work, we investigate the evolution and genomic consequences of an autosomal, multilocus, male meiotic drive system, Segregation Distorter (SD) in the fruit fly, Drosophila melanogaster. In African populations, the predominant SD chromosome variant, SD-Mal, is characterized by two overlapping, paracentric inversions on chromosome arm 2R and nearly perfect (~100%) transmission. We study the SD-Mal system in detail, exploring its components, chromosomal structure, and evolutionary history. Our findings reveal a recent chromosome-scale selective sweep mediated by strong epistatic selection for haplotypes carrying Sd, the main driving allele, and one or more factors within the double inversion. While most SD-Mal chromosomes are homozygous lethal, SD-Mal haplotypes can recombine with other, complementing haplotypes via crossing over, and with wildtype chromosomes via gene conversion. SD-Mal chromosomes have nevertheless accumulated lethal mutations, excess non-synonymous mutations, and excess transposable element insertions. Therefore, SD-Mal haplotypes evolve as a small, semi-isolated subpopulation with a history of strong selection. These results may explain the evolutionary turnover of SD haplotypes in different populations around the world and have implications for supergene evolution broadly.

Integrated Taxonomy Revealed Genetic Differences in Morphologically Similar and Non-Sympatric Scoliodon macrorhynchos and S. laticaudus

Simple Summary

In this study, the species identities of similar-looking coastal spadenose sharks from different areas were clarified by adding new molecular markers and more individual body measurements, including animals from the Malaysian Peninsula that had not been examined previously. Collective evidence showed that there are two genetically distinct species that do not overlap in their spatial occurrence. The Malacca Strait acts as a boundary delineating the distribution range of the Pacific spadenose shark Scoliodon macrorhynchos to the east and, of the Northern Indian Ocean, S. laticaudus to the west. In addition, the need to determine the species status of Scoliodon animals from Indonesian waters was identified. The present study reinforced the need to rely on comprehensive genetic information in addition to external characteristics to assess the species identities and distribution range for small sharks and rays that have apparent contiguous coastal distribution with limited dispersal abilities.

Abstract

Previous examination of the mitochondrial NADH2 gene and morphological characteristics led to the resurrection of Scoliodon macrorhynchos as a second valid species in the genus, in addition to S. laticaudus. This study applied an integrated taxonomic approach to revisit the classification of the genus Scoliodon based on new materials from the Malaysian Peninsula, Malaysian Borneo and Eastern Bay of Bengal. Mitochondrial DNA data suggested the possibility of three species of Scoliodon in the Indo-West Pacific, while the nuclear DNA data showed partially concordant results with a monophyletic clade of S. macrorhynchos and paraphyletic clades of S. laticaudus and S. cf. laticaudus from the Malacca Strait. Morphological, meristic and dental characteristics overlapped between the three putative species. Collective molecular and morphological evidence suggested that the differences that exist among the non-sympatric species of Scoliodon are consistent with isolation by distance, and Scoliodon macrorhynchos remains as a valid species, while S. cf. laticaudus is assigned as S. laticaudus. The Malacca Strait acts as a spatial delineator in separating the Pacific S. macrorhynchos (including South China Sea) from the Northern Indian Ocean S. laticaudus. Future taxonomic work should focus on clarifying the taxonomic status of Scoliodon from the Indonesian waters.

Global Patterns of Recombination across Human Viruses

Viral recombination is a major evolutionary mechanism driving adaptation processes, such as the ability of host-switching. Understanding global patterns of recombination could help to identify underlying mechanisms and to evaluate the potential risks of rapid adaptation. Conventional approaches (e.g., those based on linkage disequilibrium) are computationally demanding or even intractable when sequence alignments include hundreds of sequences, common in viral data sets. We present a comprehensive analysis of recombination across 30 genomic alignments from viruses infecting humans. In order to scale the analysis and avoid the computational limitations of conventional approaches, we apply newly developed topological data analysis methods able to infer recombination rates for large data sets. We show that viruses, such as ZEBOV and MARV, consistently displayed low levels of recombination, whereas high levels of recombination were observed in Sarbecoviruses, HBV, HEV, Rhinovirus A, and HIV. We observe that recombination is more common in positive single-stranded RNA viruses than in negatively single-stranded RNA ones. Interestingly, the comparison across multiple viruses suggests an inverse correlation between genome length and recombination rate. Positional analyses of recombination breakpoints along viral genomes, combined with our approach, detected at least 39 nonuniform patterns of recombination (i.e., cold or hotspots) in 18 viral groups. Among these, noteworthy hotspots are found in MERS-CoV and Sarbecoviruses (at spike, Nucleocapsid and ORF8). In summary, we have developed a fast pipeline to measure recombination that, combined with other approaches, has allowed us to find both common and lineage-specific patterns of recombination among viruses with potential relevance in viral adaptation.

Evidence for multi-copy Mega-NUMTs in the human genome

The maternal mode of mitochondrial DNA (mtDNA) inheritance is central to human genetics. Recently, evidence for bi-parental inheritance of mtDNA was claimed for individuals of three pedigrees that suffered mitochondrial disorders. We sequenced mtDNA using both direct Sanger and Massively Parallel Sequencing in several tissues of eleven maternally related and other affiliated healthy individuals of a family pedigree and observed mixed mitotypes in eight individuals. Cells without nuclear DNA, i.e. thrombocytes and hair shafts, only showed the mitotype of haplogroup (hg) V. Skin biopsies were prepared to generate ρ° cells void of mtDNA, sequencing of which resulted in a hg U4c1 mitotype. The position of the Mega-NUMT sequence was determined by fluorescence in situ hybridization and two different quantitative PCR assays were used to determine the number of contributing mtDNA copies. Thus, evidence for the presence of repetitive, full mitogenome Mega-NUMTs matching haplogroup U4c1 in various tissues of eight maternally related individuals was provided. Multi-copy Mega-NUMTs mimic mixtures of mtDNA that cannot be experimentally avoided and thus may appear in diverse fields of mtDNA research and diagnostics. We demonstrate that hair shaft mtDNA sequencing provides a simple but reliable approach to exclude NUMTs as source of misleading results.

Genomic signatures of recombination in a natural population of the bdelloid rotifer Adineta vaga

Sexual reproduction is almost ubiquitous among extant eukaryotes. As most asexual lineages are short-lived, abandoning sex is commonly regarded as an evolutionary dead end. Still, putative anciently asexual lineages challenge this view. One of the most striking examples are bdelloid rotifers, microscopic freshwater invertebrates believed to have completely abandoned sexual reproduction tens of Myr ago. Here, we compare whole genomes of 11 wild-caught individuals of the bdelloid rotifer Adineta vaga and present evidence that some patterns in its genetic variation are incompatible with strict clonality and lack of genetic exchange. These patterns include genotype proportions close to Hardy-Weinberg expectations within loci, lack of linkage disequilibrium between distant loci, incongruent haplotype phylogenies across the genome, and evidence for hybridization between divergent lineages. Analysis of triallelic sites independently corroborates these findings. Our results provide evidence for interindividual genetic exchange and recombination in A. vaga, a species previously thought to be anciently asexual.

Nuclear-mitochondrial DNA segments resemble paternally inherited mitochondrial DNA in humans

Several strands of evidence question the dogma that human mitochondrial DNA (mtDNA) is inherited exclusively down the maternal line, most recently in three families where several individuals harbored a 'heteroplasmic haplotype' consistent with biparental transmission. Here we report a similar genetic signature in 7 of 11,035 trios, with allelic fractions of 5-25%, implying biparental inheritance of mtDNA in 0.06% of offspring. However, analysing the nuclear whole genome sequence, we observe likely large rare or unique nuclear-mitochondrial DNA segments (mega-NUMTs) transmitted from the father in all 7 families. Independently detecting mega-NUMTs in 0.13% of fathers, we see autosomal transmission of the haplotype. Finally, we show the haplotype allele fraction can be explained by complex concatenated mtDNA-derived sequences rearranged within the nuclear genome. We conclude that rare cryptic mega-NUMTs can resemble paternally mtDNA heteroplasmy, but find no evidence of paternal transmission of mtDNA in humans.

Identity determination in diagnostic surgical pathology

From a technical perspective, specimen identity determination in surgical pathology over the last several decades has primarily focused on analysis of repetitive DNA sequences, specifically microsatellite repeats. However, a number of techniques have recently been developed that have similar, if not greater, utility in surgical pathology, most notably analysis of single nucleotide polymorphism (SNPs) and gene panels by next generation sequencing (NGS). For cases with an extremely limited sample or a degraded sample, sequence analysis of mitochondrial DNA continues to be the method of choice. From a diagnostic perspective, interest in identity determination in surgical pathology is usually centered on resolving issues of specimen provenance due to specimen labeling/accessioning deficiencies and possible contamination, but is also frequently performed in cases for which the patient's clinical course following definitive therapy is remarkably atypical, in cases of an unexpected diagnosis, and by patient request for "peace of mind". However, the methods used for identity determination have a much broader range of applications in surgical pathology beyond tissue provenance analysis. The methods can be used to provide ancillary information for cases in which the histomorphology is not definitively diagnostic, as for example for tumors that have a virtually identical microscopic appearance but for which the differential diagnosis includes synchronous/metachronous tumors versus a metastasis, and for the diagnosis of hydropic early gestations versus hydatidiform molar pregnancies. The methods also have utility in several other clinical settings, for example to rule out a donor-transmitted malignancy in a transplant recipient, to monitor bone marrow transplant engraftment, and to evaluate natural chimerism.

Clonorchis sinensis and Clonorchiasis: The Relevance of Exploring Genetic Variation

Parasitic trematodes (flukes) cause substantial mortality and morbidity in humans. The Chinese liver fluke, Clonorchis sinensis, is one of the most destructive parasitic worms in humans in China, Vietnam, Korea and the Russian Far East. Although C. sinensis infection can be controlled relatively well using anthelmintics, the worm is carcinogenic, inducing cholangiocarcinoma and causing major suffering in ~15 million people in Asia. This chapter provides an account of C. sinensis and clonorchiasis research-covering aspects of biology, epidemiology, pathogenesis and immunity, diagnosis, treatment and control, genetics and genomics. It also describes progress in the area of molecular biology (genetics, genomics, transcriptomics and proteomics) and highlights challenges associated with comparative genomics and population genetics. It then reviews recent advances in the sequencing and characterisation of the mitochondrial and nuclear genomes for a Korean isolate of C. sinensis and summarises salient comparative genomic work and the implications thereof. The chapter concludes by considering how advances in genomic and informatics will enable research on the genetics of C. sinensis and related parasites, as well as the discovery of new fluke-specific intervention targets.

An earlier revolution: genetic and genomic analyses reveal pre-existing cultural differences leading to Neolithization

Archaeological evidence shows that, in the long run, Neolitization (the transition from foraging to food production) was associated with demographic growth. We used two methods (patterns of linkage disequilibrium from whole-genome SNPs and MSMC estimates on genomes) to reconstruct the demographic profiles for respectively 64 and 24 modern-day populations with contrasting lifestyles across the Old World (sub-Saharan Africa, south-eastern Asia, Siberia). Surprisingly, in all regions, food producers had larger effective population sizes (N_e) than foragers already 20 k years ago, well before the Neolithic revolution. As expected, this difference further increased ~12–10 k years ago, around or just before the onset of food production. Using paleoclimate reconstructions, we show that the early difference in N_e cannot be explained by food producers inhabiting more favorable regions. A number of mechanisms, including ancestral differences in census size, sedentism, exploitation of the natural resources, social stratification or connectivity between groups, might have led to the early differences in Ne detected in our analyses. Irrespective of the specific mechanisms involved, our results provide further evidence that long term cultural differences among populations of Palaeolithic hunter-gatherers are likely to have played an important role in the later Neolithization process.

Pandora's Box: mitochondrial defects in ischaemic heart disease and stroke

Ischaemic heart disease and stroke are vascular events with serious health consequences worldwide. Recent genetic and epigenetic techniques have revealed many genetic determinants of these vascular events and simplified the approaches to research focused on ischaemic heart disease and stroke. The pathogenetic mechanisms of ischaemic heart disease and stroke are complex, with mitochondrial involvement (partially or entirely) recently gaining substantial support. Not only can mitochondrial reactive oxygen species give rise to ischaemic heart disease and stroke by production of oxidised low-density lipoprotein and induction of apoptosis, but the impact on pericytes contributes directly to the pathogenesis. Over the past two decades, publications implicate the causative role of nuclear genes in the development of ischaemic heart disease and stroke, in contrast to the potential role of mitochondrial DNA (mtDNA) in the pathophysiology of the disorders, which is much less understood, although recent studies do demonstrate that the involvement of mitochondria and mtDNA in the development of ischaemic heart disease and stroke is likely to be larger than originally thought, with the novel discovery of links among mitochondria, mtDNA and vascular events. Here we explore the molecular events and mtDNA alterations in relation to the role of mitochondria in ischaemic heart disease and stroke.

Mixed transmission modes and dynamic genome evolution in an obligate animal-bacterial symbiosis

Reliable transmission of symbionts between host generations facilitates the evolution of beneficial and pathogenic associations. Although transmission mode is typically characterized as either vertical or horizontal, the prevalence of intermediate transmission modes, and their impact on symbiont genome evolution, are understudied. Here, we use population genomics to explore mixed transmission modes of chemosynthetic bacterial symbionts in the bivalve Solemya velum. Despite strong evidence for symbiont inheritance through host oocytes, whole-genome analyses revealed signatures of frequent horizontal transmission, including discordant mitochondrial-symbiont genealogies, widespread recombination and a dynamic symbiont genome structure consistent with evolutionary patterns of horizontally transmitted associations. Population-level analyses thus provide a tractable means of ascertaining the fidelity of vertical versus horizontal transmission. Our data support the strong influence horizontal transmission can have on symbiont genome evolution, and shed light on the dynamic evolutionary pressures shaping symbiotic bacterial genomes.

Transcriptional quiescence of paternal mtDNA in cyprinid fish embryos

Mitochondrial homoplasmy signifies the existence of identical copies of mitochondrial DNA (mtDNA) and is essential for normal development, as heteroplasmy causes abnormal development and diseases in human. Homoplasmy in many organisms is ensured by maternal mtDNA inheritance through either absence of paternal mtDNA delivery or early elimination of paternal mtDNA. However, whether paternal mtDNA is transcribed has remained unknown. Here we report that paternal mtDNA shows late elimination and transcriptional quiescence in cyprinid fishes. Paternal mtDNA was present in zygotes but absent in larvae and adult organs of goldfish and blunt-snout bream, demonstrating paternal mtDNA delivery and elimination for maternal mtDNA inheritance. Surprisingly, paternal mtDNA remained detectable up to the heartbeat stage, suggesting its late elimination leading to embryonic heteroplasmy up to advanced embryogenesis. Most importantly, we never detected the cytb RNA of paternal mtDNA at all stages when paternal mtDNA was easily detectable, which reveals that paternal mtDNA is transcriptionally quiescent and thus excludes its effect on the development of heteroplasmic embryos. Therefore, paternal mtDNA in cyprinids shows late elimination and transcriptional quiescence. Clearly, transcriptional quiescence of paternal mtDNA represents a new mechanism for maternal mtDNA inheritance and provides implications for treating mitochondrion-associated diseases by mitochondrial transfer or replacement.

Similar Efficacies of Selection Shape Mitochondrial and Nuclear Genes in Both Drosophila melanogaster and Homo sapiens

Deleterious mutations contribute to polymorphism even when selection effectively prevents their fixation. The efficacy of selection in removing deleterious mitochondrial mutations from populations depends on the effective population size (N_e) of the mitochondrial DNA and the degree to which a lack of recombination magnifies the effects of linked selection. Using complete mitochondrial genomes from Drosophila melanogaster and nuclear data available from the same samples, we reexamine the hypothesis that nonrecombining animal mitochondrial DNA harbor an excess of deleterious polymorphisms relative to the nuclear genome. We find no evidence of recombination in the mitochondrial genome, and the much-reduced level of mitochondrial synonymous polymorphism relative to nuclear genes is consistent with a reduction in N_e. Nevertheless, we find that the neutrality index, a measure of the excess of nonsynonymous polymorphism relative to the neutral expectation, is only weakly significantly different between mitochondrial and nuclear loci. This difference is likely the result of the larger proportion of beneficial mutations in X-linked relative to autosomal loci, and we find little to no difference between mitochondrial and autosomal neutrality indices. Reanalysis of published data from Homo sapiens reveals a similar lack of a difference between the two genomes, although previous studies have suggested a strong difference in both species. Thus, despite a smaller N_e, mitochondrial loci of both flies and humans appear to experience similar efficacies of purifying selection as do loci in the recombining nuclear genome.

Extreme-Depth Re-sequencing of Mitochondrial DNA Finds No Evidence of Paternal Transmission in Humans

Recent reports have questioned the accepted dogma that mammalian mitochondrial DNA (mtDNA) is strictly maternally inherited. In humans, the argument hinges on detecting a signature of inter-molecular recombination in mtDNA sequences sampled at the population level, inferring a paternal source for the mixed haplotypes. However, interpreting these data is fraught with difficulty, and direct experimental evidence is lacking. Using extreme-high depth mtDNA re-sequencing up to ~1.2 million-fold coverage, we find no evidence that paternal mtDNA haplotypes are transmitted to offspring in humans, thus excluding a simple dilution mechanism for uniparental transmission of mtDNA present in all healthy individuals. Our findings indicate that an active mechanism eliminates paternal mtDNA which likely acts at the molecular level.

Emerging evidence raises the possibility that human mitochondrial DNA (mtDNA) is not strictly maternally inherited, but it has not been technically possible to test this hypothesis directly. We identified trios with discordant mtDNA haplotypes, parent-offspring trios were validated using polymorphic microsatellites, and then used extreme-high depth mtDNA re-sequencing to look for paternally transmitted mtDNA. Despite having up to ~1.2 million-fold coverage of mtDNA, we find no evidence that paternal mtDNA haplotypes are transmitted to offspring in humans. Our findings exclude a simple dilution mechanism for uniparental transmission of mtDNA present in all healthy individuals.

Mitochondrial data are not suitable for resolving placental mammal phylogeny

Mitochondrial data have traditionally been used in reconstructing a variety of species phylogenies. The low rates of recombination and thorough characterization of mitochondrial data across vertebrate species make it a particularly attractive phylogenetic marker. The relatively low number of fully sequenced mammal genomes and the lack of extensive sampling within Superorders have posed a serious problem for reaching agreement on the placement mammal species. The use of mitochondrial data sequences from large numbers of mammals could serve to circumvent the taxon-sampling deficit. Here we assess the suitability of mitochondrial data as a phylogenetic marker in mammal phylogenetics. MtDNA datasets of mammal origin have been filtered as follows: (i) we have sampled sparsely across the phylogenetic tree, (ii) we have constrained our sampling to genes with high taxon coverage, (iii) we have categorised rates across sites in a phylogeny independent manner and have removed fast evolving sites, and (iv), we have sampled from very shallow divergence times to reduce phylogenetic conflict. However, topologies obtained using these filters are not consistent with previous studies and are discordant across different genes. Individual mitochondrial genes, and indeed all mitochondrial genes analysed as a supermatrix, resulted in poor resolution of the species phylogeny. Overall, our study highlights the limitations of mitochondrial data, not only for resolving deep divergences and but also for shallow divergences in the mammal phylogeny.

Evidence for weak genetic recombination at the PTP2 locus of Nosema ceranae

The microsporidian Nosema ceranae is an emergent pathogen that threatens the health of honeybees and other pollinators all over the world. Its recent rapid spread across a wide variety of host species and environments demonstrated an enhanced ability of adaptation, which seems to contradict the lack of evidence for genetic recombination and the absence of a sexual stage in its life cycle. Here we retrieved fresh data of the patterns of genetic variation at the PTP2 locus in naturally infected Apis mellifera colonies, by means of single genome amplification. This technique, designed to prevent the formation of chimeric haplotypes during polymerase chain reaction (PCR), provides more reliable estimates of the diversity levels and haplotype structure than standard PCR-cloning methods. Our results are consistent with low but significant rates of recombination in the history of the haplotypes detected: estimates of the population recombination rate are of the order of 30 and support recent evidence for unexpectedly high levels of variation of the parasites within honeybee colonies. These observations suggest the existence of a diploid stage at some point in the life cycle of this parasite and are relevant for our understanding of the dynamics of its expanding population.

Coexistence of minicircular and a highly rearranged mtDNA molecule suggests that recombination shapes mitochondrial genome organization

Recombination has been proposed as a possible mechanism to explain mitochondrial (mt) gene rearrangements, although the issue of whether mtDNA recombination occurs in animals has been controversial. In this study, we sequenced the entire mt genome of the megaspilid wasp Conostigmus sp., which possessed a highly rearranged mt genome. The sequence of the A+T-rich region contained a number of different types of repeats, similar to those reported previously in the nematode Meloidogyne javanica, in which recombination was discovered. In Conostigmus, we detected the end products of recombination: a range of minicircles. However, using isolated (cloned) fragments of the A+T-rich region, we established that some of these minicircles were found to be polymerase chain reaction (PCR) artifacts. It appears that regions with repeats are prone to PCR template switching or PCR jumping. Nevertheless, there is strong evidence that one minicircle is real, as amplification primers that straddle the putative breakpoint junction produce a single strong amplicon from genomic DNA but not from the cloned A+T-rich region. The results provide support for the direct link between recombination and mt gene rearrangement. Furthermore, we developed a model of recombination which is important for our understanding of mtDNA evolution.

Rapid and sustained autologous neutralizing response leading to early spontaneous recovery after HCV infection

After HCV infection, the association between the humoral response and viral sequence evolution remains unclear. We investigated the mechanisms leading to early HCV clearance and spontaneous recovery in two patients. The early evolution of the HCV envelope glycoproteins, and the infectivity spectrum of variants were explored using retroviral pseudoparticles bearing HCV envelopes. Ability of the autologous neutralizing response to control these variants was analyzed. For the first case, the maximum neutralizing activity was for serum collected between two and three months post ALT peak, this activity was still detectable after 30 months. For the second case, autologous neutralizing activity against the variant isolated at the ALT peak was detected in every serum collected between 4 days and 13 months after. The neutralizing response was sustained beyond the time at which the virus was cleared. This raise interesting questions about the role of such antibodies in case of re-exposure.

How good are indirect tests at detecting recombination in human mtDNA?

Empirical proof of human mitochondrial DNA (mtDNA) recombination in somatic tissues was obtained in 2004; however, a lack of irrefutable evidence exists for recombination in human mtDNA at the population level. Our inability to demonstrate convincingly a signal of recombination in population data sets of human mtDNA sequence may be due, in part, to the ineffectiveness of current indirect tests. Previously, we tested some well-established indirect tests of recombination (linkage disequilibrium vs. distance using D' and r(2), Homoplasy Test, Pairwise Homoplasy Index, Neighborhood Similarity Score, and Max χ(2)) on sequence data derived from the only empirically confirmed case of human mtDNA recombination thus far and demonstrated that some methods were unable to detect recombination. Here, we assess the performance of these six well-established tests and explore what characteristics specific to human mtDNA sequence may affect their efficacy by simulating sequence under various parameters with levels of recombination (ρ) that vary around an empirically derived estimate for human mtDNA (population parameter ρ = 5.492). No test performed infallibly under any of our scenarios, and error rates varied across tests, whereas detection rates increased substantially with ρ values > 5.492. Under a model of evolution that incorporates parameters specific to human mtDNA, including rate heterogeneity, population expansion, and ρ = 5.492, successful detection rates are limited to a range of 7-70% across tests with an acceptable level of false-positive results: the neighborhood similarity score incompatibility test performed best overall under these parameters. Population growth seems to have the greatest impact on recombination detection probabilities across all models tested, likely due to its impact on sequence diversity. The implications of our findings on our current understanding of mtDNA recombination in humans are discussed.

Within-host evolution of Staphylococcus aureus during asymptomatic carriage

BACKGROUND

Staphylococcus aureus is a major cause of healthcare associated mortality, but like many important bacterial pathogens, it is a common constituent of the normal human body flora. Around a third of healthy adults are carriers. Recent evidence suggests that evolution of S. aureus during nasal carriage may be associated with progression to invasive disease. However, a more detailed understanding of within-host evolution under natural conditions is required to appreciate the evolutionary and mechanistic reasons why commensal bacteria such as S. aureus cause disease. Therefore we examined in detail the evolutionary dynamics of normal, asymptomatic carriage. Sequencing a total of 131 genomes across 13 singly colonized hosts using the Illumina platform, we investigated diversity, selection, population dynamics and transmission during the short-term evolution of S. aureus.

PRINCIPAL FINDINGS

We characterized the processes by which the raw material for evolution is generated: micro-mutation (point mutation and small insertions/deletions), macro-mutation (large insertions/deletions) and the loss or acquisition of mobile elements (plasmids and bacteriophages). Through an analysis of synonymous, non-synonymous and intergenic mutations we discovered a fitness landscape dominated by purifying selection, with rare examples of adaptive change in genes encoding surface-anchored proteins and an enterotoxin. We found evidence for dramatic, hundred-fold fluctuations in the size of the within-host population over time, which we related to the cycle of colonization and clearance. Using a newly-developed population genetics approach to detect recent transmission among hosts, we revealed evidence for recent transmission between some of our subjects, including a husband and wife both carrying populations of methicillin-resistant S. aureus (MRSA).

SIGNIFICANCE

This investigation begins to paint a picture of the within-host evolution of an important bacterial pathogen during its prevailing natural state, asymptomatic carriage. These results also have wider significance as a benchmark for future systematic studies of evolution during invasive S. aureus disease.

Mitochondria, oxidative DNA damage, and aging

Protection from reactive oxygen species (ROS) and from mitochondrial oxidative damage is well known to be necessary to longevity. The relevance of mitochondrial DNA (mtDNA) to aging is suggested by the fact that the two most commonly measured forms of mtDNA damage, deletions and the oxidatively induced lesion 8-oxo-dG, increase with age. The rate of increase is species-specific and correlates with maximum lifespan. It is less clear that failure or inadequacies in the protection from reactive oxygen species (ROS) and from mitochondrial oxidative damage are sufficient to explain senescence. DNA containing 8-oxo-dG is repaired by mitochondria, and the high ratio of mitochondrial to nuclear levels of 8-oxo-dG previously reported are now suspected to be due to methodological difficulties. Furthermore, MnSOD -/+ mice incur higher than wild type levels of oxidative damage, but do not display an aging phenotype. Together, these findings suggest that oxidative damage to mitochondria is lower than previously thought, and that higher levels can be tolerated without physiological consequence. A great deal of work remains before it will be known whether mitochondrial oxidative damage is a "clock" which controls the rate of aging. The increased level of 8-oxo-dG seen with age in isolated mitochondria needs explanation. It could be that a subset of cells lose the ability to protect or repair mitochondria, resulting in their incurring disproportionate levels of damage. Such an uneven distribution could exceed the reserve capacity of these cells and have serious physiological consequences. Measurements of damage need to focus more on distribution, both within tissues and within cells. In addition, study must be given to the incidence and repair of other DNA lesions, and to the possibility that repair varies from species to species, tissue to tissue, and young to old.

High mitochondrial mutation rates estimated from deep-rooting Costa Rican pedigrees

Estimates of mutation rates for the noncoding hypervariable Region I (HVR-I) of mitochondrial DNA vary widely, depending on whether they are inferred from phylogenies (assuming that molecular evolution is clock-like) or directly from pedigrees. All pedigree-based studies so far were conducted on populations of European origin. In this article, we analyzed 19 deep-rooting pedigrees in a population of mixed origin in Costa Rica. We calculated two estimates of the HVR-I mutation rate, one considering all apparent mutations, and one disregarding changes at sites known to be mutational hot spots and eliminating genealogy branches which might be suspected to include errors, or unrecognized adoptions along the female lines. At the end of this procedure, we still observed a mutation rate equal to 1.24 × 10(-6) , per site per year, i.e., at least threefold as high as estimates derived from phylogenies. Our results confirm that mutation rates observed in pedigrees are much higher than estimated assuming a neutral model of long-term HVRI evolution. We argue that until the cause of these discrepancies will be fully understood, both lower estimates (i.e., those derived from phylogenetic comparisons) and higher, direct estimates such as those obtained in this study, should be considered when modeling evolutionary and demographic processes.

Evolutionary genetics of MHC class II beta genes in the brown hare, Lepus europaeus

The genes of the major histocompatibility complex (MHC) are attractive candidates for investigating the link between adaptive variation and individual fitness. High levels of diversity at the MHC are thought to be the result of parasite-mediated selection and there is growing evidence to support this theory. Most studies, however, target just a single gene within the MHC and infer any evidence of selection to be representative of the entire gene region. Here we present data from three MHC class II beta genes (DPB, DQB, and DRB) for brown hares in two geographic regions and compare them against previous results from a class II alpha-chain gene (DQA). We report moderate levels of diversity and high levels of population differentiation in the DQB and DRB genes (Na = 11, D_est = 0.071 and Na = 15, D_est = 0.409, respectively), but not for the DPB gene (Na = 4, D_est = 0.00). We also detected evidence of positive selection within the peptide binding region of the DQB and DRB genes (95% CI, ω > 1.0) but found no signature of selection for DPB. Mutation and recombination were both found to be important processes shaping the evolution of the class II genes. Our findings suggest that while diversifying selection is a significant contributor to the generally high levels of MHC diversity, it does not act in a uniform manner across the entire MHC class II region. The beta-chain genes that we have characterized provide a valuable set of MHC class II markers for future studies of the evolution of adaptive variation in Leporids.

"Patchy-tachy" leads to false positives for recombination

Indirect tests have detected recombination in mitochondrial DNA (mtDNA) from many animal lineages, including mammals. However, it is possible that features of the molecular evolutionary process without recombination could be incorrectly inferred by indirect tests as being due to recombination. We have identified one such example, which we call "patchy-tachy" (PT), where different partitions of sequences evolve at different rates, that leads to an excess of false positives for recombination inferred by indirect tests. To explore this phenomena, we characterized the false positive rates of six widely used indirect tests for recombination using simulations of general models for mtDNA evolution with PT but without recombination. All tests produced 30-99% false positives for recombination, although the conditions that produced the maximal level of false positives differed between the tests. To evaluate the degree to which conditions that exacerbate false positives are found in published sequence data, we turned to 20 animal mtDNA data sets in which recombination is suggested by indirect tests. Using a model where different regions of the sequences were free to evolve at different rates in different lineages, we demonstrated that PT is prevalent in many data sets in which recombination was previously inferred using indirect tests. Taken together, our results argue that PT without recombination is a viable alternative explanation for detection of widespread recombination in animal mtDNA using indirect tests.

Potential efficacy of mitochondrial genes for animal DNA barcoding: a case study using eutherian mammals

BACKGROUND

A well-informed choice of genetic locus is central to the efficacy of DNA barcoding. Current DNA barcoding in animals involves the use of the 5' half of the mitochondrial cytochrome oxidase 1 gene (CO1) to diagnose and delimit species. However, there is no compelling a priori reason for the exclusive focus on this region, and it has been shown that it performs poorly for certain animal groups. To explore alternative mitochondrial barcoding regions, we compared the efficacy of the universal CO1 barcoding region with the other mitochondrial protein-coding genes in eutherian mammals. Four criteria were used for this comparison: the number of recovered species, sequence variability within and between species, resolution to taxonomic levels above that of species, and the degree of mutational saturation.

RESULTS

Based on 1,179 mitochondrial genomes of eutherians, we found that the universal CO1 barcoding region is a good representative of mitochondrial genes as a whole because the high species-recovery rate (> 90%) was similar to that of other mitochondrial genes, and there were no significant differences in intra- or interspecific variability among genes. However, an overlap between intra- and interspecific variability was still problematic for all mitochondrial genes. Our results also demonstrated that any choice of mitochondrial gene for DNA barcoding failed to offer significant resolution at higher taxonomic levels.

CONCLUSIONS

We suggest that the CO1 barcoding region, the universal DNA barcode, is preferred among the mitochondrial protein-coding genes as a molecular diagnostic at least for eutherian species identification. Nevertheless, DNA barcoding with this marker may still be problematic for certain eutherian taxa and our approach can be used to test potential barcoding loci for such groups.

Molecular phylogenetic study of Culex quinquefasciatus mosquito from different geographical regions of India using 16S rRNA gene sequences

Culex quinquefasciatus is a major vector of filariasis and various encephalitis in India and worldwide. Vector control remains the most successful strategy for the suppression of mosquito borne diseases. The genetic structure of vector populations in terms of insecticide resistance and susceptibility or refractoriness to infection may possibly vary. To exploit the genetic variability in vector population could pave the path for the alternative strategies in vector management. The sequences of ribosomal RNA molecules have been widely used for such studies. Here, we examined the molecular phylogenetic relationship among the Cx. quinquefasciatus collected from different geographical regions of India, using 16S ribosomal RNA (16S rRNA) gene nucleotide sequences. The distances among the species were measured using Pearson correlation; the Neighbor-Joining (NJ) method was used for the clustering with appropriate bootstrap values using Data Analysis in Molecular Biology and Evolution (DAMBE) software. The results revealed that the populations are genetically diverse. Based on the distance values and the tree topology on the basis of 16S rRNA sequences reflected the clear biogeographical and geoclimatic pattern among the different geographical populations from India.

Reducing chimera formation during PCR amplification to ensure accurate genotyping

Measurements of population diversity are fundamental to the reconstruction of the evolutionary and epidemiological history of organisms. Commonly used protocols to measure population diversity using the polymerase chain reaction (PCR) are prone to the introduction of artificial chimeras. These are often difficult to detect and can confound the correct interpretation of results due to the false generation of recombinants when the underlying DNA sample contains multiple distinct templates. This study presents a standardised procedure to suppress the formation of artificial chimeras during PCR amplification. The solution is based on the accurate determination of the efficiency and end point of the log-linear phase of a PCR. This procedure will facilitate the generation of data sets that more accurately reflect the underlying population diversity rather than artifacts introduced by the process itself.

Mitochondrial recombination increases with age in Podospora anserina

With uniparental inheritance of mitochondria, there seems little reason for homologous recombination in mitochondria, but the machinery for mitochondrial recombination is quite well-conserved in many eukaryote species. In fungi and yeasts heteroplasmons may be formed when strains fuse and transfer of organelles takes place, making it possible to study mitochondrial recombination when introduced mitochondria contain different markers. A survey of wild-type isolates from a local population of the filamentous fungus Podospora anserina for the presence of seven optional mitochondrial introns indicated that mitochondrial recombination does take place in nature. Moreover the recombination frequency appeared to be correlated with age: the more rapidly ageing fraction of the population had a significantly lower linkage disequilibrium indicating more recombination. Direct confrontation experiments with heterokaryon incompatible strains with different mitochondrial markers at different (relative) age confirmed that mitochondrial recombination increases with age. We propose that with increasing mitochondrial damage over time, mitochondrial recombination - even within a homoplasmic population of mitochondria - is a mechanism that may restore mitochondrial function.

Biparental inheritance of plastidial and mitochondrial DNA and hybrid variegation in Pelargonium

Plastidial (pt) and mitochondrial (mt) genes usually show maternal inheritance. Non-Mendelian, biparental inheritance of plastids was first described by Baur (Z Indukt Abstamm Vererbungslehre 1:330–351, 1909) for crosses between Pelargonium cultivars. We have analyzed the inheritance of pt and mtDNA by examining the progeny from reciprocal crosses of Pelargoniumzonale and P. inquinans using nucleotide sequence polymorphisms of selected pt and mt genes. Sequence analysis of the progeny revealed biparental inheritance of both pt and mtDNA. Hybrid plants exhibited variegation: our data demonstrate that the inquinans chloroplasts, but not the zonale chloroplasts bleach out, presumably due to incompatibility of the former with the hybrid nuclear genome. Different distribution of maternal and paternal sequences could be observed in different sectors of the same leaf, in different leaves of the same plant, and in different plants indicating random segregation and sorting-out of maternal and paternal plastids and mitochondria in the hybrids. The substantial transmission of both maternal and paternal mitochondria to the progeny turns Pelargonium into a particular interesting subject for studies on the inheritance, segregation and recombination of mt genes.

Mitochondrial functional complementation in mitochondrial DNA-based diseases

Mitochondria exist in networks that are continuously remodeled through fusion and fission. Why do individual mitochondria in living cells fuse and divide continuously? Protein machinery and molecular mechanism for the dynamic nature of mitochondria have been almost clarified. However, the biological significance of the mitochondrial fusion and fission events has been poorly understood, although there is a possibility that mitochondrial fusion and fission are concerned with quality controls of mitochondria. trans-mitochondrial cell and mouse models possessing heteroplasmic populations of mitochondrial DNA (mtDNA) haplotypes are quite efficient for answering this question, and one of the answers is "mitochondrial functional complementation" that is able to regulate respiratory function of individual mitochondria according to "one for all, all for one" principle. In this review, we summarize the observations about mitochondrial functional complementation in mammals and discuss its biological significance in pathogeneses of mtDNA-based diseases.

Using classical population genetics tools with heterochroneous data: time matters!

BACKGROUND

New polymorphism datasets from heterochroneous data have arisen thanks to recent advances in experimental and microbial molecular evolution, and the sequencing of ancient DNA (aDNA). However, classical tools for population genetics analyses do not take into account heterochrony between subsets, despite potential bias on neutrality and population structure tests. Here, we characterize the extent of such possible biases using serial coalescent simulations.

METHODOLOGY/PRINCIPAL FINDINGS

We first use a coalescent framework to generate datasets assuming no or different levels of heterochrony and contrast most classical population genetic statistics. We show that even weak levels of heterochrony ( approximately 10% of the average depth of a standard population tree) affect the distribution of polymorphism substantially, leading to overestimate the level of polymorphism theta, to star like trees, with an excess of rare mutations and a deficit of linkage disequilibrium, which are the hallmark of e.g. population expansion (possibly after a drastic bottleneck). Substantial departures of the tests are detected in the opposite direction for more heterochroneous and equilibrated datasets, with balanced trees mimicking in particular population contraction, balancing selection, and population differentiation. We therefore introduce simple corrections to classical estimators of polymorphism and of the genetic distance between populations, in order to remove heterochrony-driven bias. Finally, we show that these effects do occur on real aDNA datasets, taking advantage of the currently available sequence data for Cave Bears (Ursus spelaeus), for which large mtDNA haplotypes have been reported over a substantial time period (22-130 thousand years ago (KYA)).

CONCLUSIONS/SIGNIFICANCE

Considering serial sampling changed the conclusion of several tests, indicating that neglecting heterochrony could provide significant support for false past history of populations and inappropriate conservation decisions. We therefore argue for systematically considering heterochroneous models when analyzing heterochroneous samples covering a large time scale.

Can indirect tests detect a known recombination event in human mtDNA?

Whether human mitochondrial DNA (mtDNA) recombines sufficiently to influence its evolution, evolutionary analysis, and disease etiology, remains equivocal. Overall, evidence from indirect studies of population genetic data suggests that recombination is not occurring at detectable levels. This may be explained by no, or low, recombination or, alternatively, current indirect tests may be incapable of detecting recombination in human mtDNA. To investigate the latter, we have tested whether six well-established indirect tests of recombination could detect recombination in a human mtDNA data set, in which its occurrence had been empirically confirmed. Three showed statistical evidence for recombination (r(2) vs. distance, the Homoplasy test, Neighborhood Similarity Score), and three did not (D' vs. distance, Max Chi Squared, Pairwise Homoplasy Index). Possible reasons for detection failure are discussed. Further, evidence from earlier studies suggesting a lack of recombination in mtDNA in humans is reconsidered, taking into account the appropriateness of the tests used, based on our new findings.

Amino acid covariation in a functionally important human immunodeficiency virus type 1 protein region is associated with population subdivision

The frequently reported amino acid covariation of the highly polymorphic human immunodeficiency virus type 1 (HIV-1) exterior envelope glycoprotein V3 region has been assumed to reflect fitness epistasis between residues. However, nonrandom association of amino acids, or linkage disequilibrium, has many possible causes, including population subdivision. If the amino acids at a set of sequence sites differ in frequencies between subpopulations, then analysis of the whole population may reveal linkage disequilibrium even if it does not exist in any subpopulation. HIV-1 has a complex population structure, and the effects of this structure on linkage disequilibrium were investigated by estimating within- and among-subpopulation components of variance in linkage disequilibrium. The amino acid covariation previously reported is explained by differences in amino acid frequencies among virus subpopulations in different patients and by nonsystematic disequilibrium among patients. Disequilibrium within patients appears to be entirely due to differences in amino acid frequencies among sampling time points and among chemokine coreceptor usage phenotypes of virus particles, but not source tissues. Positive selection explains differences in allele frequencies among time points and phenotypes, indicating that these differences are adaptive rather than due to genetic drift. However, the absence of a correlation between linkage disequilibrium and phenotype suggests that fitness epistasis is an unlikely cause of disequilibrium. Indeed, when population structure is removed by analyzing sequences from a single time point and phenotype, no disequilibrium is detectable within patients. These results caution against interpreting amino acid covariation and coevolution as evidence for fitness epistasis.

Doubly uniparental inheritance: two mitochondrial genomes, one precious model for organelle DNA inheritance and evolution

Eukaryotes have exploited several mechanisms for organelle uniparental inheritance, so this feature arose and evolved independently many times in their history. Metazoans' mitochondria commonly experience strict maternal inheritance; that is, they are only transmitted by females. However, the most noteworthy exception comes from some bivalve mollusks, in which two mitochondrial lineages (together with their genomes) are inherited: one through females (F) and the other through males (M). M and F genomes show up to 30% sequence divergence. This inheritance mechanism is known as doubly uniparental inheritance (DUI), because both sexes inherit uniparentally their mitochondria. Here, we review what we know about this unusual system, and we propose a model for evolution of DUI that might account for its origin as sex determination mechanism. Moreover, we propose DUI as a choice model to address many aspects that should be of interest to a wide range of biological subfields, such as mitochondrial inheritance, mtDNA evolution and recombination, genomic conflicts, evolution of sex, and developmental biology. Actually, as research proceeds, mitochondria appear to have acquired a central role in many fundamental processes of life, which are not only in their metabolic activity as cellular power plants, such as cell signaling, fertilization, development, differentiation, ageing, apoptosis, and sex determination. A function of mitochondria in the origin and maintenance of sex has been also proposed.

Mitochondrial complementation preventing respiratory dysfunction caused by mutant mtDNA

The mitochondrial theory of aging is the idea that age-associated mitochondrial dysfunction is caused by accumulation of somatic mutations in mitochondrial DNA (mtDNA). However, mitochondria are considered to be a dynamic organelle that repeats fusion and fission. Through fusion and fission, there is an extensive and continuous exchange of mtDNA and its products between mitochondria. This mitochondrial complementation prevents individuals from expression of respiratory dysfunction caused by pathogenic mutant mtDNAs. Thus, the presence of mitochondrial complementation does not support the mitochondrial theory of aging. Moreover, the presence of mitochondrial complementation enables gene therapy for mitochondrial diseases using nuclear transplantation of zygotes. (c) 2009 International Union of Biochemistry and Molecular Biology, Inc.

Revealing the hidden complexities of mtDNA inheritance

Mitochondrial DNA (mtDNA) is a pivotal tool in molecular ecology, evolutionary and population genetics. The power of mtDNA analyses derives from a relatively high mutation rate and the apparent simplicity of mitochondrial inheritance (maternal, without recombination), which has simplified modelling population history compared to the analysis of nuclear DNA. However, in biology things are seldom simple, and advances in DNA sequencing and polymorphism detection technology have documented a growing list of exceptions to the central tenets of mitochondrial inheritance, with paternal leakage, heteroplasmy and recombination now all documented in multiple systems. The presence of paternal leakage, recombination and heteroplasmy can have substantial impact on analyses based on mtDNA, affecting phylogenetic and population genetic analyses, estimates of the coalescent and the myriad of other parameters that are dependent on such estimates. Here, we review our understanding of mtDNA inheritance, discuss how recent findings mean that established ideas may need to be re-evaluated, and we assess the implications of these new-found complications for molecular ecologists who have relied for decades on the assumption of a simpler mode of inheritance. We show how it is possible to account for recombination and heteroplasmy in evolutionary and population analyses, but that accurate estimates of the frequencies of biparental inheritance and recombination are needed. We also suggest how nonclonal inheritance of mtDNA could be exploited, to increase the ways in which mtDNA can be used in analyses.