Genetic Conservation of CBS Domain Containing Protein Family in Oryza Species and Their Association with Abiotic Stress Responses

Crop Wild Relatives (CWRs) form a comprehensive gene pool that can answer the queries related to plant domestication, speciation, and ecological adaptation. The genus ‘Oryza’ comprises about 27 species, of which two are cultivated, while the remaining are wild. Here, we have attempted to understand the conservation and diversification of the genes encoding Cystathionine β-synthase (CBS) domain-containing proteins (CDCPs) in domesticated and CWRs of rice. Few members of CDCPs were previously identified to be stress-responsive and associated with multiple stress tolerance in rice. Through genome-wide analysis of eleven rice genomes, we identified a total of 36 genes encoding CDCPs in O. longistaminata, 38 in O. glaberrima, 39 each in O. rufipogon, O. glumaepatula, O. brachyantha, O. punctata, and O. sativa subsp. japonica, 40 each in O. barthii and O. meridionalis, 41 in O. nivara, and 42 in O. sativa subsp. indica. Gene duplication analysis as well as non-synonymous and synonymous substitutions in the duplicated gene pairs indicated that this family is shaped majorly by the negative or purifying selection pressure through the long-term evolution process. We identified the presence of two additional hetero-domains, namely TerCH and CoatomerE (specifically in O. sativa subsp. indica), which were not reported previously in plant CDCPs. The in silico expression analysis revealed some of the members to be responsive to various abiotic stresses. Furthermore, the qRT-PCR based analysis identified some members to be highly inducive specifically in salt-tolerant genotype in response to salinity. The cis-regulatory element analysis predicted the presence of numerous stress as well as a few phytohormone-responsive elements in their promoter region. The data presented in this study would be helpful in the characterization of these CDCPs from rice, particularly in relation to abiotic stress tolerance.     

Characterisation of Indica Special Protein (ISP), a Marker Protein for the Differentiation of Oryza sativa Subspecies indica and japonica

Based on both morphological and physiological traits, Asian cultivated rice (Oryza sativa L.) can be classified into two distinct subspecies, indica and japonica. To better understand the differences between the two subspecies, a proteomic approach was used to profile proteins present in the yellow seedling stage of 10 indica and 10 japonica rice varieties. We report the discovery of a new protein, Indica Special Protein (ISP), which was only detected in yellow seedlings of indica varieties, and was absent from japonica varieties. Hence, ISP may represent a key gene for the differentiation of indica and japonica subspecies.     

Genome-Wide Identification and Analysis of the Metallothionein Genes in Oryza Genus

Metallothionein (MT) proteins are low molecular mass, cysteine-rich, and metal-binding proteins that play an important role in maintaining metal homeostasis and stress response. However, the evolutionary relationships and functional differentiation of MT in the Oryza genus remain unclear. Here we identified 53 MT genes from six Oryza genera, including O. sativa ssp. japonica, O. rufipogon, O. sativa ssp. indica, O. nivara, O. glumaepatula, and O. barthii. The MT genes were clustered into four groups based on phylogenetic analysis. MT genes are unevenly distributed on chromosomes; almost half of the MT genes were clustered on chromosome 12, which may result from a fragment duplication containing the MT genes on chromosome 12. Five pairs of segmental duplication events and ten pairs of tandem duplication events were found in the rice MT family. The Ka/Ks values of the fifteen duplicated MT genes indicated that the duplicated MT genes were under a strong negative selection during evolution. Next, combining the promoter activity assay with gene expression analysis revealed different expression patterns of MT genes. In addition, the expression of OsMT genes was induced under different stresses, including NaCl, CdCl₂, ABA, and MeJ treatments. Additionally, we found that OsMT genes were mainly located in chloroplasts. These results imply that OsMT genes play different roles in response to these stresses. All results provide important insights into the evolution of the MT gene family in the Oryza genus, and will be helpful to further study the function of MT genes.     

Evolutionary Significance of NHX Family and NHX1 in Salinity Stress Adaptation in the Genus Oryza

