Transcriptomics and co-expression networks reveal tissue-specific responses and regulatory hubs under mild and severe drought in papaya (Carica papaya L.)

MAD-FC: A fold change visualization with readability, proportionality, and symmetry

We propose a fold change transform that demonstrates a combination of visualization properties exhibited by log and linear plots of fold change. A fold change visualization should ideally exhibit: (1) readability, where fold change values are recoverable from datapoint position; (2) proportionality, where fold change values of the same direction are proportionally distant from the point of no change; (3) symmetry, where positive and negative fold changes of the same magnitude are equidistant to the point of no change; and (4) high dynamic range, where datapoint values are distinguishable across orders of magnitude within a fixed plot area and pixel resolution. A linear visualization has readability and partial proportionality but lacks high dynamic range and symmetry (because negative direction fold changes are bound between [0, 1] while positive are between (1, ∞)). Log plots of fold change have partial readability, high dynamic range, and symmetry, but lack proportionality because of the log transform. We outline a new transform, named mirrored axis distortion of fold change (MAD-FC), that extends a linear visualization of fold change data to exhibit readability, proportionality, and symmetry (but still has the limited dynamic range of linear plots). We illustrate the use of MAD-FC with biomedical data using various fold change plots. We argue that MAD plots may be a more useful visualization than log or linear plots for applications that do not require a high dynamic range (less than 8 units in log2 space).

Multiomics analysis reveals a substantial decrease in nanoplastics uptake and associated impacts by nano zinc oxide in fragrant rice (Oryza sativa L.)

Nanoplastics (NPs) have emerged as global environmental pollutants with concerning implications for sustainable agriculture. Understanding the underlying mechanisms of NPs toxicity and devising strategies to mitigate their impact is crucial for crop growth and development. Here, we investigated the nanoparticles of zinc oxide (nZnO) to mitigate the adverse effects of 80 nm NPs on fragrant rice. Our results showed that optimized nZnO (25 mg L-1) concentration rescued root length and structural deficits by improving oxidative stress response, antioxidant defense mechanism and balanced nutrient levels, compared to seedlings subjected only to NPs stress (50 mg L-1). Consequently, microscopy observations, Zeta potential and Fourier transform infrared (FTIR) results revealed that NPs were mainly accumulated on the initiation joints of secondary roots and between cortical cells that blocks the nutrients uptake, while the supplementation of nZnO led to the formation of aggregates with NPs, which effectively impedes the uptake of NPs by the roots of fragrant rice. Transcriptomic analysis identified a total of 3973, 3513 and 3380 differentially expressed genes (DEGs) in response to NPs, nZnO and NPs+nZnO, respectively, compared to the control. Moreover, DEGs were significantly enriched in multiple pathways including biosynthesis of secondary metabolite, phenylpropanoid biosynthesis, amino sugar and nucleotide sugar metabolism, carotenoid biosynthesis, plant-pathogen interactions, MAPK signaling pathway, starch and sucrose metabolism, and plant hormone signal transduction. These pathways could play a significant role in alleviating NPs toxicity and restoring fragrant rice roots. Furthermore, metabolomic analysis demonstrated that nZnO application restored 2-acetyl-1-pyrroline (2-AP) pathways genes expression, enzymatic activities, and the content of essential precursors related to 2-AP biosynthesis under NPs toxicity, which ultimately led to the restoration of 2-AP content in the leaves. In conclusion, this study shows that optimized nZnO application effectively alleviates NPs toxic effects and restores both root structure and aroma production in fragrant rice leaves. This research offers a sustainable and practical strategy to enhance crop production under NPs toxicity while emphasizing the pivotal role of essential micronutrient nanomaterials in agriculture.

Bioinformatics and meta-analysis of expression data to investigate transcriptomic response of wheat root to abiotic stresses

Abiotic stresses are predominant and main causes of the losses in the crop yield. A complexity of systems biology and involvement of numerous genes in the response to abiotic factors have challenged efforts to create tolerant cultivars with sustainable production. The root is the main organ of the plant and determines a plant tolerance under stressful conditions. In this study, we carried out a meta-analysis of expression datasets from wheat root to identify differentially expressed genes, followed by the weighted gene co-expression network analysis (WGCNA) to construct the weighted gene co-expression network. The aim was to identify consensus differentially expressed genes with regulatory functions, gene networks, and biological pathways involved in response of wheat root to a set of abiotic stresses. The meta-analysis using Fisher method (FDR<0.05) identified consensus 526 DEGs from 55,367 probe sets. Although the annotated expression data are limited for wheat, the functional analysis based on the data from model plants could identify the up-regulated seven regulatory genes involved in chromosome organization and response to oxygen-containing compounds. WGCNA identified four gene modules that were mostly associated with the ribosome biogenesis and polypeptide synthesis. This study's findings enhance our understanding of key players and gene networks related to wheat root response to multiple abiotic stresses.

