Mitochondrial genome recombination in the zone of contact between two hybridizing conifers

A Systematic Review and Developmental Perspective on Origin of CMS Genes in Crops

Cytoplasmic male sterility (CMS) arises from the incompatibility between the nucleus and cytoplasm as typical representatives of the chimeric structures in the mitochondrial genome (mitogenome), which has been extensively applied for hybrid seed production in various crops. The frequent occurrence of chimeric mitochondrial genes leading to CMS is consistent with the mitochondrial DNA (mtDNA) evolution. The sequence conservation resulting from faithfully maternal inheritance and the chimeric structure caused by frequent sequence recombination have been defined as two major features of the mitogenome. However, when and how these chimeric mitochondrial genes appear in the context of the highly conserved reproduction of mitochondria is an enigma. This review, therefore, presents the critical view of the research on CMS in plants to elucidate the mechanisms of this phenomenon. Generally, distant hybridization is the main mechanism to generate an original CMS source in natural populations and in breeding. Mitochondria and mitogenomes show pleomorphic and dynamic changes at key stages of the life cycle. The promitochondria in dry seeds develop into fully functioning mitochondria during seed imbibition, followed by massive mitochondria or mitogenome fusion and fission in the germination stage along with changes in the mtDNA structure and quantity. The mitogenome stability is controlled by nuclear loci, such as the nuclear gene Msh1. Its suppression leads to the rearrangement of mtDNA and the production of heritable CMS genes. An abundant recombination of mtDNA is also often found in distant hybrids and somatic/cybrid hybrids. Since mtDNA recombination is ubiquitous in distant hybridization, we put forward a hypothesis that the original CMS genes originated from mtDNA recombination during the germination of the hybrid seeds produced from distant hybridizations to solve the nucleo-cytoplasmic incompatibility resulting from the allogenic nuclear genome during seed germination.

Phylogenetic Relationships and Next-Generation Barcodes in the Genus Torreya Reveal a High Proportion of Misidentified Cultivated Plants

Accurate species identification is key to conservation and phylogenetic inference. Living plant collections from botanical gardens/arboretum are important resources for the purpose of scientific research, but the proportion of cultivated plant misidentification are un-tested using DNA barcodes. Here, we assembled the next-generation barcode (complete plastid genome and complete nrDNA cistron) and mitochondrial genes from genome skimming data of Torreya species with multiple accessions for each species to test the species discrimination and the misidentification proportion of cultivated plants used in Torreya studies. A total of 38 accessions were included for analyses, representing all nine recognized species of genus Torreya. The plastid phylogeny showed that all 21 wild samples formed species-specific clades, except T. jiulongshanensis. Disregarding this putative hybrid, seven recognized species sampled here were successfully discriminated by the plastid genome. Only the T. nucifera accessions grouped into two grades. The species identification rate of the nrDNA cistron was 62.5%. The Skmer analysis based on nuclear reads from genome skims showed promise for species identification with seven species discriminated. The proportion of misidentified cultivated plants from arboreta/botanical gardens was relatively high with four accessions (23.5%) representing three species. Interspecific relationships within Torreya were fully resolved with maximum support by plastomes, where Torreya jackii was on the earliest diverging branch, though sister to T. grandis in the nrDNA cistron tree, suggesting that this is likely a hybrid species between T. grandis and an extinct Torreya ancestor lineage. The findings here provide quantitative insights into the usage of cultivated samples for phylogenetic study.

Control of plastid inheritance by environmental and genetic factors

The genomes of cytoplasmic organelles (mitochondria and plastids) are maternally inherited in most eukaryotes, thus excluding organellar genomes from the benefits of sexual reproduction and recombination. The mechanisms underlying maternal inheritance are largely unknown. Here we demonstrate that two independently acting mechanisms ensure maternal inheritance of the plastid (chloroplast) genome. Conducting large-scale genetic screens for paternal plastid transmission, we discovered that mild chilling stress during male gametogenesis leads to increased entry of paternal plastids into sperm cells and strongly increased paternal plastid transmission. We further show that the inheritance of paternal plastid genomes is controlled by the activity of a genome-degrading exonuclease during pollen maturation. Our data reveal that (1) maternal inheritance breaks down under specific environmental conditions, (2) an organelle exclusion mechanism and a genome degradation mechanism act in concert to prevent paternal transmission of plastid genes and (3) plastid inheritance is determined by complex gene-environment interactions.

Recombinant Mitochondrial Genomes Reveal Recent Interspecific Hybridization between Invasive Salangid Fishes

The interspecific recombination of the mitochondrial (mt) genome, if not an experimental artifact, may result from interbreeding of species with broken reproductive barriers, which, in turn, is a frequent consequence of human activities including species translocations, habitat modifications, and climate change. This issue, however, has not been addressed for Protosalanx chinensis and other commercially important and, simultaneously, invasive salangid fishes that were the product of successful aquaculture in China. To assess the probability of interspecific hybridization, we analyzed the patterns of diversity and recombination in the complete mitochondrial (mt) genomes of these fishes using the GenBank resources. A sliding window analysis revealed a non-uniform distribution of the intraspecific differences in P. chinensis with four highly pronounced peaks of divergence centered at the COI, ND4L-ND4, and ND5 genes, and also at the control region. The corresponding divergent regions in P. chinensis show a high sequence similarity (99−100%) to the related salangid fishes, Neosalanx tangkahkeii and N. anderssoni. This observation suggests that the divergent regions of P. chinensis may represent a recombinant mitochondrial DNA (mtDNA) containing mt genome fragments belonging to different salangid species. Indeed, four, highly significant (pairwise homoplasy index test, P < 0.00001) signals of recombination have been revealed at coordinates closely corresponding to the divergent regions. The recombinant fragments are, however, not fixed, and different mt genomes of P. chinensis are mosaic, containing different numbers of recombinant events. These facts, along with the high similarity or full identity of the recombinant fragments between the donor and the recipient sequences, indicate a recent interspecific hybridization between P. chinensis and two Neosalanx species. Alternative hypotheses, including taxonomical misidentifications, sequence misalignments, DNA contamination, and/or artificial PCR recombinants, are not supported by the data. The recombinant fragments revealed in our study represent diagnostic genetic markers for the identification and distinguishing of hybrids, which can be used to control the invasive dynamics of hybrid salangid fishes.

