Utility of 17 chloroplast genes for inferring the phylogeny of the basal angiosperms

The chloroplast genome sequences of Ipomoea alba and I. obscura (Convolvulaceae): genome comparison and phylogenetic analysis

Ipomoea species have diverse uses as ornamentals, food, and medicine. However, their genomic information is limited; I. alba and I. obscura were sequenced and assembled. Their chloroplast genomes were 161,353 bp and 159,691 bp, respectively. Both genomes exhibited a quadripartite structure, consisting of a pair of inverted repeat (IR) regions, which are separated by the large single-copy (LSC) and small single-copy (SSC) regions. The overall GC content was 37.5% for both genomes. A total of 104 and 93 simple sequence repeats, 50 large repeats, and 30 and 22 short tandem repeats were identified in the two chloroplast genomes, respectively. G and T were more preferred than C and A at the third base position based on the Parity Rule 2 plot analysis, and the neutrality plot revealed correlation coefficients of 0.126 and 0.105, indicating the influence of natural selection in shaping the codon usage bias in most protein-coding genes (CDS). Genome comparative analyses using 31 selected Ipomoea taxa from Thailand showed that their chloroplast genomes are rather conserved, but the presence of expansion or contraction of the IR region was identified in some of these Ipomoea taxa. A total of five highly divergent regions were identified, including the CDS genes accD, ndhA, and ndhF, as well as the intergenic spacer regions psbI-atpA and rpl32-ccsA. Phylogenetic analysis based on both the complete chloroplast genome sequence and CDS datasets of 31 Ipomoea taxa showed that I. alba is resolved as a group member for series (ser.) Quamoclit, which contains seven other taxa, including I. hederacea, I. imperati, I. indica, I. nil, I. purpurea, I. quamoclit, and I. × sloteri, while I. obscura is grouped with I. tiliifolia, both of which are under ser. Obscura, and is closely related to I. biflora of ser. Pes-tigridis. Divergence time estimation using the complete chloroplast genome sequence dataset indicated that the mean age of the divergence for Ipomoeeae, Argyreiinae, and Astripomoeinae, was approximately 29.99 Mya, 19.81 Mya, and 13.40 Mya, respectively. The node indicating the divergence of I. alba from the other members of Ipomoea was around 10.06 Mya, and the split between I. obscura and I. tiliifolia is thought to have happened around 17.13 Mya. The split between the I. obscura accessions from Thailand and Taiwan is thought to have taken place around 0.86 Mya.

The first two whole mitochondrial genomes for the genus Dactylis species: assembly and comparative genomics analysis

BACKGROUND

Orchardgrass (Dactylis glomerata L.), a perennial forage, has the advantages of rich leaves, high yield, and good quality and is one of the most significant forage for grassland animal husbandry and ecological management in southwest China. Mitochondrial (mt) genome is one of the major genetic systems in plants. Studying the mt genome of the genus Dactylis could provide more genetic information in addition to the nuclear genome project of the genus.

RESULTS

In this study, we sequenced and assembled two mitochondrial genomes of Dactylis species of D. glomerata (597, 281 bp) and D. aschersoniana (613, 769 bp), based on a combination of PacBio and Illumina. The gene content in the mitochondrial genome of D. aschersoniana is almost identical to the mitochondrial genome of D. glomerata, which contains 22-23 protein-coding genes (PCGs), 8 ribosomal RNAs (rRNAs) and 30 transfer RNAs (tRNAs), while D. glomerata lacks the gene encoding the Ribosomal protein (rps1) and D. aschersoniana contains one pseudo gene (atp8). Twenty-three introns were found among eight of the 30 protein-coding genes, and introns of three genes (nad 1, nad2, and nad5) were trans-spliced in Dactylis aschersoniana. Further, our mitochondrial genome characteristics investigation of the genus Dactylis included codon usage, sequences repeats, RNA editing and selective pressure. The results showed that a large number of short repetitive sequences existed in the mitochondrial genome of D. aschersoniana, the size variation of two mitochondrial genomes is due largely to the presence of a large number of short repetitive sequences. We also identified 52-53 large fragments that were transferred from the chloroplast genome to the mitochondrial genome, and found that the similarity was more than 70%. ML and BI methods used in phylogenetic analysis revealed that the evolutionary status of the genus Dactylis.

CONCLUSIONS

Thus, this study reveals the significant rearrangements in the mt genomes of Pooideae species. The sequenced Dactylis mt genome can provide more genetic information and improve our evolutionary understanding of the mt genomes of gramineous plants.

Nuclear phylogenomics of angiosperms and insights into their relationships and evolution

Angiosperms (flowering plants) are by far the most diverse land plant group with over 300,000 species. The sudden appearance of diverse angiosperms in the fossil record was referred to by Darwin as the "abominable mystery," hence contributing to the heightened interest in angiosperm evolution. Angiosperms display wide ranges of morphological, physiological, and ecological characters, some of which have probably influenced their species richness. The evolutionary analyses of these characteristics help to address questions of angiosperm diversification and require well resolved phylogeny. Following the great successes of phylogenetic analyses using plastid sequences, dozens to thousands of nuclear genes from next-generation sequencing have been used in angiosperm phylogenomic analyses, providing well resolved phylogenies and new insights into the evolution of angiosperms. In this review we focus on recent nuclear phylogenomic analyses of large angiosperm clades, orders, families, and subdivisions of some families and provide a summarized Nuclear Phylogenetic Tree of Angiosperm Families. The newly established nuclear phylogenetic relationships are highlighted and compared with previous phylogenetic results. The sequenced genomes of Amborella, Nymphaea, Chloranthus, Ceratophyllum, and species of monocots, Magnoliids, and basal eudicots, have facilitated the phylogenomics of relationships among five major angiosperms clades. All but one of the 64 angiosperm orders were included in nuclear phylogenomics with well resolved relationships except the placements of several orders. Most families have been included with robust and highly supported placements, especially for relationships within several large and important orders and families. Additionally, we examine the divergence time estimation and biogeographic analyses of angiosperm on the basis of the nuclear phylogenomic frameworks and discuss the differences compared with previous analyses. Furthermore, we discuss the implications of nuclear phylogenomic analyses on ancestral reconstruction of morphological, physiological, and ecological characters of angiosperm groups, limitations of current nuclear phylogenomic studies, and the taxa that require future attention.

Mitochondrial phylogenomics of Cuscuta (Convolvulaceae) reveals a potentially functional horizontal gene transfer from the host

Horizontal gene transfers (HGTs) from host or other organisms have been reported in mitochondrial genomes of parasitic plants. Genes transferred in this fashion have usually been found non-functional. Several examples of HGT from the mitochondrial genome of parasitic Cuscuta (Convolvulaceae) to its hosts have been reported, but not vice versa. Here we used 31 protein-coding mitochondrial genes to infer the phylogeny of Cuscuta, and compared it with previous nuclear and plastid phylogenetic estimates. We also investigated the presence of HGTs within these lineages. Unlike in plastid genomes, we did not find extensive gene loss in their mitochondrial counterparts. Our results reveal the first example of organellar HGT from host to Cuscuta. Mitochondrial atp1 genes of South African subgenus Pachystigma were inferred to be transferred from Lamiales, with high support. Moreover, the horizontally transferred atp1 gene has functionally replaced the native, vertically transmitted copy, has an intact open reading frame, and is under strong purifying selection, all of which suggests that this xenolog remains functional. The mitochondrial phylogeny of Cuscuta is generally consistent with previous plastid and nuclear phylogenies, except for the misplacement of Pachystigma when atp1 is included. This incongruence may be caused by the HGT mentioned above. No example of HGT was found within mitochondrial genes of three other, independently evolved parasitic lineages we sampled: Cassytha/Laurales, Krameria/Zygophyllales, and Lennooideae/Boraginales.

Molecular Diversity and Phylogeny Reconstruction of Genus Colobanthus (Caryophyllaceae) Based on Mitochondrial Gene Sequences