Analysis of Carica papaya Informs Lineage-Specific Evolution of the Aquaporin (AQP) Family in Brassicales

Aquaporins (AQPs), a type of intrinsic membrane proteins that transport water and small solutes across biological membranes, play crucial roles in plant growth and development. This study presents a first genome-wide identification and comparative analysis of the AQP gene family in papaya (Carica papaya L.), an economically and nutritionally important fruit tree of tropical and subtropical regions. A total of 29 CpAQP genes were identified, which represent five subfamilies, i.e., nine plasma intrinsic membrane proteins (PIPs), eight tonoplast intrinsic proteins (TIPs), seven NOD26-like intrinsic proteins (NIPs), two X intrinsic proteins (XIPs), and three small basic intrinsic proteins (SIPs). Although the family is smaller than the 35 members reported in Arabidopsis, it is highly diverse, and the presence of CpXIP genes as well as orthologs in Moringa oleifera and Bretschneidera sinensis implies that the complete loss of the XIP subfamily in Arabidopsis is lineage-specific, sometime after its split with papaya but before Brassicaceae-Cleomaceae divergence. Reciprocal best hit-based sequence comparison of 530 AQPs and synteny analyses revealed that CpAQP genes belong to 29 out of 61 identified orthogroups, and lineage-specific evolution was frequently observed in Brassicales. Significantly, the well-characterized NIP3 group was completely lost; lineage-specific loss of the NIP8 group in Brassicaceae occurred sometime before the divergence with Cleomaceae, and lineage-specific loss of NIP2 and SIP3 groups in Brassicaceae occurred sometime after the split with Cleomaceae. In contrast to a predominant role of recent whole-genome duplications (WGDs) on the family expansion in B. sinensis, Tarenaya hassleriana, and Brassicaceae plants, no recent AQP repeats were identified in papaya, and ancient WGD repeats are mainly confined to the PIP subfamily. Subfamily even group-specific evolution was uncovered via comparing exon-intron structures, conserved motifs, the aromatic/arginine selectivity filter, and gene expression profiles. Moreover, down-regulation during fruit ripening and expression divergence of duplicated CpAQP genes were frequently observed in papaya. These findings will not only improve our knowledge on lineage-specific family evolution in Brassicales, but also provide valuable information for further studies of AQP genes in papaya and species beyond.

Identification and characterization of RuvBL DNA helicase genes for tolerance against abiotic stresses in bread wheat (Triticum aestivum L.) and related species

Recombination UVB (sensitivity) like (RuvBL) helicase genes represent a conserved family of genes, which are known to be involved in providing tolerance against abiotic stresses like heat and drought. We identified nine wheat RuvBL genes, one each on nine different chromosomes, belonging to homoeologous groups 2, 3, and 4. The lengths of genes ranged from 1647 to 2197 bp and exhibited synteny with corresponding genes in related species including Ae. tauschii, Z. mays, O. sativa, H. vulgare, and B. distachyon. The gene sequences were associated with regulatory cis-elements and transposable elements. Two genes, namely TaRuvBL1a-4A and TaRuvBL1a-4B, also carried targets for a widely known miRNA, tae-miR164. Gene ontology revealed that these genes were closely associated with ATP-dependent formation of histone acetyltransferase complex. Analysis of the structure and function of RuvBL proteins revealed that the proteins were localized mainly in the cytoplasm. A representative gene, namely TaRuvBL1a-4A, was also shown to be involved in protein-protein interactions with ten other proteins. On the basis of phylogeny, RuvBL proteins were placed in two sub-divisions, namely RuvBL1 and RuvBL2, which were further classified into clusters and sub-clusters. In silico studies suggested that these genes were differentially expressed under heat/drought. The qRT-PCR analysis confirmed that expression of TaRuvBL genes differed among wheat cultivars, which differed in the level of thermotolerance. The present study advances our understanding of the biological role of wheat RuvBL genes and should help in planning future studies on RuvBL genes in wheat including use of RuvBL genes in breeding thermotolerant wheat cultivars.

Network Biology Analyses and Dynamic Modeling of Gene Regulatory Networks under Drought Stress Reveal Major Transcriptional Regulators in Arabidopsis

Drought is one of the most serious abiotic stressors in the environment, restricting agricultural production by reducing plant growth, development, and productivity. To investigate such a complex and multifaceted stressor and its effects on plants, a systems biology-based approach is necessitated, entailing the generation of co-expression networks, identification of high-priority transcription factors (TFs), dynamic mathematical modeling, and computational simulations. Here, we studied a high-resolution drought transcriptome of Arabidopsis. We identified distinct temporal transcriptional signatures and demonstrated the involvement of specific biological pathways. Generation of a large-scale co-expression network followed by network centrality analyses identified 117 TFs that possess critical properties of hubs, bottlenecks, and high clustering coefficient nodes. Dynamic transcriptional regulatory modeling of integrated TF targets and transcriptome datasets uncovered major transcriptional events during the course of drought stress. Mathematical transcriptional simulations allowed us to ascertain the activation status of major TFs, as well as the transcriptional intensity and amplitude of their target genes. Finally, we validated our predictions by providing experimental evidence of gene expression under drought stress for a set of four TFs and their major target genes using qRT-PCR. Taken together, we provided a systems-level perspective on the dynamic transcriptional regulation during drought stress in Arabidopsis and uncovered numerous novel TFs that could potentially be used in future genetic crop engineering programs.

