Genome-wide identification of long non-coding RNAs suggests a potential association with effector gene transcription in Phytophthora sojae

Underground communication: Long non-coding RNA signaling in the plant rhizosphere

Long non-coding RNAs (lncRNAs) have emerged as integral gene-expression regulators underlying plant growth, development, and adaptation. To adapt to the heterogeneous and dynamic rhizosphere, plants use interconnected regulatory mechanisms to optimally fine-tune gene-expression-governing interactions with soil biota, as well as nutrient acquisition and heavy metal tolerance. Recently, high-throughput sequencing has enabled the identification of plant lncRNAs responsive to rhizosphere biotic and abiotic cues. Here, we examine lncRNA biogenesis, classification, and mode of action, highlighting the functions of lncRNAs in mediating plant adaptation to diverse rhizosphere factors. We then discuss studies that reveal the significance and target genes of lncRNAs during developmental plasticity and stress responses at the rhizobium interface. A comprehensive understanding of specific lncRNAs, their regulatory targets, and the intricacies of their functional interaction networks will provide crucial insights into how these transcriptomic switches fine-tune responses to shifting rhizosphere signals. Looking ahead, we foresee that single-cell dissection of cell-type-specific lncRNA regulatory dynamics will enhance our understanding of the precise developmental modulation mechanisms that enable plant rhizosphere adaptation. Overcoming future challenges through multi-omics and genetic approaches will more fully reveal the integral roles of lncRNAs in governing plant adaptation to the belowground environment.

Long noncoding RNAs emerge from transposon-derived antisense sequences and may contribute to infection stage-specific transposon regulation in a fungal phytopathogen

BACKGROUND

The genome of the obligate biotrophic phytopathogenic barley powdery mildew fungus Blumeria hordei is inflated due to highly abundant and possibly active transposable elements (TEs). In the absence of the otherwise common repeat-induced point mutation transposon defense mechanism, noncoding RNAs could be key for regulating the activity of TEs and coding genes during the pathogenic life cycle.

RESULTS

We performed time-course whole-transcriptome shotgun sequencing (RNA-seq) of total RNA derived from infected barley leaf epidermis at various stages of fungal pathogenesis and observed significant transcript accumulation and time point-dependent regulation of TEs in B. hordei. Using a manually curated consensus database of 344 TEs, we discovered phased small RNAs mapping to 104 consensus transposons, suggesting that RNA interference contributes significantly to their regulation. Further, we identified 5,127 long noncoding RNAs (lncRNAs) genome-wide in B. hordei, of which 823 originated from the antisense strand of a TE. Co-expression network analysis of lncRNAs, TEs, and coding genes throughout the asexual life cycle of B. hordei points at extensive positive and negative co-regulation of lncRNAs, subsets of TEs and coding genes.

CONCLUSIONS

Our work suggests that similar to mammals and plants, fungal lncRNAs support the dynamic modulation of transcript levels, including TEs, during pivotal stages of host infection. The lncRNAs may support transcriptional diversity and plasticity amid loss of coding genes in powdery mildew fungi and may give rise to novel regulatory elements and virulence peptides, thus representing key drivers of rapid evolutionary adaptation to promote pathogenicity and overcome host defense.

Genome-wide identification and characterization of Puccinia striiformis-responsive lncRNAs in Triticum aestivum

