Resolution of the ordinal phylogeny of mosses using targeted exons from organellar and nuclear genomes

SynGenes: a Python class for standardizing nomenclatures of mitochondrial and chloroplast genes and a web form for enhancing searches for evolutionary analyses

BACKGROUND

The reconstruction of the evolutionary history of organisms has been greatly influenced by the advent of molecular techniques, leading to a significant increase in studies utilizing genomic data from different species. However, the lack of standardization in gene nomenclature poses a challenge in database searches and evolutionary analyses, impacting the accuracy of results obtained.

RESULTS

To address this issue, a Python class for standardizing gene nomenclatures, SynGenes, has been developed. It automatically recognizes and converts different nomenclature variations into a standardized form, facilitating comprehensive and accurate searches. Additionally, SynGenes offers a web form for individual searches using different names associated with the same gene. The SynGenes database contains a total of 545 gene name variations for mitochondrial and 2485 for chloroplasts genes, providing a valuable resource for researchers.

CONCLUSIONS

The SynGenes platform offers a solution for standardizing gene nomenclatures of mitochondrial and chloroplast genes and providing a standardized search solution for specific markers in GenBank. Evaluation of SynGenes effectiveness through research conducted on GenBank and PubMedCentral demonstrated its ability to yield a greater number of outcomes compared to conventional searches, ensuring more comprehensive and accurate results. This tool is crucial for accurate database searches, and consequently, evolutionary analyses, addressing the challenges posed by non-standardized gene nomenclature.

Characterization of the fiber-like cortical cells in moss gametophytes

MAIN CONCLUSION

Fiber-like cells with thickened cell walls of specific structure and polymer composition that includes (1 → 4)-β-galactans develop in the outer stem cortex of several moss species gametophytes. The early land plants evolved several specialized cell types and tissues that did not exist in their aquatic ancestors. Of these, water-conducting elements and reproductive organs have received most of the research attention. The evolution of tissues specialized to fulfill a mechanical function is by far less studied despite their wide distribution in land plants. For vascular plants following a homoiohydric trajectory, the evolutionary emergence of mechanical tissues is mainly discussed starting with the fern-like plants with their hypodermal sterome or sclerified fibers that have xylan and lignin-based cell walls. However, mechanical challenges were also faced by bryophytes, which lack lignified cell-walls. To characterize mechanical tissues in the bryophyte lineage, following a poikilohydric trajectory, we used six wild moss species (Polytrichum juniperinum, Dicranum sp., Rhodobryum roseum, Eurhynchiadelphus sp., Climacium dendroides, and Hylocomium splendens) and analyzed the structure and composition of their cell walls. In all of them, the outer stem cortex of the leafy gametophytic generation had fiber-like cells with a thickened but non-lignified cell wall. Such cells have a spindle-like shape with pointed tips. The additional thick cell wall layer in those fiber-like cells is composed of sublayers with structural evidence for different cellulose microfibril orientation, and with specific polymer composition that includes (1 → 4)-β-galactans. Thus, the basic cellular characters of the cells that provide mechanical support in vascular plant taxa (elongated cell shape, location at the periphery of a primary organ, the thickened cell wall and its peculiar composition and structure) also exist in mosses.

Mitochondrial genome structure and composition in 70 fishes: a key resource for fisheries management in the South Atlantic

BACKGROUND

Phylogenetic gaps of public databases of reference sequences are a major obstacle for comparative genomics and management of marine resources, particularly in the Global South, where economically important fisheries and conservation flagship species often lack closely-related references. We applied target-enrichment to obtain complete mitochondrial genomes of marine ichthyofauna from the Brazilian coast selected based on economic significance, conservation status and lack of phylogenetically-close references. These included sardines (Dorosomatidae, Alosidae), mackerels (Scombridae) croakers (Sciaenidae), groupers (Epinephelidae) and snappers (Lutjanidae).

RESULTS

Custom baits were designed to enrich mitochondrial DNA across a broad phylogenetic range of fishes. Sequencing generated approximately 100k reads per sample, which were assembled in a total of 70 complete mitochondrial genomes and include fifty-two new additions to GenBank, including five species with no previous mitochondrial data. Departures from the typical gene content and order occurred in only three taxa and mostly involved tRNA gene duplications. Start-codons for all genes, except Cytochrome C Oxidase subunit I (COI), were consistently ATG, whilst a wide range of stop-codons deviated from the prevailing TAA. Phylogenetic analysis confirmed assembly accuracy and revealed signs of cryptic diversification within the Mullus genus. Lineage delimitation methods using Sardinella aurita and S. brasiliensis mitochondrial genomes support a single Operational Taxonomic Unit.

CONCLUSIONS

Target enrichment was highly efficient, providing complete novel mitochondrial genomes with little sequencing effort. These sequences are deposited in public databases to enable subsequent studies in population genetics and adaptation of Latin American fish species and serve as a vital resource for conservation and management programs that rely on molecular data for species and genus-level identification.

Compositae-ParaLoss-1272: A complementary sunflower-specific probe set reduces paralogs in phylogenomic analyses of complex systems

Premise

A family-specific probe set for sunflowers, Compositae-1061, enables family-wide phylogenomic studies and investigations at lower taxonomic levels, but may lack resolution at genus to species levels, especially in groups complicated by polyploidy and hybridization.

Methods

We developed a Hyb-Seq probe set, Compositae-ParaLoss-1272, that targets orthologous loci in Asteraceae. We tested its efficiency across the family by simulating target enrichment sequencing in silico. Additionally, we tested its effectiveness at lower taxonomic levels in the historically complex genus Packera. We performed Hyb-Seq with Compositae-ParaLoss-1272 for 19 Packera taxa that were previously studied using Compositae-1061. The resulting sequences from each probe set, plus a combination of both, were used to generate phylogenies, compare topologies, and assess node support.

Results

We report that Compositae-ParaLoss-1272 captured loci across all tested Asteraceae members, had less gene tree discordance, and retained longer loci than Compositae-1061. Most notably, Compositae-ParaLoss-1272 recovered substantially fewer paralogous sequences than Compositae-1061, with only ~5% of the recovered loci reporting as paralogous, compared to ~59% with Compositae-1061.

Discussion

Given the complexity of plant evolutionary histories, assigning orthology for phylogenomic analyses will continue to be challenging. However, we anticipate Compositae-ParaLoss-1272 will provide improved resolution and utility for studies of complex groups and lower taxonomic levels in the sunflower family.

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.

A target enrichment probe set for resolving phylogenetic relationships in the coffee family, Rubiaceae

Premise

Rubiaceae is among the most species-rich plant families, as well as one of the most morphologically and geographically diverse. Currently available phylogenies have mostly relied on few genomic and plastid loci, as opposed to large-scale genomic data. Target enrichment provides the ability to generate sequence data for hundreds to thousands of phylogenetically informative, single-copy loci, which often leads to improved phylogenetic resolution at both shallow and deep taxonomic scales; however, a publicly accessible Rubiaceae-specific probe set that allows for comparable phylogenetic inference across clades is lacking.

Methods

Here, we use publicly accessible genomic resources to identify putatively single-copy nuclear loci for target enrichment in two Rubiaceae groups: tribe Hillieae (Cinchonoideae) and tribal complex Palicoureeae+Psychotrieae (Rubioideae). We sequenced 2270 exonic regions corresponding to 1059 loci in our target clades and generated in silico target enrichment sequences for other Rubiaceae taxa using our designed probe set. To test the utility of our probe set for phylogenetic inference across Rubiaceae, we performed a coalescent-aware phylogenetic analysis using a subset of 27 Rubiaceae taxa from 10 different tribes and three subfamilies, and one outgroup in Apocynaceae.

Results

We recovered an average of 75% and 84% of targeted exons and loci, respectively, per Rubiaceae sample. Probes designed using genomic resources from a particular subfamily were most efficient at targeting sequences from taxa in that subfamily. The number of paralogs recovered during assembly varied for each clade. Phylogenetic inference of Rubiaceae with our target regions resolves relationships at various scales. Relationships are largely consistent with previous studies of relationships in the family with high support (≥0.98 local posterior probability) at nearly all nodes and evidence of gene tree discordance.

Discussion

Our probe set, which we call Rubiaceae2270x, was effective for targeting loci in species across and even outside of Rubiaceae. This probe set will facilitate phylogenomic studies in Rubiaceae and advance systematics and macroevolutionary studies in the family.

Chromosome-level genome assembly of Niphotrichum japonicum provides new insights into heat stress responses in mosses