Mitochondrial genomes have become an interesting object of evolutionary and systematic study both for animals and plants, including angiosperms. Although the framework of the angiosperm phylogeny was built on the information derived from chloroplast and nuclear genes, mitochondrial sequences also revealed their usefulness in solving the phylogenetic issues at different levels of plant systematics. Here, we report for the first time the complete sequences of 26 protein-coding genes of eight Colobanthus species (Caryophyllaceae). Of these, 23 of them represented core mitochondrial genes, which are directly associated with the primary function of that organelle, and the remaining three genes represented a facultative set of mitochondrial genes. Comparative analysis of the identified genes revealed a generally high degree of sequence conservation. The Ka/Ks ratio was <1 for most of the genes, which indicated purifying selection. Only for rps12 was Ka/Ks > 1 in all studied species, suggesting positive selection. We identified 146−165 potential RNA editing sites in genes of the studied species, which is lower than in most angiosperms. The reconstructed phylogeny based on mitochondrial genes was consistent with the taxonomic position of the studied species, showing the separate character of the family Caryophyllaceae and close relationships between all studied Colobanthus species, with C. lycopodioides sharing less similarity.

Plastid phylogenomic analyses of the Selaginella sanguinolenta group (Selaginellaceae) reveal conflict signatures resulting from sequence types, outlier genes, and pervasive RNA editing

Different from the generally conserved plastomes (plastid genomes) of most land plants, the Selaginellaceae plastomes exhibit dynamic structure, high GC content and high substitution rates. Previous plastome analyses identified strong conflict on several clades in Selaginella, however the factors causing the conflictions and the impact on the phylogenetic inference have not been sufficiently investigated. Here, we dissect the distribution of phylogenetic signals and conflicts in Selaginella sanguinolenta group, the plastome of which is DR (direct repeats) structure and with genome-wide RNA editing. We analyzed the data sets including 22 plastomes representing all species of the S. sanguinolenta group, covering the entire geographical distribution from the Himalayas to Siberia and the Russian Far East regions. We recovered four different topologies by applying multispecies coalescent (ASTRAL) and concatenation methods (IQ-TREE and RAxML) on four data sets of PC (protein-coding genes), NC (non-coding sequences), PCN (the concatenated PC and NC), and RC (predicted RNA editing sites "C" were corrected by "T"), respectively. Six monophyletic clades, S. nummularifolia clade, S. rossii clade, S. sajanensis clade, S. sanguinolenta I clade, S. sanguinolenta II clade, and S. sanguinolenta III clade, were consistently resolved and supported by the characteristics of GC content, RNA editing frequency, and gene content. However, the relationships among these clades varied across the four topologies. To explore the underlying causes of the uncertainty, we compared the phylogenetic signals of the four topologies. We identified that the sequence types (coding versus non-coding), outlier genes (genes with extremely high |ΔGLS| values), and C-to-U RNA editing frequency in the protein-coding genes were responsible for the unstable phylogenomic relationship. We further revealed a significant positive correlation between the |ΔGLS| values and the variation coefficient of the RNA editing number. Our results demonstrated that the coalescent method performed better than the concatenation method in overcoming the problems caused by outlier genes and extreme RNA editing events. Our study particularly focused on the importance of exploring the plastid phylogenomic conflicts and suggested conducting concatenated analyses cautiously when adopting organelle genome data.

Plastid phylogenomics of the Gynoxoid group (Senecioneae, Asteraceae) highlights the importance of motif-based sequence alignment amid low genetic distances

PREMISE

The genus Gynoxys and relatives form a species-rich lineage of Andean shrubs and trees with low genetic distances within the sunflower subtribe Tussilaginineae. Previous molecular phylogenetic investigations of the Tussilaginineae have included few, if any, representatives of this Gynoxoid group or reconstructed ambiguous patterns of relationships for it.

METHODS

We sequenced complete plastid genomes of 21 species of the Gynoxoid group and related Tussilaginineae and conducted detailed comparisons of the phylogenetic relationships supported by the gene, intron, and intergenic spacer partitions of these genomes. We also evaluated the impact of manual, motif-based adjustments of automatic DNA sequence alignments on phylogenetic tree inference.

RESULTS

Our results indicate that the inclusion of all plastid genome partitions is needed to infer well-supported phylogenetic trees of the Gynoxoid group. Whole plastome-based tree inference suggests that the genera Gynoxys and Nordenstamia are polyphyletic and form the core clade of the Gynoxoid group. This clade is sister to a clade of Aequatorium and Paragynoxys and also includes some but not all representatives of Paracalia.

CONCLUSIONS

The concatenation and combined analysis of all plastid genome partitions and the construction of manually-curated, motif-based DNA sequence alignments are found to be instrumental in the recovery of well-supported relationships of the Gynoxoid group. We demonstrate that the correct assessment of homology in genome-level plastid sequence data sets is crucial for subsequent phylogeny reconstruction and that the manual post-processing of multiple sequence alignments improves the reliability of such reconstructions amid low genetic distances between taxa.

Plastid phylogenomic insights into relationships of all flowering plant families

BACKGROUND

Flowering plants (angiosperms) are dominant components of global terrestrial ecosystems, but phylogenetic relationships at the familial level and above remain only partially resolved, greatly impeding our full understanding of their evolution and early diversification. The plastome, typically mapped as a circular genome, has been the most important molecular data source for plant phylogeny reconstruction for decades.

RESULTS

Here, we assembled by far the largest plastid dataset of angiosperms, composed of 80 genes from 4792 plastomes of 4660 species in 2024 genera representing all currently recognized families. Our phylogenetic tree (PPA II) is essentially congruent with those of previous plastid phylogenomic analyses but generally provides greater clade support. In the PPA II tree, 75% of nodes at or above the ordinal level and 78% at or above the familial level were resolved with high bootstrap support (BP ≥ 90). We obtained strong support for many interordinal and interfamilial relationships that were poorly resolved previously within the core eudicots, such as Dilleniales, Saxifragales, and Vitales being resolved as successive sisters to the remaining rosids, and Santalales, Berberidopsidales, and Caryophyllales as successive sisters to the asterids. However, the placement of magnoliids, although resolved as sister to all other Mesangiospermae, is not well supported and disagrees with topologies inferred from nuclear data. Relationships among the five major clades of Mesangiospermae remain intractable despite increased sampling, probably due to an ancient rapid radiation.

CONCLUSIONS

We provide the most comprehensive dataset of plastomes to date and a well-resolved phylogenetic tree, which together provide a strong foundation for future evolutionary studies of flowering plants.

Discordant Phylogenomic Placement of Hydnoraceae and Lactoridaceae Within Piperales Using Data From All Three Genomes

Phylogenetic relationships within the magnoliid order Piperales have been studied extensively, yet the relationships of the monotypic family Lactoridaceae and the holoparasitic Hydnoraceae to the remainder of the order remain a matter of debate. Since the first confident molecular phylogenetic placement of Hydnoraceae among Piperales, different studies have recovered various contradictory topologies. Most phylogenetic hypotheses were inferred using only a few loci and have had incomplete taxon sampling at the genus level. Based on these results and an online survey of taxonomic opinion, the Angiosperm Phylogeny Group lumped both Hydnoraceae and Lactoridaceae in Aristolochiaceae; however, the latter family continues to have unclear relationships to the aforementioned taxa. Here we present extensive phylogenomic tree reconstructions based on up to 137 loci from all three subcellular genomes for all genera of Piperales. We infer relationships based on a variety of phylogenetic methods, explore instances of phylogenomic discordance between the subcellular genomes, and test alternative topologies. Consistent with these phylogenomic results and a consideration of the principles of phylogenetic classification, we propose to exclude Hydnoraceae and Lactoridaceae from the broad circumscription of Aristolochiaceae, and instead favor recognition of four monophyletic and morphologically well circumscribed families in the perianth-bearing Piperales: Aristolochiaceae, Asaraceae, Hydnoraceae, and Lactoridaceae, with a total of six families in the order.

siRNA biogenesis and advances in topically applied dsRNA for controlling virus infections in tomato plants

