RNA-seq reveals the involvement of key genes for aerobic adaptation in rice

Genetic mapping of regions associated with root system architecture in rice using MutMap QTL-seq

The root system architecture is an important complex trait in rice. With changing climatic conditions and soil nutrient deficiencies, there is an immediate need to breed nutrient-use-efficient rice varieties with robust root system architectural (RSA) traits. To map the genomic regions associated with crucial component traits of RSA viz. root length and root volume, a biparental F2 mapping population was developed using TI-128, an Ethyl Methane Sulphonate (EMS) mutant of a mega variety BPT-5204 having high root length (RL) and root volume (RV) with wild type BPT-5204. Extreme bulks having high RL and RV and low RL and RV were the whole genome re-sequenced along with parents. Genetic mapping using the MutMap QTL-Seq approach elucidated two genomic intervals on Chr.12 (3.14-3.74 Mb, 18.11-20.85 Mb), and on Chr.2 (23.18-23.68 Mb) as potential regions associated with both RL and RV. The Kompetitive Allele Specific PCR (KASP) assays for SNPs with delta SNP index near 1 were associated with higher RL and RV in the panel of sixty-two genotypes varying in root length and volume. The KASP_SNPs viz. Chr12_S4 (C→T; Chr12:3243938), located in the 3' UTR region of LOC_Os12g06670 encoding a protein kinase domain-containing protein and Chr2_S6 (C→T; Chr2:23181622) present upstream in the regulator of chromosomal condensation protein LOC_Os2g38350. Validation of these genes using qRT-PCR and in-silico studies using various online tools and databases revealed higher expression in TI-128 as compared to BPT- 5204 at the seedling and panicle initiation stages implying the functional role in enhancing RL and RV.

Upland rice genomic signatures of adaptation to drought resistance and navigation to molecular design breeding

Upland rice is a distinctive drought-aerobic ecotype of cultivated rice highly resistant to drought stress. However, the genetic and genomic basis for the drought-aerobic adaptation of upland rice remains largely unclear due to the lack of genomic resources. In this study, we identified 25 typical upland rice accessions and assembled a high-quality genome of one of the typical upland rice varieties, IRAT109, comprising 384 Mb with a contig N50 of 19.6 Mb. Phylogenetic analysis revealed upland and lowland rice have distinct ecotype differentiation within the japonica subgroup. Comparative genomic analyses revealed that adaptive differentiation of lowland and upland rice is likely attributable to the natural variation of many genes in promoter regions, formation of specific genes in upland rice, and expansion of gene families. We revealed differentiated gene expression patterns in the leaves and roots of the two ecotypes and found that lignin synthesis mediated by the phenylpropane pathway plays an important role in the adaptive differentiation of upland and lowland rice. We identified 28 selective sweeps that occurred during domestication and validated that the qRT9 gene in selective regions can positively regulate drought resistance in rice. Eighty key genes closely associated with drought resistance were appraised for their appreciable potential in drought resistance breeding. Our study enhances the understanding of the adaptation of upland rice and provides a genome navigation map of drought resistance breeding, which will facilitate the breeding of drought-resistant rice and the "blue revolution" in agriculture.

Long non-coding RNA-mediated epigenetic response for abiotic stress tolerance in plants

Plants perceive environmental fluctuations as stress and confront several stresses throughout their life cycle individually or in combination. Plants have evolved their sensing and signaling mechanisms to perceive and respond to a variety of stresses. Epigenetic regulation plays a critical role in the regulation of genes, spatiotemporal expression of genes under stress conditions and imparts a stress memory to encounter future stress responses. It is quintessential to integrate our understanding of genetics and epigenetics to maintain plant fitness, achieve desired genetic gains with no trade-offs, and durable long-term stress tolerance. The long non-coding RNA >200 nts having no coding potential (or very low) play several roles in epigenetic memory, contributing to the regulation of gene expression and the maintenance of cellular identity which include chromatin remodeling, imprinting (dosage compensation), stable silencing, facilitating nuclear organization, regulation of enhancer-promoter interactions, response to environmental signals and epigenetic switching. The lncRNAs are involved in a myriad of stress responses by activation or repression of target genes and hence are potential candidates for deploying in climate-resilient breeding programs. This review puts forward the significant roles of long non-coding RNA as an epigenetic response during abiotic stresses in plants and the prospects of deploying lncRNAs for designing climate-resilient plants.

