Inheritance of chloroplast DNA is not strictly maternal in Silene vulgaris (Caryophyllaceae): evidence from experimental crosses and natural populations

Paternal leakage of plastids rescues inter-lineage hybrids in Silene nutans

BACKGROUND AND AIMS

Organelle genomes are usually maternally inherited in angiosperms. However, biparental inheritance has been observed, especially in hybrids resulting from crosses between divergent genetic lineages. When it concerns the plastid genome, this exceptional mode of inheritance might rescue inter-lineage hybrids suffering from plastid-nuclear incompatibilities. Genetically differentiated lineages of Silene nutans exhibit strong postzygotic isolation owing to plastid-nuclear incompatibilities, highlighted by inter-lineage hybrid chlorosis and mortality. Surviving hybrids can exhibit variegated leaves, which might indicate paternal leakage of the plastid genome. We tested whether the surviving hybrids inherited the paternal plastid genome and survived thanks to paternal leakage.

METHODS

We characterized the leaf phenotype (fully green, variegated or white) of 504 surviving inter-lineage hybrids obtained from a reciprocal cross experiment among populations of four genetic lineages (W1, W2, W3 and E1) of S. nutans from Western Europe and genotyped 560 leaf samples (both green and white leaves for variegated hybrids) using six lineage-specific plastid single nucleotide polymorphisms.

KEY RESULTS

A high proportion of the surviving hybrids (≤98 %) inherited the paternal plastid genome, indicating paternal leakage. The level of paternal leakage depended on cross type and cross direction. The E1 and W2 lineages as maternal lineages led to the highest hybrid mortality and to the highest paternal leakage from W1 and W3 lineages in the few surviving hybrids. This was consistent with E1 and W2 lineages, which contained the most divergent plastid genomes. When W3 was the mother, more hybrids survived, and no paternal leakage was detected.

CONCLUSIONS

By providing a plastid genome potentially more compatible with the hybrid nuclear background, paternal leakage has the potential to rescue inter-lineage hybrids from plastid-nuclear incompatibilities. This phenomenon might slow down the speciation process, provided hybrid survival and reproduction can occur in the wild.

Phylogenomic conflict analyses of the plastid and mitochondrial genomes via deep genome skimming highlight their independent evolutionary histories: A case study in the cinquefoil genus Potentilla sensu lato (Potentilleae, Rosaceae)

Phylogenomic conflicts are widespread among genomic data, with most previous studies primarily focusing on nuclear datasets instead of organellar genomes. In this study, we investigate phylogenetic conflict analyses within and between plastid and mitochondrial genomes using Potentilla as a case study. We generated three plastid datasets (coding, noncoding, and all-region) and one mitochondrial dataset (coding regions) to infer phylogenies based on concatenated and multispecies coalescent (MSC) methods. Conflict analyses were then performed using PhyParts and Quartet Sampling (QS). Both plastid and mitochondrial genomes divided the Potentilla into eight highly supported clades, two of which were newly identified in this study. While most organellar loci were uninformative for the majority of nodes (bootstrap value < 70%), PhyParts and QS detected conflicting signals within the two organellar genomes. Regression analyses revealed that conflict signals mainly occurred among shorter loci, whereas longer loci tended to be more concordant with the species tree. In addition, two significant disagreements between the two organellar genomes were detected, likely attributed to hybridization and/or incomplete lineage sorting. Our results demonstrate that mitochondrial genes can fully resolve the phylogenetic relationships among eight major clades of Potentilla and are not always linked with plastome in evolutionary history. Stochastic inferences appear to be the primary source of observed conflicts among the gene trees. We recommend that the loci with short sequence length or containing limited informative sites should be used cautiously in MSC analysis, and suggest the joint application of concatenated and MSC methods for phylogenetic inference using organellar genomes.

Complete Plastid Genomes of Nine Species of Ranunculeae (Ranunculaceae) and Their Phylogenetic Inferences

The tribe Ranunculeae, Ranunculaceae, comprising 19 genera widely distributed all over the world. Although a large number of Sanger sequencing-based molecular phylogenetic studies have been published, very few studies have been performed on using genomic data to infer phylogenetic relationships within Ranunculeae. In this study, the complete plastid genomes of nine species (eleven samples) from Ceratocephala, Halerpestes, and Ranunculus were de novo assembled using a next-generation sequencing method. Previously published plastomes of Oxygraphis and other related genera of the family were downloaded from GenBank for comparative analysis. The complete plastome of each Ranunculeae species has 112 genes in total, including 78 protein-coding genes, 30 transfer RNA genes, and four ribosomal RNA genes. The plastome structure of Ranunculeae samples is conserved in gene order and arrangement. There are no inverted repeat (IR) region expansions and only one IR contraction was found in the tested samples. This study also compared plastome sequences across all the samples in gene collinearity, codon usage, RNA editing sites, nucleotide variability, simple sequence repeats, and positive selection sites. Phylogeny of the available Ranunculeae species was inferred by the plastome data using maximum-likelihood and Bayesian inference methods, and data partitioning strategies were tested. The phylogenetic relationships were better resolved compared to previous studies based on Sanger sequencing methods, showing the potential value of the plastome data in inferring the phylogeny of the tribe.

Untying the Gordian knot of plastid phylogenomic conflict: A case from ferns

Phylogenomic studies based on plastid genome have resolved recalcitrant relationships among various plants, yet the phylogeny of Dennstaedtiaceae at the level of family and genera remains unresolved due to conflicting plastid genes, limited molecular data and incomplete taxon sampling of previous studies. The present study generated 30 new plastid genomes of Dennstaedtiaceae (9 genera, 29 species), which were combined with 42 publicly available plastid genomes (including 24 families, 27 genera, 42 species) to explore the evolution of Dennstaedtiaceae. In order to minimize the impact of systematic errors on the resolution of phylogenetic inference, we applied six strategies to generate 30 datasets based on CDS, intergenic spacers, and whole plastome, and two tree inference methods (maximum-likelihood, ML; and multispecies coalescent, MSC) to comprehensively analyze the plastome-scale data. Besides, the phylogenetic signal among all loci was quantified for controversial nodes using ML framework, and different topologies hypotheses among all datasets were tested. The species trees based on different datasets and methods revealed obvious conflicts at the base of the polypody ferns. The topology of the "CDS-codon-align-rm3" (CDS with the removal of the third codon) matrix was selected as the primary reference or summary tree. The final phylogenetic tree supported Dennstaedtiaceae as the sister group to eupolypods, and Dennstaedtioideae was divided into four clades with full support. This robust reconstructed phylogenetic backbone establishes a framework for future studies on Dennstaedtiaceae classification, evolution and diversification. The present study suggests considering plastid phylogenomic conflict when using plastid genomes. From our results, reducing saturated genes or sites can effectively mitigate tree conflicts for distantly related taxa. Moreover, phylogenetic trees based on amino acid sequences can be used as a comparison to verify the confidence of nucleotide-based trees.

