Development of a mobile, high-throughput, and low-cost image-based plant growth phenotyping system.

Nondestructive plant phenotyping forms a key technique for unraveling molecular processes underlying plant development and response to the environment. While the emergence of high-throughput phenotyping facilities can further our understanding of plant development and stress responses, their high costs greatly hinder scientific progress. To democratize high-throughput plant phenotyping, we developed sets of low-cost image- and weight-based devices to monitor plant shoot growth and evapotranspiration. We paired these devices to a suite of computational pipelines for integrated and straightforward data analysis. The developed tools were validated for their suitability for large genetic screens by evaluating a cowpea (Vigna unguiculata) diversity panel for responses to drought stress. The observed natural variation was used as an input for a genome-wide association study, from which we identified nine genetic loci that might contribute to cowpea drought resilience during early vegetative development. The homologs of the candidate genes were identified in Arabidopsis (Arabidopsis thaliana) and subsequently evaluated for their involvement in drought stress by using available T-DNA insertion mutant lines. These results demonstrate the varied applicability of this low-cost phenotyping system. In the future, we foresee these setups facilitating the identification of genetic components of growth, plant architecture, and stress tolerance across a wide variety of plant species.

Syntrichia ruralis: emerging model moss genome reveals a conserved and previously unknown regulator of desiccation in flowering plants.

Water scarcity, resulting from climate change, poses a significant threat to ecosystems. Syntrichia ruralis, a dryland desiccation-tolerant moss, provides valuable insights into survival of water-limited conditions. We sequenced the genome of S. ruralis, conducted transcriptomic analyses, and performed comparative genomic and transcriptomic analyses with existing genomes and transcriptomes, including with the close relative S. caninervis. We took a genetic approach to characterize the role of an S. ruralis transcription factor, identified in transcriptomic analyses, in Arabidopsis thaliana. The genome was assembled into 12 chromosomes encompassing 21 169 protein-coding genes. Comparative analysis revealed copy number and transcript abundance differences in known desiccation-associated gene families, and highlighted genome-level variation among species that may reflect adaptation to different habitats. A significant number of abscisic acid (ABA)-responsive genes were found to be negatively regulated by a MYB transcription factor (MYB55) that was upstream of the S. ruralis ortholog of ABA-insensitive 3 (ABI3). We determined that this conserved MYB transcription factor, uncharacterized in Arabidopsis, acts as a negative regulator of an ABA-dependent stress response in Arabidopsis. The new genomic resources from this emerging model moss offer novel insights into how plants regulate their responses to water deprivation.

Composition and function of stress granules and P-bodies in plants.

Stress Granules (SGs) and Processing-bodies (P-bodies) are biomolecular condensates formed in the cell with the highly conserved purpose of maintaining balance between storage, translation, and degradation of mRNA. This balance is particularly important when cells are exposed to different environmental conditions and adjustments have to be made in order for plants to respond to and tolerate stressful conditions. While P-bodies are constitutively present in the cell, SG formation is a stress-induced event. Typically thought of as protein-RNA aggregates, SGs and P-bodies are formed by a process called liquid-liquid phase separation (LLPS), and both their function and composition are very dynamic. Both foci are known to contain proteins involved in translation, protein folding, and ATPase activity, alluding to their roles in regulating mRNA and protein expression levels. From an RNA perspective, SGs and P-bodies primarily consist of mRNAs, though long non-coding RNAs (lncRNAs) have also been observed, and more focus is now being placed on the specific RNAs associated with these aggregates. Recently, metabolites such as nucleotides and amino acids have been reported in purified plant SGs with implications for the energetic dynamics of these condensates. Thus, even though the field of plant SGs and P-bodies is relatively nascent, significant progress has been made in understanding their composition and biological role in stress responses. In this review, we discuss the most recent discoveries centered around SG and P-body function and composition in plants.

Linking discoveries, mechanisms, and technologies to develop a clearer perspective on plant long noncoding RNAs.

