Genome-wide analysis of alternative splicing landscapes modulated during plant-virus interactions in Brachypodium distachyon

Nuclear proteome of virus-infected and healthy potato leaves

BACKGROUND

Infection of plants by viruses interferes with expression and subcellular localization of plant proteins. Potyviruses comprise the largest and most economically damaging group of plant-infecting RNA viruses. In virus-infected cells, at least two potyviral proteins localize to nucleus but reasons remain partly unknown.

RESULTS

In this study, we examined changes in the nuclear proteome of leaf cells from a diploid potato line (Solanum tuberosum L.) after infection with potato virus A (PVA; genus Potyvirus; Potyviridae) and compared the data with that acquired for healthy leaves. Gel-free liquid chromatography-coupled to tandem mass spectrometry was used to identify 807 nuclear proteins in the potato line v2-108; of these proteins, 370 were detected in at least two samples of healthy leaves. A total of 313 proteins were common in at least two samples of healthy and PVA-infected leaves; of these proteins, 8 showed differential accumulation. Sixteen proteins were detected exclusively in the samples from PVA-infected leaves, whereas other 16 proteins were unique to healthy leaves. The protein Dnajc14 was only detected in healthy leaves, whereas different ribosomal proteins, ribosome-biogenesis proteins, and RNA splicing-related proteins were over-represented in the nuclei of PVA-infected leaves. Two virus-encoded proteins were identified in the samples of PVA-infected leaves.

CONCLUSIONS

Our results show that PVA infection alters especially ribosomes and splicing-related proteins in the nucleus of potato leaves. The data increase our understanding of potyvirus infection and the role of nucleus in infection. To our knowledge, this is the first study of the nuclear proteome of potato leaves and one of the few studies of changes occurring in nuclear proteomes in response to plant virus infection.

Multi-omics sequencing provides insight into floral transition in Catalpa bungei. C.A. Mey

BACKGROUND

Floral transition plays an important role in development, and proper time is necessary to improve the value of valuable ornamental trees. The molecular mechanisms of floral transition remain unknown in perennial woody plants. "Bairihua" is a type of C. bungei that can undergo floral transition in the first planting year.

RESULTS

Here, we combined short-read next-generation sequencing (NGS) and single-molecule real-time (SMRT) sequencing to provide a more complete view of transcriptome regulation during floral transition in C. bungei. The circadian rhythm-plant pathway may be the critical pathway during floral transition in early flowering (EF) C. bungei, according to horizontal and vertical analysis in EF and normal flowering (NF) C. bungei. SBP and MIKC-MADS-box were seemingly involved in EF during floral transition. A total of 61 hub genes were associated with floral transition in the MEturquoise model with Weighted Gene Co-expression Network Analysis (WGCNA). The results reveal that ten hub genes had a close connection with the GASA homologue gene (Cbu.gene.18280), and the ten co-expressed genes belong to five flowering-related pathways. Furthermore, our study provides new insights into the complexity and regulation of alternative splicing (AS). The ratio or number of isoforms of some floral transition-related genes is different in different periods or in different sub-genomes.

CONCLUSIONS

Our results will be a useful reference for the study of floral transition in other perennial woody plants. Further molecular investigations are needed to verify our sequencing data.

Differential Alternative Splicing Genes and Isoform Regulation Networks of Rapeseed (Brassica napus L.) Infected with Sclerotinia sclerotiorum

Alternative splicing (AS) is a post-transcriptional level of gene expression regulation that increases transcriptome and proteome diversity. How the AS landscape of rapeseed (Brassica napus L.) changes in response to the fungal pathogen Sclerotinia sclerotiorum is unknown. Here, we analyzed 18 RNA-seq libraries of mock-inoculated and S. sclerotiorum-inoculated susceptible and tolerant B. napus plants. We found that infection increased AS, with intron retention being the main AS event. To determine the key genes functioning in the AS response, we performed a differential AS (DAS) analysis. We identified 79 DAS genes, including those encoding splicing factors, defense response proteins, crucial transcription factors and enzymes. We generated coexpression networks based on the splicing isoforms, rather than the genes, to explore the genes’ diverse functions. Using this weighted gene coexpression network analysis alongside a gene ontology enrichment analysis, we identified 11 modules putatively involved in the pathogen defense response. Within these regulatory modules, six DAS genes (ascorbate peroxidase 1, ser/arg-rich protein 34a, unknown function 1138, nitrilase 2, v-atpase f, and amino acid transporter 1) were considered to encode key isoforms involved in the defense response. This study provides insight into the post-transcriptional response of B. napus to S. sclerotiorum infection.

ORF Ι of Mycovirus SsNSRV-1 is Associated with Debilitating Symptoms of Sclerotinia sclerotiorum

We previously identified Sclerotinia sclerotiorum negative-stranded virus 1 (SsNSRV-1), the first (-) ssRNA mycovirus, associated with hypovirulence of its fungal host Sclerotinia sclerotiorum. In this study, functional analysis of Open Reading Frame Ι (ORF Ι) of SsNSRV-1 was performed. The integration and expression of ORF Ι led to defects in hyphal tips, vegetative growth, and virulence of the mutant strains of S. sclerotiorum. Further, differentially expressed genes (DEGs) responding to the expression of ORF Ι were identified by transcriptome analysis. In all, 686 DEGs consisted of 267 up-regulated genes and 419 down-regulated genes. DEGs reprogramed by ORF Ι were relevant to secretory proteins, pathogenicity, transcription, transmembrane transport, protein biosynthesis, modification, and metabolism. Alternative splicing was also detected in all mutant strains, but not in hypovirulent strain AH98, which was co-infected by SsNSRV-1 and Sclerotinia sclerotiorum hypovirus 1 (SsHV-1). Thus, the integrity of SsNSRV-1 genome may be necessary to protect viral mRNA from splicing and inactivation by the host. Taken together, the results suggested that protein ORF Ι could regulate the transcription, translation, and modification of host genes in order to facilitate viral proliferation and reduce the virulence of the host. Therefore, ORF Ι may be a potential gene used for the prevention of S. sclerotiorum.

The Emerging Role of Long Non-Coding RNAs in Plant Defense Against Fungal Stress

Growing interest and recent evidence have identified long non-coding RNA (lncRNA) as the potential regulatory elements for eukaryotes. LncRNAs can activate various transcriptional and post-transcriptional events that impact cellular functions though multiple regulatory functions. Recently, a large number of lncRNAs have also been identified in higher plants, and an understanding of their functional role in plant resistance to infection is just emerging. Here, we focus on their identification in crop plant, and discuss their potential regulatory functions and lncRNA-miRNA-mRNA network in plant pathogen stress responses, referring to possible examples in a model plant. The knowledge gained from a deeper understanding of this colossal special group of plant lncRNAs will help in the biotechnological improvement of crops.

Genome-wide analysis of mammary gland shows modulation of transcriptome landscape with alternative splice variants in Staphylococcus aureus mastitis in mice

Epidemiological mapping shows Staphylococcus aureus to be the leading mastitis causing pathogen in India with diverse genetic lineages circulating in the dairy cattle population. We previously reported that endemic clonal strains of S. aureus isolated from subclinical mastitis lead to specific alteration of epigenetic modulators resulting in deviating immune response in intramammary infection mouse model. However, the extent of transcriptome modulation and associated alternative splicing in S. aureus mastitis is poorly understood. Hence, to gain a deeper insight of the extent of modulation of transcriptome landscape, we expanded the study here using high throughput, paired-end RNA sequencing analysis of the mouse mammary gland inoculated with three strains of S. aureus (SA1, SA2, and SA3) possessing specific genotype, virulence and enterotoxin traits. Overall, we detected 35,878 transcripts in S. aureus inoculated mammary gland, 23% more than those annotated in the reference genome. Expression of 20,756 transcripts was > 1 fragment per kilobase of transcript per million mapped fragments and 25.95% of multi-exonic genes were alternatively spliced. We noted Alternative Splicing (AS) events for > 100 immune-related genes. S. aureus infection quantitatively altered AS events in mice mammary gland. Collectively, the majority of differentially expressed significant genes clustered into immune-associated, cell adhesion and metabolic process categories. We observed AS events for 379 transcripts of genes putatively encoding several splicing associated proteins and transcription factors besides inflammatory mediators. The present analysis provides new insights into global transcriptome landscape and AS events in host-defense related genes in response to S. aureus intramammary infection, suggesting the need for studies focusing on multi-target and/or network therapeutics approach to combat mastitis.

