Identification of shared single copy nuclear genes in Arabidopsis, Populus, Vitis and Oryza and their phylogenetic utility across various taxonomic levels

Cytonuclear evolution in fully heterotrophic plants: lifestyles and gene function determine scenarios

BACKGROUND

Evidence shows that full mycoheterotrophs and holoparasites often have reduced plastid genomes with rampant gene loss, elevated substitution rates, and deeply altered to conventional evolution in mitochondrial genomes, but mechanisms of cytonuclear evolution is unknown. Endoparasitic Sapria himalayana and mycoheterotrophic Gastrodia and Platanthera guangdongensis represent different heterotrophic types, providing a basis to illustrate cytonuclear evolution. Here, we focused on nuclear-encoded plastid / mitochondrial (N-pt / mt) -targeting protein complexes, including caseinolytic protease (ClpP), ribulose-1,5-bisphosphate carboxylase/oxygenase (RuBisCo), oxidative phosphorylation system (OXPHOS), DNA recombination, replication, and repair (DNA-RRR) system, and pentatricopeptide repeat (PPR) proteins, to identify evolutionary drivers for cytonuclear interaction.

RESULTS

The severity of gene loss of N-pt PPR and pt-RRR genes was positively associated with increased degree of heterotrophy in full mycoheterotrophs and S. himalayana, while N-mt PPR and mt-RRR genes were retained. Substitution rates of organellar and nuclear genes encoding N-pt/mt subunits in protein complexes were evaluated, cytonuclear coevolution was identified in S. himalayana, whereas disproportionate rates of evolution were observed in the OXPHOS complex in full mycoheterotrophs, only slight accelerations in substitution rates were identified in N-mt genes of full mycoheterotrophs.

CONCLUSIONS

Nuclear compensatory evolution was identified in protein complexes encoded by plastid and N-pt genes. Selection shaping codon preferences, functional constraint, mt-RRR gene regulation, and post-transcriptional regulation of PPR genes all facilitate mito-nuclear evolution. Our study enriches our understanding of genomic coevolution scenarios in fully heterotrophic plants.

Developing Asparagaceae1726: An Asparagaceae-specific probe set targeting 1726 loci for Hyb-Seq and phylogenomics in the family

Premise

Target sequence capture (Hyb-Seq) is a cost-effective sequencing strategy that employs RNA probes to enrich for specific genomic sequences. By targeting conserved low-copy orthologs, Hyb-Seq enables efficient phylogenomic investigations. Here, we present Asparagaceae1726-a Hyb-Seq probe set targeting 1726 low-copy nuclear genes for phylogenomics in the angiosperm family Asparagaceae-which will aid the often-challenging delineation and resolution of evolutionary relationships within Asparagaceae.

Methods