Sequence Capture of Mitochondrial Genome with PCR-Generated Baits Provides New Insights into the Biogeography of the Genus Abies Mill

Mitochondrial DNA (mtDNA), being maternally inherited in plants of the family Pinaceae, is an important source of phylogeographic information. However, its use is hindered by a low mutation rate and frequent structure rearrangements. In the present study, we tested the method of genomic libraries enrichment with mtDNA via the sequence capture method yielding mtDNA data which were further used to reconstruct the phylogenetic tree of the genus Abies. The baits for hybrid capture were obtained by long-range PCR using primers designed on the basis of the assembly of Abies sibirica Ledeb. mitochondrial genome. Mitochondrial genomes of Picea sitchensis (Bong.) Carr., Larix sibirica Ledeb., and Keteleeria davidiana (Bertrand) Beissn. were used as an outgroup. The resulting phylogenetic tree consists of two sister branches, including the Eurasian and American species, respectively, with some exceptions. The subclade of A. sachalinensis (F. Schmidt) Mast. and A. veitchii Lindl. (Japan and Sakhalin islands) occupies a basal position in the branch of American firs, probably due to the complex history of fir migrations from North America to Eurasia. The tree has high support for majority of clades. For species represented by more than one sample an intraspecific variability was found which is suitable to design mtDNA markers for phylogeographic and population studies.

Complete Mitochondrial Genome of a Gymnosperm, Sitka Spruce (Picea sitchensis), Indicates a Complex Physical Structure

Plant mitochondrial genomes vary widely in size. Although many plant mitochondrial genomes have been sequenced and assembled, the vast majority are of angiosperms, and few are of gymnosperms. Most plant mitochondrial genomes are smaller than a megabase, with a few notable exceptions. We have sequenced and assembled the complete 5.5-Mb mitochondrial genome of Sitka spruce (Picea sitchensis), to date, one of the largest mitochondrial genomes of a gymnosperm. We sequenced the whole genome using Oxford Nanopore MinION, and then identified contigs of mitochondrial origin assembled from these long reads based on sequence homology to the white spruce mitochondrial genome. The assembly graph shows a multipartite genome structure, composed of one smaller 168-kb circular segment of DNA, and a larger 5.4-Mb single component with a branching structure. The assembly graph gives insight into a putative complex physical genome structure, and its branching points may represent active sites of recombination.

Heteroplasmy and Patterns of Cytonuclear Linkage Disequilibrium in Wild Carrot

Organellar genomes are considered to be strictly uniparentally-inherited. Uniparental inheritance allows for cytonuclear coevolution and the development of highly coordinated cytonuclear interactions. Yet, instances of biparental inheritance have been documented across eukaryotes. Biparental inheritance in otherwise uniparentally-inherited organelles is termed leakage (maternal or paternal) and allows for the presence of multiple variants of the same organellar genome within an individual, called heteroplasmy. It is unclear what, if any, evolutionary consequences are placed on nuclear and/or organellar genomes due to heteroplasmy. One way of accessing cytonuclear interactions and potential coevolution is through calculating cytonuclear linkage disequilibrium (cnLD), or the non-random association of alleles between nuclear and organellar genomes. Patterns of cnLD can indicate positive or negative cytonuclear selection, coevolution between the nuclear and organellar genomes, non-traditional organellar inheritance, or instances of ancestral heteroplasmy. In plants, cytonuclear interactions have been shown to play a role in cytoplasmic male sterility which occurs in gynodioecious species and is associated with leakage. We used the gynodioecious species, Daucus carota L. spp. carota, or wild carrot, to investigate cnLD. We genotyped a total of 265 individuals from two regions of the USA at 15 nuclear microsatellites, the mitochondrial genes cox1 and atp9, and an intergenic region between trnS and trnG (StoG) in the plastid genome to calculate nuclear–nuclear LD (nucLD), cnLD, and organellar LD (i.e., within the mtDNA and between mtDNA and ptDNA) within the two regions. We were further able to identify cox1 and StoG heteroplasmy and calculate some of the same LD measures within heteroplasmic and homoplasmic (non-heteroplasmic) datasets. We used a Z-transformation test to demonstrate that heteroplasmic individuals display significantly higher levels of cnLD within both regions. In spite of this, within and between organellar LD is low to moderate. Given these patterns of LD in two regions of the USA in which gene flow has been shown to occur between crop and wild carrot, we suggest that heteroplasmy is an evolutionary mechanism which permits the maintenance of cnLD while also acting to disrupt organellar LD.

Patterns of mtDNA variation reveal complex evolutionary history of relict and endangered peat bog pine (Pinus uliginosa)