Plastaumatic: Automating plastome assembly and annotation

Plastome sequence data is most often extracted from plant whole genome sequencing data and need to be assembled and annotated separately from the nuclear genome sequence. In projects comprising multiple genomes, it is labour intense to individually process the plastomes as it requires many steps and software. This study developed Plastaumatic - an automated pipeline for both assembly and annotation of plastomes, with the scope of the researcher being able to load whole genome sequence data with minimal manual input, and therefore a faster runtime. The main structure of the current automated pipeline includes trimming of adaptor and low-quality sequences using fastp, de novo plastome assembly using NOVOPlasty, standardization and quality checking of the assembled genomes through a custom script utilizing BLAST+ and SAMtools, annotation of the assembled genomes using AnnoPlast, and finally generating the required files for NCBI GenBank submissions. The pipeline is demonstrated with 12 potato accessions and three soybean accessions.

Phylogeographic structure of common sage (Salvia officinalis L.) reveals microrefugia throughout the Balkans and colonizations of the Apennines

Studying the population-genetic and phylogeographic structures of a representative species of a particular geographical region can not only provide us with information regarding its evolutionary history, but also improve our understanding of the evolutionary processes underlying the patterns of species diversity in that area. By analysing eight highly polymorphic microsatellite loci and two chloroplast DNA regions, we have investigated the influence of Pleistocene climate fluctuations on the evolutionary history of Salvia officinalis L. (common sage). The populations with the highest genetic diversity were located in the central parts of the Balkan distribution range. A large group of closely related haplotypes was distributed throughout the Balkans and the central Apennines, while the private lineage occupied the southern Apennines. In addition, two highly differentiated lineages were scattered only over the Balkans. The results suggest that a single refugium of the studied species from the last glacial period was located in the central part of the range in the Balkans. Numerous microrefugia, probably spanning several glaciation cycles, were scattered across the Balkans, while colonisation of the Apennines from the Balkans occurred at least on two occasions.

Comparative chloroplast genomes and phylogenetic analyses of Pinellia

Background

Pinellia Tenore (Araceae) is a genus of perennial herbaceous plants, all of which have medicinal value. The chloroplast (cp) genome data of Pinellia are scarce, and the phylogenetic relationship and gene evolution remain unclear.

Methods and results

We sequenced and annotated the Pinellia pedatisecta cp genome and combined it with previously published genomes for other Pinellia species. We used bioinformatics methods to analyse the genomic structure, repetitive sequences, interspecific variation, divergence hotspots, phylogenetic relationships, divergence time estimation and selective pressure of four Pinellia plastomes. Results showed that the cp genomes of Pinellia varied in length between 168,178 (P. pedatisecta MN046890) and 164,013 bp (P. ternata KR270823). A total of 68–111 SSR loci were identified as candidate molecular markers for further genetic diversity study. Eight mutational hotspot regions were determined, including psbI-trnG-UCC, psbM-rpoB, ndhJ-trnT-UGU, trnP-UGG-trnW-CCA, ndhF-trnN-GUU, ndhG-ndhE, ycf1-rps15 and trnR-ycf1. Gene selection pressure suggested that four genes were subjected to positive selection. Phylogenetic inferences based on the complete cp genomes revealed a sister relationship between Pinellia and Arisaema plants whose divergence was estimated to occur around 22.48 million years ago. All Pinellia species formed a monophyletic evolutionary clade in which P. peltata, rather than P. pedatisecta, earlier diverged, indicating that P. pedatisecta is not the basal taxon of Pinellia but P. peltata may be.

Conclusions

The cp genomes of Pinellia will provide valuable information for species classification, identification, molecular breeding and evolutionary exploration of the genus Pinellia.

Supplementary Information

The online version of this article (10.1007/s11033-022-07617-5) contains supplementary material, which is available to authorized users.

Disentangling Complex Inheritance Patterns of Plant Organellar Genomes: An Example From Carrot

Plant mitochondria and plastids display an array of inheritance patterns and varying levels of heteroplasmy, where individuals harbor more than 1 version of a mitochondrial or plastid genome. Organelle inheritance in plants has the potential to be quite complex and can vary with plant growth, development, and reproduction. Few studies have sought to investigate these complicated patterns of within-individual variation and inheritance using experimental crosses in plants. We carried out crosses in carrot, Daucus carota L. (Apiaceae), which has previously been shown to exhibit organellar heteroplasmy. We used mitochondrial and plastid markers to begin to disentangle the patterns of organellar inheritance and the fate of heteroplasmic variation, with special focus on cases where the mother displayed heteroplasmy. We also investigated heteroplasmy across the plant, assaying leaf samples at different development stages and ages. Mitochondrial and plastid paternal leakage was rare and offspring received remarkably similar heteroplasmic mixtures to their heteroplasmic mothers, indicating that heteroplasmy is maintained over the course of maternal inheritance. When offspring did differ from their mother, they were likely to exhibit a loss of the genetic variation that was present in their mother. Finally, we found that mitochondrial variation did not vary significantly over plant development, indicating that substantial vegetative sorting did not occur. Our study is one of the first to quantitatively investigate inheritance patterns and heteroplasmy in plants using controlled crosses, and we look forward to future studies making use of whole genome information to study the complex evolutionary dynamics of plant organellar genomes.

Phylogenomic framework of the IRLC legumes (Leguminosae subfamily Papilionoideae) and intercontinental biogeography of tribe Wisterieae

