Search

Lettuce (Lactuca sativa) is a major crop and a member of the large, highly successful Compositae family of flowering plants. Here we present a reference assembly for the species and family. This was generated using whole-genome shotgun Illumina reads plus in vitro proximity ligation data to create large superscaffolds; it was validated genetically and superscaffolds were oriented in genetic bins ordered along nine chromosomal pseudomolecules. We identify several genomic features that may have contributed to the success of the family, including genes encoding Cycloidea-like transcription factors, kinases, enzymes involved in rubber biosynthesis and disease resistance proteins that are expanded in the genome. We characterize 21 novel microRNAs, one of which may trigger phasiRNAs from numerous kinase transcripts. We provide evidence for a whole-genome triplication event specific but basal to the Compositae. We detect 26% of the genome in triplicated regions containing 30% of all genes that are enriched for regulatory sequences and depleted for genes involved in defence.

The genus Phytophthora consists of many notorious pathogens of crops and forestry trees. At present, battling Phytophthora diseases is challenging due to a lack of understanding of their pathogenesis. We investigated the role of small RNAs in regulating soybean defense in response to infection by Phytophthora sojae, the second most destructive pathogen of soybean. Small RNAs, including microRNAs (miRNAs) and small interfering RNAs (siRNAs), are universal regulators that repress target gene expression in eukaryotes. We identified known and novel small RNAs that differentially accumulated during P. sojae infection in soybean roots. Among them, miR393 and miR166 were induced by heat-inactivated P. sojae hyphae, indicating that they may be involved in soybean basal defense. Indeed, knocking down the level of mature miR393 led to enhanced susceptibility of soybean to P. sojae; furthermore, the expression of isoflavonoid biosynthetic genes was drastically reduced in miR393 knockdown roots. These data suggest that miR393 promotes soybean defense against P. sojae. In addition to miRNAs, P. sojae infection also resulted in increased accumulation of phased siRNAs (phasiRNAs) that are predominantly generated from canonical resistance genes encoding nucleotide binding-leucine rich repeat proteins and genes encoding pentatricopeptide repeat-containing proteins. This work identifies specific miRNAs and phasiRNAs that regulate defense-associated genes in soybean during Phytophthora infection.

In plants, heterochromatin is maintained by a small RNA-based gene silencing mechanism known as RNA-directed DNA methylation (RdDM). RdDM requires the non-redundant functions of two plant-specific DNA-dependent RNA polymerases (RNAP), RNAP IV and RNAP V. RNAP IV plays a major role in siRNA biogenesis, while RNAP V may recruit DNA methylation machinery to target endogenous loci for silencing. Although small RNA-generating regions that are dependent on both RNAP IV and RNAP V have been identified previously, the genomic loci targeted by RNAP V for siRNA accumulation and silencing have not been described extensively. To characterize the RNAP V-dependent, heterochromatic siRNA-generating regions in the Arabidopsis genome, we deeply sequenced the small RNA populations of wild-type and RNAP V null mutant (nrpe1) plants. Our results showed that RNAP V-dependent siRNA-generating loci are associated predominately with short repetitive sequences in intergenic regions. Suppression of small RNA production from short repetitive sequences was also prominent in RdDM mutants including dms4, drd1, dms3 and rdm1, reflecting the known association of these RdDM effectors with RNAP V. The genomic regions targeted by RNAP V were small, with an estimated average length of 238 bp. Our results suggest that RNAP V affects siRNA production from genomic loci with features dissimilar to known RNAP IV-dependent loci. RNAP V, along with RNAP IV and DRM1/2, may target and silence a set of small, intergenic transposable elements located in dispersed genomic regions for silencing. Silencing at these loci may be actively reinforced by RdDM.

Most of our understanding of plant genome structure and evolution has come from the careful annotation of small (e.g., 100 kb) sequenced genomic regions or from automated annotation of complete genome sequences. Here, we sequenced and carefully annotated a contiguous 22 Mb region of maize chromosome 4 using an improved pseudomolecule for annotation. The sequence segment was comprehensively ordered, oriented, and confirmed using the maize optical map. Nearly 84% of the sequence is composed of transposable elements (TEs) that are mostly nested within each other, of which most families are low-copy. We identified 544 gene models using multiple levels of evidence, as well as five miRNA genes. Gene fragments, many captured by TEs, are prevalent within this region. Elimination of gene redundancy from a tetraploid maize ancestor that originated a few million years ago is responsible in this region for most disruptions of synteny with sorghum and rice. Consistent with other sub-genomic analyses in maize, small RNA mapping showed that many small RNAs match TEs and that most TEs match small RNAs. These results, performed on approximately 1% of the maize genome, demonstrate the feasibility of refining the B73 RefGen_v1 genome assembly by incorporating optical map, high-resolution genetic map, and comparative genomic data sets. Such improvements, along with those of gene and repeat annotation, will serve to promote future functional genomic and phylogenomic research in maize and other grasses.

