Analysing the rice young panicle transcriptome reveals the gene regulatory network controlled by TRIANGULAR HULL1

Characterization and Expression Analysis of the ALOG Gene Family in Rice (Oryza sativa L.)

ALOG (Arabidopsis LSH1 and Oryza G1) proteins constitute a plant-specific family of transcription factors that play crucial roles in lateral organ development across land plants. Initially identified through forward genetic studies of Arabidopsis LSH1 and rice G1 proteins, ALOG family members have since been functionally characterized in various plant species. However, research focusing on the characteristics and expression patterns of all ALOG family members in rice remains relatively limited. In this study, we systematically characterized OsALOG family genes in rice. Compared to other genes in rice and Arabidopsis, the ALOG family genes have a relatively simple structure. The alignment of OsALOG amino acid sequences and analysis of disorder predictions reveal that all members possess conserved ALOG domains, while the conservation of intrinsically disordered regions (IDRs) is relatively low. Four amino acids-alanine, glycine, proline, and serine-are significantly enriched in the IDRs of each ALOG protein. Synteny analysis indicates that most OsALOG genes have undergone considerable divergence compared to their counterparts in Arabidopsis. Bioinformatic analysis of cis-regulatory elements predicts that OsALOG family genes contain elements responsive to ABA, light, and methyl jasmonate, although the abundance and composition of these elements vary among different members. The expression patterns associated with the rice floral development of OsALOG genes can be broadly categorized into two types; however, even within the same type, differences in expression levels, as well as the initiation time and duration of expression, were observed. These results provide a comprehensive understanding of the structural characteristics and expression patterns of OsALOG members in rice.

Comprehensive identification and analysis of DUF640 genes associated with rice growth

Protein domains with conserved amino acid sequences and uncharacterized functions are called domains of unknown function (DUF). The DUF640 gene family plays a crucial role in plant growth, particularly in light regulation, floral organ development, and fruit development. However, there exists a lack of systematic understanding of the evolutionary relationships and functional differentiation of DUF640 within the Oryza genus. In this study, 61 DUF640 genes were identified in the Oryza genus. The expression of DUF640s is induced by multiple hormonal stressors including abscisic acid (ABA), cytokinin (CK), ethylene (ETH), and indole-3-acetic acid (IAA). Specifically, OiDUF640-10 expression significantly increased after ETH treatment. Transgenic experiments showed that overexpressing OiDUF640-10 lines were sensitive to ETH, and seedling length was obstructed. Evolutionary analysis revealed differentiation of the OiDUF640-10 gene in O. sativa ssp. indica and japonica varieties, likely driven by natural selection during the domestication of cultivated rice. These results indicate that OiDUF640-10 plays a vital role in the regulation of rice seedling length.

The ALOG family members OsG1L1 and OsG1L2 regulate inflorescence branching in rice

The architecture of the rice inflorescence is an important determinant of crop yield. The length of the inflorescence and the number of branches are among the key factors determining the number of spikelets, and thus grains, that a plant will develop. In particular, the timing of the identity transition from indeterminate branch meristem to determinate spikelet meristem governs the complexity of the inflorescence. In this context, the ALOG gene TAWAWA1 (TAW1) has been shown to delay the transition to determinate spikelet development in Oryza sativa (rice). Recently, by combining precise laser microdissection of inflorescence meristems with RNA-seq, we observed that two ALOG genes, OsG1-like 1 (OsG1L1) and OsG1L2, have expression profiles similar to that of TAW1. Here, we report that osg1l1 and osg1l2 loss-of-function CRISPR mutants have similar phenotypes to the phenotype of the previously published taw1 mutant, suggesting that these genes might act on related pathways during inflorescence development. Transcriptome analysis of the osg1l2 mutant suggested interactions of OsG1L2 with other known inflorescence architecture regulators and the data sets were used for the construction of a gene regulatory network (GRN), proposing interactions among genes potentially involved in controlling inflorescence development in rice. In this GRN, we selected the homeodomain-leucine zipper transcription factor encoding the gene OsHOX14 for further characterization. The spatiotemporal expression profiling and phenotypical analysis of CRISPR loss-of-function mutants of OsHOX14 suggests that the proposed GRN indeed serves as a valuable resource for the identification of new proteins involved in rice inflorescence development.

Transcriptome and Proteome Analysis Revealed That Hormone and Reactive Oxygen Species Synergetically Regulate Dormancy of Introgression Line in Rice (Oryza sativa L.)

