The Physcomitrella patens chromosome-scale assembly reveals moss genome structure and evolution

The rDNA Loci-Intersections of Replication, Transcription, and Repair Pathways

Genes encoding ribosomal RNA (rDNA) are essential for cell survival and are particularly sensitive to factors leading to genomic instability. Their repetitive character makes them prone to inappropriate recombinational events arising from collision of transcriptional and replication machineries, resulting in unstable rDNA copy numbers. In this review, we summarize current knowledge on the structure and organization of rDNA, its role in sensing changes in the genome, and its linkage to aging. We also review recent findings on the main factors involved in chromatin assembly and DNA repair in the maintenance of rDNA stability in the model plants Arabidopsis thaliana and the moss Physcomitrella patens, providing a view across the plant evolutionary tree.

Mutations in the Physcomitrium patens gene encoding Aminodeoxychorismate Synthase confer auxotrophic phenotypes

To facilitate genetic mapping of developmental mutants of Physcomitrium patens, we produced a genetic marker that combines recessive auxotrophy with dominant positive selection. We first identified the gene affected by the pabB4 auxotrophic mutation and then replaced it with a cassette that confers antibiotic resistance. This strain may be used to produce bi-parental somatic hybrids with nearly any other strain.

ViralRecall-A Flexible Command-Line Tool for the Detection of Giant Virus Signatures in 'Omic Data

Giant viruses are widespread in the biosphere and play important roles in biogeochemical cycling and host genome evolution. Also known as nucleo-cytoplasmic large DNA viruses (NCLDVs), these eukaryotic viruses harbor the largest and most complex viral genomes known. Studies have shown that NCLDVs are frequently abundant in metagenomic datasets, and that sequences derived from these viruses can also be found endogenized in diverse eukaryotic genomes. The accurate detection of sequences derived from NCLDVs is therefore of great importance, but this task is challenging owing to both the high level of sequence divergence between NCLDV families and the extraordinarily high diversity of genes encoded in their genomes, including some encoding for metabolic or translation-related functions that are typically found only in cellular lineages. Here, we present ViralRecall, a bioinformatic tool for the identification of NCLDV signatures in ‘omic data. This tool leverages a library of giant virus orthologous groups (GVOGs) to identify sequences that bear signatures of NCLDVs. We demonstrate that this tool can effectively identify NCLDV sequences with high sensitivity and specificity. Moreover, we show that it can be useful both for removing contaminating sequences in metagenome-assembled viral genomes as well as the identification of eukaryotic genomic loci that derived from NCLDV. ViralRecall is written in Python 3.5 and is freely available on GitHub: https://github.com/faylward/viralrecall.

A Synthetic Approach to Reconstruct the Evolutionary and Functional Innovations of the Plant Histone Variant H2A.W

Diversification of histone variants is marked by the acquisition of distinct motifs and functional properties through convergent evolution.1-4 H2A variants are distinguished by specific C-terminal motifs and tend to be segregated within defined domains of the genome.5,6 Whether evolution of these motifs pre-dated the evolution of segregation mechanisms or vice versa has remained unclear. A suitable model to address this question is the variant H2A.W, which evolved in plants through acquisition of a KSPK motif7 and is tightly associated with heterochromatin.4 We used fission yeast, where chromatin is naturally devoid of H2A.W, to study the impact of engineered chimeras combining yeast H2A with the KSPK motif. Biochemical assays showed that the KSPK motif conferred nucleosomes with specific properties. Despite uniform incorporation of the engineered H2A chimeras in the yeast genome, the KSPK motif specifically affected heterochromatin composition and function. We conclude that the KSPK motif promotes chromatin properties in yeast that are comparable to the properties and function of H2A.W in plant heterochromatin. We propose that the selection of functional motifs confer histone variants with properties that impact primarily a specific chromatin state. The association between a new histone variant and a preferred chromatin state can thus provide a setting for the evolution of mechanisms that segregate the new variant to this state, thereby enhancing the impact of the selected properties of the variant on genome activity.

Genome-wide characterization and analysis of the CCT motif family genes in soybean (Glycine max)

MAIN CONCLUSION

Soybean possesses 19 CMF genes which mainly arose from duplication events. Their features and motifs are highly conserved but transcriptional data indicated functional diversity in metabolism and stress responses. CCT [for CONSTANS, CONSTANS-like (CO-like), and timing of CAB expression1 (TOC1)] domain-containing genes play important roles in regulating flowering, plant growth, and grain yield and are also involved in stress responses. The CMF (CCT motif family) genes, included in the CCT family, contain a single CCT domain as the only identifiable domain in their predicted protein sequence and are interesting targets for breeding programs. In this study, we identified 19 putative GmCMF genes, based on the latest soybean (Glycine max) genome annotation. The predicted GmCMF proteins were characterized based on conserved structural features, and a phylogenetic tree was constructed including all CMF proteins from rice and Arabidopsis as representative examples of the monocotyledonous (monocot) and dicotyledonous (dicot) plants, respectively. High similarities in the conserved motifs of the protein sequences and the gene structures were found. In addition, by analyzing the CMF gene family in soybean, we identified seven pairs of genes that originated from segmental chromosomal duplication events attributable to the most recent whole-genome duplication (WGD) event in the Glycine lineage. Expression analysis of GmCMF genes in various tissues and after specific treatments demonstrated tissue and stress-response specific differential expression. Gene expression analysis was complemented by the identification of putative cis-elements present in the promoter regions of the genes through a bioinformatics approach, using the existing soybean reference genome sequence and gene models. Co-functional networks inferred from distinct types of genomics data-including microarrays and RNA-seq samples from soybean-revealed that GmCMF genes might play crucial roles in metabolism and transport processes. The results of this study, the first systematic analysis of the soybean CCT gene family, can serve as a strong foundation for further elucidation of their physiological functions and biological roles.

Axenic in vitro cultivation of 19 peat moss (Sphagnum L.) species as a resource for basic biology, biotechnology, and paludiculture

Sphagnum farming can substitute peat with renewable biomass and thus help mitigate climate change. Large volumes of the required founder material can only be supplied sustainably by axenic cultivation in bioreactors. We established axenic in vitro cultures from sporophytes of 19 Sphagnum species collected in Austria, Germany, Latvia, the Netherlands, Russia, and Sweden: S. angustifolium, S. balticum, S. capillifolium, S. centrale, S. compactum, S. cuspidatum, S. fallax, S. fimbriatum, S. fuscum, S. lindbergii, S. medium/divinum, S. palustre, S. papillosum, S. rubellum, S. russowii, S. squarrosum, S. subnitens, S. subfulvum and S. warnstorfii. These species cover five of the six European Sphagnum subgenera; namely, Acutifolia, Cuspidata, Rigida, Sphagnum and Squarrosa. Their growth was measured in suspension cultures, whereas their ploidy was determined by flow cytometry and compared with the genome size of Physcomitrella patens. We identified haploid and diploid Sphagnum species, found that their cells are predominantly arrested in the G1 phase of the cell cycle, and did not find a correlation between plant productivity and ploidy. DNA barcoding was achieved by sequencing introns of the BRK1 genes. With this collection, high-quality founder material for diverse large-scale applications, but also for basic Sphagnum research, is available from the International Moss Stock Center.