Orchard Management and Incorporation of Biochemical and Molecular Strategies for Improving Drought Tolerance in Fruit Tree Crops

Water scarcity is one of the greatest concerns for agronomy worldwide. In recent years, many water resources have been depleted due to multiple factors, especially mismanagement. Water resource shortages lead to cropland expansion, which likely influences climate change and affects global agriculture, especially horticultural crops. Fruit yield is the final aim in commercial orchards; however, drought can slow tree growth and/or decrease fruit yield and quality. It is therefore necessary to find approaches to solve this problem. The main objective of this review is to discuss the most recent horticultural, biochemical, and molecular strategies adopted to improve the response of temperate fruit crops to water stress. We also address the viability of cultivating fruit trees in dry areas and provide precise protection methods for planting fruit trees in arid lands. We review the main factors involved in planting fruit trees in dry areas, including plant material selection, regulated deficit irrigation (DI) strategies, rainwater harvesting (RWH), and anti-water stress materials. We also provide a detailed analysis of the molecular strategies developed to combat drought, such as Clustered Regularly Interspaced Short Palindromic Repeat (CRISPR) through gene overexpression or gene silencing. Finally, we look at the molecular mechanisms associated with the contribution of the microbiome to improving plant responses to drought.

Transcriptome analysis reveals the molecular mechanisms of Phragmites australis tolerance to CuO-nanoparticles and/or flood stress induced by arbuscular mycorrhizal fungi

The molecular mechanism of arbuscular mycorrhizal fungi (AMF) in vertical flow constructed wetlands (VFCWs) for the purification of copper oxide nanoparticles (CuO-NPs) contaminated wastewater remains unclear. In this study, transcriptome analysis was used to explore the effect of AMF inoculation on the gene expression profile of Phragmites australis roots under different concentrations of CuO-NPs and/or flood stress. 551, 429 and 2281 differentially expressed genes (DEGs) were specially regulated by AMF under combined stresses of CuO-NPs and flood, single CuO-NPs stress and single flood stress, respectively. Based on the results of DEG function annotation and enrichment analyses, AMF inoculation under CuO-NPs and/or flood stress up-regulated the expression of a number of genes involved in antioxidant defense systems, cell wall biosynthesis and transporter protein, which may contribute to plant tolerance. The expression of 30 transcription factors (TFs) was up-regulated by AMF inoculation under combined stresses of CuO-NPs and flood, and 44 and 44 TFs were up-regulated under single CuO-NPs or flood condition, respectively, which may contribute to the alleviating effect of symbiosis on CuO-NPs and/or flood stress. These results provided a theoretical basis for enhancing the ecological restoration function of wetland plants for metallic nanoparticles (MNPs) by mycorrhizal technology in the future.

Tissue-specific mRNA profiling of the Brassica napus-Sclerotinia sclerotiorum interaction uncovers novel regulators of plant immunity

White mold is caused by the fungal pathogen Sclerotinia sclerotiorum and leads to rapid and significant loss in plant yield. Among its many brassicaceous hosts, including Brassica napus (canola) and Arabidopsis, the response of individual tissue layers directly at the site of infection has yet to be explored. Using laser microdissection coupled with RNA sequencing, we profiled the epidermis, mesophyll, and vascular leaf tissue layers of B. napus in response to S. sclerotiorum. High-throughput tissue-specific mRNA sequencing increased the total number of detected transcripts compared with whole-leaf assessments and provided novel insight into the conserved and specific roles of ontogenetically distinct leaf tissue layers in response to infection. When subjected to pathogen infection, the epidermis, mesophyll, and vasculature activate both specific and shared gene sets. Putative defense genes identified through transcription factor network analysis were then screened for susceptibility against necrotrophic, hemi-biotrophic, and biotrophic pathogens. Arabidopsis deficient in PR5-like RECEPTOR KINASE (PR5K) mRNA levels were universally susceptible to all pathogens tested and were further characterized to identify putative interacting partners involved in the PR5K signaling pathway. Together, these data provide insight into the complexity of the plant defense response directly at the site of infection.

Genomic Analysis of LEA Genes in Carica papaya and Insight into Lineage-Specific Family Evolution in Brassicales