Wheat stripe rust (yellow rust) caused by Puccinia striiformis f. sp. tritici (Pst) is a serious biotic stress factor limiting wheat production worldwide. Emerging evidence demonstrates that long non-coding RNAs (lncRNAs) participate in various developmental processes in plants via post-transcription regulation. In this study, RNA sequencing (RNA-seq) was performed on a pair of near-isogenic lines-rust resistance line FLW29 and rust susceptible line PBW343-which differed only in the rust susceptibility trait. A total of 6,807 lncRNA transcripts were identified using bioinformatics analyses, among which 10 lncRNAs were found to be differentially expressed between resistance and susceptible lines. In order to find the target genes of the identified lncRNAs, their interactions with wheat microRNA (miRNAs) were predicted. A total of 199 lncRNAs showed interactions with 65 miRNAs, which further target 757 distinct mRNA transcripts. Moreover, detailed functional annotations of the target genes were used to identify the candidate genes, pathways, domains, families, and transcription factors that may be related to stripe rust resistance response in wheat plants. The NAC domain protein, disease resistance proteins RPP13 and RPM1, At1g58400, monodehydroascorbate reductase, NBS-LRR-like protein, rust resistance kinase Lr10-like, LRR receptor, serine/threonine-protein kinase, and cysteine proteinase are among the identified targets that are crucial for wheat stripe rust resistance. Semiquantitative PCR analysis of some of the differentially expressed lncRNAs revealed variations in expression profiles of two lncRNAs between the Pst-resistant and Pst-susceptible genotypes at least under one condition. Additionally, simple sequence repeats (SSRs) were also identified from wheat lncRNA sequences, which may be very useful for conducting targeted gene mapping studies of stripe rust resistance in wheat. These findings improved our understanding of the molecular mechanism responsible for the stripe rust disease that can be further utilized to develop wheat varieties with durable resistance to this disease.

Study on the Applications and Regulatory Mechanisms of Grafting on Vegetables

Grafting can overcome problems with soil sensitivity, enhance plant stress tolerance, improve product quality, and increase crop yield and value. This paper reviews the various mechanisms of vegetable grafting, the graft survival process and its influencing factors, the practical applications of grafting, and the molecular regulation of grafting in vegetables. The importance of germplasm and rootstock interactions, the mechanization of vegetable grafting, and future aspects, including intelligence and digitalization, are discussed.

Identification of Long Intergenic Noncoding RNAs in Rhizoctonia cerealis following Inoculation of Wheat

Wheat sharp eyespot caused by Rhizoctonia cerealis is primarily a severe threat to worldwide wheat production. Currently, there are no resistant wheat cultivars, and the use of fungicides is the primary method for controlling this disease. Elucidating the mechanisms of R. cerealis pathogenicity can accelerate the pace of the control of this disease. Long intergenic noncoding RNAs (lincRNAs) that function in plant-pathogen interactions might provide a new perspective. We systematically analyzed lincRNAs and identified a total of 1,319 lincRNAs in R. cerealis. We found that lincRNAs are involved in various biological processes, as shown by differential expression analysis and weighted correlation network analysis (WGCNA). Next, one of nine hub lincRNAs in the blue module that was related to infection and growth processes, MSTRG.4380.1, was verified to reduce R. cerealis virulence on wheat by a host-induced gene silencing (HIGS) assay. Following that, RNA sequencing (RNA-Seq) analysis revealed that the significantly downregulated genes in the MSTRG.4380.1 knockdown lines were associated mainly with infection-related processes, including hydrolase, transmembrane transporter, and energy metabolism activities. Additionally, 23 novel microRNAs (miRNAs) were discovered during small RNA (sRNA) sequencing (sRNA-Seq) analysis of MSTRG.4380.1 knockdown, and target prediction of miRNAs suggested that MSTRG.4380.1 does not act as a competitive endogenous RNA (ceRNA). This study performed the first genome-wide identification of R. cerealis lincRNAs and miRNAs. It confirmed the involvement of a lincRNA in the infection process, providing new insights into the mechanism of R. cerealis infection and offering a new approach for protecting wheat from R. cerealis. IMPORTANCE Rhizoctonia cerealis, the primary causal agent of wheat sharp eyespot, has caused significant losses in worldwide wheat production. Since no resistant wheat cultivars exist, chemical control is the primary method. However, this approach is environmentally unfriendly and costly. RNA interference (RNAi)-mediated pathogenicity gene silencing has been proven to reduce the growth of Rhizoctonia and provides a new perspective for disease control. Recent studies have shown that lincRNAs are involved in various biological processes across species, such as biotic and abiotic stresses. Therefore, verifying the function of lincRNAs in R. cerealis is beneficial for understanding the infection mechanism. In this study, we reveal that lincRNAs could contribute to the virulence of R. cerealis, which provides new insights into controlling this pathogen.

