Natural variation in Brachypodium distachyon responses to combined abiotic stresses.

The demand for agricultural production is becoming more challenging as climate change increases global temperature and the frequency of extreme weather events. This study examines the phenotypic variation of 149 accessions of Brachypodium distachyon under drought, heat, and the combination of stresses. Heat alone causes the largest amounts of tissue damage while the combination of stresses causes the largest decrease in biomass compared to other treatments. Notably, Bd21-0, the reference line for B. distachyon, did not have robust growth under stress conditions, especially the heat and combined drought and heat treatments. The climate of origin was significantly associated with B. distachyon responses to the assessed stress conditions. Additionally, a GWAS found loci associated with changes in plant height and the amount of damaged tissue under stress. Some of these SNPs were closely located to genes known to be involved in responses to abiotic stresses and point to potential causative loci in plant stress response. However, SNPs found to be significantly associated with a response to heat or drought individually are not also significantly associated with the combination of stresses. This, with the phenotypic data, suggests that the effects of these abiotic stresses are not simply additive, and the responses to the combined stresses differ from drought and heat alone.

Toward a data infrastructure for the Plant Cell Atlas.

We review how a data infrastructure for the Plant Cell Atlas might be built using existing infrastructure and platforms. The Human Cell Atlas has developed an extensive infrastructure for human and mouse single cell data, while the European Bioinformatics Institute has developed a Single Cell Expression Atlas, that currently houses several plant data sets. We discuss issues related to appropriate ontologies for describing a plant single cell experiment. We imagine how such an infrastructure will enable biologists and data scientists to glean new insights into plant biology in the coming decades, as long as such data are made accessible to the community in an open manner.

The plant response to high CO2 levels is heritable and orchestrated by DNA methylation.

Plant responses to abiotic environmental challenges are known to have lasting effects on the plant beyond the initial stress exposure. Some of these lasting effects are transgenerational, affecting the next generation. The plant response to elevated carbon dioxide (CO2 ) levels has been well studied. However, these investigations are typically limited to plants grown for a single generation in a high CO2 environment while transgenerational studies are rare. We aimed to determine transgenerational growth responses in plants after exposure to high CO2 by investigating the direct progeny when returned to baseline CO2 levels. We found that both the flowering plant Arabidopsis thaliana and seedless nonvascular plant Physcomitrium patens continue to display accelerated growth rates in the progeny of plants exposed to high CO2 . We used the model species Arabidopsis to dissect the molecular mechanism and found that DNA methylation pathways are necessary for heritability of this growth response. More specifically, the pathway of RNA-directed DNA methylation is required to initiate methylation and the proteins CMT2 and CMT3 are needed for the transgenerational propagation of this DNA methylation to the progeny plants. Together, these two DNA methylation pathways establish and then maintain a cellular memory to high CO2 exposure.

Specificity of ARGONAUTE7-miR390 interaction and dual functionality in TAS3 trans-acting siRNA formation.

Trans-acting siRNA form through a refined RNAi mechanism in plants. miRNA-guided cleavage triggers entry of precursor transcripts into an RNA-DEPENDENT RNA POLYMERASE6 pathway, and sets the register for phased tasiRNA formation by DICER-LIKE4. Here, we show that miR390-ARGONAUTE7 complexes function in distinct cleavage or noncleavage modes at two target sites in TAS3a transcripts. The AGO7 cleavage, but not the noncleavage, function could be provided by AGO1, the dominant miRNA-associated AGO, but only when AGO1 was guided to a modified target site through an alternate miRNA. AGO7 was highly selective for interaction with miR390, and miR390 in turn was excluded from association with AGO1 due entirely to an incompatible 5' adenosine. Analysis of AGO1, AGO2, and AGO7 revealed a potent 5' nucleotide discrimination function for some, although not all, ARGONAUTEs. miR390 and AGO7, therefore, evolved as a highly specific miRNA guide/effector protein pair to function at two distinct tasiRNA biogenesis steps.

Large structural variations in the haplotype-resolved African cassava genome.