Estimates of genetic differentiation at intra- and interspecific level are often hindered by the lack of suitable molecular markers. Low phylogeographic resolution limits development of appropriate conservation strategies especially in case of endangered forest tree species with small and disjunct distribution. In this study, we assessed fine-scale genetic structure of relict and endangered peat bog pine (Pinus uliginosa) and two other closely related European pine species (Pinus mugo and Pinus uncinata) using a set of 15 newly developed maternally inherited and seed-mediated mitochondrial DNA (mtDNA) markers and two previously known polymorphic mtDNA regions (nad1, nad7). Three main groups, corresponding in general to three investigated species were revealed in the haplotype network analysis. However, only P. uncinata was clearly distinct at all levels of analysis, whereas great genetic similarity and haplotype sharing was observed between P. uliginosa and P. mugo. Strong phylogeographic structure was found in P. uliginosa that showed high differentiation at relatively short geographical distance among populations and the existence of mitochondrial lineages of different evolutionary history. Hybridization with other pine species has likely contributed to genetic differentiation of P. uliginosa as indicated by contemporary distribution of mtDNA haplotypes. The research emphasizes the importance of accurate assessments of genetic structure of endangered species with complex evolutionary history for development of efficient conservation strategies.

Mitochondrial introgression and interspecies recombination in the Fusarium fujikuroi species complex

The Fusarium fujikuroi species complex (FFSC) is an economically important monophyletic lineage in the genus Fusarium. Incongruence observed among mitochondrial gene trees, as well as the multiple non-orthologous copies of the internal transcribed spacer region of the ribosomal RNA genes, suggests that the origin and history of this complex likely involved interspecies gene flow. Based on this hypothesis, the mitochondrial genomes of non-conspecific species should harbour signatures of introgression or introgressive hybridization. The aim of this study was therefore to search for recombination between the mitochondrial genomes of different species in the FFSC. Using methods based on mt genome sequence similarity, five significant recombinant regions in both gene and intergenic regions were detected. Using coalescent-based methods and the sequences for individual mt genes, various ancestral recombination events between different lineages of the FFSC were also detected. These findings suggest that interspecies gene flow and introgression are likely to have played key roles in the evolution of the FFSC at both ancient and more recent time scales.

Mitochondrial Recombination and Introgression during Speciation by Hybridization

Genome recombination is a major source of genotypic diversity and contributes to adaptation and speciation following interspecies hybridization. The contribution of recombination in these processes has been thought to be largely limited to the nuclear genome because organelles are mostly uniparentally inherited in animals and plants, which prevents recombination. Unicellular eukaryotes such as budding yeasts do, however, transmit mitochondria biparentally, suggesting that during hybridization, both parents could provide alleles that contribute to mitochondrial functions such as respiration and metabolism in hybrid populations or hybrid species. We examined the dynamics of mitochondrial genome transmission and evolution during speciation by hybridization in the natural budding yeast Saccharomyces paradoxus. Using population-scale mitochondrial genome sequencing in two endemic North American incipient species SpB and SpC and their hybrid species SpC*, we found that both parental species contributed to the hybrid mitochondrial genome through recombination. We support our findings by showing that mitochondrial recombination between parental types is frequent in experimental crosses that recreate the early step of this speciation event. In these artificial hybrids, we observed that mitochondrial genome recombination enhances phenotypic variation among diploid hybrids, suggesting that it could play a role in the phenotypic differentiation of hybrid species. Like the nuclear genome, the mitochondrial genome can, therefore, also play a role in hybrid speciation.

Interspecific Plastome Recombination Reflects Ancient Reticulate Evolution in Picea (Pinaceae)

Plastid sequences are a cornerstone in plant systematic studies and key aspects of their evolution, such as uniparental inheritance and absent recombination, are often treated as axioms. While exceptions to these assumptions can profoundly influence evolutionary inference, detecting them can require extensive sampling, abundant sequence data, and detailed testing. Using advancements in high-throughput sequencing, we analyzed the whole plastomes of 65 accessions of Picea, a genus of ∼35 coniferous forest tree species, to test for deviations from canonical plastome evolution. Using complementary hypothesis and data-driven tests, we found evidence for chimeric plastomes generated by interspecific hybridization and recombination in the clade comprising Norway spruce (P. abies) and 10 other species. Support for interspecific recombination remained after controlling for sequence saturation, positive selection, and potential alignment artifacts. These results reconcile previous conflicting plastid-based phylogenies and strengthen the mounting evidence of reticulate evolution in Picea. Given the relatively high frequency of hybridization and biparental plastid inheritance in plants, we suggest interspecific plastome recombination may be more widespread than currently appreciated and could underlie reported cases of discordant plastid phylogenies.

Complete Mitochondrial Genomes of the Cherskii's Sculpin Cottus czerskii and Siberian Taimen Hucho taimen Reveal GenBank Entry Errors: Incorrect Species Identification and Recombinant Mitochondrial Genome

The complete mitochondrial (mt) genome is sequenced in 2 individuals of the Cherskii's sculpin Cottus czerskii. A surprisingly high level of sequence divergence (10.3%) has been detected between the 2 genomes of C czerskii studied here and the GenBank mt genome of C czerskii (KJ956027). At the same time, a surprisingly low level of divergence (1.4%) has been detected between the GenBank C czerskii (KJ956027) and the Amur sculpin Cottus szanaga (KX762049, KX762050). We argue that the observed discrepancies are due to incorrect taxonomic identification so that the GenBank accession number KJ956027 represents actually the mt genome of C szanaga erroneously identified as C czerskii. Our results are of consequence concerning the GenBank database quality, highlighting the potential negative consequences of entry errors, which once they are introduced tend to be propagated among databases and subsequent publications. We illustrate the premise with the data on recombinant mt genome of the Siberian taimen Hucho taimen (NCBI Reference Sequence Database NC_016426.1; GenBank accession number HQ897271.1), bearing 2 introgressed fragments (≈0.9 kb [kilobase]) from 2 lenok subspecies, Brachymystax lenok and Brachymystax lenok tsinlingensis, submitted to GenBank on June 12, 2011. Since the time of submission, the H taimen recombinant mt genome leading to incorrect phylogenetic inferences was propagated in multiple subsequent publications despite the fact that nonrecombinant H taimen genomes were also available (submitted to GenBank on August 2, 2014; KJ711549, KJ711550). Other examples of recombinant sequences persisting in GenBank are also considered. A GenBank Entry Error Depositary is urgently needed to monitor and avoid a progressive accumulation of wrong biological information.