A non-transgenic approach based on RNA interference was employed to induce protection against tomato mosaic virus (ToMV) infection in tomato plants. dsRNA molecules targeting the cp gene of ToMV were topically applied on plants prior to virus inoculation. Protection was dose-dependent and sequence-specific. While no protection was achieved when 0-16 µg dsRNA were used, maximum rates of resistance (60 and 63%) were observed in doses of 200 and 400 µg/plant, respectively. Similar rates were also obtained against potato virus Y when targeting its cp gene. The protection was quickly activated upon dsRNA application and lasted for up to 4 days. In contrast, no detectable antiviral response was triggered by the dsRNA from a begomovirus genome, suggesting the method is not effective against phloem-limited DNA viruses. Deep sequencing was performed to analyze the biogenesis of siRNA populations. Although long-dsRNA remained in the treated leaves for at least 10 days, its systemic movement was not observed. Conversely, dsRNA-derived siRNA populations (mainly 21- and 22-nt) were detected in non-treated leaves, which indicates endogenous processing and transport through the plant. Altogether, this study provides critical information for the development of novel tools against plant viruses; strengths and limitations inherent to the systems are discussed.

Comprehensive study of the genes involved in chlorophyll synthesis and degradation pathways in some monocot and dicot plant species

Chlorophyll (Chl) biosynthesis is one of the most important cellular processes essential for plant photosynthesis. Chl degradation pathway is also important catabolic process occurs during leaf senescence, fruit ripening and under biotic or abiotic stress conditions. Here we have systematically investigated the molecular evolution, gene structure, compositional analysis along with ENc plot, correspondence analysis and codon usage bias of the proteins and encoded genes involved in Chl metabolism from monocots and dicots. The gene and species specific phylogenetic trees using amino acid sequences showed clear clustering formation of the selected species based on monocots and dicots but not supported by 18S rRNA. Nucleotide composition of the encoding genes showed that average GC%, GC1%, GC2% and GC3% were higher in monocots. RSCU analysis depicts that genes from monocots for both pathways and genes for synthesis pathway from dicots only biased to G/C-ending synonymous codons but in degradation pathway most optimal codons (except UUG) in dicots biased to A/U-ending synonymous codons. We found strong evidence of episodic diversifying selection at several amino acid sites in all genes investigated. Conserved domain and gene structures were observed for the genes with varying lengths of introns and exons, involved in Chl metabolism along with some intronless genes within synthesis pathway. ENc and correspondence analyses suggested the mutational or selection constraint on the genes to shape the codon usage. These comprehensive studies may be helpful in further research in molecular phylogenetics and genomics and to better understand the evolutionary dynamics of Chl metabolic pathway.Communicated by Ramaswamy H. Sarma.

Implications of plastome evolution in the true lilies (monocot order Liliales)

The families of the monocot order Liliales exhibit highly contrasting characteristic of photosynthetic and mycoheterotrophic life histories. Although previous phylogenetic and morphological studies of Liliales have been conducted, they have not examined molecular evolution associated with this contrasting phenomenon. Here, we conduct the first comparative plastome study of all ten families of Liliales using 29 newly sequenced plastid genomes analyzed together with previously published data. We also present a phylogenetic analysis for Liliales of 78 plastid genes combined with 22 genes from all three genomes (nuclear 18S rDNA and phyC; 17 plastid genes; and mitochondrial matR, atpA, and cob). Within the newly generated phylogenetic tree of Liliales, we evaluate the ancestral state changes of selected morphological traits in the order. There are no significant differences in plastid genome features among species that show divergent characteristics correlated with family circumscriptions. However, the results clearly differentiate between photosynthetic and mycoheterotrophic taxa of Liliales in terms of genome structure, and gene content and order. The newly sequenced plastid genomes and combined three-genome data revealed Smilacaceae as sister to Liliaceae instead of Philesiaceae and Ripogonaceae. Additionally, we propose a revised familial classification system of Liliales that consists of nine families, considering Ripogonaceae a synonym of Philesiaceae. The ancestral state reconstruction indicated synapomorphies for each family of Liliales, except Liliaceae, Melanthiaceae and Colchicaceae. A taxonomic key for all nine families of Liliales is also provided.

Expanded inverted repeat region with large scale inversion in the first complete plastid genome sequence of Plantago ovata

Plantago ovata (Plantaginaceae) is an economically and medicinally important species, however, least is known about its genomics and evolution. Here, we report the first complete plastome genome of P. ovata and comparison with previously published genomes of related species from Plantaginaceae. The results revealed that P. ovata plastome size was 162,116 bp and that it had typical quadripartite structure containing a large single copy region of 82,084 bp and small single copy region of 5,272 bp. The genome has a markedly higher inverted repeat (IR) size of 37.4 kb, suggesting large-scale inversion of 13.8 kb within the expanded IR regions. In addition, the P. ovata plastome contains 149 different genes, including 43 tRNA, 8 rRNA, and 98 protein-coding genes. The analysis revealed 139 microsatellites, of which 71 were in the non-coding regions. Approximately 32 forward, 34 tandem, and 17 palindromic repeats were detected. The complete genome sequences, 72 shared genes, matK gene, and rbcL gene from related species generated the same phylogenetic signals, and phylogenetic analysis revealed that P. ovata formed a single clade with P. maritima and P. media. The divergence time estimation as employed in BEAST revealed that P. ovata diverged from P. maritima and P. media about 11.0 million years ago (Mya; 95% highest posterior density, 10.06-12.25 Mya). In conclusion, P. ovata had significant variation in the IR region, suggesting a more stable P. ovata plastome genome than that of other Plantaginaceae species.

Mycoheterotrophic Epirixanthes (Polygalaceae) has a typical angiosperm mitogenome but unorthodox plastid genomes

BACKGROUND AND AIMS

Fully mycoheterotrophic plants derive carbon and other nutrients from root-associated fungi and have lost the ability to photosynthesize. While mycoheterotroph plastomes are often degraded compared with green plants, the effect of this unusual symbiosis on mitochondrial genome evolution is unknown. By providing the first complete organelle genome data from Polygalaceae, one of only three eudicot families that developed mycoheterotrophy, we explore how both organellar genomes evolved after loss of photosynthesis.

METHODS

We sequenced and assembled four complete plastid genomes and a mitochondrial genome from species of Polygalaceae, focusing on non-photosynthetic Epirixanthes. We compared these genomes with those of other mycoheterotroph and parasitic plant lineages, and assessed whether organelle genes in Epirixanthes experienced relaxed or intensified selection compared with autotrophic relatives.

KEY RESULTS

Plastomes of two species of Epirixanthes have become substantially degraded compared with that of autotrophic Polygala. Although the lack of photosynthesis is presumably homologous in the genus, the surveyed Epirixanthes species have marked differences in terms of plastome size, structural rearrangements, gene content and substitution rates. Remarkably, both apparently replaced a canonical plastid inverted repeat with large directly repeated sequences. The mitogenome of E. elongata incorporated a considerable number of fossilized plastid genes, by intracellular transfer from an ancestor with a less degraded plastome. Both plastid and mitochondrial genes in E. elongata have increased substitution rates, but the plastid genes of E. pallida do not. Despite this, both species have similar selection patterns operating on plastid housekeeping genes.

CONCLUSIONS

Plastome evolution largely fits with patterns of gene degradation seen in other heterotrophic plants, but includes highly unusual directly duplicated regions. The causes of rate elevation in the sequenced Epirixanthes mitogenome and of rate differences in plastomes of related mycoheterotrophic species are not currently understood.

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.

Genomic technologies for Hevea breeding

The commercial production of high quality natural rubber (NR) solely depends on Hevea brasiliensis Muell. Arg, (Para rubber tree) and accounts for >98% of total production worldwide. NR with its unique properties is an essential commodity for the automobile industry and its synthetic counterparts are in no way substitute to it. The rubber tree genome is very complex and plays an important role in delivering the unique properties of Hevea. But a lack of knowledge on the molecular mechanisms of rubber biosynthesis, disease resistance, etc., in elite clones of rubber still persists. Marker-assisted selection and transgenic techniques were proved to be advantageous in improving the breeding efficiency for latex yield, disease resistance, etc. The suppression subtractive hybridization (SSH), in the form of subtracted cDNA libraries and microarrays, can assist in searching the functions of expressed genes (candidate gene approach). Expressed sequence tags (ESTs) related to various metabolic aspects are well utilized to create EST banks that broadly represent the genes expressed in one tissue, such as latex cells, that assists in the study of gene function and regulation. Transcriptome analysis and gene mapping have been accomplished in Hevea at various stages. However, a selection criterion to delineate high yielding genotypes at the juvenile stage has not been accomplished so far. This is the main pit fall for rubber breeding apart from stock-scion interactions leading to yield differences among a clonally multiplied population. At least four draft genome sequences have been published on Hevea rubber, and all give different genome size and contig lengths-a comprehensive and acceptable genomic map remains unfulfilled. The progress made in molecular markers, latex biosynthesis genes, transcriptome analysis, chloroplast and mitochondrial DNA diversity, paternity identification through Breeding without Breeding (BwB), stimulated latex production and its molecular intricacies, molecular biology of tapping panel dryness, genomics for changed climates and genome mapping are discussed in this review. These information can be utilized to improvise the molecular breeding programs of Hevea in future.