Cassava (Manihot esculenta Crantz, 2n = 36) is a global food security crop. It has a highly heterozygous genome, high genetic load, and genotype-dependent asynchronous flowering. It is typically propagated by stem cuttings and any genetic variation between haplotypes, including large structural variations, is preserved by such clonal propagation. Traditional genome assembly approaches generate a collapsed haplotype representation of the genome. In highly heterozygous plants, this results in artifacts and an oversimplification of heterozygous regions. We used a combination of Pacific Biosciences (PacBio), Illumina, and Hi-C to resolve each haplotype of the genome of a farmer-preferred cassava line, TME7 (Oko-iyawo). PacBio reads were assembled using the FALCON suite. Phase switch errors were corrected using FALCON-Phase and Hi-C read data. The ultralong-range information from Hi-C sequencing was also used for scaffolding. Comparison of the two phases revealed >5000 large haplotype-specific structural variants affecting over 8 Mb, including insertions and deletions spanning thousands of base pairs. The potential of these variants to affect allele-specific expression was further explored. RNA-sequencing data from 11 different tissue types were mapped against the scaffolded haploid assembly and gene expression data are incorporated into our existing easy-to-use web-based interface to facilitate use by the broader plant science community. These two assemblies provide an excellent means to study the effects of heterozygosity, haplotype-specific structural variation, gene hemizygosity, and allele-specific gene expression contributing to important agricultural traits and further our understanding of the genetics and domestication of cassava.

Heating quinoa shoots results in yield loss by inhibiting fruit production and delaying maturity.

Increasing global temperatures and a growing world population create the need to develop crop varieties that provide higher yields in warmer climates. There is growing interest in expanding quinoa cultivation, because of the ability of quinoa to produce nutritious grain in poor soils, with little water and at high salinity. The main limitation to expanding quinoa cultivation, however, is the susceptibility of quinoa to temperatures above approximately 32°C. This study investigates the phenotypes, genes and mechanisms that may affect quinoa seed yield at high temperatures. Using a differential heating system where only roots or only shoots were heated, quinoa yield losses were attributed to shoot heating. Plants with heated shoots lost 60-85% yield as compared with control plants. Yield losses were the result of lower fruit production, which lowered the number of seeds produced per plant. Furthermore, plants with heated shoots had delayed maturity and greater non-reproductive shoot biomass, whereas plants with both heated roots and heated shoots produced higher yields from the panicles that had escaped the heat, compared with the control. This suggests that quinoa uses a type of avoidance strategy to survive heat. Gene expression analysis identified transcription factors differentially expressed in plants with heated shoots and low yield that had been previously associated with flower development and flower opening. Interestingly, in plants with heated shoots, flowers stayed closed during the day while the control flowers were open. Although a closed flower may protect the floral structures, this could also cause yield losses by limiting pollen dispersal, which is necessary to produce fruit in the mostly female flowers of quinoa.

Antiviral ARGONAUTEs Against Turnip Crinkle Virus Revealed by Image-Based Trait Analysis.

RNA-based silencing functions as an important antiviral immunity mechanism in plants. Plant viruses evolved to encode viral suppressors of RNA silencing (VSRs) that interfere with the function of key components in the silencing pathway. As effectors in the RNA silencing pathway, ARGONAUTE (AGO) proteins are targeted by some VSRs, such as that encoded by Turnip crinkle virus (TCV). A VSR-deficient TCV mutant was used to identify AGO proteins with antiviral activities during infection. A quantitative phenotyping protocol using an image-based color trait analysis pipeline on the PlantCV platform, with temporal red, green, and blue imaging and a computational segmentation algorithm, was used to measure plant disease after TCV inoculation. This process captured and analyzed growth and leaf color of Arabidopsis (Arabidopsis thaliana) plants in response to virus infection over time. By combining this quantitative phenotypic data with molecular assays to detect local and systemic virus accumulation, AGO2, AGO3, and AGO7 were shown to play antiviral roles during TCV infection. In leaves, AGO2 and AGO7 functioned as prominent nonadditive, anti-TCV effectors, whereas AGO3 played a minor role. Other AGOs were required to protect inflorescence tissues against TCV. Overall, these results indicate that distinct AGO proteins have specialized, modular roles in antiviral defense across different tissues, and demonstrate the effectiveness of image-based phenotyping to quantify disease progression.

Genome-wide association study (GWAS) of leaf cuticular wax components in Camelina sativa identifies genetic loci related to intracellular wax transport.