Resolving deep relationships of PACMAD grasses: a phylogenomic approach

BACKGROUND

Plastome sequences for 18 species of the PACMAD grasses (subfamilies Panicoideae, Aristidoideae, Chloridoideae, Micrairoideae, Arundinoideae, Danthonioideae) were analyzed phylogenomically. Next generation sequencing methods were used to provide complete plastome sequences for 12 species. Sanger sequencing was performed to determine the plastome of one species, Hakonechloa macra, to provide a reference for annotation. These analyses were conducted to resolve deep subfamilial relationships within the clade. Divergence estimates were assessed to determine potential factors that led to the rapid radiation of this lineage and its dominance of warmer open habitats.

RESULTS

New plastomes were completely sequenced and characterized for 13 PACMAD species. An autapomorphic ~1140 bp deletion was found in Hakonechloa macra putatively pseudogenizing rpl14 and eliminating rpl16 from this plastome. Phylogenomic analyses support Panicoideae as the sister group to the ACMAD clade. Complete plastome sequences provide greater support at deep nodes within the PACMAD clade. The initial diversification of PACMAD subfamilies was estimated to occur at 32.4 mya.

CONCLUSIONS

Phylogenomic analyses of complete plastomes provides resolution for deep relationships of PACMAD grasses. The divergence estimate of 32.4 mya at the crown node of the PACMAD clade coincides with the Eocene-Oligocene Transition (EOT). The Eocene was a period of global cooling and drying, which led to forest fragmentation and the expansion of open habitats now dominated by these grasses. Understanding how these grasses are related and determining a cause for their rapid radiation allows for future predictions of grassland distribution in the face of a changing global climate.

On the adaptive value of cytoplasmic genomes in plants

Is DNA variation maintained in organelle genomes selectively neutral? The answer to this question has important implications for many aspects of ecology and evolution. While traditionally the answer has been 'yes', recent studies in animals have shown that, on the contrary, mitochondrial DNA polymorphism is frequently adaptive. In plants, however, the neutrality assumption has not been strongly challenged. Here, we begin with a critical evaluation of arguments in favour of this long-held view. We then discuss the latest empirical evidence for the opposing prediction that sequence variation in plant cytoplasmic genomes is frequently adaptive. While outstanding research progress is being made towards understanding this fundamental topic, we highlight the need for studies that combine information ranging from field experiments to physiology to molecular evolutionary biology. Such an interdisciplinary approach provides a means for determining the frequency, drivers and evolutionary significance of adaptive organelle DNA variation.

Mitochondrial DNA capture and divergence in Pinus provide new insights into the evolution of the genus

The evolution of the mitochondrial (mt) genome is far from being fully understood. Systematic investigations into the modes of inheritance, rates and patterns of recombination, nucleotide substitution, and structural changes in the mt genome are still lacking in many groups of plants. In this study, we sequenced >11kbp mtDNA segments from multiple accessions of 36 pine species to characterize the evolutionary patterns of mtDNA in the genus Pinus. We found extremely low substitution rates and complex repetitive sequences scattered across different genome regions, as well as chimeric structures that were probably generated by multiple intergenomic recombinations. The mtDNA-based phylogeny of the genus differed from that based on chloroplast and nuclear DNA in the placement of several groups of species. Such discordances suggest a series of mtDNA capture events during past range shifts of the pine species and that both vertical and horizontal inheritance are implicated in the evolution of mtDNA in Pinus. MtDNA dating revealed that most extant lineages of the genus originated during Oligocene-Miocene radiation and subgenus Strobus diversified earlier than subgenus Pinus. Our findings illustrate a reticular evolutionary pathway for the mt genome through capture and recombination in the genus Pinus, and provide new insights into the evolution of the genus.

Population genetic evidence for speciation pattern and gene flow between Picea wilsonii, P. morrisonicola and P. neoveitchii

BACKGROUND AND AIMS

Genetic drift due to geographical isolation, gene flow and mutation rates together make it difficult to determine the evolutionary relationships of present-day species. In this study, population genetic data were used to model and decipher interspecific relationships, speciation patterns and gene flow between three species of spruce with similar morphology, Picea wilsonii, P. neoveitchii and P. morrisonicola. Picea wilsonii and P. neoveitchii occur from central to north-west China, where they have overlapping distributions. Picea morrisonicola, however, is restricted solely to the island of Taiwan and is isolated from the other two species by a long distance.

METHODS

Sequence variations were examined in 18 DNA fragments for 22 populations, including three fragments from the chloroplast (cp) genome, two from the mitochondrial (mt) genome and 13 from the nuclear genome.

KEY RESULTS

In both the cpDNA and the mtDNA, P. morrisonicola accumulated more species-specific mutations than the other two species. However, most nuclear haplotypes of P. morrisonicola were shared by P. wilsonii, or derived from the dominant haplotypes found in that species. Modelling of population genetic data supported the hypothesis that P. morrisonicola derived from P. wilsonii within the more recent past, most probably indicating progenitor-derivative speciation with a distinct bottleneck, although further gene flow from the progenitor to the derivative continued. In addition, the occurrence was detected of an obvious mtDNA introgression from P. neoveitchii to P. wilsonii despite their early divergence.

CONCLUSIONS

The extent of mutation, introgression and lineage sorting taking place during interspecific divergence and demographic changes in the three species had varied greatly between the three genomes. The findings highlight the complex evolutionary histories of these three Asian spruce species.

