Characterization of the transcriptional divergence between the subspecies of cultivated rice (Oryza sativa)

Inferring Historical Introgression with Deep Learning

Resolving phylogenetic relationships among taxa remains a challenge in the era of big data due to the presence of genetic admixture in a wide range of organisms. Rapidly developing sequencing technologies and statistical tests enable evolutionary relationships to be disentangled at a genome-wide level, yet many of these tests are computationally intensive and rely on phased genotypes, large sample sizes, restricted phylogenetic topologies, or hypothesis testing. To overcome these difficulties, we developed a deep learning-based approach, named ERICA, for inferring genome-wide evolutionary relationships and local introgressed regions from sequence data. ERICA accepts sequence alignments of both population genomic data and multiple genome assemblies, and efficiently identifies discordant genealogy patterns and exchanged regions across genomes when compared with other methods. We further tested ERICA using real population genomic data from Heliconius butterflies that have undergone adaptive radiation and frequent hybridization. Finally, we applied ERICA to characterize hybridization and introgression in wild and cultivated rice, revealing the important role of introgression in rice domestication and adaptation. Taken together, our findings demonstrate that ERICA provides an effective method for teasing apart evolutionary relationships using whole genome data, which can ultimately facilitate evolutionary studies on hybridization and introgression.

Genomic variations combined with epigenetic modifications rewire open chromatin in rice

Cis-regulatory elements (CREs) fine-tune gene transcription in eukaryotes. CREs with sequence variations play vital roles in driving plant or crop domestication. However, how global sequence and structural variations (SVs) are responsible for multilevel changes between indica and japonica rice (Oryza sativa) is still not fully elucidated. To address this, we conducted multiomic studies using MNase hypersensitivity sequencing (MH-seq) in combination with RNA sequencing (RNA-seq), chromatin immunoprecipitation sequencing (ChIP-seq), and bisulfite sequencing (BS-seq) between the japonica rice variety Nipponbare (NIP) and indica rice variety 93-11. We found that differential MNase hypersensitive sites (MHSs) exhibited some distinct intrinsic genomic sequence features between NIP and 93-11. Notably, through MHS-genome-wide association studies (GWAS) integration, we found that key sequence variations may be associated with differences of agronomic traits between NIP and 93-11, which is partly achieved by MHSs harboring CREs. In addition, SV-derived differential MHSs caused by transposable element (TE) insertion, especially by noncommon TEs among rice varieties, were associated with genes with distinct functions, indicating that TE-driven gene neo- or subfunctionalization is mediated by changes of chromatin openness. This study thus provides insights into how sequence and genomic SVs control agronomic traits of NIP and 93-11; it also provides genome-editing targets for molecular breeding aiming at improving favorable agronomic properties.

Haploid Induction in Indica Rice: Exploring New Opportunities

Haploid plants are of significant interest to crop breeders due to their ability to expedite the development of inbred lines. Chromosome-doubling of haploids, produced by either in vitro or in vivo methods, results in fully homozygous doubled haploids. For nearly five decades, in vitro methods of anther and microspore culture have been attempted in many crops. In rice, in vitro methods are used with some success in japonica cultivars, although indica types have remained recalcitrant to a large extent. This review aims to explore the reasons for the lack of success of in vitro methods in indica rice and discuss new advancements in in vivo haploid induction protocols in other cereals and their relevance to rice. In particular, the current level of understanding of in vivo haploid inducer systems that utilize MTL and CENH3 mutants is analyzed in detail. One notable advantage of in vivo haploid induction systems is that they do not require tissue culture competence. This makes these methods more accessible and potentially transformative for research, offering a pragmatic approach to improving indica rice cultivars. By embracing these in vivo methods and harnessing the power of gene editing technologies like CRISPR/Cas9 systems, breeders can reshape their approach to indica rice improvement.

Cis- and Trans-variations of Stearoyl-CoA Desaturase Provide New Insights into the Mechanisms of Diverged Pattern of Phenotypic Plasticity for Temperature Adaptation in Two Congeneric Oyster Species

