New gene models and alternative splicing in the maize pathogen Colletotrichum graminicola revealed by RNA-Seq analysis

Highly efficient gene knockout system in the maize pathogen Colletotrichum graminicola using Agrobacterium tumefaciens-mediated transformation (ATMT)

Colletotrichum graminicola, a hemibiotrophic pathogenic fungus, is the causal agent of anthracnose of maize, which causes significant yield losses worldwide, especially in warm and humid maize production regions. An efficient targeted genes knockout protocol is crucial to explore molecular mechanisms of fungal virulence to the host. In this study, we established a gene knockout transformation system by employing Agrobacterium tumefaciens-mediated transformation to knockout genes in M 1.001 strain of C. graminicola. The conidia germination status, induction medium type, and ratio of Agrobacterium cell and conidia suspension were optimized for the knockout of CgBRN1(OR352905), a gene relating to the fungal melanin biosynthesis pathway. Additionally, CgPKS18 (OR352906) and CgCDC25 (OR352903) were knocked out to test the applicability of the gene knockout transformation system. In this established system, transformation efficiency was 176 transformants per 1 × 105 conidia and the homologous recombination efficiency was 53.3 to 75%. Furthermore, disease index, lesion number and lesion size caused by the three above-mentioned mutant strains were found to be reduced significantly compared to the wild-type strain, which indicated reduction in fungal virulence due to the lack of those genes.

Comprehensive transcriptional analysis of pig facial skin development

Background

Skin development is a complex process that is influenced by many factors. Pig skin is used as an ideal material for xenografts because it is more anatomically and physiologically similar to human skin. It has been shown that the skin development of different pig breeds is different, and some Chinese pig breeds have the characteristics of skin thickness and facial skin folds, but the specific regulatory mechanism of this skin development is not yet clear.

Methods

In this study, the facial skin of Chenghua sows in the four developmental stages of postnatal Day 3 (D3) , Day 90 (D90) , Day 180 (D180), and Year 3 (Y3) were used as experimental materials, and RNA sequencing (RNA-seq) analysis was used to explore the changes in RNA expression in skin development at the four developmental stages, determine the differentially expressed messenger RNAs (mRNAs), long noncoding RNAs (lncRNAs), microRNAs (miRNAs), and circular RNAs (circRNAs), and perform functional analysis of related genes by Gene Ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses.

Results

A pairwise comparison of the four developmental stages identified several differentially expressed genes (DEGs) and found that the number of differentially expressed RNAs (DE RNAs) increased with increasing developmental time intervals. Elastin (ELN) is an important component of the skin. Its content affects the relaxation of the epidermis and dermal connection, and its expression is continuously downregulated during the four developmental stages. The functions of DEGs at different developmental stages were examined by performing GO and KEGG analyses, and the GO terms and enrichment pathways of mRNAs, lncRNAs, miRNAs, and circRNAs highly overlapped, among which the PPAR signaling pathway, a classical pathway for skin development, was enriched by DEGs of D3 vs. D180, D90 vs. D180 and D180 vs. Y3. In addition, we constructed lncRNA-miRNA-mRNA and circRNA-miRNA interaction networks and found genes that may be associated with skin development, but their interactions need further study.

Conclusions

We identified a number of genes associated with skin development, performed functional analyses on some important DEGs and constructed interaction networks that facilitate further studies of skin development.

Comparative transcriptomics and genomic analyses reveal differential gene expression related to Colletotrichum brevisporum resistance in papaya (Carica papaya L.)