Long noncoding RNAs (lncRNAs) are a large and diverse class of genes in eukaryotic genomes that contribute to a variety of regulatory processes. Functionally characterized lncRNAs play critical roles in plants, ranging from regulating flowering to controlling lateral root formation. However, findings from the past decade have revealed that thousands of lncRNAs are present in plant transcriptomes, and characterization has lagged far behind identification. In this setting, distinguishing function from noise is challenging. However, the plant community has been at the forefront of discovery in lncRNA biology, providing many functional and mechanistic insights that have increased our understanding of this gene class. In this review, we examine the key discoveries and insights made in plant lncRNA biology over the past two and a half decades. We describe how discoveries made in the pregenomics era have informed efforts to identify and functionally characterize lncRNAs in the subsequent decades. We provide an overview of the functional archetypes into which characterized plant lncRNAs fit and speculate on new avenues of research that may uncover yet more archetypes. Finally, this review discusses the challenges facing the field and some exciting new molecular and computational approaches that may help inform lncRNA comparative and functional analyses.

Telomerase RNA in Hymenoptera (Insecta) switched to plant/ciliate-like biogenesis.

In contrast to the catalytic subunit of telomerase, its RNA subunit (TR) is highly divergent in size, sequence and biogenesis pathways across eukaryotes. Current views on TR evolution assume a common origin of TRs transcribed with RNA polymerase II in Opisthokonta (the supergroup including Animalia and Fungi) and Trypanosomida on one hand, and TRs transcribed with RNA polymerase III under the control of type 3 promoter, found in TSAR and Archaeplastida supergroups (including e.g. ciliates and Viridiplantae taxa, respectively). Here, we focus on unknown TRs in one of the largest Animalia order - Hymenoptera (Arthropoda) with more than 300 available representative genomes. Using a combination of bioinformatic and experimental approaches, we identify their TRs. In contrast to the presumed type of TRs (H/ACA box snoRNAs transcribed with RNA Polymerase II) corresponding to their phylogenetic position, we find here short TRs of the snRNA type, likely transcribed with RNA polymerase III under the control of the type 3 promoter. The newly described insect TRs thus question the hitherto assumed monophyletic origin of TRs across Animalia and point to an evolutionary switch in TR type and biogenesis that was associated with the divergence of Arthropods.

Evolutionary analysis of the LORELEI gene family in plants reveals regulatory subfunctionalization.

A signaling complex comprising members of the LORELEI (LRE)-LIKE GPI-anchored protein (LLG) and Catharanthus roseus RECEPTOR-LIKE KINASE 1-LIKE (CrRLK1L) families perceive RAPID ALKALINIZATION FACTOR (RALF) peptides and regulate growth, reproduction, immunity, and stress responses in Arabidopsis (Arabidopsis thaliana). Genes encoding these proteins are members of multigene families in most angiosperms and could generate thousands of signaling complex variants. However, the links between expansion of these gene families and the functional diversification of this critical signaling complex as well as the evolutionary factors underlying the maintenance of gene duplicates remain unknown. Here, we investigated LLG gene family evolution by sampling land plant genomes and explored the function and expression of angiosperm LLGs. We found that LLG diversity within major land plant lineages is primarily due to lineage-specific duplication events, and that these duplications occurred both early in the history of these lineages and more recently. Our complementation and expression analyses showed that expression divergence (i.e. regulatory subfunctionalization), rather than functional divergence, explains the retention of LLG paralogs. Interestingly, all but one monocot and all eudicot species examined had an LLG copy with preferential expression in male reproductive tissues, while the other duplicate copies showed highest levels of expression in female or vegetative tissues. The single LLG copy in Amborella trichopoda is expressed vastly higher in male compared to in female reproductive or vegetative tissues. We propose that expression divergence plays an important role in retention of LLG duplicates in angiosperms.

High-Throughput Evolutionary Comparative Analysis of Long Intergenic Noncoding RNAs in Multiple Organisms.