Paternal leakage, heteroplasmy, and the evolution of plant mitochondrial genomes

Plant mitochondrial genomes are usually transmitted to the progeny from the maternal parent. However, cases of paternal transmission are known and are perhaps more common than once thought. This review will consider recent evidence, both direct and indirect, of paternal transmission (leakage) of the mitochondrial genome of seed plants, especially in natural populations, and how this can result in offspring that carry a mixture of maternally and paternally derived copies of the genome; a type of heteroplasmy. It will further consider how this heteroplasmy facilitates recombination between genetically distinct partners; a process that can enhance mitochondrial genotypic diversity. This will then form the basis for a discussion of five evolutionary questions that arise from these observations. Questions include how plant mitochondrial genome evolution can be placed on a sexual to asexual continuum, whether cytoplasmic male sterility (CMS) facilitates the evolution of paternal leakage, whether paternal leakage is more likely in populations undergoing admixture, how leakage influences patterns of gene flow, and whether heteroplasmy occurs in natural populations at a frequency greater than predicted by crossing experiments. It is proposed that each of these questions offers fertile ground for future research on a diversity of plant species.

Mitochondrial DNA variation and introgression in Siberian taimen Hucho taimen

Siberian taimen Hucho taimen is the largest representative of the family Salmonidae inhabiting rivers of northern Eurasia. The species is under intensive aquaculture activity. To monitor natural taimen populations we have sequenced a portion (8,141 bp) of the mitochondrial (mt) genome in 28 specimens of H. taimen from six localities in the Amur River basin. Nucleotide variability is low (π = 0.0010), but structured in two divergent haplotype groups. A comparison of the data with the GenBank H. taimen mt genome (HQ897271) reveals significant differences between them in spite of the fact that the fish specimens come from neighboring geographical areas. The distribution of divergence is non-uniform with two highly pronounced divergent regions centered on two genes, ND3 and ND6. To clarify the pattern of divergence we sequenced the corresponding portion of the mt genome of lenok Brachymystax tumensis and analyzed the GenBank complete mt genomes of related species. We have found that the first and second divergent regions are identical between the GenBank H. taimen and two lenok subspecies, B. lenok and B. lenok tsinlingensis, respectively. Consequently, both divergent regions represent introgressed mtDNA resulting from intergeneric hybridization between the two lenok subspecies and H. taimen. Introgression is, however, not detected in our specimens. This plus the precise identity of the introgressed fragments between the donor and the recipient GenBank sequence suggests that the introgression is local and very recent, probably due to artificial manipulations involving taimen – lenok intergeneric hybridization. Human-mediated hybridization may become a major threat to aquatic biodiversity. Consequently we suggest that due attention needs to be given to this threat by means of responsible breeding program management, so as to prevent a potential spread of hybrid fishes that could jeopardize the resilience of locally adapted gene pools of the native H. taimen populations.

Biparental inheritance of organelles in Pelargonium: evidence for intergenomic recombination of mitochondrial DNA

While uniparental transmission of mtDNA is widespread and dominating in eukaryotes leaving mutation as the major source of genotypic diversity, recently, biparental inheritance of mitochondrial genes has been demonstrated in reciprocal crosses of Pelargonium zonale and P. inquinans. The thereby arising heteroplasmy carries the potential for recombination between mtDNAs of different descent, i.e. between the parental mitochondrial genomes. We have analyzed these Pelargonium hybrids for mitochondrial intergenomic recombination events by examining differences in DNA blot hybridization patterns of the mitochondrial genes atp1 and cob. Further investigation of these genes and their flanking regions using nucleotide sequence polymorphisms and PCR revealed DNA segments in the progeny, which contained both P. zonale and P. inquinans sequences suggesting an intergenomic recombination in hybrids of Pelargonium. This turns Pelargonium into an interesting subject for studies of recombination and evolutionary dynamics of mitochondrial genomes.

Nucleotide variation in the mitochondrial genome provides evidence for dual routes of postglacial recolonization and genetic recombination in the northeastern brook trout (Salvelinus fontinalis)

Levels and patterns of mitochondrial DNA (mtDNA) variation were examined to investigate the population structure and possible routes of postglacial recolonization of the world's northernmost native populations of brook trout (Salvelinus fontinalis), which are found in Labrador, Canada. We analyzed the sequence diversity of a 1960-bp portion of the mitochondrial genome (NADH dehydrogenase 1 gene and part of cytochrome oxidase 1) of 126 fish from 32 lakes distributed throughout seven regions of northeastern Canada. These populations were found to have low levels of mtDNA diversity, a characteristic trait of populations at northern extremes, with significant structuring at the level of the watershed. Upon comparison of northeastern brook trout sequences to the publicly available brook trout whole mitochondrial genome (GenBank AF154850), we infer that the GenBank sequence is from a fish whose mtDNA has recombined with that of Arctic charr (S. alpinus). The haplotype distribution provides evidence of two different postglacial founding groups contributing to present-day brook trout populations in the northernmost part of their range; the evolution of the majority of the haplotypes coincides with the timing of glacier retreat from Labrador. Our results exemplify the strong influence that historical processes such as glaciations have had on shaping the current genetic structure of northern species such as the brook trout.

Large-scale asymmetric introgression of cytoplasmic DNA reveals Holocene range displacement in a North American boreal pine complex