Rice (Oryza sativa), a staple crop for a substantial part of the world’s population, is highly sensitive to soil salinity; however, some wild Oryza relatives can survive in highly saline environments. Sodium/hydrogen antiporter (NHX) family members contribute to Na⁺ homeostasis in plants and play a major role in conferring salinity tolerance. In this study, we analyzed the evolution of NHX family members using phylogeny, conserved domains, tertiary structures, expression patterns, and physiology of cultivated and wild Oryza species to decipher the role of NHXs in salt tolerance in Oryza. Phylogenetic analysis showed that the NHX family can be classified into three subfamilies directly related to their subcellular localization: endomembrane, plasma membrane, and tonoplast (vacuolar subfamily, vNHX1). Phylogenetic and structural analysis showed that vNHX1s have evolved from streptophyte algae (e.g., Klebsormidium nitens) and are abundant and highly conserved in all major land plant lineages, including Oryza. Moreover, we showed that tissue tolerance is a crucial trait conferring tolerance to salinity in wild rice species. Higher Na⁺ accumulation and reduced Na⁺ effluxes in leaf mesophyll were observed in the salt-tolerant wild rice species O. alta, O. latifolia, and O. coarctata. Among the key genes affecting tissue tolerance, expression of NHX1 and SOS1/NHX7 exhibited significant correlation with salt tolerance among the rice species and cultivars. This study provides insights into the evolutionary origin of plant NHXs and their role in tissue tolerance of Oryza species and facilitates the inclusion of this trait during the development of salinity-tolerant rice cultivars.     

Improved Transformation and Regeneration of Indica Rice: Disruption of SUB1A as a Test Case via CRISPR-Cas9

Gene editing by use of clustered regularly interspaced short palindromic repeats (CRISPR) has become a powerful tool for crop improvement. However, a common bottleneck in the application of this approach to grain crops, including rice (Oryza sativa), is efficient vector delivery and calli regeneration, which can be hampered by genotype-dependent requirements for plant regeneration. Here, methods for Agrobacterium-mediated and biolistic transformation and regeneration of indica rice were optimized using CRISPR-Cas9 gene-editing of the submergence tolerance regulator SUBMERGENCE 1A-1 gene of the cultivar Ciherang-Sub1. Callus induction and plantlet regeneration methods were optimized for embryogenic calli derived from immature embryos and mature seed-derived calli. Optimized regeneration (95%) and maximal editing efficiency (100%) were obtained from the immature embryo-derived calli. Phenotyping of T₁ seeds derived from the edited T₀ plants under submergence stress demonstrated inferior phenotype compared to their controls, which phenotypically validates the disruption of SUB1A-1 function. The methods pave the way for rapid CRISPR-Cas9 gene editing of recalcitrant indica rice cultivars.     

Haplotype Analysis of BADH1 by Next-Generation Sequencing Reveals Association with Salt Tolerance in Rice during Domestication

Betaine aldehyde dehydrogenase 1 (BADH1), a paralog of the fragrance gene BADH2, is known to be associated with salt stress through the accumulation of synthesized glycine betaine (GB), which is involved in the response to abiotic stresses. Despite the unclear association between BADH1 and salt stress, we observed the responses of eight phenotypic characteristics (germination percentage (GP), germination energy (GE), germination index (GI), mean germination time (MGT), germination rate (GR), shoot length (SL), root length (RL), and total dry weight (TDW)) to salt stress during the germination stage of 475 rice accessions to investigate their association with BADH1 haplotypes. We found a total of 116 SNPs and 77 InDels in the whole BADH1 gene region, representing 39 haplotypes. Twenty-nine haplotypes representing 27 mutated alleles (two InDels and 25 SNPs) were highly (p < 0.05) associated with salt stress, including the five SNPs that have been previously reported to be associated with salt tolerance. We observed three predominant haplotypes associated with salt tolerance, Hap_2, Hap_18, and Hap_23, which were Indica specific, indicating a comparatively high number of rice accessions among the associated haplotypes. Eight plant parameters (phenotypes) also showed clear responses to salt stress, and except for MGT (mean germination time), all were positively correlated with each other. Different signatures of domestication for BADH1 were detected in cultivated rice by identifying the highest and lowest Tajima’s D values of two major cultivated ecotypes (Temperate Japonica and Indica). Our findings on these significant associations and BADH1 evolution to plant traits can be useful for future research development related to its gene expression.     

Programmed Editing of Rice (Oryza sativa L.) OsSPL16 Gene Using CRISPR/Cas9 Improves Grain Yield by Modulating the Expression of Pyruvate Enzymes and Cell Cycle Proteins