Combined strategy employing MutMap and RNA-seq reveals genomic regions and genes associated with complete panicle exsertion in rice

Complete panicle exsertion (CPE) in rice is an important determinant of yield and a desirable trait in breeding. However, the genetic basis of CPE in rice still remains to be completely characterized. An ethyl methane sulfonate (EMS) mutant line of an elite cultivar Samba Mahsuri (BPT 5204), displaying stable and consistent CPE, was identified and named as CPE-110. MutMap and RNA-seq were deployed for unraveling the genomic regions, genes, and markers associated with CPE. Two major genomic intervals, on chromosome 8 (25668481-25750456) and on chromosome 11 (20147154-20190400), were identified to be linked to CPE through MutMap. A non-synonymous SNP (G/A; Chr8:25683828) in the gene LOC_Os08g40570 encoding pyridoxamine 5'-phosphate oxidase with the SNP index 1 was converted to Kompetitive allele-specific PCR (KASP) marker. This SNP (KASP 8-1) exhibited significant association with CPE and further validated through assay in the F2 mapping population, released varieties and CPE exhibiting BPT 5204 mutant lines. RNA-seq of the flag leaves at the booting stage, 1100 genes were upregulated and 1305 downregulated differentially in CPE-110 and BPT 5204. Metabolic pathway analysis indicated an enrichment of genes involved in photosynthesis, glyoxylate, dicarboxylate, porphyrin, pyruvate, chlorophyll, carotenoid, and carbon metabolism. Further molecular and functional studies of the candidate genes could reveal the mechanistic aspects of CPE.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-023-01412-1.

Silicate solubilizing and plant growth promoting bacteria interact with biogenic silica to impart heat stress tolerance in rice by modulating physiology and gene expression

Heat stress caused due to increasing warming climate has become a severe threat to global food production including rice. Silicon plays a major role in improving growth and productivity of rice by aiding in alleviating heat stress in rice. Soil silicon is only sparingly available to the crops can be made available by silicate solubilizing and plant-growth-promoting bacteria that possess the capacity to solubilize insoluble silicates can increase the availability of soluble silicates in the soil. In addition, plant growth promoting bacteria are known to enhance the tolerance to abiotic stresses of plants, by affecting the biochemical and physiological characteristics of plants. The present study is intended to understand the role of beneficial bacteria viz. Rhizobium sp. IIRR N1 a silicate solublizer and Gluconacetobacter diazotrophicus, a plant growth promoting bacteria and their interaction with insoluble silicate sources on morpho-physiological and molecular attributes of rice (Oryza sativa L.) seedlings after exposure to heat stress in a controlled hydroponic system. Joint inoculation of silicates and both the bacteria increased silicon content in rice tissue, root and shoot biomass, significantly increased the antioxidant enzyme activities (viz. superoxidase dismutase, catalase and ascorbate peroxidase) compared to other treatments with sole application of either silicon or bacteria. The physiological traits (viz. chlorophyll content, relative water content) were also found to be significantly enhanced in presence of silicates and both the bacteria after exposure to heat stress conditions. Expression profiling of shoot and root tissues of rice seedlings revealed that seedlings grown in the presence of silicates and both the bacteria exhibited higher expression of heat shock proteins (HSPs viz., OsHsp90, OsHsp100 and 60 kDa chaperonin), hormone-related genes (OsIAA6) and silicon transporters (OsLsi1 and OsLsi2) as compared to seedlings treated with either silicates or with the bacteria alone. The results thus reveal the interactive effect of combined application of silicates along with bacteria Rhizobium sp. IIRR N1, G. diazotrophicus inoculation not only led to augmented silicon uptake by rice seedlings but also influenced the plant biomass and elicited higher expression of HSPs, hormone-related and silicon transporter genes leading to improved tolerance of seedling to heat stress.

