Fungal pathogens can have devastating effects on global crop production, leading to annual economic losses ranging from 10% to 23%. In light of climate change-related challenges, researchers anticipate an increase in fungal infections as a result of shifting environmental conditions. However, plants have developed intricate molecular mechanisms for effective defense against fungal attacks. Understanding these mechanisms is essential to the development of new strategies for protecting crops from multiple fungi threats. Public omics databases provide valuable resources for research on plant–pathogen interactions; however, integrating data from different studies can be challenging due to experimental variation. In this study, we aimed to identify the core genes that defend against the pathogenic fungi Colletotrichum higginsianum and Botrytis cinerea in Arabidopsis thaliana. Using a custom framework to control batch effects and construct Gene Co-expression Networks in publicly available RNA-seq dataset from infected A. thaliana plants, we successfully identified a gene module that was responsive to both pathogens. We also performed gene annotation to reveal the roles of previously unknown protein-coding genes in plant defenses against fungal infections. This research demonstrates the potential of publicly available RNA-seq data for identifying the core genes involved in defending against multiple fungal pathogens.

BACKGROUND:Pentatricopeptide repeat (PPR) proteins play an essential role in the post-transcriptional regulation of genes in plastid genomes. Although important advances have been made in understanding the functions of these genes, there is little information available on chloroplastic PPR genes in non-model plants and less in plants without chloroplasts. In the present study, a comprehensive and multifactorial bioinformatic strategy was applied to search for putative PPR genes in the foliar and meristematic tissues of green and albino plantlets of the non-model plant Agave angustifolia Haw.RESULTS:A total of 1581 PPR transcripts were identified, of which 282 were chloroplastic. Leaf tissue in the albino plantlets showed the highest levels of expression of chloroplastic PPRs. The search for hypothetical targets of 12 PPR sequences in the chloroplast genes of A. angustifolia revealed their action on transcripts related to ribosomes and translation, photosystems, ATP synthase, plastid-encoded RNA polymerase and RuBisCO.CONCLUSIONS:Our results suggest that the expression of PPR genes depends on the state of cell differentiation and plastid development. In the case of the albino leaf tissue, which lacks functional chloroplasts, it is possible that anterograde and retrograde signaling networks are severely compromised, leading to a compensatory anterograde response characterized by an increase in the expression of PPR genes.

Banana is the most popular fruit in the world, with a relevant role in food security for more than 400 million people. However, fungal diseases cause substantial losses every year. A better understanding of the banana immune system should facilitate the development of new disease-resistant cultivars. In this study, we performed a genome-wide analysis of the leucine-rich repeat receptor-like protein (LRR-RLP) disease resistance gene family in a wild banana. We identified 78 LRR-RLP genes in the banana genome. Remarkably, seven MaLRR-RLPs formed a gene cluster in the distal part of chromosome 10, where resistance to Fusarium wilt caused by Foc race 1 has been previously mapped. Hence, we proposed these seven MaLRR-RLPs as resistance gene candidates (RGCs) for Fusarium wilt. We also identified seven other banana RGCs based on their close phylogenetic relationships with known LRR-RLP proteins. Moreover, phylogenetic analysis of the banana, rice, and Arabidopsis LRR-RLP families revealed five major phylogenetic clades shared by these plant species. Finally, transcriptomic analysis of the MaLRR-RLP gene family in plants treated with Foc race 1 or Foc TR4 showed the expression of several members of this family, and some of them were upregulated in response to these Foc races. Our study provides novel insights into the structure, distribution, evolution, and expression of the LRR-RLP gene family in bananas as well as valuable RGCs that will facilitate the identification of disease resistance genes for the genetic improvement of this crop.