Colletotrichum brevisporum is an important causal pathogen of anthracnose that seriously affects the fruit quality and yield of papaya (Carica papaya L.). Although many genes and biological processes involved in anthracnose resistance have been reported in other species, the molecular mechanisms involved in the response or resistance to anthracnose in post-harvest papaya fruits remain unclear. In this study, we compared transcriptome changes in the post-harvest fruits of the anthracnose-susceptible papaya cultivar Y61 and the anthracnose-resistant cultivar G20 following C. brevisporum inoculation. More differentially expressed genes (DEGs) and differentially expressed long non-coding RNAs (DElnRNAs) were identified in G20 than in Y61, especially at 24 h post-inoculation (hpi), suggesting a prompt activation of defense responses in G20 in the first 24 h after C. brevisporum inoculation. These DEGs were mainly enriched in plant-pathogen interaction, phenylpropanoid biosynthesis/metabolism, and peroxisome and flavonoid biosynthesis pathways in both cultivars. However, in the first 24 hpi, the number of DEGs related to anthracnose resistance was greater in G20 than in Y61, and changes in their expression levels were faster in G20 than in Y61. We also identified a candidate anthracnose-resistant gene cluster, which consisted of 12 genes, 11 in G20 and Y61, in response to C. brevisporum inoculation. Moreover, 529 resistance gene analogs were identified in papaya genome, most of which responded to C. brevisporum inoculation and were genetically different between papaya cultivars and wild-type populations. The total expression dose of the resistance gene analogs may help papaya resist C. brevisporum infection. This study revealed the mechanisms underlying different anthracnose resistance between the anthracnose-resistant and anthracnose-susceptible cultivars based on gene expression, and identified some potential anthracnose resistance-related candidate genes/major regulatory factors. Our findings provided potential targets for developing novel genetic strategies to overcome anthracnose in papaya.

Maize Anthracnose Stalk Rot in the Genomic Era

Anthracnose stalk rot (ASR) of maize results in millions of dollars in losses annually in the United States. ASR, together with anthracnose leaf blight and anthracnose top dieback, is caused by the fungus Colletotrichum graminicola. Current ASR management recommendations emphasize host resistance and reduction of plant stressors (e.g., drought, heat, low fertility, or soil acidity). Stress reduction may be more difficult to achieve in the future due to more high-intensity production protocols and climate change. Moreover, cultural and chemical management practices may conflict with other important goals, including environmental sustainability and maximization of yield potential. Thus, future ASR management may rely more heavily on host resistance, for which there are relatively few highly effective sources. The last comprehensive review of C. graminicola and maize anthracnose was written over two decades ago. The genomic age has brought important new insights into mechanisms governing the host-pathogen interaction from the application of molecular and cytological technologies. This review provides a summary of our current model of maize anthracnose etiology, including how increased knowledge of molecular and cellular events could contribute to better ASR management. Improved understanding of C. graminicola taxonomy has confirmed that the fungus is specific to Zea mays, and that it colonizes living maize tissues via a critical biotrophic phase. Successful biotrophic establishment relies on an array of secreted protein effectors and secondary metabolites produced at different stages of infection and dispersed to multiple locations. These molecules could provide therapeutic targets for the next generation of transgenic or gene-edited ASR-resistant hybrids.

Alternative splicing diversifies the transcriptome and proteome of the rice blast fungus during host infection

Alternative splicing (AS) contributes to diversifying and regulating cellular responses to environmental conditions and developmental cues by differentially producing multiple mRNA and protein isoforms from a single gene. Previous studies on AS in pathogenic fungi focused on profiling AS isoforms under a limited number of conditions. We analysed AS profiles in the rice blast fungus Magnaporthe oryzae, a global threat to rice production, using high-quality transcriptome data representing its vegetative growth (mycelia) and multiple host infection stages. We identified 4,270 AS isoforms derived from 2,413 genes, including 499 genes presumably regulated by infection-specific AS. AS appears to increase during infection, with 32.7% of the AS isoforms being produced during infection but absent in mycelia. Analysis of the isoforms observed at each infection stage showed that 636 AS isoforms were more abundant than corresponding annotated mRNAs, especially after initial hyphal penetration into host cell. Many such dominant isoforms were predicted to encode regulatory proteins such as transcription factors and phospho-transferases. We also identified the genes encoding distinct proteins via AS and confirmed the translation of some isoforms via a proteomic analysis, suggesting potential AS-mediated neo-functionalization of some genes during infection. Comprehensive profiling of the pattern of genome-wide AS during multiple stages of rice-M. oryzae interaction established a foundational resource that will help investigate the role and regulation of AS during rice infection.

Full-Length Transcriptome Analysis Reveals Candidate Genes Involved in Terpenoid Biosynthesis in Artemisia argyi