With a diversity of approximately 22,000 species, bryophytes (hornworts, liverworts, and mosses) represent a major and diverse lineage of land plants. Bryophytes can thrive in many extreme environments as they can endure the stresses of drought, heat, and cold. The moss Niphotrichum japonicum (Grimmiaceae, Grimmiales) can subsist for extended periods under heat and drought conditions, providing a good candidate for studying the genetic basis underlying such high resilience. Here, we de novo assembled the genome of N. japonicum using Nanopore long reads combined with Hi-C scaffolding technology to anchor the 191.61 Mb assembly into 14 pseudochromosomes. The genome structure of N. japonicum's autosomes is mostly conserved and highly syntenic, in contrast to the sparse and disordered genes present in its sex chromosome. Comparative genomic analysis revealed the presence of 10,019 genes exclusively in N. japonicum. These genes may contribute to the species-specific resilience, as demonstrated by the gene ontology (GO) enrichment. Transcriptome analysis showed that 37.44% (including 3,107 unique genes) of the total annotated genes (26,898) exhibited differential expression as a result of heat-induced stress, and the mechanisms that respond to heat stress are generally conserved across plants. These include the upregulation of HSPs, LEAs, and reactive oxygen species (ROS) scavenging genes, and the downregulation of PPR genes. N. japonicum also appears to have distinctive thermal mechanisms, including species-specific expansion and upregulation of the Self-incomp_S1 gene family, functional divergence of duplicated genes, structural clusters of upregulated genes, and expression piggybacking of hub genes. Overall, our study highlights both shared and species-specific heat tolerance strategies in N. japonicum, providing valuable insights into the heat tolerance mechanism and the evolution of resilient plants.

Symbiosis between cyanobacteria and plants: from molecular studies to agronomic applications

Nitrogen-fixing cyanobacteria from the order Nostocales are able to establish symbiotic relationships with diverse plant species. They are promiscuous symbionts, as the same strain of cyanobacterium is able to form symbiotic biological nitrogen-fixing relationships with different plants species. This review will focus on the different types of cyanobacterial-plant associations, both endophytic and epiphytic, and provide insights from a structural viewpoint, as well as our current understanding of the mechanisms involved in the symbiotic crosstalk. In all these symbioses, the benefit for the plant is clear; it obtains from the cyanobacterium fixed nitrogen and other bioactive compounds, such as phytohormones, polysaccharides, siderophores, or vitamins, leading to enhanced plant growth and productivity. Additionally, there is increasing use of different cyanobacterial species as bio-inoculants for biological nitrogen fixation to improve soil fertility and crop production, thus providing an eco-friendly, alternative, and sustainable approach to reduce the over-reliance on synthetic chemical fertilizers.

Phylogenomics of weevils revisited: data curation and modelling compositional heterogeneity

Weevils represent one of the most prolific radiations of beetles and the most diverse group of herbivores on land. The phylogeny of weevils (Curculionoidea) has received extensive attention, and a largely satisfactory framework for their interfamilial relationships has been established. However, a recent phylogenomic study of Curculionoidea based on anchored hybrid enrichment (AHE) data yielded an abnormal placement for the family Belidae (strongly supported as sister to Nemonychidae + Anthribidae). Here we reanalyse the genome-scale AHE data for Curculionoidea using various models of molecular evolution and data filtering methods to mitigate anticipated systematic errors and reduce compositional heterogeneity. When analysed with the infinite mixture model CAT-GTR or using appropriately filtered datasets, Belidae are always recovered as sister to the clade (Attelabidae, (Caridae, (Brentidae, Curculionidae))), which is congruent with studies based on morphology and other sources of molecular data. Although the relationships of the 'higher Curculionidae' remain challenging to resolve, we provide a consistent and robust backbone phylogeny of weevils. Our extensive analyses emphasize the significance of data curation and modelling across-site compositional heterogeneity in phylogenomic studies.

Magnetically controlled assembly: a new approach to organic integrated photonics

Hierarchical self-assembly of organic molecules or assemblies is of great importance for organic photonics to move from fundamental research to integrated and practical applications. Magnetic fields with the advantages of high controllability, non-contact manipulation, and instantaneous response have emerged as an elegant way to prepare organic hierarchical nanostructures. In this perspective, we outline the development history of organic photonic materials and highlight the importance of organic hierarchical nanostructures for a wide range of applications, including microlasers, optical displays, information encoding, sensing, and beyond. Then, we will discuss recent advances in magnetically controlled assembly for creating organic hierarchical nanostructures, with a particular focus on their potential for enabling the development of integrated photonic devices with unprecedented functionality and performance. Finally, we present several perspectives on the further development of magnetically controlled assembly strategies from the perspective of performance optimization and functional design of organic integrated photonics.

hybpiper-nf and paragone-nf: Containerization and additional options for target capture assembly and paralog resolution

Premise

The HybPiper pipeline has become one of the most widely used tools for the assembly of target capture data for phylogenomic analysis. After the production of locus sequences and before phylogenetic analysis, the identification of paralogs is a critical step for ensuring the accurate inference of evolutionary relationships. Algorithmic approaches using gene tree topologies for the inference of ortholog groups are computationally efficient and broadly applicable to non-model organisms, especially in the absence of a known species tree.

Methods and Results

We containerized and expanded the functionality of both HybPiper and a pipeline for the inference of ortholog groups, providing novel options for the treatment of target capture sequence data, and allowing seamless use of the outputs of the former as inputs for the latter. The Singularity container presented here includes all dependencies, and the corresponding pipelines (hybpiper-nf and paragone-nf, respectively) are implemented via two Nextflow scripts for easier deployment and to vastly reduce the number of commands required for their use.

Conclusions

The hybpiper-nf and paragone-nf pipelines are easily installed and provide a user-friendly experience and robust results to the phylogenetic community. They are used by the Australian Angiosperm Tree of Life project. The pipelines are available at https://github.com/chrisjackson-pellicle/hybpiper-nf and https://github.com/chrisjackson-pellicle/paragone-nf.

Nano-Strainer: A workflow for the identification of single-copy nuclear loci for plant systematic studies, using target capture kits and Oxford Nanopore long reads

In modern plant systematics, target enrichment enables simultaneous analysis of hundreds of genes. However, when dealing with reticulate or polyploidization histories, few markers may suffice, but often are required to be single-copy, a condition that is not necessarily met with commercial capture kits. Also, large genome sizes can render target capture ineffective, so that amplicon sequencing would be preferable; however, knowledge about suitable loci is often missing. Here, we present a comprehensive workflow for the identification of putative single-copy nuclear markers in a genus of interest, by mining a small dataset from target capture using a few representative taxa. The proposed pipeline assesses sequence variability contained in the data from targeted loci and assigns reads to their respective genes, via a combined BLAST/clustering procedure. Cluster consensus sequences are then examined based on four pre-defined criteria presumably indicative for absence of paralogy. This is done by calculating four specialized indices; loci are ranked according to their performance in these indices, and top-scoring loci are considered putatively single- or low copy. The approach can be applied to any probe set. As it relies on long reads, the present contribution also provides template workflows for processing Nanopore-based target capture data. Obtained markers are further tested and then entered into amplicon sequencing. For the detection of possibly remaining paralogy in these data, which might occur in groups with rampant paralogy, we also employ the long-read assembly tool canu. In diploid representatives of the young Compositae genus Leucanthemum, characterized by high levels of polyploidy, our approach resulted in successful amplification of 13 loci. Modifications to remove traces of paralogy were made in seven of these. A species tree from the markers correctly reproduced main relationships in the genus, however, at low resolution. The presented workflow has the potential to valuably support phylogenetic research, for example in polyploid plant groups.

The complete mitochondrial genome of an Antarctic moss, Andreaea regularis Müll. Hal. 1890 (Andreaeaceae)

In the present study, we determined the complete mitochondrial genome of Andreaea regularis Müll. Hal. 1890, a lantern moss of the genus Andreaea Hedw. (Andreaeaceae). The A. regularis mitochondrial genome, with a total length of 118,833 bp, consists of 40 protein-coding genes, 3 ribosomal RNA genes, and 24 transfer RNA genes. A phylogenetic tree constructed with 19 complete mitochondrial genomes composed of liverworts, hornworts, and 15 mosses showed that Andreaeales formed the closest sister to Sphagnales before divergence of the remaining moss groups, indicating A. regularis being one of the earliest mosses. Our findings could be beneficial to investigate the bryophyte evolution.

