Re-programming of gene expression in the CS8 rice line over-expressing ADPglucose pyrophosphorylase induces a suppressor of starch biosynthesis

High Daytime Temperature Responsive MicroRNA Profiles in Developing Grains of Rice Varieties with Contrasting Chalkiness

High temperature impairs starch biosynthesis in developing rice grains and thereby increases chalkiness, affecting the grain quality. Genome encoded microRNAs (miRNAs) fine-tune target transcript abundances in a spatio-temporal specific manner, and this mode of gene regulation is critical for a myriad of developmental processes as well as stress responses. However, the role of miRNAs in maintaining rice grain quality/chalkiness during high daytime temperature (HDT) stress is relatively unknown. To uncover the role of miRNAs in this process, we used five contrasting rice genotypes (low chalky lines Cyp, Ben, and KB and high chalky lines LaGrue and NB) and compared the miRNA profiles in the R6 stage caryopsis samples from plants subjected to prolonged HDT (from the onset of fertilization through R6 stage of caryopsis development). Our small RNA analysis has identified approximately 744 miRNAs that can be grouped into 291 families. Of these, 186 miRNAs belonging to 103 families are differentially regulated under HDT. Only two miRNAs, Osa-miR444f and Osa-miR1866-5p, were upregulated in all genotypes, implying that the regulations greatly varied between the genotypes. Furthermore, not even a single miRNA was commonly up/down regulated specifically in the three tolerant genotypes. However, three miRNAs (Osa-miR1866-3p, Osa-miR5150-3p and canH-miR9774a,b-3p) were commonly upregulated and onemiRNA (Osa-miR393b-5p) was commonly downregulated specifically in the sensitive genotypes (LaGrue and NB). These observations suggest that few similarities exist within the low chalky or high chalky genotypes, possibly due to high genetic variation. Among the five genotypes used, Cypress and LaGrue are genetically closely related, but exhibit contrasting chalkiness under HDT, and thus, a comparison between them is most relevant. This comparison revealed a general tendency for Cypress to display miRNA regulations that could decrease chalkiness under HDT compared with LaGrue. This study suggests that miRNAs could play an important role in maintaining grain quality in HDT-stressed rice.

Cytological, transcriptome and miRNome temporal landscapes decode enhancement of rice grain size

BACKGROUND

Rice grain size (GS) is an essential agronomic trait. Though several genes and miRNA modules influencing GS are known and seed development transcriptomes analyzed, a comprehensive compendium connecting all possible players is lacking. This study utilizes two contrasting GS indica rice genotypes (small-grained SN and large-grained LGR). Rice seed development involves five stages (S1-S5). Comparative transcriptome and miRNome atlases, substantiated with morphological and cytological studies, from S1-S5 stages and flag leaf have been analyzed to identify GS proponents.

RESULTS

Histology shows prolonged endosperm development and cell enlargement in LGR. Stand-alone and comparative RNAseq analyses manifest S3 (5-10 days after pollination) stage as crucial for GS enhancement, coherently with cell cycle, endoreduplication, and programmed cell death participating genes. Seed storage protein and carbohydrate accumulation, cytologically and by RNAseq, is shown to be delayed in LGR. Fourteen transcription factor families influence GS. Pathway genes for four phytohormones display opposite patterns of higher expression. A total of 186 genes generated from the transcriptome analyses are located within GS trait-related QTLs deciphered by a cross between SN and LGR. Fourteen miRNA families express specifically in SN or LGR seeds. Eight miRNA-target modules display contrasting expressions amongst SN and LGR, while 26 (SN) and 43 (LGR) modules are differentially expressed in all stages.

CONCLUSIONS