How to Grow a Tree: Plant Voltage-Dependent Cation Channels in the Spotlight of Evolution

Phylogenetic analysis can be a powerful tool for generating hypotheses regarding the evolution of physiological processes. Here, we provide an updated view of the evolution of the main cation channels in plant electrical signalling: the Shaker family of voltage-gated potassium channels and the two-pore cation (K⁺) channel (TPC1) family. Strikingly, the TPC1 family followed the same conservative evolutionary path as one particular subfamily of Shaker channels (K_out) and remained highly invariant after terrestrialisation, suggesting that electrical signalling was, and remains, key to survival on land. We note that phylogenetic analyses can have pitfalls, which may lead to erroneous conclusions. To avoid these in the future, we suggest guidelines for analyses of ion channel evolution in plants.

TE-greedy-nester: structure-based detection of LTR retrotransposons and their nesting

Motivation

Transposable elements (TEs) in eukaryotes often get inserted into one another, forming sequences that become a complex mixture of full-length elements and their fragments. The reconstruction of full-length elements and the order in which they have been inserted is important for genome and transposon evolution studies. However, the accumulation of mutations and genome rearrangements over evolutionary time makes this process error-prone and decreases the efficiency of software aiming to recover all nested full-length TEs.

Results

We created software that uses a greedy recursive algorithm to mine increasingly fragmented copies of full-length LTR retrotransposons in assembled genomes and other sequence data. The software called TE-greedy-nester considers not only sequence similarity but also the structure of elements. This new tool was tested on a set of natural and synthetic sequences and its accuracy was compared to similar software. We found TE-greedy-nester to be superior in a number of parameters, namely computation time and full-length TE recovery in highly nested regions.

Availability and implementation

http://gitlab.fi.muni.cz/lexa/nested.

Supplementary information

Supplementary data are available at Bioinformatics online.

DNA methylation mutants in Physcomitrella patens elucidate individual roles of CG and non-CG methylation in genome regulation

Cytosine (DNA) methylation in plants regulates the expression of genes and transposons. While methylation in plant genomes occurs at CG, CHG, and CHH sequence contexts, the comparative roles of the individual methylation contexts remain elusive. Here, we present Physcomitrella patens as the second plant system, besides Arabidopsis thaliana, with viable mutants with an essentially complete loss of methylation in the CG and non-CG contexts. In contrast to A. thaliana, P. patens has more robust CHH methylation, similar CG and CHG methylation levels, and minimal cross-talk between CG and non-CG methylation, making it possible to study context-specific effects independently. Our data found CHH methylation to act in redundancy with symmetric methylation in silencing transposons and to regulate the expression of CG/CHG-depleted transposons. Specific elimination of CG methylation did not dysregulate transposons or genes. In contrast, exclusive removal of non-CG methylation massively up-regulated transposons and genes. In addition, comparing two exclusively but equally CG- or CHG-methylated genomes, we show that CHG methylation acts as a greater transcriptional regulator than CG methylation. These results disentangle the transcriptional roles of CG and non-CG, as well as symmetric and asymmetric methylation in a plant genome, and point to the crucial role of non-CG methylation in genome regulation.

Ties between Stress and Lipid Droplets Pre-date Seeds

Seeds were a key evolutionary innovation. These durable structures provide a concerted solution to two challenges on land: dispersal and stress. Lipid droplets (LDs) that act as nutrient storage reservoirs are one of the main cell-biological reasons for seed endurance. Although LDs are key structures in spermatophytes and are especially abundant in seeds, they are found across plants and algae, and increase during stress. Further, the proteins that underpin their form and function often have deep homologs. We propose an evolutionary scenario in which (i) the generation of LDs arose as a mechanism to mediate general drought and desiccation resilience, and (ii) the required protein framework was co-opted by spermatophytes for a seed-specific program.

A blind and independent benchmark study for detecting differeally methylated regions in plants

MOTIVATION

Bisulfite sequencing (BS-seq) is a state-of-the-art technique for investigating methylation of the DNA to gain insights into the epigenetic regulation. Several algorithms have been published for identification of differentially methylated regions (DMRs). However, the performances of the individual methods remain unclear and it is difficult to optimally select an algorithm in application settings.

RESULTS

We analyzed BS-seq data from four plants covering three taxonomic groups. We first characterized the data using multiple summary statistics describing methylation levels, coverage, noise, as well as frequencies, magnitudes and lengths of methylated regions. Then, simulated datasets with most similar characteristics to real experimental data were created. Seven different algorithms (metilene, methylKit, MOABS, DMRcate, Defiant, BSmooth, MethylSig) for DMR identification were applied and their performances were assessed. A blind and independent study design was chosen to reduce bias and to derive practical method selection guidelines. Overall, metilene had superior performance in most settings. Data attributes, such as coverage and spread of the DMR lengths, were found to be useful for selecting the best method for DMR detection. A decision tree to select the optimal approach based on these data attributes is provided. The presented procedure might serve as a general strategy for deriving algorithm selection rules tailored to demands in specific application settings.

AVAILABILITY AND IMPLEMENTATION

Scripts that were used for the analyses and that can be used for prediction of the optimal algorithm are provided at https://github.com/kreutz-lab/DMR-DecisionTree. Simulated and experimental data are available at https://doi.org/10.6084/m9.figshare.11619045.

CONTACT

ckreutz@imbi.uni-freiburg.de.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Evolutionary aspects of the Viridiplantae nitroreductases

Background

Nitroreductases are a family of evolutionarily related proteins catalyzing the reduction of nitro-substituted compounds. Nitroreductases are widespread enzymes, but nearly all modern research and practical application have been concentrated on the bacterial proteins, mainly nitroreductases of Escherichia coli. The main aim of this study is to describe the phylogenic distribution of the nitroreductases in the photosynthetic eukaryotes (Viridiplantae) to highlight their structural similarity and areas for future research and application.

Results

This study suggests that homologs of nitroreductase proteins are widely presented also in Viridiplantae. Maximum likelihood phylogenetic tree reconstruction method and comparison of the structural models suggest close evolutional relation between cyanobacterial and Viridiplantae nitroreductases.

Conclusions

This study provides the first attempt to understand the evolution of nitroreductase protein family in Viridiplantae. Our phylogeny estimation and preservation of the chloroplasts/mitochondrial localization indicate the evolutional origin of the plant nitroreductases from the cyanobacterial endosymbiont. A defined high level of the similarity on the structural level suggests conservancy also for the functions. Directions for the future research and industrial application of the Viridiplantae nitroreductases are discussed.

Tracking the evolutionary innovations of plant terrestrialization