Roles of long non-coding RNAs in plant immunity

Robust plant immune systems are fine-tuned by both protein-coding genes and non-coding RNAs. Long non-coding RNAs (lncRNAs) refer to RNAs with a length of more than 200 nt and usually do not have protein-coding function and do not belong to any other well-known non-coding RNA types. The non-protein-coding, low expression, and non-conservative characteristics of lncRNAs restrict their recognition. Although studies of lncRNAs in plants are in the early stage, emerging studies have shown that plants employ lncRNAs to regulate plant immunity. Moreover, in response to stresses, numerous lncRNAs are differentially expressed, which manifests the actions of low-expressed lncRNAs and makes plant-microbe/insect interactions a convenient system to study the functions of lncRNAs. Here, we summarize the current advances in plant lncRNAs, discuss their regulatory effects in different stages of plant immunity, and highlight their roles in diverse plant-microbe/insect interactions. These insights will not only strengthen our understanding of the roles and actions of lncRNAs in plant-microbe/insect interactions but also provide novel insight into plant immune responses and a basis for further research in this field.

The SET domain protein PsKMT3 regulates histone H3K36 trimethylation and modulates effector gene expression in the soybean pathogen Phytophthora sojae

Plant pathogens secrete effector proteins to overcome host immunity and promote colonization. In oomycete plant pathogens, the expression of many effector genes is altered upon infection; however, the regulatory mechanisms are unclear. In this study, we identified a su(var)3-9, enhancer of zeste, and trithorax (SET) domain protein-encoding gene, PsKMT3, that was highly induced at early infection stages in Phytophthora sojae. Deletion of PsKMT3 led to asexual development and pathogenicity defects. Chromatin immunoprecipitation followed by sequencing (ChIP-seq) and western blot analyses demonstrated that histone H3K36 trimethylation (H3K36me3) was significantly reduced genome-wide in mutants. RNA-seq analysis identified 374 genes encoding secreted proteins that were differentially expressed in pskmt3 at the mycelium stage. The significantly altered genes encompassed the RxLR (Arg-x-Lys-Arg) effector gene family, including the essential effector genes Avh23, Avh181, Avh240, and Avh241. Transcriptome analysis at early infection stages showed misregulation of effector gene expression waves in pskmt3. H3K36me3 was directly and indirectly associated with RxLR effector gene activation. Our results reveal a role of a SET domain protein in regulating effector gene expression and modulating histone methylation in P. sojae.

Systematic Characterization and Regulatory Role of lncRNAs in Asian Honey Bees Responding to Microsporidian Infestation