BACKGROUND: It is important to explore renewable alternatives (e.g. biofuels) that can produce energy sources to help reduce reliance on fossil oils, and reduce greenhouse gases and waste solids resulted from fossil oils consumption. Camelina sativa is an oilseed crop which has received increasing attention due to its short life cycle, broader adaptation regions, high oil content, high level of omega-3 unsaturated fatty acids, and low-input requirements in agriculture practices. To expand its Camelina production areas into arid regions, there is a need to breed for new drought-tolerant cultivars. Leaf cuticular wax is known to facilitate plant development and growth under water-limited conditions. Dissecting the genetic loci underlying leaf cuticular waxes is important to breed for cultivars with improved drought tolerance. RESULTS: Here we combined phenotypic data and single nucleotide polymorphism (SNP) data from a spring C. sativa diversity panel using genotyping-by-sequencing (GBS) technology, to perform a large-scale genome-wide association study (GWAS) on leaf wax compositions. A total of 42 SNP markers were significantly associated with 15 leaf wax traits including major wax components such as total primary alcohols, total alkanes, and total wax esters as well as their constituents. The vast majority of significant SNPs were associated with long-chain carbon monomers (carbon chain length longer than C28), indicating the important effects of long-chain carbon monomers on leaf total wax biosynthesis. These SNP markers are located on genes directly or indirectly related to wax biosynthesis such as maintaining endoplasmic reticulum (ER) morphology and enabling normal wax secretion from ER to plasma membrane or Golgi network-mediated transport. CONCLUSIONS: These loci could potentially serve as candidates for the genetic control involved in intracellular wax transport that might directly or indirectly facilitate leaf wax accumulation in C. sativa and can be used in future marker-assisted selection (MAS) to breed for the cultivars with high wax content to improve drought tolerance.

Components of Water Use Efficiency Have Unique Genetic Signatures in the Model C4 Grass Setaria.

Plant growth and water use are interrelated processes influenced by genetically controlled morphological and biochemical characteristics. Improving plant water use efficiency (WUE) to sustain growth in different environments is an important breeding objective that can improve crop yields and enhance agricultural sustainability. However, genetic improvement of WUE using traditional methods has proven difficult due to the low throughput and environmental heterogeneity of field settings. To overcome these limitations, this study utilizes a high-throughput phenotyping platform to quantify plant size and water use of an interspecific Setaria italica × Setaria viridis recombinant inbred line population at daily intervals in both well-watered and water-limited conditions. Our findings indicate that measurements of plant size and water use are correlated strongly in this system; therefore, a linear modeling approach was used to partition this relationship into predicted values of plant size given water use and deviations from this relationship at the genotype level. The resulting traits describing plant size, water use, and WUE all were heritable and responsive to soil water availability, allowing for a genetic dissection of the components of plant WUE under different watering treatments. Linkage mapping identified major loci underlying two different pleiotropic components of WUE. This study indicates that alleles controlling WUE derived from both wild and domesticated accessions can be utilized to predictably modulate trait values given a specified precipitation regime in the model C4 genus Setaria.

CG gene body DNA methylation changes and evolution of duplicated genes in cassava.

DNA methylation is important for the regulation of gene expression and the silencing of transposons in plants. Here we present genome-wide methylation patterns at single-base pair resolution for cassava (Manihot esculenta, cultivar TME 7), a crop with a substantial impact in the agriculture of subtropical and tropical regions. On average, DNA methylation levels were higher in all three DNA sequence contexts (CG, CHG, and CHH, where H equals A, T, or C) than those of the most well-studied model plant Arabidopsis thaliana. As in other plants, DNA methylation was found both on transposons and in the transcribed regions (bodies) of many genes. Consistent with these patterns, at least one cassava gene copy of all of the known components of Arabidopsis DNA methylation pathways was identified. Methylation of LTR transposons (GYPSY and COPIA) was found to be unusually high compared with other types of transposons, suggesting that the control of the activity of these two types of transposons may be especially important. Analysis of duplicated gene pairs resulting from whole-genome duplication showed that gene body DNA methylation and gene expression levels have coevolved over short evolutionary time scales, reinforcing the positive relationship between gene body methylation and high levels of gene expression. Duplicated genes with the most divergent gene body methylation and expression patterns were found to have distinct biological functions and may have been under natural or human selection for cassava traits.