Rice (Oryza sativa L.) is one of the major crops in the world and significant increase in grain yield is constant demand for breeders to meet the needs of a rapidly growing population. The size of grains is one of major components determining rice yield and a vital trait for domestication and breeding. To increase the grain size in rice, OsSPL16/qGW8 was mutagenized through CRISPR/Cas9, and proteomic analysis was performed to reveal variations triggered by mutations. More specifically, mutants were generated with two separate guide RNAs targeting recognition sites on opposite strands and genomic insertions and deletions were characterized. Mutations followed Mendelian inheritance and homozygous and heterozygous mutants lacking any T-DNA and off-target effects were screened. The mutant lines showed a significant increase in grain yield without any change in other agronomic traits in T₀, T₁, and T₂ generations. Proteomic screening found a total of 44 differentially expressed proteins (DEPs), out of which 33 and 11 were up and downregulated, respectively. Most of the DEPs related to pyruvate kinase, pyruvate dehydrogenase, and cell division and proliferation were upregulated in the mutant plants. Pathway analysis revealed that DEPs were enriched in the biosynthesis of secondary metabolites, pyruvate metabolism, glycolysis/gluconeogenesis, carbon metabolism, ubiquinone and other terpenoid-quinone biosynthesis, and citrate cycle. Gene Ontology (GO) analysis presented that most of the DEPs were involved in the pyruvate metabolic process and pyruvate dehydrogenase complex. Proteins related to pyruvate dehydrogenase E1 component subunit alpha-1 displayed higher interaction in the protein-protein interaction (PPI) network. Thus, the overall results revealed that CRISPR/Cas9-guided OsSPL16 mutations have the potential to boost the grain yield of rice. Additionally, global proteome analysis has broad applications for discovering molecular components and dynamic regulation underlying the targeted gene mutations.     

Conserved Imprinted Genes between Intra-Subspecies and Inter-Subspecies Are Involved in Energy Metabolism and Seed Development in Rice

Genomic imprinting is an epigenetic phenomenon in which a subset of genes express dependent on the origin of their parents. In plants, it is unclear whether imprinted genes are conserved between subspecies in rice. Here we identified imprinted genes from embryo and endosperm 5–7 days after pollination from three pairs of reciprocal hybrids, including inter-subspecies, japonica intra-subspecies, and indica intra-subspecies reciprocal hybrids. A total of 914 imprinted genes, including 546 in inter-subspecies hybrids, 211 in japonica intra-subspecies hybrids, and 286 in indica intra-subspecies hybrids. In general, the number of maternally expressed genes (MEGs) is more than paternally expressed genes (PEGs). Moreover, imprinted genes tend to be in mini clusters. The number of shared genes by R9N (reciprocal crosses between 9311 and Nipponbare) and R9Z (reciprocal crosses between 9311 and Zhenshan 97), R9N and RZN (reciprocal crosses between Zhonghua11 and Nipponbare), R9Z and RZN was 72, 46, and 16. These genes frequently involved in energy metabolism and seed development. Five imprinted genes (Os01g0151700, Os07g0103100, Os10g0340600, Os11g0679700, and Os12g0632800) are commonly detected in all three pairs of reciprocal hybrids and were validated by RT-PCR sequencing. Gene editing of two imprinted genes revealed that both genes conferred grain filling. Moreover, 15 and 27 imprinted genes with diverse functions in rice were shared with Arabidopsis and maize, respectively. This study provided valuable resources for identification of imprinting genes in rice or even in cereals.     

γ-Oryzanols of North American Wild Rice (Zizania palustris)

γ-Oryzanol, a natural mixture of ferulic acid esters of triterpene alcohols and sterols, are an important bioactive components present in rice bran oil. In light of the recent increase in the popularity of wild rice among consumers, and the possibility of a direct relationship between γ-oryzanol composition and its bioactivity, the oryzanol profile of major wild rice (Zizania palustris) grown in North America was studied and compared to regular brown rice (Oryza sativa L.). A total of twenty-three γ-oryzanol components were separated, identified and quantified by HPLC coupled to an Orbitrap MS. The distribution of individual γ-oryzanols was similar for all the wild rice but significantly different from those of the regular brown rice. Unlike in the regular brown rice, a significant amount of steryl caffeate and cinnamate were found in the wild rice samples. Generally, the amounts of γ-oryzanol in the wild rice were higher compared to the regular brown rice, 1,352 vs. 688 μg/g. The results from this study showed that wild rice had a more diverse γ-oryzanol composition and the higher amounts compared to the regular brown rice.     

Effect of steam distillate extracts of selected resistant cultivated and wild rices on behavior of leaffolderCnaphalocrocis medinalis (Guenée) (Lepidoptera: Pyralidae)

Rice plant volatiles extracted as steam distillates significantly affected the behavior of the rice leaffolder,Cnaphalocrocis medinalis (Guenée). Leaffolder moths laid significantly fewer eggs on TN1 rice plants treated with extracts of resistant wild species of rice,Oryza officinalis andOryza punctata, than on TN1 plants treated with extracts of resistant cultivated rices. Extracts ofO. officinalis andO. punctata adversely affected egg hatchability.O. officinalis andO. punctata extracts were more toxic to first-instar larvae than extracts of other resistant varieties. Leaf area consumed by the larvae was reduced on TN1 plants treated withO. officinalis andO. punctata extracts compared with that on plants treated with extracts of resistant cultivated rice varieties.     