The gradual transition of the algal ancestor from the freshwater to land has always attracted evolutionary biologists. The recent report of high-quality reference genomes of five Charophyta algae (Spirogloea muscicola, Mesotaenium endlicherianum, Mesostigma viride, Chlorokybus atmophyticus and Penium margaritaceum) and one hornwort (Anthoceros angustus) species sheds light on this fascinating transition. These early diverging plants and algae could have gained new genes from soil bacteria and fungi through horizontal gene transfer (HGT), which was so common during plant terrestrialization and may outrun our expectations. Through reviewing and critical thinking about the advancements on these plant genomes, here, I propose three prospective research directions that need to address in the future: (i) due to the ubiquitous nature of viruses that is similar to soil bacteria and fungi, there is less attention to viruses that probably also play an important role in the genome evolution of plants via HGT; (ii) multicellularity has occurred many times independently, but we still know a little about the biological and ecological mechanisms leading to multi-cellularity in Streptophyta; (iii) and most importantly, the quantitative relationships between genetic innovations and environmental variables such as temperature, precipitation and solar radiation, need pioneering research collaborated by biological evolutionists, computer scientists, and ecologists, which are crucial for understanding the macroevolution of plants and could also be used to simulate the evolution of plants under future climate change.

Fibrillarin evolution through the Tree of Life: Comparative genomics and microsynteny network analyses provide new insights into the evolutionary history of Fibrillarin

Fibrillarin (FIB), a methyltransferase essential for life in the vast majority of eukaryotes, is involved in methylation of rRNA required for proper ribosome assembly, as well as methylation of histone H2A of promoter regions of rRNA genes. RNA viral progression that affects both plants and animals requires FIB proteins. Despite the importance and high conservation of fibrillarins, there little is known about the evolutionary dynamics of this small gene family. We applied a phylogenomic microsynteny-network approach to elucidate the evolutionary history of FIB proteins across the Tree of Life. We identified 1063 non-redundant FIB sequences across 1049 completely sequenced genomes from Viruses, Bacteria, Archaea, and Eukarya. FIB is a highly conserved single-copy gene through Archaea and Eukarya lineages, except for plants, which have a gene family expansion due to paleopolyploidy and tandem duplications. We found a high conservation of the FIB genomic context during plant evolution. Surprisingly, FIB in mammals duplicated after the Eutheria split (e.g., ruminants, felines, primates) from therian mammals (e.g., marsupials) to form two main groups of sequences, the FIB and FIB-like groups. The FIB-like group transposed to another genomic context and remained syntenic in all the eutherian mammals. This transposition correlates with differences in the expression patterns of FIB-like proteins and with elevated Ks values potentially due to reduced evolutionary constraints of the duplicated copy. Our results point to a unique evolutionary event in mammals, between FIB and FIB-like genes, that led to non-redundant roles of the vital processes in which this protein is involved.

Metabolome analysis using multiple data mining approaches suggests luteolin biosynthesis in Physcomitrella patens

The model land plant Physcomitrella patens synthesizes flavonoids which may act as protectant of ultraviolet-B radiation. We aimed to uncover its flavonoid profile, for which metabolome analysis using liquid chromatography coupled with Ion trap/Orbitrap mass spectrometry was performed. From the 80% methanol extracts, 661 valid peaks were detected. Prediction of the elemental compositions within a mass accuracy of 2 ppm indicated that 217 peaks had single elemental composition. A compound database search revealed 47 peaks to be annotated as secondary metabolites based on the compound database search. Comprehensive substituent search by ShiftedIonsFinder showed there were 13 peaks of potential flavonoid derivatives. Interestingly, a peak having m/z 287.0551, corresponding to that of luteolin, was detected, even though flavone synthase has never been identified in P. patens. Using P. patens labeled with stable isotopes (¹³C-, ¹⁵N-, ¹⁸O-, and ³⁴S), we confirmed the elemental composition of the peak as C₁₅H₁₀O₆. By a comparison of MS/MS spectra with that of authentic standard, the peak was identified as luteolin or related flavone isomers. This is the first report of luteolin or related flavones synthesis and the possibility of the existence of an unknown enzyme with flavone synthase activity in P. patens.

High gene space divergence contrasts with frozen vegetative architecture in the moss family Funariaceae

A new paradigm has slowly emerged regarding the diversification of bryophytes, with inferences from molecular data highlighting a dynamic evolution of their genome. However, comparative studies of expressed genes among closely related taxa is so far missing. Here we contrast the dimensions of the vegetative transcriptome of Funaria hygrometrica and Physcomitrium pyriforme against the genome of their relative, Physcomitrium (Physcomitrella) patens. These three species of Funariaceae share highly conserved vegetative bodies, and are partially sympatric, growing on mineral soil in mostly temperate regions. We analyzed the vegetative gametophytic transcriptome of F. hygrometrica and P. pyriforme and mapped short reads, transcripts, and proteins to the genome and gene space of P. patens. Only about half of the transcripts of F. hygrometrica map to their ortholog in P. patens, whereas at least 90% of those of P. pyriforme align to loci in P. patens. Such divergence is unexpected given the high morphological similarity of the gametophyte but reflects the estimated times of divergence of F. hygrometrica and P. pyriforme from P. patens, namely 55 and 20 mya, respectively. The newly sampled transcriptomes bear signatures of at least one, rather ancient, whole genome duplication (WGD), which may be shared with one reported for P. patens. The transcriptomes of F. hygrometrica and P. pyriforme reveal significant contractions or expansions of different gene families. While transcriptomes offer only an incomplete estimate of the gene space, the high number of transcripts obtained suggest a significant divergence in gene sequences, and gene number among the three species, indicative of a rather strong, dynamic genome evolution, shaped in part by whole, partial or localized genome duplication. The gene ontology of their specific and rapidly-evolving protein families, suggests that the evolution of the Funariaceae may have been driven by the diversification of metabolic genes that may optimize the adaptations to environmental conditions, a hypothesis well in line with ecological patterns in the genetic diversity and structure in seed plants.

Reconstructing trait evolution in plant evo-devo studies

Our planet is teeming with an astounding diversity of plants. In a mere single group of closely related species, tremendous diversity can be observed in their form and function - the colour of petals in flowering plants, the shape of the fronds in ferns, and the branching pattern of the gametophyte in mosses. Diversity can also be found in subtler traits, such as the resistance to pathogens or the ability to recruit symbiotic microbes from the environment. Plant traits can also be highly conserved - at the cellular and metabolic levels, entire biosynthetic pathways are present in all plant groups, and morphological characteristics such as vascular tissues have been conserved for hundreds of millions of years. The research community that seeks to understand these traits - both the diverse and the conserved - by taking an evolutionary point-of-view on plant biology is growing. Here, we summarize a subset of the different aspects of plant evolutionary biology, provide a guide for structuring comparative biology approaches and discuss the pitfalls that (plant) researchers should avoid when embarking on such studies.

DNA damage triggers reprogramming of differentiated cells into stem cells in Physcomitrella

DNA damage can result from intrinsic cellular processes and from exposure to stressful environments. Such DNA damage generally threatens genome integrity and cell viability¹. However, here we report that the transient induction of DNA strand breaks (single-strand breaks, double-strand breaks or both) in the moss Physcomitrella patens can trigger the reprogramming of differentiated leaf cells into stem cells without cell death. After intact leafy shoots (gametophores) were exposed to zeocin, an inducer of DNA strand breaks, the STEM CELL-INDUCING FACTOR 1 (STEMIN1)² promoter was activated in some leaf cells. These cells subsequently initiated tip growth and underwent asymmetric cell divisions to form chloronema apical stem cells, which are in an earlier phase of the life cycle than leaf cells and have the ability to form new gametophores. This DNA-strand-break-induced reprogramming required the DNA damage sensor ATR kinase, but not ATM kinase, together with STEMIN1 and closely related proteins. ATR was also indispensable for the induction of STEMIN1 by DNA strand breaks. Our findings indicate that DNA strand breaks, which are usually considered to pose a severe threat to cells, trigger cellular reprogramming towards stem cells via the activity of ATR and STEMINs.