Fatty Acid Profiles of Some Siberian Bryophytes and Prospects of Their Use in Chemotaxonomy

The composition of fatty acids (FAs) in gametophyte samples of 20 Siberian bryophyte species from four orders of mosses and four orders of liverworts collected in relatively cold months (April and/or October) was examined. FA profiles were obtained using gas chromatography. Thirty-seven FAs were found, from 12:0 to 26:0; they included mono-, polyunsaturated (PUFAs) and rare FAs, such as 22:5n-3 and two acetylenic FAs, 6a,9,12-18:3 and 6a,9,12,15-18:4 (dicranin). Acetylenic FAs were found in all examined species of the Bryales and Dicranales orders, dicranin being the predominant FA. The role of particular PUFAs in mosses and liverworts is discussed. Multivariate discriminant analysis (MDA) was performed to determine whether FAs can be used in the chemotaxonomy of bryophytes. Based on the MDA results, FA composition is related to the taxonomic status of species. Thus, several individual FAs were identified as chemotaxonomic markers at the level of bryophyte orders. These were 18:3n-3; 18:4n-3; 6a,9,12-18:3; 6a,9,12,15-18:4; 20:4n-3 and EPA in mosses and 16:3n-3; 16:2n-6; 18:2n-6; 18:3n-3 and EPA in liverworts. These findings indicate that further research into bryophyte FA profiles can shed light on phylogenetic relationships within this group of plants and the evolution of their metabolic pathways.

Complete chloroplast and mitochondrial genomes of Ditrichum rhynchostegium Kindb. (Ditrichaceae, Bryophyta)

The moss family Pottiaceae is one of the most diverse lineages of the subclass Dicranidae (haplolepideous mosses). Nevertheless, the phylogenetic relationships of Pottiaceae with other Dicranidae families remain unclear. To better understand the ancestral genomic structure and evolution of the Pottiaceae, herein, we present the chloroplast and mitochondrial genomes of Ditrichum rhynchostegium (Ditrichaceae, Bryophyta). The chloroplast genome is 124,628 bp long and displayed a circular structure composed of a large single-copy region, a small single-copy region, and a pair of inverted repeats. It has 118 genes, including 82 protein-coding genes, 32 tRNA genes, and four rRNA genes. The mitochondrial genome is 106,246 bp long and has a circular structure. It contains 67 genes, including 40 protein-coding genes, 24 tRNA genes, and three rRNA genes. Phylogenetic trees based on the coding sequences strongly support the sister relationship of D. rhynchostegium with all Pottiaceous accessions, and the dextrosely arranged operculum cells suggest its affinity for Pottiaceae. This study also demonstrates that long-read sequencing employing the Nanopore platform facilitates the repair of unassembled or misassembled organellar genomic regions.

Phylotranscriptomics of liverworts: revisiting the backbone phylogeny and ancestral gene duplications

BACKGROUND AND AIMS

With some 7300 extant species, liverworts (Marchantiophyta) represent one of the major land plant lineages. The backbone relationships, such as the phylogenetic position of Ptilidiales, and the occurrence and timing of whole-genome duplications, are still contentious.

METHODS

Based on analyses of the newly generated transcriptome data for 38 liverworts and complemented with those publicly available, we reconstructed the evolutionary history of liverworts and inferred gene duplication events along the 55 taxon liverwort species tree.

KEY RESULTS

Our phylogenomic study provided an ordinal-level liverwort nuclear phylogeny and identified extensive gene tree conflicts and cyto-nuclear incongruences. Gene duplication analyses based on integrated phylogenomics and Ks distributions indicated no evidence of whole-genome duplication events along the backbone phylogeny of liverworts.

CONCLUSIONS

With a broadened sampling of liverwort transcriptomes, we re-evaluated the backbone phylogeny of liverworts, and provided evidence for ancient hybridizations followed by incomplete lineage sorting that shaped the deep evolutionary history of liverworts. The lack of whole-genome duplication during the deep evolution of liverworts indicates that liverworts might represent one of the few major embryophyte lineages whose evolution was not driven by whole-genome duplications.

Multiplexed Target Enrichment Enables Efficient and In-Depth Analysis of Antimicrobial Resistome in Metagenomes

Antibiotic resistance genes (ARGs) pose a serious threat to public health and ecological security in the 21st century. However, the resistome only accounts for a tiny fraction of metagenomic content, which makes it difficult to investigate low-abundance ARGs in various environmental settings. Thus, a highly sensitive, accurate, and comprehensive method is needed to describe ARG profiles in complex metagenomic samples. In this study, we established a high-throughput sequencing method based on targeted amplification, which could simultaneously detect ARGs (n = 251), mobile genetic element genes (n = 8), and metal resistance genes (n = 19) in metagenomes. The performance of amplicon sequencing was compared with traditional metagenomic shotgun sequencing (MetaSeq). A total of 1421 primer pairs were designed, achieving extremely high coverage of target genes. The amplicon sequencing significantly improved the recovery of target ARGs (~9 × 104-fold), with higher sensitivity and diversity, less cost, and computation burden. Furthermore, targeted enrichment allows deep scanning of single nucleotide polymorphisms (SNPs), and elevated SNPs detection was shown in this study. We further performed this approach for 48 environmental samples (37 feces, 20 soils, and 7 sewage) and 16 clinical samples. All samples tested in this study showed high diversity and recovery of targeted genes. Our results demonstrated that the approach could be applied to various metagenomic samples and served as an efficient tool in the surveillance and evolution assessment of ARGs. Access to the resistome using the enrichment method validated in this study enabled the capture of low-abundance resistomes while being less costly and time-consuming, which can greatly advance our understanding of local and global resistome dynamics. IMPORTANCE ARGs, an increasing global threat to human health, can be transferred into health-related microorganisms in the environment by horizontal gene transfer, posing a serious threat to public health. Advancing profiling methods are needed for monitoring and predicting the potential risks of ARGs in metagenomes. Our study described a customized amplicon sequencing assay that could enable a high-throughput, targeted, in-depth analysis of ARGs and detect a low-abundance portion of resistomes. This method could serve as an efficient tool to assess the variation and evolution of specific ARGs in the clinical and natural environment.

Codonoblepharonteae, a New Major Lineage among Orthotrichoideae (Orthotrichaceae, Bryophyta)

Orthotrichoideae aggregates epiphytic mosses widespread throughout temperate regions and high tropical mountains of the world. Recently, important advances have been made in elucidating its phylogenetic relationships and evolutionary patterns. Fourteen genera are currently recognized within the subfamily, which are spread over two main tribes: Orthotricheae, comprising Orthotrichinae and Lewinskyinae, and Zygodonteae. Despite the progress, some groups have received little attention, as is the case of genus Codonoblepharon. Recent studies have suggested that this genus may represent a separate lineage from Zygodonteae, in which it traditionally has been considered. Although, none of the studies were conclusive as they did not include a representative sampling of the Codonoblepharon species. This work aims to evaluate the taxonomic position of Codonoblepharon and its phylogenetic relationships within Orthotrichoideae. For this purpose, we present an updated phylogenetic tree based on four different loci, one belonging to the nuclear genome (ITS2) and the rest to the plastid genome (rps4, trnG and trnL-F). The phylogenetic reconstruction recovers all samples of Codonoblepharon in a monophyletic group, sister to the rest of the subfamily, constituting a lineage independent of the two currently recognized tribes. For this reason, we propose the new tribe Codonoblepharonteae to accommodate Codonoblepharon.

Chemotaxonomic and evolutionary perspectives of Bryophyta based on multivariate analysis of fatty acid fingerprints of Eastern Himalayan mosses

