Rice LGD1 containing RNA binding activity affects growth and development through alternative promoters

A meta-analysis of rice phosphoproteomics data to understand variation in cell signalling across the rice pan-genome

Phosphorylation is the most studied post-translational modification, and has multiple biological functions. In this study, we have re-analysed publicly available mass spectrometry proteomics datasets enriched for phosphopeptides from Asian rice (Oryza sativa). In total we identified 15,522 phosphosites on serine, threonine and tyrosine residues on rice proteins. We identified sequence motifs for phosphosites, and link motifs to enrichment of different biological processes, indicating different downstream regulation likely caused by different kinase groups. We cross-referenced phosphosites against the rice 3,000 genomes, to identify single amino acid variations (SAAVs) within or proximal to phosphosites that could cause loss of a site in a given rice variety. The data was clustered to identify groups of sites with similar patterns across rice family groups, for example those highly conserved in Japonica, but mostly absent in Aus type rice varieties - known to have different responses to drought. These resources can assist rice researchers to discover alleles with significantly different functional effects across rice varieties. The data has been loaded into UniProt Knowledge-Base - enabling researchers to visualise sites alongside other data on rice proteins e.g. structural models from AlphaFold2, PeptideAtlas and the PRIDE database - enabling visualisation of source evidence, including scores and supporting mass spectra.

Combined genome-wide association study and epistasis analysis reveal multifaceted genetic architectures of plant height in Asian cultivated rice

Plant height (PH) in rice (Oryza sativa) is an important trait for its adaptation and agricultural performance. Discovery of the semi-dwarf1 (SD1) mutation initiated the Green Revolution, boosting rice yield and fitness, but the underlying genetic regulation of PH in rice remains largely unknown. Here, we performed genome-wide association study (GWAS) and identified 12 non-repetitive QTL/genes regulating PH variation in 619 Asian cultivated rice accessions. One of these was an SD1 structural variant, not normally detected in standard GWAS analyses. Given the strong effect of SD1 on PH, we also divided 619 accessions into subgroups harbouring distinct SD1 haplotypes, and found a further 85 QTL/genes for PH, revealing genetic heterogeneity that may be missed by analysing a broad, diverse population. Moreover, we uncovered two epistatic interaction networks of PH-associated QTL/genes in the japonica (Geng)-dominant SD1^NIP subgroup. In one of them, the hub QTL/gene qphSN1.4/GAMYB interacted with qphSN3.1/OsINO80, qphSN3.4/HD16/EL1, qphSN6.2/LOC_Os06g11130, and qphSN10.2/MADS56. Sequence variations in GAMYB and MADS56 were associated with their expression levels and PH variations, and MADS56 was shown to physically interact with MADS57 to coregulate expression of gibberellin (GA) metabolic genes OsGA2ox3 and Elongated Uppermost Internode1 (EUI1). Our study uncovered the multifaceted genetic architectures of rice PH, and provided novel and abundant genetic resources for breeding semi-dwarf rice and new candidates for further mechanistic studies on regulation of PH in rice.

A wheat integrative regulatory network from large-scale complementary functional datasets enables trait-associated gene discovery for crop improvement