P-SAMS: a web site for plant artificial microRNA and synthetic trans-acting small interfering RNA design.

SUMMARY: The Plant Small RNA Maker Site (P-SAMS) is a web tool for the simple and automated design of artificial miRNAs (amiRNAs) and synthetic trans-acting small interfering RNAs (syn-tasiRNAs) for efficient and specific targeted gene silencing in plants. P-SAMS includes two applications, P-SAMS amiRNA Designer and P-SAMS syn-tasiRNA Designer. The navigation through both applications is wizard-assisted, and the job runtime is relatively short. Both applications output the sequence of designed small RNA(s), and the sequence of the two oligonucleotides required for cloning into 'B/c' compatible vectors. AVAILABILITY AND IMPLEMENTATION: The P-SAMS website is available at http://p-sams.carringtonlab.org. CONTACT: acarbonell@ibmcp.upv.es or nfahlgren@danforthcenter.org.

Roles and programming of Arabidopsis ARGONAUTE proteins during Turnip mosaic virus infection.

In eukaryotes, ARGONAUTE proteins (AGOs) associate with microRNAs (miRNAs), short interfering RNAs (siRNAs), and other classes of small RNAs to regulate target RNA or target loci. Viral infection in plants induces a potent and highly specific antiviral RNA silencing response characterized by the formation of virus-derived siRNAs. Arabidopsis thaliana has ten AGO genes of which AGO1, AGO2, and AGO7 have been shown to play roles in antiviral defense. A genetic analysis was used to identify and characterize the roles of AGO proteins in antiviral defense against Turnip mosaic virus (TuMV) in Arabidopsis. AGO1, AGO2 and AGO10 promoted anti-TuMV defense in a modular way in various organs, with AGO2 providing a prominent antiviral role in leaves. AGO5, AGO7 and AGO10 had minor effects in leaves. AGO1 and AGO10 had overlapping antiviral functions in inflorescence tissues after systemic movement of the virus, although the roles of AGO1 and AGO10 accounted for only a minor amount of the overall antiviral activity. By combining AGO protein immunoprecipitation with high-throughput sequencing of associated small RNAs, AGO2, AGO10, and to a lesser extent AGO1 were shown to associate with siRNAs derived from silencing suppressor (HC-Pro)-deficient TuMV-AS9, but not with siRNAs derived from wild-type TuMV. Co-immunoprecipitation and small RNA sequencing revealed that viral siRNAs broadly associated with wild-type HC-Pro during TuMV infection. These results support the hypothesis that suppression of antiviral silencing during TuMV infection, at least in part, occurs through sequestration of virus-derived siRNAs away from antiviral AGO proteins by HC-Pro. These findings indicate that distinct AGO proteins function as antiviral modules, and provide a molecular explanation for the silencing suppressor activity of HC-Pro.

Distinct argonaute-mediated 22G-RNA pathways direct genome surveillance in the C. elegans germline.

Endogenous small RNAs (endo-siRNAs) interact with Argonaute (AGO) proteins to mediate sequence-specific regulation of diverse biological processes. Here, we combine deep-sequencing and genetic approaches to explore the biogenesis and function of endo-siRNAs in C. elegans. We describe conditional alleles of the Dicer-related helicase, drh-3, that abrogate both RNA interference and the biogenesis of endo-siRNAs, called 22G-RNAs. DRH-3 is a core component of RNA-dependent RNA polymerase (RdRP) complexes essential for several distinct 22G-RNA systems. We show that, in the germline, one system is dependent on worm-specific AGOs, including WAGO-1, which localizes to germline nuage structures called P granules. WAGO-1 silences certain genes, transposons, pseudogenes, and cryptic loci. Finally, we demonstrate that components of the nonsense-mediated decay pathway function in at least one WAGO-mediated surveillance pathway. These findings broaden our understanding of the biogenesis and diversity of 22G-RNAs and suggest additional regulatory functions for small RNAs.

Highly specific gene silencing in a monocot species by artificial microRNAs derived from chimeric miRNA precursors.