Jack pine (Pinus banksiana) and lodgepole pine (Pinus contorta var. latifolia) are two North American boreal hard pines that hybridize in their zone of contact in western Canada. The main objective of this study was to characterize their patterns of introgression resulting from past and recent gene flow, using cytoplasmic markers having maternal or paternal inheritance. Mitochondrial DNA (mtDNA) and chloroplast DNA (cpDNA) diversity was assessed in allopatric populations of each species and in stands from the current zone of contact containing morphological hybrids. Cluster analyses were used to identify genetic discontinuities among groups of populations. A canonical analysis was also conducted to detect putative associations among cytoplasmic DNA variation, tree morphology, and site ecological features. MtDNA introgression was extensive and asymmetric: it was detected in P. banksiana populations from the hybrid zone and from allopatric areas, but not in P. contorta populations. Very weak cpDNA introgression was observed, and only in P. banksiana populations. The mtDNA introgression pattern indicated that central Canada was first colonized by migrants from a P. contorta glacial population located west of the Rocky Mountains, before being replaced by P. banksiana migrating westward during the Holocene. In contrast, extensive pollen gene flow would have erased the cpDNA traces of this ancient presence of P. contorta. Additional evidence for this process was provided by the results of canonical analysis, which indicated that the current cpDNA background of trees reflected recent pollen gene flow from the surrounding dominant species rather than historical events that took place during the postglacial colonization.

Genetic divergence, range expansion and possible homoploid hybrid speciation among pine species in Northeast China

Although homoploid hybrid speciation in plants is probably more common than previously realized, there are few well-documented cases of homoploid hybrid origin in conifers. We examined genetic divergence between two currently widespread pines in Northeast China, Pinus sylvestris var. mongolica and Pinus densiflora, and also whether two narrowly distributed pines in the same region, Pinus funebris and Pinus takahasii, might have originated from the two widespread species by homoploid hybrid speciation. Our results, based on population genetic analysis of chloroplast (cp), mitochondrial (mt) DNA, and nuclear gene sequence variation, showed that the two widespread species were divergent for both cp- and mtDNA variation, and also for haplotype variation at two of eight nuclear gene loci surveyed. Our analysis further indicated that P. sylvestris var. mongolica and P. densiflora remained allopatric during the most severe Quaternary glacial period that occurred in Northeast China, but subsequently exhibited rapid range expansions. P. funebris and P. takahasii, were found to contain a mixture of chlorotypes and nuclear haplotypes that distinguish P. sylvestris var. mongolica and P. densiflora, in support of the hypothesis that they possibly originated via homoploid hybrid speciation following secondary contact and hybridization between P. sylvestris var. mongolica and P. densiflora.

Numts help to reconstruct the demographic history of the ocellated lizard (Lacerta lepida) in a secondary contact zone

In northwestern Iberia, two largely allopatric Lacerta lepida mitochondrial lineages occur, L5 occurring to the south of Douro River and L3 to the north, with a zone of putative secondary contact in the region of the Douro River valley. Cytochrome b sequence chromatograms with polymorphisms at nucleotide sites diagnostic for the two lineages were detected in individuals in the region of the Douro River and further north within the range of L3. We show that these polymorphisms are caused by the presence of four different numts (I-IV) co-occurring with the L3 genome, together with low levels of heteroplasmy. Two of the numts (I and II) are similar to the mitochondrial genome of L5 but are quite divergent from the mitochondrial genome of L3 where they occur. We show that these numts are derived from the mitochondrial genome of L5 and were incorporated in L3 through hybridization at the time of secondary contact between the lineages. The additional incidence of these numts to the north of the putative contact zone is consistent with an earlier postglacial northward range expansion of L5, preceding that of L3. We show that genetic exchange between the lineages responsible for the origin of these numts in L3 after secondary contact occurred prior to, or coincident with, the northward expansion of L3. This study shows that, in the context of phylogeographic analysis, numts can provide evidence for past demographic events and can be useful tools for the reconstruction of complex evolutionary histories.

Tracing the recombination and colonization history of hybrid species in space and time

Hybrid speciation has long fascinated evolutionary biologists and laymen alike, presumably because it challenges our classical view of evolution as a 'one-way street' leading to strictly tree-like patterns of ancestry and descent. Homoploid hybrid speciation (HHS) has been a particularly interesting puzzle, as it appears to occur extremely rapidly, perhaps within less than 50 generations (McCarthy et al. 1995; Buerkle et al. 2000). Nevertheless, HHS may sometimes involve extended or repeated periods of recombination and gene exchange between populations subject to strong divergent natural selection (Buerkle & Rieseberg 2008). Thus, HHS provides a highly interesting setting for understanding the drivers and tempo of adaptive divergence and speciation in the face of gene flow (Arnold 2006; Rieseberg & Willis 2007; Nolte & Tautz 2009). In the present issue of Molecular Ecology, Wang et al. (2011) explore a particularly challenging issue connected to HHS: they attempt to trace the colonization and recombination history of an ancient (several MYA) hybrid species, from admixture and recombination in the ancestral hybrid zone to subsequent range shifts triggered by tectonic events (uplift of the Tibetan plateau) and climatic shifts (Pleistocene ice ages). This work is important because it addresses key issues related to the origin of the standing genetic variation available for adaptive responses (e.g. to climate change) and speciation in temperate species, which are topics of great current interest (Rieseberg et al. 2003; Barrett & Schluter 2008; de Carvalho et al. 2010).

Colonization of the Tibetan Plateau by the homoploid hybrid pine Pinus densata