The inverted repeat-lacking clade (IRLC) is one of the most derived clades within the subfamily Papilionoideae of the legume family, and includes various economically important plants, e.g., chickpeas, peas, liquorice, and the largest genus of angiosperms, Astragalus. Tribe Wisterieae is one of the earliest diverged groups of the IRLC, and its generic delimitation and spatiotemporal diversification needs further clarifications. Based on genome skimming data, we herein reconstruct the phylogenomic framework of the IRLC, and infer the inter-generic relationships and historical biogeography of Wisterieae. We redefine tribe Caraganeae to contain Caragana only, and tribe Astragaleae is reduced to the Erophaca-Astragalean clade. The chloroplast capture scenario was hypothesized as the most plausible explanation of the topological incongruences between the chloroplast CDSs and nuclear ribosomal DNA trees in both the Glycyrrhizinae-Adinobotrys-Wisterieae clade and the Chesneyeae-Caraganeae-Hedysareae clade. A new name, Caragana lidou L. Duan & Z.Y. Chang, is proposed within Caraganeae. Thirteen genera are herein supported within Wisterieae, including a new genus, Villosocallerya L. Duan, J. Compton & Schrire, segregated from Callerya. Our biogeographic analyses suggest that Wisterieae originated in the late Eocene and its most recent common ancestor (MRCA) was distributed in continental southeastern Asia. Lineages of Wisterieae remained in the ancestral area from the early Oligocene to the early Miocene. By the middle Miocene, Whitfordiodendron and the MRCA of Callerya-Kanburia-Villosocallerya Clade became disjunct between the Sunda area and continental southeastern Asia, respectively; the MRCA of Wisteria migrated to North America via the Bering land bridge. The ancestor of Austrocallerya and Padbruggea migrated to the Wallacea-Oceania area, which split in the early Pliocene. In the Pleistocene, Wisteria brachybotrys, W. floribunda and Wisteriopsis japonica reached Japan, and Callerya cinerea dispersed to South Asia. This study provides a solid phylogenomic for further evolutionary/biogeographic/systematic investigations on the ecologically diverse and economically important IRLC legumes.

Clade-Specific Plastid Inheritance Patterns Including Frequent Biparental Inheritance in Passiflora Interspecific Crosses

Plastid inheritance in angiosperms is presumed to be largely maternal, with the potential to inherit plastids biparentally estimated for about 20% of species. In Passiflora, maternal, paternal and biparental inheritance has been reported; however, these studies were limited in the number of crosses and progeny examined. To improve the understanding of plastid transmission in Passiflora, the progeny of 45 interspecific crosses were analyzed in the three subgenera: Passiflora, Decaloba and Astrophea. Plastid types were assessed following restriction digestion of PCR amplified plastid DNA in hybrid embryos, cotyledons and leaves at different developmental stages. Clade-specific patterns of inheritance were detected such that hybrid progeny from subgenera Passiflora and Astrophea predominantly inherited paternal plastids with occasional incidences of maternal inheritance, whereas subgenus Decaloba showed predominantly maternal and biparental inheritance. Biparental plastid inheritance was also detected in some hybrids from subgenus Passiflora. Heteroplasmy due to biparental inheritance was restricted to hybrid cotyledons and first leaves with a single parental plastid type detectable in mature plants. This indicates that in Passiflora, plastid retention at later stages of plant development may not reflect the plastid inheritance patterns in embryos. Passiflora exhibits diverse patterns of plastid inheritance, providing an excellent system to investigate underlying mechanisms in angiosperms.

Effective double-digest RAD sequencing and genotyping despite large genome size

Obtaining informative data is the ambition of any genomic project, but in nonmodel species with very large genomes, pursuing such a goal requires surmounting a series of analytical challenges. Double-digest RAD sequencing is routinely used in nonmodel organisms and offers some control over the volume of data obtained. However, the volume of data recovered is not always an indication of the reliability of data sets, and quality checks are necessary to ensure that true and artefactual information is set apart. In the present study, we aim to fill the gap existing between the known applicability of RAD sequencing methods in plants with large genomes and the use of the retrieved loci for population genetic inference. By analysing two populations of Cypripedium calceolus, a nonmodel orchid species with a large genome size (1C ~ 31.6 Gbp), we provide a complete workflow from library preparation to bioinformatic filtering and inference of genetic diversity and differentiation. We show how filtering strategies to dismiss potentially misleading data need to be explored and adapted to data set-specific features. Moreover, we suggest that the occurrence of organellar sequences in libraries should not be neglected when planning the experiment and analysing the results. Finally, we explain how, in the absence of prior information about the genome of the species, seeking high standards of quality during library preparation and sequencing can provide an insurance against unpredicted technical or biological constraints.

Guidelines for the Choice of Sequences for Molecular Plant Taxonomy

This chapter presents an overview of the major plant DNA sequences and molecular methods available for plant taxonomy. Guidelines are provided for the choice of sequences and methods to be used, based on the DNA compartment (nuclear, chloroplastic, mitochondrial), evolutionary mechanisms, and the level of taxonomic differentiation of the plants under survey.

Complete Chloroplast Genome Characterization of Oxalis Corniculata and Its Comparison with Related Species from Family Oxalidaceae

Oxalis corniculata L. (family Oxalidaceae) is a small creeper wood sorrel plant that grows well in moist climates. Despite being medicinally important, little is known about the genomics of this species. Here, we determined the complete chloroplast genome sequence of O. corniculata for the first time and compared it with other members of family Oxalidaceae. The genome was 152,189 bp in size and comprised of a pair of 25,387 bp inverted repeats (IR) that separated a large 83,427 bp single copy region (LSC) and a small 16,990 bp single copy region (SSC). The chloroplast genome of O. corniculata contains 131 genes with 83 protein coding genes, 40 tRNA genes, and 8 rRNA genes. The analysis revealed 46 microsatellites, of which 6 were present in coding sequences (CDS) regions, 34 in the LSC, 8 in the SSC, and 2 in the single IR region. Twelve palindromic repeats, 30 forward repeats, and 32 tandem repeats were also detected. Chloroplast genome comparisons revealed an overall high degree of sequence similarity between O. corniculata and O. drummondii and some divergence in the intergenic spacers of related species in Oxalidaceae. Furthermore, the seven most divergent genes (ccsA, clpP, rps8, rps15, rpl22, matK, and ycf1) among genomes were observed. Phylogenomic characterization on the basis of 60 shared genes revealed that O. corniculata is closely related to O. drummondii. The complete O. corniculata genome sequenced in the present study is a valuable resource for investigating the population and evolutionary genetics of family Oxalidaceae and can be used to identify related species.