BackgroundSeveral high-throughput technologies can measure in parallel the abundance of many mRNA transcripts within a sample. These include the widely-used microarray as well as the more recently developed methods based on sequence tag abundances such as the Massively Parallel Signature Sequencing (MPSS) technology. A comparison of microarray and MPSS technologies can help to establish the metrics for data comparisons across these technology platforms and determine some of the factors affecting the measurement of mRNA abundances using different platforms.ResultsWe compared transcript abundance (gene expression) measurement data obtained using Affymetrix and Agilent microarrays with MPSS data. All three technologies were used to analyze the same set of mRNA samples; these samples were extracted from various wild type Arabidopsis thaliana tissues and floral mutants. We calculated correlations and used clustering methodology to compare the normalized expression data and expression ratios across samples and technologies. Abundance expression measurements were more similar between different samples measured by the same technology than between the same sample measured by different technologies. However, when expression ratios were employed, samples measured by different technologies were found to cluster together more frequently than with abundance expression levels.Furthermore, the two microarray technologies were more consistent with each other than with MPSS. We also investigated probe-position effects on Affymetrix data and tag-position effects in MPSS. We found a similar impact on Affymetrix and MPSS measurements, which suggests that these effects were more likely a characteristic of the RNA sample rather than technology-specific biases.ConclusionComparisons of transcript expression ratios showed greater consistency across platforms than measurements of transcript abundance. In addition, for measurements based on abundances, technology differences can mask the impact of biological differences between samples and tissues.

BackgroundNucleotide binding site-leucine rich repeat (NBS-LRR)-encoding genes comprise the largest class of plant disease resistance genes. The 149 NBS-LRR-encoding genes and the 58 related genes that do not encode LRRs represent approximately 0.8% of all ORFs so far annotated in Arabidopsis ecotype Col-0. Despite their prevalence in the genome and functional importance, there was little information regarding expression of these genes.ResultsWe analyzed the expression patterns of approximately 170 NBS-LRR-encoding and related genes in Arabidopsis Col-0 using multiple analytical approaches: expressed sequenced tag (EST) representation, massively parallel signature sequencing (MPSS), microarray analysis, rapid amplification of cDNA ends (RACE) PCR, and gene trap lines. Most of these genes were expressed at low levels with a variety of tissue specificities. Expression was detected by at least one approach for all but 10 of these genes. The expression of some but not the majority of NBS-LRR-encoding and related genes was affected by salicylic acid (SA) treatment; the response to SA varied among different accessions. An analysis of previously published microarray data indicated that ten NBS-LRR-encoding and related genes exhibited increased expression in wild-type Landsberg erecta (Ler) after flagellin treatment. Several of these ten genes also showed altered expression after SA treatment, consistent with the regulation of R gene expression during defense responses and overlap between the basal defense response and salicylic acid signaling pathways. Enhancer trap analysis indicated that neither jasmonic acid nor benzothiadiazole (BTH), a salicylic acid analog, induced detectable expression of the five NBS-LRR-encoding genes and one TIR-NBS-encoding gene tested; however, BTH did induce detectable expression of the other TIR-NBS-encoding gene analyzed. Evidence for alternative mRNA polyadenylation sites was observed for many of the tested genes. Evidence for alternative splicing was found for at least 12 genes, 11 of which encode TIR-NBS-LRR proteins. There was no obvious correlation between expression pattern, phylogenetic relationship or genomic location of the NBS-LRR-encoding and related genes studied.ConclusionTranscripts of many NBS-LRR-encoding and related genes were defined. Most were present at low levels and exhibited tissue-specific expression patterns. Expression data are consistent with most Arabidopsis NBS-LRR-encoding and related genes functioning in plant defense responses but do not preclude other biological roles.