Role of Long Noncoding RNAs ZlMSTRG.11348 and UeMSTRG.02678 in Temperature-Dependent Culm Swelling in Zizania latifolia

Temperature influences the physiological processes and ecology of both hosts and endophytes; however, it remains unclear how long noncoding RNAs (lncRNAs) modulate the consequences of temperature-dependent changes in host–pathogen interactions. To explore the role of lncRNAs in culm gall formation induced by the smut fungus Ustilago esculenta in Zizania latifolia, we employed RNA sequencing to identify lncRNAs and their potential cis-targets in Z. latifolia and U. esculenta under different temperatures. In Z. latifolia and U. esculenta, we identified 3194 and 173 lncRNAs as well as 126 and four potential target genes for differentially expressed lncRNAs, respectively. Further function and expression analysis revealed that lncRNA ZlMSTRG.11348 regulates amino acid metabolism in Z. latifolia and lncRNA UeMSTRG.02678 regulates amino acid transport in U. esculenta. The plant defence response was also found to be regulated by lncRNAs and suppressed in Z. latifolia infected with U. esculenta grown at 25 °C, which may result from the expression of effector genes in U. esculenta. Moreover, in Z. latifolia infected with U. esculenta, the expression of genes related to phytohormones was altered under different temperatures. Our results demonstrate that lncRNAs are important components of the regulatory networks in plant-microbe-environment interactions, and may play a part in regulating culm swelling in Z. latifolia plants.     

Genomic Architecture of Phenotypic Plasticity in Response to Water Stress in Tetraploid Wheat

Phenotypic plasticity is one of the main mechanisms of adaptation to abiotic stresses via changes in critical developmental stages. Altering flowering phenology is a key evolutionary strategy of plant adaptation to abiotic stresses, to achieve the maximum possible reproduction. The current study is the first to apply the linear regression residuals as drought plasticity scores while considering the variation in flowering phenology and traits under non-stress conditions. We characterized the genomic architecture of 17 complex traits and their drought plasticity scores for quantitative trait loci (QTL) mapping, using a mapping population derived from a cross between durum wheat (Triticum turgidum ssp. durum) and wild emmer wheat (T. turgidum ssp. dicoccoides). We identified 79 QTLs affected observed traits and their plasticity scores, of which 33 reflected plasticity in response to water stress and exhibited epistatic interactions and/or pleiotropy between the observed and plasticity traits. Vrn-B3 (TaTF1) residing within an interval of a major drought-escape QTL was proposed as a candidate gene. The favorable alleles for most of the plasticity QTLs were contributed by wild emmer wheat, demonstrating its high potential for wheat improvement. Our study presents a new approach for the quantification of plant adaptation to various stresses and provides new insights into the genetic basis of wheat complex traits under water-deficit stress.     

Glycoalkaloids of Wild and Cultivated Solanum: Effects on Specialist and Generalist Insect Herbivores

Plant domestication by selective breeding may reduce plant chemical defense in favor of growth. However, few studies have simultaneously studied the defensive chemistry of cultivated plants and their wild congeners in connection to herbivore susceptibility. We compared the constitutive glycoalkaloids (GAs) of cultivated potato, Solanum tuberosum, and a wild congener, S. commersonii, by liquid chromatography coupled to mass spectrometry. We also determined the major herbivores present on the two species in field plots, and tested their preference for the plants and their isolated GAs in two-choice bioassays. Solanum commersonii had a different GA profile and higher concentrations than S. tuberosum. In the field, S. tuberosum was mostly attacked by the generalist aphids Myzus persicae and Macrosiphum euphorbiae, and by the specialist flea beetle Epitrix argentinensis. In contrast, the most common herbivore on S. commersonii was the specialist sawfly Tequus sp. Defoliation levels were higher on the wild species, probably due to the chewing feeding behavior of Tequus sp. As seen in the field, M. persicae and E. argentinensis preferred leaf disks of the cultivated plant, while Tequus sp. preferred those of the wild one. Congruently, GAs from S. commersonii were avoided by M. persicae and preferred by Tequus sp. The potato aphid performed well on both species and was not deterred by S. commersonii GAs. These observations suggest that different GA profiles explain the feeding preferences of the different herbivores, and that domestication has altered the defensive capacity of S. tuberosum. However, the wild relative is still subject to severe defoliation by a specialist herbivore that may cue on the GAs.