GenFam: A web application and database for gene family-based classification and functional enrichment analysis

Genome-scale studies using high-throughput sequencing (HTS) technologies generate substantial lists of differentially expressed genes under different experimental conditions. These gene lists need to be further mined to narrow down biologically relevant genes and associated functions in order to guide downstream functional genetic analyses. A popular approach is to determine statistically overrepresented genes in a user-defined list through enrichment analysis tools, which rely on functional annotations of genes based on Gene Ontology (GO) terms. Here, we propose a new computational approach, GenFam, which allows annotation, classification, and enrichment of genes based on their gene family, thus simplifying identification of candidate gene families and associated genes that may be relevant to the query. GenFam and its integrated database comprises of three hundred and eighty-four unique gene families and supports gene family analyses for sixty plant genomes. Four comparative case studies with plant species belonging to different clades and families were performed using GenFam which demonstrated its robustness and comprehensiveness over preexisting functional enrichment tools. To make it readily accessible for plant biologists, GenFam is available as a web-based application where users can input gene IDs and export enrichment results in both tabular and graphical formats. Users can also customize analysis parameters by choosing from the various statistical enrichment tests and multiple testing correction methods. Additionally, the web-based application, source code, and database are freely available to use and download. Website: http://mandadilab.webfactional.com/home/. Source code and database: http://mandadilab.webfactional.com/home/dload/.

Alternative splicing is a Sorghum bicolor defense response to fungal infection

This study provides new insights that alternative splicing participates with transcriptional control in defense responses to Colletotrichum sublineola in sorghum In eukaryotic organisms, alternative splicing (AS) is an important post-transcriptional mechanism to generate multiple transcript isoforms from a single gene. Protein variants translated from splicing isoforms may have altered molecular characteristics in signal transduction and metabolic activities. However, which transcript isoforms will be translated into proteins and the biological functions of the resulting proteoforms are yet to be identified. Sorghum is one of the five major cereal crops, but its production is severely affected by fungal diseases. For example, sorghum anthracnose caused by Colletotrichum sublineola greatly reduces grain yield and biomass production. In this study, next-generation sequencing technology was used to analyze C. sublineola-inoculated sorghum seedlings compared with mock-inoculated control. It was identified that AS regulation may be as important as traditional transcriptional control during defense responses to fungal infection. Moreover, several genes involved in flavonoid and phenylpropanoid biosynthetic pathways were found to undergo multiple AS modifications. Further analysis demonstrated that non-conventional targets of both 5'- and 3'-splice sites were alternatively used in response to C. sublineola infection. Splicing factors were also affected at both transcriptional and post-transcriptional levels. As the first transcriptome report on C. sublineola infected sorghum, our work also suggested that AS plays crucial functions in defense responses to fungal invasion.

Transcriptome-wide alternative splicing modulation during plant-pathogen interactions in wheat

Alternative splicing (AS) enhances the diversities of both transcripts and proteins in eukaryotes, which contribute to stress adaptation. To catalog wheat (Triticum aestivum L.) AS genes, we characterized 45 RNA-seq libraries from wheat seedlings infected by powdery mildew, Blumeria graminis f. sp. tritici (Bgt) or stripe rust fungus, Puccinia striiformis f. sp. tritici (Pst). We discovered that 11.2% and 10.4% of the multiexon genes had AS transcripts during Bgt and Pst infections, respectively. In response to fungal infection, wheat modulated AS not only in disease resistance proteins, but also in splicing related factors. Apart from the stress induced or activated splicing variants by pathogen, the differential expression profiles were fold increased through changing the ratio of full spliced transcripts versus intron retention (IR) transcripts. Comparing AS transcripts produced by the same gene in Bgt with Pst stress, the spliced terminal exons and the stranded introns are independent and different. This demonstrated that differential induction of specific splice variants were activated against two fungal pathogens. The specific induced AS genes in the Pst-resistant plants were enriched in improving the membrane permeability and protein modification ability, whereas gene expression involved in protein translation and transport were strengthened in Pst-susceptible plants.

siRNAs Derived from Cymbidium Mosaic Virus and Odontoglossum Ringspot Virus Down-modulated the Expression Levels of Endogenous Genes in Phalaenopsis equestris

Interplay between Cymbidium mosaic virus (CymMV)/Odontoglossum ringspot virus (ORSV) and its host plant Phalaenopsis equestris remain largely unknown, which led to deficiency of effective measures to control disease of P. equestris caused by infecting viruses. In this study, for the first time, we characterized viral small interfering RNAs (vsiRNAs) profiles in P. equestris co-infected with CymMV and ORSV through small RNA sequencing technology. CymMV and ORSV small interfering RNAs (siRNAs) demonstrated several general and specific/new characteristics. vsiRNAs, with A/U bias at the first nucleotide, were predominantly 21-nt long and they were derived predominantly (90%) from viral positive-strand RNA. 21-nt siRNA duplexes with 0-nt overhangs were the most abundant 21-nt duplexes, followed by 2-nt overhangs and then 1-nt overhangs 21-nt duplexes in infected P. equestris. Continuous but heterogeneous distribution and secondary structures prediction implied that vsiRNAs originate predominantly by direct Dicer-like enzymes cleavage of imperfect duplexes in the most folded regions of the positive strand of both viruses RNA molecular. Furthermore, we totally predicted 54 target genes by vsiRNAs with psRNATarget server, including disease/stress response-related genes, RNA interference core components, cytoskeleton-related genes, photosynthesis or energy supply related genes. Gene Ontology classification showed that a majority of the predicted targets were related to cellular components and cellular processes and performed a certain function. All target genes were down-regulated with different degree by vsiRNAs as shown by real-time reverse transcription polymerase chain reaction. Taken together, CymMV and ORSV siRNAs played important roles in interplay with P. equestris by down modulating the expression levels of endogenous genes in host plant.

Alternative Splicing in the Regulation of Plant-Microbe Interactions

As sessile organisms, plants are continuously exposed to a wide range of biotic interactions. While some biotic interactions are beneficial or even essential for the plant (e.g. rhizobia and mycorrhiza), others such as pathogens are detrimental and require fast adaptation. Plants partially achieve this growth and developmental plasticity by modulating the repertoire of genes they express. In the past few years, high-throughput transcriptome sequencing have revealed that, in addition to transcriptional control of gene expression, post-transcriptional processes, notably alternative splicing (AS), emerged as a key mechanism for gene regulation during plant adaptation to the environment. AS not only can increase proteome diversity by generating multiple transcripts from a single gene but also can reduce gene expression by yielding isoforms degraded by mechanisms such as nonsense-mediated mRNA decay. In this review, we will summarize recent discoveries detailing the contribution of AS to the regulation of plant-microbe interactions, with an emphasis on the modulation of immunity receptor function and other components of the signaling pathways that deal with pathogen responses. We will also discuss emerging evidences that AS could contribute to dynamic reprogramming of the plant transcriptome during beneficial interactions, such as the legume-symbiotic interaction.

Panicum Mosaic Virus and Its Satellites Acquire RNA Modifications Associated with Host-Mediated Antiviral Degradation