Mitochondrial genome evolution in parasitic plants

Background

Parasitic plants rely on their host to cover their nutritional requirements either for their entire life or a smaller part of it. Depending on the level of parasitism, a proportional reduction on the plastid genome has been found. However, knowledge on gene loss and evolution of the mitogenome of parasitic plants is only available for four hemiparasitic Viscum species (Viscaceae), which lack many of the mitochondrial genes, while the remaining genes exhibit very fast molecular evolution rates. In this study, we include another genus, Phoradendron, from the Viscaceae, as well as 10 other hemiparasitic or holoparasitic taxa from across the phylogeny of the angiosperms to investigate how fast molecular evolution works on their mitogenomes, and the extent of gene loss.

Results

Our observations from Viscum were replicated in Phoradendron liga, whereas the remaining parasitic plants in the study have a complete set of the core mitochondrial genes and exhibit moderate or only slightly raised substitution rates compared to most autotrophic taxa, without any statistically significant difference between the different groups (autotrophs, hemiparasites and holoparasites). Additionally, further evidence is provided for the placement of Balanophoraceae within the order Santalales, while the exact placement of Cynomoriaceae still remains elusive.

Conclusions

We examine the mitochondrial gene content of 11 hemiparasitic and holoparasitic plants and confirm previous observations in Viscaceae. We show that the remaining parasitic plants do not have significantly higher substitution rates than autotrophic plants in their mitochondrial genes. We provide further evidence for the placement of Balanophoraceae in the Santalales.

Electronic supplementary material

The online version of this article (10.1186/s12862-019-1401-8) contains supplementary material, which is available to authorized users.

The Complete Chloroplast Genomes of Echinacanthus Species (Acanthaceae): Phylogenetic Relationships, Adaptive Evolution, and Screening of Molecular Markers

Among the four species of Echinacanthus (Acanthaceae), one distributed in the West Himalayan region and three restricted to the Sino-Vietnamese karst region. Because of its ecological significance, molecular markers are necessary for proper assessment of its genetic diversity and phylogenetic relationships. Herein, the complete chloroplast genomes of four Echinacanthus species were determined for the first time. The results indicated that all the chloroplast genomes were mapped as a circular structure and each genomes included 113 unique genes, of which 80 were protein-coding, 29 were tRNAs, and 4 were rRNAs. However, the four cp genomes ranged from 151,333 to 152,672 bp in length. Comparison of the four cp genomes showed that the divergence level was greater between geographic groups. We also analyzed IR expansion or contraction in the four cp genomes and the fifth type of the large single copy/inverted repeat region in Lamiales was suggested. Furthermore, based on the analyses of comparison and nucleotide variability, six most divergent sequences (rrn16, ycf1, ndhA, rps16-trnQ-UUG, trnS-GCU-trnG-UCC, and psaA-ycf3) were identified. A total of 37-45 simple sequence repeats were discovered in the four species and 22 SSRs were identified as candidate effective molecular markers for detecting interspecies polymorphisms. These SSRs and hotspot regions could be used as potential molecular markers for future study. Phylogenetic analysis based on Bayesian and parsimony methods did not support the monophyly of Echinacanthus. The phylogenetic relationships among the four species were clearly resolved and the results supported the recognition of the Sino-Vietnamese Echinacanthus species as a new genus. Based on the protein sequence evolution analysis, 12 genes (rpl14, rpl16, rps4, rps15, rps18, rps19, psbK, psbN, ndhC, ndhJ, rpoB, and infA) were detected under positive selection in branch of Sino-Vietnamese Echinacanthus species. These genes will lead to understanding the adaptation of Echinacanthus species to karst environment. The study will help to resolve the phylogenetic relationship and understand the adaptive evolution of Echinacanthus. It will also provide genomic resources and potential markers suitable for future species identification and speciation studies of the genus.

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.

Key questions and challenges in angiosperm macroevolution

Contents Summary 1170 I. Introduction 1170 II. Six key questions 1172 III. Three key challenges 1177 IV. Conclusions 1181 Acknowledgements 1182 References 1183 SUMMARY: The origin and rapid diversification of angiosperms (flowering plants) represent one of the most intriguing topics in evolutionary biology. Despite considerable progress made in complementary fields over the last two decades (paleobotany, phylogenetics, ecology, evo-devo, genomics), many important questions remain. For instance, what has been the impact of mass extinctions on angiosperm diversification? Are the angiosperms an adaptive radiation? Has morphological evolution in angiosperms been gradual or pulsed? We propose that the recent and ongoing revolution in macroevolutionary methods provides an unprecedented opportunity to explore long-standing questions that probably hold important clues to understand present-day biodiversity. We present six key questions that explore the origin and diversification of angiosperms. We also identify three key challenges to address these questions: (1) the development of new integrative models that include diversification, multiple intrinsic and environmental traits, biogeography and the fossil record all at once, whilst accounting for sampling bias and heterogeneity of macroevolutionary processes through time and among lineages; (2) the need for large and standardized synthetic databases of morphological variation; and (3) continuous effort on sampling the fossil record, but with a revolution in current paleobotanical practice.

Molecular dating of phylogenetic divergence between Urochloa species based on complete chloroplast genomes

BACKGROUND

Forage species of Urochloa are planted in millions of hectares of tropical and subtropical pastures in South America. Most of the planted area is covered with four species (U. ruziziensis, U. brizantha, U. decumbens and U. humidicola). Breeding programs rely on interspecific hybridizations to increase genetic diversity and introgress traits of agronomic importance. Knowledge of phylogenetic relationships is important to optimize compatible hybridizations in Urochloa, where phylogeny has been subject of some controversy. We used next-generation sequencing to assemble the chloroplast genomes of four Urochloa species to investigate their phylogenetic relationships, compute their times of divergence and identify chloroplast DNA markers (microsatellites, SNPs and InDels).

RESULTS

Whole plastid genome sizes were 138,765 bp in U. ruziziensis, 138,945 bp in U. decumbens, 138,946 bp in U. brizantha and 138,976 bp in U. humidicola. Each Urochloa chloroplast genome contained 130 predicted coding regions and structural features that are typical of Panicoid grasses. U. brizantha and U. decumbens chloroplast sequences are highly similar and show reduced SNP, InDel and SSR polymorphism as compared to U. ruziziensis and U. humidicola. Most of the structural and sequence polymorphisms were located in intergenic regions, and reflected phylogenetic distances between species. Divergence of U. humidicola from a common ancestor with the three other Urochloa species was estimated at 9.46 mya. U. ruziziensis, U. decumbens, and U. brizantha formed a clade where the U. ruziziensis lineage would have diverged by 5.67 mya, followed by a recent divergence event between U. decumbens and U. brizantha around 1.6 mya.

CONCLUSION

Low-coverage Illumina sequencing allowed the successful sequence analysis of plastid genomes in four species of Urochloa used as forages in the tropics. Pairwise sequence comparisons detected multiple microsatellite, SNP and InDel sites prone to be used as molecular markers in genetic analysis of Urochloa. Our results placed the origin of U. humidicola and U. ruziziensis divergence in the Miocene-Pliocene boundary, and the split between U. brizantha and U. decumbens in the Pleistocene.

The phylogeny of Heliconia (Heliconiaceae) and the evolution of floral presentation

