Organelle Phylogenomics and Extensive Conflicting Phylogenetic Signals in the Monocot Order Poales

New insights into the phylogenetic relationships within the Lauraceae from mitogenomes

BACKGROUND

The family Lauraceae is subdivided into six main lineages: Caryodaphnopsideae, Cassytheae, Cryptocaryeae, Hypodaphnideae, Laureae, and Neocinnamomeae. However, phylogenetic relationships among these lineages have been debatable due to incongruence between trees constructed using nuclear ribosomal DNA (nrDNA) sequences and chloroplast (cp) genomes. As with cp DNA, the mitochondrial (mt) DNA of most flowering plants is maternally inherited, so the phylogenetic relationships recovered with mt genomes are expected to be consistent with that from cp genomes, rather than nrDNA sequences.

RESULTS

The mitogenome of Machilus yunnanensis, with a length of 735,392 bp, has a very different genome size and gene linear order from previously published magnoliid mitogenomes. Phylogenomic reconstructions based on 41 mt genes from 92 Lauraceae mitogenomes resulted in highly supported relationships: sisterhood of the Laureae and a group containing Neocinnamomeae and Caryodaphnopsideae, with Cassytheae being the next sister group, followed by Cryptocaryeae. However, we found significant incongruence among the mitochondrial, chloroplast, and nuclear phylogenies, especially for the species within the Caryodaphnopsideae and Neocinnamomeae lineages. Time-calibrated phylogenetic analyses showed that the split between Caryodaphnopsideae and Neocinnamomeae dated to the later Eocene, around 38.5 Ma, Laureae originated in the Late Cretaceous, around 84.9 Ma, Cassytheae originated in the mid-Cretaceous around 102 Ma, and Cryptocaryeae originated in the Early Cretaceous around 116 Ma. From the Late Cretaceous to the Paleocene, net diversification rates significantly increased across extant clades of major lineages, and both speciation rates and net diversification rates continued steady growth towards the present.

CONCLUSIONS

The topology obtained here for the first time shows that mt genes can be used to support relationships among lineages of Lauraceae. Our results highlight that both Caryodaphnopsideae and Neocinnamomeae lineages are younger than previously thought, likely first diversifying in the Eocene, and species in the other extant lineages of Lauraceae dates in a long-time span from the Early Cretaceous to the Eocene, and the climate of a period of about 90 million years was relatively warm, while the extant species of Lauraceae then continuously diversified with global cooling from the Eocene to the present day.

The complete chloroplast genome of Typha angustifolia and its phylogenetic position in Typhaceae

Typha angustifolia, commonly known as narrowleaf cattail, is a marginal, semi-aquatic, herbaceous perennial species with both ecological and edible values. In this study, the complete chloroplast (cp) genome of T. angustifolia was assembled using the next-generation sequencing technology. The whole cp genome was 161,597 bp in length, consisting of a large single copy (LSC, 89,119 bp) and a small single copy (SSC, 18,550 bp) separated by two copies of inverted region (IR, 26,964 bp). The genome encoded 113 unique genes, including 79 protein-coding genes, 30 tRNA genes, four rRNA genes, with 19 duplicated genes in the IR regions. Phylogenetic analysis showed that T. angustifolia is sister to Typha orientalis in the family Typhaceae. The cp genome of T. angustifolia is reported for the first time, which will provide essential and important genetic resources for future phylogenetic investigation within the genus Typha.

The complete chloroplast genome and phylogenetic analysis of Cyperus brevifolius (Rottb.) Hassk. 1844 (Cyperaceae)

The medicinal plant Cyperus brevifolius (Rottb.) Hassk. 1844 has a long history of use in traditional Chinese medicine. In this study, we determined and systematically analyzed the complete chloroplast (cp) genome of C. brevifolius. The genome is 183,717 bp in length with a GC content of 33.24%. It comprises four distinct regions: a large-single copy (LSC) region of 101,190 bp, a small-single copy (SSC) region of 10,366 bp, and two inverted repeat (IR) regions of 36,079 bp each. A total of 137 genes are present in the genome including 89 protein-coding genes, 40 tRNA genes, and eight rRNA genes. Phylogenetic analysis reveals that C. brevifolius belongs to the Cyperus genus. This newly sequenced cp genome provides valuable insights for future genetic and genomic studies on Cyperus.

Unprecedented variation pattern of plastid genomes and the potential role in adaptive evolution in Poales

BACKGROUND

The plastid is the photosynthetic organelle in plant cell, and the plastid genomes (plastomes) are generally conserved in evolution. As one of the most economically and ecologically important order of angiosperms, Poales was previously documented to exhibit great plastomic variation as an order of photoautotrophic plants.

RESULTS