Bryophyta comprises one of the earliest lineages of land plants that had implemented remarkable innovations to their lipid metabolic systems for successful adaptation to terrestrial habitat. This study presents a comprehensive investigation of fatty acid profiles of mosses from Eastern Himalayas with an aim to trace their chemotaxonomic and evolutionary implications. Fatty acid compositions of 40 random mosses belonging to major families of Bryophyta were explored by gas chromatographic analysis. A diverse array of saturated, monounsaturated and polyunsaturated fatty acids including rare acetylenic fatty acids were detected. Hexadecanoic acid (C16:0), 9,12 (Z,Z)-octadecadienoic acid (C18:2n6) and 9,12,15 (Z,Z,Z)-octadecatrienoic acid (C18:3n3) were the predominant fatty acids in all the mosses. However, quantitative variation of C20 polyunsaturated fatty acids (PUFAs), specifically 5,8,11,14 (Z,Z,Z,Z)-eicosatetraenoic acid (C20:4n6), among the investigated mosses was the most prominent outcome. The diplolepidous members of Bryidae, especially the mosses of Hypnales, Bryales and Bartramiales contained higher amount of C20 PUFAs compared with the haplolepidous orders. Principal component analyses based on individual fatty acids and other related parameters validated C20:4n6 content and the ratio of C20:4n6/C18:2n6 as the apparent chemotaxonomic discriminants. The prevalent notion of considering 9,12,15-octadecatrien-6-ynoic acid (C18:4a) as the chemomarker of Dicranaceae has also been challenged, since the compound was detected not only in different families of Dicranales, but also in a Pottiales member, Leptodontium viticulosoides. Therefore, an ensemble of fatty acids instead of a single one can be considered as the chemical signature for taxonomic interpretation which may also be vital from an evolutionary standpoint.

An overview of bioinformatics, genomics, and transcriptomics resources for bryophytes

Bryophytes are useful models for the study of plant evolution, development, plant-fungal symbiosis, stress responses, and gametogenesis. Additionally, their dominant haploid gametophytic phase makes them great models for functional genomics research, allowing straightforward genome editing and gene knockout via CRISPR or homologous recombination. Until 2016, however, the only bryophyte genome sequence published was that of Physcomitrium patens. Throughout recent years, several other bryophyte genomes and transcriptome datasets became available, enabling better comparative genomics in evolutionary studies. The increase in the number of bryophyte genome and transcriptome resources available has yielded a plethora of annotations, databases, and bioinformatics tools to access the new data, which covers the large diversity of this clade and whose biology comprises features such as association with arbuscular mycorrhiza fungi, sex chromosomes, low gene redundancy, or loss of RNA editing genes for organellar transcripts. Here we provide a guide to resources available for bryophytes with regards to genome and transcriptome databases and bioinformatics tools.

A scoping review of bryophyte microbiota: diverse microbial communities in small plant packages

Plant health depends not only on the condition of the plant itself but also on its diverse community of microbes, or microbiota. Just like the better-studied angiosperms, bryophytes (mosses, liverworts, and hornworts) harbor diverse communities of bacteria, archaea, fungi, and other microbial eukaryotes. Bryophytes are increasingly recognized as important model systems for understanding plant evolution, development, physiology, and symbiotic interactions. Much of the work on bryophyte microbiota in the past focused on specific symbiont types for each bryophyte group, but more recent studies are taking a broader view acknowledging the coexistence of diverse microbial communities in bryophytes. Therefore, this review integrates studies of bryophyte microbes from both perspectives to provide a holistic view of the existing research for each bryophyte group and on key themes. The systematic search also reveals the taxonomic and geographic biases in this field, including a severe under-representation of the tropics, very few studies on viruses or eukaryotic microbes beyond fungi, and a focus on mycorrhizal fungi studies in liverworts. Such gaps may have led to errors in conclusions about evolutionary patterns in symbiosis. This analysis points to a wealth of future research directions that promise to reveal how the distinct life cycles and physiology of bryophytes interact with their microbiota.

Diversity, phylogeny, and adaptation of bryophytes: insights from genomic and transcriptomic data

Bryophytes including mosses, liverworts, and hornworts are among the earliest land plants, and occupy a crucial phylogenetic position to aid in the understanding of plant terrestrialization. Despite their small size and simple structure, bryophytes are the second largest group of extant land plants. They live ubiquitously in various habitats and are highly diversified, with adaptive strategies to modern ecosystems on Earth. More and more genomes and transcriptomes have been assembled to address fundamental questions in plant biology. Here, we review recent advances in bryophytes associated with diversity, phylogeny, and ecological adaptation. Phylogenomic studies have provided increasing supports for the monophyly of bryophytes, with hornworts sister to the Setaphyta clade including liverworts and mosses. Further comparative genomic analyses revealed that multiple whole-genome duplications might have contributed to the species richness and morphological diversity in mosses. We highlight that the biological changes through gene gain or neofunctionalization that primarily evolved in bryophytes have facilitated the adaptation to early land environments; among the strategies to adapt to modern ecosystems in bryophytes, desiccation tolerance is the most remarkable. More genomic information for bryophytes would shed light on key mechanisms for the ecological success of these 'dwarfs' in the plant kingdom.

Unraveling host-microbe interactions and ecosystem functions in moss-bacteria symbioses

Mosses are non-vascular plants usually found in moist and shaded areas, with great ecological importance in several ecosystems. This is especially true in northern latitudes, where mosses are responsible for up to 100% of primary production in some ecosystems. Mosses establish symbiotic associations with unique bacteria that play key roles in the carbon and nitrogen cycles. For instance, in boreal environments, more than 35% of the nitrogen fixed by diazotrophic symbionts in peatlands is transferred to mosses, directly affecting carbon fixation by the hosts, while moss-associated methanotrophic bacteria contribute 10-30% of moss carbon. Further, half of ecosystem N input may derive from moss-cyanobacteria associations in pristine ecosystems. Moss-bacteria interactions have consequences on a global scale since northern environments sequester 20% of all the carbon generated by forests in the world and stock at least 32% of global terrestrial carbon. Different moss hosts influence bacteria in distinct ways, which suggests that threats to mosses also threaten unique microbial communities with important ecological and biogeochemical consequences. Since their origin ~500 Ma, mosses have interacted with bacteria, making these associations ideal models for understanding the evolution of plant-microbe associations and their contribution to biogeochemical cycles.

Phylogenomic Analysis Reconstructed the Order Matoniales from Paleopolyploidy Veil

Phylogenetic conflicts limit our understanding of the evolution of terrestrial life under multiple whole genome duplication events, and the phylogeny of early terrestrial plants remains full of controversy. Although much incongruence has been solved with so-called robust topology based on single or lower copy genes, the evolutionary mechanisms behind phylogenetic conflicts such as polyploidization remain poorly understood. Here, through decreasing the effects of polyploidization and increasing the samples of species, which represent all four orders and eight families that comprise early leptosporangiate ferns, we have reconstructed a robust phylogenetic tree and network with 1125 1-to-1 orthologs based on both coalescent and concatenation methods. Our data consistently suggest that Matoniales, as a monophyletic lineage including Matoniaceae and Dipteridaceae, should be redefined as an ordinal rank. Furthermore, we have identified and located at least 11 whole-genome duplication events within the evolutionary history of four leptosporangiates lineages, and associated polyploidization with higher speciation rates and mass extinction events. We hypothesize that paleopolyploidization may have enabled leptosporangiate ferns to survive during mass extinction events at the end Permian period and then flourish throughout the Mesozoic era, which is supported by extensive fossil records. Our results highlight how ancient polyploidy can result in rapid species radiation, thus causing phylogenetic conflicts yet allowing plants to survive and thrive during mass extinction events.

An NGS-Based Phylogeny of Orthotricheae (Orthotrichaceae, Bryophyta) With the Proposal of the New Genus Rehubryum From Zealandia

Phylogenomic data increase the possibilities of resolving the evolutionary and systematic relationships among taxa. This is especially valuable in groups with few and homoplasious morphological characters, in which systematic and taxonomical delimitations have been traditionally difficult. Such is the case of several lineages within Bryophyta, like Orthotrichaceae, the second most diverse family of mosses. Members of tribe Orthotricheae are common in temperate and cold regions, as well as in high tropical mountains. In extratropical areas, they represent one of the main components of epiphytic communities, both in dry and oceanic or hyperoceanic conditions. The epiphytic environment is considered a hostile one for plant development, mainly due to its low capacity of moisture retention. Thus, the diversification of the Orthotrichaceae in this environment could be seen as striking. Over the last two decades, great taxonomic and systematic progresses have led to a rearrangement at the generic level in this tribe, providing a new framework to link environment to patterns of diversification. Here, we use nuclear loci targeted with the GoFlag 408 enrichment probe set to generate a well-sampled phylogeny with well-supported suprageneric taxa and increasing the phylogenetic resolution within the two recognized subtribes. Specifically, we show that several genera with Ulota-like morphology jointly constitute an independent lineage. Within this lineage, the recently described Atlantichella from Macaronesia and Western Europe appears as the sister group of Ulota bellii from Zealandia. This latter species is here segregated in the new genus Rehubryum. Assessment of the ecological and biogeographical affinities of the species within the phylogenetic framework suggests that niche adaptation (including climate and substrate) may be a key evolutionary driver that shaped the high diversification of Orthotricheae.