The evolution of phenotypic plasticity plays an essential role in adaptive responses to climate change; however, its regulatory mechanisms in marine organisms which exhibit high phenotypic plasticity still remain poorly understood. The temperature-responsive trait oleic acid content and its major gene stearoyl-CoA desaturase (Scd) expression have diverged in two allopatric congeneric oyster species, cold-adapted Crassostrea gigas and warm-adapted Crassostrea angulata. In this study, genetic and molecular methods were used to characterize fatty acid desaturation and membrane fluidity regulated by oyster Scd. Sixteen causative single-nucleotide polymorphisms (SNPs) were identified in the promoter/cis-region of the Scd between wild C. gigas and C. angulata. Further functional experiments showed that an SNP (g.-333C [C. gigas allele] >T [C. angulata allele]) may influence Scd transcription by creating/disrupting the binding motif of the positive trans-factor Y-box factor in C. gigas/C. angulata, which mediates the higher/lower constitutive expression of Scd in C. gigas/C. angulata. Additionally, the positive trans-factor sterol-regulatory element-binding proteins (Srebp) were identified to specifically bind to the promoter of Scd in both species, and were downregulated during cold stress in C. gigas compared to upregulated in C. angulata. This partly explains the relatively lower environmental sensitivity (plasticity) of Scd in C. gigas. This study serves as an experimental case to reveal that both cis- and trans-variations shape the diverged pattern of phenotypic plasticity, which provides new insights into the formation of adaptive traits and the prediction of the adaptive potential of marine organisms to future climate change.

Genome-wide association mapping of sodium and potassium concentration in rice grains and shoots under alternate wetting and drying and continuously flooded irrigation

Identification of a large number of QTL and candidate genes for sodium accumulation in a field grown population of rice derived from the aus subpopulation. Rice (Oryza sativa L.) is a globally important cereal crop. Sodium (Na⁺) and potassium (K⁺) are the major monovalent ions which affect rice growth, and exploring their uptake mechanisms will be useful for understanding rice biology. Since the balance of Na⁺ and K⁺ plays a significant role in adaptation of rice to salinity, that biology might inform the search for tolerance. In this study, the Na⁺ and K⁺ concentration and Na⁺/K⁺ ratio in grains and shoots were analyzed in the Bengal and Assam Aus Panel grown in field conditions under continuously flooded (CF) and alternate wetting and drying (AWD) irrigation. Overall, AWD irrigation significantly reduced the Na⁺ concentration and increased the K⁺ concentration in shoots and grains compared to the plants grown under CF. Genome-wide association mapping was conducted on Na⁺, K⁺ concentration and Na⁺/K⁺ ratio with 2 million SNPs using an efficient mixed model. Only QTLs which contained more than two significant SNPs (p < 0.0001) and where at least one of these significant SNPs passed a 10% false discovery rate were reported. A total of 106 QTLs were identified as being associated with Na⁺ concentration and Na⁺/K⁺ ratio across all traits and field conditions, with 48 QTLs found in multiple traits and/or water conditions. Four notable QTLs (one each on chromosomes 1 and 11, two on chromosome 2) and the haplotype variants of four candidate genes (OsHKT1;5, OsNHX2, LOC_Os02g32490 and OsFAD2_1) are discussed. The QTLs/candidate genes identified here could be useful for breeding rice that accumulates lower concentrations of sodium.

Cell Type-Specific Differentiation Between Indica and Japonica Rice Root Tip Responses to Different Environments Based on Single-Cell RNA Sequencing

Background: As Oryza sativa ssp. indica and Oryza sativa ssp. japonica are the two major subspecies of Asian cultivated rice, the adaptative evolution of these varieties in divergent environments is an important topic in both theoretical and practical studies. However, the cell type-specific differentiation between indica and japonica rice varieties in response to divergent habitat environments, which facilitates an understanding of the genetic basis underlying differentiation and environmental adaptation between rice subspecies at the cellular level, is little known.

Methods: We analyzed a published single-cell RNA sequencing dataset to explore the differentially expressed genes between indica and japonica rice varieties in each cell type. To estimate the relationship between cell type-specific differentiation and environmental adaptation, we focused on genes in the WRKY, NAC, and BZIP transcription factor families, which are closely related to abiotic stress responses. In addition, we integrated five bulk RNA sequencing datasets obtained under conditions of abiotic stress, including cold, drought and salinity, in this study. Furthermore, we analyzed quiescent center cells in rice root tips based on orthologous markers in Arabidopsis.

Results: We found differentially expressed genes between indica and japonica rice varieties with cell type-specific patterns, which were enriched in the pathways related to root development and stress reposes. Some of these genes were members of the WRKY, NAC, and BZIP transcription factor families and were differentially expressed under cold, drought or salinity stress. In addition, LOC_Os01g16810, LOC_Os01g18670, LOC_Os04g52960, and LOC_Os08g09350 may be potential markers of quiescent center cells in rice root tips.