Heliconia (Heliconiaceae, order Zingiberales) is among the showiest plants of the Neotropical rainforest and represent a spectacular co-evolutionary radiation with hummingbirds. Despite the attractiveness and ecological importance of many Heliconia, the genus has been the subject of limited molecular phylogenetic studies. We sample seven markers from the plastid and nuclear genomes for 202 samples of Heliconia. This represents ca. 75% of accepted species and includes coverage of all taxonomic subgenera and sections. We date this phylogeny using fossils associated with other families in the Zingiberales; in particular we review and evaluate the Eocene fossil Ensete oregonense. We use this dated phylogenetic framework to evaluate the evolution of two components of flower orientation that are hypothesized to be important for modulating pollinator discrimination and pollen placement: resupination and erect versus pendant inflorescence habit. Our phylogenetic results suggest that the monophyletic Melanesian subgenus Heliconiopsis and a small clade of Ecuadorian species are together the sister group to the rest of Heliconia. Extant diversity of Heliconia originated in the Late Eocene (39Ma) with rapid diversification through the Early Miocene, making it the oldest known clade of hummingbird-pollinated plants. Most described subgenera and sections are not monophyletic, though closely related groups of species, often defined by shared geography, mirror earlier morphological cladistic analyses. Evaluation of changes in resupination and inflorescence habit suggests that these characters are more homoplasious than expected, and this largely explains the non-monophyly of previously circumscribed subgenera, which were based on these characters. We also find strong evidence for the correlated evolution of resupination and inflorescence habit. The correlated model suggests that the most recent common ancestor of all extant Heliconia had resupinate flowers and erect inflorescences. Finally, we note our nearly complete species sampling and dated phylogeny allow for an assessment of taxonomic history in terms of phylogenetic diversity. We find approximately half of the currently recognized species, corresponding to half of the phylogenetic diversity, have been described since 1975, highlighting the continued importance of basic taxonomic research and conservation initiatives to preserve both described and undiscovered species of Heliconia.

Complete plastid genome sequence of Primula sinensis (Primulaceae): structure comparison, sequence variation and evidence for accD transfer to nucleus

Species-rich genus Primula L. is a typical plant group with which to understand genetic variance between species in different levels of relationships. Chloroplast genome sequences are used to be the information resource for quantifying this difference and reconstructing evolutionary history. In this study, we reported the complete chloroplast genome sequence of Primula sinensis and compared it with other related species. This genome of chloroplast showed a typical circular quadripartite structure with 150,859 bp in sequence length consisting of 37.2% GC base. Two inverted repeated regions (25,535 bp) were separated by a large single-copy region (82,064 bp) and a small single-copy region (17,725 bp). The genome consists of 112 genes, including 78 protein-coding genes, 30 tRNA genes and four rRNA genes. Among them, seven coding genes, seven tRNA genes and four rRNA genes have two copies due to their locations in the IR regions. The accD and infA genes lacking intact open reading frames (ORF) were identified as pseudogenes. SSR and sequence variation analyses were also performed on the plastome of Primula sinensis, comparing with another available plastome of P. poissonii. The four most variable regions, rpl36–rps8, rps16–trnQ, trnH–psbA and ndhC–trnV, were identified. Phylogenetic relationship estimates using three sub-datasets extracted from a matrix of 57 protein-coding gene sequences showed the identical result that was consistent with previous studies. A transcript found from P. sinensis transcriptome showed a high similarity to plastid accD functional region and was identified as a putative plastid transit peptide at the N-terminal region. The result strongly suggested that plastid accD has been functionally transferred to the nucleus in P. sinensis.

A few-gene plastid phylogenetic framework for mycoheterotrophic monocots

PREMISE OF THE STUDY

Few-gene studies with broad taxon sampling have provided major insights into phylogeny and underpin plant classification. However, they have typically excluded heterotrophic plants because of loss, pseudogenization, or rapid evolution of plastid genes. Here we performed a phylogenetic survey of three commonly retained plastid genes to assess their utility in placing mycoheterotrophs.

METHODS

We surveyed accD, clpP, and matK for 34 taxa from seven monocot families that include full mycoheterotrophs and a broad sampling of photosynthetic taxa. After screening for weak contaminants, we conducted phylogenetic analyses and characterized among-lineage rate variation.

KEY RESULTS

Likelihood analyses strongly supported local placements of fully mycoheterotrophic taxa for Corsiaceae, Iridaceae, Orchidaceae, and Petrosaviaceae, in positions consistent with other studies. Depression of likelihood bootstrap support values near mycoheterotrophic clades was alleviated when each mycoheterotrophic family was considered separately. Triuridaceae (Sciaphila) monophyly was recovered in a partitioned likelihood analysis, and the family then placed as sister to Cyclanthaceae-Pandanaceae. Burmanniaceae placed in Dioscoreales with weak to strong support depending on analysis details, and we inferred a plastid-based phylogeny for the family. Thismiaceae species may retain a plastid genome, based on accD retention. The inferred position of Thismiaceae is unstable, but was close to Taccaceae (Dioscoreales) in some analyses.

CONCLUSIONS

Long branches/elevated substitution rates, missing genes, and occasional contaminants are challenges for plastid-based phylogenetic inference with full mycoheterotrophs. However, most mycoheterotrophs can be readily integrated into the broad picture of plant phylogeny using several plastid genes and broad taxonomic sampling.

A Phylogenomic Approach Based on PCR Target Enrichment and High Throughput Sequencing: Resolving the Diversity within the South American Species of Bartsia L. (Orobanchaceae)

Advances in high-throughput sequencing (HTS) have allowed researchers to obtain large amounts of biological sequence information at speeds and costs unimaginable only a decade ago. Phylogenetics, and the study of evolution in general, is quickly migrating towards using HTS to generate larger and more complex molecular datasets. In this paper, we present a method that utilizes microfluidic PCR and HTS to generate large amounts of sequence data suitable for phylogenetic analyses. The approach uses the Fluidigm Access Array System (Fluidigm, San Francisco, CA, USA) and two sets of PCR primers to simultaneously amplify 48 target regions across 48 samples, incorporating sample-specific barcodes and HTS adapters (2,304 unique amplicons per Access Array). The final product is a pooled set of amplicons ready to be sequenced, and thus, there is no need to construct separate, costly genomic libraries for each sample. Further, we present a bioinformatics pipeline to process the raw HTS reads to either generate consensus sequences (with or without ambiguities) for every locus in every sample or—more importantly—recover the separate alleles from heterozygous target regions in each sample. This is important because it adds allelic information that is well suited for coalescent-based phylogenetic analyses that are becoming very common in conservation and evolutionary biology. To test our approach and bioinformatics pipeline, we sequenced 576 samples across 96 target regions belonging to the South American clade of the genus Bartsia L. in the plant family Orobanchaceae. After sequencing cleanup and alignment, the experiment resulted in ~25,300bp across 486 samples for a set of 48 primer pairs targeting the plastome, and ~13,500bp for 363 samples for a set of primers targeting regions in the nuclear genome. Finally, we constructed a combined concatenated matrix from all 96 primer combinations, resulting in a combined aligned length of ~40,500bp for 349 samples.

Five major shifts of diversification through the long evolutionary history of Magnoliidae (angiosperms)

Background

With 10,000 species, Magnoliidae are the largest clade of flowering plants outside monocots and eudicots. Despite an ancient and rich fossil history, the tempo and mode of diversification of Magnoliidae remain poorly known. Using a molecular data set of 12 markers and 220 species (representing >75% of genera in Magnoliidae) and six robust, internal fossil age constraints, we estimate divergence times and significant shifts of diversification across the clade. In addition, we test the sensitivity of magnoliid divergence times to the choice of relaxed clock model and various maximum age constraints for the angiosperms.

Results

Compared with previous work, our study tends to push back in time the age of the crown node of Magnoliidae (178.78-126.82 million years, Myr), and of the four orders, Canellales (143.18-125.90 Myr), Piperales (158.11-88.15 Myr), Laurales (165.62-112.05 Myr), and Magnoliales (164.09-114.75 Myr). Although families vary in crown ages, Magnoliidae appear to have diversified into most extant families by the end of the Cretaceous. The strongly imbalanced distribution of extant diversity within Magnoliidae appears to be best explained by models of diversification with 6 to 13 shifts in net diversification rates. Significant increases are inferred within Piperaceae and Annonaceae, while the low species richness of Calycanthaceae, Degeneriaceae, and Himantandraceae appears to be the result of decreases in both speciation and extinction rates.