Incomplete lineage sorting and local extinction shaped the complex evolutionary history of the Paleogene relict conifer genus, Chamaecyparis (Cupressaceae)

Inferring accurate biogeographic history of plant taxa with an East Asia (EA)-North America (NA) is usually hindered by conflicting phylogenies and a poor fossil record. The current distribution of Chamaecyparis (false cypress; Cupressaceae) with four species in EA, and one each in western and eastern NA, and its relatively rich fossil record, make it an excellent model for studying the EA-NA disjunction. Here we reconstruct phylogenomic relationships within Chamaecyparis using > 1400 homologous nuclear and 61 plastid genes. Our phylogenomic analyses using concatenated and coalescent approaches revealed strong cytonuclear discordance and conflicting topologies between nuclear gene trees. Incomplete lineage sorting (ILS) and hybridization are possible explanations of conflict; however, our coalescent analyses and simulations suggest that ILS is the major contributor to the observed phylogenetic discrepancies. Based on a well-resolved species tree and four fossil calibrations, the crown lineage of Chamaecyparis is estimated to have originated in the upper Cretaceous, followed by diversification events in the early and middle Paleogene. Ancestral area reconstructions suggest that Chamaecyparis had an ancestral range spanning both EA and NA. Fossil records further indicate that this genus is a relict of the "boreotropical" flora, and that local extinctions of European species were caused by global cooling. Overall, our results unravel a complex evolutionary history of a Paleogene relict conifer genus, which may have involved ILS, hybridization and the extinction of local species.

Breaking the limits - multichromosomal structure of an early eudicot Pulsatilla patens mitogenome reveals extensive RNA-editing, longest repeats and chloroplast derived regions among sequenced land plant mitogenomes

BACKGROUND

The mitogenomes of vascular plants are one of the most structurally diverse molecules. In the present study we characterize mitogenomes of a rare and endangered species Pulsatilla patens. We investigated the gene content and its RNA editing potential, repeats distribution and plastid derived sequences.

RESULTS

The mitogenome structure of early divergent eudicot, endangered Pulsatilla patens does not support the master chromosome hypothesis, revealing the presence of three linear chromosomes of total length 986 613 bp. The molecules are shaped by the presence of extremely long, exceeding 87 kbp repeats and multiple chloroplast-derived regions including nearly complete inverted repeat. Since the plastid IR content of Ranunculales is very characteristic, the incorporation into mitogenome could be explained rather by intracellular transfer than mitochondrial HGT. The mitogenome contains almost a complete set of genes known from other vascular plants with exception of rps10 and sdh3, the latter being present but pseudogenized. Analysis of long ORFs enabled the identification of genes which are rarely present in plant mitogenomes, including RNA and DNA polymerases, albeit their presence even at species level is variable. Mitochondrial transcripts of P. patens were edited with a high frequency, which exceeded the level known in other analyzed angiosperms, despite the strict qualification criteria of counting the editing events and taking into analysis generally less frequently edited leaf transcriptome. The total number of edited sites was 902 and nad4 was identified as the most heavily edited gene with 65 C to U changes. Non-canonical, reverse U to C editing was not detected. Comparative analysis of mitochondrial genes of three Pulsatilla species revealed a level of variation comparable to chloroplast CDS dataset and much higher infrageneric differentiation than in other known angiosperm genera. The variation found in CDS of mitochondrial genes is comparable to values found among Pulsatilla plastomes. Despite the complicated mitogenome structure, 14 single copy regions of 329 kbp, not splitted by repeats or plastid-derived sequences (MTPT), revealed the potential for phylogenetic, phylogeographic and population genetics studies by revealing intra- and interspecific collinearity.

CONCLUSIONS

This study provides valuable new information about mitochondrial genome of early divergent eudicots, Pulsatilla patens, revealed multi-chromosomal structure and shed new light on mitogenomics of early eudicots.

Taming the beast: a revised classification of Cortinariaceae based on genomic data

Family Cortinariaceae currently includes only one genus, Cortinarius, which is the largest Agaricales genus, with thousands of species worldwide. The species are important ectomycorrhizal fungi and form associations with many vascular plant genera from tropicals to arctic regions. Genus Cortinarius contains a lot of morphological variation, and its complexity has led many taxonomists to specialize in particular on infrageneric groups. The previous attempts to divide Cortinarius have been shown to be unnatural and the phylogenetic studies done to date have not been able to resolve the higher-level classification of the group above section level. Genomic approaches have revolutionized our view on fungal relationships and provide a way to tackle difficult groups. We used both targeted capture sequencing and shallow whole genome sequencing to produce data and to perform phylogenomic analyses of 75 single-copy genes from 19 species. In addition, a wider 5-locus analysis of 245 species, from the Northern and Southern Hemispheres, was also done. Based on our results, a classification of the family Cortinariaceae into ten genera—Cortinarius, Phlegmacium, Thaxterogaster, Calonarius, Aureonarius, Cystinarius, Volvanarius, Hygronarius, Mystinarius, and Austrocortinarius—is proposed. Seven genera, 10 subgenera, and four sections are described as new to science and five subgenera are introduced as new combinations in a new rank. In addition, 41 section names and 514 species names are combined in new genera and four lecto- and epitypes designated. The position of Stephanopus in suborder Agaricineae remains to be studied. Targeted capture sequencing is used for the first time in fungal taxonomy in Basidiomycetes. It provides a cost-efficient way to produce -omics data in species-rich groups. The -omics data was produced from fungarium specimens up to 21 years old, demonstrating the value of museum specimens in the study of the fungal tree of life. This study is the first family revision in Agaricales based on genomics data and hopefully many others will soon follow.

Stress-associated developmental reprogramming in moss protonemata by synthetic activation of the common symbiosis pathway

Symbioses between angiosperms and rhizobia or arbuscular mycorrhizal fungi are controlled through a conserved signaling pathway. Microbe-derived, chitin-based elicitors activate plant cell surface receptors and trigger nuclear calcium oscillations, which are decoded by a calcium/calmodulin-dependent protein kinase (CCaMK) and its target transcription factor interacting protein of DMI3 (IPD3). Genes encoding CCaMK and IPD3 have been lost in multiple non-mycorrhizal plant lineages yet retained among non-mycorrhizal mosses. Here, we demonstrated that the moss Physcomitrium is equipped with a bona fide CCaMK that can functionally complement a Medicago loss-of-function mutant. Conservation of regulatory phosphosites allowed us to generate predicted hyperactive forms of Physcomitrium CCaMK and IPD3. Overexpression of synthetically activated CCaMK or IPD3 in Physcomitrium led to abscisic acid (ABA) accumulation and ectopic development of brood cells, which are asexual propagules that facilitate escape from local abiotic stresses. We therefore propose a functional role for Physcomitrium CCaMK-IPD3 in stress-associated developmental reprogramming.

The complete chloroplast and mitochondrial genomes of Scopelophila cataractae (Mitt.) Broth. (Pottiaceae, Bryophyta)

The complete chloroplast and mitochondrial genome sequences of Scopelophila cataractae (Pottiaceae, Bryophyta) are determined. The chloroplast genome is 122,290 bp with 118 genes and the mitochondrial genome is 105,607 bp with 67 genes, both genomes are circular. This study showed the S. cataractae plastome contains the smallest genome size, and a functional trnP^GGG gene, relative to other pottiaceous species. Phylogenetic inferences support the sister relationship of S. cataractae to all other pottiaceous accessions. 

Analysis of Paralogs in Target Enrichment Data Pinpoints Multiple Ancient Polyploidy Events in Alchemilla s.l. (Rosaceae)

Target enrichment is becoming increasingly popular for phylogenomic studies. Although baits for enrichment are typically designed to target single-copy genes, paralogs are often recovered with increased sequencing depth, sometimes from a significant proportion of loci, especially in groups experiencing whole-genome duplication (WGD) events. Common approaches for processing paralogs in target enrichment data sets include random selection, manual pruning, and mainly, the removal of entire genes that show any evidence of paralogy. These approaches are prone to errors in orthology inference or removing large numbers of genes. By removing entire genes, valuable information that could be used to detect and place WGD events is discarded. Here, we used an automated approach for orthology inference in a target enrichment data set of 68 species of Alchemilla s.l. (Rosaceae), a widely distributed clade of plants primarily from temperate climate regions. Previous molecular phylogenetic studies and chromosome numbers both suggested ancient WGDs in the group. However, both the phylogenetic location and putative parental lineages of these WGD events remain unknown. By taking paralogs into consideration and inferring orthologs from target enrichment data, we identified four nodes in the backbone of Alchemilla s.l. with an elevated proportion of gene duplication. Furthermore, using a gene-tree reconciliation approach, we established the autopolyploid origin of the entire Alchemilla s.l. and the nested allopolyploid origin of four major clades within the group. Here, we showed the utility of automated tree-based orthology inference methods, previously designed for genomic or transcriptomic data sets, to study complex scenarios of polyploidy and reticulate evolution from target enrichment data sets.[Alchemilla; allopolyploidy; autopolyploidy; gene tree discordance; orthology inference; paralogs; Rosaceae; target enrichment; whole genome duplication.].