Dormancy is a complex agronomy phenotype controlled by multiple signaling and a key trait repressing pre-harvest sprouting (PHS). However, the signaling network of dormancy remains unclear. In this study, we used Zhonghua11 (ZH11) with a weak dormancy, and Introgression line (IL) with a strong dormancy to study the mechanism of hormones and reactive oxygen species (ROS) crosstalk regulating rice dormancy. The germination experiment showed that the germination rate of ZH11 was 76.86%, while that of IL was only 1.25%. Transcriptome analysis showed that there were 1658 differentially expressed genes (DEGs) between IL and ZH11, of which 577 were up-regulated and 1081 were down-regulated. Additionally, DEGs were mainly enriched in oxidoreductase activity, cell periphery, and plant hormone signal transduction pathways. Tandem mass tags (TMT) quantitative proteomics analysis showed 275 differentially expressed proteins (DEPs) between IL and ZH11, of which 176 proteins were up-regulated, 99 were down-regulated, and the DEPs were mainly enriched in the metabolic process and oxidation-reduction process. The comprehensive transcriptome and proteome analysis showed that their correlation was very low, and only 56 genes were co-expressed. Hormone content detection showed that IL had significantly lower abscisic acid (ABA) contents than the ZH11 while having significantly higher jasmonic acid (JA) contents than the ZH11. ROS content measurement showed that the hydrogen peroxide (H₂O₂) content of IL was significantly lower than the ZH11, while the production rate of superoxide anion (O₂^.-) was significantly higher than the ZH11. These results indicate that hormones and ROS crosstalk to regulate rice dormancy. In particular, this study has deepened our mechanism of ROS and JA crosstalk regulating rice dormancy and is conducive to our precise inhibition of PHS.

Combined Analysis of BSA-Seq Based Mapping, RNA-Seq, and Metabolomic Unraveled Candidate Genes Associated with Panicle Grain Number in Rice (Oryza sativa L.)

Rice grain yield is a complex and highly variable quantitative trait consisting of several key components, including the grain weight, the effective panicles per unit area, and the grain number per panicle (GNPP). The GNPP is a significant contributor to grain yield controlled by multiple genes (QTL) and is crucial for improvement. Attempts have been made to find genes for this trait, which has always been a challenging and arduous task through conventional methods. We combined a BSA analysis, RNA profiling, and a metabolome analysis in the present study to identify new candidate genes involved in the GNPP. The F₂ population from crossing R4233 (high GNPP) and Ce679 (low GNPP) revealed a frequency distribution fitting two segregated genes. Three pools, including low, middle, and high GNPP, were constructed and a BSA analysis revealed six candidate regions spanning 5.38 Mb, containing 739 annotated genes. Further, a conjunctive analysis of BSA-Seq and RNA-Seq showed 31 differentially expressed genes (DEGs) in the candidate intervals. Subsequently, a metabolome analysis showed 1024 metabolites, with 71 significantly enriched, including 44 up and 27 downregulated in Ce679 vs. R4233. A KEGG enrichment analysis of these 31 DEGs and 71 differentially enriched metabolites (DEMs) showed two genes, Os12g0102100 and Os01g0580500, significantly enriched in the metabolic pathways’ biosynthesis of secondary metabolites, cysteine and methionine metabolism, and fatty acid biosynthesis. Os12g0102100, which encodes for the alcohol dehydrogenase superfamily and a zinc-containing protein, is a novel gene whose contribution to the GNPP is not yet elucidated. This gene coding for mitochondrial trans-2-enoyl-CoA reductase is involved in the biosynthesis of myristic acid, also known as tetradecanoic acid. The Os01g0580500 coding for the enzyme 1-aminoclopropane-1-carboxylate oxidase (OsACO7) is responsible for the final step of the ethylene biosynthesis pathway through the conversion of 1-aminocyclopropane-1-carboxylic acid (ACC) into ethylene. Unlike Os12g0102100, this gene was significantly upregulated in R4233, downregulated in Ce679, and significantly enriched in two of the three metabolite pathways. This result pointed out that these two genes are responsible for the difference in the GNPP in the two cultivars, which has never been identified. Further validation studies may disclose the physiological mechanisms through which they regulate the GNPP in rice.

Comprehensive Transcriptome Analysis of GS3 Near-Isogenic Lines During Panicle Development in Rice (Oryza sativa L.)

Panicle architecture is an important agronomic trait in rice that affects rice yields and quality. The GRAIN SIZE 3 (GS3) locus has been identified as a major quantitative trait locus (QTL) affecting grain length and weight. The current understanding of the function of the GS3 gene, especially concerning the regulatory mechanism of panicle development, is still in its infancy. In this study, we generated GS3 near-isogenic lines (NILs) by successive crossing and backcrossing of TD70 (large grain) with Kasalath (small grain), using Kasalath as the recurrent parent. To identify potential transcription dynamic changes in rice panicle formation and grain shape, we deeply analyzed transcriptional profiles for the NILs (NIL-GS3 and NIL-gs3) at three different panicle developmental stages (S, M, and L). A total of 887, 1,768, and 1,478 differentially expressed genes (DEGs) were identified at stages S, M, and L, respectively. We also found 542 differential expressed long non-coding RNAs (lncRNAs). Co-expression analysis further revealed significant clusters associated with different development periods in NIL-gs3 lines. Gene Ontology and KEGG enrichment analysis revealed G-protein signaling and hormones pathway were successively activated at the M and L stages of NIL-gs3, which indicated activation of the G-protein signaling pathway might trigger the down-streaming hormone signaling transduction. we found that other hormones such ABA, Auxin, CK were significantly enriched in the L stage in the NIL-gs3. We highlighted the synergistic interplay of G-protein and multiple hormones signaling pathways and their essential roles in regulating rice panicle formation and the grain shape. Our study provides an invaluable resource for further molecular mechanistic studies that affect rice grain size and provide new insight for directed selection by marker-assisted backcross breeding.