Pterocarpan synthase (PTS) structures suggest a common quinone methide-stabilizing function in dirigent proteins and proteins with dirigent-like domains

The biochemical activities of dirigent proteins (DPs) give rise to distinct complex classes of plant phenolics. DPs apparently began to emerge during the aquatic-to-land transition, with phylogenetic analyses revealing the presence of numerous DP subfamilies in the plant kingdom. The vast majority (>95%) of DPs in these large multigene families still await discovery of their biochemical functions. Here, we elucidated the 3D structures of two pterocarpan-forming proteins with dirigent-like domains. Both proteins stereospecifically convert distinct diastereomeric chiral isoflavonoid precursors to the chiral pterocarpans, (-)- and (+)-medicarpin, respectively. Their 3D structures enabled comparisons with stereoselective lignan- and aromatic terpenoid-forming DP orthologs. Each protein provides entry into diverse plant natural products classes, and our experiments suggest a common biochemical mechanism in binding and stabilizing distinct plant phenol-derived mono- and bis-quinone methide intermediates during different C-C and C-O bond-forming processes. These observations provide key insights into both their appearance and functional diversification of DPs during land plant evolution/adaptation. The proposed biochemical mechanisms based on our findings provide important clues to how additional physiological roles for DPs and proteins harboring dirigent-like domains can now be rationally and systematically identified.

Contiguous and stochastic CHH methylation patterns of plant DRM2 and CMT2 revealed by single-read methylome analysis

Cytosine methylome data is commonly generated through next-generation sequencing, with analyses averaging methylation states of individual reads. We propose an alternative method of analysing single-read methylome data. Using this method, we identify patterns relating to the mechanism of two plant non-CG-methylating enzymes, CMT2 and DRM2. CMT2-methylated regions show higher stochasticity, while DRM2-methylated regions have higher variation among cells. Based on these patterns, we develop a classifier that predicts enzyme activity in different species and tissues. To facilitate further single-read analyses, we develop a genome browser, SRBrowse, optimised for visualising and analysing sequencing data at single-read resolution.

Heat stress response in the closest algal relatives of land plants reveals conserved stress signaling circuits

All land plants (embryophytes) share a common ancestor that likely evolved from a filamentous freshwater alga. Elucidating the transition from algae to embryophytes - and the eventual conquering of Earth's surface - is one of the most fundamental questions in plant evolutionary biology. Here, we investigated one of the organismal properties that might have enabled this transition: resistance to drastic temperature shifts. We explored the effect of heat stress in Mougeotia and Spirogyra, two representatives of Zygnematophyceae - the closest known algal sister lineage to land plants. Heat stress induced pronounced phenotypic alterations in their plastids, and high-performance liquid chromatography-tandem mass spectroscopy-based profiling of 565 transitions for the analysis of main central metabolites revealed significant shifts in 43 compounds. We also analyzed the global differential gene expression responses triggered by heat, generating 92.8 Gbp of sequence data and assembling a combined set of 8905 well-expressed genes. Each organism had its own distinct gene expression profile; less than one-half of their shared genes showed concordant gene expression trends. We nevertheless detected common signature responses to heat such as elevated transcript levels for molecular chaperones, thylakoid components, and - corroborating our metabolomic data - amino acid metabolism. We also uncovered the heat-stress responsiveness of genes for phosphorelay-based signal transduction that links environmental cues, calcium signatures and plastid biology. Our data allow us to infer the molecular heat stress response that the earliest land plants might have used when facing the rapidly shifting temperature conditions of the terrestrial habitat.

The identification of differentially expressed genes in male and female gametophytes of simple thalloid liverwort Pellia endiviifolia sp. B using an RNA-seq approach

MAIN CONCLUSION

This study shows differences in gene expression between male and female gametophytes of the simple thalloid liverwort with a distinction between the vegetative and reproductive phases of growth. Pellia endiviifolia is a simple thalloid liverwort that, together with hornworts and mosses, represents the oldest living land plants. The limited taxon sampling for genomic and functional studies hampers our understanding of processes governing evolution of these plants. RNA sequencing represents an attractive way to elucidate the molecular mechanisms of non-model species development. In the present study, RNA-seq was used to profile the differences in gene expression between P. endiviifolia male and female gametophytes, with a distinction between the vegetative and reproductive phases of growth. By comparison of the gene expression profiles from individuals producing sex organs with the remaining thalli types, we have determined a set of genes whose expression might be important for the development of P. endiviifolia reproductive organs. The selected differentially expressed genes (DEGs) were categorized into five main pathways: metabolism, genetic information processing, environmental information processing, cellular processes, and organismal systems. A comparison of the obtained data with the Marchantia polymorpha transcriptome resulted in the identification of genes exhibiting a similar expression pattern during the reproductive phase of growth between members of the two distinct liverwort classes. The common expression profile of  87 selected genes suggests a common mechanism governing sex organ development in both liverwort species. The obtained RNA-seq results were confirmed by RT-qPCR for the DEGs with the highest differences in expression level. Five Pellia-female-specific and two Pellia-male-specific DEGs showed enriched expression in archegonia and antheridia, respectively. The identified genes are promising candidates for functional studies of their involvement in liverwort sexual reproduction.

Single Nucleotide Polymorphism Charting of P. patens Reveals Accumulation of Somatic Mutations During in vitro Culture on the Scale of Natural Variation by Selfing

Introduction: Physcomitrium patens (Hedw.) Mitten (previously known as Physcomitrella patens) was collected by H.L.K. Whitehouse in Gransden Wood (Huntingdonshire, United Kingdom) in 1962 and distributed across the globe starting in 1974. Hence, the Gransden accession has been cultured in vitro in laboratories for half a century. Today, there are more than 13 different pedigrees derived from the original accession. Additionally, accessions from other sites worldwide were collected during the last decades. Methods and Results: In this study, 250 high throughput RNA sequencing (RNA-seq) samples and 25 gDNA samples were used to detect single nucleotide polymorphisms (SNPs). Analyses were performed using five different P. patens accessions and 13 different Gransden pedigrees. SNPs were overlaid with metadata and known phenotypic variations. Unique SNPs defining Gransden pedigrees and accessions were identified and experimentally confirmed. They can be successfully employed for PCR-based identification. Conclusion: We show independent mutations in different Gransden laboratory pedigrees, demonstrating that somatic mutations occur and accumulate during in vitro culture. The frequency of such mutations is similar to those observed in naturally occurring populations. We present evidence that vegetative propagation leads to accumulation of deleterious mutations, and that sexual reproduction purges those. Unique SNP sets for five different P. patens accessions were isolated and can be used to determine individual accessions as well as Gransden pedigrees. Based on that, laboratory methods to easily determine P. patens accessions and Gransden pedigrees are presented.

The Importance of ATM and ATR in Physcomitrella patens DNA Damage Repair, Development, and Gene Targeting