Long noncoding RNAs (lncRNAs) are pivotal regulators in gene expression and diverse biological processes, such as immune defense and host-pathogen interactions. However, little is known about the roles of lncRNAs in the response of the Asian honey bee (Apis cerana) to microsporidian infestation. Based on our previously obtained high-quality transcriptome datasets from the midgut tissues of Apis cerana cerana workers at 7 days post inoculation (dpi) and 10 dpi with Nosema ceranae (AcT7 and AcT10 groups) and the corresponding un-inoculated midgut tissues (AcCK7 and AcCK10 groups), the transcriptome-wide identification and structural characterization of lncRNAs were conducted, and the differential expression pattern of lncRNAs was then analyzed, followed by investigation of the regulatory roles of differentially expressed lncRNAs (DElncRNAs) in host response. Here, 2365, 2322, 2487, and 1986 lncRNAs were, respectively, identified in the AcCK7, AcT7, AcCK7, and AcT10 groups. After removing redundant ones, a total of 3496 A. c. cerana lncRNAs were identified, which shared similar structural characteristics with those discovered in other animals and plants, such as shorter exons and introns than mRNAs. Additionally, 79 and 73 DElncRNAs were screened from the workers' midguts at 7 dpi and 10 dpi, respectively, indicating the alteration of the overall expression pattern of lncRNAs in host midguts after N. ceranae infestation. These DElncRNAs could, respectively, regulate 87 and 73 upstream and downstream genes, involving a suite of functional terms and pathways, such as metabolic process and Hippo signaling pathway. Additionally, 235 and 209 genes co-expressed with DElncRNAs were found to enrich in 29 and 27 terms, as well as 112 and 123 pathways, such as ABC transporters and the cAMP signaling pathway. Further, it was detected that 79 (73) DElncRNAs in the host midguts at 7 (10) dpi could target 321 (313) DEmiRNAs and further target 3631 (3130) DEmRNAs. TCONS_00024312 and XR_001765805.1 were potential precursors for ame-miR-315 and ame-miR-927, while TCONS_00006120 was the putative precursor for both ame-miR-87-1 and ame-miR-87-2. These results together suggested that DElncRNAs are likely to play regulatory roles in the host response to N. ceranae infestation through the regulation of neighboring genes via a cis-acting effect, modulation of co-expressed mRNAs via trans-acting effect, and control of downstream target genes' expression via competing endogenous RNA networks. Our findings provide a basis for disclosing the mechanism underlying DElncRNA-mediated host N. ceranae response and a new perspective into the interaction between A. c. cerana and N. ceranae.

An Oomycete-Specific Leucine-Rich Repeat-Containing Protein Is Involved in Zoospore Flagellum Development in Phytophthora sojae

A leucine-rich repeat (LRR) is a widespread structural motif of 20 to 30 amino acids with characteristic repetitive sequences rich in leucine. LRR-containing proteins are critical for ligand recognition and binding, participating in plant development and defense. Like plants, oomycetes also harbor genes encoding LRR-containing proteins, but their functions remain largely unknown. We identified a zoospore-upregulated gene from Phytophthora sojae with LRRs and an extra structural maintenance of chromosomes-like domain. We generated knockout and complemented knockout strains of this LRR protein and found that its deletion resulted in a pronounced reduction in zoospore mobility and chemotaxis, cyst germination, and virulence. Interestingly, micro-examination of zoospores under a scanning electron microscope revealed irregularly shaped zoospores without flagella in these deletion mutants. In addition, the reintroduction of this LRR protein into the knockout mutant reversed all the deficiencies. Our data demonstrate a critical role for the Phytophthora LRR protein in modulating zoospore development, which impairs migration to the host soybean and affects the spread of Phytophthora pathogens.

Multi-Omics Techniques for Soybean Molecular Breeding

Soybean is a major crop that provides essential protein and oil for food and feed. Since its origin in China over 5000 years ago, soybean has spread throughout the world, becoming the second most important vegetable oil crop and the primary source of plant protein for global consumption. From early domestication and artificial selection through hybridization and ultimately molecular breeding, the history of soybean breeding parallels major advances in plant science throughout the centuries. Now, rapid progress in plant omics is ushering in a new era of precision design breeding, exemplified by the engineering of elite soybean varieties with specific oil compositions to meet various end-use targets. The assembly of soybean reference genomes, made possible by the development of genome sequencing technology and bioinformatics over the past 20 years, was a great step forward in soybean research. It facilitated advances in soybean transcriptomics, proteomics, metabolomics, and phenomics, all of which paved the way for an integrated approach to molecular breeding in soybean. In this review, we summarize the latest progress in omics research, highlight novel findings made possible by omics techniques, note current drawbacks and areas for further research, and suggest that an efficient multi-omics approach may accelerate soybean breeding in the future. This review will be of interest not only to soybean breeders but also to researchers interested in the use of cutting-edge omics technologies for crop research and improvement.