SUMMARY: Signals originating within plastids modulate organelle differentiation by transcriptionally regulating nuclear‐encoded genes. These retrograde signals are also integral regulators of plant development, including leaf morphology. The clb5 mutant displays severe leaf morphology defects due to Apocarotenoid Signal 1 (ACS1) accumulation in the developmentally arrested plastid. Transcriptomic analysis of clb5 validates that ACS1 accumulation deregulates hundreds of nuclear genes, including the suppression of most genes encoding plastid ribosomal proteins. Herein, we order the molecular events causing the leaf phenotype associated with the accumulation of ACS1, which includes two consecutive retrograde signaling cascades. Firstly, ACS1 originating in the plastid drives inhibition of plastid translation (IPT) via nuclear transcriptome remodeling of chlororibosomal proteins, requiring light as an essential component. Subsequently, IPT results in leaf morphological defects via a GUN1‐dependent pathway shared with seedlings undergoing chemical IPT treatments and is restricted to an early window of the leaf development. Collectively, this work advances our understanding of the complexity within plastid retrograde signaling exemplified by sequential signal exchange and consequences that in a particular temporal and spatial context contribute to the modulation of leaf development. Significance Statement: Accumulation of an undefined apocarotenoid, ACS1, derived from cis‐carotenoids leads to important defects in leaf development via a previously unknown molecular mechanism. In this work we show that accumulation of ACS1 initiates a retrograde signal cascade resulting in the transcriptional modulation of diverse genes and define the molecular events involved in ACS1's regulation of leaf development. Our work illustrates how two sequential plastidial retrograde signaling cascades modulate leaf development during a particular developmental window. These results will be instrumental in efforts to biochemically define ACS1. [ABSTRACT FROM AUTHOR]

Summary: •Plant disease resistance genes (R genes) encode proteins that function to monitor signals indicating pathogenic infection, thus playing a critical role in the plant's defense system. Although many studies have been performed to explore the functional details of these important genes, their origin and evolutionary history remain unclear.•In this study, focusing on the largest group of R genes, the nucleotide‐binding site–leucine‐rich repeat (NBS‐LRR) genes, we conducted an extensive genome‐wide survey of 38 representative model organisms and obtained insights into the evolutionary stage and timing of NBS‐LRR genes.•Our data show that the two major domains, NBS and LRR, existed before the split of prokaryotes and eukaryotes but their fusion was observed only in land plant lineages. The Toll/interleukin‐1 receptor (TIR) class of NBS‐LRR genes probably had an earlier origin than its nonTIR counterpart. The similarities of the innate immune systems of plants and animals are likely to have been shaped by convergent evolution after their independent origins.•Our findings start to unravel the evolutionary history of these important genes from the perspective of comparative genomics and also highlight the important role of reorganizing pre‐existing building blocks in generating evolutionary novelties. [ABSTRACT FROM AUTHOR]

Background: Small non-coding RNAs (21 to 24 nucleotides) regulate a number of developmental processes in plants and animals by silencing genes using multiple mechanisms. Among these, the most conserved classes are microRNAs (miRNAs) and small interfering RNAs (siRNAs), both of which are produced by RNase III-like enzymes called Dicers. Many plant miRNAs play critical roles in nutrient homeostasis, developmental processes, abiotic stress and pathogen responses. Currently, only 70 miRNA have been identified in soybean. Methods: We utilized Illumina's SBS sequencing technology to generate high-quality small RNA (sRNA) data from four soybean (Glycine max) tissues, including root, seed, flower, and nodules, to expand the collection of currently known soybean miRNAs. We developed a bioinformatics pipeline using in-house scripts and publicly available structure prediction tools to differentiate the authentic mature miRNA sequences from other sRNAs and short RNA fragments represented in the public sequencing data. Results: The combined sequencing and bioinformatics analyses identified 129 miRNAs based on hairpin secondary structure features in the predicted precursors. Out of these, 42 miRNAs matched known miRNAs in soybean or other species, while 87 novel miRNAs were identified. We also predicted the putative target genes of all identified miRNAs with computational methods and verified the predicted cleavage sites in vivo for a subset of these targets using the 5' RACE method. Finally, we also studied the relationship between the abundance of miRNA and that of the respective target genes by comparison to Solexa cDNA sequencing data. Conclusion: Our study significantly increased the number of miRNAs known to be expressed in soybean. The bioinformatics analysis provided insight on regulation patterns between the miRNAs and their predicted target genes expression. We also deposited the data in a soybean genome browser based on the UCSC Genome Browser architecture. Using the browser, we annotated the soybean data with miRNA sequences from four tissues and cDNA sequencing data. Overlaying these two datasets in the browser allows researchers to analyze the miRNA expression levels relative to that of the associated target genes. The browser can be accessed at http://digbio. missouri.edu/soybean_mirna/. [ABSTRACT FROM AUTHOR]