Late embryogenesis abundant (LEA) proteins comprise a diverse superfamily involved in plant development and stress responses. This study presents a first genome-wide analysis of LEA genes in papaya (Carica papaya L., Caricaceae), an economically important tree fruit crop widely cultivated in the tropics and subtropics. A total of 28 members were identified from the papaya genome, which belong to eight families with defined Pfam domains, i.e., LEA_1 (3), LEA_2 (4), LEA_3 (5), LEA_4 (5), LEA_5 (2), LEA_6 (2), DHN (4), and SMP (3). The family numbers are comparable to those present in Ricinus communis (Euphorbiaceae, 28) and Moringa oleifera (Moringaceae, 29), but relatively less than that found in Moringa oleifera (Cleomaceae, 39) and Arabidopsis thaliana (Brassicaceae, 51), implying lineage-specific evolution in Brassicales. Indeed, best-reciprocal-hit-based sequence comparison and synteny analysis revealed the presence of 29 orthogroups, and significant gene expansion in Tarenaya and Arabidopsis was mainly contributed by whole-genome duplications that occurred sometime after their split with the papaya. Though a role of transposed duplication was also observed, tandem duplication was shown to be a key contributor in gene expansion of most species examined. Further comparative analyses of exon-intron structures and protein motifs supported fast evolution of this special superfamily, especially in Arabidopsis. Transcriptional profiling revealed diverse expression patterns of CpLEA genes over various tissues and different stages of developmental fruit. Moreover, the transcript level of most genes appeared to be significantly regulated by drought, cold, and salt stresses, corresponding to the presence of cis-acting elements associated with stress response in their promoter regions. These findings not only improve our knowledge on lineage-specific family evolution in Brassicales, but also provide valuable information for further functional analysis of LEA genes in papaya.

Myco-Synergism Boosts Herbivory-Induced Maize Defense by Triggering Antioxidants and Phytohormone Signaling

BACKGROUND

Biocontrol strategies are the best possible and eco-friendly solution to develop resistance against O furnacalis and improve the maize yield. However, the knowledge about underlying molecular mechanisms, metabolic shifts, and hormonal signaling is limited.

METHODS

Here, we used an axenic and a consortium of entomopathogenic Beauveria bassiana OFDH1-5 and a pathogen-antagonistic Trichoderma asperellum GDFS1009 in maize and observed that consortium applications resulted in higher chlorophyll contents and antioxidants activities [superoxide dismutase (SOD), peroxidase (POD), proline, protease, and polyphenol oxidase (PPO)] with a decrease in O. furnacalis survival. We performed a comprehensive transcriptome and an untargeted metabolome profiling for the first time at a vegetative stage in fungal inoculated maize leaves at 0-, 12-, 24-, 48-, and 72-h post insect infestation.

RESULTS

The consortium of B. bassiana and T. asperellum leads to 80-95% of O. furnacalis mortality. A total of 13,156 differentially expressed genes were used for weighted gene coexpression network analysis. We identified the six significant modules containing thirteen candidate genes [protein kinase (GRMZM2G025459), acyl-CoA dehydrogenase (GRMZM5G864319), thioredoxin gene (GRMZM2G091481), glutathione S-transferase (GRMZM2G116273), patatin-like phospholipase gene (GRMZM2G154523), cytochrome P450 (GRMZM2G139874), protease inhibitor (GRMZM2G004466), (AC233926.1_FG002), chitinase (GRMZM2G453805), defensin (GRMZM2G392863), peroxidase (GRMZM2G144153), GDSL- like lipase (AC212068.4_FG005), and Beta-glucosidase (GRMZM2G031660)], which are not previously reported that are highly correlated with Jasmonic acid - Ethylene (JA-ET) signaling pathway and antioxidants. We detected a total of 130 negative and 491 positive metabolomic features using a ultrahigh-performance liquid chromatography ion trap time-of-flight mass spectrometry (UHPLC-QTOF-MS). Intramodular significance and real time-quantitative polymerase chain reaction (RT-qPCR) expressions showed that these genes are the true candidate genes. Consortium treated maize had higher jasmonic acid (JA), salicylic acid (SA), and ethylene (ET) levels.

CONCLUSION

Our results provide insights into the genetics, biochemicals, and metabolic diversity and are useful for future biocontrol strategies against ACB attacks.

Phylogenomic and Microsynteny Analysis Provides Evidence of Genome Arrangements of High-Affinity Nitrate Transporter Gene Families of Plants

Nitrate transporter 2 (NRT2) and NRT3 or nitrate-assimilation-related 2 (NAR2) proteins families form a two-component, high-affinity nitrate transport system, which is essential for the acquisition of nitrate from soils with low N availability. An extensive phylogenomic analysis across land plants for these families has not been performed. In this study, we performed a microsynteny and orthology analysis on the NRT2 and NRT3 genes families across 132 plants (Sensu lato) to decipher their evolutionary history. We identified significant differences in the number of sequences per taxonomic group and different genomic contexts within the NRT2 family that might have contributed to N acquisition by the plants. We hypothesized that the greater losses of NRT2 sequences correlate with specialized ecological adaptations, such as aquatic, epiphytic, and carnivory lifestyles. We also detected expansion on the NRT2 family in specific lineages that could be a source of key innovations for colonizing contrasting niches in N availability. Microsyntenic analysis on NRT3 family showed a deep conservation on land plants, suggesting a high evolutionary constraint to preserve their function. Our study provides novel information that could be used as guide for functional characterization of these gene families across plant lineages.