ATAC-Seq Reveals the Landscape of Open Chromatin and cis-Regulatory Elements in the Phytophthora sojae Genome

Nucleosome-free open chromatin often harbors transcription factor (TF)-binding sites that are associated with active cis-regulatory elements. However, analysis of open chromatin regions has rarely been applied to oomycete or fungal plant pathogens. In this study, we performed the assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq) to identify open chromatin and cis-regulatory elements in Phytophthora sojae at the mycelial stage. We identified 10,389 peaks representing nucleosome-free regions (NFRs). The peaks were enriched in gene-promoter regions and associated with 40% of P. sojae genes; transcription levels were higher for genes with multiple peaks than genes with a single peak and were higher for genes with a single peak than genes without peak. Chromatin accessibility was positively correlated with gene transcription level. Through motif discovery based on NFR peaks in core promoter regions, 25 candidate cis-regulatory motifs with evidence of TF-binding footprints were identified. These motifs exhibited various preferences for location in the promoter region and associations with the transcription level of their target genes, which included some putative pathogenicity-related genes. As the first study revealing the landscape of open chromatin and the correlation between chromatin accessibility and gene transcription level in oomycetes, the results provide a technical reference and data resources for future studies on the regulatory mechanisms of gene transcription.[Formula: see text] Copyright © 2022 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Function identification of miR394 in tomato resistance to Phytophthora infestans

MiR394 plays a negative role in tomato resistance to late blight. The lncRNA40787 severing as an eTM for miR394 to regulate LCR and exerting functions in tomato resistance. Tomato (Solanum lycopersicum), which was used as model species for studying the mechanism of plant disease defense, is susceptible to multiple pathogens. Non-coding RNA (ncRNA) has a pivotal role in plants response to biological stresses. It has previously been observed that the expression level of miR394 changed significantly after the infection of various pathogens. However, there has been no detailed investigation of the accumulated or suppressed mechanism of miR394. Our previous study predicted three lncRNAs (lncRNA40787, lncRNA27177, and lncRNA42566) that contain miR394 endogenous target mimics (eTM), which may exist as the competitive endogenous RNAs (ceRNAs) of miR394. In our study, the transcription levels of these three lncRNAs were strongly up-regulated in tomato upon infection with P. infestans. In contrast with the three lncRNAs, the accumulation of miR394 was significantly suppressed. Based on the expression pattern, and value of minimum free energy (mfes) that represents the binding ability between lncRNA and miRNA, lncRNA40787 was chosen for further investigation. Results showed that overexpression of lncRNA40787 reduced the expression of miR394 along with decreased lesion area and enhanced disease resistance. Overexpression of miR394, however, decreased the expression of its target gene Leaf Curling Responsiveness (LCR), and suppressed the synthesis components genes of jasmonic acid (JA), depressing the resistance of tomato to P. infestans infection. Taken together, our findings indicated that miR394 can be decoyed by lncRNA40787, and negatively regulated the expression of LCR to enhance tomato susceptibility under P. infestans infection. Our study provided detailed information on the lncRNA40787-miR394-LCR regulatory network and serves as a reference for future research.

Transcriptomic profiling identifies novel transcripts, isomorphs, and noncoding RNAs in Paracoccidioides brasiliensis

This paper describes a transcriptomic profiling of Paracoccidioides brasiliensis (Pb) performed with the aid of an RNA-seq-based approach, aimed at characterizing the general transcriptome in this human pathogenic fungus, responsible for paracoccidioidomycosis (PCM). Results confirm that ∼75% of the genes currently annotated in the P. brasiliensis genome are, in fact, transcribed in vivo and that ∼19% of them may display alternative isomorphs. Moreover, we identified 627 transcripts that do not match any gene currently mapped in the genome, represented by 114 coding transcripts (probably derived from previously unmapped protein-coding genes) and 513 noncoding RNAs (ncRNAs), including 203 long-noncoding RNAs (lncRNAs).