Integration of all analyses concludes in a "Domino effect" model for GS regulation highlighting chronology and fruition of each event. This study delineates the essence of GS regulation, providing scope for future exploits. The rice grain development database (RGDD) ( www.nipgr.ac.in/RGDD/index.php ; https://doi.org/10.5281/zenodo.7762870 ) has been developed for easy access of data generated in this paper.

Co-expression pan-network reveals genes involved in complex traits within maize pan-genome

BACKGROUND

With the advances in the high throughput next generation sequencing technologies, genome-wide association studies (GWAS) have identified a large set of variants associated with complex phenotypic traits at a very fine scale. Despite the progress in GWAS, identification of genotype-phenotype relationship remains challenging in maize due to its nature with dozens of variants controlling the same trait. As the causal variations results in the change in expression, gene expression analyses carry a pivotal role in unraveling the transcriptional regulatory mechanisms behind the phenotypes.

RESULTS

To address these challenges, we incorporated the gene expression and GWAS-driven traits to extend the knowledge of genotype-phenotype relationships and transcriptional regulatory mechanisms behind the phenotypes. We constructed a large collection of gene co-expression networks and identified more than 2 million co-expressing gene pairs in the GWAS-driven pan-network which contains all the gene-pairs in individual genomes of the nested association mapping (NAM) population. We defined four sub-categories for the pan-network: (1) core-network contains the highest represented ~ 1% of the gene-pairs, (2) near-core network contains the next highest represented 1-5% of the gene-pairs, (3) private-network contains ~ 50% of the gene pairs that are unique to individual genomes, and (4) the dispensable-network contains the remaining 50-95% of the gene-pairs in the maize pan-genome. Strikingly, the private-network contained almost all the genes in the pan-network but lacked half of the interactions. We performed gene ontology (GO) enrichment analysis for the pan-, core-, and private- networks and compared the contributions of variants overlapping with genes and promoters to the GWAS-driven pan-network.

CONCLUSIONS

Gene co-expression networks revealed meaningful information about groups of co-regulated genes that play a central role in regulatory processes. Pan-network approach enabled us to visualize the global view of the gene regulatory network for the studied system that could not be well inferred by the core-network alone.

Starch synthases SSIIa and GBSSI control starch structure but do not determine starch granule morphology in the absence of SSIIIa and SSIVb

High levels of two major starch synthases, SSIIa and GBSSI, in ss3a ss4b double mutant rice alter the starch structure but fail to recover the polygonal starch granule morphology. The endosperm starch granule is polygonal in wild-type rice but spherical in double mutant japonica rice lacking genes encoding two of the five major Starch synthase (SS) isozymes expressed in endosperm, SSIIIa and SSIVb. Japonica rice naturally has low levels of SSIIa and Granule-bound SSI (GBSSI). Therefore, introduction of active SSIIa allele and/or high-expressing GBSSI allele from indica rice into the japonica rice mutant lacking SS isozymes can help elucidate the compensatory roles of SS isozymes in starch biosynthesis. In this study, we crossed the ss3a ss4a double mutant japonica rice with the indica rice to generate three new rice lines with high and/or low SSIIa and GBSSI levels, and examined their starch structure, physicochemical properties, and levels of other starch biosynthetic enzymes. Lines with high SSIIa levels showed more SSI and SSIIa bound to starch granule, reduced levels of short amylopectin chains (7 ≤ DP ≤ 12), increased levels of amylopectin chains with DP > 13, and consequently higher gelatinization temperature. Lines with high GBSSI levels showed an increase in amylose content. The ADP-glucose content of the crude extract was high in lines with low or high SSIIa and low GBSSI levels, but was low in lines with high GBSSI. Addition of high SSIIa and GBSSI altered the starch structure and physicochemical properties but did not affect the starch granule morphology, confirming that SSIIIa and SSIVb are key enzymes affecting starch granule morphology in rice. The relationship among SS isozymes and its effect on the amount of substrate (ADP-glucose) is discussed.

Evolutionary expansion and functional divergence of sugar transporters in Saccharum (S. spontaneum and S. officinarum)

The sugar transporter (ST) family is considered to be the most important gene family for sugar accumulation, but limited information about the ST family in the important sugar-yielding crop Saccharum is available due to its complex genetic background. Here, 105 ST genes were identified and clustered into eight subfamilies in Saccharum spontaneum. Comparative genomics revealed that tandem duplication events contributed to ST gene expansions of two subfamilies, PLT and STP, in S. spontaneum, indicating an early evolutionary step towards high sugar content in Saccharum. The analyses of expression patterns were based on four large datasets with a total of 226 RNA sequencing samples from S. spontaneum and Saccharum officinarum. The results clearly demonstrated 50 ST genes had different spatiotemporal expression patterns in leaf tissues, 10 STs were specifically expressed in the stem, and 10 STs responded to the diurnal rhythm. Heterologous expression experiments in the defective yeast strain EBY.VW4000 indicated STP13, pGlcT2, VGT3, and TMT4 are the STs with most affinity for glucose/fructose and SUT1_T1 has the highest affinity to sucrose. Furthermore, metabolomics analysis suggested STP7 is a sugar starvation-induced gene and STP13 has a function in retrieving sugar in senescent tissues. PLT11, PLT11_T1, TMT3, and TMT4 contributed to breaking the limitations of the storage sink. SUT1, SUT1_T1, PLT11, TMT4, pGlcT2, and VGT3 responded for different functions in these two Saccharum species. This study demonstrated the evolutionary expansion and functional divergence of the ST gene family and will enable the further investigation of the molecular mechanism of sugar metabolism in Saccharum.

The NAC transcription factor NAC019-A1 is a negative regulator of starch synthesis in wheat developing endosperm

Starch is a major component of wheat (Triticum aestivum L.) endosperm and is an important part of the human diet. The functions of many starch synthesis genes have been elucidated. However, little is known about their regulatory mechanisms in wheat. Here, we identified a novel NAC transcription factor, TaNAC019-A1 (TraesCS3A02G077900), that negatively regulates starch synthesis in wheat and rice (Oryza sativa L.) endosperms. TaNAC019-A1 was highly expressed in the endosperm of developing grains and encoded a nucleus-localized transcriptional repressor. Overexpression of TaNAC019-A1 in rice and wheat led to significantly reduced starch content, kernel weight, and kernel width. The TaNAC019-A1-overexpression wheat lines had smaller A-type starch granules and fewer B-type starch granules than wild-type. Moreover, TaNAC019-A1 could directly bind to the 'ACGCAG' motif in the promoter regions of ADP-glucose pyrophosphorylase small subunit 1 (TaAGPS1-A1, TraesCS7A02G287400) and TaAGPS1-B1 (TraesCS7B02G183300) and repress their expression, thereby inhibiting starch synthesis in wheat endosperm. One haplotype of TaNAC019-B1 (TaNAC019-B1-Hap2, TraesCS3B02G092800) was positively associated with thousand-kernel weight and underwent positive selection during the Chinese wheat breeding process. Our data demonstrate that TaNAC019-A1 is a negative regulator of starch synthesis in wheat endosperm and provide novel insight into wheat yield improvement.

Identifying genes for resistant starch, slowly digestible starch, and rapidly digestible starch in rice using genome-wide association studies

BACKGROUND

The digestibility of starch is important for the nutritive value of staple food. Although several genes are responsible for resistant starch (RS) and slowly digestible starch (SDS), gaps persist concerning the molecular basis of RS and SDS formation due to the complex genetic mechanisms of starch digestibility.

OBJECTIVES

The objective of this study was to identify new genes for starch digestibility in rice and interprete the genetic mechanisms of RS and SDS by GWAS.

METHODS

Genome-wide association studies were conducted by associating the RS and SDS phenotypes of 104 re-sequenced rice lines to an SNP dataset of 2,288,867 sites using a compressed mixed linear model. Candidate genes were identified according to the position of the SNPs based on data from the MSU Rice Genome Annotation Project.

RESULTS

Seven quantitative trait loci (QTLs) were detected to be associated with the RS content, among which the SNP 6 m1765761 was located on Waxy. Starch branching enzymes IIa (BEIIa) close to QTL qRS-I4 was detected and further identified as a specific candidate gene for RS in INDICA. Two QTLs were associated with SDS, and the LOC_Os09g09360 encoding lipase was identified as a causal gene for SDS.

CONCLUSIONS

GWAS is a valid strategy to genetically dissect the formation of starch digestion properties in rice. RS formation in grains is dependent on the rice type; lipid might also contribute to starch digestibility and should be an alternative factor to improve rice starch digestibility.

Deep-reaching acceleration of global mean ocean circulation over the past two decades

Ocean circulation redistributes Earth's energy and water masses and influences global climate. Under historical greenhouse warming, regional ocean currents show diverse tendencies, but whether there is an emerging trend of the global mean ocean circulation system is not yet clear. Here, we show a statistically significant increasing trend in the globally integrated oceanic kinetic energy since the early 1990s, indicating a substantial acceleration of global mean ocean circulation. The increasing trend in kinetic energy is particularly prominent in the global tropical oceans, reaching depths of thousands of meters. The deep-reaching acceleration of the ocean circulation is mainly induced by a planetary intensification of surface winds since the early 1990s. Although possibly influenced by wind changes associated with the onset of a negative Pacific decadal oscillation since the late 1990s, the recent acceleration is far larger than that associated with natural variability, suggesting that it is principally part of a long-term trend.