Genomic Approaches for Improvement of Tropical Fruits: Fruit Quality, Shelf Life and Nutrient Content

The breeding of tropical fruit trees for improving fruit traits is complicated, due to the long juvenile phase, generation cycle, parthenocarpy, polyploidy, polyembryony, heterozygosity and biotic and abiotic factors, as well as a lack of good genomic resources. Many molecular techniques have recently evolved to assist and hasten conventional breeding efforts. Molecular markers linked to fruit development and fruit quality traits such as fruit shape, size, texture, aroma, peel and pulp colour were identified in tropical fruit crops, facilitating Marker-assisted breeding (MAB). An increase in the availability of genome sequences of tropical fruits further aided in the discovery of SNP variants/Indels, QTLs and genes that can ascertain the genetic determinants of fruit characters. Through multi-omics approaches such as genomics, transcriptomics, metabolomics and proteomics, the identification and quantification of transcripts, including non-coding RNAs, involved in sugar metabolism, fruit development and ripening, shelf life, and the biotic and abiotic stress that impacts fruit quality were made possible. Utilizing genomic assisted breeding methods such as genome wide association (GWAS), genomic selection (GS) and genetic modifications using CRISPR/Cas9 and transgenics has paved the way to studying gene function and developing cultivars with desirable fruit traits by overcoming long breeding cycles. Such comprehensive multi-omics approaches related to fruit characters in tropical fruits and their applications in breeding strategies and crop improvement are reviewed, discussed and presented here.

Transcriptome-Based Identification and Functional Characterization of NAC Transcription Factors Responsive to Drought Stress in Capsicum annuum L

Capsicum annuum L. is one of the most cultivated Solanaceae species, and in the open field, water limitation leading to drought stress affects its fruit quality, fruit setting, fruit size and ultimately yield. We identified stage-specific and a common core set of differentially expressed genes, following RNA-seq transcriptome analyses of a breeding line subjected to acute drought stress followed by recovery (rewatering), at three stages of plant development. Among them, two NAC transcription factor (TF) genes, i.e., CaNAC072 and CaNAC104, were always upregulated after drought stress and downregulated after recovery. The two TF proteins were observed to be localized in the nucleus following their transient expression in Nicotiana benthamiana leaves. The expression of the two NACs was also induced by NaCl, polyethylene glycol (PEG) and abscisic acid (ABA) treatments, suggesting that CaNAC072 is an early, while CaNAC104 is a late abiotic stress-responsive gene. Virus-induced gene silencing (VIGS) of CaNAC104 did not affect the pepper plantlet’s tolerance to drought stress, while VIGS of CaNAC072 increased drought tolerance. Heterologous expression of CaNAC072 in Arabidopsis thaliana as well as in plants mutated for its homolog ANAC072 did not increase drought stress tolerance. This highlights a different role of the two NAC homologs in the two species. Here, we discuss the complex role of NACs as transcriptional switches in the response to drought stress in bell pepper.

Gene-Metabolite Network Analysis Revealed Tissue-Specific Accumulation of Therapeutic Metabolites in Mallotus japonicus

Mallotus japonicus is a valuable traditional medicinal plant in East Asia for applications as a gastrointestinal drug. However, the molecular components involved in the biosynthesis of bioactive metabolites have not yet been explored, primarily due to a lack of omics resources. In this study, we established metabolome and transcriptome resources for M. japonicus to capture the diverse metabolite constituents and active transcripts involved in its biosynthesis and regulation. A combination of untargeted metabolite profiling with data-dependent metabolite fragmentation and metabolite annotation through manual curation and feature-based molecular networking established an overall metabospace of M. japonicus represented by 2129 metabolite features. M. japonicus de novo transcriptome assembly showed 96.9% transcriptome completeness, representing 226,250 active transcripts across seven tissues. We identified specialized metabolites biosynthesis in a tissue-specific manner, with a strong correlation between transcripts expression and metabolite accumulations in M. japonicus. The correlation- and network-based integration of metabolome and transcriptome datasets identified candidate genes involved in the biosynthesis of key specialized metabolites of M. japonicus. We further used phylogenetic analysis to identify 13 C-glycosyltransferases and 11 methyltransferases coding candidate genes involved in the biosynthesis of medicinally important bergenin. This study provides comprehensive, high-quality multi-omics resources to further investigate biological properties of specialized metabolites biosynthesis in M. japonicus.

Comparative Transcriptome and Weighted Gene Co-expression Network Analysis Identify Key Transcription Factors of Rosa chinensis 'Old Blush' After Exposure to a Gradual Drought Stress Followed by Recovery