Positive-sense RNA viruses in the Tombusviridae family have genomes lacking a 5' cap structure and prototypical 3' polyadenylation sequence. Instead, these viruses utilize an extensive network of intramolecular RNA-RNA interactions to direct viral replication and gene expression. Here we demonstrate that the genomic RNAs of Panicum mosaic virus (PMV) and its satellites undergo sequence modifications at their 3' ends upon infection of host cells. Changes to the viral and subviral genomes arise de novo within Brachypodium distachyon (herein called Brachypodium) and proso millet, two alternative hosts of PMV, and exist in the infections of a native host, St. Augustinegrass. These modifications are defined by polyadenylation [poly(A)] events and significant truncations of the helper virus 3' untranslated region-a region containing satellite RNA recombination motifs and conserved viral translational enhancer elements. The genomes of PMV and its satellite virus (SPMV) were reconstructed from multiple poly(A)-selected Brachypodium transcriptome data sets. Moreover, the polyadenylated forms of PMV and SPMV RNAs copurify with their respective mature icosahedral virions. The changes to viral and subviral genomes upon infection are discussed in the context of a previously understudied poly(A)-mediated antiviral RNA degradation pathway and the potential impact on virus evolution.IMPORTANCE The genomes of positive-sense RNA viruses have an intrinsic capacity to serve directly as mRNAs upon viral entry into a host cell. These RNAs often lack a 5' cap structure and 3' polyadenylation sequence, requiring unconventional strategies for cap-independent translation and subversion of the cellular RNA degradation machinery. For tombusviruses, critical translational regulatory elements are encoded within the 3' untranslated region of the viral genomes. Here we describe RNA modifications occurring within the genomes of Panicum mosaic virus (PMV), a prototypical tombusvirus, and its satellite agents (i.e., satellite virus and noncoding satellite RNAs), all of which depend on the PMV-encoded RNA polymerase for replication. The atypical RNAs are defined by terminal polyadenylation and truncation events within the 3' untranslated region of the PMV genome. These modifications are reminiscent of host-mediated RNA degradation strategies and likely represent a previously underappreciated defense mechanism against invasive nucleic acids.

Genome-wide analysis of alternative splicing divergences between Brassica hexaploid and its parents

Compared with its parents, Brassica hexaploid underwent significant AS changes, which may provide diversified gene expression regulation patterns and could enhance its adaptability during evolution Polyploidization is considered a significant evolution force that promotes species formation. Alternative splicing (AS) plays a crucial role in multiple biological processes during plant growth and development. To explore the effects of allopolyploidization on the AS patterns of genes, a genome-wide AS analysis was performed by RNA-seq in Brassica hexaploid and its parents. In total, we found 7913 (27540 AS events), 14447 (70179 AS events), and 13205 (60804 AS events) AS genes in Brassica rapa, Brassica carinata, and Brassica hexaploid, respectively. A total of 920 new AS genes were discovered in Brassica hexaploid. There were 56 differently spliced genes between Brassica hexaploid and its parents. In addition, most of the alternative 5' splice sites were located 4 bp upstream of the dominant 5' splice sites, and most of the alternative 3' splice sites were located 3 bp downstream of the dominant 3' splice sites in Brassica hexapliod, which was similar to B. carinata. Furthermore, we cloned and sequenced all amplicons from the RT-PCR products of GRP7/8, namely, Bol045859, Bol016025 and Bol02880. The three genes were found to produce AS transcripts in a new way. The AS patterns of genes were diverse between Brassica hexaploid and its parents, including the loss and gain of AS events. Allopolyploidization changed alternative splicing sites of pre-mRNAs in Brassica hexaploid, which brought about alterations in the sequences of transcripts. Our study provided novel insights into the AS patterns of genes in allopolyploid plants, which may provide a reference for the study of polyploidy adaptability.

Identification, expression, alternative splicing and functional analysis of pepper WRKY gene family in response to biotic and abiotic stresses

WRKY proteins are a large group of plant transcription factors that are involved in various biological processes, including biotic and abiotic stress responses, hormone response, plant development, and metabolism. WRKY proteins have been identified in several plants, but only a few have been identified in Capsicum annuum. Here, we identified a total of 62 WRKY genes in the latest pepper genome. These genes were classified into three groups (Groups 1–3) based on the structural features of their proteins. The structures of the encoded proteins, evolution, and expression under normal growth conditions were analyzed and 35 putative miRNA target sites were predicted in 20 CaWRKY genes. Moreover, the response to cold or CMV treatments of selected WRKY genes were examined to validate the roles under stresses. And alternative splicing (AS) events of some CaWRKYs were also identified under CMV infection. Promoter analysis confirmed that CaWRKY genes are involved in growth, development, and biotic or abiotic stress responses in hot pepper. The comprehensive analysis provides fundamental information for better understanding of the signaling pathways involved in the WRKY-mediated regulation of developmental processes, as well as biotic and abiotic stress responses.

The Arabidopsis splicing regulator SR45 confers salt tolerance in a splice isoform-dependent manner

Functions of most splice isoforms that are generated by alternative splicing are unknown. We show that two splice variants that encode proteins differing in only eight amino acids have distinct functions in a stress response. Serine/arginine-rich (SR) and SR-like proteins, a conserved family of RNA binding proteins across eukaryotes, play important roles in pre-mRNA splicing and other post-transcriptional processes. Pre-mRNAs of SR and SR-like proteins undergo extensive alternative splicing in response to diverse stresses and produce multiple splice isoforms. However, the functions of most splice isoforms remain elusive. Alternative splicing of pre-mRNA of Arabidopsis SR45, which encodes an SR-like splicing regulator, generates two isoforms (long-SR45.1 and short-SR45.2). The proteins encoded by these two isoforms differ in eight amino acids. Here, we investigated the role of SR45 and its splice variants in salt stress tolerance. The loss of SR45 resulted in enhanced sensitivity to salt stress and changes in expression and splicing of genes involved in regulating salt stress response. Interestingly, only the long isoform (SR45.1) rescued the salt-sensitive phenotype as well as the altered gene expression and splicing patterns in the mutant. These results suggest that SR45 positively regulates salt tolerance. Furthermore, only the long isoform is required for SR45-mediated salt tolerance.

Genome-Wide Analysis of Serine/Arginine-Rich Protein Family in Wheat and Brachypodium distachyon

By regulating the pre-mRNA splicing of other genes and themselves, plant serine/arginine-rich (SR) proteins play important roles in development and in response to abiotic stresses. Presently, the functions of most plant SR protein genes remain unclear. Wheat (Triticumaestivum) and Brachypodiumdistachyon are closely related species. In this study, 40 TaSR and 18 BdSR proteins were identified respectively, and they were classified into seven subfamilies: SR, RS, SCL, RSZ, RS2Z, SC35, and SR45. Similar to Arabidopsis and rice SR protein genes, most TaSR and BdSR protein genes are expressed extensively. Surprisingly, real-time polymerase chain reaction (RT-PCR) analyses showed that no alternative splicing event was found in TaSR protein genes, and only six BdSR protein genes are alternatively spliced genes. The detected alternatively spliced BdSR protein genes and transcripts are much fewer than in Arabidopsis, rice, maize, and sorghum. In the promoter regions, 92 development-related, stress-related, and hormone-related cis-elements were detected, indicating their functions in development and in response to environmental stresses. Meanwhile, 19 TaSR and 16 BdSR proteins were predicted to interact with other SR proteins or non-SR proteins, implying that they are involved in other functions in addition to modulating pre-mRNA splicing as essential components of the spliceosome. These results lay a foundation for further analyses of these genes.

Genome-wide alternative splicing landscapes modulated by biotrophic sugarcane smut pathogen

Alternative splicing (AS) promotes transcriptome and proteome diversity during growth, development, and stress responses in eukaryotes. Genome-wide studies of AS in sugarcane (Saccharum spp.) are lacking, mainly due to the absence of a high-quality sequenced reference genome, sugarcane's large, complex genome, and the variable chromosome numbers and polyploidy of sugarcane cultivars. Here, we analyzed changes in the sugarcane isoform-level transcriptome and AS landscape during infection with the smut fungus (Sporisorium scitamineum) using a hybrid approach involving Sorghum bicolor reference-based and Trinity de novo mapping tools. In total, this analysis detected 16,039 and 15,379 transcripts (≥2 FPKM) at 5 and 200 days after infection, respectively. A conservative estimate of isoform-level expression suggested that approximately 5,000 (14%) sugarcane genes undergo AS. Differential expression analysis of the alternatively spliced genes in healthy and smut-infected sugarcane revealed 896 AS events modulated at different stages of infection. Gene family and gene ontology functional enrichment analysis of the differentially spliced genes revealed overrepresentation of functional categories related to the cell wall, defense, and redox homeostasis pathways. Our study provides novel insight into the AS landscape of sugarcane during smut disease interactions.