Optimization and application of non-native Phragmites australis transcriptome assemblies

Phragmites australis (common reed) has a cosmopolitan distribution and has been suggested as a model organism for the study of invasive plant species. In North America, the non-native subspecies (ssp. australis) is widely distributed across the contiguous 48 states in the United States and large parts of Canada. Even though millions of dollars are spent annually on Phragmites management, insufficient knowledge of P. australis impeded the efficiency of management. To solve this problem, transcriptomic information generated from multiple types of tissue could be a valuable resource for future studies. Here, we constructed forty-nine P. australis transcriptomes assemblies via different assembly tools and multiple parameter settings. The optimal transcriptome assembly for functional annotation and downstream analyses was selected among these transcriptome assemblies by comprehensive assessments. For a total of 422,589 transcripts assembled in this transcriptome assembly, 319,046 transcripts (75.5%) have at least one functional annotation. Within the transcriptome assembly, we further identified 1,495 transcripts showing tissue-specific expression pattern, 10,828 putative transcription factors, and 72,165 candidates for simple sequence repeats markers. The identification and analyses of predicted transcripts related to herbicide- and salinity-resistant genes were shown as two applications of the transcriptomic information to facilitate further research on P. australis. Transcriptome assembly and selection would be important for the transcriptome annotation. With this optimal transcriptome assembly and all relative information from downstream analyses, we have helped to establish foundations for future studies on the mechanisms underlying the invasiveness of non-native P. australis subspecies.

Genomics and transcriptomics to protect rice (Oryza sativa. L.) from abiotic stressors: -pathways to achieving zero hunger

More over half of the world's population depends on rice as a major food crop. Rice (Oryza sativa L.) is vulnerable to abiotic challenges including drought, cold, and salinity since it grown in semi-aquatic, tropical, or subtropical settings. Abiotic stress resistance has bred into rice plants since the earliest rice cultivation techniques. Prior to the discovery of the genome, abiotic stress-related genes were identified using forward genetic methods, and abiotic stress-tolerant lines have developed using traditional breeding methods. Dynamic transcriptome expression represents the degree of gene expression in a specific cell, tissue, or organ of an individual organism at a specific point in its growth and development. Transcriptomics can reveal the expression at the entire genome level during stressful conditions from the entire transcriptional level, which can be helpful in understanding the intricate regulatory network relating to the stress tolerance and adaptability of plants. Rice (Oryza sativa L.) gene families found comparatively using the reference genome sequences of other plant species, allowing for genome-wide identification. Transcriptomics via gene expression profiling which have recently dominated by RNA-seq complements genomic techniques. The identification of numerous important qtl,s genes, promoter elements, transcription factors and miRNAs involved in rice response to abiotic stress was made possible by all of these genomic and transcriptomic techniques. The use of several genomes and transcriptome methodologies to comprehend rice (Oryza sativa, L.) ability to withstand abiotic stress have been discussed in this review.

The mechanism of Megalobrama amblycephala muscle injury repair based on RNA-seq

Skeletal muscle myogenesis and injury-induced muscle regeneration contribute to muscle formation. Skeletal muscle stem cells, termed satellite cells (SCs), proliferate to repair injured muscle. To identify the molecular mechanism of regeneration after muscle injury as well as the genes related to muscle development in fish, in this study, the immunohistochemistry and the high-throughput RNA sequencing (RNA-seq) analysis were performed after Megalobrama amblycephala muscle was injured by needle stab. The results showed that paired box7-positive (Pax7⁺) SCs increased, and peaked at 96 to 144 h-post injury (hpi). The 6729 differentially expressed genes (DEGs), including 2125 up-regulated and 4604 down-regulated genes were found. GO terms significantly enriched by DEGs contained intercellular connections, signaling transduction and enzyme activity. KEGG enrichment analysis showed that most of the pathways were related to immunity, metabolism and cells related molecules, including actin skeleton regulation, Epstein Barr virus infection and plaque adhesion. The WGCNA results revealed that actin cytoskeleton and lipid metabolism related genes probably played crucial roles during repair after muscle injury. Collectively, all these results will provide new insights into the molecular mechanisms underlying muscle injury repair of fish.