Conclusions

This study provides a new time scale for the evolutionary history of an important, but underexplored, part of the tree of angiosperms. The ages of the main clades of Magnoliidae (above the family level) are older than previously thought, and in several lineages, there were significant increases and decreases in net diversification rates. This study is a new robust framework for future investigations of trait evolution and of factors influencing diversification in this group as well as angiosperms as a whole.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0320-6) contains supplementary material, which is available to authorized users.

Chloroplast DNA sequence utility for the lowest phylogenetic and phylogeographic inferences in angiosperms: the tortoise and the hare IV

PREMISE OF THE STUDY

Noncoding chloroplast DNA (NC-cpDNA) sequences are the staple data source of low-level phylogeographic and phylogenetic studies of angiosperms. We followed up on previous papers (tortoise and hare II and III) that sought to identify the most consistently variable regions of NC-cpDNA. We used an exhaustive literature review and newly available whole plastome data to assess applicability of previous conclusions at low taxonomic levels.

METHODS

We aligned complete plastomes of 25 species pairs from across angiosperms, comparing the number of genetic differences found in 107 NC-cpDNA regions and matK. We surveyed Web of Science for the plant phylogeographic literature between 2007 and 2013 to assess how NC-cpDNA has been used at the intraspecific level.

KEY RESULTS

Several regions are consistently the most variable across angiosperm lineages: ndhF-rpl32, rpl32-trnL((UAG)), ndhC-trnV((UAC)), 5'rps16-trnQ((UUG)), psbE-petL, trnT((GGU))-psbD, petA-psbJ, and rpl16 intron. However, there is no universally best region. The average number of regions applied to low-level studies is ∼2.5, which may be too little to access the full discriminating power of this genome.

CONCLUSIONS

Plastome sequences have been used successfully at lower and lower taxonomic levels. Our findings corroborate earlier works, suggesting that there are regions that are most likely to be the most variable. However, while NC-cpDNA sequences are commonly used in plant phylogeographic studies, few of the most variable regions are applied in that context. Furthermore, it appears that in most studies too few NC-cpDNAs are used to access the discriminating power of the cpDNA genome.

Arabinogalactan proteins mark stigmatic receptivity in the protogynous flowers of Magnolia virginiana (Magnoliaceae)

PREMISE OF THE STUDY

Factors affecting floral receptivity in angiosperms remain opaque, but recent studies suggest that the acquisition of stigmatic receptivity associated with cell-wall-related arabinogalactan proteins (AGPs) may be a widespread feature of flowering plants. Here, the time during which a stigma is receptive is evaluated and related to the secretion of AGPs in Magnolia virginiana, a protogynous member of an early-divergent angiosperm clade (magnoliids) with a clearly discernible female receptive phase.

METHODS

Magnolia virginiana flower phenology was documented, and histochemical changes in the stigma before and after pollination were examined. Stigmatic receptivity was evaluated in relation to the secretion of AGPs detected in whole mounts and immunolocalized in sectioned stigmas.

KEY RESULTS

Protogynous Magnolia flowers had a precise window of stigmatic receptivity, which is concomitant with the secretion of two AGPs labeled for different epitopes. After pollen germination and tube growth, these two AGPs could no longer be detected in the stigmas, suggesting that these AGPs interact with the growing male gametophytes and could be markers of stigmatic receptivity.

CONCLUSIONS

These results show that the period of stigmatic receptivity is finely coordinated with the secretion of two arabinogalactan proteins on stigmas of flowers of M. virginiana. This first report of AGP presence in stigmatic tissues in a member of the magnoliids, together with recently described similar patterns in eudicots, monocots, and members of early-divergent lineages of flowering plants, suggests an ancient and widespread role for AGPs on stigmatic receptivity in angiosperms.

Analysis of variation in chloroplast DNA sequences

This chapter introduces and reviews methods for analyzing variation in chloroplast DNA, mainly by polymerase chain reaction (PCR) and subsequent revelation of polymorphisms. Sources for chloroplast primers are discussed, as well as methods such as Sanger sequencing, PCR followed by restriction fragment length polymorphism (RFLP), gel electrophoresis, fragment analysis on automated DNA sequencers, denaturing high-performance liquid chromatography (dHPLC), and next-generation sequencing (NGS). A special section deals with peculiarities of chloroplast DNA variation, such as tandem repeats and mini- and microsatellites.

Studying genetic variability of pomegranate (Punica granatum L.) based on chloroplast DNA and barcode genes

Chloroplast DNA has been used extensively to analyze plant phylogenies at different taxonomic levels because of its size, organization and sequence conservation. In the present research, two chloroplastic regions, petA–psaJ, trnC–trnD and four DNA barcodes (trnH–psbA, ITS, rbcL, matK), were used to introduce suitable regions for the assessment of genetic diversity among P. granatum L. genotypes. Analysis of psbE–petL in petA–psaJ region revealed 1,300 nucleotides with 4.29 % genetic diversity among genotypes, while trnC–petN in trnC–trnD region showed 1.8 % genetic diversity. Therefore, despite the results obtained from the study of other plants, the trnC–trnD region had a low potential for the evaluation of diversity among pomegranate genotypes. Analysis of DNA barcodes in pomegranate showed that trnH–psbA (genetic diversity 2.91 %) provides the highest intra-species variation, followed by ITS (genetic diversity 0.44 %). Eighteen genotypes from different geographical origins of Iran were used to investigate psbE–petL and trnH–psbA potential as novel barcodes to determine genetic polymorphism and characterize pomegranate genotypes. The results suggested that two regions, psbE–petL and trnH–psbA, were more suitable for determining intra-species relationships of pomegranate.

Impact of gene molecular evolution on phylogenetic reconstruction: a case study in the rosids (Superorder Rosanae, Angiosperms)

Rate of substitution of genomic regions is among the most debated intrinsic features that impact phylogenetic informativeness. However, this variable is also coupled with rates of nonsynonymous substitutions that underscore the nature and degree of selection on the selected genes. To empirically address these variables, we constructed four completely overlapping data sets of plastid matK, atpB, rbcL, and mitochondrial matR genes and used the rosid lineage (angiosperms) as a working platform. The genes differ in combinations of overall rates of nucleotide and amino acid substitutions. Tree robustness, homoplasy, accuracy in contrast to a reference tree, and phylogenetic informativeness are evaluated. The rapidly evolving/unconstrained matK faired best, whereas remaining genes varied in degrees of contribution to rosid phylogenetics across the lineage's 108 million years evolutionary history. Phylogenetic accuracy was low with the slowly evolving/unconstrained matR despite least amount of homoplasy. Third codon positions contributed the highest amount of parsimony informative sites, resolution and informativeness, but magnitude varied with gene mode of evolution. These findings are in clear contrast with the views that rapidly evolving regions and the 3rd codon position have inevitable negative impact on phylogenetic reconstruction at deep historic level due to accumulation of multiple hits and subsequent elevation in homoplasy and saturation. Relaxed evolutionary constraint in rapidly evolving genes distributes substitutions across codon positions, an evolutionary mode expected to reduce the frequency of multiple hits. These findings should be tested at deeper evolutionary histories.

Reconstructing the age and historical biogeography of the ancient flowering-plant family Hydatellaceae (Nymphaeales)

Background

The aquatic flowering-plant family Hydatellaceae has a classic Gondwanan distribution, as it is found in Australia, India and New Zealand. To shed light on the biogeographic history of this apparently ancient branch of angiosperm phylogeny, we dated the family in the context of other seed-plant divergences, and evaluated its biogeography using parsimony and likelihood methods. We also explicitly tested the effect of different extinction rates on biogeographic inferences.

Results

We infer that the stem lineage of Hydatellaceae originated in the Lower Cretaceous; in contrast, its crown originated much more recently, in the early Miocene, with the bulk of its diversification after the onset of the Pliocene. Biogeographic reconstructions predict a mix of dispersal and vicariance events, but considerations of geological history preclude most vicariance events, besides a split at the root of the family between southern and northern clades. High extinction rates are plausible in the family, and when these are taken into account there is greater uncertainty in biogeographic inferences.