Artemisia argyi is an important medicinal plant widely utilized for moxibustion heat therapy in China. The terpenoid biosynthesis process in A. argyi is speculated to play a key role in conferring its medicinal value. However, the molecular mechanism underlying terpenoid biosynthesis remains unclear, in part because the reference genome of A. argyi is unavailable. Moreover, the full-length transcriptome of A. argyi has not yet been sequenced. Therefore, in this study, de novo transcriptome sequencing of A. argyi's root, stem, and leaf tissues was performed to obtain those candidate genes related to terpenoid biosynthesis, by combining the PacBio single-molecule real-time (SMRT) and Illumina sequencing NGS platforms. And more than 55.4 Gb of sequencing data and 108,846 full-length reads (non-chimeric) were generated by the Illumina and PacBio platform, respectively. Then, 53,043 consensus isoforms were clustered and used to represent 36,820 non-redundant transcripts, of which 34,839 (94.62%) were annotated in public databases. In the comparison sets of leaves vs roots, and leaves vs stems, 13,850 (7,566 up-regulated, 6,284 down-regulated) and 9,502 (5,284 up-regulated, 4,218 down-regulated) differentially expressed transcripts (DETs) were obtained, respectively. Specifically, the expression profile and KEGG functional enrichment analysis of these DETs indicated that they were significantly enriched in the biosynthesis of amino acids, carotenoids, diterpenoids and flavonoids, as well as the metabolism processes of glycine, serine and threonine. Moreover, multiple genes encoding significant enzymes or transcription factors related to diterpenoid biosynthesis were highly expressed in the A. argyi leaves. Additionally, several transcription factor families, such as RLK-Pelle_LRR-L-1 and RLK-Pelle_DLSV, were also identified. In conclusion, this study offers a valuable resource for transcriptome information, and provides a functional genomic foundation for further research on molecular mechanisms underlying the medicinal use of A. argyi leaves.

Transcriptome profiling of Capsicum annuum using Illumina- and PacBio SMRT-based RNA-Seq for in-depth understanding of genes involved in trichome formation

Trichomes, specialized epidermal cells located in aerial parts of plants, play indispensable roles in resisting abiotic and biotic stresses. However, the regulatory genes essential for multicellular trichrome development in Capsicum annuum L. (pepper) remain unclear. In this study, the transcript profiles of peppers GZZY-23 (hairy) and PI246331 (hairless) were investigated to gain insights into the genes responsible for the formation of multicellular trichomes. A total of 40,079 genes, including 4743 novel genes and 13,568 differentially expressed genes (DEGs), were obtained. Functional enrichment analysis revealed that the most noticeable pathways were transcription factor activity, sequence-specific DNA binding, and plant hormone signal transduction, which might be critical for multicellular trichome formation in hairy plants. We screened 11 DEGs related to trichome development; 151 DEGs involved in plant hormone signal transduction; 312 DEGs belonging to the MYB, bHLH, HD-Zip, and zinc finger transcription factor families; and 1629 DEGs predicted as plant resistance genes (PRGs). Most of these DEGs were highly expressed in GZZY-23 or trichomes. Several homologs of trichome regulators, such as SlCycB2, SlCycB3, and H, were considerably upregulated in GZZY-23, especially in the trichomes. The transcriptomic data generated in this study provide a basis for future characterization of trichome formation in pepper.

Comprehensive transcriptome profiling reveals abundant long non-coding RNAs associated with development of the rice false smut fungus, Ustilaginoidea virens

Long non-coding RNAs (lncRNAs) play an important role in biological processes but regulation and function of lncRNAs remain largely unelucidated, especially in fungi. Ustilaginoidea virens is an economically important fungus causing a devastating disease of rice. By combining microscopic and RNA-seq analyses, we comprehensively characterized lncRNAs of this fungus in infection and developmental processes and defined four serial typical stages. RNA-seq analyses revealed 1724 lncRNAs in U. virens, including 1084 long intergenic non-coding RNAs (lincRNAs), 51 intronic RNAs (incRNAs), 566 natural antisense transcripts (lncNATs) and 23 sense transcripts. Gene Ontology enrichment of differentially expressed lincRNAs and lncNATs demonstrated that these were mainly involved in transport-related regulation. Functional studies of transport-related lncRNAs revealed that UvlncNAT-MFS, a cytoplasm localized lncNAT of a putative MFS transporter gene, UvMFS, could form an RNA duplex with UvMFS and was required for regulation of growth, conidiation and various stress responses. Our results were the first to elucidate the lncRNA profiles during infection and development of this important phytopathogen U. virens. The functional discovery of the novel lncRNA, UvlncNAT-MFS, revealed the potential of lncRNAs in regulation of life processes in fungi.