Coordinated by ataxia-telangiectasia-mutated (ATM) and ATM and Rad3-related (ATR), two highly conserved kinases, DNA damage repair ensures genome integrity and survival in all organisms. The Arabidopsis thaliana (A. thaliana) orthologues are well characterized and exhibit typical mammalian characteristics. We mutated the Physcomitrellapatens (P. patens) PpATM and PpATR genes by deleting functionally important domains using gene targeting. Both mutants showed growth abnormalities, indicating that these genes, particularly PpATR, are important for normal vegetative development. ATR was also required for repair of both direct and replication-coupled double-strand breaks (DSBs) and dominated the transcriptional response to direct DSBs, whereas ATM was far less important, as shown by assays assessing resistance to DSB induction and SuperSAGE-based transcriptomics focused on DNA damage repair genes. These characteristics differed significantly from the A. thaliana genes but resembled those in yeast (Saccharomyces cerevisiae). PpATR was not important for gene targeting, pointing to differences in the regulation of gene targeting and direct DSB repair. Our analysis suggests that ATM and ATR functions can be substantially diverged between plants. The differences in ATM and ATR reflect the differences in DSB repair pathway choices between A. thaliana and P. patens, suggesting that they represent adaptations to different demands for the maintenance of genome stability.

The Ceratodon purpureus genome uncovers structurally complex, gene rich sex chromosomes.. [PMID: 0 DOI: 10.1101/2020.07.03.163634]

Non-recombining sex chromosomes, like the mammalian Y, often lose genes and accumulate transposable elements, a process termed degeneration1,2. The correlation between suppressed recombination and degeneration is clear in animal XY systems1,2, but the absence of recombination is confounded with other asymmetries between the X and Y. In contrast, UV sex chromosomes, like those found in bryophytes, experience symmetrical population genetic conditions3,4. Here we test for degeneration in the bryophyte UV sex chromosome system through genomic comparisons with new female and male chromosome-scale reference genomes of the moss Ceratodon purpureus. We show that the moss sex chromosomes evolved over 300 million years ago and expanded via two chromosomal fusions. Although the sex chromosomes show signs of weaker purifying selection than autosomes, we find suppressed recombination alone is insufficient to drive gene loss on sex-specific chromosomes. Instead, the U and V sex chromosomes harbor thousands of broadly-expressed genes, including numerous key regulators of sexual development across land plants.

PHYTOCHROME INTERACTING FACTORs in the moss Physcomitrella patens regulate light-controlled gene expression

Phytochromes are red and far-red light receptors in plants that control growth and development in response to changes in the environment. Light-activated phytochromes enter the nucleus and act on a set of downstream signalling components to regulate gene expression. PHYTOCHROME INTERACTING FACTORs (PIFs) belong to the basic helix-loop-helix family of transcription factors and directly bind to light-activated phytochromes. Potential homologues of PIFs have been identified in ferns, bryophytes and streptophyte algae, and it has been shown that the potential PIF homologues from Physcomitrella patens, PIF1 to PIF4, have PIF function when expressed in Arabidopsis. However, their function in Physcomitrella is still unknown. Seed plant PIFs bind to G-box-containing promoters and, therefore, we searched the Physcomitrella genome for genes that contain G-boxes in their promoter. Here, we show that Physcomitrella PIFs activate these promoters in a G-box-dependent manner, suggesting that they could be direct PIF targets. Furthermore, we generated Physcomitrella pif1, pif2, pif3 and pif4 knock out mutant lines and quantified the expression of potential PIF direct target genes. The expression of these genes was generally reduced in pif mutants compared to the wildtype, but for several genes, the relative induction upon a short light treatment was higher in pif mutants than the wildtype. In contrast, expression of these genes was strongly repressed in continuous light, and pif mutants showed partial downregulation of these genes in the dark. Thus, the overall function of PIFs in light-regulated gene expression might be an ancient property of PIFs.

Characterisation of evolutionarily conserved key players affecting eukaryotic flagellar motility and fertility using a moss model

Defects in flagella/cilia are often associated with infertility and disease. Motile male gametes (sperm cells) are an ancestral eukaryotic trait that has been lost in several lineages like flowering plants. Here, we made use of a phenotypic male fertility difference between two moss (Physcomitrella patens) ecotypes to explore spermatozoid function. We compare genetic and epigenetic variation as well as expression profiles between the Gransden and Reute ecotype to identify a set of candidate genes associated with moss male infertility. We generated a loss-of-function mutant of a coiled-coil domain containing 39 (ccdc39) gene that is part of the flagellar hydin network. Defects in mammal and algal homologues of this gene coincide with a loss of fertility, demonstrating the evolutionary conservation of flagellar function related to male fertility across kingdoms. The Ppccdc39 mutant resembles the Gransden phenotype in terms of male fertility. Potentially, several somatic (epi-)mutations occurred during prolonged vegetative propagation of Gransden, causing regulatory differences of for example the homeodomain transcription factor BELL1. Probably these somatic changes are causative for the observed male fertility defect. We propose that moss spermatozoids might be employed as an easily accessible system to study male infertility of humans and animals in terms of flagellar structure and movement.

The evolving role of abscisic acid in cell function and plant development over geological time

Abscisic acid (ABA) is found in a wide diversity of organisms, yet we know most about the hormonal action of this compound in the ecologically dominant and economically important angiosperms. In angiosperms, ABA regulates a suite of critical responses from desiccation tolerance through to seed dormancy and stomatal closure. Work exploring the function of key genes in the ABA signalling pathway of angiosperms has revealed that this signal transduction pathway is ancient, yet considerable change in the physiological roles of this hormone have occurred over geological time. With recent advances in our capacity to characterise gene function in non-angiosperms we are on the cusp of revealing the origins of this critical hormonal signalling pathway in plants, and understanding how a simple hormone may have shaped land plant diversity, ecology and adaptation over the past 500 million years.

Step by step evolution of Indeterminate Domain (IDD) transcriptional regulators: from algae to angiosperms

INTRODUCTION

The Indeterminate Domain (IDD) proteins are a plant-specific subclass of C2H2 Zinc Finger transcription factors. Some of these transcription factors play roles in diverse aspects of plant metabolism and development, but the function of most of IDD genes is unknown and the molecular evolution of the subfamily has not been explored in detail.

METHODS

In this study, we mined available genome sequences of green plants (Viridiplantae) to reconstruct the phylogeny and then described the motifs/expression patterns of IDD genes.

KEY RESULTS

We identified the complete set of IDD genes of 16 Streptophyta genomes. We found that IDD and its sister clade STOP arose by a duplication at the base of Streptophyta. Once on land, the IDD genes duplicated extensively, giving rise to at least ten lineages. Some of these lineages were lost in extant non-vascular plants and gymnosperms, but all of them were retained in angiosperms, duplicating profoundly in dicots and monocots and acquiring, at the same time, surprising heterogeneity in their C-terminal regions and expression patterns.

CONCLUSIONS

IDDs were present in the last common ancestor of Streptophyta. On land, IDDs duplicated extensively, leading to ten lineages. Later, IDDs were recruited by angiosperms where they diversified greatly in number, C-terminal and expression patterns. Interestingly, such diversification occurred during the evolution of novel traits of the plant body. This study provides a solid framework of the orthology relationships of green land plant IDD transcription factors, thus increasing the accuracy of orthologue identification in model and non-model species and facilitating the identification of agronomically important genes related to plant metabolism and development.