Progress in epigenetics has revealed mechanisms that can heritably regulate gene function independent of genetic alterations. Nevertheless, little is known about the role of epigenetics in evolution. This is due in part to scant data on epigenetic variation among natural populations. In plants, small interfering RNA (siRNA) is involved in both the initiation and maintenance of gene silencing by directing DNA methylation and/or histone methylation. Here, we report that, in the model plant Arabidopsis thaliana, a cluster of ∼24 nt siRNAs found at high levels in the ecotype Landsberg erecta (Ler) could direct DNA methylation and heterochromatinization at a hAT element adjacent to the promoter of FLOWERING LOCUS C (FLC), a major repressor of flowering, whereas the same hAT element in ecotype Columbia (Col) with almost identical DNA sequence, generates a set of low abundance siRNAs that do not direct these activities. We have called this hAT element MPF for Methylated region near Promoter of FLC, although de novo methylation triggered by an inverted repeat transgene at this region in Col does not alter its FLC expression. DNA methylation of the Ler allele MPF is dependent on genes in known silencing pathways, and such methylation is transmissible to Col by genetic crosses, although with varying degrees of penetrance. A genome-wide comparison of Ler and Col small RNAs identified at least 68 loci matched by a significant level of ∼24 nt siRNAs present specifically in Ler but not Col, where nearly half of the loci are related to repeat or TE sequences. Methylation analysis revealed that 88% of the examined loci (37 out of 42) were specifically methylated in Ler but not Col, suggesting that small RNA can direct epigenetic differences between two closely related Arabidopsis ecotypes. [ABSTRACT FROM AUTHOR]

Small RNAs (21-24 nt) are involved in gene regulation through translation inhibition, mRNA cleavage, or directing chromatin modifications. In rice, currently ≈240 microRNAs (miRNAs) have been annotated. We sequenced more than four million small RNAs from rice and identified another 24 miRNA genes. Among these, we `found a unique class of miRNAs that derive from natural cis-antisense transcript pairs. This configuration generates miRNAs that can perfectly match their targets. We provide evidence that `the miRNAs function by inducing mRNA cleavage in the middle of their complementary site. Their production requires Dicer-like 1 (DCL1) activity, which is essential for canonical miRNA biogenesis., All of the natural antisense miRNAs (nat-miRNAs) identified in this study have large introns in their precursors that appear critical for nat-miRNA evolution and for the formation of functional mIRNA loci. These findings suggest that other natural cis-antisense loci with similar exon-intron arrangements could be another source of miRNA genes. [ABSTRACT FROM AUTHOR]