Conclusions

A stem origin for Hydatellaceae in the Lower Cretaceous is consistent with the initial appearance of fossils attributed to its sister clade, the water lilies. In contrast, the crown clade is young, indicating that vicariant explanations for species outside Australia are improbable. Although long-distance dispersal is likely the primary driver of biogeographic distribution in Hydatellaceae, we infer that the recent drying out of central Australia divided the family into tropical vs. subtropical/temperate clades around the beginning of the Miocene.

A comparison of 16 DNA regions for use as phylogenetic markers in the pleurocarpous moss genus Plagiothecium (Hypnales)

PREMISE OF THE STUDY

Within the Hypnales-the most derived and species-rich order of pleurocarpous mosses - phylogenies at or below the family level often show poor resolution. In preparation for a phylogeny of the genus Plagiothecium, we wished to identify the DNA markers best suited for evolutionary reconstruction in this group of hypnalean pleurocarps.

METHODS

For each of 25 collections of Plagiothecium and associated taxa, 16 DNA regions were sequenced: nuclear ITS and 26S, and plastid rps4, rps4-trnL, trnL-F, trnK (matK)-psbA, psbA-trnH, trnM-V, trnD-T, rbcL, atpB-rbcL, psbT-H, rpoC1 exon 2 (partial), the trnG intron, the rpl16 intron and the plastid ribosomal spacer DNA (cpITS). Each region was evaluated on the basis of its ability to resolve clades, the amount of homoplasy present in the data set, and the relative ease of obtaining the data. Descriptive statistics for each region are given.

KEY RESULTS

Under-utilized plastid markers for bryophytes such as trnK-psbA, rps4-trnL, and trnD-T outperformed more traditional markers such as trnL-F and rps4. Individual plastid topologies were similar, suggesting that only a limited amount of plastid data are needed to recover a backbone phylogeny. Adding a small amount of nuclear ribosomal data to a large plastid matrix restructured the recovered topology, emphasizing the importance of sampling multiple genomes and the need for new low-copy nuclear markers in bryophyte systematics.

CONCLUSIONS

Future genus-level phylogenies of pleurocarpous mosses should target under-utilized plastid markers such as trnK-psbA and rps4-trnL in conjunction with low-copy nuclear markers.

A long PCR-based approach for DNA enrichment prior to next-generation sequencing for systematic studies

PREMISE OF THE STUDY

We present an alternative approach for molecular systematic studies that combines long PCR and next-generation sequencing. Our approach can be used to generate templates from any DNA source for next-generation sequencing. Here we test our approach by amplifying complete chloroplast genomes, and we present a set of 58 potentially universal primers for angiosperms to do so. Additionally, this approach is likely to be particularly useful for nuclear and mitochondrial regions. •

METHODS AND RESULTS

Chloroplast genomes of 30 species across angiosperms were amplified to test our approach. Amplification success varied depending on whether PCR conditions were optimized for a given taxon. To further test our approach, some amplicons were sequenced on an Illumina HiSeq 2000. •

CONCLUSIONS

Although here we tested this approach by sequencing plastomes, long PCR amplicons could be generated using DNA from any genome, expanding the possibilities of this approach for molecular systematic studies.

Land plant evolutionary timeline: gene effects are secondary to fossil constraints in relaxed clock estimation of age and substitution rates

PREMISE OF THE STUDY

Land plants play an essential role in the evolution of terrestrial life. Their time of origin and diversification is fundamental to understanding the evolution of life on land. We investigated the timing and the rate of molecular evolution of land plants, evaluating the effects of different types of molecular data, including temporal information from fossils, and using different molecular clock methods. •

METHODS

Ages and absolute rates were estimated independently with two substitutionally different data sets: a highly conserved 4-gene data set and matK, a fast-evolving gene. The vascular plant backbone and the crown nodes of all major lineages were calibrated with fossil-derived ages. Dates and absolute rates were estimated while including or excluding the calibrations and using two relaxed clocks that differ in their implementation of temporal autocorrelation. •

KEY RESULTS

Land plants diverged from streptophyte alga 912 (870-962) million years ago (Mya) but diversified into living lineages 475 (471-480) Mya. Ages estimated for all major land-plant lineages agree with their fossil record, except for angiosperms. Different genes estimated very similar ages and correlated absolute rates across the tree. Excluding calibrations resulted in the greatest age differences. Different relaxed clocks provided similar ages, but different and uncorrelated absolute rates. •

CONCLUSIONS

Whole-genome rate accelerations or decelerations may underlie the similar ages and correlated absolute rates estimated with different genes. We suggest that pronounced substitution rate changes around the angiosperm crown node may represent a challenge for relaxed clocks to model adequately.

Re-mind the gap! Insertion - deletion data reveal neglected phylogenetic potential of the nuclear ribosomal internal transcribed spacer (ITS) of fungi

Rapidly evolving, indel-rich phylogenetic markers play a pivotal role in our understanding of the relationships at multiple levels of the tree of life. There is extensive evidence that indels provide conserved phylogenetic signal, however, the range of phylogenetic depths for which gaps retain tree signal has not been investigated in detail. Here we address this question using the fungal internal transcribed spacer (ITS), which is central in many phylogenetic studies, molecular ecology, detection and identification of pathogenic and non-pathogenic species. ITS is repeatedly criticized for indel-induced alignment problems and the lack of phylogenetic resolution above species level, although these have not been critically investigated. In this study, we examined whether the inclusion of gap characters in the analyses shifts the phylogenetic utility of ITS alignments towards earlier divergences. By re-analyzing 115 published fungal ITS alignments, we found that indels are slightly more conserved than nucleotide substitutions, and when included in phylogenetic analyses, improved the resolution and branch support of phylogenies across an array of taxonomic ranges and extended the resolving power of ITS towards earlier nodes of phylogenetic trees. Our results reconcile previous contradicting evidence for the effects of data exclusion: in the case of more sophisticated indel placement, the exclusion of indel-rich regions from the analyses results in a loss of tree resolution, whereas in the case of simpler alignment methods, the exclusion of gapped sites improves it. Although the empirical datasets do not provide to measure alignment accuracy objectively, our results for the ITS region are consistent with previous simulations studies alignment algorithms. We suggest that sophisticated alignment algorithms and the inclusion of indels make the ITS region and potentially other rapidly evolving indel-rich loci valuable sources of phylogenetic information, which can be exploited at multiple taxonomic levels.

Discrimination of Phellodendron amurense and P. chinense based on DNA analysis and the simultaneous analysis of alkaloids

Phellodendri Cortex is the bark of the stems of Phellodendron amurense Ruprecht or P. chinense Schneider (Rutaceae), which is orginated from periderm. The internal transcribed spacer sequences of 20 originated plants and identified samples were analyzed. The result showed that the 99% of the base sequences of P. amurense were identical to that of P. chinense, but the differentiation of P. amurense and P. chinense was difficult. In addition, the ribulose-1, 5-bisphospate carboxylase large subunit (rbcL) intergenic spacer sequences of specific parts produced the same result. However, when the analysis was carried out by using the RAPD (randomly amplification polymorphism DNA) analysis method, which utilizes 48 randomly primers, it allowed us to confirm the polymorphism of P. amurense and P. chinense in 12 primers. A high-performance liquid chromatographic (HPLC) method was developed and validated for the simultaneous quantitation of berberine, palmatine and jatrorrhizine in a traditional herbal drug, Phellodendri Cortex. The HPLC method was applied successfully to the quantification of three constituents in the extract of twenty Phellodendri Cortex. The results indicated that the established HPLC and RAPD methods are suitable for the quantitative analysis and the quality control multi-simultaneous discrimination in Phellodendri Cortex.

Plastid genome evolution in mycoheterotrophic Ericaceae