Identification and genetic analysis of alternative splicing of long non-coding RNAs in tomato initial flowering stage

Alternative splicing (AS) is a key modulator of development in many eukaryotic organisms. Long non-coding RNAs (lncRNAs) are a class of non-coding RNAs that play essential regulatory roles in various developmental processes and stress responses. However, the functions of AS lncRNAs during the initial flowering of tomato are largely unknown. This study was designed to investigate the AS pattern of lncRNAs in tomato flower, leaf, and root tissues at the initial flowering stage. Using RNA-Seq, we found that 72.55% of lncRNAs underwent AS in these tissues, yielding a total of 16,995 AS events. Among them, the main type of AS event is alternative first exon (AFE), followed by retained intron (RI). We performed candidate target genes analysis on tissue-specific AS lncRNA, and the results indicated that the candidate target genes of these lncRNAs may be involved in the regulation of circadian rhythm, plant immunity, cellulose synthesis and phosphate-containing compound metabolic process. Moreover, a total of 73,085 putative SNPs and 15,679 InDels were detected, and the potential relationship between the AS of lncRNAs and interesting SNP and InDel loci, as well as their numbers, revealed their effects on tomato genetic diversity and genomic stability. Our data provide new insights into the complexity of the transcriptome and the regulation of AS.

Comparative transcriptomics analysis uncovers alternative splicing events and molecular markers in cabbage (Brassica oleracea L.)

Alternative splicing (AS) events were identified and verified in cabbage by comparative transcriptome analysis. The corresponding markers were developed and the germplasm resources were identified. Alternative splicing (AS) is a central regulatory mechanism that greatly contributes to plant gene expression and transcriptome diversity. A large body of evidence has shown that AS complexity is relevant for plant development, evolution, complexity, and adaptation. Both insertion/deletion (InDel) and single nucleotide polymorphism (SNP) are typically co-dominant inheritance markers and have abundant polymorphisms. These have been widely used for marker-assisted selection, genetic mapping, and germplasm identification in plants. However, little is known about the molecular mechanisms underlying AS events and the development of markers including SNP and InDel from the cabbage transcriptome. In this study, three cabbage transcriptome datasets were collected and aligned to the cabbage reference genome to analyze AS events and marker development. 31,524 AS events were identified from three cabbage genotypes, accounting for 20.8% of the total cabbage genes. Alternative 3' splice site donor (A3SS) was the most frequent type of the four main AS events in cabbage. 70,475 InDels and 706,269 SNPs were identified with average frequencies of 1 InDel/6.9 kb and 1 SNP/0.7 kb, respectively. 71,942 potential SSRs were identified in 53,129 assembled unigenes with a density of 1 SSR/6.8 kb. The ratio of SNPs with synonymous/non-synonymous mutations was 1:0.65. 142 InDels and 36 SNPs were randomly selected and validated via Sanger sequencing and polymorphism was found among 66.2% of the InDels and 78.6% of the SNPs. Furthermore, 35 informative InDel markers were successfully used for genetic diversity analysis on 36 cabbage accessions. These results facilitate understanding of the molecular regulation mechanism underlying AS events in cabbage. They also provide molecular marker resource data for genetic mapping construction and germplasm identification, and facilitate the genetic improvement of cabbage via breeding.

Alternative splicing coupled to nonsense-mediated mRNA decay contributes to the high-altitude adaptation of maca (Lepidium meyenii)

Alpine plants remain the least studied plant communities in terrestrial ecosystems. However, how they adapt to high-altitude environments is far from clear. Here, we used RNA-seq to investigate a typical alpine plant maca (Lepidium meyenii) to understand its high-altitude adaptation at transcriptional and post-transcriptional level. At transcriptional level, we found that maca root significantly up-regulated plant immunity genes in day-time comparing to night-time, and up-regulated abiotic (cold/osmotic) stress response genes in Nov and Dec comparing to Oct. In addition, 17 positively selected genes were identified, which could be involved in mitochondrion. At post-transcriptional level, we found that maca had species-specific characterized alternative splicing (AS) profile which could be influenced by stress environments. For example, the alternative 3' splice site events (A3SS, 39.62%) were predominate AS events in maca, rather than intron retention (IR, 23.17%). Interestingly, besides serine/arginine-rich (SR) proteins and long non-coding RNAs (lncRNAs), a lot of components in nonsense-mediated mRNA decay (NMD) were identified under differential alternative splicing (DAS), supporting AS coupled to NMD as essential mechanisms for maca's stress responses and high-altitude adaptation. Taken together, we first attempted to unveil maca's high-altitude adaptation mechanisms based on transcriptome and post-transcriptome evidence. Our data provided valuable insights to understand the high-altitude adaptation of alpine plants.

Viral modulation of cellular RNA alternative splicing: A new key player in virus-host interactions?

Upon viral infection, a tug of war is triggered between host cells and viruses to maintain/gain control of vital cellular functions, the result of which will ultimately dictate the fate of the host cell. Among these essential cellular functions, alternative splicing (AS) is an important RNA maturation step that allows exons, or parts of exons, and introns to be retained in mature transcripts, thereby expanding proteome diversity and function. AS is widespread in higher eukaryotes, as it is estimated that nearly all genes in humans are alternatively spliced. Recent evidence has shown that upon infection by numerous viruses, the AS landscape of host‐cells is affected. In this review, we summarize recent advances in our understanding of how virus infection impacts the AS of cellular transcripts. We also present various molecular mechanisms allowing viruses to modulate cellular AS. Finally, the functional consequences of these changes in the RNA splicing signatures during virus–host interactions are discussed.

This article is categorized under:

RNA in Disease and Development > RNA in Disease

RNA Processing > Splicing Regulation/Alternative Splicing

Potential use of the Pteris vittata arsenic hyperaccumulation-regulation network for phytoremediation

Arsenic accumulation in soil is a global problem typically addressed using phytoremediation methods. Pteris vittata, a model arsenic hyperaccumulator, has great potential as a genetically engineered plant for phytoremediation. However, the lack of omic information on this species has severely limited the identification and application of its arsenic hyperaccumulation and regulation components. In this study, we used an optimized single-molecular real-time (SMRT) strategy to create a de novo full-length transcriptomic-tonoplast proteomic database for this unsequenced fern and to determine the genetic components underlying its arsenic hyperaccumulation-regulation mechanisms. We established a comprehensive network consisting of six major transporter families, two novel resistance pathways, and a regulatory system by examining alternative splicing (AS) and long non-coding RNA (lncRNA) in different tissues following As(III) and As(V) treatment. The database and network established in this study will deepen our understanding of the unique hyperaccumulation and regulation mechanisms of P. vittata, ultimately providing a valuable resource for futher research on phytoremediation of arsenic-contaminated soil.

Brassica napus Infected with Leptosphaeria maculans

Alternative splicing (AS) is a post-transcriptional regulatory process that enhances transcriptome diversity, thereby affecting plant growth, development, and stress responses. To identify the new transcripts and changes in the isoform-level AS landscape of rapeseed (Brassica napus) infected with the fungal pathogen Leptosphaeria maculans, we compared eight RNA-seq libraries prepared from mock-inoculated and inoculated B. napus cotyledons and stems. The AS events that occurred in stems were almost the same as those in cotyledons, with intron retention representing the most common AS pattern. We identified 1892 differentially spliced genes between inoculated and uninoculated plants. We performed a weighted gene co-expression network analysis (WGCNA) to identify eight co-expression modules and their Hub genes, which are the genes most connected with other genes within each module. There are nine Hub genes, encoding nine transcription factors, which represent key regulators of each module, including members of the NAC, WRKY, TRAF, AP2/ERF-ERF, C2H2, C2C2-GATA, HMG, bHLH, and C2C2-CO-like families. Finally, 52 and 117 alternatively spliced genes in cotyledons and stems were also differentially expressed between mock-infected and infected materials, such as HMG and C2C2-Dof; which have dual regulatory mechanisms in response to L. maculans. The splicing of the candidate genes identified in this study could be exploited to improve resistance to L. maculans.

Spaceflight-induced alternative splicing during seedling development in Arabidopsis thaliana