Genome-wide alternative splicing profiling in the fungal plant pathogen Sclerotinia sclerotiorum during the colonization of diverse host families

Sclerotinia sclerotiorum is a notorious generalist plant pathogen that threatens more than 600 host plants, including wild and cultivated species. The molecular bases underlying the broad compatibility of S. sclerotiorum with its hosts is not fully elucidated. In contrast to higher plants and animals, alternative splicing (AS) is not well studied in plant-pathogenic fungi. AS is a common regulated cellular process that increases cell protein and RNA diversity. In this study, we annotated spliceosome genes in the genome of S. sclerotiorum and characterized their expression in vitro and during the colonization of six host species. Several spliceosome genes were differentially expressed in planta, suggesting that AS was altered during infection. Using stringent parameters, we identified 1,487 S. sclerotiorum genes differentially expressed in planta and exhibiting alternative transcripts. The most common AS events during the colonization of all plants were retained introns and the alternative 3' receiver site. We identified S. sclerotiorum genes expressed in planta for which (a) the relative accumulation of alternative transcripts varies according to the host being colonized and (b) alternative transcripts harbour distinct protein domains. This notably included 42 genes encoding predicted secreted proteins showing high-confidence AS events. This study indicates that AS events are taking place in the plant pathogenic fungus S. sclerotiorum during the colonization of host plants and could generate functional diversity in the repertoire of proteins secreted by S. sclerotiorum during infection.

Physiological RNA dynamics in RNA-Seq analysis

Physiological RNA dynamics cause problems in transcriptome analysis. Physiological RNA accumulation affects the analysis of RNA quantification, and physiological RNA degradation affects the analysis of the RNA sequence length, feature site and quantification. In the present article, we review the effects of physiological degradation and accumulation of RNA on analysing RNA sequencing data. Physiological RNA accumulation and degradation probably led to such phenomena as incorrect estimations of transcription quantification, differential expressions, co-expressions, RNA decay rates, alternative splicing, boundaries of transcription, novel genes, new single-nucleotide polymorphisms, small RNAs and gene fusion. Thus, the transcriptomic data obtained up to date warrant further scrutiny. New and improved techniques and bioinformatics software are needed to produce accurate data in transcriptome research.

Transcriptomic Analysis Reveals the Involvement of lncRNA-miRNA-mRNA Networks in Hair Follicle Induction in Aohan Fine Wool Sheep Skin

Long non-coding RNAs (lncRNA) and microRNAs (miRNA) are new found classes of non-coding RNAs (ncRNAs) that are not translated into proteins but regulate various cellular and biological processes. In this study, we conducted a transcriptomic analysis of ncRNA and mRNA expression in Aohan fine wool sheep (AFWS) at different growth stages (embryonic day 90, embryonic day 120, and the day of birth), and explored their relationship with wool follicle growth. In total, 461 lncRNAs, 106 miRNAs, and 1,009 mRNAs were found to be differentially expressed during the three stages of wool follicle development. Gene Ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis were performed to clarify the roles of the differentially expressed lncRNA, miRNA, and mRNA in the different stages of wool follicle development. Quantitative real-time PCR (qRT-PCR) was used to validate the results of RNA-seq analysis. lncRNA (MSTRG.223165) was found to act as a competing endogenous RNA (ceRNA) and may participate in wool follicle development by acting as an miR-21 sponge. Network prediction implicated the MSTRG.223165-miR-21-SOX6 axis in the wool follicle development. The targeting relationships of miR-21 with SOX6 and MSTRG.223165 were validated in dual-luciferase assays. This is the first report indicating the association of the lncRNA-miRNA-mRNA network with wool follicle development in AFWS. This study provides new insights into the regulation of the wool follicle growth and represents a solid foundation for wool sheep breeding programs.

The Genetic Regulation of Alternative Splicing in Populus deltoides