Identification of all expressed transcripts in a sequenced genome is essential both for genome analysis and for realization of the goals of systems biology. We used the transcriptional profiling technology called 'massively parallel signature sequencing' to develop a comprehensive expression atlas of rice (Oryza sativa cv Nipponbare). We sequenced 46,971,553 mRNA transcripts from 22 libraries, and 2,953,855 small RNAs from 3 libraries. The data demonstrate widespread transcription throughout the genome, including sense expression of at least 25,500 annotated genes and antisense expression of nearly 9,000 annotated genes. An additional set of ∼15,000 mRNA signatures mapped to unannotated genomic regions. The majority of the small RNA data represented lower abundance short interfering RNAs that match repetitive sequences, intergenic regions and genes. Among these, numerous clusters of highly regulated small RNAs were readily observed. We developed a genome browser (http://mpss.udel.edu/rice) for public access to the transcriptional profiling data for this important crop. [ABSTRACT FROM AUTHOR]

Background: Rice blast, caused by the fungal pathogen Magnaporthe grisea, is a devastating disease causing tremendous yield loss in rice production. The public availability of the complete genome sequence of M. grisea provides ample opportunities to understand the molecular mechanism of its pathogenesis on rice plants at the transcriptome level. To identify all the expressed genes encoded in the fungal genome, we have analyzed the mycelium and appressorium transcriptomes using massively parallel signature sequencing (MPSS), robust-long serial analysis of gene expression (RLSAGE) and oligoarray methods. Results: The MPSS analyses identified 12,531 and 12,927 distinct significant tags from mycelia and appressoria, respectively, while the RL-SAGE analysis identified 16,580 distinct significant tags from the mycelial library. When matching these 12,531 mycelial and 12,927 appressorial significant tags to the annotated CDS, 500 bp upstream and 500 bp downstream of CDS, 6,735 unique genes in mycelia and 7,686 unique genes in appressoria were identified. A total of 7,135 mycelium-specific and 7,531 appressorium-specific significant MPSS tags were identified, which correspond to 2,088 and 1,784 annotated genes, respectively, when matching to the same set of reference sequences. Nearly 85% of the significant MPSS tags from mycelia and appressoria and 65% of the significant tags from the RL-SAGE mycelium library matched to the M. grisea genome. MPSS and RL-SAGE methods supported the expression of more than 9,000 genes, representing over 80% of the predicted genes in M. grisea. About 40% of the MPSS tags and 55% of the RL-SAGE tags represent novel transcripts since they had no matches in the existing M. grisea EST collections. Over 19% of the annotated genes were found to produce both sense and antisense tags in the protein-coding region. The oligoarray analysis identified the expression of 3,793 mycelium-specific and 4,652 appressorium-specific genes. A total of 2,430 mycelial genes and 1,886 appressorial genes were identified by both MPSS and oligoarray. Conclusion: The comprehensive and deep transcriptome analysis by MPSS and RL-SAGE methods identified many novel sense and antisense transcripts in the M. grisea genome at two important growth stages. The differentially expressed transcripts that were identified, especially those specifically expressed in appressoria, represent a genomic resource useful for gaining a better understanding of the molecular basis of M. grisea pathogenicity. Further analysis of the novel antisense transcripts will provide new insights into the regulation and function of these genes in fungal growth, development and pathogenesis in the host plants. [ABSTRACT FROM AUTHOR]

Large-scale sequencing of short mRNA-derived tags can establish the qualitative and quantitative characteristics of a complex transcriptome. We sequenced 12,304,362 tags from five diverse libraries of Arabidopsis thaliana using massively parallel signature sequencing (MPSS). A total of 48,572 distinct signatures, each representing a different transcript, were expressed at significant levels. These signatures were compared to the annotation of the A. thaliana genomic sequence; in the five libraries, this comparison yielded between 17,353 and 18,361 genes with sense expression, and between 5,487 and 8,729 genes with antisense expression. An additional 6,691 MPSS signatures mapped to unannotated regions of the genome. Expression was demonstrated for 1,168 genes for which expression data were previously unknown. Alternative polyadenylation was observed for more than 25% of A. thaliana genes transcribed in these libraries. The MPSS expression data suggest that the A. thaliana transcriptome is complex and contains many as-yet uncharacterized variants of normal coding transcripts. [ABSTRACT FROM AUTHOR]

Provides a correction to the article 'Genome-Wide Analysis of NBS-LLR-encoding Genes in Arabidopsis by Blake C. Meyers, Alexander Kozik, Alyssa Griego, Hanhui Kuang and Richard W. Michelmore, published in the 2003 issue of the journal 'Plant Cell.'

Research on a few model plant–pathogen systems has benefitted from years of tool and resource development. This is not the case for the vast majority of economically and nutritionally important plants, creating a crop improvement bottleneck. Cassava bacterial blight (CBB), caused by Xanthomonas axonopodis pv. manihotis (Xam), is an important disease in all regions where cassava (Manihot esculenta Crantz) is grown. Here, we describe the development of cassava that can be used to visualize one of the initial steps of CBB infection in vivo. Using CRISPR-mediated homology-directed repair (HDR), we generated plants containing scarless insertion of GFP at the 3’ end of CBB susceptibility (S) gene MeSWEET10a. Activation of MeSWEET10a-GFP by the transcription activator-like (TAL) effector TAL20 was subsequently visualized at transcriptional and translational levels. To our knowledge, this is the first such demonstration of HDR via gene editing in cassava. [ABSTRACT FROM AUTHOR]

Polyadenylation sites mark the ends of mRNA transcripts. Alternative polyadenylation (APA) may alter sequence elements and/or the coding capacity of transcripts, a mechanism that has been demonstrated to regulate gene expression and transcriptome diversity. To study the role of APA in transcriptome dynamics, we analyzed a large-scale data set of RNA "tags" that signify poly(A) sites and expression levels of mRNA. These tags were derived from a wide range of tissues and developmental stages that were mutated or exposed to environmental treatments, and generated using digital gene expression (DGE)-based protocols of the massively parallel signature sequencing (MPSS-DGE) and the Illumina sequencing-by-synthesis (SBS-DGE) sequencing platforms. The data offer a global view of APA and how it contributes to transcriptome dynamics. Upon analysis of these data, we found that 60% of Arabidopsis genes have multiple poly(A) sites. Likewise, 47% and 82% of rice genes use APA, supported by MPSS-DGE and SBS-DGE tags, respectively. In both species, 49%-66% of APA events were mapped upstream of annotated stop codons. Interestingly, 10% of the transcriptomes are made up of APA transcripts that are differentially distributed among developmental stages and in tissues responding to environmental stresses, providing an additional level of transcriptome dynamics. Examples of pollen-specific APA switching and salicylic acid treatment-specific APA clearly demonstrated such dynamics. The significance of these APAs is more evident in the 3034 genes that have conserved APA events between rice and Arabidopsis. [ABSTRACT FROM AUTHOR]