RNA-Seq Analysis of Magnaporthe grisea Transcriptome Reveals the High Potential of ZnO Nanoparticles as a Nanofungicide

Magnaporthe grisea is one of the most destructive pathogen that encounters a challenge to rice production around the worldwide. The unique properties of ZnO nanoparticles (NPs), have high attractiveness as nanofungicide. In the present study, the response of fungi to ZnO NPs was evaluated using RNA sequencing (RNA-seq). Two different aligners (STAR and Hisat2) were used for aligning the clean reads, and the DEseq2 package was used to identify the differentially expressed genes (DEGs). In total, 1,438 and 761 fungal genes were significantly up- and down-regulated in response to ZnO NPs, respectively. The DEGs were subjected to functional enrichment analysis to identify significantly enriched biological pathways. Functional enrichment analysis revealed that "cell membrane components," "ion (calcium) transmembrane transporter activity," "steroid biosynthesis pathway" and "catalytic activity" were the contributed terms to fungal response mechanisms. The genes involved in aflatoxin efflux pumps and ribosome maturation were among the genes showing significant up- and down-regulation after ZnO NPs application. To confirm the obtained RNA-seq results, the expression of six randomly selected genes were evaluated using q-RT-PCR. Overall, the RNA-seq results suggest that ZnO NPs primarily act on the fungal cell membrane, but accumulation of ROS inside the cell induces oxidative stress, the fungal catalytic system is disrupted, resulting into the inhibition of ROS scavenging and eventually, to the death of fungal cells. Our findings provide novel insights into the effect of ZnO NPs as a promising nanofungicide for effective control of rice blast disease.

Transcriptome Analysis of the Cf-13-Mediated Hypersensitive Response of Tomato to Cladosporium fulvum Infection

Tomato leaf mold disease caused by Cladosporium fulvum (C. fulvum) is one of the most common diseases affecting greenhouse tomato production. Cf proteins can recognize corresponding AVR proteins produced by C. fulvum, and Cf genes are associated with leaf mold resistance. Given that there are many physiological races of C. fulvum and that these races rapidly mutate, resistance to common Cf genes (such as Cf-2, Cf-4, Cf-5, and Cf-9) has decreased. In the field, Ont7813 plants (carrying the Cf-13 gene) show effective resistance to C. fulvum; thus, these plants could be used as new, disease-resistant materials. To explore the mechanism of the Cf-13-mediated resistance response, transcriptome sequencing was performed on three replicates each of Ont7813 (Cf-13) and Moneymaker (MM; carrying the Cf-0 gene) at 0, 9, and 15 days after inoculation (dai) for a total of 18 samples. In total, 943 genes were differentially expressed, specifically in the Ont7813 response process as compared to the Moneymaker response process. Gene ontology (GO) classification of these 943 differentially expressed genes (DEGs) showed that GO terms, including "hydrogen peroxide metabolic process (GO_Process)", "secondary active transmembrane transporter activity (GO_Function)", and "mismatch repair complex (GO_Component)", which were the same as 11 other GO terms, were significantly enriched. An analysis of the Kyoto Encyclopedia of Genes and Genomes (KEGG) revealed that many key regulatory genes of the Cf-13-mediated resistance response processes were involved in the "plant hormone signal transduction" pathway, the "plant-pathogen interaction" pathway, and the "MAPK signaling pathway-plant" pathway. Moreover, during C. fulvum infection, jasmonic acid (JA) and salicylic acid (SA) contents significantly increased in Ont7813 at the early stage. These results lay a vital foundation for further understanding the molecular mechanism of the Cf-13 gene in response to C. fulvum infection.

Insights into the Host-Pathogen Interaction Pathways through RNA-Seq Analysis of Lens culinaris Medik. in Response to Rhizoctonia bataticola Infection