Alternative splicing (AS) is a mechanism of regulation of the proteome via enabling the production of multiple mRNAs from a single gene. To date, the dynamics of AS and its effects on the protein sequences of individuals in a large and genetically unrelated population of trees have not been investigated. Here we describe the diversity of AS events within a previously genotyped population of 268 individuals of Populus deltoides and their putative downstream functional effects. Using a robust bioinformatics pipeline, the AS events and resulting transcript isoforms were discovered and quantified for each individual in the population. Analysis of the AS revealed that, as expected, most AS isoforms are conserved. However, we also identified a substantial collection of new, unannotated splice junctions and transcript isoforms. Heritability estimates for the expression of transcript isoforms showed that approximately half of the isoforms are heritable. The genetic regulators of these AS isoforms and splice junction usage were then identified using a genome-wide association analysis. The expression of AS isoforms was predominately cis regulated while splice junction usage was generally regulated in trans. Additionally, we identified 696 genes encoding alternatively spliced isoforms that changed putative protein domains relative to the longest protein coding isoform of the gene, and 859 genes exhibiting this same phenomenon relative to the most highly expressed isoform. Finally, we found that 748 genes gained or lost micro-RNA binding sites relative to the longest protein coding isoform of a given gene, while 940 gained or lost micro-RNA binding sites relative to the most highly expressed isoform. These results indicate that a significant fraction of AS events are genetically regulated and that this isoform usage can result in protein domain architecture changes.

Choice of library size normalization and statistical methods for differential gene expression analysis in balanced two-group comparisons for RNA-seq studies

Background

High-throughput RNA sequencing (RNA-seq) has evolved as an important analytical tool in molecular biology. Although the utility and importance of this technique have grown, uncertainties regarding the proper analysis of RNA-seq data remain. Of primary concern, there is no consensus regarding which normalization and statistical methods are the most appropriate for analyzing this data. The lack of standardized analytical methods leads to uncertainties in data interpretation and study reproducibility, especially with studies reporting high false discovery rates. In this study, we compared a recently developed normalization method, UQ-pgQ2, with three of the most frequently used alternatives including RLE (relative log estimate), TMM (Trimmed-mean M values) and UQ (upper quartile normalization) in the analysis of RNA-seq data. We evaluated the performance of these methods for gene-level differential expression analysis by considering the factors, including: 1) normalization combined with the choice of a Wald test from DESeq2 and an exact test/QL (Quasi-likelihood) F-Test from edgeR; 2) sample sizes in two balanced two-group comparisons; and 3) sequencing read depths.

Results

Using the MAQC RNA-seq datasets with small sample replicates, we found that UQ-pgQ2 normalization combined with an exact test can achieve better performance in term of power and specificity in differential gene expression analysis. However, using an intra-group analysis of false positives from real and simulated data, we found that a Wald test performs better than an exact test when the number of sample replicates is large and that a QL F-test performs the best given sample sizes of 5, 10 and 15 for any normalization. The RLE, TMM and UQ methods performed similarly given a desired sample size.

Conclusion

We found the UQ-pgQ2 method combined with an exact test/QL F-test is the best choice in order to control false positives when the sample size is small. When the sample size is large, UQ-pgQ2 with a QL F-test is a better choice for the type I error control in an intra-group analysis. We observed read depths have a minimal impact for differential gene expression analysis based on the simulated data.

Two genes in a pathogenicity gene cluster encoding secreted proteins are required for appressorial penetration and infection of the maize anthracnose fungus Colletotrichum graminicola

To avoid pathogen-associated molecular pattern recognition, the hemibiotrophic maize pathogen Colletotrichum graminicola secretes proteins mediating the establishment of biotrophy. Targeted deletion of 26 individual candidate genes and seven gene clusters comprising 32 genes of C. graminicola identified a pathogenicity cluster (CLU5) of five co-linear genes, all of which, with the exception of CLU5b, encode secreted proteins. Targeted deletion of all genes of CLU5 revealed that CLU5a and CLU5d are required for full appressorial penetration competence, with virulence deficiencies independent of the host genotype and organ inoculated. Cytorrhysis experiments and microscopy showed that Δclu5a mutants form pressurized appressoria, but they are hampered in forming penetration pores and fail to differentiate a penetration peg. Whereas Δclu5d mutants elicited WT-like papillae, albeit at increased frequencies, papillae induced by Δclu5a mutants were much smaller than those elicited by the WT. Synteny of CLU5 is not only conserved in Colletotrichum spp. but also in additional species of Sordariomycetes including insect pathogens and saprophytes suggesting importance of CLU5 for fungal biology. Since CLU5a and CLU5d also occur in non-pathogenic fungi and since they are expressed prior to plant invasion and even in vegetative hyphae, the encoded proteins probably do not act primarily as effectors.

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.

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.