Resilience of growing in arid and semiarid regions and a high capacity of accumulating sugar-rich biomass with low lignin percentages have placed Agave species as an emerging bioenergy crop. Although transcriptome sequencing of fiber-producing agave species has been explored, molecular bases that control wall cell biogenesis and metabolism in agave species are still poorly understood. Here, through RNAseq data mining, we reconstructed the cellulose biosynthesis pathway and the phenylpropanoid route producing lignin monomers in A. tequilana, and evaluated their expression patterns in silico and experimentally. Most of the orthologs retrieved showed differential expression levels when they were analyzed in different tissues with contrasting cellulose and lignin accumulation. Phylogenetic and structural motif analyses of putative CESA and CAD proteins allowed to identify those potentially involved with secondary cell wall formation. RT-qPCR assays revealed enhanced expression levels of AtqCAD5 and AtqCESA7 in parenchyma cells associated with extraxylary fibers, suggesting a mechanism of formation of sclerenchyma fibers in Agave similar to that reported for xylem cells in model eudicots. Overall, our results provide a framework for understanding molecular bases underlying cell wall biogenesis in Agave species studying mechanisms involving in leaf fiber development in monocots.

Together with their undeniable role in the ecology of arid and semiarid ecosystems, Agave species are emerging as a model to dissect the relationships between crassulacean acid metabolism and high efficiency of light and water use, and as an energy crop for bioethanol production. Transcriptome resources from economically valuable Agaves species, such as Agave tequilana and A. salmiana, as well as hybrids for fibers, are now available, and multiple gene expression landscape analyses have been reported. Key components in molecular mechanisms underlying drought tolerance could be uncovered by analyzing gene expression patterns of roots. This study describes an efficient protocol for high-quality total RNA isolation from phenolic compounds-rich Agave roots. Our methodology involves suitable root handling and collecting in the field and using saving-time commercial kits available. RNA isolated from roots free of lignified out-layers and clean cortex showed high values of quality and integrity according to electrophoresis and microfluidics-based platform. Synthesis of long full-length cDNAs and PCR amplification tested the suitability for downstream applications of extracted RNA. The protocol was applied successfully to A. tequilana roots but can be used for other Agave species that also develop lignified epidermis/exodermis in roots.

PAPhys are phytases that belong to the family of purple acid phosphatases (PAPs). Currently, there is limited information about how to distinguish a PAPhy sequence from a PAP without phytase activity. PAPhys have only been identified in angiosperm plants; therefore, it is unknown how prevalent they are in other photosynthetic organisms and what their phylogenetic relationships are. Therefore, the purpose of this research is to discover potential PAPhys using novel PAPhy-associated motifs and provide new information on their phylogeny. Experimentally validated PAPhys were identified and analyzed using the neighbor-joining phylogeny, and from each clade, a consensus sequence (CS) was generated. CSs were analyzed with the Multiple Em for Motif Elicitation (MEME) tool to identify characteristic motifs; with these motifs, potential PAPhys were searched in a broad range of photosynthetic organisms using the Motif Alignment & Search Tool (MAST). Genes encoding potential PAPhys were analyzed by maximum-likelihood phylogeny. Twenty-one PAPhys were identified in the literature and resolved into five clades. The MEME analysis generated four PAPhy-associated motifs (M1 to M4). A total of 160 potential PAPhys were identified with MAST. The phylogeny of these potential PAPhys revealed 11 major clades, among which six clades do not possess any experimentally studied PAPhy. The PAPhy motifs identified in this work were useful for identifying potential PAPhys. Phylogenetic analysis suggested that PAPhys are an ancient class of proteins in photosynthetic organisms. These new insights have provided valuable information on the evolution of this class of phytases and may assist in the cloning of novel PAPhys, which could lead to biotechnological applications.

Most of the commercial papaya genotypes show susceptibility to water deficit stress and require high volumes of irrigation water to yield properly. To tackle this problem, we have collected wild native genotypes of Carica papaya that have proved to show better physiological performance under water deficit stress than the commercial cultivar grown in Mexico. In the present study, plants from a wild Carica papaya genotype and a commercial genotype were subjected to water deficit stress (WDS), and their response was characterized in physiological and molecular terms. The physiological parameters measured (water potential, photosynthesis, Fv/Fm and electrolyte leakage) confirmed that the papaya wild genotype showed better physiological responses than the commercial one when exposed to WDS. Subsequently, RNA-Seq was performed for 4 cDNA libraries in both genotypes (susceptible and tolerant) under well-watered conditions, and when they were subjected to WDS for 14 days. Consistently, differential expression analysis revealed that after 14 days of WDS, the wild tolerant genotype had a higher number of up-regulated genes, and a higher number of transcription factors (TF) that were differentially expressed in response to WDS, than the commercial genotype. Thus, six TF genes (CpHSF, CpMYB, CpNAC, CpNFY-A, CpERF and CpWRKY) were selected for further qRT-PCR analysis as they were highly expressed in response to WDS in the wild papaya genotype. qRT-PCR results confirmed that the wild genotype had higher expression levels (REL) in all 6 TF genes than the commercial genotype. Our transcriptomic analysis should help to unravel candidate genes that may be useful in the development of new drought-tolerant cultivars of this important tropical crop.