Dry root rot (Rhizoctonia bataticola) is an important disease of lentils (Lens culinaris Medik.).To gain an insight into the molecular aspects of host-pathogen interactions, the RNA-seq approach was used in lentils following inoculation with R.bataticola. The RNA-Seq has generated >450 million high-quality reads (HQRs) and nearly 96.97% were properly aligned to the reference genome. Very high similarity in FPKM (fragments per kilobase of exon per million mapped fragments) values (R > 0.9) among biological replicates showed the consistency of the RNA-Seq results. The study revealed various DEGs (differentially expressed genes) that were associated with changes in phenolic compounds, transcription factors (TFs), antioxidants, receptor kinases, hormone signals which corresponded to the cell wall modification enzymes, defense-related metabolites, and jasmonic acid (JA)/ethylene (ET) pathways. Gene ontology (GO) categorization also showed similar kinds of significantly enriched similar GO terms. Interestingly, of the total unigenes (42,606), 12,648 got assembled and showed significant hit with Rhizoctonia species. String analysis also revealed the role of various disease responsive proteins viz., LRR family proteins, LRR-RLKs, protein kinases, etc. in the host-pathogen interaction. Insilico validation analysis was performed using Genevestigator^® and DEGs belonging to six major defense-response groups viz., defense-related enzymes, disease responsive genes, hormones, kinases, PR (pathogenesis related) proteins, and TFs were validated. For the first time some key miRNA targets viz. miR156, miR159, miR167, miR169, and miR482 were identified from the studied transcriptome, which may have some vital role in Rhizoctonia-based responses in lentils. The study has revealed the molecular mechanisms of the lentil/R.bataticola interactions and also provided a theoretical approach for the development of lentil genotypes resistant to R.bataticola.

Stress-Induced Changes in Alternative Splicing Landscape in Rice: Functional Significance of Splice Isoforms in Stress Tolerance

Improvements in yield and quality of rice are crucial for global food security. However, global rice production is substantially hindered by various biotic and abiotic stresses. Making further improvements in rice yield is a major challenge to the rice research community, which can be accomplished through developing abiotic stress-resilient rice varieties and engineering durable agrochemical-independent pathogen resistance in high-yielding elite rice varieties. This, in turn, needs increased understanding of the mechanisms by which stresses affect rice growth and development. Alternative splicing (AS), a post-transcriptional gene regulatory mechanism, allows rapid changes in the transcriptome and can generate novel regulatory mechanisms to confer plasticity to plant growth and development. Mounting evidence indicates that AS has a prominent role in regulating rice growth and development under stress conditions. Several regulatory and structural genes and splicing factors of rice undergo different types of stress-induced AS events, and the functional significance of some of them in stress tolerance has been defined. Both rice and its pathogens use this complex regulatory mechanism to devise strategies against each other. This review covers the current understanding and evidence for the involvement of AS in biotic and abiotic stress-responsive genes, and its relevance to rice growth and development. Furthermore, we discuss implications of AS for the virulence of different rice pathogens and highlight the areas of further research and potential future avenues to develop climate-smart and disease-resistant rice varieties.

Probing into the Effects of Grapevine Leafroll-Associated Viruses on the Physiology, Fruit Quality and Gene Expression of Grapes

Grapevine leafroll is one of the most widespread and highly destructive grapevine diseases that is responsible for great economic losses to the grape and wine industries throughout the world. Six distinct viruses have been implicated in this disease complex. They belong to three genera, all in the family Closteroviridae. For the sake of convenience, these viruses are named as grapevine leafroll-associated viruses (GLRaV-1, -2, -3, -4, -7, and -13). However, their etiological role in the disease has yet to be established. Furthermore, how infections with each GLRaV induce the characteristic disease symptoms remains unresolved. Here, we first provide a brief overview on each of these GLRaVs with a focus on genome structure, expression strategies and gene functions, where available. We then provide a review on the effects of GLRaV infection on the physiology, fruit quality, fruit chemical composition, and gene expression of grapevine based on the limited information so far reported in the literature. We outline key methodologies that have been used to study how GLRaV infections alter gene expression in the grapevine host at the transcriptomic level. Finally, we present a working model as an initial attempt to explain how infections with GLRaVs lead to the characteristic symptoms of grapevine leafroll disease: leaf discoloration and downward rolling. It is our hope that this review will serve as a starting point for grapevine virology and the related research community to tackle this vastly important and yet virtually uncharted territory in virus-host interactions involving woody and perennial fruit crops.