Unlike parasitic plants, which are linked to their hosts directly through haustoria, mycoheterotrophic (MHT) plants derive all or part of their water and nutrients from autothrophs via fungal mycorrhizal intermediaries. Ericaceae, the heather family, are a large and diverse group of plants known to form elaborate symbiotic relationships with mycorrhizal fungi. Using PHYA sequence data, we first investigated relationships among mycoheterotrophic Ericaceae and their close autotrophic relatives. Phylogenetic results suggest a minimum of two independent origins of MHT within this family. Additionally, a comparative investigation of plastid genomes (plastomes) grounded within this phylogenetic framework was conducted using a slot-blot Southern hybridization approach. This survey encompassed numerous lineages of Ericaceae with different life histories and trophic levels, including multiple representatives from mixotrophic Pyroleae and fully heterotrophic Monotropeae and Pterosporeae. Fifty-four probes derived from all categories of protein coding genes typically found within the plastomes of flowering plants were used. Our results indicate that the holo-mycoheterotrophic Ericaceae exhibit extensive loss of genes relating to photosynthetic function and expression of the plastome but retain genes with possible functions outside photosynthesis. Mixotrophic taxa tend to retain most genes relating to photosynthetic functions but are varied regarding the plastid ndh gene content. This investigation extends previous inferences that the loss of the NDH complex occurs prior to becoming holo-heterotrophic and it shows that the pattern of gene losses among mycoheterotrophic Ericaceae is similar to that of haustorial parasites. Additionally, we identify the most desirable candidate species for entire plastome sequencing.

Molecular phylogenetics of Hydatellaceae (Nymphaeales): sexual-system homoplasy and a new sectional classification

PREMISE OF THE STUDY

Species relationships are unknown in Hydatellaceae, a small family of dwarf aquatics related to water lilies that arose near the base of angiosperm phylogeny. Here we use molecular evidence to infer a species tree for the family and apply this to reconstructing major transitions in morphology and sexual system in this early branch of angiosperms.

METHODS

We assembled plastid (atpB, matK, ndhF, rbcL) and nuclear (ribosomal ITS) data for 50 samples (including outgroups) and estimated a species tree for Hydatellaceae using a Bayesian multispecies coalescent approach. We reconstructed the evolution of several morphological characters, then tested for associations between sexual system and reproductive morphology using phylogenetic ANOVA.

KEY RESULTS

Dioecious species of Hydatellaceae have significantly greater stamen number and anther length than do cosexual species, suggesting changes in male function. The perennial habit that defines one subclade likely represents a reversion from annuality. Species relationships do not fall along traditional morphological divisions, but new sections proposed here are supported by fruit and seed synapomorphies. The earliest split in the family is reflected in geography and climate (i.e., tropical vs. subtropical/temperate clades). We found limited evidence of incongruence between plastid and nuclear trees, with one exception involving gene-tree nonmonophyly for two close relatives (Trithuria submersa, T. bibracteata).

CONCLUSIONS

While the direction of sexual-system evolution is ambiguous, transitions are significantly associated with changes in involucral phyllome length and proxies of pollen production. We propose a new sectional circumscription based on fruit, seed, and DNA evidence.

Reconciling gene and genome duplication events: using multiple nuclear gene families to infer the phylogeny of the aquatic plant family Pontederiaceae

Most plant phylogenetic inference has used DNA sequence data from the plastid genome. This genome represents a single genealogical sample with no recombination among genes, potentially limiting the resolution of evolutionary relationships in some contexts. In contrast, nuclear DNA is inherently more difficult to employ for phylogeny reconstruction because major mutational events in the genome, including polyploidization, gene duplication, and gene extinction can result in homologous gene copies that are difficult to identify as orthologs or paralogs. Gene tree parsimony (GTP) can be used to infer the rooted species tree by fitting gene genealogies to species trees while simultaneously minimizing the estimated number of duplications needed to reconcile conflicts among them. Here, we use GTP for five nuclear gene families and a previously published plastid data set to reconstruct the phylogenetic backbone of the aquatic plant family Pontederiaceae. Plastid-based phylogenetic studies strongly supported extensive paraphyly of Eichhornia (one of the four major genera) but also depicted considerable ambiguity concerning the true root placement for the family. Our results indicate that species trees inferred from the nuclear genes (alone and in combination with the plastid data) are highly congruent with gene trees inferred from plastid data alone. Consideration of optimal and suboptimal gene tree reconciliations place the root of the family at (or near) a branch leading to the rare and locally restricted E. meyeri. We also explore methods to incorporate uncertainty in individual gene trees during reconciliation by considering their individual bootstrap profiles and relate inferred excesses of gene duplication events on individual branches to whole-genome duplication events inferred for the same branches. Our study improves understanding of the phylogenetic history of Pontederiaceae and also demonstrates the utility of GTP for phylogenetic analysis.

A multi gene sequence-based phylogeny of the Musaceae (banana) family

Background

The classification of the Musaceae (banana) family species and their phylogenetic inter-relationships remain controversial, in part due to limited nucleotide information to complement the morphological and physiological characters. In this work the evolutionary relationships within the Musaceae family were studied using 13 species and DNA sequences obtained from a set of 19 unlinked nuclear genes.

Results

The 19 gene sequences represented a sample of ~16 kb of genome sequence (~73% intronic). The sequence data were also used to obtain estimates for the divergence times of the Musaceae genera and Musa sections. Nucleotide variation within the sample confirmed the close relationship of Australimusa and Callimusa sections and showed that Eumusa and Rhodochlamys sections are not reciprocally monophyletic, which supports the previous claims for the merger between the two latter sections. Divergence time analysis supported the previous dating of the Musaceae crown age to the Cretaceous/Tertiary boundary (~ 69 Mya), and the evolution of Musa to ~50 Mya. The first estimates for the divergence times of the four Musa sections were also obtained.

Conclusions

The gene sequence-based phylogeny presented here provides a substantial insight into the course of speciation within the Musaceae. An understanding of the main phylogenetic relationships between banana species will help to fine-tune the taxonomy of Musaceae.

The evolution of chloroplast genome structure in ferns

The plastid genome (plastome) is a rich source of phylogenetic and other comparative data in plants. Most land plants possess a plastome of similar structure. However, in a major group of plants, the ferns, a unique plastome structure has evolved. The gene order in ferns has been explained by a series of genomic inversions relative to the plastome organization of seed plants. Here, we examine for the first time the structure of the plastome across fern phylogeny. We used a PCR-based strategy to map and partially sequence plastomes. We found that a pair of partially overlapping inversions in the region of the inverted repeat occurred in the common ancestor of most ferns. However, the ancestral (seed plant) structure is still found in early diverging branches leading to the osmundoid and filmy fern lineages. We found that a second pair of overlapping inversions occurred on a branch leading to the core leptosporangiates. We also found that the unique placement of the gene matK in ferns (lacking a flanking intron) is not a result of a large-scale inversion, as previously thought. This is because the intron loss maps to an earlier point on the phylogeny than the nearby inversion. We speculate on why inversions may occur in pairs and what this may mean for the dynamics of plastome evolution.

Characterization of the complete chloroplast genome of Hevea brasiliensis reveals genome rearrangement, RNA editing sites and phylogenetic relationships

Rubber tree (Hevea brasiliensis) is an economical plant and widely grown for natural rubber production. However, genomic research of rubber tree has lagged behind other species in the Euphorbiaceae family. We report the complete chloroplast genome sequence of rubber tree as being 161,191 bp in length including a pair of inverted repeats of 26,810 bp separated by a small single copy region of 18,362 bp and a large single copy region of 89,209 bp. The chloroplast genome contains 112 unique genes, 16 of which are duplicated in the inverted repeat. Of the 112 unique genes, 78 are predicted protein-coding genes, 4 are ribosomal RNA genes and 30 are tRNA genes. Relative to other plant chloroplast genomes, we observed a unique rearrangement in the rubber tree chloroplast genome: a 30-kb inversion between the trnE(UUC)-trnS(GCU) and the trnT(GGU)-trnR(UCU). A comparison between the rubber tree chloroplast genes and cDNA sequences revealed 51 RNA editing sites in which most (48 sites) were located in 26 protein coding genes and the other 3 sites were in introns. Phylogenetic analysis based on chloroplast genes demonstrated a close relationship between Hevea and Manihot in Euphorbiaceae and provided a strong support for a monophyletic group of the eurosid I.