Conclusion: These results identified cell type-specific differentially expressed genes between indica-japonica rice varieties that were related to various environmental stresses and provided putative markers of quiescent center cells. This study provides new clues for understanding the development and physiology of plants during the process of adaptative divergence, in addition to identifying potential target genes for the improvement of stress tolerance in rice breeding applications.

Genome-wide discovery of natural variation in pre-mRNA splicing and prioritising causal alternative splicing to salt stress response in rice

Pre-mRNA splicing is an essential step for the regulation of gene expression. In order to specifically capture splicing variants in plants for genome-wide association studies (GWAS), we developed a software tool to quantify and visualise Variations of Splicing in Population (VaSP). VaSP can quantify splicing variants from short-read RNA-seq datasets and discover genotype-specific splicing (GSS) events, which can be used to prioritise causal pre-mRNA splicing events in GWAS. We applied our method to an RNA-seq dataset with 328 samples from 82 genotypes from a rice diversity panel exposed to optimal and saline growing conditions. In total, 764 significant GSS events were identified in salt stress conditions. GSS events were used as markers for a GWAS with the shoot Na⁺ accumulation, which identified six GSS events in five genes significantly associated with the shoot Na⁺ content. Two of these genes, OsNUC1 and OsRAD23 emerged as top candidate genes with splice variants that exhibited significant divergence between the variants for shoot growth under salt stress conditions. VaSP is a versatile tool for alternative splicing analysis in plants and a powerful tool for prioritising candidate causal pre-mRNA splicing and corresponding genomic variations in GWAS.

Harnessing Knowledge from Maize and Rice Domestication for New Crop Breeding

Crop domestication has fundamentally altered the course of human history, causing a shift from hunter-gatherer to agricultural societies and stimulating the rise of modern civilization. A greater understanding of crop domestication would provide a theoretical basis for how we could improve current crops and develop new crops to deal with environmental challenges in a sustainable manner. Here, we provide a comprehensive summary of the similarities and differences in the domestication processes of maize and rice, two major staple food crops that feed the world. We propose that maize and rice might have evolved distinct genetic solutions toward domestication. Maize and rice domestication appears to be associated with distinct regulatory and evolutionary mechanisms. Rice domestication tended to select de novo, loss-of-function, coding variation, while maize domestication more frequently favored standing, gain-of-function, regulatory variation. At the gene network level, distinct genetic paths were used to acquire convergent phenotypes in maize and rice domestication, during which different central genes were utilized, orthologous genes played different evolutionary roles, and unique genes or regulatory modules were acquired for establishing new traits. Finally, we discuss how the knowledge gained from past domestication processes, together with emerging technologies, could be exploited to improve modern crop breeding and domesticate new crops to meet increasing human demands.

Genome-wide association mapping for grain manganese in rice (Oryza sativa L.) using a multi-experiment approach

Manganese (Mn) is an essential trace element for plants and commonly contributes to human health; however, the understanding of the genes controlling natural variation in Mn in crop plants is limited. Here, the integration of two of genome-wide association study approaches was used to increase the identification of valuable quantitative trait loci (QTL) and candidate genes responsible for the concentration of grain Mn across 389 diverse rice cultivars grown in Arkansas and Texas, USA, in multiple years. Single-trait analysis was initially performed using three different SNP datasets. As a result, significant loci could be detected using the high-density SNP dataset. Based on the 5.2 M SNP dataset, major QTLs were located on chromosomes 3 and 7 for Mn containing six candidate genes. In addition, the phenotypic data of grain Mn concentration were combined from three flooded-field experiments from the two sites and 3 years using multi-experiment analysis based on the 5.2 M SNP dataset. Two previous QTLs on chromosome 3 were identified across experiments, whereas new Mn QTLs were identified that were not found in individual experiments, on chromosomes 3, 4, 9 and 11. OsMTP8.1 was identified in both approaches and is a good candidate gene that could be controlling grain Mn concentration. This work demonstrates the utilisation of multi-experiment analysis to identify constitutive QTLs and candidate genes associated with the grain Mn concentration. Hence, the approach should be advantageous to facilitate genomic breeding programmes in rice and other crops considering QTLs and genes associated with complex traits in natural populations.