Plants grown in spaceflight experience novel environmental signals, including those associated with microgravity and ionizing radiation. Spaceflight triggers a response involving transcriptional re-programming and altered cell morphology, though many aspects of this response remain uncharacterized. We analyzed the spaceflight-induced transcriptome with a focus on genes that undergo alternative splicing to examine differential splicing associated with spaceflight—an unstudied characteristic of the molecular response to spaceflight exposure. RNA sequence data obtained during the APEX03 spaceflight experiment that was collected from two Arabidopsis thaliana ecotypes at two seedling stages grown onboard the International Space Station, or as ground controls at Kennedy Space Center, were re-examined to detect alternative splicing differences induced by spaceflight. Presence/absence variation analysis was used to identify putative expression-level differences in alternatively spliced isoforms between spaceflight and ground controls and was followed by analysis of significant differential alternative splicing. This study provides the first evidence of a role for alternative splicing in the molecular processes of physiological adaptation to the spaceflight environment.

A Two-Amino Acid Difference in the Coat Protein of Satellite panicum mosaic virus Isolates Is Responsible for Differential Synergistic Interactions with Panicum mosaic virus

Panicum mosaic virus (PMV) (genus Panicovirus, family Tombusviridae) and its molecular parasite, Satellite panicum mosaic virus (SPMV), synergistically interact in coinfected proso and pearl millet (Panicum miliaceum L.) plants resulting in a severe symptom phenotype. In this study, we examined synergistic interactions between the isolates of PMV and SPMV by using PMV-NE, PMV85, SPMV-KS, and SPMV-Type as interacting partner viruses in different combinations. Coinfection of proso millet plants by PMV-NE and SPMV-KS elicited severe mosaic, chlorosis, stunting, and eventual plant death compared with moderate mosaic, chlorotic streaks, and stunting by PMV85 and SPMV-Type. In reciprocal combinations, coinfection of proso millet by either isolate of PMV with SPMV-KS but not with SPMV-Type elicited severe disease synergism, suggesting that SPMV-KS was the main contributor for efficient synergistic interaction with PMV isolates. Coinfection of proso millet plants by either isolate of PMV and SPMV-KS or SPMV-Type caused increased accumulation of coat protein (CP) and genomic RNA copies of PMV, compared with infections by individual PMV isolates. Additionally, CP and genomic RNA copies of SPMV-KS accumulated at substantially higher levels, compared with SMPV-Type in coinfected proso millet plants with either isolate of PMV. Hybrid viruses between SPMV-KS and SPMV-Type revealed that SPMV isolates harboring a CP fragment with four differing amino acids at positions 18, 35, 59, and 98 were responsible for differential synergistic interactions with PMV in proso millet plants. Mutation of amino acid residues at these positions in different combinations in SPMV-KS, similar to those as in SPMV-Type or vice-versa, revealed that A35 and R98 in SPMV-KS CP play critical roles in enhanced synergistic interactions with PMV isolates. Taken together, these data suggest that the two distinct amino acids at positions 35 and 98 in the CP of SPMV-KS and SPMV-Type are involved in the differential synergistic interactions with the helper viruses.

The Tug-of-War between Plants and Viruses: Great Progress and Many Remaining Questions

Plants are persistently challenged by various phytopathogens. To protect themselves, plants have evolved multilayered surveillance against all pathogens. For intracellular parasitic viruses, plants have developed innate immunity, RNA silencing, translation repression, ubiquitination-mediated and autophagy-mediated protein degradation, and other dominant resistance gene-mediated defenses. Plant viruses have also acquired diverse strategies to suppress and even exploit host defense machinery to ensure their survival. A better understanding of the defense and counter-defense between plants and viruses will obviously benefit from the development of efficient and broad-spectrum virus resistance for sustainable agriculture. In this review, we summarize the cutting edge of knowledge concerning the defense and counter-defense between plants and viruses, and highlight the unexploited areas that are especially worth investigating in the near future.

Comprehensive identification of the full-length transcripts and alternative splicing related to the secondary metabolism pathways in the tea plant (Camellia sinensis)

Flavonoids, theanine and caffeine are the main secondary metabolites of the tea plant (Camellia sinensis), which account for the tea's unique flavor quality and health benefits. The biosynthesis pathways of these metabolites have been extensively studied at the transcriptional level, but the regulatory mechanisms are still unclear. In this study, to explore the transcriptome diversity and complexity of tea plant, PacBio Iso-Seq and RNA-seq analysis were combined to obtain full-length transcripts and to profile the changes in gene expression during the leaf development. A total of 1,388,066 reads of insert (ROI) were generated with an average length of 1,762 bp, and more than 54% (755,716) of the ROIs were full-length non-chimeric (FLNC) reads. The Benchmarking Universal Single-Copy Orthologue (BUSCO) completeness was 92.7%. A total of 93,883 non-redundant transcripts were obtained, and 87,395 (93.1%) were new alternatively spliced isoforms. Meanwhile, 7,650 differential expression transcripts (DETs) were identified. A total of 28,980 alternative splicing (AS) events were predicted, including 1,297 differential AS (DAS) events. The transcript isoforms of the key genes involved in the flavonoid, theanine and caffeine biosynthesis pathways were characterized. Additionally, 5,777 fusion transcripts and 9,052 long non-coding RNAs (lncRNAs) were also predicted. Our results revealed that AS potentially plays a crucial role in the regulation of the secondary metabolism of the tea plant. These findings enhanced our understanding of the complexity of the secondary metabolic regulation of tea plants and provided a basis for the subsequent exploration of the regulatory mechanisms of flavonoid, theanine and caffeine biosynthesis in tea plants.

Computational analysis of alternative splicing in plant genomes

Computational analyses play crucial roles in characterizing splicing isoforms in plant genomes. In this review, we provide a survey of computational tools used in recently published, genome-scale splicing analyses in plants. We summarize the commonly used software and pipelines for read mapping, isoform reconstruction, isoform quantification, and differential expression analysis. We also discuss methods for analyzing long reads and the strategies to combine long and short reads in identifying splicing isoforms. We review several tools for characterizing local splicing events, splicing graphs, coding potential, and visualizing splicing isoforms. We further discuss the procedures for identifying conserved splicing isoforms across plant species. Finally, we discuss the outlook of integrating other genomic data with splicing analyses to identify regulatory mechanisms of AS on genome-wide scale.

Cis-regulated alternative splicing divergence and its potential contribution to environmental responses in Arabidopsis

Alternative splicing (AS) plays key roles in plant development and the responses of plants to environmental changes. However, the mechanisms underlying AS divergence (differential expression of transcript isoforms resulting from AS) in plant accessions and its contribution to responses to environmental stimuli remain unclear. In this study, we investigated genome-wide variation of AS in Arabidopsis thaliana accessions Col-0, Bur-0, C24, Kro-0 and Ler-1, as well as their F₁ hybrids, and characterized the regulatory mechanisms for AS divergence by RNA sequencing. We found that most of the divergent AS events in Arabidopsis accessions were cis-regulated by sequence variation, including those in core splice site and splicing motifs. Many genes that differed in AS between Col-0 and Bur-0 were involved in stimulus responses. Further genome-wide association analyses of 22 environmental variables showed that single nucleotide polymorphisms influencing known splice site strength were also associated with environmental stress responses. These results demonstrate that cis-variation in genomic sequences among Arabidopsis accessions was the dominant contributor to AS divergence, and it may contribute to differences in environmental responses among Arabidopsis accessions.

Proteomic and Phosphoproteomic Analysis in Tobacco Mosaic Virus-Infected Tobacco (Nicotiana tabacum)