Disequilibrium evolution of the Fructose-1,6-bisphosphatase gene family leads to their functional biodiversity in Gossypium species

Background

Fructose-1,6-bisphosphatase (FBP) is a key enzyme in the plant sucrose synthesis pathway, in the Calvin cycle, and plays an important role in photosynthesis regulation in green plants. However, no systemic analysis of FBPs has been reported in Gossypium species.

Results

A total of 41 FBP genes from four Gossypium species were identified and analyzed. These FBP genes were sorted into two groups and seven subgroups. Results revealed that FBP family genes were under purifying selection pressure that rendered FBP family members as being conserved evolutionarily, and there was no tandem or fragmental DNA duplication in FBP family genes. Collinearity analysis revealed that a FBP gene was located in a translocated DNA fragment and the whole FBP gene family was under disequilibrium evolution that led to a faster evolutionary progress of the members in G. barbadense and in A_t subgenome than those in other Gossypium species and in the D_t subgenome, respectively, in this study. Through RNA-seq analyses and qRT-PCR verification, different FBP genes had diversified biological functions in cotton fiber development (two genes in 0 DPA and 1DPA ovules and four genes in 20–25 DPA fibers), in plant responses to Verticillium wilt onset (two genes) and to salt stress (eight genes).

Conclusion

The FBP gene family displayed a disequilibrium evolution pattern in Gossypium species, which led to diversified functions affecting not only fiber development, but also responses to Verticillium wilt and salt stress. All of these findings provide the foundation for further study of the function of FBP genes in cotton fiber development and in environmental adaptability.

Phylogenomic synteny network analyses reveal ancestral transpositions of auxin response factor genes in plants

BACKGROUND

Auxin response factors (ARFs) have long been a research focus and represent a class of key regulators of plant growth and development. Integrated phylogenomic synteny network analyses were able to provide novel insights into the evolution of the ARF gene family.

RESULTS

Here, more than 3500 ARFs collected from plant genomes and transcriptomes covering major streptophyte lineages were used to reconstruct the broad-scale family phylogeny, where the early origin and diversification of ARF in charophytes was delineated. Based on the family phylogeny, we proposed a unified six-group classification system for angiosperm ARFs. Phylogenomic synteny network analyses revealed the deeply conserved genomic syntenies within each of the six ARF groups and the interlocking syntenic relationships connecting distinct groups. Recurrent duplication events, such as those that occurred in seed plants, angiosperms, core eudicots and grasses contributed to the expansion of ARF genes which facilitated functional diversification. Ancestral transposition activities in important plant families, including crucifers, legumes and grasses, were unveiled by synteny network analyses. Ancestral gene duplications along with transpositions have profound evolutionary significance which may have accelerated the functional diversification process of paralogues.

CONCLUSIONS

The broad-scale family phylogeny in combination with the state-of-art phylogenomic synteny network analyses not only allowed us to infer the evolutionary trajectory of ARF genes across distinct plant lineages, but also facilitated to generate a more robust classification regime for this transcription factor family. Our study provides insights into the evolution of ARFs which will enhance our current understanding of this important transcription factor family.

DEK1 displays a strong subcellular polarity during Physcomitrella patens 3D growth

Defective Kernel 1 (DEK1) is genetically at the nexus of the 3D morphogenesis of land plants. We aimed to localize DEK1 in the moss Physcomitrella patens to decipher its function during this process. To detect DEK1 in vivo, we inserted the tdTomato fluorophore into PpDEK1 gene locus. Confocal microscopy coupled with the use of time-gating allowed the precise DEK1 subcellular localization during 3D morphogenesis. DEK1 localization displays a strong polarized signal, as it is restricted to the plasma membrane domain between recently divided cells during the early steps of 3D growth development as well as during the subsequent vegetative growth. The signal furthermore displays a clear developmental pattern because it is only detectable in recently divided and elongating cells. Additionally, DEK1 localization appears to be independent of its calpain domain proteolytic activity. The DEK1 polar subcellular distribution in 3D tissue developing cells defines a functional cellular framework to explain its role in this developmental phase. Also, the observation of DEK1 during spermatogenesis suggests another biological function for this protein in plants. Finally the DEK1-tagged strain generated here provides a biological platform upon which further investigations into 3D developmental processes can be performed.

Evolution of Plant Hormone Response Pathways

This review focuses on the evolution of plant hormone signaling pathways. Like the chemical nature of the hormones themselves, the signaling pathways are diverse. Therefore, we focus on a group of hormones whose primary perception mechanism involves an Skp1/Cullin/F-box-type ubiquitin ligase: auxin, jasmonic acid, gibberellic acid, and strigolactone. We begin with a comparison of the core signaling pathways of these four hormones, which have been established through studies conducted in model organisms in the Angiosperms. With the advent of next-generation sequencing and advanced tools for genetic manipulation, the door to understanding the origins of hormone signaling mechanisms in plants beyond these few model systems has opened. For example, in-depth phylogenetic analyses of hormone signaling components are now being complemented by genetic studies in early diverging land plants. Here we discuss recent investigations of how basal land plants make and sense hormones. Finally, we propose connections between the emergence of hormone signaling complexity and major developmental transitions in plant evolution.

Origin and Diversity of Plant Receptor-Like Kinases

Because of their high level of diversity and complex evolutionary histories, most studies on plant receptor-like kinase subfamilies have focused on their kinase domains. With the large amount of genome sequence data available today, particularly on basal land plants and Charophyta, more attention should be paid to primary events that shaped the diversity of the RLK gene family. We thus focus on the motifs and domains found in association with kinase domains to illustrate their origin, organization, and evolutionary dynamics. We discuss when these different domain associations first occurred and how they evolved, based on a literature review complemented by some of our unpublished results.

PEATmoss (Physcomitrella Expression Atlas Tool): a unified gene expression atlas for the model plant Physcomitrella patens

Physcomitrella patens is a bryophyte model plant that is often used to study plant evolution and development. Its resources are of great importance for comparative genomics and evo-devo approaches. However, expression data from Physcomitrella patens were so far generated using different gene annotation versions and three different platforms: CombiMatrix and NimbleGen expression microarrays and RNA sequencing. The currently available P. patens expression data are distributed across three tools with different visualization methods to access the data. Here, we introduce an interactive expression atlas, Physcomitrella Expression Atlas Tool (PEATmoss), that unifies publicly available expression data for P. patens and provides multiple visualization methods to query the data in a single web-based tool. Moreover, PEATmoss includes 35 expression experiments not previously available in any other expression atlas. To facilitate gene expression queries across different gene annotation versions, and to access P. patens annotations and related resources, a lookup database and web tool linked to PEATmoss was implemented. PEATmoss can be accessed at https://peatmoss.online.uni-marburg.de.

Biological and Cellular Functions of the Microdomain-Associated FWL/CNR Protein Family in Plants