Evolution and Regulation of a Large Effector Family of Pyricularia oryzae

Plant pathogen effectors play important roles in parasitism, including countering plant immunity. However, investigations of the emergence and diversification of fungal effectors across host-adapted populations has been limited. We previously identified a gene encoding a suppressor of plant cell death in Pyricularia oryzae (syn. Magnaporthe oryzae). Here, we report the gene is one of a 21-member gene family and we characterize sequence diversity in different populations. Within the rice pathogen population, nucleotide diversity is low, however; the majority of gene family members display presence-absence polymorphism or other null alleles. Gene family allelic diversity is greater between host-adapted populations and, thus, we named them host-adapted genes (HAGs). Multiple copies of HAGs were found in some genome assemblies and sequence divergence between the alleles in two cases suggested they were the result of repeat-induced point mutagenesis. Transfer of family members between populations and novel HAG haplotypes resulting from apparent recombination were observed. HAG family transcripts were induced in planta and a subset of HAGs are dependent on a key regulator of pathogenesis, PMK1. We also found differential intron splicing for some HAGs that would prevent ex planta protein expression. For some genes, spliced transcript was expressed in antiphase with an overlapping antisense transcript. Characterization of HAG expression patterns and allelic diversity reveal novel mechanisms for HAG regulation and mechanisms generating sequence diversity and novel allele combinations. This evidence of strong in planta-specific expression and selection operating on the HAG family is suggestive of a role in parasitism.[Formula: see text] Copyright © 2021 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Secretory Peptides as Bullets: Effector Peptides from Pathogens against Antimicrobial Peptides from Soybean

Soybean is an important crop as both human food and animal feed. However, the yield of soybean is heavily impacted by biotic stresses including insect attack and pathogen infection. Insect bites usually make the plants vulnerable to pathogen infection, which causes diseases. Fungi, oomycetes, bacteria, viruses, and nematodes are major soybean pathogens. The infection by pathogens and the defenses mounted by soybean are an interactive and dynamic process. Using fungi, oomycetes, and bacteria as examples, we will discuss the recognition of pathogens by soybean at the molecular level. In this review, we will discuss both the secretory peptides for soybean plant infection and those for pathogen inhibition. Pathogenic secretory peptides and peptides secreted by soybean and its associated microbes will be included. We will also explore the possible use of externally applied antimicrobial peptides identical to those secreted by soybean and its associated microbes as biopesticides.

Sl-lncRNA15492 interacts with Sl-miR482a and affects Solanum lycopersicum immunity against Phytophthora infestans

Long non-coding RNAs (lncRNAs) are involved in the resistance of plants to infection by pathogens via interactions with microRNAs (miRNAs). Long non-coding RNAs are cleaved by miRNAs to produce phased small interfering RNAs (phasiRNAs), which, as competing endogenous RNAs (ceRNAs), function as decoys for mature miRNAs, thus inhibiting their expression, and contain pre-miRNA sequences to produce mature miRNAs. However, whether lncRNAs and miRNAs mediate other molecular mechanisms during plant resistance to pathogens is unknown. In this study, as a positive regulator, Sl-lncRNA15492 from tomato (Solanum lycopersicum Zaofen No. 2) plants affected tomato resistance to Phytophthora infestans. Gain- and loss-of-function experiments and RNA ligase-mediated 5'-amplification of cDNA ends (RLM-5' RACE) also revealed that Sl-miR482a was negatively involved in tomato resistance by targeting Sl-NBS-LRR genes and that silencing of Sl-NBS-LRR1 decreased tomato resistance. Sl-lncRNA15492 inhibited the expression of mature Sl-miR482a, whose precursor was located within the antisense sequence of Sl-lncRNA15492. Further degradome analysis and additional RLM-5' RACE experiments verified that mature Sl-miR482a could also cleave Sl-lncRNA15492. These results provide a mechanism by which lncRNAs might inhibit precursor miRNA expression through antisense strands of lncRNAs, and demonstrate that Sl-lncRNA15492 and Sl-miR482a mutually inhibit the maintenance of Sl-NBS-LRR1 homeostasis during tomato resistance to P. infestans.