Diverse Roles of MAX1 Homologues in Rice

Cytochrome P450 enzymes encoded by MORE AXILLARY GROWTH1 (MAX1)-like genes produce most of the structural diversity of strigolactones during the final steps of strigolactone biosynthesis. The diverse copies of MAX1 in Oryza sativa provide a resource to investigate why plants produce such a wide range of strigolactones. Here we performed in silico analyses of transcription factors and microRNAs that may regulate each rice MAX1, and compared the results with available data about MAX1 expression profiles and genes co-expressed with MAX1 genes. Data suggest that distinct mechanisms regulate the expression of each MAX1. Moreover, there may be novel functions for MAX1 homologues, such as the regulation of flower development or responses to heavy metals. In addition, individual MAX1s could be involved in specific functions, such as the regulation of seed development or wax synthesis in rice. Our analysis reveals potential new avenues of strigolactone research that may otherwise not be obvious.

Identification and Validation of Aerobic Adaptation QTLs in Upland Rice

The aerobic adaptation of upland rice is considered as the key genetic difference between upland rice and lowland rice. Genetic dissection of the aerobic adaptation is important as the basis for improving drought tolerance and terrestrial adaptation by using the upland rice. We raised BC₁-BC₃ introgression lines (ILs) in lowland rice Minghui 63 (MH63) background. The QTLs of yield and yield-related traits were detected based on ILs under the aerobic and lowland environments, and then the yield-related QTLs were identified in a backcrossed inbred population of BC₄F₅ under aerobic condition. We further verified phenotypes of QTL near-isogenic lines. Finally, three QTLs responsible for increasing yield in aerobic environment were detected by multiple locations and generations, which were designated as qAER1, qAER3, and qAER9 (QTL of aerobic adaptation). The qAER1 and qAER9 were fine-mapped. We found that qAER1 and qAER9 controlled plant height and heading date, respectively; while both of them increased yields simultaneously by suitable plant height and heading date without delay in the aerobic environment. The phenotypic differences between lowland rice and upland rice in the aerobic environment further supported the above results. We pyramided the two QTLs as corresponding molecular modules in the irrigated lowland rice MH63 background, and successfully developed a new upland rice variety named as Zhongkexilu 2. This study will lay the foundation for using aerobic adaptation QTLs in rice breeding programs and for further cloning the key genes involved in aerobic adaptation.

Comparative transcriptome analysis reveals the response mechanism of Cf-16-mediated resistance to Cladosporium fulvum infection in tomato

BACKGROUND

Leaf mold disease caused by Cladosporium fulvum is a serious threat affecting the global production of tomato. Cf genes are associated with leaf mold resistance, including Cf-16, which confers effective resistance to leaf mold in tomato. However, the molecular mechanism of the Cf-16-mediated resistance response is largely unknown.

RESULTS

We performed a comparative transcriptome analysis of C. fulvum-resistant (cv. Ontario7816) and C. fulvum-susceptible (cv. Moneymaker) tomato cultivars to identify differentially expressed genes (DEGs) at 4 and 8 days post inoculation (dpi) with C. fulvum. In total, 1588 and 939 more DEGs were found in Cf-16 tomato than in Moneymaker at 4 and 8 dpi, respectively. Additionally, 1350 DEGs were shared between the 4- and 8-dpi Cf-16 groups, suggesting the existence of common core DEGs in response to C. fulvum infection. The up-regulated DEGs in Cf-16 tomato were primarily associated with defense processes and phytohormone signaling, including salicylic acid (SA) and jasmonic acid (JA). Moreover, SA and JA levels were significantly increased in Cf-16 tomato at the early stages of C. fulvum infection. Contrary to the previous study, the number of up-regulated genes in Cf-16 compared to Cf-10 and Cf-12 tomatoes was significantly higher at the early stages of C. fulvum infection.

CONCLUSION

Our results provide new insight into the Cf-mediated mechanism of resistance to C. fulvum, especially the unique characteristics of Cf-16 tomato in response to this fungus.