Tobacco mosaic virus (TMV) is a common source of biological stress that significantly affects plant growth and development. It is also useful as a model in studies designed to clarify the mechanisms involved in plant viral disease. Plant responses to abiotic stress were recently reported to be regulated by complex mechanisms at the post-translational modification (PTM) level. Protein phosphorylation is one of the most widespread and major PTMs in organisms. Using immobilized metal ion affinity chromatography (IMAC) enrichment, high-pH C18 chromatography fraction, and high-accuracy mass spectrometry (MS), a set of proteins and phosphopeptides in both TMV-infected tobacco and control tobacco were identified. A total of 4905 proteins and 3998 phosphopeptides with 3063 phosphorylation sites were identified. These 3998 phosphopeptides were assigned to 1311 phosphoproteins, as some proteins carried multiple phosphorylation sites. Among them, 530 proteins and 337 phosphopeptides corresponding to 277 phosphoproteins differed between the two groups. There were 43 upregulated phosphoproteins, including phosphoglycerate kinase, pyruvate phosphate dikinase, protein phosphatase 2C, and serine/threonine protein kinase. To the best of our knowledge, this is the first phosphoproteomic analysis of leaves from a tobacco cultivar, K326. The results of this study advance our understanding of tobacco development and TMV action at the protein phosphorylation level.

Translatomics: The Global View of Translation

In all kingdoms of life, proteins are synthesized by ribosomes in a process referred to as translation. The amplitude of translational regulation exceeds the sum of transcription, mRNA degradation and protein degradation. Therefore, it is essential to investigate translation in a global scale. Like the other “omics”-methods, translatomics investigates the totality of the components in the translation process, including but not limited to translating mRNAs, ribosomes, tRNAs, regulatory RNAs and nascent polypeptide chains. Technical advances in recent years have brought breakthroughs in the investigation of these components at global scale, both for their composition and dynamics. These methods have been applied in a rapidly increasing number of studies to reveal multifaceted aspects of translation control. The process of translation is not restricted to the conversion of mRNA coding sequences into polypeptide chains, it also controls the composition of the proteome in a delicate and responsive way. Therefore, translatomics has extended its unique and innovative power to many fields including proteomics, cancer research, bacterial stress response, biological rhythmicity and plant biology. Rational design in translation can enhance recombinant protein production for thousands of times. This brief review summarizes the main state-of-the-art methods of translatomics, highlights recent discoveries made in this field and introduces applications of translatomics on basic biological and biomedical research.

Alternative Splicing and Protein Diversity: Plants Versus Animals

Plants, unlike animals, exhibit a very high degree of plasticity in their growth and development and employ diverse strategies to cope with the variations during diurnal cycles and stressful conditions. Plants and animals, despite their remarkable morphological and physiological differences, share many basic cellular processes and regulatory mechanisms. Alternative splicing (AS) is one such gene regulatory mechanism that modulates gene expression in multiple ways. It is now well established that AS is prevalent in all multicellular eukaryotes including plants and humans. Emerging evidence indicates that in plants, as in animals, transcription and splicing are coupled. Here, we reviewed recent evidence in support of co-transcriptional splicing in plants and highlighted similarities and differences between plants and humans. An unsettled question in the field of AS is the extent to which splice isoforms contribute to protein diversity. To take a critical look at this question, we presented a comprehensive summary of the current status of research in this area in both plants and humans, discussed limitations with the currently used approaches and suggested improvements to current methods and alternative approaches. We end with a discussion on the potential role of epigenetic modifications and chromatin state in splicing memory in plants primed with stresses.

Segmental and Tandem Duplications Driving the Recent NBS-LRR Gene Expansion in the Asparagus Genome

Garden asparagus is an important horticultural plant worldwide. It is, however, susceptible to a variety of diseases, which can affect the potential yield, spear quality, and lifespan of production fields. Screening studies have identified resistant germplasm. The genetic resistance is usually complex, and the genes underlying that resistance are still unknown. Most often, disease resistance is determined by resistance genes (R). The most predominant R-genes contain nucleotide binding site and leucine-rich repeat (NBS-LRR) domains. Using bioinformatics and data mining approaches, we identified and characterized 68 NBS predicted proteins encoded by 49 different loci in the asparagus genome. The NBS-encoding genes were grouped into seven distinct classes based on their domain architecture. The NBS genes are unevenly distributed through the genome and nearly 50% of the genes are present in clusters. Chromosome 6 is significantly NBS-enriched and one single cluster hosts 10% of the genes. Phylogenetic analysis points to their diversification into three families during their evolution. Recent duplications are likely to have dominated the NBS expansion with both tandem genes and duplication events across multiple chromosomes. Transcriptome sequencing data provided evidence for their transcription and tissue-specific expression. The total number of cis-regulatory elements as well as their relative positions within the NBS promoters suggests a complex transcriptional network regulating defense responses. Our study provides a strong groundwork for the isolation of candidate R-genes in garden asparagus.

Chromosome-level reference genome and alternative splicing atlas of moso bamboo (Phyllostachys edulis)

Background

Bamboo is one of the most important nontimber forestry products worldwide. However, a chromosome-level reference genome is lacking, and an evolutionary view of alternative splicing (AS) in bamboo remains unclear despite emerging omics data and improved technologies.

Results

Here, we provide a chromosome-level de novo genome assembly of moso bamboo (Phyllostachys edulis) using additional abundance sequencing data and a Hi-C scaffolding strategy. The significantly improved genome is a scaffold N50 of 79.90 Mb, approximately 243 times longer than the previous version. A total of 51,074 high-quality protein-coding loci with intact structures were identified using single-molecule real-time sequencing and manual verification. Moreover, we provide a comprehensive AS profile based on the identification of 266,711 unique AS events in 25,225 AS genes by large-scale transcriptomic sequencing of 26 representative bamboo tissues using both the Illumina and Pacific Biosciences sequencing platforms. Through comparisons with orthologous genes in related plant species, we observed that the AS genes are concentrated among more conserved genes that tend to accumulate higher transcript levels and share less tissue specificity. Furthermore, gene family expansion, abundant AS, and positive selection were identified in crucial genes involved in the lignin biosynthetic pathway of moso bamboo.

Conclusions

These fundamental studies provide useful information for future in-depth analyses of comparative genome and AS features. Additionally, our results highlight a global perspective of AS during evolution and diversification in bamboo.

Potato Spindle Tuber Viroid Modulates Its Replication through a Direct Interaction with a Splicing Regulator

Viroids are circular noncoding RNAs (ncRNAs) that infect plants. Despite differences in the genetic makeup and biogenesis, viroids and various long ncRNAs all rely on RNA structure-based interactions with cellular factors for function. Viroids replicating in the nucleus utilize DNA-dependent RNA polymerase II for transcription, a process that involves a unique splicing form of transcription factor IIIA (TFIIIA-7ZF). Here, we provide evidence showing that potato spindle tuber viroid (PSTVd) interacts with a TFIIIA splicing regulator (ribosomal protein L5 [RPL5]) in vitro and in vivo PSTVd infection compromises the regulatory role of RPL5 over splicing of TFIIIA transcripts, while ectopic expression of RPL5 reduces TFIIIA-7ZF expression and attenuates PSTVd accumulation. Furthermore, we illustrate that the RPL5 binding site on the PSTVd genome resides in the central conserved region critical for replication. Together, our data suggest that viroids can regulate their own replication and modulate specific regulatory factors leading to splicing changes in only one or a few genes. This study also has implications for understanding the functional mechanisms of ncRNAs and elucidating the global splicing changes in various host-pathogen interactions.IMPORTANCE Viroids are the smallest replicons among all living entities. As circular noncoding RNAs, viroids can replicate and spread in plants, often resulting in disease symptoms. Potato spindle tuber viroid (PSTVd), the type species of nuclear-replicating viroids, requires a unique splicing form of transcription factor IIIA (TFIIIA-7ZF) for its propagation. Here, we provide evidence showing that PSTVd directly interacts with a splicing regulator, RPL5, to favor the expression of TFIIIA-7ZF, thereby promoting viroid replication. This finding provides new insights to better understand viroid biology and sheds light on the noncoding RNA-based regulation of splicing. Our discovery also establishes RPL5 as a novel negative factor regulating viroid replication in the nucleus and highlights a potential means for viroid control.

Brachypodium: A Monocot Grass Model Genus for Plant Biology