Maize stalk rot caused by Fusarium graminearum can lead to severe losses and accumulation of mycotoxins with detrimental effects on livestock health. Because few management strategies are available, the development of resistant varieties is considered the most cost-effective way to control the disease. However, the stalk-tissue-specific mechanisms underlying resistance to F. graminearum remain poorly understood, although it is believed to be strongly influenced by environmental factors. In this study, we performed a temporal transcriptome analysis of two maize inbred lines with contrasting responses to stalk rot using gene expression profiling. We observed differential downregulation of gene expression during the first 2 weeks in a resistant inbred line inoculated with F. graminearum. Time-course gene ontology enrichment analysis suggests that resistance may be caused by a modulation of gene expression associated with redox homeostasis, hormone biosynthesis, cytoskeleton activity, and cell wall remodeling. We validated our gene expression profiling data by measuring the expression of 10 differentially expressed genes using quantitative reverse-transcription PCR. Our analyses also revealed the effect of two environmental conditions with contrasting temperatures and relative humidity on the resistant phenotype and gene expression. This research expands our knowledge of molecular events underlying resistance to stalk rot and the effect of environmental conditions on the disease interaction. Our findings can be exploited for the development of resistant varieties. [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 .

Currently, there is a need of non-computationally-intensive bioinformatics tools to cope with the increase of large datasets produced by Next Generation Sequencing technologies. We present a simple and robust bioinformatics pipeline to search for novel enzymes in metagenomic sequences. The strategy is based on pattern searching using as reference conserved motifs coded as regular expressions. As a case study, we applied this scheme to search for novel proteases S8A in a publicly available metagenome. Briefly, (1) the metagenome was assembled and translated into amino acids; (2) patterns were matched using regular expressions; (3) retrieved sequences were annotated; and (4) diversity analyses were conducted. Following this pipeline, we were able to identify nine sequences containing an S8 catalytic triad, starting from a metagenome containing 9,921,136 Illumina reads. Identity of these nine sequences was confirmed by BLASTp against databases at NCBI and MEROPS. Identities ranged from 62 to 89% to their respective nearest ortholog, which belonged to phyla Proteobacteria, Actinobacteria, Planctomycetes, Bacterioidetes, and Cyanobacteria, consistent with the most abundant phyla reported for this metagenome. All these results support the idea that they all are novel S8 sequences and strongly suggest that our methodology is robust and suitable to detect novel enzymes.

The diversity and community structure of arbuscular mycorrhizal fungi (AMF) associated with coconut (Cocos nucifera) roots was evaluated by next generation sequencing (NGS) using partial sequences of the 18S rDNA gene and by spore isolation and morphological identification from rhizosphere soil. Root samples from six different Green Dwarf coconut plantations and from one organic plantation surrounded by tropical dry forest along the coastal sand dunes in Yucatan, Mexico, were collected during the rainy and dry seasons. In total, 14 root samples were sequenced with the Illumina MiSeq platform. Additionally, soil samples from the dry season were collected to identify AMF glomerospores. Based on a 95–97% similarity, a total of 36 virtual taxa (VT) belonging to nine genera were identified including one new genus-like clade. Glomus was the most abundant genus, both in number of VT and sequences. The comparison of dry and rainy season samples revealed differences in the richness and composition of AMF communities colonizing coconut roots. Our study shows that the main AMF genera associated with coconut tree roots in all samples were Glomus, Sclerocystis, Rhizophagus, Redeckera, and Diversispora. Based on glomerospore morphology, 22 morphospecies were recorded among which 14 were identified to species. Sclerocystis sinuosa, Sclerocystis rubiformis, Glomus microaggregatum, and Acaulospora scrobiculata were dominant in field rhizosphere samples. This is the first assessment of the composition of AMF communities colonizing coconut roots in rainy and dry seasons. It is of importance for selection of AMF species to investigate for their potential application in sustainable agriculture of coconut.

Background Somatic embryogenesis (SE) is a useful biotechnological tool to study the morpho-physiological, biochemical and molecular processes during the development of Coffea canephora. Plant growth regulators (PGR) play a key role during cell differentiation in SE. The Auxin-response-factor (ARF) and Auxin/Indole-3-acetic acid (Aux/IAA) are fundamental components involved in the signaling of the IAA. The IAA signaling pathway activates or represses the expression of genes responsive to auxins during the embryogenic transition of the somatic cells. The growing development of new generation sequencing technologies (NGS), as well as bioinformatics tools, has allowed us to broaden the landscape of SE study of various plant species and identify the genes directly involved. Methods Analysis of transcriptome expression profiles of the C. canephora genome and the identification of a particular set of differentially expressed genes (DEG) during SE are described in this study. Results A total of eight ARF and seven Aux/IAA differentially expressed genes were identified during the different stages of the SE induction process. The quantitative expression analysis showed that ARF18 and ARF5 genes are highly expressed after 21 days of the SE induction, while Aux/IAA7 and Aux/IAA12 genes are repressed. Discussion The results of this study allow a better understanding of the genes involved in the auxin signaling pathway as well as their expression profiles during the SE process.

Both suberin and its associated waxes contribute to the formation of apoplastic barriers that protect plants from the environment. Some transcription factors have emerged as regulators of the suberization process. The potato StNAC103 gene was reported as a repressor of suberin polyester and suberin-associated waxes deposition because its RNAi-mediated downregulation (StNAC103-RNAi) over-accumulated suberin and associated waxes in the tuber phellem concomitantly with the induction of representative biosynthetic genes. Here, to explore if other genes of the large NAC gene family participate to this repressive function, we extended the silencing to other NAC members by targeting the conserved NAC domain of StNAC103 (StNAC103-RNAi-c). Transcript profile of the StNAC103-RNAi-c phellem indicated that StNAC101 gene was an additional potential target. In comparison with StNAC103-RNAi, the silencing with StNAC103-RNAi-c construct resulted in a similar effect in suberin but yielded an increased load of associated waxes in tuber phellem, mainly alkanes and feruloyl esters. Globally, the chemical effects in both silenced lines are supported by the transcript accumulation profile of genes involved in the biosynthesis, transport and regulation of apoplastic lipids. In contrast, the genes of polyamine biosynthesis were downregulated. Altogether these results point out to StNAC101 as a candidate to repress the suberin-associated waxes.

Pseudoperonospora humuli is an obligate oomycete pathogen of hop (Humulus lupulus) that causes downy mildew, an important disease in most production regions in the Northern Hemisphere. The pathogen can cause a systemic infection in hop, overwinter in the root system, and infect propagation material. Substantial yield loss may occur owing to P. humuli infection of strobiles (seed cones), shoots, and cone-bearing branches. Fungicide application and cultural practices are the primary methods to manage hop downy mildew. However, effective, sustainable, and cost-effective management of downy mildew can be improved by developing early detection systems to inform on disease risk and timely fungicide application. However, no species-specific diagnostic assays or genomic resources are available for P. humuli. The genome of the P. humuli OR502AA isolate was partially sequenced using Illumina technology and assembled with ABySS. The assembly had a minimum scaffold length of 500 bp and an N50 (median scaffold length of the assembled genome) of 19.2 kbp. A total number of 18,656 genes were identified using MAKER standard gene predictions. Additionally, transcriptome assemblies were generated using RNA-seq and Trinity for seven additional P. humuli isolates. Bioinformatics analyses of next generation sequencing reads of P. humuli and P. cubensis (a closely related sister species) identified 242 candidate species-specific P. humuli genes that could be used as diagnostic molecular markers. These candidate genes were validated using polymerase chain reaction against a diverse collection of isolates from P. humuli, P. cubensis, and other oomycetes. Overall, four diagnostic markers were found to be uniquely present in P. humuli. These candidate markers identified through comparative genomics can be used for pathogen diagnostics in propagation material, such as rhizomes and vegetative cuttings, or adapted for biosurveillance of airborne sporangia, an important source of inoculum in hop downy mildew epidemics.

RNA sequencing (RNA-seq) coupled to DNA methylation strategies enables the detection and characterization of genes which expression levels might be mediated by DNA methylation. Here we describe a bioinformatics protocol to analyze gene expression levels using RNA-seq data that allow us to identify candidate genes to be tested by bisulfite assays. The candidate methylated genes are usually those that are low expressed in a particular condition or developmental stage.

Drought is the main single factor affecting crop production worldwide. Water-deficit stress causes severe impairment of cell homeostasis, resulting from osmotic and oxidative stress. However, plants have evolved complex mechanisms to deal with this stress. In the case of papaya (Carica papaya), commercial cultivars show susceptibility to drought stress. Given that Central America and the south of Mexico appear to be the center of origin of this species, the possibility exists that native populations collected in Yucatan may show better responses to drought than their commercial counterparts. In fact, our genotypes collected from undisturbed regions of Yucatan have proved to be physiologically more tolerant to water-deficit stress than those commercial cultivars commonly grown in Mexico. In order to obtain insights into the genetic mechanisms that C. papaya may use to deal with water-deficit stress; a transcriptome (RNA-Seq) approach was used with C. papaya leaves. Total RNA was obtained from two different genotypes, one that is susceptible to water deficit and another genotype that shows tolerance to water-deficit stress. Four cDNA libraries were constructed, from leaf tissue from plants from the two genotypes grown either under irrigation or under water deficit for 14 days. The quality of those libraries was confirmed by bioanalyzer (Agilent 2100). The cDNA libraries were then sequenced using MiSeq (Illumina platform). A total of 20,140,717 reads were obtained; these reads were pre-processed to remove adaptors and low-quality reads, and 11,804,864 clean sequences were assembled into 135,459 contigs. After removing isoforms and short contigs (<300 bp), a total of 29,070 sequences were identified as genes that represent the transcriptome of C. papaya. Sequence similarity searches against the Arabidopsis proteome revealed that 28,878 (99%) of 29,070 genes had significant matches (TAIR 10; bit-score >= 90), confirming the good quality of the transcriptome. Thus, these data provide a valuable source of information for improving our understanding of the molecular mechanisms by which C. papaya overcomes water-deficit stress.

Somatic embryogenesis (SE) is one of the most studied developmental processes due to its applications, such as plant micropropagation, transformation, and germplasm conservation. The use of massive techniques of sequencing, as well as the use of subtractive hybridization and macroarrays, has led to the identification of hundreds of genes involved in the SE process. These have been important developments to study the molecular aspects of the progress of SE. With the advent of the new massive techniques for sequencing RNA, it has been possible to see a more complete picture of whole processes. In this chapter we present a technique to handle the elaboration of the transcriptome from the extraction of RNA until the assembly of the complete transcriptome.

Tocopherols, tocotrienols, and plastochromanols (collectively termed tocochromanols) are lipid-soluble antioxidants synthesized by all plants. Their dietary intake, primarily from seed oils, provides vitamin E and other health benefits. Tocochromanol biosynthesis has been dissected in the dicot Arabidopsis thaliana, which has green, photosynthetic seeds, but our understanding of tocochromanol accumulation in major crops, whose seeds are nonphotosynthetic, remains limited. To understand the genetic control of tocochromanols in grain, we conducted a joint linkage and genome-wide association study in the 5000-line U.S. maize (Zea mays) nested association mapping panel. Fifty-two quantitative trait loci for individual and total tocochromanols were identified, and of the 14 resolved to individual genes, six encode novel activities affecting tocochromanols in plants. These include two chlorophyll biosynthetic enzymes that explain the majority of tocopherol variation, which was not predicted given that, like most major cereal crops, maize grain is nonphotosynthetic. This comprehensive assessment of natural variation in vitamin E levels in maize establishes the foundation for improving tocochromanol and vitamin E content in seeds of maize and other major cereal crops.