The complete plastid genome sequence of the enigmatic moss, Takakia lepidozioides (Takakiopsida, Bryophyta): evolutionary perspectives on the largest collection of genes in mosses and the intensive RNA editing

Complete chloroplast genome sequence of a moss, Takakia lepidozioides (Takakiopsida) is reported. The largest collection of genes in mosses and the intensive RNA editing were discussed from evolutionary perspectives. We assembled the entire plastid genome sequence of Takakia lepidozioides (Takakiopsida), emerging from the first phylogenetic split among extant mosses. The genome sequences were assembled into a circular molecule 149,016 bp in length, with a quadripartite structure comprising a large and a small single-copy region separated by inverted repeats. It contained 88 genes coding for proteins, 32 for tRNA, four for rRNA, two open reading frames, and at least one pseudogene (tufA). This is the largest number of genes of all sequenced plastid genomes in mosses and Takakia is the only moss that retains the seven coding genes ccsA, cysA, cysT, petN rpoA, rps16 and trnP^GGG. Parsimonious interpretation of gene loss suggests that the last common ancestor of bryophytes had all seven genes and that mosses lost at least three of them during their diversification. Analyses of the plastid transcriptome identified the extraordinary frequency of RNA editing with more than 1100 sites. We indicated a close correlation between the monoplastidy of vegetative tissue and the intensive RNA editing sites in the plastid genome in land plant lineages. Here, we proposed a hypothesis that the small population size of plastids in each vegetative cell of some early diverging land plants, including Takakia, might cause the frequent fixation of mutations in plastid genome through the intracellular genetic drift and that deleterious mutations might be continuously compensated by RNA editing during or following transcription.

Insight into early diversification of leucine-rich repeat receptor-like kinases provided by the sequenced moss and hornwort genomes

Identification of the subfamily X leucine-rich repeat receptor-like kinases in the recently sequenced moss and hornwort genomes points to their diversification into distinct groups during early evolution of land plants. Signal transduction mediated through receptor-ligand interactions plays key roles in controlling developmental and physiological processes of multicellular organisms, and plants employ diverse receptors in signaling. Leucine-rich repeat receptor-like kinases (LRR-RLKs) represent one of the largest receptor classes in plants and are structurally classified into subfamilies. LRR-RLKs of the subfamily X are unique in the variety of their signaling roles; they include receptors for steroid or peptide hormones as well as negative regulators of signaling through binding to other LRR-RLKs, raising a question as to how they diversified. However, our understanding of diversification processes of LRR-RLKs has been hindered by the paucity of genomic data in non-seed plants and limited taxa sampling in previous phylogenetic analyses. Here we analyzed the phylogeny of LRR-RLK X sequences collected from all major land plant lineages and show that this subfamily diversified into six major clades before the divergence between bryophytes and vascular plants. Notably, we have identified homologues of the brassinosteroid receptor, BRASSINOSTEROID INSENSITIVE 1 (BRI1), in the genomes of Sphagnum mosses, hornworts, and ferns, contrary to earlier reports that postulate the origin of BRI1-like LRR-RLKs in the seed plant lineage. The phylogenetic distribution of major clades illustrates that the current receptor repertoire was shaped through lineage-specific gene family expansion and independent gene losses, highlighting dynamic changes in the evolution of LRR-RLKs.

Complete genus-level plastid phylogenomics of Alismataceae with revisited historical biogeography

Alismataceae, an ancient lineage of monocots, has attracted attention due to its complex evolutionary history, ornamental value, and ecological role. However, the phylogenetic relationships and evolutionary history of the family have not been conclusively resolved. Here, we constructed the first complete genus-level plastid phylogeny of Alismataceae by using 78 genes and updated the historical biogeography based on the phylogenomic tree. Our results divide the Alismataceae into three major clades with robust support values; one clade comprises the former Limnocharitaceae, and the second clade includes the mainly temperate genera Alisma, Baldellia, Damasonium and Luronium, and the monotypic African genus Burnatia as a sister of the temperate genera. The remaining genera are either tropical or have some temperate species in addition to tropical ones, and they constitute the third major clade. Molecular dating and biogeographic analyses suggest that Alismataceae arose in Neotropical, West Palearctic, and Afrotropical regions during the Cretaceous, followed by the split into three main clades due to a combination of vicariance and dispersal events. Unlike earlier studies, we inferred that the mainly temperate clade likely originated from Afrotropical and West Palearctic regions during the Eocene. The most recent common ancestor of the other two clades lived in the Neotropical area during the Late Cretaceous. Long-distance dispersal and vicariance together seem to contribute to the transoceanic distribution of this family.

Effects of extraction parameters on lipid profiling of mosses using UPLC-ESI-QTOF-MS and multivariate data analysis

INTRODUCTION

Non-target lipid profiling by using ultra-performance liquid chromatography coupled to electrospray ionization quadrupole time-of-flight mass spectrometry (UPLC-ESI-QTOF-MS) has been used extensively in the past decades in plant studies. However, the lipidomes of bryophytes have only been scarcely studied, although they are the second largest group in plant kingdom.

OBJECTIVES

We evaluated the effects of different cell disruption methods (no disruption, shake, ultrasound, and bead beating), and storage conditions (air-dried, freeze-dried, and fresh frozen) of five moss species (including Racomitrium lanuginosum B and D, Philonotis fontana, Sphagnum teres, and Hylocomium splendens).

METHODS

The lipid profiling results of each extraction parameter were analyzed by using multivariate data analysis including unsupervised principal component analysis and supervised orthogonal projections to latent structures discriminant analysis.

RESULTS

The results showed that extraction with bead beating resulted in the highest lipid content and the most detected features, but these were caused by the contamination from plastic tubes. Minor lipid metabolite changes were found in shaking and ultrasonication methods when compared with no disruption method. Significant amounts of phosphatidylcholine, diacylglyceryltrimethylhomoserine and their lyso lipids were observed in air-dried moss tissues, whereas diacylglycerol, triacylglycerol and ceramide were mostly exclusively detected when fresh frozen tissues were used for extraction.

CONCLUSION

We concluded that lipid extraction using fresh frozen samples with ultrasound assistance provide the most original lipid composition and gave a relatively high lipid content.

High-Throughput Genomic Data Reveal Complex Phylogenetic Relationships in Stylosanthes Sw (Leguminosae)