Rose is one of the most fundamental ornamental crops, but its yield and quality are highly limited by drought. The key transcription factors (TFs) and co-expression networks during rose’s response to drought stress and recovery after drought stress are still limited. In this study, the transcriptomes of leaves of 2-year-old cutting seedlings of Rosa chinensis ‘Old Blush’ from three continuous droughted stages (30, 60, 90 days after full watering) and rewatering were analyzed using RNA sequencing. Weighted gene co-expression network analysis (WGCNA) was used to construct a co-expression network, which was associated with the physiological traits of drought response to discovering the hub TFs involved in drought response. More than 45 million high-quality clean reads were generated from the sample and used for comparison with the rose reference genome. A total of 46433 differentially expressed genes (DEGs) were identified. Gene Ontology (GO) term enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis indicated that drought stress caused significant changes in signal transduction, plant hormones including ABA, auxin, brassinosteroid (BR), cytokinin, ethylene (ET), jasmonic acid (JA) and salicylic acid (SA), primary and secondary metabolism, and a certain degree of recovery after rewatering. Gene co-expression analysis identified 18 modules, in which four modules showed a high degree of correlation with physiological traits. In addition, 42 TFs including members of NACs, WRKYs, MYBs, AP2/ERFs, ARFs, and bHLHs with high connectivity in navajowhite1 and blue modules were screened. This study provides the transcriptome sequencing report of R. chinensis ‘Old Blush’ during drought stress and rewatering process. The study also identifies the response of candidate TFs to drought stress, providing guidelines for improving the drought tolerance of the rose through molecular breeding in the future.

An efficient papaya leaf distortion mosaic potyvirus vector for virus-induced gene silencing in papaya

Papaya (Carica papaya L.) is regarded as an excellent model for genomic studies of tropical trees because of its short generation time and its small genome that has been sequenced. However, functional genomic studies in papaya depend on laborious genetic transformations because no rapid tools exist for this species. Here, we developed a highly efficient virus-induced gene silencing (VIGS) vector for use in papaya by modifying an artificially attenuated infectious clone of papaya leaf distortion mosaic virus (PLDMV; genus: Potyvirus), PLDMV-E, into a stable Nimble Cloning (NC)-based PLDMV vector, pPLDMV-NC, in Escherichia coli. The target fragments for gene silencing can easily be cloned into pPLDMV-NC without multiple digestion and ligation steps. Using this PLDMV VIGS system, we silenced and characterized five endogenous genes in papaya, including two common VIGS marker genes, namely, phytoene desaturase, Mg-chelatase H subunit, putative GIBBERELLIN (GA)-INSENSITIVE DWARF1A and 1B encoding GA receptors; and the cytochrome P450 gene CYP83B1, which encodes a key enzyme involved in benzylglucosinolate biosynthesis. The results demonstrate that our newly developed PLDMV VIGS vector is a rapid and convenient tool for functional genomic studies in papaya.

Systematic analysis of differentially expressed ZmMYB genes related to drought stress in maize

MYB transcription factors play pivotal roles in hormone conduction signaling and abiotic stress response. In this study, 54 differentially expressed ZmMYB genes were identified and comprehensive analyses were conducted including gene's structure, chromosomal localization, phylogenetic tree, motif prediction, cis-elements and expression patterns. The results showed that 54 genes were unevenly distributed on 10 chromosomes and classified into eleven main subgroups by phylogenetic analysis, supported by motif and exon/intron analyses. The mainly stress-related cis-elements were ABRE, ARE, MBS and DRE-core. In addition, 8 core ZmMYB genes were identified by co-expression network. qRT-PCR results showed that the 8 ZmMYB genes exhibited different expression levels under different abiotic stresses, indicating that they were responsive to various abiotic stress. These results will provide insight for further functional investigation of ZmMYB genes.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-01013-2.

Phloem exudate metabolic content reflects the response to water-deficit stress in pea plants (Pisum sativum L.)

Drought stress impacts the quality and yield of Pisum sativum. Here, we show how short periods of limited water availability during the vegetative stage of pea alters phloem sap content and how these changes are connected to strategies used by plants to cope with water deficit. We have investigated the metabolic content of phloem sap exudates and explored how this reflects P. sativum physiological and developmental responses to drought. Our data show that drought is accompanied by phloem-mediated redirection of the components that are necessary for cellular respiration and the proper maintenance of carbon/nitrogen balance during stress. The metabolic content of phloem sap reveals a shift from anabolic to catabolic processes as well as the developmental plasticity of P. sativum plants subjected to drought. Our study underlines the importance of phloem-mediated transport for plant adaptation to unfavourable environmental conditions. We also show that phloem exudate analysis can be used as a useful proxy to study stress responses in plants. We propose that the decrease in oleic acid content within phloem sap could be considered as a potential marker of early signalling events mediating drought response.