Transcriptional Analysis Allows Genome Reannotation and Reveals that Cryptococcus gattii VGII Undergoes Nutrient Restriction during Infection

Cryptococcus gattii is a human and animal pathogen that infects healthy hosts and caused the Pacific Northwest outbreak of cryptococcosis. The inhalation of infectious propagules can lead to internalization of cryptococcal cells by alveolar macrophages, a niche in which C. gattii cells can survive and proliferate. Although the nutrient composition of macrophages is relatively unknown, the high induction of amino acid transporter genes inside the phagosome indicates a preference for amino acid uptake instead of synthesis. However, the presence of countable errors in the R265 genome annotation indicates significant inhibition of transcriptomic analysis in this hypervirulent strain. Thus, we analyzed RNA-Seq data from in vivo and in vitro cultures of C. gattii R265 to perform the reannotation of the genome. In addition, based on in vivo transcriptomic data, we identified highly expressed genes and pathways of amino acid metabolism that would enable C. gattii to survive and proliferate in vivo. Importantly, we identified high expression in three APC amino acid transporters as well as the GABA permease. The use of amino acids as carbon and nitrogen sources, releasing ammonium and generating carbohydrate metabolism intermediaries, also explains the high expression of components of several degradative pathways, since glucose starvation is an important host defense mechanism.

Transcriptome analysis of smut fungi reveals widespread intergenic transcription and conserved antisense transcript expression

Background

Biotrophic fungal plant pathogens cause billions of dollars in losses to North American crops annually. The model for functional investigation of these fungi is Ustilago maydis. Its 20.5 Mb annotated genome sequence has been an excellent resource for investigating biotrophic plant pathogenesis. Expressed-sequence tag libraries and microarray hybridizations have provided insight regarding the type of transcripts produced by U. maydis but these analyses were not comprehensive and there were insufficient data for transcriptome comparison to other smut fungi. To improve transcriptome annotation and enable comparative analyses, comprehensive strand-specific RNA-seq was performed on cell-types of three related smut species: U. maydis (common smut of corn), Ustilago hordei (covered smut of barley), and Sporisorium reilianum (head smut of corn).

Results

In total, >1 billion paired-end sequence reads were obtained from haploid cell, dikaryon and teliospore RNA of U. maydis, haploid cell RNA of U. hordei, and haploid and dikaryon cell RNA of S. reilianum. The sequences were assembled into transfrags using Trinity, and updated gene models were created using PASA and categorized with Cufflinks Cuffcompare. Representative genes that were predicted for the first time with these RNA-seq analyses and genes with novel annotation features were independently assessed by reverse transcriptase PCR. The analyses indicate hundreds more predicted proteins, relative to the previous genome annotation, could be produced by U. maydis from altered transcript forms, and that the number of non-coding RNAs produced, including transcribed intergenic sequences and natural antisense transcripts, approximately equals the number of mRNAs. This high representation of non-coding RNAs appears to be a conserved feature of the smut fungi regardless of whether they have RNA interference machinery. Approximately 50% of the identified NATs were conserved among the smut fungi.

Conclusions

Overall, these analyses revealed: 1) smut genomes encode a number of transcriptional units that is twice the number of annotated protein-coding genes, 2) a small number of intergenic transcripts may encode proteins with characteristics of fungal effectors, 3) the vast majority of intergenic and antisense transcripts do not contain ORFs, 4) a large proportion of the identified antisense transcripts were detected at orthologous loci among the smut fungi, and 5) there is an enrichment of functional categories among orthologous loci that suggests antisense RNAs could have a genome-wide, non-RNAi-mediated, influence on gene expression in smut fungi.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3720-8) contains supplementary material, which is available to authorized users.