Conserved Subgroups of the Plant-Specific RWP-RK Transcription Factor Family Are Present in Oomycete Pathogens

Nitrogen is a major constituent of proteins, chlorophyll, nucleotides, and hormones and has profound effects on plant growth and productivity. RWP-RK family transcription factors (TFs) are key regulators that bind to cis-acting elements in the promoter regions of nitrogen use efficiency-related genes and genes responsible for gametogenesis and embryogenesis. The proteins share a conserved RWPxRK motif; have been found in all vascular plants, green algae, and slime molds; and are considered to be a plant-specific TF family. In this study, we show that RWP-RK proteins are also widely present in the Stramenopila kingdom, particularly among the oomycetes, with 12-15 members per species. These proteins form three distinct phylogenetic subgroups, two of which are relatively closely related to the nodule inception (NIN)-like protein (NLP) or the RWP-RK domain protein (RKD) subfamilies of plant RWP-RK proteins. The donor for horizontal gene transfer of RWP-RK domains to slime molds is likely to have been among the Stramenopila, predating the divide between brown algae and oomycetes. The RWP-RK domain has secondary structures that are conserved across plants and oomycetes, but several amino acids that may affect DNA-binding affinity differ. The transcriptional activities of orthologous RWP-RK genes were found to be conserved in oomycetes. Our results demonstrate that RWP-RK family TF genes are present in the oomycetes and form specific subgroups with functions that are likely conserved. Our results provide new insights for further understanding the evolution and function of this TF family in specific eukaryotic organisms.

XRN1-associated long non-coding RNAs may contribute to fungal virulence and sexual development in entomopathogenic fungus Cordyceps militaris

BACKGROUND

Numerous long non-coding RNAs (lncRNAs) identified and characterized in mammals, plants, and fungi have been found to play critical regulatory roles in biological processes. However, little is known about the role of lncRNAs in insect pathogenic fungi.

RESULTS

By profiling the transcriptomes of sexual and asexual development in the insect-pathogenic fungus Cordyceps militaris, 4140 lncRNAs were identified and found to be dynamically expressed during fungal development. The lncRNAs had shorter transcript lengths and lower numbers of exons compared to protein-coding genes. The expressed target genes (neighboring and cis-regulated) of various expressed lncRNAs were predicted, and these genes showed significant enrichment in energy metabolism and signaling pathways, such as 'Glycolysis/Gluconeogenesis' and "MAPK signaling pathway". To better understand how lncRNAs function in the fungus, xrn1, the final gene of the NMD pathway, which determines the fate of lncRNAs, was disrupted. The Δxrn1 deletion mutant displayed significant (P < 0.05) attenuation of virulence and a lower growth rate in C. militaris. Quantitative RT-PCR results revealed 10 lncRNAs with significantly higher expression, while 8 of these 10 lncRNA target genes (virulence- and sexual development-related) showed significantly lower expression in Δxrn1 compared to in the wild-type, suggesting that lncRNA expression regulates fungal virulence and sexual development by affecting gene expression.

CONCLUSION

These findings suggest that lncRNAs in C. militaris play important roles in the fungal infection progress and fruiting body production, providing a broad repertoire and resource for further studies of lncRNAs. © 2019 Society of Chemical Industry.

Whole Genome Re-sequencing Reveals Natural Variation and Adaptive Evolution of Phytophthora sojae