We acquired 93 plastomes, representing all the 16 families and 5 major clades of Poales to reveal the extent of their variation and evolutionary pattern. Extensive variation including the largest one in monocots with 225,293 bp in size, heterogeneous GC content, and a wide variety of gene duplication and loss were revealed. Moreover, rare occurrences of three inverted repeat (IR) copies in angiosperms and one IR loss were observed, accompanied by short IR (sIR) and small direct repeat (DR). Widespread structural heteroplasmy, diversified inversions, and unusual genomic rearrangements all appeared in Poales, occasionally within a single species. Extensive repeats in the plastomes were found to be positively correlated with the observed inversions and rearrangements. The variation all showed a "small-large-moderate" trend along the evolution of Poales, as well as for the sequence substitution rate. Finally, we found some positively selected genes, mainly in C4 lineages, while the closely related lineages of those experiencing gene loss tended to have undergone more relaxed purifying selection.

CONCLUSIONS

The variation of plastomes in Poales may be related to its successful diversification into diverse habitats and multiple photosynthetic pathway transitions. Our order-scale analyses revealed unusual evolutionary scenarios for plastomes in the photoautotrophic order of Poales and provided new insights into the plastome evolution in angiosperms as a whole.

The low level of plastome differentiation observed in some lineages of Poales hinders molecular species identification

Chloroplast genomes are a source of information successfully used in various fields of plant genetics, including molecular species identification. However, recent studies indicate an extremely low level of interspecific variability in the plastomes of some taxonomic groups of plants, including the genus Stipa L., which is a representative of the grass family. In this study we aimed to analyze the level of chloroplast genome diversity within particular genera as well as the effectiveness of identifying plant species in the Poaceae family and the other representatives of Poales order. Analysis of complete plastid genome alignments created for 96 genera comprising 793 species and 1707 specimens obtained from the GenBank database allowed defining and categorizing molecular diagnostic characters distinguishing the analyzed species from the other representatives of the genus. The results also demonstrate which species do not have any species-specific mutations, thereby they cannot be identified on the basis of differences between the complete chloroplast genomes. Our research showed a huge diversity of the analyzed species in terms of the number of molecular diagnostic characters and indicated which genera pose a particular challenge in terms of molecular species identification. The results show that a very low level of genetic diversity between plastomes is not uncommon in Poales. This is the first extensive research on super-barcoding that tests this method on a large data set and illustrates its effectiveness against the background of phylogenetic relationships.

Intrageneric structural variation in organelle genomes from the genus Dystaenia (Apiaceae): genome rearrangement and mitochondrion-to-plastid DNA transfer

Introduction

During plant evolution, intracellular DNA transfer (IDT) occurs not only from organelles to the nucleus but also between organelles. To further comprehend these events, both organelle genomes and transcriptomes are needed.

Methods

In this study, we constructed organelle genomes and transcriptomes for two Dystaenia species and described their dynamic IDTs between their nuclear and mitochondrial genomes, or plastid and mitochondrial genomes (plastome and mitogenome).

Results and Discussion

We identified the putative functional transfers of the mitochondrial genes 5' rpl2, rps10, rps14, rps19, and sdh3 to the nucleus in both Dystaenia species and detected two transcripts for the rpl2 and sdh3 genes. Additional transcriptomes from the Apicaceae species also provided evidence for the transfers and duplications of these mitochondrial genes, showing lineage-specific patterns. Intrageneric variations of the IDT were found between the Dystaenia organelle genomes. Recurrent plastid-to-mitochondrion DNA transfer events were only identified in the D. takeshimana mitogenome, and a pair of mitochondrial DNAs of plastid origin (MIPTs) may generate minor alternative isoforms. We only found a mitochondrion-to-plastid DNA transfer event in the D. ibukiensis plastome. This event may be linked to inverted repeat boundary shifts in its plastome. We inferred that the insertion region involved an MIPT that had already acquired a plastid sequence in its mitogenome via IDT. We propose that the MIPT acts as a homologous region pairing between the donor and recipient sequences. Our results provide insight into the evolution of organelle genomes across the family Apiaceae.

Phylogenomics and the flowering plant tree of life

The advances accelerated by next-generation sequencing and long-read sequencing technologies continue to provide an impetus for plant phylogenetic study. In the past decade, a large number of phylogenetic studies adopting hundreds to thousands of genes across a wealth of clades have emerged and ushered plant phylogenetics and evolution into a new era. In the meantime, a roadmap for researchers when making decisions across different approaches for their phylogenomic research design is imminent. This review focuses on the utility of genomic data (from organelle genomes, to both reduced representation sequencing and whole-genome sequencing) in phylogenetic and evolutionary investigations, describes the baseline methodology of experimental and analytical procedures, and summarizes recent progress in flowering plant phylogenomics at the ordinal, familial, tribal, and lower levels. We also discuss the challenges, such as the adverse impact on orthology inference and phylogenetic reconstruction raised from systematic errors, and underlying biological factors, such as whole-genome duplication, hybridization/introgression, and incomplete lineage sorting, together suggesting that a bifurcating tree may not be the best model for the tree of life. Finally, we discuss promising avenues for future plant phylogenomic studies.