Genetic Characterization of an Endangered Chilean Endemic Species, Prosopis burkartii Muñoz, Reveals its Hybrids Parentage

The hybridization of Prosopis burkartii, a critically endangered endemic species, and the identification of its paternal species has not been genetically studied before. In this study we aimed to genetically confirm the origin of this species. To resolve the parental status of P. burkartii, inter-simple sequence repeat (ISSR), simple sequence repeats (SSR) and intron trnL molecular markers were used, and compared with Chilean species from the Algarobia and Strombocarpa sections. Out of seven ISSRs, a total of 70 polymorphic bands were produced in four species of the Strombocarpa section. An Multi-dimensional scaling (MDS) and Bayasian (STRUCTURE) analysis showed signs of introgression of genetic material in P. burkartii. Unweighted pair group method with arithmetic average (UPGMA) cluster analysis showed three clusters, and placed the P. burkartii cluster nested within the P. tamarugo group. Sequencing of the trnL intron showed a fragment of 535 bp and 529 bp in the species of the Algarobia and Strombocarpa sections, respectively. Using maximum parsimony (MP) and maximum likelihood (ML) trees with the trnL intron, revealed four clusters. A species-specific diagnostic method was performed, using the trnL intron Single Nucleotide Polymorphism (SNP). This method identified if individuals of P. burkartii inherited their maternal DNA from P. tamarugo or from P. strombulifera. We deduced that P. tamarugo and P. strombulifera are involved in the formation of P. burkartii.

Intra-individual heteroplasmy in the Gentiana tongolensis plastid genome (Gentianaceae)

Chloroplasts are typically inherited from the female parent and are haploid in most angiosperms, but rare intra-individual heteroplasmy in plastid genomes has been reported in plants. Here, we report an example of plastome heteroplasmy and its characteristics in Gentiana tongolensis (Gentianaceae). The plastid genome of G. tongolensis is 145,757 bp in size and is missing parts of petD gene when compared with other Gentiana species. A total of 112 single nucleotide polymorphisms (SNPs) and 31 indels with frequencies of more than 2% were detected in the plastid genome, and most were located in protein coding regions. Most sites with SNP frequencies of more than 10% were located in six genes in the LSC region. After verification via cloning and Sanger sequencing at three loci, heteroplasmy was identified in different individuals. The cause of heteroplasmy at the nucleotide level in plastome of G. tongolensis is unclear from the present data, although biparental plastid inheritance and transfer of plastid DNA seem to be most likely. This study implies that botanists should reconsider the heredity and evolution of chloroplasts and be cautious with using chloroplasts as genetic markers, especially in Gentiana.

Characterizing gene tree conflict in plastome-inferred phylogenies

Evolutionary relationships among plants have been inferred primarily using chloroplast data. To date, no study has comprehensively examined the plastome for gene tree conflict. Using a broad sampling of angiosperm plastomes, we characterize gene tree conflict among plastid genes at various time scales and explore correlates to conflict (e.g., evolutionary rate, gene length, molecule type). We uncover notable gene tree conflict against a backdrop of largely uninformative genes. We find alignment length and tree length are strong predictors of concordance, and that nucleotides outperform amino acids. Of the most commonly used markers, matK, greatly outperforms rbcL; however, the rarely used gene rpoC2 is the top-performing gene in every analysis. We find that rpoC2 reconstructs angiosperm phylogeny as well as the entire concatenated set of protein-coding chloroplast genes. Our results suggest that longer genes are superior for phylogeny reconstruction. The alleviation of some conflict through the use of nucleotides suggests that stochastic and systematic error is likely the root of most of the observed conflict, but further research on biological conflict within plastome is warranted given documented cases of heteroplasmic recombination. We suggest that researchers should filter genes for topological concordance when performing downstream comparative analyses on phylogenetic data, even when using chloroplast genomes.

Complete chloroplast genomes of medicinally important Teucrium species and comparative analyses with related species from Lamiaceae

Teucrium is one of the most economically and ecologically important genera in the Lamiaceae family; however, it is currently the least well understood at the plastome level. In the current study, we sequenced the complete chloroplast (cp) genomes of T. stocksianum subsp. stenophyllum R.A.King (TSS), T. stocksianum subsp. stocksianum Boiss. (TS) and T. mascatense Boiss. (TM) through next-generation sequencing and compared them with the cp genomes of related species in Lamiaceae (Ajuga reptans L., Caryopteris mongholica Bunge, Lamium album L., Lamium galeobdolon (L.) Crantz, and Stachys byzantina K.Koch). The results revealed that the TSS, TS and TM cp genomes have sizes of 150,087, 150,076 and 150,499 bp, respectively. Similarly, the large single-copy (LSC) regions of TSS, TS and TM had sizes of 81,707, 81,682 and 82,075 bp, respectively. The gene contents and orders of these genomes were similar to those of other angiosperm species. However, various differences were observed at the inverted repeat (IR) junctions, and the extent of the IR expansion into ψrps19 was 58 bp, 23 bp and 61 bp in TSS, TS and TM, respectively. Similarly, in all genomes, the pbsA gene was present in the LSC at varying distances from the JLA (IRa-LSC) junction. Furthermore, 89, 72, and 92 repeats were identified in the TSS, TM and TS cp genomes, respectively. The highest number of simple sequence repeats was found in TSS (128), followed by TS (127) and TM (121). Pairwise alignments of the TSS cp genome with related cp genomes showed a high degree of synteny. However, relatively lower sequence identity was observed when various coding regions were compared to those of related cp genomes. The average pairwise divergence among the complete cp genomes showed that TSS was more divergent from TM (0.018) than from TS (0.006). The current study provides valuable genomic insight into the genus Teucrium and its subspecies that may be applied to a more comprehensive study.