Artificial microRNAs (amiRNAs) are used for selective gene silencing in plants. However, current methods to produce amiRNA constructs for silencing transcripts in monocot species are not suitable for simple, cost-effective and large-scale synthesis. Here, a series of expression vectors based on Oryza sativa MIR390 (OsMIR390) precursor was developed for high-throughput cloning and high expression of amiRNAs in monocots. Four different amiRNA sequences designed to target specifically endogenous genes and expressed from OsMIR390-based vectors were validated in transgenic Brachypodium distachyon plants. Surprisingly, amiRNAs accumulated to higher levels and were processed more accurately when expressed from chimeric OsMIR390-based precursors that include distal stem-loop sequences from Arabidopsis thaliana MIR390a (AtMIR390a). In all cases, transgenic plants displayed the predicted phenotypes induced by target gene repression, and accumulated high levels of amiRNAs and low levels of the corresponding target transcripts. Genome-wide transcriptome profiling combined with 5'-RLM-RACE analysis in transgenic plants confirmed that amiRNAs were highly specific.

PRG-1 and 21U-RNAs interact to form the piRNA complex required for fertility in C. elegans.

In metazoans, Piwi-related Argonaute proteins have been linked to germline maintenance, and to a class of germline-enriched small RNAs termed piRNAs. Here we show that an abundant class of 21 nucleotide small RNAs (21U-RNAs) are expressed in the C. elegans germline, interact with the C. elegans Piwi family member PRG-1, and depend on PRG-1 activity for their accumulation. The PRG-1 protein is expressed throughout development and localizes to nuage-like structures called P granules. Although 21U-RNA loci share a conserved upstream sequence motif, the mature 21U-RNAs are not conserved and, with few exceptions, fail to exhibit complementarity or evidence for direct regulation of other expressed sequences. Our findings demonstrate that 21U-RNAs are the piRNAs of C. elegans and link this class of small RNAs and their associated Piwi Argonaute to the maintenance of temperature-dependent fertility.

New generation of artificial MicroRNA and synthetic trans-acting small interfering RNA vectors for efficient gene silencing in Arabidopsis.

Artificial microRNAs (amiRNAs) and synthetic trans-acting small interfering RNAs (syn-tasiRNAs) are used for small RNA-based, specific gene silencing or knockdown in plants. Current methods to generate amiRNA or syn-tasiRNA constructs are not well adapted for cost-effective, large-scale production or for multiplexing to specifically suppress multiple targets. Here, we describe simple, fast, and cost-effective methods with high-throughput capability to generate amiRNA and multiplexed syn-tasiRNA constructs for efficient gene silencing in Arabidopsis (Arabidopsis thaliana) and other plant species. amiRNA or syn-tasiRNA inserts resulting from the annealing of two overlapping and partially complementary oligonucleotides are ligated directionally into a zero background BsaI/ccdB-based expression vector. BsaI/ccdB vectors for amiRNA or syn-tasiRNA cloning and expression contain a modified version of Arabidopsis MIR390a or TAS1c precursors, respectively, in which a fragment of the endogenous sequence was substituted by a ccdB cassette flanked by two BsaI sites. Several amiRNA and syn-tasiRNA sequences designed to target one or more endogenous genes were validated in transgenic plants that (1) exhibited the expected phenotypes predicted by loss of target gene function, (2) accumulated high levels of accurately processed amiRNAs or syn-tasiRNAs, and (3) had reduced levels of the corresponding target RNAs.

Functional analysis of three Arabidopsis ARGONAUTES using slicer-defective mutants.

In RNA-directed silencing pathways, ternary complexes result from small RNA-guided ARGONAUTE (AGO) associating with target transcripts. Target transcripts are often silenced through direct cleavage (slicing), destabilization through slicer-independent turnover mechanisms, and translational repression. Here, wild-type and active-site defective forms of several Arabidopsis thaliana AGO proteins involved in posttranscriptional silencing were used to examine several AGO functions, including small RNA binding, interaction with target RNA, slicing or destabilization of target RNA, secondary small interfering RNA formation, and antiviral activity. Complementation analyses in ago mutant plants revealed that the catalytic residues of AGO1, AGO2, and AGO7 are required to restore the defects of Arabidopsis ago1-25, ago2-1, and zip-1 (AGO7-defective) mutants, respectively. AGO2 had slicer activity in transient assays but could not trigger secondary small interfering RNA biogenesis, and catalytically active AGO2 was necessary for local and systemic antiviral activity against Turnip mosaic virus. Slicer-defective AGOs associated with miRNAs and stabilized AGO-miRNA-target RNA ternary complexes in individual target coimmunoprecipitation assays. In genome-wide AGO-miRNA-target RNA coimmunoprecipitation experiments, slicer-defective AGO1-miRNA associated with target RNA more effectively than did wild-type AGO1-miRNA. These data not only reveal functional roles for AGO1, AGO2, and AGO7 slicer activity, but also indicate an approach to capture ternary complexes more efficiently for genome-wide analyses.