Membrane microdomains/nanodomains are sub-compartments of the plasma membrane enriched in sphingolipids and characterized by their unique protein composition. They play important roles in regulating plant development and plant-microbe interactions including mutualistic symbiotic interactions. Several protein families are associated with the microdomain fraction of biological membranes such as flotillins, prohibitins, and remorins. More recently, GmFWL1, a FWL/CNR protein exclusively expressed in the soybean nodule, was functionally characterized as a new microdomain-associated protein. Interestingly, GmFWL1 is homologous to the tomato FW2-2 protein, a major regulator of tomato fruit development. In this review, we summarize the knowledge gained about the biological, cellular, and physiological functions of members of the FWL/CNR family across various plant species. The role of the FWL/CNR proteins is also discussed within the scope of their evolution and transcriptional regulation.

Mosses: Versatile plants for biotechnological applications

Mosses have long been recognized as powerful experimental tools for the elucidation of complex processes in plant biology. Recent increases in the availability of sequenced genomes and mutant collections, the establishment of novel technologies for targeted mutagenesis, and the development of viable protocols for large-scale production in bioreactors are now transforming mosses into one of the most versatile tools for biotechnological applications. In the present review, we highlight the astonishing biotechnological potential of mosses and how these plants are being exploited for industrial, pharmaceutical, and environmental applications. We focus on the biological features that support their use as model organisms for basic and applied research, and how these are being leveraged to explore the biotechnological potential in an increasing number of species. Finally, we also provide an overview of the available moss cultivation protocols from an industrial perspective, offering insights into batch operations that are not yet well established or do not even exist in the literature. Our goal is to bolster the use of mosses as factories for the biosynthesis of molecules of interest and to show how these species can be harnessed for the generation of novel and commercially useful bioproducts.

Two ANGUSTIFOLIA genes regulate gametophore and sporophyte development in Physcomitrella patens

In Arabidopsis thaliana the ANGUSTIFOLIA (AN) gene regulates the width of leaves by controlling the diffuse growth of leaf cells in the medio-lateral direction. In the genome of the moss Physcomitrella patens, we found two normal ANs (PpAN1-1 and 1-2). Both PpAN1 genes complemented the A. thaliana an-1 mutant phenotypes. An analysis of spatiotemporal promoter activity of each PpAN1 gene, using transgenic lines that contained each PpAN1-promoter- uidA (GUS) gene, showed that both promoters are mainly active in the stems of haploid gametophores and in the middle to basal region of the young sporophyte that develops into the seta and foot. Analyses of the knockout lines for PpAN1-1 and PpAN1-2 genes suggested that these genes have partially redundant functions and regulate gametophore height by controlling diffuse cell growth in gametophore stems. In addition, the seta and foot were shorter and thicker in diploid sporophytes, suggesting that cell elongation was reduced in the longitudinal direction, whereas no defects were detected in tip-growing protonemata. These results indicate that both PpAN1 genes in P. patens function in diffuse growth of the haploid and diploid generations but not in tip growth. To visualize microtubule distribution in gametophore cells of P. patens, transformed lines expressing P. patens α-tubulin fused to sGFP were generated. Contrary to expectations, the orientation of microtubules in the tips of gametophores in the PpAN1-1/1-2 double-knockout lines was unchanged. The relationships among diffuse cell growth, cortical microtubules and AN proteins are discussed.

A pseudomolecule-scale genome assembly of the liverwort Marchantia polymorpha

Marchantia polymorpha has recently become a prime model for cellular, evo-devo, synthetic biological, and evolutionary investigations. We present a pseudomolecule-scale assembly of the M. polymorpha genome, making comparative genome structure analysis and classical genetic mapping approaches feasible. We anchored 88% of the M. polymorpha draft genome to a high-density linkage map resulting in eight pseudomolecules. We found that the overall genome structure of M. polymorpha is in some respects different from that of the model moss Physcomitrella patens. Specifically, genome collinearity between the two bryophyte genomes and vascular plants is limited, suggesting extensive rearrangements since divergence. Furthermore, recombination rates are greatest in the middle of the chromosome arms in M. polymorpha like in most vascular plant genomes, which is in contrast with P. patens where recombination rates are evenly distributed along the chromosomes. Nevertheless, some other properties of the genome are shared with P. patens. As in P. patens, DNA methylation in M. polymorpha is spread evenly along the chromosomes, which is in stark contrast with the angiosperm model Arabidopsis thaliana, where DNA methylation is strongly enriched at the centromeres. Nevertheless, DNA methylation and recombination rate are anticorrelated in all three species. Finally, M. polymorpha and P. patens centromeres are of similar structure and marked by high abundance of retroelements unlike in vascular plants. Taken together, the highly contiguous genome assembly we present opens unexplored avenues for M. polymorpha research by linking the physical and genetic maps, making novel genomic and genetic analyses, including map-based cloning, feasible.

Anthoceros genomes illuminate the origin of land plants and the unique biology of hornworts

Hornworts comprise a bryophyte lineage that diverged from other extant land plants >400 million years ago and bears unique biological features, including a distinct sporophyte architecture, cyanobacterial symbiosis and a pyrenoid-based carbon-concentrating mechanism (CCM). Here, we provide three high-quality genomes of Anthoceros hornworts. Phylogenomic analyses place hornworts as a sister clade to liverworts plus mosses with high support. The Anthoceros genomes lack repeat-dense centromeres as well as whole-genome duplication, and contain a limited transcription factor repertoire. Several genes involved in angiosperm meristem and stomatal function are conserved in Anthoceros and upregulated during sporophyte development, suggesting possible homologies at the genetic level. We identified candidate genes involved in cyanobacterial symbiosis and found that LCIB, a Chlamydomonas CCM gene, is present in hornworts but absent in other plant lineages, implying a possible conserved role in CCM function. We anticipate that these hornwort genomes will serve as essential references for future hornwort research and comparative studies across land plants.

Chromatin Organization in Early Land Plants Reveals an Ancestral Association between H3K27me3, Transposons, and Constitutive Heterochromatin

Genome packaging by nucleosomes is a hallmark of eukaryotes. Histones and the pathways that deposit, remove, and read histone modifications are deeply conserved. Yet, we lack information regarding chromatin landscapes in extant representatives of ancestors of the main groups of eukaryotes, and our knowledge of the evolution of chromatin-related processes is limited. We used the bryophyte Marchantia polymorpha, which diverged from vascular plants circa 400 mya, to obtain a whole chromosome genome assembly and explore the chromatin landscape and three-dimensional genome organization in an early diverging land plant lineage. Based on genomic profiles of ten chromatin marks, we conclude that the relationship between active marks and gene expression is conserved across land plants. In contrast, we observed distinctive features of transposons and other repetitive sequences in Marchantia compared with flowering plants. Silenced transposons and repeats did not accumulate around centromeres. Although a large fraction of constitutive heterochromatin was marked by H3K9 methylation as in flowering plants, a significant proportion of transposons were marked by H3K27me3, which is otherwise dedicated to the transcriptional repression of protein-coding genes in flowering plants. Chromatin compartmentalization analyses of Hi-C data revealed that repressed B compartments were densely decorated with H3K27me3 but not H3K9 or DNA methylation as reported in flowering plants. We conclude that, in early plants, H3K27me3 played an essential role in heterochromatin function, suggesting an ancestral role of this mark in transposon silencing.

Dynamics of the Transcriptome Response to Heat in the Moss, Physcomitrella patens