Pinus densata is an intriguingly successful homoploid hybrid species that occupies vast areas of the southeastern Tibetan Plateau in which neither of its parental species are present, but the colonization processes involved are poorly understood. To shed light on how this species colonized and became established on the plateau, we surveyed paternally inherited chloroplast (cp) and maternally inherited mitochondrial (mt) DNA variation within and among 54 populations of P. densata and its putative parental species throughout their respective ranges. Strong spatial genetic structure of both cp and mtDNA were detected in P. densata populations. Mitotypes specific to P. densata were likely generated by complex recombination events. A putative ancestral hybrid zone in the northeastern periphery of P. densata was identified, and we propose that the species then colonized the plateau by migrating westwards. Along the colonization route, consecutive bottlenecks and surfing of rare alleles caused a significant reduction in genetic diversity and strong population differentiation. The direction and intensity of introgression from parental species varied among geographic regions. In western parts of its range, the species seems to have been isolated from seed and pollen flow from its parent species for a long time. The observed spatial distribution of genetic diversity in P. densata also appears to reflect the persistence of this species on the plateau during the last glaciation. Our results indicate that both ancient and contemporary population dynamics have contributed to the spatial distribution of genetic diversity in P. densata, which accordingly reflects its evolutionary history.

A test of seven candidate barcode regions from the plastome in Picea (Pinaceae)

DNA barcoding, as a tool for species discrimination, has been used efficiently in animals, algae and fungi, but there are still debates on which DNA region(s) can be used as the standard barcode(s) for land plants. Gymnosperms, especially conifers, are important components of forests, and there is an urgent need for them to be identified through DNA barcoding because of their high frequency of collection in the field. However, the feasibility of DNA barcoding in gymnosperms has not been examined based on a dense species sampling. Here we selected seven candidate DNA barcodes from the plastome (matK, rbcL, rpoB, rpoC1, atpF-atpH, psbA-trnH, and psbK-psbI) to evaluate their suitability in Picea (spruce). The results showed that none of them or their different combinations has sufficient resolution for spruce species, although matK+rbcL might be used as a two-locus barcode. The low efficiency of these candidate barcodes in Picea might be caused by the paternal inheritance of the chloroplast genome, long generation time, recent radiation, and frequent inter-specific hybridization aided by wind pollination. Some of these factors could also be responsible for the difficulties in barcoding other plant groups. Furthermore, the potential of the nuclear LEAFY gene as a land plant barcode was discussed.

Gene flow and species delimitation: a case study of two pine species with overlapping distributions in southeast China

Species delimitation detected by molecular markers is complicated by introgression and incomplete lineage sorting between species. Recent modeling suggests that fixed genetic differences between species are highly related to rates of intraspecific gene flow. However, it remains unclear whether such differences are due to high levels of intraspecific gene flow overriding the spread of introgressed alleles or favoring rapid lineage sorting between species. In pines, chloroplast (cp) and mitochondrial (mt) DNAs are normally paternally and maternally inherited, respectively, and thus their relative rates of intraspecific gene flow are expected to be high and low, respectively. In this study, we used two pine species with overlapping geographical distributions in southeast China, P. massoniana and P. hwangshanensis, as a model system to examine the association between organelle gene flow and variation within and between species. We found that cpDNA variation across these two pine species is more species specific than mtDNA variation and almost delimits taxonomic boundaries. The shared mt/cp DNA genetic variation between species shows no bias in regard to parapatric versus allopatric species' distributions. Our results therefore support the hypothesis that high intraspecific gene flow has accelerated cpDNA lineage sorting between these two pine species.

Unraveling cryptic reticulate relationships and the origin of orphan hybrid disjunct populations in Narcissus

Evolutionary consequences of natural hybridization between species may vary so drastically depending on spatial, genetic, and ecological factors that multiple approaches are required to uncover them. To unravel the evolutionary history of a controversial hybrid (Narcissus x perezlarae), here we use four approaches: DNA sequences from five regions (four organellar, one nuclear), cytological studies (chromosome counts and genome size), crossing experiments, and niche modeling. We conclude that (1) it actually consists of two different hybrid taxa, N.xperezlarae s.s. (N. cavanillesii x N. miniatus) and N.xalentejanus (N. cavanillesii x N. serotinus); (2) both have been formed several times independently, that is, polytopically; (3) N. cavanillesii was the mother progenitor in most hybridization events. We also address the origin of orphan hybrid populations of N.xperezlarae in eastern Spain, hundreds of kilometers away from N. cavanillesii. Although long-distance dispersal of already formed hybrids cannot be completely rejected, extirpation of N. cavanillesii via demographic competition is a more likely explanation. Low-reproductive barriers to fertilization by foreign pollen in N. cavanillesii, molecular footprints of the former presence of this species in the area, active asexual propagation by bulbs in N.xperezlarae, and overlapping ecological niches are consistent with the extirpation scenario.

OrgConv: detection of gene conversion using consensus sequences and its application in plant mitochondrial and chloroplast homologs

Background

The ancestry of mitochondria and chloroplasts traces back to separate endosymbioses of once free-living bacteria. The highly reduced genomes of these two organelles therefore contain very distant homologs that only recently have been shown to recombine inside the mitochondrial genome. Detection of gene conversion between mitochondrial and chloroplast homologs was previously impossible due to the lack of suitable computer programs. Recently, I developed a novel method and have, for the first time, discovered recurrent gene conversion between chloroplast mitochondrial genes. The method will further our understanding of plant organellar genome evolution and help identify and remove gene regions with incongruent phylogenetic signals for several genes widely used in plant systematics. Here, I implement such a method that is available in a user friendly web interface.

Results

OrgConv (Organellar Conversion) is a computer package developed for detection of gene conversion between mitochondrial and chloroplast homologous genes. OrgConv is available in two forms; source code can be installed and run on a Linux platform and a web interface is available on multiple operating systems. The input files of the feature program are two multiple sequence alignments from different organellar compartments in FASTA format. The program compares every examined sequence against the consensus sequence of each sequence alignment rather than exhaustively examining every possible combination. Making use of consensus sequences significantly reduces the number of comparisons and therefore reduces overall computational time, which allows for analysis of very large datasets. Most importantly, with the significantly reduced number of comparisons, the statistical power remains high in the face of correction for multiple tests.