Gene regulation is central to all aspects of organism growth, and understanding it using large-scale functional datasets can provide a whole view of biological processes controlling complex phenotypic traits in crops. However, the connection between massive functional datasets and trait-associated gene discovery for crop improvement is still lacking. In this study, we constructed a wheat integrative gene regulatory network (wGRN) by combining an updated genome annotation and diverse complementary functional datasets, including gene expression, sequence motif, transcription factor (TF) binding, chromatin accessibility, and evolutionarily conserved regulation. wGRN contains 7.2 million genome-wide interactions covering 5947 TFs and 127 439 target genes, which were further verified using known regulatory relationships, condition-specific expression, gene functional information, and experiments. We used wGRN to assign genome-wide genes to 3891 specific biological pathways and accurately prioritize candidate genes associated with complex phenotypic traits in genome-wide association studies. In addition, wGRN was used to enhance the interpretation of a spike temporal transcriptome dataset to construct high-resolution networks. We further unveiled novel regulators that enhance the power of spike phenotypic trait prediction using machine learning and contribute to the spike phenotypic differences among modern wheat accessions. Finally, we developed an interactive webserver, wGRN (http://wheat.cau.edu.cn/wGRN), for the community to explore gene regulation and discover trait-associated genes. Collectively, this community resource establishes the foundation for using large-scale functional datasets to guide trait-associated gene discovery for crop improvement.

AirMeasurer: open-source software to quantify static and dynamic traits derived from multiseason aerial phenotyping to empower genetic mapping studies in rice

Low-altitude aerial imaging, an approach that can collect large-scale plant imagery, has grown in popularity recently. Amongst many phenotyping approaches, unmanned aerial vehicles (UAVs) possess unique advantages as a consequence of their mobility, flexibility and affordability. Nevertheless, how to extract biologically relevant information effectively has remained challenging. Here, we present AirMeasurer, an open-source and expandable platform that combines automated image analysis, machine learning and original algorithms to perform trait analysis using 2D/3D aerial imagery acquired by low-cost UAVs in rice (Oryza sativa) trials. We applied the platform to study hundreds of rice landraces and recombinant inbred lines at two sites, from 2019 to 2021. A range of static and dynamic traits were quantified, including crop height, canopy coverage, vegetative indices and their growth rates. After verifying the reliability of AirMeasurer-derived traits, we identified genetic variants associated with selected growth-related traits using genome-wide association study and quantitative trait loci mapping. We found that the AirMeasurer-derived traits had led to reliable loci, some matched with published work, and others helped us to explore new candidate genes. Hence, we believe that our work demonstrates valuable advances in aerial phenotyping and automated 2D/3D trait analysis, providing high-quality phenotypic information to empower genetic mapping for crop improvement.

OsCKX2 regulates phosphate deficiency tolerance by modulating cytokinin in rice

Cytokinin oxidase/dehydrogenases (CKXs) are key enzymes that degrade cytokinins (CTKs) and play an essential role in plant growth and development. The present study analyzed the phenotypic and physiological characteristics of OsCKX2 overexpressing (OE) and knockout (KO) rice plants after exposure to phosphate (Pi) deficiency and the transcriptome and metabolome to investigate the function of OsCKX2 in response to Pi deficiency. OsCKX2 KO plants demonstrated higher endogenous CTK levels than wild-type (WT) under Pi deficiency. Further analysis indicated more robust tolerance of OsCKX2 KO plants to Pi deficiency, which exhibited higher phosphorus concentration, larger shoot biomass, and lesser leaf yellowing under Pi deficiency; whereas the opposite was observed for OsCKX2 OE plants. Transcriptome and metabolome analyses revealed that overexpression of OsCKX2 downregulated the transcriptional levels of genes related to Pi transporters, membrane lipid metabolism, and glycolysis, and reduced the consumption of metabolites in membrane lipid metabolism and glycolysis. On the contrary, knockout of OsCKX2 upregulated the expression of Pi transporters, and increased the consumption of metabolites in membrane lipid metabolism and glycolysis. These results indicated that OsCKX2 impacted Pi uptake, recycling, and plant growth via Pi transporters, phospholipid hydrolysis, and glycolysis under Pi deficiency. Overall, OsCKX2 negatively regulated Pi deficiency tolerance by modulating CTKs in rice.

Uncovering the genetic mechanisms regulating panicle architecture in rice with GPWAS and GWAS

Background

The number of panicles per plant, number of grains per panicle, and 1000-grain weight are important factors contributing to the grain yield per plant in rice. The Rice Diversity Panel 1 (RDP1) contains a total of 421 purified, homozygous rice accessions representing diverse genetic variations within O. sativa. The release of High-Density Rice Array (HDRA, 700 k SNPs) dataset provides a new opportunity to discover the genetic variants of panicle architectures in rice.

Results

In this report, a new method genome-phenome wide association study (GPWAS) was performed with 391 individuals and 27 traits derived from RDP1 to scan the relationship between the genes and multi-traits. A total of 1985 gene models were linked to phenomic variation with a p-value cutoff of 4.49E-18. Besides, 406 accessions derived from RDP1 with 411,066 SNPs were used to identify QTLs associated with the total spikelets number per panicle (TSNP), grain number per panicle (GNP), empty grain number per panicle (EGNP), primary branch number (PBN), panicle length (PL), and panicle number per plant (PN) by GLM, MLM, FarmCPU, and BLINK models for genome-wide association study (GWAS) analyses. A total of 18, 21, 18, 17, 15, and 17 QTLs were identified tightly linked with TSNP, GNP, EGNP, PBN, PL, and PN, respectively. Then, a total of 23 candidate genes were mapped simultaneously using both GWAS and GPWAS methods, composed of 6, 4, 5, 4, and 4 for TSNP, GNP, EGNP, PBN, and PL. Notably, one overlapped gene (Os01g0140100) were further investigated based on the haplotype and gene expression profile, indicating this gene might regulate the TSNP or panicle architecture in rice.

Conclusions

Nearly 30 % (30/106) QTLs co-located with the previous published genes or QTLs, indicating the power of GWAS. Besides, GPWAS is a new method to discover the relationship between genes and traits, especially the pleiotropy genes. Through comparing the results from GWAS and GPWAS, we identified 23 candidate genes related to panicle architectures in rice. This comprehensive study provides new insights into the genetic basis controlling panicle architectures in rice, which lays a foundation in rice improvement.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07391-x.

Genome-wide association study of agronomic traits in rice cultivated in temperate regions

BACKGROUND

Rice plants are sensitive to the agro-climate conditions, being photoperiod one of main factor contributing to their adaptation to the region where they are grown. Dissecting the genetic bases underlying diversity in rice populations adapted to specific environmental conditions is a fundamental resource for breeding. In this study we have analysed a collection of japonica varieties adapted to temperate regions to perform association studies with traits of high agronomical interest such as heading date, plant height, number of panicles, panicle length and number of grains per panicle.

RESULTS

We have performed a genome wide association study using a panel of 1713 SNPs that, based on previous linkage disequilibrium estimations, provides a full coverage of the whole genome. We have found a total of 43 SNPs associated with variations in the different traits. The identified SNPs were distributed across the genome except in chromosome 12, where no associated SNPs were found. The inspection of the vicinity of these markers also revealed a set of genes associated with physiological functions strongly linked to agronomic traits. Of special relevance are two genes involved in gibberellin homeostasis that are associated with plant height and panicle length. We also detected novel associated sites with heading date, panicle length and number of grain per panicle.

CONCLUSION

We have identified loci associated with important agronomic traits among cultivars adapted to temperate conditions. Some of these markers co-localized with already known genes or QTLs, but the association also provided novel molecular markers that can be of help to elucidate the complicated genetic mechanism controlling important agronomic traits, as flowering regulation in the non-dependent photoperiod pathway. The detected associated markers may provide important tools for the genetic improvement of rice cultivars in temperate regions.

Meta-expression analysis of unannotated genes in rice and approaches for network construction to suggest the probable roles

This work suggests 2020 potential candidates in rice for the functional annotation of unannotated genes using meta-analysis of anatomical samples derived from microarray and RNA-seq technologies and this information will be useful to identify novel morphological agronomic traits. Although the genome of rice (Oryza sativa) has been sequenced, 14,365 genes are considered unannotated because they lack putative annotation information. According to the Rice Genome Annotation Project Database ( http://rice.plantbiology.msu.edu/ ), the proportion of functionally characterized unannotated genes (0.35%) is quite limited when compared with the approximately 3.9% of annotated genes with assigned putative functions. Researchers require additional information to help them investigate the molecular mechanisms associated with those unannotated genes. To determine which of them might regulate morphological or physiological traits in the rice genome, we conducted a meta-analysis of expression data that covered a wide range of tissue/organ samples. Overall, 2020 genes showed cultivar-, tissue-, or organ-preferential patterns of expression. Representative candidates from featured groups were validated by RT-PCR, and the GUS reporter system was used to validate the expression of genes that were clustered according to their leaf or root preference. Taking a molecular and genetics approach, we examined meta-expression data and found that 127 genes were differentially expressed between japonica and indica rice cultivars. This is potentially significant for future agronomic applications. We also used a T-DNA insertional mutant and performed a co-expression network analysis of Sword shape dwarf1 (SSD1), a gene that regulates cell division. This network was refined via RT-PCR analysis. Our results suggested that SSD1 represses the expression of four genes related to the processes of DNA replication or cell division and provides insight into possible molecular mechanisms. Together, these strategies present a valuable tool for in-depth characterization of currently unannotated genes.

Genome-wide identification and analysis of rice genes preferentially expressed in pollen at an early developmental stage

Microspore production using endogenous developmental programs has not been well studied. The main limitation is the difficulty in identifying genes preferentially expressed in pollen grains at early stages. To overcome this limitation, we collected transcriptome data from anthers and microspore/pollen and performed meta-expression analysis. Subsequently, we identified 410 genes showing preferential expression patterns in early developing pollen samples of both japonica and indica cultivars. The expression patterns of these genes are distinguishable from genes showing pollen mother cell or tapetum-preferred expression patterns. Gene Ontology enrichment and MapMan analyses indicated that microspores in rice are closely linked with protein degradation, nucleotide metabolism, and DNA biosynthesis and regulation, while the pollen mother cell or tapetum are strongly associated with cell wall metabolism, lipid metabolism, secondary metabolism, and RNA biosynthesis and regulation. We also generated transgenic lines under the control of the promoters of eight microspore-preferred genes and confirmed the preferred expression patterns in plants using the GUS reporting system. Furthermore, cis-regulatory element analysis revealed that pollen specific elements such as POLLEN1LELAT52, and 5659BOXLELAT5659 were commonly identified in the promoter regions of eight rice genes with more frequency than estimation. Our study will provide new sights on early pollen development in rice, a model crop plant.

Genetic resources offer efficient tools for rice functional genomics research

Rice is an important crop and major model plant for monocot functional genomics studies. With the establishment of various genetic resources for rice genomics, the next challenge is to systematically assign functions to predicted genes in the rice genome. Compared with the robustness of genome sequencing and bioinformatics techniques, progress in understanding the function of rice genes has lagged, hampering the utilization of rice genes for cereal crop improvement. The use of transfer DNA (T-DNA) insertional mutagenesis offers the advantage of uniform distribution throughout the rice genome, but preferentially in gene-rich regions, resulting in direct gene knockout or activation of genes within 20-30 kb up- and downstream of the T-DNA insertion site and high gene tagging efficiency. Here, we summarize the recent progress in functional genomics using the T-DNA-tagged rice mutant population. We also discuss important features of T-DNA activation- and knockout-tagging and promoter-trapping of the rice genome in relation to mutant and candidate gene characterizations and how to more efficiently utilize rice mutant populations and datasets for high-throughput functional genomics and phenomics studies by forward and reverse genetics approaches. These studies may facilitate the translation of rice functional genomics research to improvements of rice and other cereal crops.

Differential gene expression and alternative splicing between diploid and tetraploid watermelon

The exploitation of synthetic polyploids for producing seedless fruits is well known in watermelon. Tetraploid progenitors of triploid watermelon plants, compared with their diploid counterparts, exhibit wide phenotypic differences. Although many factors modulate alternative splicing (AS) in plants, the effects of autopolyploidization on AS are still unknown. In this study, we used tissues of leaf, stem, and fruit of diploid and tetraploid sweet watermelon to understand changes in gene expression and the occurrence of AS. RNA-sequencing analysis was performed along with reverse transcription quantitative PCR and rapid amplification of cDNA ends (RACE)-PCR to demonstrate changes in expression and splicing. All vegetative tissues except fruit showed an increased level of AS in the tetraploid watermelon throughout the growth period. The ploidy levels of diploids and the tetraploid were confirmed using a ploidy analyser. We identified 5362 and 1288 genes that were up- and downregulated, respectively, in tetraploid as compared with diploid plants. We further confirmed that 22 genes underwent AS events across tissues, indicating possibilities of generating different protein isoforms with altered functions of important transcription factors and transporters. Arginine biosynthesis, chlorophyllide synthesis, GDP mannose biosynthesis, trehalose biosynthesis, and starch and sucrose degradation pathways were upregulated in autotetraploids. Phloem protein 2, chloroplastic PGR5-like protein, zinc-finger protein, fructokinase-like 2, MYB transcription factor, and nodulin MtN21 showed AS in fruit tissues. These results should help in developing high-quality seedless watermelon and provide additional transcriptomic information related to other cucurbits.

Arabidopsis mTERF15 is required for mitochondrial nad2 intron 3 splicing and functional complex I activity

Mitochondria play a pivotal role in most eukaryotic cells, as they are responsible for the generation of energy and diverse metabolic intermediates for many cellular events. During endosymbiosis, approximately 99% of the genes encoded by the mitochondrial genome were transferred into the host nucleus, and mitochondria import more than 1000 nuclear-encoded proteins from the cytosol to maintain structural integrity and fundamental functions, including DNA replication, mRNA transcription and RNA metabolism of dozens of mitochondrial genes. In metazoans, a family of nuclear-encoded proteins called the mitochondrial transcription termination factors (mTERFs) regulates mitochondrial transcription, including transcriptional termination and initiation, via their DNA-binding activities, and the dysfunction of individual mTERF members causes severe developmental defects. Arabidopsis thaliana and Oryza sativa contain 35 and 48 mTERFs, respectively, but the biological functions of only a few of these proteins have been explored. Here, we investigated the biological role and molecular mechanism of Arabidopsis mTERF15 in plant organelle metabolism using molecular genetics, cytological and biochemical approaches. The null homozygous T-DNA mutant of mTERF15, mterf15, was found to result in substantial retardation of both vegetative and reproductive development, which was fully complemented by the wild-type genomic sequence. Surprisingly, mitochondria-localized mTERF15 lacks obvious DNA-binding activity but processes mitochondrial nad2 intron 3 splicing through its RNA-binding ability. Impairment of this splicing event not only disrupted mitochondrial structure but also abolished the activity of mitochondrial respiratory chain complex I. These effects are in agreement with the severe phenotype of the mterf15 homozygous mutant. Our study suggests that Arabidopsis mTERF15 functions as a splicing factor for nad2 intron 3 splicing in mitochondria, which is essential for normal plant growth and development.

Genetic analysis of vegetative branching in sorghum

We identified quantitative trait loci influencing plant architecture that may be valuable in breeding of optimized genotypes for sustainable food and/or cellulosic biomass production, and advancing resilience to changing climates. We describe a 3-year study to identify quantitative trait loci (QTLs) for vegetative branching of sorghum in a recombinant inbred line population of 161 genotypes derived from two morphologically distinct parents, S. bicolor × S. propinquum. We quantify vegetative branching based on morphological position and physiological status. Different sets of QTLs for different levels of branching were identified. QTLs discovered on chromosomes 1, 3, 7 and 8 affect multiple vegetative branching variables, suggesting that these regions may contain genes that control general axillary meristem initiation. Other regions that only influence one vegetative branching trait could contain genes that influence developmental processes contributing to divergent patterns of plant architecture. We investigate the relationship between vegetative branching patterns and dry biomass, and conclude that tillers with mature panicles and immature secondary branches each show consistent positive correlation with dry biomass. Among 19 branching-related genes from rice, eight sorghum homologs of seven rice genes are in syntenic blocks within branching-related QTL likelihood intervals. Five of these eight genes are within 700 kb of SNPs significantly associated with differences in branching in genome-wide association study of a diversity panel of 377 sorghum accessions, and three contain striking allelic variations between S. bicolor and S. propinquum that are likely to impact gene functions. Unraveling genetic determinants for vegetative branching may contribute to deterministic breeding of optimized genotypes for sustainable food and cellulosic biomass production in both optimal and marginal conditions, which are resilient to future climates that are more volatile and more stressful.

Identification and validation of promoters and cis-acting regulatory elements

Studies of promoters that largely regulate gene expression at the transcriptional level are crucial for improving our basic understanding of gene regulation and will expand the toolbox of available promoters for use in plant biotechnology. In this review, we present a comprehensive analysis of promoters and their underlying mechanisms in transcriptional regulation, including epigenetic marks and chromatin-based regulation. Large-scale prediction of promoter sequences and their contributing cis-acting elements has become routine due to recent advances in transcriptomic technologies and genome sequencing of several plants. However, predicted regulatory sequences may or may not be functional and demonstration of the contribution of the element to promoter activity is essential for confirmation of regulatory sequences. Synthetic promoters and introns provide useful approaches for functional validation of promoter sequences. The development and improvement of gene expression tools for rapid, efficient, predictable, and high-throughput analysis of promoter components will be critical for confirmation of the functional regulatory element sequences identified through transcriptomic and genomic analyses.

Proteins involved in distinct phases of cold hardening process in frost resistant winter barley (Hordeum vulgare L.) cv Luxor

Winter barley is an economically important cereal crop grown in higher latitudes and altitudes where low temperatures represent an important environmental constraint limiting crop productivity. In this study changes in proteome of leaves and crowns in a frost tolerant winter barley cv. Luxor in relation to short and long term periods of cold followed by a brief frost treatment were studied in order to disclose proteins responsible for the cold hardening process in distinct plant tissues. The mentioned changes have been monitored using two dimensional difference gel electrophoresis (2D-DIGE) with subsequent peptide-mapping protein identification. Regarding approximately 600–700 distinct protein spots detected on 2D gels, there has been found at least a two-fold change after exposure to low temperatures in about 10% of proteins in leaves and 13% of proteins in crowns. Protein and nitrogen metabolic processes have been influenced by low temperature to a similar extent in both tissues while catabolism, carbohydrate metabolism and proteins involved in stress response have been more affected in crowns than in leaves. The range of changes in protein abundance was generally higher in leaves and chloroplast proteins were frequently affected which suggests a priority to protect photosynthetic apparatus. Overall, our data proved existence of slightly different response strategies to low temperature stress in crowns and leaves, i.e., tissues with different biological role. Moreover, there have been found several proteins with large increase in accumulation, e.g., 33 kDa oxygen evolving protein of photosystem II in leaves and “enhanced disease susceptibility 1” in crowns; these proteins might have potential to indicate an enhanced level of frost tolerance in barley.