Thermal stress negatively impacts crop yields, and as the overall temperature of the earth's atmosphere is gradually increasing, the identification of the temperature transduction pathway of the heat signal is essential in developing new strategies in order to adapt plant breeding to warmer climates. Heat stress damages the molecular structures and physiological processes in plants in proportion to the level and duration of the stress, which leads to different types of responses. In general, plants respond more efficiently when they are first subjected to a moderate temperature increase before being subjected to a higher temperature stress. This adaptive response is called the acclimation period and has been investigated in several plant species. However, there is a lack of information on the dynamic of the Heat Shock Response (HSR) over a continuous period of temperature rise without an acclimation period. In this paper, we investigated the effects of mild (30 °C) and high (37 °C) continuous heat stress over a 24-h period. Through RNA-Seq analysis, we assessed the remodeling of the transcriptome in the moss Physcomitrella patens. Our results showed that the 30 °C treatment particularly affected the expression of a few genes at 1 and 24 h, suggesting a biphasic response. Up-regulated genes at 1 h encode mainly HSR proteins (protein folding and endoplasmic reticulum stress), indicating an early heat response; while the up-regulated genes at 24 h belong to the thiamine biosynthesis pathway. In contrast, the genes involved in photosynthesis and carbon partitioning were repressed by this treatment. Under a higher temperature stress (37 °C), the induction of the HSR occurred rapidly (1 h) and was then attenuated throughout the time points investigated. A network approach (Weighted Gene Correlation Network Analysis, WGCNA) was used to identify the groups of genes expressing similar profiles, highlighting a HsfA1E binding motif within the promoters of some unrelated genes which displayed rapid and transient heat-activation. Therefore, it could be suggested that these genes could be direct targets of activation by a HsfA1E transcription factors.

The hornwort genome and early land plant evolution

Hornworts, liverworts and mosses are three early diverging clades of land plants, and together comprise the bryophytes. Here, we report the draft genome sequence of the hornwort Anthoceros angustus. Phylogenomic inferences confirm the monophyly of bryophytes, with hornworts sister to liverworts and mosses. The simple morphology of hornworts correlates with low genetic redundancy in plant body plan, while the basic transcriptional regulation toolkit for plant development has already been established in this early land plant lineage. Although the Anthoceros genome is small and characterized by minimal redundancy, expansions are observed in gene families related to RNA editing, UV protection and desiccation tolerance. The genome of A. angustus bears the signatures of horizontally transferred genes from bacteria and fungi, in particular of genes operating in stress-response and metabolic pathways. Our study provides insight into the unique features of hornworts and their molecular adaptations to live on land.

Genome-Wide Analysis of the MYB Transcription Factor Superfamily in Physcomitrella patens

MYB transcription factors (TFs) are one of the largest TF families in plants to regulate numerous biological processes. However, our knowledge of the MYB family in Physcomitrella patens is limited. We identified 116 MYB genes in the P. patens genome, which were classified into the R2R3-MYB, R1R2R3-MYB, 4R-MYB, and MYB-related subfamilies. Most R2R3 genes contain 3 exons and 2 introns, whereas R1R2R3 MYB genes contain 10 exons and 9 introns. N3R-MYB (novel 3RMYB) and NR-MYBs (novel RMYBs) with complicated gene structures appear to be novel MYB proteins. In addition, we found that the diversity of the MYB domain was mainly contributed by domain shuffling and gene duplication. RNA-seq analysis suggested that MYBs exhibited differential expression to heat and might play important roles in heat stress responses, whereas CCA1-like MYB genes might confer greater flexibility to the circadian clock. Some R2R3-MYB and CCA1-like MYB genes are preferentially expressed in the archegonium and during the transition from the chloronema to caulonema stage, suggesting their roles in development. Compared with that of algae, the numbers of MYBs have significantly increased, thus our study lays the foundation for further exploring the potential roles of MYBs in the transition from aquatic to terrestrial environments.

Factors Affecting Organelle Genome Stability in Physcomitrella patens

Organelle genomes are essential for plants; however, the mechanisms underlying the maintenance of organelle genomes are incompletely understood. Using the basal land plant Physcomitrella patens as a model, nuclear-encoded homologs of bacterial-type homologous recombination repair (HRR) factors have been shown to play an important role in the maintenance of organelle genome stability by suppressing recombination between short dispersed repeats. In this review, I summarize the factors and pathways involved in the maintenance of genome stability, as well as the repeats that cause genomic instability in organelles in P. patens, and compare them with findings in other plant species. I also discuss the relationship between HRR factors and organelle genome structure from the evolutionary standpoint.

Current Challenges in Studying Alternative Splicing in Plants: The Case of Physcomitrella patens SR Proteins

To colonize different terrestrial habitats, early land plants had to overcome the challenge of coping with harsh new environments. Alternative splicing - an RNA processing mechanism through which splice sites are differentially recognized, originating multiple transcripts and potentially different proteins from the same gene - can be key for plant stress tolerance. Serine/arginine-rich (SR) proteins constitute an evolutionarily conserved family of major alternative splicing regulators that in plants subdivides into six subfamilies. Despite being well studied in animals and a few plant species, such as the model angiosperm Arabidopsis thaliana and the crop Oryza sativa, little is known of these splicing factors in early land plants. Establishing the whole complement of SR proteins in different species is essential to understand the functional and evolutionary significance of alternative splicing. An in silico search for SR proteins in the extant moss Physcomitrella patens revealed inconsistencies both in the published data and available databases, likely arising from automatic annotation lacking adequate manual curation. These misannotations interfere with the description not only of the number and subfamily classification of Physcomitrella SR proteins but also of their domain architecture, potentially hindering the elucidation of their molecular functions. We therefore advise caution when looking into P. patens genomic resources. Our systematic survey nonetheless confidently identified 16 P. patens SR proteins that fall into the six described subfamilies and represent counterparts of well-established members in Arabidopsis and rice. Intensified research efforts should disclose whether SR proteins were already determining alternative splicing modulation and stress tolerance in early land plants.

What Can Long Terminal Repeats Tell Us About the Age of LTR Retrotransposons, Gene Conversion and Ectopic Recombination?

LTR retrotransposons constitute a significant part of plant genomes and their evolutionary dynamics play an important role in genome size changes. Current methods of LTR retrotransposon age estimation are based only on LTR (long terminal repeat) divergence. This has prompted us to analyze sequence similarity of LTRs in 25,144 LTR retrotransposons from fifteen plant species as well as formation of solo LTRs. We found that approximately one fourth of nested retrotransposons showed a higher LTR divergence than the pre-existing retrotransposons into which they had been inserted. Moreover, LTR similarity was correlated with LTR length. We propose that gene conversion can contribute to this phenomenon. Gene conversion prediction in LTRs showed potential converted regions in 25% of LTR pairs. Gene conversion was higher in species with smaller genomes while the proportion of solo LTRs did not change with genome size in analyzed species. The negative correlation between the extent of gene conversion and the abundance of solo LTRs suggests interference between gene conversion and ectopic recombination. Since such phenomena limit the traditional methods of LTR retrotransposon age estimation, we recommend an improved approach based on the exclusion of regions affected by gene conversion.