Selfish Mitonuclear Conflict

Mitochondria, a nearly ubiquitous feature of eukaryotes, are derived from an ancient symbiosis. Despite billions of years of cooperative coevolution - in what is arguably the most important mutualism in the history of life - the persistence of mitochondrial genomes also creates conditions for genetic conflict with the nucleus. Because mitochondrial genomes are present in numerous copies per cell, they are subject to both within- and among-organism levels of selection. Accordingly, 'selfish' genotypes that increase their own proliferation can rise to high frequencies even if they decrease organismal fitness. It has been argued that uniparental (often maternal) inheritance of cytoplasmic genomes evolved to curtail such selfish replication by minimizing within-individual variation and, hence, within-individual selection. However, uniparental inheritance creates conditions for cytonuclear conflict over sex determination and sex ratio, as well as conditions for sexual antagonism when mitochondrial variants increase transmission by enhancing maternal fitness but have the side-effect of being harmful to males (i.e., 'mother's curse'). Here, we review recent advances in understanding selfish replication and sexual antagonism in the evolution of mitochondrial genomes and the mechanisms that suppress selfish interactions, drawing parallels and contrasts with other organelles (plastids) and bacterial endosymbionts that arose more recently. Although cytonuclear conflict is widespread across eukaryotes, it can be cryptic due to nuclear suppression, highly variable, and lineage-specific, reflecting the diverse biology of eukaryotes and the varying architectures of their cytoplasmic genomes.

Molecular analysis indicates high levels of carabid weed seed consumption in cereal fields across Central Europe

Carabid beetles are abundant in temperate agroecosystems and can play a pivotal role as biocontrol agents. While there is good knowledge regarding their effects on invertebrate pests in some systems, comparably little is known on the rate of seed feeding under field conditions. Molecular approaches are ideally suited for investigating carabid feeding interactions; to date, however, they have only been applied to animal prey. We sampled adult carabid beetles in organic cereal fields in three regions along a Central European transect. Regurgitates from populations of the three most common species, Poecilus cupreus, Pseudoophonus rufipes and Pterostichus melanarius, were screened for plant DNA, cereal aphids, collembolans and earthworms. The frequency of carabid individuals positive for plant DNA was high (> 70%) and independent of carabid species, sex, region and the time point of sampling. Detections for non-pest and pest prey were comparably lower, with 21.6% for collembolans, 18.1% for earthworms and 4.2% for aphids, respectively. Despite the prolonged detection period of plant DNA in carabid guts, as compared to animal prey, these first results suggest that weed seeds form an important part of the adult carabid diet. It would also lend support to the hypothesis that seed-feeding carabids are biocontrol agents of weeds, with effects of regulation on the weed seedbank that depend on behavioural and contextual factors including carabid species preferences for weed seed species, their life stage and tillage practices.

Characterizing gene tree conflict in plastome-inferred phylogenies

Evolutionary relationships among plants have been inferred primarily using chloroplast data. To date, no study has comprehensively examined the plastome for gene tree conflict. Using a broad sampling of angiosperm plastomes, we characterize gene tree conflict among plastid genes at various time scales and explore correlates to conflict (e.g., evolutionary rate, gene length, molecule type). We uncover notable gene tree conflict against a backdrop of largely uninformative genes. We find alignment length and tree length are strong predictors of concordance, and that nucleotides outperform amino acids. Of the most commonly used markers, matK, greatly outperforms rbcL; however, the rarely used gene rpoC2 is the top-performing gene in every analysis. We find that rpoC2 reconstructs angiosperm phylogeny as well as the entire concatenated set of protein-coding chloroplast genes. Our results suggest that longer genes are superior for phylogeny reconstruction. The alleviation of some conflict through the use of nucleotides suggests that stochastic and systematic error is likely the root of most of the observed conflict, but further research on biological conflict within plastome is warranted given documented cases of heteroplasmic recombination. We suggest that researchers should filter genes for topological concordance when performing downstream comparative analyses on phylogenetic data, even when using chloroplast genomes.

Phylogeographic and phylogenetic analysis for Tripterygium species delimitation

Tripterygium wilfordii (Celastraceae) is a traditional Chinese medicine; and the dried root and rhizome constitute the main officinal parts. Tripterygium wilfordii has been identified as a potential candidate for the treatment of systemic lupus erythematosus, rheumatoid arthritis, nephritis, asthma, leprosy, and cancer. The phylogenetic relationships within the Tripterygium genus are ambiguous; thus, our aim is to clarify the relationships within this genus using phylogeographic and phylogenetic analyses. Here, we first sequenced three plastid DNA regions (i.e., psbA‐trnH, rpl32‐trnL, and trnL‐trnF) and found that Tripterygium hypoglaucum and T. wilfordii were clustered together based on the strength of the topology in the phylogenetic analysis: T. hypoglaucum is polyphyletic, and T. wilfordii is paraphyletic. A spatial analysis of molecular variance showed that the best group value is 4, and the groups were almost consistent with the topology of in the phylogenetic analysis. The Mantel analyses of Tripterygium using IBD web showed statistically significant relationships between genetic and geographical distance distributions (r = .3479, p < .0001). The molecular dating using Fossil calibration indicated that the divergence in Tripterygium was approximately 8.13 Ma. Furthermore, we also analyzed four DNA regions (i.e., ITS2, psbA‐trnH, matK, and rbcL) that were obtained from the NCBI nucleotide database; these results showed that T. wilfordii and T. hypoglaucum clustered together, while Tripterygium regelii represented a separate cluster. Tripterygium hypoglaucum and T. wilfordii were never distinct lineages, and the species circumscriptions are artificial. We propose that T. wilfordii and T. hypoglaucum are conspecific, while T. regelii likely constitutes a separate species.

Complete chloroplast genomes from apomictic Taraxacum (Asteraceae): Identity and variation between three microspecies

Chloroplast DNA sequences show substantial variation between higher plant species, and less variation within species, so are typically excellent markers to investigate evolutionary, population and genetic relationships and phylogenies. We sequenced the plastomes of Taraxacum obtusifrons Markl. (O978); T. stridulum Trávniček ined. (S3); and T. amplum Markl. (A978), three apomictic triploid (2n = 3x = 24) dandelions from the T. officinale agg. We aimed to characterize the variation in plastomes, define relationships and correlations with the apomictic microspecies status, and refine placement of the microspecies in the evolutionary or phylogenetic context of the Asteraceae. The chloroplast genomes of accessions O978 and S3 were identical and 151,322 bp long (where the nuclear genes are known to show variation), while A978 was 151,349 bp long. All three genomes contained 135 unique genes, with an additional copy of the trnF-GGA gene in the LSC region and 20 duplicated genes in the IR region, along with short repeats, the typical major Inverted Repeats (IR1 and IR2, 24,431bp long), and Large and Small Single Copy regions (LSC 83,889bp and SSC 18,571bp in O978). Between the two Taraxacum plastomes types, we identified 28 SNPs. The distribution of polymorphisms suggests some parts of the Taraxacum plastome are evolving at a slower rate. There was a hemi-nested inversion in the LSC region that is common to Asteraceae, and an SSC inversion from ndhF to rps15 found only in some Asteraceae lineages. A comparative repeat analysis showed variation between Taraxacum and the phylogenetically close genus Lactuca, with many more direct repeats of 40bp or more in Lactuca (1% larger plastome than Taraxacum). When individual genes and non-coding regions were for Asteraceae phylogeny reconstruction, not all showed the same evolutionary scenario suggesting care is needed for interpretation of relationships if a limited number of markers are used. Studying genotypic diversity in plastomes is important to characterize the nature of evolutionary processes in nuclear and cytoplasmic genomes with the different selection pressures, population structures and breeding systems.

A molecular phylogeny of Caraganeae (Leguminosae, Papilionoideae) reveals insights into new generic and infrageneric delimitations

Based on sequence data of nuclear ITS and plastid matK, trnL-F and psbA-trnH markers, the phylogeny of the subtribes Caraganinae and Chesneyinae in tribe Caraganeae was inferred. The results support the monophyly of each of the subtribes. Within subtribes Caraganinae, Calophaca and Halimodendron are herein transferred into Caragana to ensure its generic monophyly. The subtribe Chesneyinae is composed of four well-supported genera: Chesneya, Chesniella, Gueldenstaedtia and Tibetia. Based on phylogenetic, morphological, distributional and habitat type evidence, the genus Chesneya was divided into three monophyletic sections: Chesneya sect. Chesneya, Chesneya sect. Pulvinatae and Chesneya sect. Spinosae. Chesneya macrantha is herein transferred into Chesniella. Spongiocarpella is polyphyletic and its generic rank is not maintained. The position of Chesneya was incongruent in the nuclear ITS and the plastid trees. A paternal chloroplast capture event via introgression is hypothesized for the origin of Chesneya, which is postulated to have involved the common ancestor of Chesniella (♂) and that of the Gueldenstaedtia - Tibetia (GUT) clade (♀) as the parents.

Species delimitation in plants using the Qinghai-Tibet Plateau endemic Orinus (Poaceae: Tridentinae) as an example

BACKGROUND AND AIMS

Accurate identification of species is essential for the majority of biological studies. However, defining species objectively and consistently remains a challenge, especially for plants distributed in remote regions where there is often a lack of sufficient previous specimens. In this study, multiple approaches and lines of evidence were used to determine species boundaries for plants occurring in the Qinghai-Tibet Plateau, using the genus Orinus (Poaceae) as a model system for an integrative approach to delimiting species.

METHODS

A total of 786 individuals from 102 populations of six previously recognized species were collected for niche, morphological and genetic analyses. Three plastid DNA regions (matK, rbcL and trnH-psbA) and one nuclear DNA region [internal transcribed space (ITS)] were sequenced.

KEY RESULTS

Whereas six species had been previously recognized, statistical analyses based on character variation, molecular data and niche differentiation identified only two well-delimited clusters, together with a third possibly originating from relatively recent hybridization between, or historical introgression from, the other two.

CONCLUSIONS

Based on a principle of integrative species delimitation to reconcile different sources of data, the results provide compelling evidence that the six previously recognized species of the genus Orinus that were examined should be reduced to two, with new circumscriptions, and a third, identified in this study, should be described as a new species. This empirical study highlights the value of applying genetic differentiation, morphometric statistics and ecological niche modelling in an integrative approach to re-circumscribing species boundaries. The results produce relatively objective, operational and unbiased taxonomic classifications of plants occurring in remote regions.

A trnI_CAU triplication event in the complete chloroplast genome of Paris verticillata M.Bieb. (Melanthiaceae, Liliales)

The chloroplast is an essential plant organelle responsible for photosynthesis. Gene duplication, relocation, and loss in the chloroplast genome (cpDNA) are useful for exploring the evolution and phylogeny of plant species. In this study, the complete chloroplast genome of Paris verticillata was sequenced using the 454 sequencing system and Sanger sequencing method to trace the evolutionary pattern in the tribe Parideae of the family Melanthiaceae (Liliales). The circular double-stranded cpDNA of P. verticillata (157,379 bp) consists of two inverted repeat regions each of 28,373 bp, a large single copy of 82,726 bp, and a small single copy of 17,907 bp. Gene content and order are generally similar to the previously reported cpDNA sequences within the order Liliales. However, we found that trnI_CAU was triplicated in P. verticillata. In addition, cemA is suspected to be a pseudogene due to the presence of internal stop codons created by poly(A) insertion and single small CA repeats. Such changes were not found in previously examined cpDNAs of the Melanthiaceae or other families of the Liliales, suggesting that such features are unique to the tribe Parideae of Melanthiaceae. The characteristics of P. verticillata cpDNA will provide useful information for uncovering the evolution within Paris and for further research of plastid genome evolution and phylogenetic studies in Liliales.

Guidelines for the choice of sequences for molecular plant taxonomy

This chapter presents an overview of the major plant DNA sequences and molecular methods available for plant taxonomy. Guidelines are provided for the choice of sequences and methods to be used, based on the DNA compartment (nuclear, chloroplastic, mitochondrial), evolutionary mechanisms, and the level of taxonomic differentiation of the plants under survey.

Diverse microbiota identified in whole intact nest chambers of the red mason bee Osmia bicornis (Linnaeus 1758)

Microbial activity is known to have profound impact on bee ecology and physiology, both by beneficial and pathogenic effects. Most information about such associations is available for colony-building organisms, and especially the honey bee. There, active manipulations through worker bees result in a restricted diversity of microbes present within the colony environment. Microbial diversity in solitary bee nests remains unstudied, although their larvae face a very different situation compared with social bees by growing up in isolated compartments. Here, we assessed the microbiota present in nests and pre-adults of Osmia bicornis, the red mason bee, by culture-independent pyrosequencing. We found high bacterial diversity not comparable with honey bee colonies. We identified a variety of bacteria potentially with positive or negative interactions for bee larvae. However, most of the other diverse bacteria present in the nests seem to originate from environmental sources through incorporated nest building material and stored pollen. This diversity of microorganisms may cause severe larval mortality and require specific physiological or symbiotic adaptations against microbial threats. They may however also profit from such a diverse environment through gain of mutualistic partners. We conclude that further studies of microbiota interaction in solitary bees will improve the understanding of fitness components and populations dynamics.

Reprint of: using nuclear gene data for plant phylogenetics: progress and prospects

The paper reviews the current state of low and single copy nuclear markers that have been applied successfully in plant phylogenetics to date, and discusses case studies highlighting the potential of massively parallel high throughput or next-generation sequencing (NGS) approaches for molecular phylogenetic and evolutionary investigations. The current state, prospects and challenges of specific single- or low-copy plant nuclear markers as well as phylogenomic case studies are presented and evaluated.

Bioclimatic, ecological, and phenotypic intermediacy and high genetic admixture in a natural hybrid of octoploid strawberries

PREMISE OF THE STUDY

Hybrid zones provide "natural laboratories" for understanding the processes of selection, reinforcement, and speciation. We sought to gain insight into the degree of introgression and the extent of ecological-phenotypic intermediacy in the natural hybrid strawberry, Fragaria × ananassa subsp. cuneifolia. •

METHODS

We used whole-plastome sequencing to identify parental species-specific (Fragaria chiloensis and F. virginiana) chloroplast single-nucleotide polymorphisms and combined the use of these with nuclear microsatellite markers to genetically characterize the hybrid zone. We assessed the potential role of selection in the observed geographic patterns by bioclimatically characterizing the niche of the hybrid populations and phenotypically characterizing hybrid individuals of known genomic constitution. •

KEY RESULTS

Significant admixture and little overall maternal bias in chloroplast or nuclear genomes suggest a high degree of interfertility among the parental and hybrid species and point to a long history of backcrossing and genetic mixing in the hybrid zone. Even though hybrids were phenotypically intermediate to the parental species, there was a discernible fingerprint of the parental genotype within hybrid individuals. Thus, although the pattern of introgression observed suggests geographic limitations to gene flow, it may be reinforced by selection for specific parental traits in the bioclimatically intermediate habitat occupied by the hybrid. •

CONCLUSIONS

This work uncovered the genetic complexity underlying the hybrid zone of the wild relatives of the cultivated strawberry. It lays the foundation for experimental dissection of the causes of genomic introgression and nuclear-cytoplasmic disassociation, and for understanding other parts of Fragaria evolutionary history.

Sex-determining chromosomes and sexual dimorphism: insights from genetic mapping of sex expression in a natural hybrid Fragaria × ananassa subsp. cuneifolia

We studied the natural hybrid (Fragaria × ananassa subsp. cuneifolia) between two sexually dimorphic octoploid strawberry species (Fragaria virginiana and Fragaria chiloensis) to gain insight into the dynamics of sex chromosomes and the genesis of sexual dimorphism. Male sterility is dominant in both the parental species and thus will be inherited maternally, but the chromosome that houses the sex-determining region differs. Thus, we asked whether (1) the cytotypic composition of hybrid populations represents one or both maternal species, (2) the sex-determining chromosome of the hybrid reflects the location of male sterility within the maternal donor species and (3) crosses from the hybrid species show less sexual dimorphism than the parental species. We found that F. × ananassa subsp. cuneifolia populations consisted of both parental cytotypes but one predominated within each population. Genetic linkage mapping of two crosses showed dominance of male sterility similar to the parental species, however, the map location of male sterility reflected the maternal donor in one cross, but not the other. Moreover, female function mapped to a single region in the first cross, but to two regions in the second cross. Aside from components of female function (fruit set and seed set), other traits that have been found to be significantly sexually dimorphic in the pure species were either not dimorphic or were dimorphic in the opposite direction to the parental species. These results suggest that hybrids experience some disruption of dimorphism in secondary sexual traits, as well as novel location and number of quantitative trait locus (QTL) affecting sex function.

De novo transcriptome assembly and polymorphism detection in the flowering plant Silene vulgaris (Caryophyllaceae)

Members of the angiosperm genus Silene are widely used in studies of ecology and evolution, but available genomic and population genetic resources within Silene remain limited. Deep transcriptome (i.e. expressed sequence tag or EST) sequencing has proven to be a rapid and cost-effective means to characterize gene content and identify polymorphic markers in non-model organisms. In this study, we report the results of 454 GS-FLX Titanium sequencing of a polyA-selected and normalized cDNA library from Silene vulgaris. The library was generated from a single pool of transcripts, combining RNA from leaf, root and floral tissue from three genetically divergent European subpopulations of S. vulgaris. A single full-plate 454 run produced 959,520 reads totalling 363.6 Mb of sequence data with an average read length of 379.0 bp after quality trimming and removal of custom library adaptors. We assembled 832,251 (86.7%) of these reads into 40,964 contigs, which have a total length of 25.4 Mb and can be organized into 18,178 graph-based clusters or 'isogroups'. Assembled sequences were annotated based on homology to genes in multiple public databases. Analysis of sequence variants identified 13,432 putative single-nucleotide polymorphisms (SNPs) and 1320 simple sequence repeats (SSRs) that are candidates for microsatellite analysis. Estimates of nucleotide diversity from 1577 contigs were used to generate genome-wide distributions that revealed several outliers with high diversity. All of these resources are publicly available through NCBI and/or our website (http://silenegenomics.biology.virginia.edu) and should provide valuable genomic and population genetic tools for the Silene research community.

Speciation genes in plants

BACKGROUND

Analyses of speciation genes--genes that contribute to the cessation of gene flow between populations--can offer clues regarding the ecological settings, evolutionary forces and molecular mechanisms that drive the divergence of populations and species. This review discusses the identities and attributes of genes that contribute to reproductive isolation (RI) in plants, compares them with animal speciation genes and investigates what these genes can tell us about speciation.

SCOPE

Forty-one candidate speciation genes were identified in the plant literature. Of these, seven contributed to pre-pollination RI, one to post-pollination, prezygotic RI, eight to hybrid inviability, and 25 to hybrid sterility. Genes, gene families and genetic pathways that were frequently found to underlie the evolution of RI in different plant groups include the anthocyanin pathway and its regulators (pollinator isolation), S RNase-SI genes (unilateral incompatibility), disease resistance genes (hybrid necrosis), chimeric mitochondrial genes (cytoplasmic male sterility), and pentatricopeptide repeat family genes (cytoplasmic male sterility).

CONCLUSIONS

The most surprising conclusion from this review is that identities of genes underlying both prezygotic and postzygotic RI are often predictable in a broad sense from the phenotype of the reproductive barrier. Regulatory changes (both cis and trans) dominate the evolution of pre-pollination RI in plants, whereas a mix of regulatory mutations and changes in protein-coding genes underlie intrinsic postzygotic barriers. Also, loss-of-function mutations and copy number variation frequently contribute to RI. Although direct evidence of positive selection on speciation genes is surprisingly scarce in plants, analyses of gene family evolution, along with theoretical considerations, imply an important role for diversifying selection and genetic conflict in the evolution of RI. Unlike in animals, however, most candidate speciation genes in plants exhibit intraspecific polymorphism, consistent with an important role for stochastic forces and/or balancing selection in development of RI in plants.

Phylogenetic analysis of mitochondrial substitution rate variation in the angiosperm tribe Sileneae

BACKGROUND

Recent phylogenetic studies have revealed that the mitochondrial genome of the angiosperm Silene noctiflora (Caryophyllaceae) has experienced a massive mutation-driven acceleration in substitution rate, placing it among the fastest evolving eukaryotic genomes ever identified. To date, it appears that other species within Silene have maintained more typical substitution rates, suggesting that the acceleration in S. noctiflora is a recent and isolated evolutionary event. This assessment, however, is based on a very limited sampling of taxa within this diverse genus.

RESULTS

We analyzed the substitution rates in 4 mitochondrial genes (atp1, atp9, cox3 and nad9) across a broad sample of 74 species within Silene and related genera in the tribe Sileneae. We found that S. noctiflora shares its history of elevated mitochondrial substitution rate with the closely related species S. turkestanica. Another section of the genus (Conoimorpha) has experienced an acceleration of comparable magnitude. The phylogenetic data remain ambiguous as to whether the accelerations in these two clades represent independent evolutionary events or a single ancestral change. Rate variation among genes was equally dramatic. Most of the genus exhibited elevated rates for atp9 such that the average tree-wide substitution rate for this gene approached the values for the fastest evolving branches in the other three genes. In addition, some species exhibited major accelerations in atp1 and/or cox3 with no correlated change in other genes. Rates of non-synonymous substitution did not increase proportionally with synonymous rates but instead remained low and relatively invariant.

CONCLUSION

The patterns of phylogenetic divergence within Sileneae suggest enormous variability in plant mitochondrial mutation rates and reveal a complex interaction of gene and species effects. The variation in rates across genomic and phylogenetic scales raises questions about the mechanisms responsible for the evolution of mutation rates in plant mitochondrial genomes.

Effects of environmental heterogeneity on the distribution of sexes within and among populations in a gynodioecious species, Geranium maculatum

Populations containing both females and hermaphrodites (dimorphic) are generally found in drier sites than those with only hermaphrodites (monomorphic). The sex-differential plasticity hypothesis (SDP) suggests that this is caused by hermaphrodites reducing allocation to seeds in harsh environments, allowing female establishment. We proposed that a similar process could explain sex distribution within populations. We compared light availability and soil moisture between sites of three monomorphic and three dimorphic populations of Geranium maculatum and between microsites occupied by females and hermaphrodites within populations. We also correlated seed production in dimorphic populations with environmental measures. We found that dimorphic and monomorphic populations occurred in sites with similar soil moisture but within two dimorphic populations females occurred in drier microsites than hermaphrodites, as predicted by the SDP hypothesis. Contrary to the predictions, hermaphrodites' seed production was not influenced by the environment. Rather, females' seed production was correlated with environmental conditions in two populations, although the direction of the correlation differed between populations. Our results suggest that in this species, the SDP hypothesis does not explain sex distribution among or within populations. However, microsite environments may influence the distribution of sexes within a population and potentially aid in maintaining gynodioecy.

Detection of rare paternal chloroplast inheritance in controlled crosses of the endangered sunflower Helianthus verticillatus

A variety of questions in population and evolutionary biology are studied using chloroplast DNA (cpDNA). The presumed maternal inheritance in angiosperms allows for certain assumptions and calculations to be made when studying plant hybridization, phylogeography, molecular systematics and seed dispersal. Further, the placement of transgenes in the chloroplast to lessen the probability of 'escape' to weedy relatives has been proposed since such genes would not move through pollen. In many studies, however, strict maternal inheritance is assumed but not tested directly, and some studies may have sample sizes too small to be able to detect rare paternal leakage. Here, we study the inheritance of cpDNA simple sequence repeats in 323 offspring derived from greenhouse crosses of the rare sunflower Helianthus verticillatus Small. We found evidence for rare chloroplast paternal leakage and heteroplasmy in 1.86% of the offspring. We address the question of whether one can extrapolate the mode of chloroplast transmission within a genus by comparing our results to the findings of another sunflower species study. The findings of occasional paternal transmission of the chloroplast genome are discussed in the framework of using these markers in studies of population and evolutionary biology both in Helianthus and other angiosperms.