Application of Allele Specific PCR in Identifying Offspring Genotypes of Bi-Allelic SbeIIb Mutant Lines in Rice

Bi-allelic mutant lines induced by clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR associated (Cas) systems are important genetic materials. It is very important to establish a rapid and cheap method in identifying homozygous mutant plants from offspring segregation populations of bi-allelic mutant lines. In this study, the offspring genotypes of rice bi-allelic starch branching enzyme IIb mutant lines were identified using the allele specific PCR (AS-PCR) method. The target sequences of two alleles were aligned from their 5′ to 3′ ends, and the first different bases were used as the 3′ ends of mismatch primers. Another mismatched base was introduced at the third nucleotide from the 3′ end of mismatch primer. The PCR reaction mixture and amplification program were optimized according to the differences of mutation target sequence and mismatch primers. The offspring plant genotypes of bi-allelic mutant lines could be accurately identified using the amplified DNA fragments by agarose gel electrophoresis. This study could provide a method reference for the rapid screening of homozygous mutant plants from offspring segregation population of heterozygous and bi-allelic mutant lines.

Transcriptional analysis for the difference in carotenoids accumulation in flesh and peel of white-fleshed loquat fruit

Loquat (Eriobotrya japonica Lindl.) is divided into yellow- and white-fleshed based on the difference in fruit color, and the variations in carotenoids accumulation are considered as the main reasons for this difference. Using RNA-seq technology, a transcriptome analysis was carried out on the flesh and peel of ‘Baiyu’ fruit during four different fruit development stages. A total of 172.53 Gb clean reads with an average of 6.33 Gb reads were detected for each library, and the percentage of Q30 was higher than 90.84%. We identified 16 carotenogenic and 13 plastid-lipid-associated protein (PAP) genes through RNA-seq. Of these, five carotenogenic and four PAP related genes exhibited remarkable differences in the expression patterns. Carotenoids biosynthetic genes, including DXS, PSY1 and VDE displayed higher expression levels in peel than that in the flesh. However, carotenoids decomposition gene, such as NCDE1, exhibited higher expression in flesh than that in the peel. Notably, all differentially expressed PAP genes showed higher expression levels in peel than flesh. We inferred that the differential accumulation of carotenoids in flesh and peel of 'Baiyu' is caused by the up- or down-regulation of the carotenogenic and PAP related genes. The functional analysis of these important genes will provide valuable information about underlying molecular mechanism of carotenoids accumulation in loquat.

Genome-wide identification and characterization of the ALOG gene family in Petunia

BACKGROUND

The ALOG (Arabidopsis LSH1 and Oryza G1) family of proteins, namely DUF640 (domain of unknown function 640) domain proteins, were found in land plants. Functional characterization of a few ALOG members in model plants such as Arabidopsis and rice suggested they play important regulatory roles in plant development. The information about its evolution, however, is largely limited, and there was no any report on the ALOG genes in Petunia, an important ornamental species.

RESULTS

The ALOG genes were identified in four species of Petunia including P. axillaris, P. inflata, P. integrifolia, and P. exserta based on the genome and/or transcriptome databases, which were further confirmed by cloning from P. hybrida 'W115' (Mitchel diploid), a popular laboratorial petunia line susceptible to genetic transformation. Phylogenetic analysis indicated that Petunia ALOG genes (named as LSHs according to their closest Arabidopsis homologs) were grouped into four clades, which can be further divided into eight groups, and similar exon-intron structure and motifs are reflected in the same group. The PhLSH genes of hybrid petunia 'W115' were mainly derived from P. axillaris. The qPCR analysis revealed distinct spatial expression patterns among them suggesting potentially functional diversification. Moreover, over-expressing PhLSH7a and PhLSH7b in Arabidopsis uncovered their functions in the development of both vegetative and reproductive organs.

CONCLUSIONS

Petunia genome includes 11 ALOG genes that can be divided into eight distinct groups, and they also show different expression patterns. Among these genes, PhLSH7b and PhLSH7a play significant roles in plant growth and development, especially in fruit development. Our results provide new insight into the evolution of ALOG gene family and have laid a good foundation for the study of petunia LSH gene in the future.