A comparison of per sample global scaling and per gene normalization methods for differential expression analysis of RNA-seq data

Normalization is an essential step with considerable impact on high-throughput RNA sequencing (RNA-seq) data analysis. Although there are numerous methods for read count normalization, it remains a challenge to choose an optimal method due to multiple factors contributing to read count variability that affects the overall sensitivity and specificity. In order to properly determine the most appropriate normalization methods, it is critical to compare the performance and shortcomings of a representative set of normalization routines based on different dataset characteristics. Therefore, we set out to evaluate the performance of the commonly used methods (DESeq, TMM-edgeR, FPKM-CuffDiff, TC, Med UQ and FQ) and two new methods we propose: Med-pgQ2 and UQ-pgQ2 (per-gene normalization after per-sample median or upper-quartile global scaling). Our per-gene normalization approach allows for comparisons between conditions based on similar count levels. Using the benchmark Microarray Quality Control Project (MAQC) and simulated datasets, we performed differential gene expression analysis to evaluate these methods. When evaluating MAQC2 with two replicates, we observed that Med-pgQ2 and UQ-pgQ2 achieved a slightly higher area under the Receiver Operating Characteristic Curve (AUC), a specificity rate > 85%, the detection power > 92% and an actual false discovery rate (FDR) under 0.06 given the nominal FDR (≤0.05). Although the top commonly used methods (DESeq and TMM-edgeR) yield a higher power (>93%) for MAQC2 data, they trade off with a reduced specificity (<70%) and a slightly higher actual FDR than our proposed methods. In addition, the results from an analysis based on the qualitative characteristics of sample distribution for MAQC2 and human breast cancer datasets show that only our gene-wise normalization methods corrected data skewed towards lower read counts. However, when we evaluated MAQC3 with less variation in five replicates, all methods performed similarly. Thus, our proposed Med-pgQ2 and UQ-pgQ2 methods perform slightly better for differential gene analysis of RNA-seq data skewed towards lowly expressed read counts with high variation by improving specificity while maintaining a good detection power with a control of the nominal FDR level.

Transcriptomics and Alternative Splicing Analyses Reveal Large Differences between Maize Lines B73 and Mo17 in Response to Aphid Rhopalosiphum padi Infestation

Maize (Zea mays L.) is a staple crop worldwide with extensive genetic variations. Various insects attack maize plants causing large yield loss. Here, we investigated the responses of maize B73, a susceptible line, and Mo17, a resistant line, to the aphid Rhopalosiphum padi on metabolite and transcriptome levels. R. padi feeding had no effect on the levels of the defensive metabolites benzoxazinoids (Bxs) in either line, and Mo17 contained substantially greater levels of Bxs than did B73. Profiling of the differentially expressed genes revealed that B73 and Mo17 responded to R. padi infestation specifically, and importantly, these two lines showed large gene expression differences even without R. padi herbivory. Correlation analysis identified four transcription factors (TFs) that might account for the high Bx levels in Mo17. Similarly, genome-wide alternative splicing (AS) analyses indicated that both B73 and Mo17 had temporally specific responses to R. padi infestation, and these two lines also exhibited large differences of AS regulation under normal condition, and 340 genes, including 10 TFs, were constantly differentially spliced. This study provides large-scale resource datasets for further studies on the mechanisms underlying maize-aphid interactions, and highlights the phenotypic divergence in defense against aphids among maize varieties.

Transcriptome Profile Reveals that Pu-Erh Tea Represses the Expression of Vitellogenin Family to Reduce Fat Accumulation in Caenorhabditis elegans