The genus Brachypodium represents a model system that is advancing our knowledge of the biology of grasses, including small grains, in the postgenomics era. The most widely used species, Brachypodium distachyon, is a C3 plant that is distributed worldwide. B. distachyon has a small genome, short life cycle, and small stature and is amenable to genetic transformation. Due to the intensive and thoughtful development of this grass as a model organism, it is well-suited for laboratory and field experimentation. The intent of this review is to introduce this model system genus and describe some key outcomes of nearly a decade of research since the first draft genome sequence of the flagship species, B. distachyon, was completed. We discuss characteristics and features of B. distachyon and its congeners that make the genus a valuable model system for studies in ecology, evolution, genetics, and genomics in the grasses, review current hot topics in Brachypodium research, and highlight the potential for future analysis using this system in the coming years.

Transcriptome analysis of alternative splicing in peanut (Arachis hypogaea L.)

BACKGROUND

Alternative splicing (AS) represents a mechanism widely used by eukaryotes for the post-transcriptional regulation of genes. The detailed exploration of AS in peanut has not been documented.

RESULTS

The strand-specific RNA-Seq technique was exploited to characterize the distribution of AS in the four samples of peanut (FH1-seed1, FH1-seed2, FH1-root and FH1-leaf). AS was detected as affecting around 37.2% of the full set of multi-exon genes. Some of these genes experienced AS throughout the plant, while in the case of others, the effect was organ-specific. Overall, AS was more frequent in the seed than in either the root or leaf. The predominant form of AS was intron retention, and AS in transcription start site and transcription terminal site were commonly identified in all the four samples. It is interesting that in genes affected by AS, the majority experienced only a single type of event. Not all of the in silico predicted transcripts appeared to be translated, implying that these are either degraded or sequestered away from the translation machinery. With respect to genes involved in fatty acid metabolism, about 61.6% were shown to experience AS.

CONCLUSION

Our report contributes significantly in AS analysis of peanut genes in general, and these results have not been mentioned before. The specific functions of different AS forms need further investigation.

Application of Brachypodium genotypes to the analysis of type II resistance to Fusarium head blight (FHB)

The resistance to Fusarium head blight (FHB) in wheat is mainly via the restrain of fungal expansion through spike rachis (type II resistance). In order to unravel the resistance mechanisms, Brachypodium distachyon 21 (Bd21), a monocotyledonous model plant, was previously proved to interact with F. graminearum, while the disease development in spike still needs to be explored in detail. Herein, it is found that the fungal spores mainly germinate on pistil of Bd21, then the hyphae rapidly extend to the bottom of floret and enter spike rachis, similar with the infection progress in wheat. However, structural difference of spike rachis was found between Brachypodium and wheat. It was found that the spread of the fungus through the rachis node of inoculated spikelets is an important index for the evaluation of type II FHB resistance in Brachypodium under optimal conditions at 28 °C and 50%-70% humidity. To verify the feasibility of this strategy, the transcription factor TaTGA2 was overexpressed in Bd21, and transgenic plants were found to show improved resistance to F. graminearum in both spikes and detached leaves, which was further supported by the increased disease severity when silencing TaTGA2 in the wheat cultivar "Sumai 3" or in tilling "Kronos" mutants. Except for Bd21, another 49 Brachypodium germplasms were further screened for FHB resistance, and three moderately susceptible germplasms, namely, PI 317418, W6-39284, and PI 254868, feasible for transformation, were determined to be better hosts than Bd21 when evaluating heterologous genes that positively regulate FHB resistance. The present study also observed variations in the levels of FHB resistance between coleoptiles and spikes or transgenic plants and natural germplasms.

Genome-Wide Association Analyses Reveal the Importance of Alternative Splicing in Diversifying Gene Function and Regulating Phenotypic Variation in Maize

Alternative splicing (AS) enhances transcriptome diversity and plays important roles in regulating plant processes. Although widespread natural variation in AS has been observed in plants, how AS is regulated and contribute to phenotypic variation is poorly understood. Here, we report a population-level transcriptome assembly and genome-wide association study to identify splicing quantitative trait loci (sQTLs) in developing maize (Zea mays) kernels from 368 inbred lines. We detected 19,554 unique sQTLs for 6570 genes. Most sQTLs showed small isoform usage changes without involving major isoform switching between genotypes. The sQTL-affected isoforms tend to display distinct protein functions. We demonstrate that nonsense-mediated mRNA decay, microRNA-mediated regulation, and small interfering peptide-mediated peptide interference are frequently involved in sQTL regulation. The natural variation in AS and overall mRNA level appears to be independently regulated with different cis-sequences preferentially used. We identified 214 putative trans-acting splicing regulators, among which ZmGRP1, encoding an hnRNP-like glycine-rich RNA binding protein, regulates the largest trans-cluster. Knockout of ZmGRP1 by CRISPR/Cas9 altered splicing of numerous downstream genes. We found that 739 sQTLs colocalized with previous marker-trait associations, most of which occurred without changes in overall mRNA level. Our findings uncover the importance of AS in diversifying gene function and regulating phenotypic variation.

Transcriptome profiling using Illumina- and SMRT-based RNA-seq of hot pepper for in-depth understanding of genes involved in CMV infection

Hot pepper (Capsicum annuum L.) is becoming an increasingly important vegetable crop in the world. Cucumber mosaic virus (CMV) is a destructive virus that can cause leaf distortion and fruit lesions, affecting pepper production. However, studies on the response to CMV infection in pepper at the transcriptional level are limited. In this study, the transcript profiles of pepper leaves after CMV infection were investigated using Illumina and single-molecule real-time (SMRT) RNA-sequencing (RNA-seq). A total of 2143 differentially expressed genes (DEGs) were identified at five different stages. Gene ontology (GO) and KEGG analysis revealed that these DEGs were involved in the response to stress, defense response and plant-pathogen interaction pathways. Among these DEGs, several key genes that consistently appeared in studies of plant-pathogen interactions had increased transcript abundance after inoculation, including chitinase, pathogenesis-related (PR) protein, TMV resistance protein, WRKY transcription factor and jasmonate ZIM-domain protein. Four of these DEGs were further validated by quantitative real-time RT-PCR (qRT-PCR). Furthermore, a total of 73, 597 alternative splicing (AS) events were identified in the pepper leaves after CMV infection, distributed in 12, 615 genes. The intron retention of WRKY33 (Capana09g001251) might be involved in the regulation of CMV infection. Taken together, our study provides a transcriptome-wide insight into the molecular basis of resistance to CMV infection in pepper leaves and potential candidate genes for improving resistance cultivars.

Transcript-level expression control of plant NLR genes

Plant NLR genes encode sensitive immune receptors that can mediate the specific recognition of pathogen avirulence effectors and activate a strong defence response, termed effector-triggered immunity. The expression of NLRs requires strict regulation, as their ability to trigger immunity is dependent on their dose, and overexpression of NLRs results in autoimmunity and massive fitness costs. An elaborate interplay of different mechanisms controlling NLR transcript levels allows plants to maximize their defence capacity, whilst limiting negative impact on their fitness. Global suppression of NLR transcripts may be a prerequisite for the fast evolution of new NLR variants and the expansion of this gene family. Here, we summarize recent progress made towards a comprehensive understanding of NLR transcript-level expression control. Multiple mechanistic steps, including transcription as well as co-/post-transcriptional processing and transcript turn-over, contribute to balanced base levels of NLR transcripts and allow for dynamic adjustments to defence situations.

Subcellular Compartmentation of Alternatively Spliced Transcripts Defines SERINE/ARGININE-RICH PROTEIN30 Expression

Alternative splicing (AS) is prevalent in higher eukaryotes, and generation of different AS variants is tightly regulated. Widespread AS occurs in response to altered light conditions and plays a critical role in seedling photomorphogenesis, but despite its frequency and effect on plant development, the functional role of AS remains unknown for most splicing variants. Here, we characterized the light-dependent AS variants of the gene encoding the splicing regulator Ser/Arg-rich protein SR30 in Arabidopsis (Arabidopsis thaliana). We demonstrated that the splicing variant SR30.2, which is predominantly produced in darkness, is enriched within the nucleus and strongly depleted from ribosomes. Light-induced AS from a downstream 3' splice site gives rise to SR30.1, which is exported to the cytosol and translated, coinciding with SR30 protein accumulation upon seedling illumination. Constitutive expression of SR30.1 and SR30.2 fused to fluorescent proteins revealed their identical subcellular localization in the nucleoplasm and nuclear speckles. Furthermore, expression of either variant shifted splicing of a genomic SR30 reporter toward SR30.2, suggesting that an autoregulatory feedback loop affects SR30 splicing. We provide evidence that SR30.2 can be further spliced and, unlike SR30.2, the resulting cassette exon variant SR30.3 is sensitive to nonsense-mediated decay. Our work delivers insight into the complex and compartmentalized RNA processing mechanisms that control the expression of the splicing regulator SR30 in a light-dependent manner.