Transcriptomic analysis reveals key transcription factors associated to drought tolerance in a wild papaya (Carica papaya) genotype

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.

Discovery of SNPs and InDels in papaya genotypes and its potential for marker assisted selection of fruit quality traits

Papaya is a tropical and climacteric fruit that is recognized for its nutritional benefits and medicinal applications. Its fruits ripen quickly and show a drastic fruit softening, leading to great post-harvest losses. To overcome this scenario, breeding programs of papaya must invest in exploring the available genetic variation to continue developing superior cultivars with improved fruit quality traits. The objective of this study was to perform a whole-genome genotyping (WGG) of papaya, predict the effects of the identified variants, and develop a list of ripening-related genes (RRGs) with linked variants. The Formosa elite lines of papaya Sekati and JS-12 were submitted to WGG with an Illumina Miseq platform. The effects of variants were predicted using the snpEff program. A total of 28,451 SNPs having Ts/Tv (Transition/Transversion) ratio of 2.45 and 1,982 small insertions/deletions (InDels) were identified. Most variant effects were predicted in non-coding regions, with only 2,104 and 138 effects placed in exons and splice site regions, respectively. A total of 106 RRGs were found to be associated with 460 variants, which may be converted into PCR markers to facilitate genetic mapping and diversity studies and to apply marker-assisted selection (MAS) for specific traits in papaya breeding programs.

What is quantitative plant biology?

Quantitative plant biology is an interdisciplinary field that builds on a long history of biomathematics and biophysics. Today, thanks to high spatiotemporal resolution tools and computational modelling, it sets a new standard in plant science. Acquired data, whether molecular, geometric or mechanical, are quantified, statistically assessed and integrated at multiple scales and across fields. They feed testable predictions that, in turn, guide further experimental tests. Quantitative features such as variability, noise, robustness, delays or feedback loops are included to account for the inner dynamics of plants and their interactions with the environment. Here, we present the main features of this ongoing revolution, through new questions around signalling networks, tissue topology, shape plasticity, biomechanics, bioenergetics, ecology and engineering. In the end, quantitative plant biology allows us to question and better understand our interactions with plants. In turn, this field opens the door to transdisciplinary projects with the society, notably through citizen science.

Exploring the complexity of soybean (Glycine max) transcriptional regulation using global gene co-expression networks

MAIN CONCLUSION

We report a soybean gene co-expression network built with data from 1284 RNA-Seq experiments, which was used to identify important regulators, modules and to elucidate the fates of gene duplicates. Soybean (Glycine max (L.) Merr.) is one of the most important crops worldwide, constituting a major source of protein and edible oil. Gene co-expression networks (GCN) have been extensively used to study transcriptional regulation and evolution of genes and genomes. Here, we report a soybean GCN using 1284 publicly available RNA-Seq samples from 15 distinct tissues. We found modules that are differentially regulated in specific tissues, comprising processes such as photosynthesis, gluconeogenesis, lignin metabolism, and response to biotic stress. We identified transcription factors among intramodular hubs, which probably integrate different pathways and shape the transcriptional landscape in different conditions. The top hubs for each module tend to encode proteins with critical roles, such as succinate dehydrogenase and RNA polymerase subunits. Importantly, gene essentiality was strongly correlated with degree centrality and essential hubs were enriched in genes involved in nucleic acids metabolism and regulation of cell replication. Using a guilt-by-association approach, we predicted functions for 93 of 106 hubs without functional description in soybean. Most of the duplicated genes had different transcriptional profiles, supporting their functional divergence, although paralogs originating from whole-genome duplications (WGD) are more often preserved in the same module than those from other mechanisms. Together, our results highlight the importance of GCN analysis in unraveling key functional aspects of the soybean genome, in particular those associated with hub genes and WGD events.

Gene co-expression network analysis to identify critical modules and candidate genes of drought-resistance in wheat

AIM

To establish a gene co-expression network for identifying principal modules and hub genes that are associated with drought resistance mechanisms, analyzing their mechanisms, and exploring candidate genes.

METHODS AND FINDINGS

42 data sets including PRJNA380841 and PRJNA369686 were used to construct the co-expression network through weighted gene co-expression network analysis (WGCNA). A total of 1,896,897,901 (284.30 Gb) clean reads and 35,021 differentially expressed genes (DEGs) were obtained from 42 samples. Functional enrichment analysis indicated that photosynthesis, DNA replication, glycolysis/gluconeogenesis, starch and sucrose metabolism, arginine and proline metabolism, and cell cycle were significantly influenced by drought stress. Furthermore, the DEGs with similar expression patterns, detected by K-means clustering, were grouped into 29 clusters. Genes involved in the modules, such as dark turquoise, yellow, and brown, were found to be appreciably linked with drought resistance. Twelve central, greatly correlated genes in stage-specific modules were subsequently confirmed and validated at the transcription levels, including TraesCS7D01G417600.1 (PP2C), TraesCS5B01G565300.1 (ERF), TraesCS4A01G068200.1 (HSP), TraesCS2D01G033200.1 (HSP90), TraesCS6B01G425300.1 (RBD), TraesCS7A01G499200.1 (P450), TraesCS4A01G118400.1 (MYB), TraesCS2B01G415500.1 (STK), TraesCS1A01G129300.1 (MYB), TraesCS2D01G326900.1 (ALDH), TraesCS3D01G227400.1 (WRKY), and TraesCS3B01G144800.1 (GT).