Comparative genomic and transcriptomic analyses can help prioritize and facilitate the functional analysis of long noncoding RNAs (lncRNAs). Evolinc-II is a bioinformatic pipeline that automates comparative analyses, searching for sequence and structural conservation for thousands of lncRNAs at once. In addition, Evolinc-II takes a phylogenetic approach to infer key evolutionary events that may have occurred during the emergence of each query lncRNA. Here, we describe how to use command line or GUI (CyVerse's Discovery Environment) versions of Evolinc-II to identify lncRNA homologs and prioritize them for functional analysis.

A Conserved Long Intergenic Non-coding RNA Containing snoRNA Sequences, lncCOBRA1, Affects Arabidopsis Germination and Development.

Long non-coding RNAs (lncRNAs) are an increasingly studied group of non-protein coding transcripts with a wide variety of molecular functions gaining attention for their roles in numerous biological processes. Nearly 6,000 lncRNAs have been identified in Arabidopsis thaliana but many have yet to be studied. Here, we examine a class of previously uncharacterized lncRNAs termed CONSERVED IN BRASSICA RAPA (lncCOBRA) transcripts that were previously identified for their high level of sequence conservation in the related crop species Brassica rapa, their nuclear-localization and protein-bound nature. In particular, we focus on lncCOBRA1 and demonstrate that its abundance is highly tissue and developmental specific, with particularly high levels early in germination. lncCOBRA1 contains two snoRNAs domains within it, making it the first sno-lincRNA example in a non-mammalian system. However, we find that it is processed differently than its mammalian counterparts. We further show that plants lacking lncCOBRA1 display patterns of delayed germination and are overall smaller than wild-type plants. Lastly, we identify the proteins that interact with lncCOBRA1 and propose a novel mechanism of lincRNA action in which it may act as a scaffold with the RACK1A protein to regulate germination and development, possibly through a role in ribosome biogenesis.

Identification and functional annotation of long intergenic non-coding RNAs in Brassicaceae.

Long intergenic noncoding RNAs (lincRNAs) are a large yet enigmatic class of eukaryotic transcripts that can have critical biological functions. The wealth of RNA-sequencing (RNA-seq) data available for plants provides the opportunity to implement a harmonized identification and annotation effort for lincRNAs that enables cross-species functional and genomic comparisons as well as prioritization of functional candidates. In this study, we processed >24 Tera base pairs of RNA-seq data from >16,000 experiments to identify ∼130,000 lincRNAs in four Brassicaceae: Arabidopsis thaliana, Camelina sativa, Brassica rapa, and Eutrema salsugineum. We used nanopore RNA-seq, transcriptome-wide structural information, peptide data, and epigenomic data to characterize these lincRNAs and identify conserved motifs. We then used comparative genomic and transcriptomic approaches to highlight lincRNAs in our data set with sequence or transcriptional conservation. Finally, we used guilt-by-association analyses to assign putative functions to lincRNAs within our data set. We tested this approach on a subset of lincRNAs associated with germination and seed development, observing germination defects for Arabidopsis lines harboring T-DNA insertions at these loci. LincRNAs with Brassicaceae-conserved putative miRNA binding motifs, small open reading frames, or abiotic-stress modulated expression are a few of the annotations that will guide functional analyses into this cryptic portion of the transcriptome.

2',3'-cAMP treatment mimics the stress molecular response in Arabidopsis thaliana.

The role of the RNA degradation product 2',3'-cyclic adenosine monophosphate (2',3'-cAMP) is poorly understood. Recent studies have identified 2',3'-cAMP in plant material and determined its role in stress signaling. The level of 2',3'-cAMP increases upon wounding, in the dark, and under heat, and 2',3'-cAMP binding to an RNA-binding protein, Rbp47b, promotes stress granule (SG) assembly. To gain further mechanistic insights into the function of 2',3'-cAMP, we used a multi-omics approach by combining transcriptomics, metabolomics, and proteomics to dissect the response of Arabidopsis (Arabidopsis thaliana) to 2',3'-cAMP treatment. We demonstrated that 2',3'-cAMP is metabolized into adenosine, suggesting that the well-known cyclic nucleotide-adenosine pathway of human cells might also exist in plants. Transcriptomics analysis revealed only minor overlap between 2',3'-cAMP- and adenosine-treated plants, suggesting that these molecules act through independent mechanisms. Treatment with 2',3'-cAMP changed the levels of hundreds of transcripts, proteins, and metabolites, many previously associated with plant stress responses, including protein and RNA degradation products, glucosinolates, chaperones, and SG components. Finally, we demonstrated that 2',3'-cAMP treatment influences the movement of processing bodies, confirming the role of 2',3'-cAMP in the formation and motility of membraneless organelles.

Evolution of plant telomerase RNAs: farther to the past, deeper to the roots.

The enormous sequence heterogeneity of telomerase RNA (TR) subunits has thus far complicated their characterization in a wider phylogenetic range. Our recent finding that land plant TRs are, similarly to known ciliate TRs, transcribed by RNA polymerase III and under the control of the type-3 promoter, allowed us to design a novel strategy to characterize TRs in early diverging Viridiplantae taxa, as well as in ciliates and other Diaphoretickes lineages. Starting with the characterization of the upstream sequence element of the type 3 promoter that is conserved in a number of small nuclear RNAs, and the expected minimum TR template region as search features, we identified candidate TRs in selected Diaphoretickes genomes. Homologous TRs were then used to build covariance models to identify TRs in more distant species. Transcripts of the identified TRs were confirmed by transcriptomic data, RT-PCR and Northern hybridization. A templating role for one of our candidates was validated in Physcomitrium patens. Analysis of secondary structure demonstrated a deep conservation of motifs (pseudoknot and template boundary element) observed in all published TRs. These results elucidate the evolution of the earliest eukaryotic TRs, linking the common origin of TRs across Diaphoretickes, and underlying evolutionary transitions in telomere repeats.

Chloroplast quality control pathways are dependent on plastid DNA synthesis and nucleotides provided by cytidine triphosphate synthase two.

Reactive oxygen species (ROS) produced in chloroplasts cause oxidative damage, but also signal to initiate chloroplast quality control pathways, cell death, and gene expression. The Arabidopsis thaliana plastid ferrochelatase two (fc2) mutant produces the ROS singlet oxygen in chloroplasts that activates such signaling pathways, but the mechanisms are largely unknown. Here we characterize one fc2 suppressor mutation and map it to CYTIDINE TRIPHOSPHATE SYNTHASE TWO (CTPS2), which encodes one of five enzymes in Arabidopsis necessary for de novo cytoplasmic CTP (and dCTP) synthesis. The ctps2 mutation reduces chloroplast transcripts and DNA content without similarly affecting mitochondria. Chloroplast nucleic acid content and singlet oxygen signaling are restored by exogenous feeding of the dCTP precursor deoxycytidine, suggesting ctps2 blocks signaling by limiting nucleotides for chloroplast genome maintenance. An investigation of CTPS orthologs in Brassicaceae showed CTPS2 is a member of an ancient lineage distinct from CTPS3. Complementation studies confirmed this analysis; CTPS3 was unable to compensate for CTPS2 function in providing nucleotides for chloroplast DNA and signaling. Our studies link cytoplasmic nucleotide metabolism with chloroplast quality control pathways. Such a connection is achieved by a conserved clade of CTPS enzymes that provide nucleotides for chloroplast function, thereby allowing stress signaling to occur.

Identification and Characterization of the Heat-Induced Plastidial Stress Granules Reveal New Insight Into Arabidopsis Stress Response.

Plants exhibit different physiological and molecular responses to adverse changes in their environment. One such molecular response is the sequestration of proteins, RNAs, and metabolites into cytoplasmic bodies called stress granules (cSGs). Here we report that, in addition to cSGs, heat stress also induces the formation of SG-like foci (cGs) in the chloroplasts of the model plant Arabidopsis thaliana. Similarly to the cSGs, (i) cpSG assemble rapidly in response to stress and disappear when the stress ceases, (ii) cpSG formation is inhibited by treatment with a translation inhibitor (lincomycin), and (iii) cpSG are composed of a stable core and a fluid outer shell. A previously published protocol for cSG extraction was successfully adapted to isolate cpSG, followed by protein, metabolite, and RNA analysis. Analogously to the cSGs, cpSG sequester proteins essential for SG formation, dynamics, and function, also including RNA-binding proteins with prion-like domain, ATPases and chaperones, and the amino acids proline and glutamic acid. However, the most intriguing observation relates to the cpSG localization of proteins, such as a complete magnesium chelatase complex, which is involved in photosynthetic acclimation to stress. These data suggest that cpSG have a role in plant stress tolerance.

Biased Gene Retention in the Face of Introgression Obscures Species Relationships.

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

N6-methyladenosine and RNA secondary structure affect transcript stability and protein abundance during systemic salt stress in Arabidopsis.

After transcription, a messenger RNA (mRNA) is further post-transcriptionally regulated by several features including RNA secondary structure and covalent RNA modifications (specifically N6-methyladenosine, m6A). Both RNA secondary structure and m6A have been demonstrated to regulate mRNA stability and translation and have been independently linked to plant responses to soil salinity levels. However, the effect of m6A on regulating RNA secondary structure and the combinatorial interplay between these two RNA features during salt stress response has yet to be studied. Here, we globally identify RNA-protein interactions and RNA secondary structure during systemic salt stress. This analysis reveals that RNA secondary structure changes significantly during salt stress, and that it is independent of global changes in RNA-protein interactions. Conversely, we find that m6A is anti-correlated with RNA secondary structure in a condition-dependent manner, with salt-specific m6A correlated with a decrease in mRNA secondary structure during salt stress. Taken together, we suggest that salt-specific m6A deposition and the associated loss of RNA secondary structure results in increases in mRNA stability for transcripts encoding abiotic stress response proteins and ultimately increases in protein levels from these stabilized transcripts. In total, our comprehensive analyses reveal important post-transcriptional regulatory mechanisms involved in plant long-term salt stress response and adaptation.

Transcriptomic and evolutionary analysis of the mechanisms by which P. argentatum, a rubber producing perennial, responds to drought.

BACKGROUND: Guayule (Parthenium argentatum Gray) is a drought tolerant, rubber producing perennial shrub native to northern Mexico and the US Southwest. Hevea brasiliensis, currently the world's only source of natural rubber, is grown as a monoculture, leaving it vulnerable to both biotic and abiotic stressors. Isolation of rubber from guayule occurs by mechanical harvesting of the entire plant. It has been reported that environmental conditions leading up to harvest have a profound impact on rubber yield. The link between rubber biosynthesis and drought, a common environmental condition in guayule's native habitat, is currently unclear. RESULTS: We took a transcriptomic and comparative genomic approach to determine how drought impacts rubber biosynthesis in guayule. We compared transcriptional profiles of stem tissue, the location of guayule rubber biosynthesis, collected from field-grown plants subjected to water-deficit (drought) and well-watered (control) conditions. Plants subjected to the imposed drought conditions displayed an increase in production of transcripts associated with defense responses and water homeostasis, and a decrease in transcripts associated with rubber biosynthesis. An evolutionary and comparative analysis of stress-response transcripts suggests that more anciently duplicated transcripts shared among the Asteraceae, rather than recently derived duplicates, are contributing to the drought response observed in guayule. In addition, we identified several deeply conserved long non-coding RNAs (lncRNAs) containing microRNA binding motifs. One lncRNA in particular, with origins at the base of Asteraceae, may be regulating the vegetative to reproductive transition observed in water-stressed guayule by acting as a miRNA sponge for miR166. CONCLUSIONS: These data represent the first genomic analyses of how guayule responds to drought like conditions in agricultural production settings. We identified an inverse relationship between stress-responsive transcripts and those associated with precursor pathways to rubber biosynthesis suggesting a physiological trade-off between maintaining homeostasis and plant productivity. We also identify a number of regulators of abiotic responses, including transcription factors and lncRNAs, that are strong candidates for future projects aimed at modulating rubber biosynthesis under water-limiting conditions common to guayules' native production environment.

Recent emergence and extinction of the protection of telomeres 1c gene in Arabidopsis thaliana.

KEY MESSAGE: Duplicate POT1 genes must rapidly diverge or be inactivated. Protection of telomeres 1 (POT1) encodes a conserved telomere binding protein implicated in both chromosome end protection and telomere length maintenance. Most organisms harbor a single POT1 gene, but in the few lineages where the POT1 family has expanded, the duplicate genes have diversified. Arabidopsis thaliana bears three POT1-like loci, POT1a, POT1b and POT1c. POT1a retains the ancestral function of telomerase regulation, while POT1b is implicated in chromosome end protection. Here we examine the function and evolution of the third POT1 paralog, POT1c. POT1c is a new gene, unique to A. thaliana, and was derived from a duplication event involving the POT1a locus and a neighboring gene encoding ribosomal protein S17. The duplicate S17 locus (dS17) is highly conserved across A. thaliana accessions, while POT1c is highly divergent, harboring multiple deletions within the gene body and two transposable elements within the promoter. The POT1c locus is transcribed at very low to non-detectable levels under standard growth conditions. In addition, no discernable molecular or developmental defects are associated with plants bearing a CRISPR mutation in the POT1c locus. However, forced expression of POT1c leads to decreased telomerase enzyme activity and shortened telomeres. Evolutionary reconstruction indicates that transposons invaded the POT1c promoter soon after the locus was formed, permanently silencing the gene. Altogether, these findings argue that POT1 dosage is critically important for viability and duplicate gene copies are retained only upon functional divergence.

Author Correction: Origin and evolution of the octoploid strawberry genome.

In the version of this article originally published, author Joshua R. Puzey was incorrectly listed as having affiliation 7 (School of Plant Sciences, University of Arizona, Tucson, AZ, USA); affiliation 6 (Department of Biology, College of William and Mary, Williamsburg, VA, USA) is the correct affiliation. The error has been corrected in the HTML and PDF versions of the article.

Origin and evolution of the octoploid strawberry genome.

Cultivated strawberry emerged from the hybridization of two wild octoploid species, both descendants from the merger of four diploid progenitor species into a single nucleus more than 1 million years ago. Here we report a near-complete chromosome-scale assembly for cultivated octoploid strawberry (Fragaria × ananassa) and uncovered the origin and evolutionary processes that shaped this complex allopolyploid. We identified the extant relatives of each diploid progenitor species and provide support for the North American origin of octoploid strawberry. We examined the dynamics among the four subgenomes in octoploid strawberry and uncovered the presence of a single dominant subgenome with significantly greater gene content, gene expression abundance, and biased exchanges between homoeologous chromosomes, as compared with the other subgenomes. Pathway analysis showed that certain metabolomic and disease-resistance traits are largely controlled by the dominant subgenome. These findings and the reference genome should serve as a powerful platform for future evolutionary studies and enable molecular breeding in strawberry.

Tail Wags the Dog? Functional Gene Classes Driving Genome-Wide GC Content in Plasmodium spp.

Plasmodium parasites are valuable models to understand how nucleotide composition affects mutation, diversification, and adaptation. No other observed eukaryotes have undergone such large changes in genomic Guanine-Cytosine (GC) content as seen in the genus Plasmodium (∼30% within 35-40 Myr). Although mutational biases are known to influence GC content in the human-infective Plasmodium vivax and Plasmodium falciparum; no study has addressed how different gene functional classes contribute to genus-wide compositional changes, or if Plasmodium GC content variation is driven by natural selection. Here, we tested the hypothesis that certain gene processes and functions drive variation in global GC content between Plasmodium species. We performed a large-scale comparative genomic analysis using the genomes and predicted genes of 17 Plasmodium species encompassing a wide genomic GC content range. Genic GC content was sorted and divided into ten equally sized quantiles that were then assessed for functional enrichment classes. In agreement that selection on gene classes may drive genomic GC content, trans-membrane proteins were enriched within extreme GC content quantiles (Q1 and Q10). Specifically, variant surface antigens, which primarily interact with vertebrate immune systems, showed skewed GC content distributions compared with other trans-membrane proteins. Although a definitive causation linking GC content, expression, and positive selection within variant surface antigens from Plasmodium vivax, Plasmodium berghei, and Plasmodium falciparum could not be established, we found that regardless of genomic nucleotide composition, genic GC content and expression were positively correlated during trophozoite stages. Overall, these data suggest that, alongside mutational biases, functional protein classes drive Plasmodium GC content change.