De novo generation of helper virus-satellite chimera RNAs results in disease attenuation and satellite sequence acquisition in a host-dependent manner

Panicum mosaic virus (PMV) is a helper RNA virus for satellite RNAs (satRNAs) and a satellite virus (SPMV). Here, we describe modifications that occur at the 3'-end of a satRNA of PMV, satS. Co-infections of PMV+satS result in attenuation of the disease symptoms induced by PMV alone in Brachypodium distachyon and proso millet. The 375 nt satS acquires ~100-200 nts from the 3'-end of PMV during infection and is associated with decreased abundance of the PMV RNA and capsid protein in millet. PMV-satS chimera RNAs were isolated from native infections of St. Augustinegrass and switchgrass. Phylogenetic analyses revealed that the chimeric RNAs clustered according to the host species from which they were isolated. Additionally, the chimera satRNAs acquired non-viral "linker" sequences in a host-specific manner. These results highlight the dynamic regulation of viral pathogenicity by satellites, and the selective host-dependent, sequence-based pressures for driving satRNA generation and genome compositions.

Integrative Analysis of Three RNA Sequencing Methods Identifies Mutually Exclusive Exons of MADS-Box Isoforms During Early Bud Development in Picea abies

Recent efforts to sequence the genomes and transcriptomes of several gymnosperm species have revealed an increased complexity in certain gene families in gymnosperms as compared to angiosperms. One example of this is the gymnosperm sister clade to angiosperm TM3-like MADS-box genes, which at least in the conifer lineage has expanded in number of genes. We have previously identified a member of this sub-clade, the conifer gene DEFICIENS AGAMOUS LIKE 19 (DAL19), as being specifically upregulated in cone-setting shoots. Here, we show through Sanger sequencing of mRNA-derived cDNA and mapping to assembled conifer genomic sequences that DAL19 produces six mature mRNA splice variants in Picea abies. These splice variants use alternate first and last exons, while their four central exons constitute a core region present in all six transcripts. Thus, they are likely to be transcript isoforms. Quantitative Real-Time PCR revealed that two mutually exclusive first DAL19 exons are differentially expressed across meristems that will form either male or female cones, or vegetative shoots. Furthermore, mRNA in situ hybridization revealed that two mutually exclusive last DAL19 exons were expressed in a cell-specific pattern within bud meristems. Based on these findings in DAL19, we developed a sensitive approach to transcript isoform assembly from short-read sequencing of mRNA. We applied this method to 42 putative MADS-box core regions in P. abies, from which we assembled 1084 putative transcripts. We manually curated these transcripts to arrive at 933 assembled transcript isoforms of 38 putative MADS-box genes. 152 of these isoforms, which we assign to 28 putative MADS-box genes, were differentially expressed across eight female, male, and vegetative buds. We further provide evidence of the expression of 16 out of the 38 putative MADS-box genes by mapping PacBio Iso-Seq circular consensus reads derived from pooled sample sequencing to assembled transcripts. In summary, our analyses reveal the use of mutually exclusive exons of MADS-box gene isoforms during early bud development in P. abies, and we find that the large number of identified MADS-box transcripts in P. abies results not only from expansion of the gene family through gene duplication events but also from the generation of numerous splice variants.

Genomic Approaches to Analyze Alternative Splicing, A Key Regulator of Transcriptome and Proteome Diversity in Brachypodium distachyon

Alternative splicing (AS) promotes transcriptome and proteome diversity in plants, which influences growth and development, and host responses to stress. Advancements in next-generation sequencing, bioinformatics, and computational biology tools have allowed biologists to investigate AS landscapes on a genome-wide scale in several plant species. Furthermore, the development of Brachypodium distachyon (Brachypodium) as a model system for grasses has facilitated comparative studies of AS within the Poaceae. These analyses revealed a plethora of genes in several biological processes that are alternatively spliced and identified conserved AS patterns among monocot and dicot plants. In this chapter, using a Brachypodium-virus pathosystem as a research template, we provide an overview of genomic and bioinformatic tools that can be used to investigate constitutive and alternative splicing in plants.

Transcriptome analysis of Brachypodium during fungal pathogen infection reveals both shared and distinct defense responses with wheat

Fusarium crown rot (FCR) of wheat and barley, predominantly caused by the fungal pathogen Fusarium pseudograminearum, is a disease of economic significance. The quantitative nature of FCR resistance within cultivated wheat germplasm has significantly limited breeding efforts to enhanced FCR resistance in wheat. In this study, we characterized the molecular responses of Brachypodium distachyon (Brachypodium hereafter) to F. pseudograminearum infection using RNA-seq to determine whether Brachypodium can be exploited as a model system towards better understanding of F. pseudograminearum-wheat interaction. The transcriptional response to infection in Brachypodium was strikingly similar to that previously reported in wheat, both in shared expression patterns of wheat homologs of Brachypodium genes and functional overlap revealed through comparative gene ontology analysis in both species. Metabolites produced by various biosynthetic pathways induced in both wheat and Brachypodium were quantified, revealing a high degree of overlap between these two species in metabolic response to infection but also showed Brachypodium does not produce certain defence-related metabolites found in wheat. Functional analyses of candidate genes identified in this study will improve our understanding of resistance mechanisms and may lead to the development of new strategies to protect cereal crops from pathogen infection.

Long noncoding RNAs in the model species Brachypodium distachyon

Eukaryotic genomes are pervasively transcribed and only a small portion of the transcribed sequences belongs to protein coding genes. High-throughput sequencing technology contributed to consolidate this perspective, allowing the identification of numerous noncoding RNAs with key roles in biological processes. Long noncoding RNAs (lncRNAs) are transcripts longer than 200 nt with limited phylogenetic conservation, expressed at low levels and characterized by tissue/organ specific expression profiles. Although a large set of lncRNAs has been identified, the functional roles of lncRNAs are only beginning to be recognized and the molecular mechanism of lncRNA-mediated gene regulation remains largely unexplored, particularly in plants where their annotation and characterization are still incomplete. Using public and proprietary poly-(A)+ RNA-seq data as well as a collection of full length ESTs from several organs, developmental stages and stress conditions in three Brachypodium distachyon inbred lines, we describe the identification and the main features of thousands lncRNAs. Here we provide a genome-wide characterization of lncRNAs, highlighting their intraspecies conservation and describing their expression patterns among several organs/tissues and stress conditions. This work represents a fundamental resource to deepen our knowledge on long noncoding RNAs in C3 cereals, allowing the Brachypodium community to exploit these results in future research programs.

Comprehensive transcriptome analyses reveal tomato plant responses to tobacco rattle virus-based gene silencing vectors

In plants, virus-induced gene silencing (VIGS) is a popular tool for functional genomic studies or rapidly assessing individual gene functions. However, molecular details regarding plant responses to viral vectors remain elusive, which may complicate experimental designs and data interpretation. To this end, we documented whole transcriptome changes of tomato elicited by the application of the most widely used tobacco rattle virus (TRV)-based vectors, using comprehensive genome-wide analyses. Our data illustrated multiple biological processes with functional implications, including (1) the enhanced activity of miR167 in guiding the cleavage of an auxin response factor; (2) reduced accumulation of phased secondary small interfering RNAs from two genomic loci; (3) altered expression of ~500 protein-coding transcripts; and (4) twenty long noncoding RNAs specifically responsive to TRV vectors. Importantly, we unraveled large-scale changes in mRNA alternative splicing patterns. These observations will facilitate future application of VIGS vectors for functional studies benefiting the plant research community and help deepen the understanding of plant-virus interactions.