Due to misbalanced energy surplus and expenditure, obesity has become a common chronic disorder that is highly associated with many metabolic diseases. Pu-erh tea, a traditional Chinese beverage, has been believed to have numerous health benefits, such as anti-obesity. However, the underlying mechanisms of its anti-obesity effect are yet to be understood. Here, we take the advantages of transcriptional profile by RNA sequencing (RNA-Seq) to view the global gene expression of Pu-erh tea. The model organism Caenorhabditis elegans was treated with different concentrations of Pu-erh tea water extract (PTE, 0 g/mL, 0.025 g/mL, and 0.05 g/mL). Compared with the control, PTE indeed decreases lipid droplets size and fat accumulation. The high-throughput RNA-Sequence technique detected 18073 and 18105 genes expressed in 0.025 g/mL and 0.05 g/mL PTE treated groups, respectively. Interestingly, the expression of the vitellogenin family (vit-1, vit-2, vit-3, vit-4 and vit-5) was significantly decreased by PTE, which was validated by qPCR analysis. Furthermore, vit-1(ok2616), vit-3(ok2348) and vit-5(ok3239) mutants are insensitive to PTE triggered fat reduction. In conclusion, our transcriptional profile by RNA-Sequence suggests that Pu-erh tea lowers the fat accumulation primarily through repression of the expression of vit(vitellogenin) family, in addition to our previously reported (sterol regulatory element binding protein) SREBP-SCD (stearoyl-CoA desaturase) axis.

The Transient Receptor Potential (TRP) Channel Family in Colletotrichum graminicola: A Molecular and Physiological Analysis

Calcium (Ca2+) is a universal second messenger in all higher organisms and centrally involved in the launch of responses to environmental stimuli. Ca2+ signals in the cytosol are initiated by the activation of Ca2+ channels in the plasma membrane and/or in endomembranes. Yeast (Saccharomyces cerevisiae) contains a Ca2+-permeable channel of the TRP family, TRPY1, which is localized in the vacuolar membrane and contributes to cytosolic free Ca2+ ([Ca2+]cyt) elevations, for example in response to osmotic upshock. A TRPY1 homologue in the rice blast fungus is known to be important for growth and pathogenicity. To determine the role of the TRP channel family in the maize pathogen Colletotrichum graminicola, proteins homologous to TRPY1 were searched. This identified not one, but four genes in the C. graminicola genome, which had putative orthologs in other fungi, and which we named CgTRPF1 through 4. The topology of the CgTRPF proteins resembled that of TRPY1, albeit with a variable number of transmembrane (TM) domains additional to the six-TM-domain core and a diverse arrangement of putatively Ca2+-binding acidic motifs. All CgTRPF genes were expressed in axenic culture and throughout the infection of maize. Like TRPY1, all TRPF proteins of C. graminicola were localized intracellularly, albeit three of them were found not in large vacuoles, but co-localized in vesicular structures. Deletion strains for the CgTRPF genes were not altered in processes thought to involve Ca2+ release from internal stores, i.e. spore germination, the utilization of complex carbon sources, and the generation of tip-focussed [Ca2+]cyt spikes. Heterologous expression of CgTRPF1 through 4 in a tryp1Δ yeast mutant revealed that none of the channels mediated the release of Ca2+ in response to osmotic upshock. Accordingly, aequorin-based [Ca2+]cyt measurements of C. graminicola showed that in this fungus, osmotic upshock-triggered [Ca2+]cyt elevations were generated entirely by influx of Ca2+ from the extracellular space. Cgtrpf mutants did not show pathogenicity defects in leaf infection assays. In summary, our study reveals major differences between different fungi in the contribution of TRP channels to Ca2+-mediated signal transduction.

Structure-function characterization reveals new catalytic diversity in the galactose oxidase and glyoxal oxidase family

Alcohol oxidases, including carbohydrate oxidases, have a long history of research that has generated fundamental biological understanding and biotechnological applications. Despite a long history of study, the galactose 6-oxidase/glyoxal oxidase family of mononuclear copper-radical oxidases, Auxiliary Activity Family 5 (AA5), is currently represented by only very few characterized members. Here we report the recombinant production and detailed structure-function analyses of two homologues from the phytopathogenic fungi Colletotrichum graminicola and C. gloeosporioides, CgrAlcOx and CglAlcOx, respectively, to explore the wider biocatalytic potential in AA5. EPR spectroscopy and crystallographic analysis confirm a common active-site structure vis-à-vis the archetypal galactose 6-oxidase from Fusarium graminearum. Strikingly, however, CgrAlcOx and CglAlcOx are essentially incapable of oxidizing galactose and galactosides, but instead efficiently catalyse the oxidation of diverse aliphatic alcohols. The results highlight the significant potential of prospecting the evolutionary diversity of AA5 to reveal novel enzyme specificities, thereby informing both biology and applications.