Conclusions

Both the source code and the web interface of OrgConv are available for free from the OrgConv website http://www.indiana.edu/~orgconv. Although OrgConv has been developed with main focus on detection of gene conversion between mitochondrial and chloroplast genes, it may also be used for detection of gene conversion between any two distinct groups of homologous sequences.

Heteroplasmy and stoichiometric complexity of plant mitochondrial genomes--though this be madness, yet there's method in't

Mitochondrial heteroplasmy is defined as the coexistence of divergent mitochondrial genotypes in a cell. The ratio of the alternative genomes may be variable, but in plants, the usually prevalent main genome is accompanied by sublimons--substoichiometric mitochondrial DNA (mtDNA) molecules. Plant mitochondrial heteroplasmy was originally viewed as being associated with pathological mutations or was found in non-natural plant populations. Currently, it is considered to be a common situation in plants. Recent years have changed the previous view on the role of homologous recombination, small-scale mutations, and paternal leakage of mtDNA in the generation of heteroplasmy. Newly developed sensitive techniques have allowed the precise estimation of mtDNA stoichiometry. Mechanisms of maintenance and transmission of heteroplasmic genomes, including DNA recombination and replication, as well as mitochondrial fusion and fission, have been studied. This review describes the high level of plant mitochondrial genome complication--the 'madness' resulting from the heteroplasmic state and explains the method hidden in this madness. Heteroplasmy is described as the evolutionary strategy of uniparentally inherited plant mitochondrial genomes which do not undergo sexual recombination. In order to compensate for this deficiency, alternative types of mtDNA are substoichiometrically accumulated as a reservoir of genetic variability and may undergo accelerated evolution. Occasionally, sublimons are selected and amplified in the process called substoichiometric shifting, to take over the role of the main genome. Alternative mitochondrial genomes may recombine, yielding new mtDNA variants, or segregate during plant growth resulting in plants with mosaic phenotypes. Two opposite roles of mitochondrial heteroplasmy with respect to acceleration or counteracting of mutation accumulation are also discussed. Finally, nuclear control of heteroplasmy and substoichiometric shifting is described.

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.

Incongruence among mitochondrial, chloroplast and nuclear gene trees in Pinus subgenus Strobus (Pinaceae)

Introgression has been considered to be one of main factors leading to phylogenetic incongruence among different datasets at lower taxonomic levels. In the plants of Pinaceae, the mtDNA, cpDNA, and nuclear DNA (nrDNA) may have different evolutionary histories through introgression because they are inherited maternally, paternally and biparentally, respectively. We compared mtDNA, cpDNA, and two low-copy nrDNA phylogenetic trees in the genus Pinus subgenus Strobus, in order to detect unknown past introgression events in this group. nrDNA trees were mostly congruent with the cpDNA tree, and supported the recent sectional and subsectional classification system. In contrast, mtDNA trees split the members of sect. Quinquefoliae into two groups that were not observed in the other gene trees. The factors constituting incongruence may be divided into the following two categories: the different splits within subsect. Strobus, and the non-monophyly of subsect. Gerardianae. The former was hypothesized to have been caused by the past introgression of cpDNA, mtDNA or both between Eurasian and North American species through Beringia. The latter was likely caused by the chimeric structure of the mtDNA sequence of P. bungeana, which might have originated through past hybridization, or through a horizontal transfer event and subsequent recombination.

Mitochondrial DNA recombination in a free-ranging Australian lizard

Mitochondrial DNA (mtDNA) is the traditional workhorse for reconstructing evolutionary events. The frequent use of mtDNA in such analyses derives from the apparent simplicity of its inheritance: maternal and lacking bi-parental recombination. However, in hybrid zones, the reproductive barriers are often not completely developed, resulting in the breakdown of male mitochondrial elimination mechanisms, leading to leakage of paternal mitochondria and transient heteroplasmy, resulting in an increased possibility of recombination. Despite the widespread occurrence of heteroplasmy and the presence of the molecular machinery necessary for recombination, we know of no documented example of recombination of mtDNA in any terrestrial wild vertebrate population. By sequencing the entire mitochondrial genome (16761bp), we present evidence for mitochondrial recombination in the hybrid zone of two mitochondrial haplotypes in the Australian frillneck lizard (Chlamydosaurus kingii).

Molecular phylogeny and biogeography of Picea (Pinaceae): Implications for phylogeographical studies using cytoplasmic haplotypes

The center of diversity is not necessarily the place of origin, as has been established by many plant molecular phylogenies. Picea is a complicated but very important genus in coniferous forests of the Northern Hemisphere, with a high species diversity in Asia. Its phylogeny and biogeography were investigated here using sequence analysis of the paternally inherited chloroplast trnC-trnD and trnT-trnF regions and the maternally inherited mitochondrial nad5 intron 1. We found that the North American P. breweriana and P. sitchensis were basal to the other spruces that were further divided into three clades in the cpDNA phylogeny, and that the New World species harbored four of five mitotypes detected, including two ancestral ones and three endemics. These results, combined with biogeographic analyses using DIVA and MacClade and fossil evidence, suggest that Picea originated in North America, and that its present distribution could stem from two times of dispersal from North America to Asia by the Beringian land bridge, and then from Asia to Europe. Most of the northeastern Asian species and the European P. abies could arise from a recent radiation given the very low interspecific genetic differentiation and pure mitotype of them. Considering that the ancestral mtDNA polymorphism can be preserved in many descendant species, even distantly related ones, we suggest that more species, at least the closely related ones, should be sampled in the phylogeographical study using cytoplasmic haplotypes if possible. In addition, we also discussed the evolution and phylogenetic utility of morphological characters in Picea.