Molecular evolution of the Pi-d2 gene conferring resistance to rice blast in Oryza

The exploitation of plant disease resistance (R) genes in breeding programs is an effective strategy for coping with pathogens. An understanding of R gene variation is the basis for this strategy. Rice blast disease, caused by the Magnaporthe oryzae fungus, is a destructive disease of rice. The rice blast resistance gene Pi-d2 represents a new class of plant R gene because of its novel extracellular domain. We investigated the nucleotide polymorphism, phylogenetic topology and evolution patterns of the Pi-d2 gene among 67 cultivated and wild rice relatives. The Pi-d2 gene originated early in the basal Poales and has remained as a single gene without expansion. The striking finding is that susceptible Pi-d2 alleles might be derived from a single nucleotide substitution of the resistant alleles after the split of Oryza subspecies. Functional pleiotropy and linkage effects are proposed for the evolution and retention of the disease-susceptible alleles in rice populations. One set of DNA primers was developed from the polymorphic position to detect the functional nucleotide polymorphism for disease resistance of the Pi-d2 gene based on conventional Polymerase Chain Reaction. The nucleotide diversity level varied between different domains of the Pi-d2 gene, which might be related to distinct functions of each domain in the disease defense response. Directional (or purifying) selection appears dominant in the molecular evolution of the Pi-d2 gene and has shaped its conserved variation pattern.

Evolutionary history of the grass gynoecium

The grass family (Poaceae) includes cereal crops that provide a key food source for the human population. The food industry uses the starch deposited in the cereal grain, which develops directly from the gynoecium. Morphological interpretation of the grass gynoecium remains controversial. We re-examine earlier hypotheses and studies of morphology and development in the context of more recent analyses of grass phylogenetics and developmental genetics. Taken in isolation, data on gynoecium development in bistigmatic grasses do not contradict its interpretation as a solitary ascidiate carpel. Nevertheless, in the context of other data, this interpretation is untenable. Broad comparative analysis in a modern phylogenetic context clearly demonstrates that the grass gynoecium is pseudomonomerous. A bistigmatic grass gynoecium has two sterile carpels, each producing a stigma, and a fertile carpel that lacks a stigma. To date, studies of grass developmental genetics and developmental morphology have failed to fully demonstrate the composite nature of the grass gynoecium be-cause its complex evolutionary history is hidden by extreme organ integration. It is problematic to interpret the gynoecium of grasses in terms of normal angiosperm gynoecium typology. Even the concept of a carpel becomes misleading in grasses; instead, we recommend the term pistil for descriptive purposes.

New plastome structural rearrangements discovered in core Tillandsioideae (Bromeliaceae) support recently adopted taxonomy

Full plastome sequences for land plants have become readily accessible thanks to the development of Next Generation Sequencing (NGS) techniques and powerful bioinformatic tools. Despite this vast amount of genomic data, some lineages remain understudied. Full plastome sequences from the highly diverse (>1,500 spp.) subfamily Tillandsioideae (Bromeliaceae, Poales) have been published for only three (i.e., Guzmania, Tillandsia, and Vriesea) out of 22 currently recognized genera. Here, we focus on core Tillandsioideae, a clade within subfamily Tillandsioideae, and explore the contribution of individual plastid markers and data categories to inform deep divergences of a plastome phylogeny. We generated 37 high quality plastome assemblies and performed a comparative analysis in terms of plastome structure, size, gene content and order, GC content, as well as number and type of repeat motifs. Using the obtained phylogenetic context, we reconstructed the evolution of these plastome attributes and assessed if significant shifts on the evolutionary traits' rates have occurred in the evolution of the core Tillandsioideae. Our results agree with previously published phylogenetic hypotheses based on plastid data, providing stronger statistical support for some recalcitrant nodes. However, phylogenetic discordance with previously published nuclear marker-based hypotheses was found. Several plastid markers that have been consistently used to address phylogenetic relationships within Tillandsioideae were highly informative for the retrieved plastome phylogeny and further loci are here identified as promising additional markers for future studies. New lineage-specific plastome rearrangements were found to support recently adopted taxonomic groups, including large inversions, as well as expansions and contractions of the inverted repeats. Evolutionary trait rate shifts associated with changes in size and GC content of the plastome regions were found across the phylogeny of core Tillandsioideae.