Here we describe and validate the Asparagaceae1726 probe set (https://github.com/bentzpc/Asparagaceae1726) in six of the seven subfamilies of Asparagaceae. We perform phylogenomic analyses with these 1726 loci and evaluate how inclusion of paralogs and bycatch plastome sequences can enhance phylogenomic inference with target-enriched data sets.

Results

We recovered at least 82% of target orthologs from all sampled taxa, and phylogenomic analyses resulted in strong support for all subfamilial relationships. Additionally, topology and branch support were congruent between analyses with and without inclusion of target paralogs, suggesting that paralogs had limited effect on phylogenomic inference.

Discussion

Asparagaceae1726 is effective across the family and enables the generation of robust data sets for phylogenomics of any Asparagaceae taxon. Asparagaceae1726 establishes a standardized set of loci for phylogenomic analysis in Asparagaceae, which we hope will be widely used for extensible and reproducible investigations of diversification in the family.

Shared single copy genes are generally reliable for inferring phylogenetic relationships among polyploid taxa

Polyploidy, or whole-genome duplication, is expected to confound the inference of species trees with phylogenetic methods for two reasons. First, the presence of retained duplicated genes requires the reconciliation of the inferred gene trees to a proposed species tree. Second, even if the analyses are restricted to shared single copy genes, the occurrence of reciprocal gene loss, where the surviving genes in different species are paralogs from the polyploidy rather than orthologs, will mean that such genes will not have evolved under the corresponding species tree and may not produce gene trees that allow inference of that species tree. Here we analyze three different ancient polyploidy events, using synteny-based inferences of orthology and paralogy to infer gene trees from nearly 17,000 sets of homologous genes. We find that the simple use of single copy genes from polyploid organisms provides reasonably robust phylogenetic signals, despite the presence of reciprocal gene losses. Such gene trees are also most often in accord with the inferred species relationships inferred from maximum likelihood models of gene loss after polyploidy: a completely distinct phylogenetic signal present in these genomes. As seen in other studies, however, we find that methods for inferring phylogenetic confidence yield high support values even in cases where the underlying data suggest meaningful conflict in the phylogenetic signals.

The wild side of grape genomics

With broad genetic diversity and as a source of key agronomic traits, wild grape species (Vitis spp.) are crucial to enhance viticulture's climatic resilience and sustainability. This review discusses how recent breakthroughs in the genome assembly and analysis of wild grape species have led to discoveries on grape evolution, from wild species' adaptation to environmental stress to grape domestication. We detail how diploid chromosome-scale genomes from wild Vitis spp. have enabled the identification of candidate disease-resistance and flower sex determination genes and the creation of the first Vitis graph-based pangenome. Finally, we explore how wild grape genomics can impact grape research and viticulture, including aspects such as data sharing, the development of functional genomics tools, and the acceleration of genetic improvement.

Evolution of Whirly1 in the angiosperms: sequence, splicing, and expression in a clade of early transitional mycoheterotrophic orchids

The plastid-targeted transcription factor Whirly1 (WHY1) has been implicated in chloroplast biogenesis, plastid genome stability, and fungal defense response, which together represent characteristics of interest for the study of autotrophic losses across the angiosperms. While gene loss in the plastid and nuclear genomes has been well studied in mycoheterotrophic plants, the evolution of the molecular mechanisms impacting genome stability is completely unknown. Here, we characterize the evolution of WHY1 in four early transitional mycoheterotrophic orchid species in the genus Corallorhiza by synthesizing the results of phylogenetic, transcriptomic, and comparative genomic analyses with WHY1 genomic sequences sampled from 21 orders of angiosperms. We found an increased number of non-canonical WHY1 isoforms assembled from all but the greenest Corallorhiza species, including intron retention in some isoforms. Within Corallorhiza, phylotranscriptomic analyses revealed the presence of tissue-specific differential expression of WHY1 in only the most photosynthetically capable species and a coincident increase in the number of non-canonical WHY1 isoforms assembled from fully mycoheterotrophic species. Gene- and codon-level tests of WHY1 selective regimes did not infer significant signal of either relaxed selection or episodic diversifying selection in Corallorhiza but did so for relaxed selection in the late-stage full mycoheterotrophic orchids Epipogium aphyllum and Gastrodia elata. Additionally, nucleotide substitutions that most likely impact the function of WHY1, such as nonsense mutations, were only observed in late-stage mycoheterotrophs. We propose that our findings suggest that splicing and expression changes may precede the selective shifts we inferred for late-stage mycoheterotrophic species, which therefore does not support a primary role for WHY1 in the transition to mycoheterotrophy in the Orchidaceae. Taken together, this study provides the most comprehensive view of WHY1 evolution across the angiosperms to date.

Development of a Target Enrichment Probe Set for Conifer (REMcon)

Conifers are an ecologically and economically important seed plant group that can provide significant insights into the evolution of land plants. Molecular phylogenetics has developed as an important approach in evolutionary studies, although there have been relatively few studies of conifers that employ large-scale data sourced from multiple nuclear genes. Target enrichment sequencing (target capture, exon capture, or Hyb-Seq) has developed as a key approach in modern phylogenomic studies. However, until now, there has been no bait set that specifically targets the entire conifer clade. REMcon is a target sequence capture probe set intended for family- and species-level phylogenetic studies of conifers that target c. 100 single-copy nuclear loci. We tested the REMcon probe set using 69 species, including 44 conifer genera across six families and four other gymnosperm taxa, to evaluate the efficiency of target capture to efficiently generate comparable DNA sequence data across conifers. The recovery of target loci was high, with, on average, 94% of the targeted regions recovered across samples with high read coverage. A phylogenetic analysis of these data produced a well-supported topology that is consistent with the current understanding of relationships among conifers. The REMcon bait set will be useful in generating relatively large-scale nuclear data sets consistently for any conifer lineage.

Validation of a simplified small-scale DNA extraction protocol from wine by quantitative real-time PCR

In the present study, we compared a simplified small-scale purification protocol to obtain DNA admixtures out of wine, with our large-scale published method. The extraction methods must provide DNA free of PCR inhibitors, that can interfere with DNA amplification. To evaluate the efficiency of grapevine's nuclear DNA extraction from wine, the new protocol was also compared in terms of purity and yield to the DNA obtained out of grapevine's (Vitis vinifera) leaf tissue, using a commercial kit. Two single-copy nuclear genes, nine-cis-epoxy carotenoid dioxygenase 2 (NCED2), and prefoldin subunit 5-like (PS5) were amplified in DNA extracted from wine and grapevine by real-time TaqMan PCR to determine the presence of inhibitors in relation to the diversity of starting biological matrix. This study showed that the small-scale, simpler method for extracting DNA from wine produced effective results in terms of inhibitor presence and purity. Furthermore, even though the initial biological matrix was more complicated, the grapevine nuclear DNA that was removed from wine was qualitatively equivalent to the DNA that was isolated from the leaves.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-024-03992-x.

Overexpression of SLIM1 transcription factor accelerates vegetative development in Arabidopsis thaliana

The transcription factor Sulfur Limitation 1 (SLIM1) belongs to the plant-specific Ethylene Insenstive3-Like transcription factor family and is known to coordinate gene expression in response to sulfur deficiency. However, the roles of SLIM1 in nutrient-sufficient conditions have not been characterized. Employing constitutive SLIM1 overexpression (35S::SLIM1) and CRISPR/Cas9 mutant plants (slim1-cr), we identified several distinct phenotypes in nutrient-sufficient conditions in Arabidopsis thaliana. Overexpression of SLIM1 results in plants with approximately twofold greater rosette area throughout vegetative development. 35S::SLIM1 plants also bolt earlier and exhibit earlier downregulation of photosynthesis-associated genes and earlier upregulation of senescence-associated genes than Col-0 and slim1-cr plants. This suggests that overexpression of SLIM1 accelerates development in A. thaliana. Genome-wide differential gene expression analysis relative to Col-0 at three time points with slim1-cr and two 35S::SLIM1 lines allowed us to identify 1,731 genes regulated directly or indirectly by SLIM1 in vivo.

Genome copy number predicts extreme evolutionary rate variation in plant mitochondrial DNA

Nuclear and organellar genomes can evolve at vastly different rates despite occupying the same cell. In most bilaterian animals, mitochondrial DNA (mtDNA) evolves faster than nuclear DNA, whereas this trend is generally reversed in plants. However, in some exceptional angiosperm clades, mtDNA substitution rates have increased up to 5,000-fold compared with closely related lineages. The mechanisms responsible for this acceleration are generally unknown. Because plants rely on homologous recombination to repair mtDNA damage, we hypothesized that mtDNA copy numbers may predict evolutionary rates, as lower copy numbers may provide fewer templates for such repair mechanisms. In support of this hypothesis, we found that copy number explains 47% of the variation in synonymous substitution rates of mtDNA across 60 diverse seed plant species representing ~300 million years of evolution. Copy number was also negatively correlated with mitogenome size, which may be a cause or consequence of mutation rate variation. Both relationships were unique to mtDNA and not observed in plastid DNA. These results suggest that homologous recombinational repair plays a role in driving mtDNA substitution rates in plants and may explain variation in mtDNA evolution more broadly across eukaryotes. Our findings also contribute to broader questions about the relationships between mutation rates, genome size, selection efficiency, and the drift-barrier hypothesis.

Genomic decoding of Theobroma grandiflorum (cupuassu) at chromosomal scale: evolutionary insights for horticultural innovation

BACKGROUND

Theobroma grandiflorum (Malvaceae), known as cupuassu, is a tree indigenous to the Amazon basin, valued for its large fruits and seed pulp, contributing notably to the Amazonian bioeconomy. The seed pulp is utilized in desserts and beverages, and its seed butter is used in cosmetics. Here, we present the sequenced telomere-to-telomere genome of cupuassu, disclosing its genomic structure, evolutionary features, and phylogenetic relationships within the Malvaceae family.

FINDINGS

The cupuassu genome spans 423 Mb, encodes 31,381 genes distributed in 10 chromosomes, and exhibits approximately 65% gene synteny with the Theobroma cacao genome, reflecting a conserved evolutionary history, albeit punctuated with unique genomic variations. The main changes are pronounced by bursts of long-terminal repeat retrotransposons at postspecies divergence, retrocopied and singleton genes, and gene families displaying distinctive patterns of expansion and contraction. Furthermore, positively selected genes are evident, particularly among retained and dispersed tandem and proximal duplicated genes associated with general fruit and seed traits and defense mechanisms, supporting the hypothesis of potential episodes of subfunctionalization and neofunctionalization following duplication, as well as impact from distinct domestication process. These genomic variations may underpin the differences observed in fruit and seed morphology, ripening, and disease resistance between cupuassu and the other Malvaceae species.

CONCLUSIONS

The cupuassu genome offers a foundational resource for both breeding improvement and conservation biology, yielding insights into the evolution and diversity within the genus Theobroma.

One hundred single-copy nuclear sequence markers for olive variety identification: a case of fingerprinting database construction in China

Olive is an ancient oil-producing tree, widely cultivated in Mediterranean countries, and now spread to other areas of the world, including China. Recently, several molecular databases were constructed in different countries and platforms for olive identification using simple sequence repeats (SSRs) or single-nucleotide polymorphisms (SNPs). However, comparing their results across laboratories was difficult. Herein, hundreds of polymorphic single-copy nuclear sequence markers were developed from the olive genome. Using the advantage of multiplex PCR amplification and high-throughput sequencing, a fingerprint database was constructed for the majority of olives cultivated in China. We used 100 high-quality sequence loci and estimated the genetic diversity and structure among all these varieties. We found that compared with that based on SSRs, the constructed fingerprint database based on these 100 sequences or a few of them, could provide a reliable olive variety identification platform in China, with high discrimination among different varieties using the principle of BLAST algorithm. An example of such identification platform based on this study was displayed on the web for the olive database in China (http://olivedb.cn/jianding). After resolving redundant genotypes, we identified 126 olive varieties with distinct genotypes in China. These varieties could be divided into two clusters, and it was revealed that the grouping of the varieties has a certain relationship with their origin. Herein, it is concluded that these single-copy orthologous nuclear sequences could be used to construct a universal fingerprint database of olives across different laboratories and platforms inexpensively. Based on such a database, variety identification can be performed easily by any laboratory, which would further facilitate olive breeding and variety exchange globally.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-023-01434-9.

Phylogenomic and syntenic data demonstrate complex evolutionary processes in early radiation of the rosids

Some of the most vexing problems of deep level relationship that remain in angiosperms involve the superrosids. The superrosid clade contains a quarter of all angiosperm species, with 18 orders in three subclades (Vitales, Saxifragales and core rosids) exhibiting remarkable morphological and ecological diversity. To help resolve deep-level relationships, we constructed a high-quality chromosome-level genome assembly for Tiarella polyphylla (Saxifragaceae) thus providing broader genomic representation of Saxifragales. Whole genome microsynteny analysis of superrosids showed that Saxifragales shared more synteny clusters with core rosids than Vitales, further supporting Saxifragales as more closely related with core rosids. To resolve the ordinal phylogeny of superrosids, we screened 122 single copy nuclear genes from genomes of 36 species, representing all 18 superrosid orders. Vitales were recovered as sister to all other superrosids (Saxifragales + core rosids). Our data suggest dramatic differences in relationships compared to earlier studies within core rosids. Fabids should be restricted to the nitrogen-fixing clade, while Picramniales, the Celastrales-Malpighiales (CM) clade, Huerteales, Oxalidales, Sapindales, Malvales and Brassicales formed an "expanded" malvid clade. The Celastrales-Oxalidales-Malpighiales (COM) clade (sensu APG IV) was not monophyletic. Crossosomatales, Geraniales, Myrtales and Zygophyllales did not belong to either of our well-supported malvids or fabids. There is strong discordance between nuclear and plastid phylogenetic hypotheses for superrosid relationships; we show that this is best explained by a combination of incomplete lineage sorting and ancient reticulation.

Functional divergence of CYP76AKs shapes the chemodiversity of abietane-type diterpenoids in genus Salvia

The genus Salvia L. (Lamiaceae) comprises myriad distinct medicinal herbs, with terpenoids as one of their major active chemical groups. Abietane-type diterpenoids (ATDs), such as tanshinones and carnosic acids, are specific to Salvia and exhibit taxonomic chemical diversity among lineages. To elucidate how ATD chemical diversity evolved, we carried out large-scale metabolic and phylogenetic analyses of 71 Salvia species, combined with enzyme function, ancestral sequence and chemical trait reconstruction, and comparative genomics experiments. This integrated approach showed that the lineage-wide ATD diversities in Salvia were induced by differences in the oxidation of the terpenoid skeleton at C-20, which was caused by the functional divergence of the cytochrome P450 subfamily CYP76AK. These findings present a unique pattern of chemical diversity in plants that was shaped by the loss of enzyme activity and associated catalytic pathways.

The genome of Acorus deciphers insights into early monocot evolution

Acorales is the sister lineage to all the other extant monocot plants. Genomic resource enhancement of this genus can help to reveal early monocot genomic architecture and evolution. Here, we assemble the genome of Acorus gramineus and reveal that it has ~45% fewer genes than the majority of monocots, although they have similar genome size. Phylogenetic analyses based on both chloroplast and nuclear genes consistently support that A. gramineus is the sister to the remaining monocots. In addition, we assemble a 2.2 Mb mitochondrial genome and observe many genes exhibit higher mutation rates than that of most angiosperms, which could be the reason leading to the controversies of nuclear genes- and mitochondrial genes-based phylogenetic trees existing in the literature. Further, Acorales did not experience tau (τ) whole-genome duplication, unlike majority of monocot clades, and no large-scale gene expansion is observed. Moreover, we identify gene contractions and expansions likely linking to plant architecture, stress resistance, light harvesting, and essential oil metabolism. These findings shed light on the evolution of early monocots and genomic footprints of wetland plant adaptations.

Transcriptome-Based Study on the Phylogeny and Hybridization of Marattialean Ferns (Marattiaceae)

Marattiaceae is a phylogenetically isolated family of tropical eusporangiate ferns including six genera with more than one-hundred species. In Marattiaceae, monophyly of genera has been well-supported phylogenetically. However, the phylogenetic relationships among them were elusive and controversial. Here, a dataset of 26 transcriptomes (including 11 newly generated) were used to assess single-copy nuclear genes and to obtain the organelle gene sequences. Through phylotranscriptomic analysis, the phylogeny and hybridization events of Marattiaceae were explored and a robust phylogenomic framework for the evolution of Marattiaceae was provided. Using both concatenation- and coalescent-based phylogenies, the gene-tree discordance, incomplete lineage sorting (ILS) simulations, and network inference were examined. Except the low support with mitochondrial genes of Marattiaceae, nuclear genes and chloroplast genes strongly supported a sister relationship between Marattiaceae and leptosporangiate ferns. At the genus level, all phylogenetic analysis based on nuclear genes datasets recovered five genera in Marattiaceae as monophyletic with strong support. Danaea and Ptisana were the first two diverged clades in turn. Christensenia was a sister clade to the clade Marattia + Angiopteris s.l. In Angiopteris s.l., three clades (Angiopteris s.s., the Archangiopteris group, and An. sparsisora) were well identified with maximum support. The Archangiopteris group was derived from Angiopteris s.s. at ca. 18 Ma. The putative hybrid species An. sparsisora between Angiopteris s.s. and the Archangiopteris group was verified by the species network analyses and the maternal plastid genes. This study will improve our understanding for using the phylotranscriptomic method to explore phylogeny and investigate hybridization events for difficult taxa in ferns.

Polyploidy and diploidization in soybean

Polyploidy is widespread and particularly common in angiosperms. The prevalence of polyploidy in the plant suggests it as a crucial driver of diversification and speciation. The paleopolyploid soybean (Glycine max) is one of the most important crops of plant protein and oil for humans and livestock. Soybean experienced two rounds of whole genome duplication around 13 and 59 million years ago. Due to the relatively slow process of post-polyploid diploidization, most genes are present in multiple copies across the soybean genome. Growing evidence suggests that polyploidization and diploidization could cause rapid and dramatic changes in genomic structure and epigenetic modifications, including gene loss, transposon amplification, and reorganization of chromatin architecture. This review is focused on recent progresses about genetic and epigenetic changes during polyploidization and diploidization of soybean and represents the challenges and potentials for application of polyploidy in soybean breeding.

Deep reticulation: the long legacy of hybridization in vascular plant evolution

Hybridization has long been recognized as a fundamental evolutionary process in plants but, until recently, our understanding of its phylogenetic distribution and biological significance across deep evolutionary scales has been largely obscure. Over the past decade, genomic and phylogenomic datasets have revealed, perhaps not surprisingly, that hybridization, often associated with polyploidy, has been common throughout the evolutionary history of plants, particularly in various lineages of flowering plants. However, phylogenomic studies have also highlighted the challenges of disentangling signals of ancient hybridization from other sources of genomic conflict (in particular, incomplete lineage sorting). Here, we provide a critical review of ancient hybridization in vascular plants, outlining well-documented cases of ancient hybridization across plant phylogeny, as well as the challenges unique to documenting ancient versus recent hybridization. We provide a definition for ancient hybridization, which, to our knowledge, has not been explicitly attempted before. Further documenting the extent of deep reticulation in plants should remain an important research focus, especially because published examples likely represent the tip of the iceberg in terms of the total extent of ancient hybridization. However, future research should increasingly explore the macroevolutionary significance of this process, in terms of its impact on evolutionary trajectories (e.g. how does hybridization influence trait evolution or the generation of biodiversity over long time scales?), as well as how life history and ecological factors shape, or have shaped, the frequency of hybridization across geologic time and plant phylogeny. Finally, we consider the implications of ubiquitous ancient hybridization for how we conceptualize, analyze, and classify plant phylogeny. Networks, as opposed to bifurcating trees, represent more accurate representations of evolutionary history in many cases, although our ability to infer, visualize, and use networks for comparative analyses is highly limited. Developing improved methods for the generation, visualization, and use of networks represents a critical future direction for plant biology. Current classification systems also do not generally allow for the recognition of reticulate lineages, and our classifications themselves are largely based on evidence from the chloroplast genome. Updating plant classification to better reflect nuclear phylogenies, as well as considering whether and how to recognize hybridization in classification systems, will represent an important challenge for the plant systematics community.

Comparative chloroplast genomics of 34 species in subtribe Swertiinae (Gentianaceae) with implications for its phylogeny

BACKGROUND

Subtribe Swertiinae, a medicinally significant and highly speciose Subtribe of family Gentianaceae. Despite previous extensive studies based on both morphology and molecular data, intergeneric and infrageneric relationships within subtribe Swertiinae remain controversial.

METHODS

Here, we employed four newly generated Swertia chloroplast genomes with thirty other published genomes to elucidate their genomic characteristics.

RESULTS

The 34 chloroplast genomes were small and ranged in size from 149,036 to 154,365 bp, each comprising two inverted repeat regions (size range 25,069-26,126 bp) that separated large single-copy (80,432-84,153 bp) and small single-copy (17,887-18,47 bp) regions, and all the chloroplast genomes showed similar gene orders, contents, and structures. These chloroplast genomes contained 129-134 genes each, including 84-89 protein-coding genes, 37 tRNAs, and 8 rRNAs. The chloroplast genomes of subtribe Swertiinae appeared to have lost some genes, such as rpl33, rpl2 and ycf15 genes. Comparative analyses revealed that two mutation hotspot regions (accD-psaI and ycf1) could serve as effective molecular markers for further phylogenetic analyses and species identification in subtribe Swertiinae. Positive selection analyses showed that two genes (ccsA and psbB) had high Ka/Ks ratios, indicating that chloroplast genes may have undergone positive selection in their evolutionary history. Phylogenetic analysis showed that the 34 subtribe Swertiinae species formed a monophyletic clade, with Veratrilla, Gentianopsis and Pterygocalyx located at the base of the phylogenetic tree. Some genera of this subtribe, however, were not monophyletic, including Swertia, Gentianopsis, Lomatogonium, Halenia, Veratrilla and Gentianopsis. In addition, our molecular phylogeny was consistent with taxonomic classification of subtribe Swertiinae in the Roate group and Tubular group. The results of molecular dating showed that the divergence between subtrib Gentianinae and subtrib Swertiinae was estimated to occur in 33.68 Ma. Roate group and Tubular group in subtribe Swertiinae approximately diverged in 25.17 Ma.

CONCLUSION

Overall, our study highlighted the taxonomic utility of chloroplast genomes in subtribe Swertiinae, and the genetic markers identified here will facilitate future studies on the evolution, conservation, population genetics, and phylogeography of subtribe Swertiinae species.

Comparative transcriptome analyses of three Gentiana species provides signals for the molecular footprints of selection effects and the phylogenetic relationships

Species in Gentiana section Cruciata are important alpine plants with a center of diversity and speciation in Qinghai-Tibet Plateau (QTP), and some of these species are sympatrically distributed in northeastern QTP. Studies on genome features and natural selection signatures of sympatric species in section Crucata have been impeded by a lack of genomic resources. Here, we showed transcript characterizations and molecular footprints of selection effects on G. straminea, G. dahurica and G. officinalis based on the comparative transcriptome. A total of 62.97 Gb clean reads were obtained with unigene numbers per species ranging from 141,819 to 236,408 after assembly. We found that these three species had similar distribution of functional categories in different databases, and key enzyme-encoding genes involved in the iridoids biosynthesis were also obtained. The selective pressure analyses indicated that most paired orthologs between these three species were subject to negative selection, and only a low proportion of the orthologs that underwent positive selection were detected. We found that some positive selected genes were involved in "catalytic activity", "metabolic process", "response to stimulus" and "response to stress". Besides, large numbers of SSR primer pairs with transferabilities were successfully designed based on the available transcriptome datasets of three Gentiana species. The phylogenetic relationships reconstructed based on 352 single-copy nuclear genes provided a rough phylogenetic framework for this genus and confirmed the monophyly of section Cruciata. Our study not only provides insights for the natural selection effects on sympatric Gentiana species, but also enhances future genetic breeding or evolutionary studies on Qinjiao species.

Phylotranscriptomic Analyses of Mycoheterotrophic Monocots Show a Continuum of Convergent Evolutionary Changes in Expressed Nuclear Genes From Three Independent Nonphotosynthetic Lineages

Mycoheterotrophy is an alternative nutritional strategy whereby plants obtain sugars and other nutrients from soil fungi. Mycoheterotrophy and associated loss of photosynthesis have evolved repeatedly in plants, particularly in monocots. Although reductive evolution of plastomes in mycoheterotrophs is well documented, the dynamics of nuclear genome evolution remains largely unknown. Transcriptome datasets were generated from four mycoheterotrophs in three families (Orchidaceae, Burmanniaceae, Triuridaceae) and related green plants and used for phylogenomic analyses to resolve relationships among the mycoheterotrophs, their relatives, and representatives across the monocots. Phylogenetic trees based on 602 genes were mostly congruent with plastome phylogenies, except for an Asparagales + Liliales clade inferred in the nuclear trees. Reduction and loss of chlorophyll synthesis and photosynthetic gene expression and relaxation of purifying selection on retained genes were progressive, with greater loss in older nonphotosynthetic lineages. One hundred seventy-four of 1375 plant benchmark universally conserved orthologous genes were undetected in any mycoheterotroph transcriptome or the genome of the mycoheterotrophic orchid Gastrodia but were expressed in green relatives, providing evidence for massively convergent gene loss in nonphotosynthetic lineages. We designate this set of deleted or undetected genes Missing in Mycoheterotrophs (MIM). MIM genes encode not only mainly photosynthetic or plastid membrane proteins but also a diverse set of plastid processes, genes of unknown function, mitochondrial, and cellular processes. Transcription of a photosystem II gene (psb29) in all lineages implies a nonphotosynthetic function for this and other genes retained in mycoheterotrophs. Nonphotosynthetic plants enable novel insights into gene function as well as gene expression shifts, gene loss, and convergence in nuclear genomes.

Genome assembly of an Australian native grass species reveals a recent whole-genome duplication and biased gene retention of genes involved in stress response

BACKGROUND

The adaptive significance of polyploidy has been extensively debated, and chromosome-level genome assemblies of polyploids can provide insight into this. The Australian grass Bothriochloa decipiens belongs to the BCD clade, a group with a complex history of hybridization and polyploid. This is the first genome assembly and annotation of a species that belongs to this fascinating yet complex group.

FINDINGS

Using Illumina short reads, 10X Genomics linked reads, and Hi-C sequencing data, we assembled a highly contiguous genome of B. decipiens, with a total length of 1,218.22 Mb and scaffold N50 of 42.637 Mb. Comparative analysis revealed that the species experienced a relatively recent whole-genome duplication. We clustered the 20 major scaffolds, representing the 20 chromosomes, into the 2 subgenomes of the parental species using unique repeat signatures. We found evidence of biased fractionation and differences in the activity of transposable elements between the subgenomes prior to hybridization. Duplicates were enriched for genes involved in transcription and response to external stimuli, supporting a biased retention of duplicated genes following whole-genome duplication.

CONCLUSIONS

Our results support the hypotheses of a biased retention of duplicated genes following polyploidy and point to differences in repeat activity associated with subgenome dominance. B. decipiens is a widespread species with the ability to establish across many soil types, making it a prime candidate for climate change- resilient ecological restoration of Australian grasslands. This reference genome is a valuable resource for future population genomic research on Australian grasses.

A pipeline for assembling low copy nuclear markers from plant genome skimming data for phylogenetic use

Background

Genome skimming is a popular method in plant phylogenomics that do not include a biased enrichment step, relying on random shallow sequencing of total genomic DNA. From these data the plastome is usually readily assembled and constitutes the bulk of phylogenetic information generated in these studies. Despite a few attempts to use genome skims to recover low copy nuclear loci for direct phylogenetic use, such endeavor remains neglected. Causes might include the trade-off between libraries with few reads and species with large genomes (i.e., missing data caused by low coverage), but also might relate to the lack of pipelines for data assembling.

Methods

A pipeline and its companion R package designed to automate the recovery of low copy nuclear markers from genome skimming libraries are presented. Additionally, a series of analyses aiming to evaluate the impact of key assembling parameters, reference selection and missing data are presented.

Results

A substantial amount of putative low copy nuclear loci was assembled and proved useful to base phylogenetic inference across the libraries tested (4 to 11 times more data than previously assembled plastomes from the same libraries).

Discussion

Critical aspects of assembling low copy nuclear markers from genome skims include the minimum coverage and depth of a sequence to be used. More stringent values of these parameters reduces the amount of assembled data and increases the relative amount of missing data, which can compromise phylogenetic inference, in turn relaxing the same parameters might increase sequence error. These issues are discussed in the text, and parameter tuning through multiple comparisons tracking their effects on support and congruence is highly recommended when using this pipeline. The skimmingLoci pipeline (https://github.com/mreginato/skimmingLoci) might stimulate the use of genome skims to recover nuclear loci for direct phylogenetic use, increasing the power of genome skimming data to resolve phylogenetic relationships, while reducing the amount of sequenced DNA that is commonly wasted.

Cytonuclear coevolution in a holoparasitic plant with highly disparate organellar genomes

Contrasting substitution rates in the organellar genomes of Lophophytum agree with the DNA repair, replication, and recombination gene content. Plastid and nuclear genes whose products form multisubunit complexes co-evolve. The organellar genomes of the holoparasitic plant Lophophytum (Balanophoraceae) show disparate evolution. In the plastid, the genome has been severely reduced and presents a > 85% AT content, while in the mitochondria most protein-coding genes have been replaced by homologs acquired by horizontal gene transfer (HGT) from their hosts (Fabaceae). Both genomes carry genes whose products form multisubunit complexes with those of nuclear genes, creating a possible hotspot of cytonuclear coevolution. In this study, we assessed the evolutionary rates of plastid, mitochondrial and nuclear genes, and their impact on cytonuclear evolution of genes involved in multisubunit complexes related to lipid biosynthesis and proteolysis in the plastid and those in charge of the oxidative phosphorylation in the mitochondria. Genes from the plastid and the mitochondria (both native and foreign) of Lophophytum showed extremely high and ordinary substitution rates, respectively. These results agree with the biased loss of plastid-targeted proteins involved in angiosperm organellar repair, replication, and recombination machinery. Consistent with the high rate of evolution of plastid genes, nuclear-encoded subunits of plastid complexes showed disproportionate increases in non-synonymous substitution rates, while those of the mitochondrial complexes did not show different rates than the control (i.e. non-organellar nuclear genes). Moreover, the increases in the nuclear-encoded subunits of plastid complexes were positively correlated with the level of physical interaction they possess with the plastid-encoded ones. Overall, these results suggest that a structurally-mediated compensatory factor may be driving plastid-nuclear coevolution in Lophophytum, and that mito-nuclear coevolution was not altered by HGT.

Genome-Wide Identification, Characterization, and Comparative Analysis of NLR Resistance Genes in Coffea spp

The largest family of disease resistance genes in plants are nucleotide-binding site leucine-rich repeat genes (NLRs). The products of these genes are responsible for recognizing avirulence proteins (Avr) of phytopathogens and triggering specific defense responses. Identifying NLRs in plant genomes with standard gene annotation software is challenging due to their multidomain nature, sequence diversity, and clustered genomic distribution. We present the results of a genome-wide scan and comparative analysis of NLR loci in three coffee species (Coffea canephora, Coffea eugenioides and their interspecific hybrid Coffea arabica). A total of 1311 non-redundant NLR loci were identified in C. arabica, 927 in C. canephora, and 1079 in C. eugenioides, of which 809, 562, and 695 are complete loci, respectively. The NLR-Annotator tool used in this study showed extremely high sensitivities and specificities (over 99%) and increased the detection of putative NLRs in the reference coffee genomes. The NLRs loci in coffee are distributed among all chromosomes and are organized mostly in clusters. The C. arabica genome presented a smaller number of NLR loci when compared to the sum of the parental genomes (C. canephora, and C. eugenioides). There are orthologous NLRs (orthogroups) shared between coffee, tomato, potato, and reference NLRs and those that are shared only among coffee species, which provides clues about the functionality and evolutionary history of these orthogroups. Phylogenetic analysis demonstrated orthologous NLRs shared between C. arabica and the parental genomes and those that were possibly lost. The NLR family members in coffee are subdivided into two main groups: TIR-NLR (TNL) and non-TNL. The non-TNLs seem to represent a repertoire of resistance genes that are important in coffee. These results will support functional studies and contribute to a more precise use of these genes for breeding disease-resistant coffee cultivars.

Phylogenomic analysis does not support a classic but controversial hypothesis of progenitor-derivative origins for the serpentine endemic Clarkia franciscana

Budding speciation involves isolation of marginal populations at the periphery of a species range and is thought to be a prominent mode of speciation in organisms with low dispersal and/or strong local adaptation among populations. Budding speciation is typically evidenced by abutting, asymmetric ranges of ecologically divergent sister species and low genetic diversity in putative budded species. Yet these indirect patterns may be unreliable, instead caused by postspeciation processes such as range or demographic shifts. Nested phylogenetic relationships provide the most conclusive evidence of budding speciation. A putative case of budding speciation in the serpentine endemic Clarkia franciscana and two closely related widespread congeners was studied by Harlan Lewis, Peter Raven, Leslie Gottlieb, and others over a 20-year period, yet the origin of C. franciscana remains controversial. Here, we reinvestigate this system with phylogenomic analyses to determine whether C. franciscana is a recently derived budded species, phylogenetically nested within one of the other two putative progenitor species. In contrast to the hypothesized pattern of relatedness among the three Clarkia species, we find no evidence for recent budding speciation. Instead, the data suggest the three species diverged simultaneously. We urge caution in using contemporary range patterns to infer geographic modes of speciation.

Target enrichment and extensive population sampling help untangle the recent, rapid radiation of Oenothera sect. Calylophus

Oenothera sect. Calylophus is a North American group of 13 recognized taxa in the evening primrose family (Onagraceae) with an evolutionary history that may include independent origins of bee pollination, edaphic endemism, and permanent translocation heterozygosity. Like other groups that radiated relatively recently and rapidly, taxon boundaries within Oenothera sect. Calylophus have remained challenging to circumscribe. In this study, we used target enrichment, flanking non-coding regions, gene tree/species tree methods, tests for gene flow modified for target-enrichment data, and morphometric analysis to reconstruct phylogenetic hypotheses, evaluate current taxon circumscriptions, and examine character evolution in Oenothera sect. Calylophus. Because sect. Calylophus comprises a clade with a relatively restricted geographic range, we were able to extensively sample across the range of geographic, edaphic and morphological diversity in the group. We found that the combination of exons and flanking non-coding regions led to improved support for species relationships. We reconstructed potential hybrid origins of some accessions and note that if processes such as hybridization are not taken into account, the number of inferred evolutionary transitions may be artificially inflated. We recovered strong evidence for multiple evolutionary origins of bee pollination from ancestral hawkmoth pollination, edaphic specialization on gypsum, and permanent translocation heterozygosity. This study applies newly emerging techniques alongside dense infraspecific sampling and morphological analyses to effectively reconstruct the recalcitrant history of a rapid radiation.

Phylotranscriptomics Resolves the Phylogeny of Pooideae and Uncovers Factors for Their Adaptive Evolution

Adaptation to cool climates has occurred several times in different angiosperm groups. Among them, Pooideae, the largest grass subfamily with ∼3,900 species including wheat and barley, have successfully occupied many temperate regions and play a prominent role in temperate ecosystems. To investigate possible factors contributing to Pooideae adaptive evolution to cooling climates, we performed phylogenetic reconstruction using five gene sets (with 1,234 nuclear genes and their subsets) from 157 transcriptomes/genomes representing all 15 tribes and 24 of 26 subtribes. Our phylogeny supports the monophyly of all tribes (except Diarrheneae) and all subtribes with at least two species, with strongly supported resolution of their relationships. Molecular dating suggests that Pooideae originated in the late Cretaceous, with subsequent divergences under cooling conditions first among many tribes from the early middle to late Eocene and again among genera in the middle Miocene and later periods. We identified a cluster of gene duplications (CGD5) shared by the core Pooideae (with 80% Pooideae species) near the Eocene–Oligocene transition, coinciding with the transition from closed to open habitat and an upshift of diversification rate. Molecular evolutionary analyses homologs of CBF for cold resistance uncovered tandem duplications during the core Pooideae history, dramatically increasing their copy number and possibly promoting adaptation to cold habitats. Moreover, duplication of AP1/FUL-like genes before the Pooideae origin might have facilitated the regulation of the vernalization pathway under cold environments. These and other results provide new insights into factors that likely have contributed to the successful adaptation of Pooideae members to temperate regions.

Comparative Genome Analysis Reveals Phylogenetic Identity of Bacillus velezensis HNA3 and Genomic Insights into Its Plant Growth Promotion and Biocontrol Effects

Bacillus velezensis HNA3, a potential plant growth promoter and biocontrol rhizobacterium, was isolated from plant rhizosphere soils in our previous work. Here, we sequenced the entire genome of the HNA3 strain and performed a comparative genome analysis. We found that HNA3 has a 3,929-kb chromosome with 46.5% GC content and 4,080 CDSs. We reclassified HNA3 as a Bacillus velezensis strain by core genome analysis between HNA3 and 74 previously defined Bacillus strains in the evolutionary tree. A comparative genomic analysis among Bacillus velezensis HNA3, Bacillus velezensis FZB42, Bacillus amyloliquefaciens DSM7, and Bacillus subtilis 168 showed that only HNA3 has one predicated secretory protein feruloyl esterase that catalyzes the hydrolysis of plant cell wall polysaccharides. The analysis of gene clusters revealed that whole biosynthetic gene clusters type Lanthipeptide was exclusively identified in HNA3 and might lead to the synthesis of new bioactive compounds. Twelve gene clusters were detected in HNA3 responsible for the synthesis of 14 secondary metabolites including Bacillaene, Fengycin, Bacillomycin D, Surfactin, Plipastatin, Mycosubtilin, Paenilarvins, Macrolactin, Difficidin, Amylocyclicin, Bacilysin, Iturin, Bacillibactin, Paenibactin, and others. HNA3 has 77 genes encoding for possible antifungal and antibacterial secreting carbohydrate active enzymes. It also contains genes involved in plant growth promotion, such as 11 putative indole acetic acid (IAA)-producing genes, spermidine and polyamine synthase genes, volatile compound producing genes, and multiple biofilm related genes. HNA3 also has 19 phosphatase genes involved in phosphorus solubilization. Our results provide insights into the genetic characteristics responsible for the bioactivities and potential application of HNA3 as plant growth-promoting strain in ecological agriculture.

IMPORTANCE This study is the primary initiative to identify Bacillus velezensis HNA3 whole genome sequence and reveal its genomic properties as an effective biocontrol agent against plant pathogens and a plant growth stimulator. HNA3 genetic profile can be used as a reference for future studies that can be applied as a highly effective biofertilizer and biofungicide inoculum to improve agriculture productivity. HNA3 reclassified in the phylogenetic tree which may be helpful for highly effective strain engineering and taxonomy. The genetic comparison among HNA3 and closely similar species B. velezensis FZB42, B. amyloliquefaciens DSM7, and B. subtilis 168 demonstrates some distinctive genetic properties of HNA3 and provides a basis for the genetic diversity of the Bacillus genus, which allows developing more effective eco-friendly resources for agriculture and separation of Bacillus velezensis as distinct species in the phylogenetic tree.

Perspectives in Myrtaceae evolution from plastomes and nuclear phylogenies

Myrtaceae is a large and species-rich family of woody eudicots, with prevalent distribution in the Southern Hemisphere. Classification and taxonomy of species belonging to this family is quite challenging, sometimes with difficulty in species identification and producing phylogenies with low support for species relationships. Most of the current knowledge comes from few molecular markers, such as plastid genes and intergenic regions, which can be difficult to handle and produce conflicting results. Based on plastid protein-coding sequences and nuclear markers, we present a topology for the phylogenetic relationships among Myrtaceae tribes. Our phylogenetic estimate offers a contrasting topology over previous analysis with fewer markers. Plastome phylogeny groups the tribes Syzygieae and Eucalypteae and individual chloroplast genes produce divergent topologies, especially among species within Myrteae tribe, but also in regard to the grouping of Syzygieae and Eucalypteae. Results are consistent and reproducible with both nuclear and organellar datasets. It confronts previous data about the deep nodes of Myrtaceae phylogeny.

When tropical and subtropical congeners met: Multiple ancient hybridization events within Eriobotrya in the Yunnan-Guizhou Plateau, a tropical-subtropical transition area in China

Global climate changes during the Miocene may have created ample opportunities for hybridization between members of tropical and subtropical biomes at the boundary between these zones. Yet, very few studies have explored this possibility. The Yunnan-Guizhou Plateau (YGP) in Southwest China is a biodiversity hotspot for vascular plants, located in a transitional area between the floristic regions of tropical Southeast Asia and subtropical East Asia. The genus Eriobotrya (Rosaceae) comprises both tropical and subtropical taxa, with 12 species recorded in the YGP, making it a suitable basis for testing the hypothesis of between-biome hybridization. Therefore, we surveyed the evolutionary history of Eriobotrya by examining three chloroplast regions and five nuclear genes for 817 individuals (47 populations) of 23 Eriobotrya species (including 19 populations of 12 species in the YGP), plus genome re-sequencing of 33 representative samples. We concluded that: (1) phylogenetic positions for 16 species exhibited strong cytonuclear conflicts, most probably due to ancient hybridization; (2) the YGP is a hotspot for hybridization, with 11 species showing clear evidence of chloroplast capture; and (3) Eriobotrya probably originated in tropical Asia during the Eocene. From the Miocene onwards, the intensification of the Eastern Asia monsoon and global cooling may have shifted the tropical-subtropical boundary and caused secondary contact between species, thus providing ample opportunity for hybridization and diversification of Eriobotrya, especially in the YGP. Our study highlights the significant role that paleoclimate changes probably played in driving hybridization and generating rich species diversity in climate transition zones.

Andean uplift, drainage basin formation, and the evolution of plants living in fast-flowing aquatic ecosystems in northern South America

Northern South America is a geologically dynamic and species-rich region. Fossil and stratigraphic data show that mountain uplift in the tropical Andes reconfigured river drainages. These landscape changes shaped the evolution of the flora in the region, yet the impacts on aquatic taxa have been overlooked. We explore the role of landscape change on the evolution of plants living strictly in rivers across drainage basins in northern South America by conducting population structure, phylogenetic inference, and divergence-dating analyses for two species in the genus Marathrum (Podostemaceae). Mountain uplift and drainage basin formation isolated populations of M. utile and M. foeniculaceum in northern South America and created barriers to gene flow across river drainages. Sympatric species hybridize and the hybrids show the phenotype of one parental line. We propose that the pattern of divergence of populations reflects the formation of river drainages, which was not complete until < 4.1 million yr ago (Ma). Our study provides a clear picture of the role of landscape change on the evolution of plants living strictly in rivers in northern South America. By shifting the focus to aquatic taxa, we provide a novel perspective on the processes shaping the evolution of the Neotropical flora.

Comparative transcriptomics provides a strategy for phylogenetic analysis and SSR marker development in Chaenomeles

The genus Chaenomeles has long been considered an important ornamental, herbal and cash crop and is widely cultivated in East Asia. Traditional studies of Chaenomeles mainly focus on evolutionary relationships at the phenotypic level. In this study, we conducted RNA-seq on 10 Chaenomeles germplasms supplemented with one outgroup species, Docynia delavayi (D. delavayi), on the Illumina HiSeq2500 platform. After de novo assemblies, we generated from 40,084 to 49,571 unigenes for each germplasm. After pairwise comparison of the orthologous sequences, 9,659 orthologues within the 11 germplasms were obtained, with 6,154 orthologous genes identified as single-copy genes. The phylogenetic tree was visualized to reveal evolutionary relationships for these 11 germplasms. GO and KEGG analyses were performed for these common single-copy genes to compare their functional similarities and differences. Selective pressure analysis based on 6,154 common single-copy genes revealed that 45 genes were under positive selection. Most of these genes are involved in building the plant disease defence system. A total of 292 genes containing simple sequence repeats (SSRs) were used to develop SSR markers and compare their functions in secondary metabolism pathways. Finally, 10 primers were chosen as SSR marker candidates for Chaenomeles germplasms by comprehensive standards. Our research provides a new methodology and reference for future related research in Chaenomeles and is also useful for improvement, breeding and selection projects in other related species.

Host-associated microbe PCR (hamPCR) enables convenient measurement of both microbial load and community composition

The ratio of microbial population size relative to the amount of host tissue, or 'microbial load', is a fundamental metric of colonization and infection, but it cannot be directly deduced from microbial amplicon data such as 16S rRNA gene counts. Because existing methods to determine load, such as serial dilution plating, quantitative PCR, and whole metagenome sequencing add substantial cost and/or experimental burden, they are only rarely paired with amplicon sequencing. We introduce host-associated microbe PCR (hamPCR), a robust strategy to both quantify microbial load and describe interkingdom microbial community composition in a single amplicon library. We demonstrate its accuracy across multiple study systems, including nematodes and major crops, and further present a cost-saving technique to reduce host overrepresentation in the library prior to sequencing. Because hamPCR provides an accessible experimental solution to the well-known limitations and statistical challenges of compositional data, it has far-reaching potential in culture-independent microbiology.

On the potential of Angiosperms353 for population genomic studies

PREMISE

The successful application of universal targeted sequencing markers, such as those developed for the Angiosperms353 probe set, within populations could reduce or eliminate the need for specific marker development, while retaining the benefits of full-gene sequences in population-level analyses. However, whether the Angiosperms353 markers provide sufficient variation within species to calculate demographic parameters is untested.

METHODS

Using herbarium specimens from a 50-year-old floristic survey in Texas, we sequenced 95 samples from 24 species using the Angiosperms353 probe set. Our data workflow calls variants within species and prepares data for population genetic analysis using standard metrics. In our case study, gene recovery was affected by genomic library concentration only at low concentrations and displayed limited phylogenetic bias.

RESULTS

We identified over 1000 segregating variants with zero missing data for 92% of species and demonstrate that Angiosperms353 markers contain sufficient variation to estimate pairwise nucleotide diversity (π)-typically between 0.002 and 0.010, with most variation found in flanking non-coding regions. In a subset of variants that were filtered to reduce linkage, we uncovered high heterozygosity in many species, suggesting that denser sampling within species should permit estimation of gene flow and population dynamics.

DISCUSSION

Angiosperms353 should benefit conservation genetic studies by providing universal repeatable markers, low missing data, and haplotype information, while permitting inclusion of decades-old herbarium specimens.

HAG1 and SWI3A/B control of male germ line development in P. patens suggests conservation of epigenetic reproductive control across land plants

KEY MESSAGE

Bryophytes as models to study the male germ line: loss-of-function mutants of epigenetic regulators HAG1 and SWI3a/b demonstrate conserved function in sexual reproduction. With the water-to-land transition, land plants evolved a peculiar haplodiplontic life cycle in which both the haploid gametophyte and the diploid sporophyte are multicellular. The switch between these phases was coined alternation of generations. Several key regulators that control the bauplan of either generation are already known. Analyses of such regulators in flowering plants are difficult due to the highly reduced gametophytic generation, and the fact that loss of function of such genes often is embryo lethal in homozygous plants. Here we set out to determine gene function and conservation via studies in bryophytes. Bryophytes are sister to vascular plants and hence allow evolutionary inferences. Moreover, embryo lethal mutants can be grown and vegetatively propagated due to the dominance of the bryophyte gametophytic generation. We determined candidates by selecting single copy orthologs that are involved in transcriptional control, and of which flowering plant mutants show defects during sexual reproduction, with a focus on the under-studied male germ line. We selected two orthologs, SWI3a/b and HAG1, and analyzed loss-of-function mutants in the moss P. patens. In both mutants, due to lack of fertile spermatozoids, fertilization and hence the switch to the diploid generation do not occur. Pphag1 additionally shows arrested male and impaired female gametangia development. We analyzed HAG1 in the dioecious liverwort M. polymorpha and found that in Mphag1 the development of gametangiophores is impaired. Taken together, we find that involvement of both regulators in sexual reproduction is conserved since the earliest divergence of land plants.

Genome-wide signatures of plastid-nuclear coevolution point to repeated perturbations of plastid proteostasis systems across angiosperms

Nuclear and plastid (chloroplast) genomes experience different mutation rates, levels of selection, and transmission modes, yet key cellular functions depend on their coordinated interactions. Functionally related proteins often show correlated changes in rates of sequence evolution across a phylogeny [evolutionary rate covariation (ERC)], offering a means to detect previously unidentified suites of coevolving and cofunctional genes. We performed phylogenomic analyses across angiosperm diversity, scanning the nuclear genome for genes that exhibit ERC with plastid genes. As expected, the strongest hits were highly enriched for genes encoding plastid-targeted proteins, providing evidence that cytonuclear interactions affect rates of molecular evolution at genome-wide scales. Many identified nuclear genes functioned in post-transcriptional regulation and the maintenance of protein homeostasis (proteostasis), including protein translation (in both the plastid and cytosol), import, quality control, and turnover. We also identified nuclear genes that exhibit strong signatures of coevolution with the plastid genome, but their encoded proteins lack organellar-targeting annotations, making them candidates for having previously undescribed roles in plastids. In sum, our genome-wide analyses reveal that plastid-nuclear coevolution extends beyond the intimate molecular interactions within chloroplast enzyme complexes and may be driven by frequent rewiring of the machinery responsible for maintenance of plastid proteostasis in angiosperms.

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.

Gene Duplication and Gain in the Trematode Atriophallophorus winterbourni Contributes to Adaptation to Parasitism

Gene duplications and novel genes have been shown to play a major role in helminth adaptation to a parasitic lifestyle because they provide the novelty necessary for adaptation to a changing environment, such as living in multiple hosts. Here we present the de novo sequenced and annotated genome of the parasitic trematode Atriophallophorus winterbourni and its comparative genomic analysis to other major parasitic trematodes. First, we reconstructed the species phylogeny, and dated the split of A. winterbourni from the Opisthorchiata suborder to approximately 237.4 Ma (±120.4 Myr). We then addressed the question of which expanded gene families and gained genes are potentially involved in adaptation to parasitism. To do this, we used hierarchical orthologous groups to reconstruct three ancestral genomes on the phylogeny leading to A. winterbourni and performed a GO (Gene Ontology) enrichment analysis of the gene composition of each ancestral genome, allowing us to characterize the subsequent genomic changes. Out of the 11,499 genes in the A. winterbourni genome, as much as 24% have arisen through duplication events since the speciation of A. winterbourni from the Opisthorchiata, and as much as 31.9% appear to be novel, that is, newly acquired. We found 13 gene families in A. winterbourni to have had more than ten genes arising through these recent duplications; all of which have functions potentially relating to host behavioral manipulation, host tissue penetration, and hiding from host immunity through antigen presentation. We identified several families with genes evolving under positive selection. Our results provide a valuable resource for future studies on the genomic basis of adaptation to parasitism and point to specific candidate genes putatively involved in antagonistic host-parasite adaptation.

Orthology and synteny analysis of receptor-like kinases "RLK" and receptor-like proteins "RLP" in legumes

BACKGROUND

Legume species are an important plant model because of their protein-rich physiology. The adaptability and productivity of legumes are limited by major biotic and abiotic stresses. Responses to these stresses directly involve plasma membrane receptor proteins known as receptor-like kinases and receptor-like proteins. Evaluating the homology relations among RLK and RLP for seven legume species, and exploring their presence among synteny blocks allow an increased understanding of evolutionary relations, physical position, and chromosomal distribution in related species and their shared roles in stress responses.

RESULTS

Typically, a high proportion of RLK and RLP legume proteins belong to orthologous clusters, which is confirmed in this study, where between 66 to 90% of the RLKs and RLPs per legume species were classified in orthologous clusters. One-third of the evaluated syntenic blocks had shared RLK/RLP genes among both legumes and non-legumes. Among the legumes, between 75 and 98% of the RLK/RLP were present in syntenic blocks. The distribution of chromosomal segments between Phaseolus vulgaris and Vigna unguiculata, two species that diverged ~ 8 mya, were highly similar. Among the RLK/RLP synteny clusters, seven experimentally validated resistance RLK/RLP genes were identified in syntenic blocks. The RLK resistant genes FLS2, BIR2, ERECTA, IOS1, and AtSERK1 from Arabidopsis and SLSERK1 from Solanum lycopersicum were present in different pairwise syntenic blocks among the legume species. Meanwhile, only the LYM1- RLP resistant gene from Arabidopsis shared a syntenic blocks with Glycine max.

CONCLUSIONS

The orthology analysis of the RLK and RLP suggests a dynamic evolution in the legume family, with between 66 to 85% of RLK and 83 to 88% of RLP belonging to orthologous clusters among the species evaluated. In fact, for the 10-species comparison, a lower number of singleton proteins were reported among RLP compared to RLK, suggesting that RLP positions are more physically conserved compared to RLK. The identification of RLK and RLP genes among the synteny blocks in legumes revealed multiple highly conserved syntenic blocks on multiple chromosomes. Additionally, the analysis suggests that P. vulgaris is an appropriate anchor species for comparative genomics among legumes.

Historical biogeography of Pomaderris (Rhamnaceae): Continental vicariance in Australia and repeated independent dispersals to New Zealand

AIM

Gondwanan biogeographic patterns include a combination of old vicariance events following the breakup of the supercontinent, and more recent long-distance dispersals across the southern landmasses. Floristic relationships between Australia and New Zealand have mostly been attributed to recent dispersal events rather than vicariance. We assessed the biogeographic history of Pomaderris (Rhamnaceae), which occurs in both Australia and New Zealand, by constructing a time-calibrated molecular phylogeny to infer (1) phylogenetic relationships and (2) the relative contributions of vicariance and dispersal events in the biogeographic history of the genus.

LOCATION

Australia and New Zealand.

METHODS

Using hybrid capture and high throughput sequencing, we generated nuclear and plastid data sets to estimate phylogenetic relationships and fossil calibrated divergence time estimates for Pomaderris. BioGeoBEARS and biogeographical stochastic mapping (BSM) were used to assess the ancestral area of the genus and the relative contributions of vicariance vs dispersal, and the directionality of dispersal events.

RESULTS

Our analyses indicate that Pomaderris originated in the Oligocene and had a widespread Australian distribution. Vicariance of western and eastern Australian clades coincides with the uplift of the Nullarbor Plain c. 14 Ma, followed by subsequent in-situ and within-biome diversification with little exchange across regions. A rapid radiation of southeastern Australian taxa beginning c. 10 Ma was the source for at least six independent long-distance dispersal events to New Zealand during the Pliocene-Pleistocene.

MAIN CONCLUSIONS

Our study demonstrates the importance of dispersal in explaining not only the current cross-Tasman distributions of Pomaderris, but for the New Zealand flora more broadly. The pattern of multiple independent long-distance dispersal events for Pomaderris, without significant radiation within New Zealand, is congruent with other lowland plant groups, suggesting that this biome has a different evolutionary history compared with the younger alpine flora of New Zealand, which exhibits extensive radiations often following single long distance dispersal events.

Phylogenomic analyses based on genome-skimming data reveal cyto-nuclear discordance in the evolutionary history of Cotoneaster (Rosaceae)

As a consequence of hybridization, polyploidization, and apomixis, the genus Cotoneaster (Rosaceae) represents one of the most complicated and controversial lineages in Rosaceae, with ca. 370 species which have been classified into two subgenera and several sections, and is notorious for its taxonomic difficulty. The infrageneric relationships and taxonomy of Cotoneaster have remained poorly understood. Previous studies have focused mainly on natural hybridization involving only several species, and phylogeny based on very limited markers. In the present study, the sequences of complete chloroplast genomes and 204 low-copy nuclear genes of 72 accessions, representing 69 species as ingroups, were used to conduct the most comprehensive phylogenetic analysis so far for Cotoneaster. Based on the sequences of complete chloroplast genomes and many nuclear genes, our analyses yield two robust phylogenetic trees respectively. Chloroplast genome and nuclear data confidently resolved relationships of this genus into two major clades which largely supported current classification based on morphological evidence. However, conflicts between the chloroplast genome and low-copy nuclear phylogenies were observed in both the species level and clade level. Cyto-nuclear discordance in the phylogeny could be caused by frequent hybridization events and incomplete sorting lineage (ILS). In addition, our divergence-time analysis revealed an evolutionary radiation of the genus from late Miocene to date.

A chromosome-scale reference genome of trifoliate orange (Poncirus trifoliata) provides insights into disease resistance, cold tolerance and genome evolution in Citrus

Trifoliate orange (Poncirus trifoliata), a deciduous close relative of evergreen Citrus, has important traits for citrus production, including tolerance/resistance to citrus greening disease (Huanglongbing, HLB) and other major diseases, and cold tolerance. It has been one of the most important rootstocks, and one of the most valuable sources of resistance and tolerance genes for citrus. Here we present a high-quality, chromosome-scale genome assembly of P. trifoliata. The 264.9-Mb assembly contains nine chromosomal pseudomolecules with 25 538 protein-coding genes, covering 97.2% of the estimated gene space. Comparative analyses of P. trifoliata and nine Citrus genomes revealed 605 species-specific genes and six rapidly evolving gene families in the P. trifoliata genome. Poncirus trifoliata has evolved specific adaptation in the C-repeat/DREB binding factor (CBF)-dependent and CBF-independent cold signaling pathways to tolerate cold. We identified candidate genes within quantitative trait loci for HLB tolerance, and at the loci for resistance to citrus tristeza virus and citrus nematode. Genetic diversity analysis of Poncirus accessions and Poncirus/Citrus hybrids shows a narrow genetic base in the US germplasm collection, and points to the importance of collecting and preserving more natural genetic variation. Two phenotypically divergent Poncirus accessions are found to be clonally related, supporting a previous conjecture that dwarf Flying Dragon originated as a mutant of a non-dwarfing type. The high-quality genome reveals features and evolutionary insights of Poncirus, and it will serve as a valuable resource for genetic, genomic and molecular research and manipulation in citrus.

Evolution of Homeologous Gene Expression in Polyploid Wheat

Polyploidization has played a prominent role in the evolutionary history of plants. Two recent and sequential allopolyploidization events have resulted in the formation of wheat species with different ploidies, and which provide a model to study the effects of polyploidization on the evolution of gene expression. In this study, we identified differentially expressed genes (DEGs) between four BBAA tetraploid wheats of three different ploidy backgrounds. DEGs were found to be unevenly distributed among functional categories and duplication modes. We observed more DEGs in the extracted tetraploid wheat (ETW) than in natural tetraploid wheats (TD and TTR13) as compared to a synthetic tetraploid (AT2). Furthermore, DEGs showed higher Ka/Ks ratios than those that did not show expression changes (non-DEGs) between genotypes, indicating DEGs and non-DEGs experienced different selection pressures. For A-B homeolog pairs with DEGs, most of them had only one differentially expressed copy, however, when both copies of a homeolog pair were DEGs, the A and B copies were more likely to be regulated to the same direction. Our results suggest that both cis- and inter-subgenome trans-regulatory changes are important drivers in the evolution of homeologous gene expression in polyploid wheat, with ploidy playing a significant role in the process.

Biased Gene Retention in the Face of Introgression Obscures Species Relationships

Phylogenomic analyses are recovering previously hidden histories of hybridization, revealing the genomic consequences of these events on the architecture of extant genomes. We applied phylogenomic techniques and several complementary statistical tests to show that introgressive hybridization appears to have occurred between close relatives of Arabidopsis, resulting in cytonuclear discordance and impacting our understanding of species relationships in the group. The composition of introgressed and retained genes indicates that selection against incompatible cytonuclear and nuclear-nuclear interactions likely acted during introgression, whereas linkage also contributed to genome composition through the retention of ancient haplotype blocks. We also applied divergence-based tests to determine the species branching order and distinguish donor from recipient lineages. Surprisingly, these analyses suggest that cytonuclear discordance arose via extensive nuclear, rather than cytoplasmic, introgression. If true, this would mean that most of the nuclear genome was displaced during introgression whereas only a small proportion of native alleles were retained.

Molecular Evidence for Hybrid Origin and Phenotypic Variation of Rosa Section Chinenses

Rosa sect. Chinenses (Rosaceae) is an important parent of modern rose that is widely distributed throughout China and plays an important role in breeding and molecular biological research. R. sect. Chinenses has variable morphological traits and mixed germplasm. However, the taxonomic status and genetic background of sect. Chinenses varieties remain unclear. In this study, we collected germplasm resources from sect. Chinenses varieties with different morphological traits. Simple sequence repeat (SSR) markers, chloroplast markers, and single copy nuclear markers were used to explore the genetic background of these germplasm resources. We described the origin of hybridization of rose germplasm resources by combining different molecular markers. The results showed that the flower and hip traits of different species in R. sect. Chinenses were significantly different. The SSR analysis showed that the two wild type varieties have different genetic backgrounds. The double petal varieties of R. sect. Chinenses could be hybrids of two wild type varieties. A phylogenetic analysis showed that the maternal inheritance of sect. Chinenses varieties had two different origins. To some extent, variation in the morphological traits of double petal species of R. sect. Chinenses reflects the influence of cultivation process. This study emphasizes that different genetic markers vary in their characteristics. Therefore, analyzing different genetic markers in could provide an insight into highly heterozygous species.

Phylogeny and multiple independent whole-genome duplication events in the Brassicales

PREMISE

Whole-genome duplications (WGDs) are prevalent throughout the evolutionary history of plants. For example, dozens of WGDs have been phylogenetically localized across the order Brassicales, specifically, within the family Brassicaceae. A WGD event has also been identified in the Cleomaceae, the sister family to Brassicaceae, yet its placement, as well as that of WGDs in other families in the order, remains unclear.

METHODS

Phylo-transcriptomic data were generated and used to infer a nuclear phylogeny for 74 Brassicales taxa. Genome survey sequencing was also performed on 66 of those taxa to infer a chloroplast phylogeny. These phylogenies were used to assess and confirm relationships among the major families of the Brassicales and within Brassicaceae. Multiple WGD inference methods were then used to assess the placement of WGDs on the nuclear phylogeny.

RESULTS

Well-supported chloroplast and nuclear phylogenies for the Brassicales and the putative placement of the Cleomaceae-specific WGD event Th-ɑ are presented. This work also provides evidence for previously hypothesized WGDs, including a well-supported event shared by at least two members of the Resedaceae family, and a possible event within the Capparaceae.

CONCLUSIONS

Phylogenetics and the placement of WGDs within highly polyploid lineages continues to be a major challenge. This study adds to the conversation on WGD inference difficulties by demonstrating that sampling is especially important for WGD identification and phylogenetic placement. Given its economic importance and genomic resources, the Brassicales continues to be an ideal group for assessing WGD inference methods.

A transcriptome-based study on the phylogeny and evolution of the taxonomically controversial subfamily Apioideae (Apiaceae)

BACKGROUND AND AIMS

A long-standing controversy in the subfamily Apioideae concerns relationships among the major lineages, which has prevented a comprehensive study of their fruits and evolutionary history. Here we use single copy genes (SCGs) generated from transcriptome datasets to generate a reliable species tree and explore the evolutionary history of Apioideae.

METHODS

In total, 3351 SCGs were generated from 27 transcriptome datasets and one genome, and further used for phylogenetic analysis using coalescent-based methods. Fruit morphology and anatomy were studied in combination with the species tree. Eleven SCGs were screened out for dating analysis with two fossils selected for calibration.

KEY RESULTS

A well-supported species tree was generated with a topology [Chamaesieae, (Bupleureae, (Pleurospermeae, (Physospermopsis Clade, (Group C, (Group A, Group B)))))] that differed from previous trees. Daucinae and Torilidinae were not in the tribe Scandiceae and existed as sister groups to the Acronema Clade. Five branches (I-V) of the species tree showed low quartet support but strong local posterior probabilities. Dating analysis suggested that Apioideae originated around 56.64 Mya (95 % highest posterior density interval, 45.18-73.53 Mya).

CONCLUSIONS

This study resolves a controversial phylogenetic relationship in Apioideae based on 3351 SCGs and coalescent-based species tree estimation methods. Gene trees that contributed to the species tree may undergoing rapid evolutionary divergence and incomplete lineage sorting. Fruits of Apioideae might have evolved in two directions, anemochorous and hydrochorous, with epizoochorous as a derived mode. Molecular and morphological evidence suggests that Daucinae and Torilidinae should be restored to the tribe level. Our results provide new insights into the morphological evolution of this subfamily, which may contribute to a better understanding of species diversification in Apioideae. Molecular dating analysis suggests that uplift of the Qinghai-Tibetan Plateau (QTP) and climate changes probably drove rapid speciation and diversification of Apioideae in the QTP region.

Target Nuclear and Off-Target Plastid Hybrid Enrichment Data Inform a Range of Evolutionary Depths in the Orchid Genus Epidendrum

Universal angiosperm enrichment probe sets designed to enrich hundreds of putatively orthologous nuclear single-copy loci are increasingly being applied to infer phylogenetic relationships of different lineages of angiosperms at a range of evolutionary depths. Studies applying such probe sets have focused on testing the universality and performance of the target nuclear loci, but they have not taken advantage of off-target data from other genome compartments generated alongside the nuclear loci. Here we do so to infer phylogenetic relationships in the orchid genus Epidendrum and closely related genera of subtribe Laeliinae. Our aims are to: 1) test the technical viability of applying the plant anchored hybrid enrichment (AHE) method (Angiosperm v.1 probe kit) to our focal group, 2) mine plastid protein coding genes from off-target reads; and 3) evaluate the performance of the target nuclear and off-target plastid loci in resolving and supporting phylogenetic relationships along a range of taxonomical depths. Phylogenetic relationships were inferred from the nuclear data set through coalescent summary and site-based methods, whereas plastid loci were analyzed in a concatenated partitioned matrix under maximum likelihood. The usefulness of target and flanking non-target nuclear regions and plastid loci was assessed through the estimation of their phylogenetic informativeness. Our study successfully applied the plant AHE probe kit to Epidendrum, supporting the universality of this kit in angiosperms. Moreover, it demonstrated the feasibility of mining plastome loci from off-target reads generated with the Angiosperm v.1 probe kit to obtain additional, uniparentally inherited sequence data at no extra sequencing cost. Our analyses detected some strongly supported incongruences between nuclear and plastid data sets at shallow divergences, an indication of potential lineage sorting, hybridization, or introgression events in the group. Lastly, we found that the per site phylogenetic informativeness of the ycf1 plastid gene surpasses that of all other plastid genes and several nuclear loci, making it an excellent candidate for assessing phylogenetic relationships at medium to low taxonomic levels in orchids.