Due to the monocultural basis of agricultural crops, mutated plant microbes with increased pathogenicity can easily spread in the field and lead to serious yield losses. As a major threat to a wide range of crop plants, oomycete pathogens continuously undergo adaptive evolution to overcome plant defense barriers. However, the genetic basis of their evolution at the molecular level remains largely unknown. Here, we investigated the nature variation and the population genomics of the soybean pathogen Phytophthora sojae by high-throughput genome re-sequencing. Genomic variation analysis revealed uneven “two-speed” evolutionary pattern with genes in gene-sparse regions (GSRs) showing higher rates of structural polymorphisms and positive selection. GSRs are enriched in effector genes and transposase-related genes. Our results also suggested that the NADH oxidase and MIP transporter gene families undergo rapid and diversifying selection. Furthermore, we demonstrated that P. sojae isolates possess varying numbers of RxLR effectors with diverse sequences, totaling 471 members. Among them, 42 core RxLR effectors are assumed to be important for infection. Finally, we observed that Avr genes exhibit abundant sequence variation in P. sojae isolates. Several novel variants lead to the evading of host resistance, including a complete deletion in Avr3c and amino acid mutations in Avr1a. Taken together, our results provide an adaptive landscape of P. sojae at single-nucleotide resolution, as well as resources for further resistance breeding and disease prevention against this important plant pathogen.

Genome-Wide Identification and Functional Prediction of Long Non-coding RNAs Involved in the Heat Stress Response in Metarhizium robertsii

Long non-coding RNAs (lncRNAs) play a significant role in stress responses. To date, only a few studies have reported the role of lncRNAs in insect-pathogenic fungi. Here, we report a genome-wide transcriptional analysis of lncRNAs produced in response to heat stress in Metarhizium robertsii, a model insect-pathogenic fungus, using strand-specific RNA sequencing. A total of 1655 lncRNAs with 1742 isoforms were identified, of which 1081 differentially expressed (DE) lncRNAs were characterized as being heat responsive. By characterizing their genomic structures and expression patterns, we found that the lncRNAs possessed shorter transcripts, fewer exons, and lower expression levels than the protein-coding genes in M. robertsii. Furthermore, target prediction analysis of the lncRNAs revealed thousands of potential DE lncRNA–messenger RNA (mRNA) pairs, among which 5381 pairs function in the cis-regulatory mode. Further pathway enrichment analysis of the corresponding cis-regulated target genes showed that the targets were significantly enriched in the following biological pathways: the Hippo signaling pathway and cell cycle. This finding suggested that these DE lncRNAs control the expression of their corresponding neighboring genes primarily through environmental information processing and cellular processes. Moreover, only 26 trans-regulated lncRNA–mRNA pairs were determined. In addition, among the targets of heat-responsive lncRNAs, two classic genes that may be involved in the response to heat stress were also identified, including hsp70 (XM_007821830 and XM_007825705). These findings expand our knowledge of lncRNAs as important regulators of the response to heat stress in filamentous fungi, including M. robertsii.

Distribution, Characteristics, and Regulatory Potential of Long Noncoding RNAs in Brown-Rot Fungi

Long noncoding RNAs have been thoroughly studied in plants, animals, and yeasts, where they play important roles as regulators of transcription. Nevertheless, almost nothing is known about their presence and characteristics in filamentous fungi, especially in basidiomycetes. In the present study, we have carried out an exhaustive annotation and characterization of lncRNAs in two lignin degrader basidiomycetes, Coniophora puteana and Serpula lacrymans. We identified 2,712 putative lncRNAs in the former and 2,242 in the latter, mainly originating from intergenic locations of transposon-sparse genomic regions. The lncRNA length, GC content, expression levels, and stability of the secondary structure differ from coding transcripts but are similar in these two species and resemble that of other eukaryotes. Nevertheless, they lack sequence conservation. Also, we found that lncRNAs are transcriptionally regulated in the same proportion as genes when the fungus actively decomposes soil organic matter. Finally, up to 7% of the upstream gene regions of Coniophora puteana and Serpula lacrymans are transcribed and produce lncRNAs. The study of expression trends in these gene-lncRNA pairs uncovered groups with similar and opposite transcriptional profiles which may be the result of cis-transcriptional regulation.