Genome sequence and analysis of the Irish potato famine pathogen Phytophthora infestans.

Phytophthora infestans is the most destructive pathogen of potato and a model organism for the oomycetes, a distinct lineage of fungus-like eukaryotes that are related to organisms such as brown algae and diatoms. As the agent of the Irish potato famine in the mid-nineteenth century, P. infestans has had a tremendous effect on human history, resulting in famine and population displacement. To this day, it affects world agriculture by causing the most destructive disease of potato, the fourth largest food crop and a critical alternative to the major cereal crops for feeding the world's population. Current annual worldwide potato crop losses due to late blight are conservatively estimated at $6.7 billion. Management of this devastating pathogen is challenged by its remarkable speed of adaptation to control strategies such as genetically resistant cultivars. Here we report the sequence of the P. infestans genome, which at approximately 240 megabases (Mb) is by far the largest and most complex genome sequenced so far in the chromalveolates. Its expansion results from a proliferation of repetitive DNA accounting for approximately 74% of the genome. Comparison with two other Phytophthora genomes showed rapid turnover and extensive expansion of specific families of secreted disease effector proteins, including many genes that are induced during infection or are predicted to have activities that alter host physiology. These fast-evolving effector genes are localized to highly dynamic and expanded regions of the P. infestans genome. This probably plays a crucial part in the rapid adaptability of the pathogen to host plants and underpins its evolutionary potential.

Evolution and functional diversification of MIRNA genes.

MicroRNAs (miRNAs) are small regulatory RNAs found in diverse eukaryotic lineages. In plants, a minority of annotated MIRNA gene families are conserved between plant families, while the majority are family- or species-specific, suggesting that most known MIRNA genes arose relatively recently in evolutionary time. Given the high proportion of young MIRNA genes in plant species, new MIRNA families are likely spawned and then lost frequently. Unlike highly conserved, ancient miRNAs, young miRNAs are often weakly expressed, processed imprecisely, lack targets, and display patterns of neutral variation, suggesting that young MIRNA loci tend to evolve neutrally. Genome-wide analyses from several plant species have revealed that variation in miRNA foldback expression, structure, processing efficiency, and miRNA size have resulted in the unique functionality of MIRNA loci and resulting miRNAs. Additionally, some miRNAs have evolved specific properties and functions that regulate other transcriptional or posttranscriptional silencing pathways. The evolution of miRNA processing and functional diversity underscores the dynamic nature of miRNA-based regulation in complex regulatory networks.

mut-16 and other mutator class genes modulate 22G and 26G siRNA pathways in Caenorhabditis elegans.

Argonaute-associated siRNAs and Piwi-associated piRNAs have overlapping roles in silencing mobile genetic elements in animals. In Caenorhabditis elegans, mutator (mut) class genes mediate siRNA-guided repression of transposons as well as exogenous RNAi, but their roles in endogenous RNA silencing pathways are not well-understood. To characterize the endogenous small RNAs dependent on mut class genes, small RNA populations from a null allele of mut-16 as well as a regulatory mut-16(mg461) allele that disables only somatic RNAi were subjected to deep sequencing. Additionally, each of the mut class genes was tested for a requirement in 26G siRNA pathways. The results indicate that mut-16 is an essential factor in multiple endogenous germline and somatic siRNA pathways involving several distinct Argonautes and RNA-dependent RNA polymerases. The results also reveal essential roles for mut-2 and mut-7 in the ERGO-1 class 26G siRNA pathway and less critical roles for mut-8, mut-14, and mut-15. We show that transposons are hypersusceptible to mut-16-dependent silencing and identify a requirement for the siRNA machinery in piRNA biogenesis from Tc1 transposons. We also show that the soma-specific mut-16(mg461) mutant allele is present in multiple C. elegans laboratory strains.