Allopolyploidy is widely present across plant lineages. Though estimating the correct phylogenetic relationships and origin of allopolyploids may sometimes become a hard task. In the genus Stylosanthes Sw. (Leguminosae), an important legume crop, allopolyploidy is a key speciation force. This makes difficult adequate species recognition and breeding efforts on the genus. Based on comparative analysis of nine high-throughput sequencing (HTS) samples, including three allopolyploids (S. capitata Vogel cv. “Campo Grande,” S. capitata “RS024” and S. scabra Vogel) and six diploids (S. hamata Taub, S. viscosa (L.) Sw., S. macrocephala M. B. Ferreira and Sousa Costa, S. guianensis (Aubl.) Sw., S. pilosa M. B. Ferreira and Sousa Costa and S. seabrana B. L. Maass & 't Mannetje) we provide a working pipeline to identify organelle and nuclear genome signatures that allowed us to trace the origin and parental genome recognition of allopolyploids. First, organelle genomes were de novo assembled and used to identify maternal genome donors by alignment-based phylogenies and synteny analysis. Second, nuclear-derived reads were subjected to repetitive DNA identification with RepeatExplorer2. Identified repeats were compared based on abundance and presence on diploids in relation to allopolyploids by comparative repeat analysis. Third, reads were extracted and grouped based on the following groups: chloroplast, mitochondrial, satellite DNA, ribosomal DNA, repeat clustered- and total genomic reads. These sets of reads were then subjected to alignment and assembly free phylogenetic analyses and were compared to classical alignment-based phylogenetic methods. Comparative analysis of shared and unique satellite repeats also allowed the tracing of allopolyploid origin in Stylosanthes, especially those with high abundance such as the StyloSat1 in the Scabra complex. This satellite was in situ mapped in the proximal region of the chromosomes and made it possible to identify its previously proposed parents. Hence, with simple genome skimming data we were able to provide evidence for the recognition of parental genomes and understand genome evolution of two Stylosanthes allopolyploids.

Taxon-specific or universal? Using target capture to study the evolutionary history of rapid radiations

Target capture has emerged as an important tool for phylogenetics and population genetics in nonmodel taxa. Whereas developing taxon‐specific capture probes requires sustained efforts, available universal kits may have a lower power to reconstruct relationships at shallow phylogenetic scales and within rapidly radiating clades. We present here a newly developed target capture set for Bromeliaceae, a large and ecologically diverse plant family with highly variable diversification rates. The set targets 1776 coding regions, including genes putatively involved in key innovations, with the aim to empower testing of a wide range of evolutionary hypotheses. We compare the relative power of this taxon‐specific set, Bromeliad1776, to the universal Angiosperms353 kit. The taxon‐specific set results in higher enrichment success across the entire family; however, the overall performance of both kits to reconstruct phylogenetic trees is relatively comparable, highlighting the vast potential of universal kits for resolving evolutionary relationships. For more detailed phylogenetic or population genetic analyses, for example the exploration of gene tree concordance, nucleotide diversity or population structure, the taxon‐specific capture set presents clear benefits. We discuss the potential lessons that this comparative study provides for future phylogenetic and population genetic investigations, in particular for the study of evolutionary radiations.

Excluding loci with substitution saturation improves inferences from phylogenomic data

The historical signal in nucleotide sequences becomes eroded over time by substitutions occurring repeatedly at the same sites. This phenomenon, known as substitution saturation, is recognized as one of the primary obstacles to deep-time phylogenetic inference using genome-scale data sets. We present a new test of substitution saturation and demonstrate its performance in simulated and empirical data. For some of the 36 empirical phylogenomic data sets that we examined, we detect substitution saturation in around 50% of loci. We found that saturation tends to be flagged as problematic in loci with highly discordant phylogenetic signals across sites. Within each data set, the loci with smaller numbers of informative sites are more likely to be flagged as containing problematic levels of saturation. The entropy saturation test proposed here is sensitive to high evolutionary rates relative to the evolutionary timeframe, while also being sensitive to several factors known to mislead phylogenetic inference, including short internal branches relative to external branches, short nucleotide sequences, and tree imbalance. Our study demonstrates that excluding loci with substitution saturation can be an effective means of mitigating the negative impact of multiple substitutions on phylogenetic inferences. [Phylogenetic model performance; phylogenomics; substitution model; substitution saturation; test statistics.]

Rates and patterns of molecular evolution in bryophyte genomes, with focus on complex thalloid liverworts, Marchantiopsida

Plants commonly referred to as "bryophytes" belong to three major lineages of non-vascular plants: the liverworts, the hornworts and the mosses. They are unique among land plants in having a dominant haploid generation and a short-lived diploid sporophytic generation. The dynamics of selection acting on a haploid genome differs from those acting on a diploid genome: new mutations are directly exposed to selection. The general aim of this paper is to investigate the diversification rateof bryophytes - measured as silent site substitution rate representing neutral evolution (mutation rate) and the nonsynonymous to synonymous substitution rate ratio (dN/dS) representing selective evolution - and compare it with earlier studies on vascular plants. Results show that the silent site substitution rate is lower for liverworts as compared to angiosperms, but not as low as for gymnosperms. The selection pressure, measured as dN/dS, isnot remarkably lower for bryophytes as compared to other diploid dominant plants as would be expected by the masking hypothesis, indicating that other factors are more important than ploidy.

Construction of DNA Tools for Hyperexpression in Marchantia Chloroplasts

Chloroplasts are attractive platforms for synthetic biology applications since they are capable of driving very high levels of transgene expression, if mRNA production and stability are properly regulated. However, plastid transformation is a slow process and currently limited to a few plant species. The liverwort Marchantia polymorpha is a simple model plant that allows rapid transformation studies; however, its potential for protein hyperexpression has not been fully exploited. This is partially due to the fact that chloroplast post-transcriptional regulation is poorly characterized in this plant. We have mapped patterns of transcription in Marchantia chloroplasts. Furthermore, we have obtained and compared sequences from 51 bryophyte species and identified putative sites for pentatricopeptide repeat protein binding that are thought to play important roles in mRNA stabilization. Candidate binding sites were tested for their ability to confer high levels of reporter gene expression in Marchantia chloroplasts, and levels of protein production and effects on growth were measured in homoplastic transformed plants. We have produced novel DNA tools for protein hyperexpression in this facile plant system that is a test-bed for chloroplast engineering.

Evolutionary Rate Variation Among Lineages in Gene Trees has a Negative Impact on Species-Tree Inference

Phylogenetic analyses of genomic data provide a powerful means of reconstructing the evolutionary relationships among organisms, yet such analyses are often hindered by conflicting phylogenetic signals among loci. Identifying the signals that are most influential to species-tree estimation can help to inform the choice of data for phylogenomic analysis. We investigated this in an analysis of 30 phylogenomic data sets. For each data set, we examined the association between several branch-length characteristics of gene trees and the distance between these gene trees and the corresponding species trees. We found that the distance of each gene tree to the species tree inferred from the full data set was positively associated with variation in root-to-tip distances and negatively associated with mean branch support. However, no such associations were found for gene-tree length, a measure of the overall substitution rate at each locus. We further explored the usefulness of the best-performing branch-based characteristics for selecting loci for phylogenomic analyses. We found that loci that yield gene trees with high variation in root-to-tip distances have a disproportionately distant signal of tree topology compared with the complete data sets. These results suggest that rate variation across lineages should be taken into consideration when exploring and even selecting loci for phylogenomic analysis.[Branch support; data filtering; nucleotide substitution model; phylogenomics; substitution rate; summary coalescent methods.]

New targets acquired: Improving locus recovery from the Angiosperms353 probe set

PREMISE

Universal target enrichment kits maximize utility across wide evolutionary breadth while minimizing the number of baits required to create a cost-efficient kit. The Angiosperms353 kit has been successfully used to capture loci throughout the angiosperms, but the default target reference file includes sequence information from only 6-18 taxa per locus. Consequently, reads sequenced from on-target DNA molecules may fail to map to references, resulting in fewer on-target reads for assembly, and reducing locus recovery.

METHODS

We expanded the Angiosperms353 target file, incorporating sequences from 566 transcriptomes to produce a 'mega353' target file, with each locus represented by 17-373 taxa. This mega353 file is a drop-in replacement for the original Angiosperms353 file in HybPiper analyses. We provide tools to subsample the file based on user-selected taxon groups, and to incorporate other transcriptome or protein-coding gene data sets.

RESULTS

Compared to the default Angiosperms353 file, the mega353 file increased the percentage of on-target reads by an average of 32%, increased locus recovery at 75% length by 49%, and increased the total length of the concatenated loci by 29%.

DISCUSSION

Increasing the phylogenetic density of the target reference file results in improved recovery of target capture loci. The mega353 file and associated scripts are available at: https://github.com/chrisjackson-pellicle/NewTargets.

Evidence of horizontal gene transfer between land plant plastids has surprising conservation implications

BACKGROUND AND AIMS

Horizontal gene transfer (HGT) is an important evolutionary mechanism because it transfers genetic material that may code for traits or functions between species or genomes. It is frequent in mitochondrial and nuclear genomes but has not been demonstrated between plastid genomes of different green land plant species.

METHODS

We Sanger-sequenced the nuclear internal transcribed spacers (ITS1 and 2) and the plastid rpl16 G2 intron (rpl16). In five individuals with foreign rpl16 we also sequenced atpB-rbcL and trnLUAA-trnFGAA.

KEY RESULTS

We discovered 14 individuals of a moss species with typical nuclear ITSs but foreign plastid rpl16 from a species of a distant lineage. None of the individuals with three plastid markers sequenced contained all foreign markers, demonstrating the transfer of plastid fragments rather than the entire plastid genome, i.e. entire plastids were not transferred. The two lineages diverged 165-185 Myr BP. The extended time interval since lineage divergence suggests that the foreign rpl16 is more likely explained by HGT than by hybridization or incomplete lineage sorting.

CONCLUSIONS

We provide the first conclusive evidence of interspecific plastid-to-plastid HGT among land plants. Two aspects are critical: it occurred at several localities during the massive colonization of recently disturbed open habitats that were created by large-scale liming as a freshwater biodiversity conservation measure; and it involved mosses whose unique life cycle includes spores that first develop a filamentous protonema phase. We hypothesize that gene transfer is facilitated when protonema filaments of different species intermix intimately when colonizing disturbed early succession habitats.

Geographic Pattern of Bryophyte Species Richness in China: The Influence of Environment and Evolutionary History

How contemporary environment interacts with macroevolutionary processes to generate the geographic pattern of bryophyte species is still unresolved. China is very rich in bryophytes, with more than 3,000 bryophytes covering 70% of the families in the world. In this study, we assessed the effects of the contemporary environment (average temperature of the coldest season TCQ, precipitation of the warmest season PWQ, and elevational range) and the recent diversification rates (estimated as mean species number per genus, MSG) on the geographical pattern of species richness for bryophytes and two groups (i.e., liverworts and mosses) in China. We compiled the provincial level distribution of bryophyte species and estimated the geographic pattern of the recent diversification rate by MSG for species in China. Univariate, multivariate regressions and path model analyses were used to assess the relationships between species richness, MSG, and their potential environmental drivers. Species richness of all bryophytes and liverworts significantly increased with the increase of MSG, either in regressions or path analyses, indicating that provinces with high bryophyte richness were mainly inhabited by species (especially liverworts) from lineages with particularly high MSG. In contrast, the species richness of mosses was insignificantly decreased with MSG in univariate regression or insignificantly increased with MSG in path analysis. Both species richness and MSG of all bryophytes and liverworts increased with the increase in energy and water availability. In contrast, for mosses, the species richness significantly increased with the increase of energy and water availability, while MSG decreased with the increase of energy and water availability. The MSG of liverworts increase with the increase of elevational range but the MSG of mosses decrease with the increase of elevational range. Our study suggests that the humid tropical and subtropical mountains in China are not only diversity hotspots for bryophytes, but also cradles for high recent diversification of liverworts, and refuges for mosses to hold many monotypic and oligotypic genera.

The complete plastome of Andreaea rupestris Hedw. (Andreaeaceae, Bryophyta)

Andreaea rupestris Hedw., one of the lantern mosses, is the lectotype of the genus Andreaea Hedw. (Andreaeaceae). Here we present its complete plastome. The plastome of A. rupestris is successfully assembled from raw reads sequenced by HiSeq X ten system. Its total length is 135,214 bp consisting of four regions: large single copy (LSC) region (92,780 bp), small single copy (SSC) region (21,102 bp), and two inverted repeat regions (IRs; 10,666 bp per each). It contains 134 genes (88 coding genes, 8 rRNAs, and 38 tRNAs). The overall GC content is 30.3% and in the LSC, SSC, and IR regions are 27.5%, 26.5%, and 46.2%, respectively. The present data will be an important sequence resource for further studies on the important early diverging lineage of mosses.

The complete plastome of Polytrichum commune Hedw. (Polytrichaceae, Bryophyta)

Polytrichum commune, one of hair-cap mosses, is the type species of the genus Polytrichum Hedw. (Polytrichaceae). Here we present its complete plastome. The plastome of P. commune is successfully assembled from raw reads sequenced by HiSeq X ten system. Its total length is 126,323 bp consisting of four regions: large single copy (LSC) region (88,070 bp), small single copy (SSC) region (16,717 bp), and inverted repeats (IRs; 9,680 bp per each). It contains 128 genes (84 coding genes, eight rRNAs, and 36 tRNAs); nine genes (four rRNAs and five tRNAs) are duplicated in IR regions. The overall GC content is 28.9% and in the LSC, SSC and IR regions is 26.1%, 25.1%, and 45.5%, respectively. This plastome is an important sequence resource for further studies on the class Polytrichopsida.

Large-Scale Phylogenomic Analyses Reveal the Monophyly of Bryophytes and Neoproterozoic Origin of Land Plants

The relationships among the four major embryophyte lineages (mosses, liverworts, hornworts, vascular plants) and the timing of the origin of land plants are enigmatic problems in plant evolution. Here, we resolve the monophyly of bryophytes by improving taxon sampling of hornworts and eliminating the effect of synonymous substitutions. We then estimate the divergence time of crown embryophytes based on three fossil calibration strategies, and reveal that maximum calibration constraints have a major effect on estimating the time of origin of land plants. Moreover, comparison of priors and posteriors provides a guide for evaluating the optimal calibration strategy. By considering the reliability of fossil calibrations and the influences of molecular data, we estimate that land plants originated in the Precambrian (980–682 Ma), much older than widely recognized. Our study highlights the important contribution of molecular data when faced with contentious fossil evidence, and that fossil calibrations used in estimating the timescale of plant evolution require critical scrutiny.

Plastid genomes and phylogenomics of liverworts (Marchantiophyta): Conserved genome structure but highest relative plastid substitution rate in land plants

With some 7300 species of small nonvascular spore-producing plants, liverworts represent one of the major lineages of land plants. Although multi-locus molecular phylogenetic studies have elucidated relationships of liverworts at different taxonomic categories, the backbone phylogeny of liverworts is still to be fully resolved, especially for the placement of Ptilidiales and the relationships within Jungermanniales and Marchantiales. Here, we provided phylogenomic inferences of liverworts based on 42 newly sequenced and 24 published liverwort plastid genomes representing all but two orders of liverworts, and characterized the evolution of the plastome in liverworts. The structure of the plastid genome is overall conserved across the phylogeny of liverworts, with only two structural variants detected from simple thalloids, besides 18 out of 43 liverwort genera showing intron variations in their plastomes. Complex thalloid liverworts maintain the most plastid genes, and seem to undergo fewer gene deletions and pseudogenization events than other liverworts. Plastid phylogenetic inferences yielded mostly robustly supported relationships, and consistently resolved Ptilidiales as the sister to Porellales. The relative ratio of silent substitutions across the three genetic compartments (i.e., 1:15:10, for mitochondrial:plastid:nuclear) suggests that liverwort plastid genes have the potential to evolve faster than their nuclear counterparts, unlike in any other major land plant lineages where the mutation rate of nuclear genes overwhelm those of their plastid and mitochondrial counterparts.

Insights Into the Evolutionary History of the Subfamily Orthotrichoideae (Orthotrichaceae, Bryophyta): New and Former Supra-Specific Taxa So Far Obscured by Prevailing Homoplasy

Mosses of the subfamily Orthotrichoideae represent one of the main components of the cryptogam epiphytic communities in temperate areas. During the last two decades, this taxonomical group has undergone an extensive revision that has led to its rearrangement at the generic level. However, their phylogenetic relationships and inferences on the evolutionary patterns that have driven the present diversity have little advanced. In this study, we present a dated molecular phylogenetic reconstruction at the subfamily level, including 130 samples that represent the 12 genera currently recognized within the subfamily, and the analysis of four molecular markers: ITS2, rps4, trnG, and trnL-F. We also analyze 13 morphological characters of systematic value to infer their origin and diagnostic utility within the subfamily. The phylogenetic reconstruction yields three main clades within the subfamily, two of which correspond to the tribe Zygodonteae, and one to Orthotricheae. Within Zygodonteae, the genus Zygodon results to be a polyphyletic artificial assembly, and we propose to separate a new genus named Australoria. Conversely, our results do not support the separation of Pentastichella and Pleurorthotrichum at the genus level and we therefore propose to include Pleurorthotrichum in Pentastichella. Regarding Orthotricheae, our analyses clearly allow the distinction of two subtribes: Orthotrichinae and Lewinskyinae. Within the latter, Ulota results a polyphyletic entity, and therefore we propose the segregation of a separate new genus named Atlantichella. Dating analyses allow us to conclude that the split of the tribes within Orthotrichoideae dates from the Middle Jurassic, while the diversification of Orthotrichum and Zygodon probably started during the Late Cretaceous. However, most of the extant genera of this subfamily seem to be younger, and apparently its highest diversification burst took place during the Oligocene. Finally, the analysis of the morphological traits reveals that most of the characters previously used to separate genera and here tested are homoplastic, which has hindered the taxonomical and systematic proposals for decades. However, even if there are no exclusive characters, all of the genera can be defined by the combination of a few characters.