CONCLUSIONS

Analyzing the response of wheat to drought stress during different growth stages, we have detected three modules and 12 hub genes that are associated with drought resistance mechanisms, and five of those genes are newly identified for drought resistance. The references provided by these modules will promote the understanding of the drought-resistance mechanism. In addition, the candidate genes can be used as a basis of transgenic or molecular marker-assisted selection for improving the drought resistance and increasing the yields of wheat.

Gene coexpression network analysis and tissue-specific profiling of gene expression in jute (Corchorus capsularis L.)

BACKGROUND

Jute (Corchorus spp.), belonging to the Malvaceae family, is an important natural fiber crop, second only to cotton, and a multipurpose economic crop. Corchorus capsularis L. is one of the only two commercially cultivated species of jute. Gene expression is spatiotemporal and is influenced by many factors. Therefore, to understand the molecular mechanisms of tissue development, it is necessary to study tissue-specific gene expression and regulation. We used weighted gene coexpression network analysis, to predict the functional roles of gene coexpression modules and individual genes, including those underlying the development of different tissue types. Although several transcriptome studies have been conducted on C. capsularis, there have not yet been any systematic and comprehensive transcriptome analyses for this species.

RESULTS

There was significant variation in gene expression between plant tissues. Comparative transcriptome analysis and weighted gene coexpression network analysis were performed for different C. capsularis tissues at different developmental stages. We identified numerous tissue-specific differentially expressed genes for each tissue, and 12 coexpression modules, comprising 126 to 4203 genes, associated with the development of various tissues. There was high consistency between the genes in modules related to tissues, and the candidate upregulated genes for each tissue. Further, a gene network including 21 genes directly regulated by transcription factor OMO55970.1 was discovered. Some of the genes, such as OMO55970.1, OMO51203.1, OMO50871.1, and OMO87663.1, directly involved in the development of stem bast tissue.

CONCLUSION

We identified genes that were differentially expressed between tissues of the same developmental stage. Some genes were consistently up- or downregulated, depending on the developmental stage of each tissue. Further, we identified numerous coexpression modules and genes associated with the development of various tissues. These findings elucidate the molecular mechanisms underlying the development of each tissue, and will promote multipurpose molecular breeding in jute and other fiber crops.

Construction, comparison and evolution of networks in life sciences and other disciplines

Network approaches have become pervasive in many research fields. They allow for a more comprehensive understanding of complex relationships between entities as well as their group-level properties and dynamics. Many networks change over time, be it within seconds or millions of years, depending on the nature of the network. Our focus will be on comparative network analyses in life sciences, where deciphering temporal network changes is a core interest of molecular, ecological, neuropsychological and evolutionary biologists. Further, we will take a journey through different disciplines, such as social sciences, finance and computational gastronomy, to present commonalities and differences in how networks change and can be analysed. Finally, we envision how borrowing ideas from these disciplines could enrich the future of life science research.

Transcriptomic Responses to Water Deficit and Nematode Infection in Mycorrhizal Tomato Roots

Climate changes include the intensification of drought in many parts of the world, increasing its frequency, severity and duration. However, under natural conditions, environmental stresses do not occur alone, and, in addition, more stressed plants may become more susceptible to attacks by pests and pathogens. Studies on the impact of the arbuscular mycorrhizal (AM) symbiosis on tomato response to water deficit showed that several drought-responsive genes are differentially regulated in AM-colonized tomato plants (roots and leaves) during water deficit. To date, global changes in mycorrhizal tomato root transcripts under water stress conditions have not been yet investigated. Here, changes in root transcriptome in the presence of an AM fungus, with or without water stress (WS) application, have been evaluated in a commercial tomato cultivar already investigated for the water stress response during AM symbiosis. Since root-knot nematodes (RKNs, Meloidogyne incognita) are obligate endoparasites and cause severe yield losses in tomato, the impact of the AM fungal colonization on RKN infection at 7 days post-inoculation was also evaluated. Results offer new information about the response to AM symbiosis, highlighting a functional redundancy for several tomato gene families, as well as on the tomato and fungal genes involved in WS response during symbiosis, underlying the role of the AM fungus. Changes in the expression of tomato genes related to nematode infection during AM symbiosis highlight a role of AM colonization in triggering defense responses against RKN in tomato. Overall, new datasets on the tomato response to an abiotic and biotic stress during AM symbiosis have been obtained, providing useful data for further researches.