Evolutionary Analyses Reveal Diverged Patterns of SQUAMOSA Promoter Binding Protein-Like (SPL) Gene Family in Oryza Genus

Genome-wide identification and expression analysis of the SPL transcription factor family and its response to abiotic stress in Pisum sativum L

Squamous promoter binding protein-like (SPL) genes encode plant-specific transcription factors (TFs) that play essential roles in modulating plant growth, development, and stress response. Pea (Pisum sativum L.) is a coarse grain crop of great importance in food production, biodiversity conservation and molecular genetic research, providing genetic information and nutritional resources for improving agricultural production and promoting human health. However, only limited researches on the structure and functions of SPL genes exist in pea (PsSPLs). In this study, we identified 22 PsSPLs and conducted a genome-wide analysis of their physical characteristics, chromosome distribution, gene structure, phylogenetic evolution and gene expression patterns. As a result, the PsSPLs were unevenly distributed on the seven chromosomes of pea and harbored the SBP domain, which is composed of approximately 76 amino acid residues. The phylogenetic analysis revealed that the PsSPLs clustered into eight subfamilies and showed high homology with SPL genes in soybean. Further analysis showed the presence of segmental duplications in the PsSPLs. The expression patterns of 22 PsSPLs at different tissues, developmental stages and under various stimulus conditions were evaluated by qRT-PCR method. It was found that the expression patterns of PsSPLs from the same subfamily were similar in different tissues, the transcripts of most PsSPLs reached the maximum peak value at 14 days after anthesis in the pod. Abiotic stresses can cause significantly up-regulated PsSPL19 expression with spatiotemporal specificity, in addition, four plant hormones can cause the up-regulated expression of most PsSPLs including PsSPL19 in a time-dependent manner. Therefore, PsSPL19 could be a key candidate gene for signal transduction during pea growth and development, pod formation, abiotic stress and plant hormone response. Our findings should provide insights for the elucidating of development regulation mechanism and breeding for resistance to abiotic stress pea.

Genome-wide identification, stress- and hormone-responsive expression characteristics, and regulatory pattern analysis of Scutellaria baicalensis SbSPLs

SQUAMOSA PROMOTER BINDING PROTEIN-LIKEs (SPLs) encode plant-specific transcription factors that regulate plant growth and development, stress response, and metabolite accumulation. However, there is limited information on Scutellaria baicalensis SPLs. In this study, 14 SbSPLs were identified and divided into 8 groups based on phylogenetic relationships. SbSPLs in the same group had similar structures. Abscisic acid-responsive (ABRE) and MYB binding site (MBS) cis-acting elements were found in the promoters of 8 and 6 SbSPLs. Segmental duplications and transposable duplications were the main causes of SbSPL expansion. Expression analysis based on transcriptional profiling showed that SbSPL1, SbSPL10, and SbSPL13 were highly expressed in roots, stems, and flowers, respectively. Expression analysis based on quantitative real-time polymerase chain reaction (RT‒qPCR) showed that most SbSPLs responded to low temperature, drought, abscisic acid (ABA) and salicylic acid (SA), among which the expression levels of SbSPL7/9/10/12 were significantly upregulated in response to abiotic stress. These results indicate that SbSPLs are involved in the growth, development and stress response of S. baicalensis. In addition, 8 Sba-miR156/157 s were identified, and SbSPL1-5 was a potential target of Sba-miR156/157 s. The results of target gene prediction and coexpression analysis together indicated that SbSPLs may be involved in the regulation of L-phenylalanine (L-Phe), lignin and jasmonic acid (JA) biosynthesis. In summary, the identification and characterization of the SbSPL gene family lays the foundation for functional research and provides a reference for improved breeding of S. baicalensis stress resistance and quality traits.

Identification and expression analysis of SBP-Box-like (SPL) gene family disclose their contribution to abiotic stress and flower budding in pigeon pea (Cajanus cajan)

The SPL gene family (for Squamosa Promoter-binding like Proteins) represents specific transcription factors that have significant roles in abiotic stress tolerance, development and the growth processes of different plants, including initiation of the leaf, branching and development of shoot and fruits. The SPL gene family has been studied in different plant species; however, its role is not yet fully explored in pigeon pea (Cajanus cajan ). In the present study, 11 members of the CcSPL gene family were identified in C. cajan . The identified SPLs were classified into nine groups based on a phylogenetic analysis involving SPL protein sequences from C. cajan , Arabidopsis thaliana , Cicer arietinum , Glycine max , Phaseolus vulgaris , Vigna unguiculata and Arachis hypogaea . Further, the identification of gene structure, motif analysis, domain analysis and presence of cis -regulatory elements in the SPL family members were studied. Based on RNA-sequencing data, gene expression analysis was performed, revealing that CcSPL2.1, 3 and 13A were significantly upregulated for salt-tolerance and CcSPL14 and 15 were upregulated in a salt-susceptible cultivar. Real-time qPCR validation indicated that CcSPL3, 4, 6 and 13A were upregulated under salt stress conditions. Therefore, molecular docking was performed against the proteins of two highly expressed genes (CcSPL3 and CcSPL14 ) with three ligands: abscisic acid, gibberellic acid and indole-3-acetic acid. Afterward, their binding affinity was obtained and three-dimensional structures were predicted. In the future, our study may open avenues for harnessing CcSPL genes in pigeon pea for enhanced abiotic stress resistance and developmental traits.

Genome-wide identification, evolution, and role of SPL gene family in beet (Beta vulgaris L.) under cold stress

BACKGROUND

SPL transcription factors play vital roles in regulating plant growth, development, and abiotic stress responses. Sugar beet (Beta vulgaris L.), one of the world's main sugar-producing crops, is a major source of edible and industrial sugars for humans. Although the SPL gene family has been extensively identified in other species, no reports on the SPL gene family in sugar beet are available.

RESULTS

Eight BvSPL genes were identified at the whole-genome level and were renamed based on their positions on the chromosome. The gene structure, SBP domain sequences, and phylogenetic relationship with Arabidopsis were analyzed for the sugar beet SPL gene family. The eight BvSPL genes were divided into six groups (II, IV, V, VI, VII, and VIII). Of the BvSPL genes, no tandem duplication events were found, but one pair of segmental duplications was present. Multiple cis-regulatory elements related to growth and development were identified in the 2000-bp region upstream of the BvSPL gene start codon (ATG). Using quantitative real-time polymerase chain reaction (qRT-PCR), the expression profiles of the eight BvSPL genes were examined under eight types of abiotic stress and during the maturation stage. BvSPL transcription factors played a vital role in abiotic stress, with BvSPL3 and BvSPL6 being particularly noteworthy.

CONCLUSION

Eight sugar beet SPL genes were identified at the whole-genome level. Phylogenetic trees, gene structures, gene duplication events, and expression profiles were investigated. The qRT-PCR analysis indicated that BvSPLs play a substantial role in the growth and development of sugar beet, potentially participating in the regulation of root expansion and sugar accumulation.

Comparative analysis of SPL transcription factors from streptophyte algae and embryophytes reveals evolutionary trajectories of SPL family in streptophytes

SQUAMOSA-PROMOTER BINDING PROTEIN-LIKE (SPL) genes encode plant-specific transcription factors which are important regulators of diverse plant developmental processes. We took advantage of available genome sequences of streptophyte algae representatives to investigate the relationships of SPL genes between freshwater green algae and land plants. Our analysis showed that streptophyte algae, hornwort and liverwort genomes encode from one to four SPL genes which is the smallest set, in comparison to other land plants studied to date. Based on the phylogenetic analysis, four major SPL phylogenetic groups were distinguished with Group 3 and 4 being sister to Group 1 and 2. Comparative motif analysis revealed conserved protein motifs within each phylogenetic group and unique bryophyte-specific motifs within Group 1 which suggests lineage-specific protein speciation processes. Moreover, the gene structure analysis also indicated the specificity of each by identifying differences in exon-intron structures between the phylogenetic groups, suggesting their evolutionary divergence. Since current understanding of SPL genes mostly arises from seed plants, the presented comparative and phylogenetic analyzes from freshwater green algae and land plants provide new insights on the evolutionary trajectories of the SPL gene family in different classes of streptophytes.

Transcriptomic responses to cold stress in Dendrobium huoshanense C.Z. Tang et S.J. Cheng

Dendrobium huoshanense C.Z. Tang et S.J. Cheng is a perennial epiphytic herb of the family Orchidaceae. The main metabolites of D. huoshanense include polysaccharides and flavonoids. Low temperature is the main environmental factor that limits the growth and development of plants. However, changes that occur at the molecular level in response to low temperatures in D. huoshanense are poorly understood. We performed a transcriptome analysis at two time points of 0 d (control group) and 7 d (cold stress group) under culture of D. huoshanense at 4 °C. A total of 37.63 Gb transcriptomic data were generated using the MGI 2000 platform. These reads were assembled into 170,754 transcripts and 23,724 differentially expressed genes (DEGs) were obtained. Pathway analysis indicated that "flavonoid biosynthesis," "anthocyanin biosynthesis," "flavone and flavonol biosynthesis," and "plant hormone signal transduction" might play a vital role in the response of D. huoshanense to cold stress. Several important pathway genes were identified to be altered under cold stress, such as genes encoding polysaccharides, flavonoids, and plant hormone-signaling transduction kinase. In addition, the content of mannose and total flavonoids increased under cold stress. Twelve DEGs related to polysaccharides, flavonoid, and hormone pathways were selected from the transcriptome data for validation with real-time quantitative PCR (RT-qPCR). Our results provide a transcriptome database and candidate genes for further study of the response of D. huoshanense to cold stress.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-023-01385-7.

Genome-Wide Association Analysis for Resistance to Coniothyrium glycines Causing Red Leaf Blotch Disease in Soybean

Soybean is a high oil and protein-rich legume with several production constraints. Globally, several fungi, viruses, nematodes, and bacteria cause significant yield losses in soybean. Coniothyrium glycines (CG), the causal pathogen for red leaf blotch disease, is the least researched and causes severe damage to soybean. The identification of resistant soybean genotypes and mapping of genomic regions associated with resistance to CG is critical for developing improved cultivars for sustainable soybean production. This study used single nucleotide polymorphism (SNP) markers generated from a Diversity Arrays Technology (DArT) platform to conduct a genome-wide association (GWAS) analysis of resistance to CG using 279 soybean genotypes grown in three environments. A total of 6395 SNPs was used to perform the GWAS applying a multilocus model Fixed and random model Circulating Probability Unification (FarmCPU) with correction of the population structure and a statistical test p-value threshold of 5%. A total of 19 significant marker-trait associations for resistance to CG were identified on chromosomes 1, 5, 6, 9, 10, 12, 13, 15, 16, 17, 19, and 20. Approximately 113 putative genes associated with significant markers for resistance to red leaf blotch disease were identified across soybean genome. Positional candidate genes associated with significant SNP loci-encoding proteins involved in plant defense responses and that could be associated with soybean defenses against CG infection were identified. The results of this study provide valuable insight for further dissection of the genetic architecture of resistance to CG in soybean. They also highlight SNP variants and genes useful for genomics-informed selection decisions in the breeding process for improving resistance traits in soybean.

Genome-Wide Identification and Evolution Analysis of the CYP76 Subfamily in Rice (Oryza sativa)

The CYP76 subfamily, a member of the CYP superfamily, plays crucial roles in the biosynthesis of phytohormones in plants, involving biosynthesis of secondary metabolites, hormone signaling, and response to environmental stresses. Here, we conducted a genome-wide analysis of the CYP76 subfamily in seven AA genome species: Oryza sativa ssp. japonica, Oryza sativa ssp. indica, Oryza rufipogon, Oryza glaberrima, Oryza meridionalis, Oryza barthii, and Oryza glumaepatula. These were identified and classified into three groups, and it was found that Group 1 contained the largest number of members. Analysis of cis-acting elements revealed a large number of elements related to jasmonic acid and light response. The gene duplication analysis revealed that the CYP76 subfamily expanded mainly in SD/WGD and tandem forms and underwent strong purifying selection during evolution. Expression pattern analysis of OsCYP76s in various developmental stages revealed that the majority of OsCYP76s exhibit relatively restricted expression patterns in leaves and roots. We further analyzed the expression of CYP76s in O. sativa, japonica, and O. sativa, indica under cold, flooding, drought, and salt abiotic stresses by qRT-PCR. We found that OsCYP76-11 showed a huge increase in relative expression after drought and salt stresses. After flooding stress, OsiCYP76-4 showed a greater increase in expression compared to other genes. CYP76 in japonica and indica showed different response patterns to the same abiotic stresses, revealing functional divergence in the gene family during evolution; these may be the key genes responsible for the differences in tolerance to indica japonica. Our results provide valuable insights into the functional diversity and evolutionary history of the CYP76 subfamily and pave the way for the development of new strategies for improving stress tolerance and agronomic traits in rice.

Zinc oxide nanoparticles (ZnONPs) influenced seed development, grain quality, and remobilization by affecting the transcription of microRNA 171 (miR171), miR156, NAM, and SUT genes in wheat (Triticum aestivum): a biological advantage and risk assessment study

Limited studies have been conducted on the role of microRNAs (miRs) and transcription factors in regulating plant cell responses to nanoparticles. This study attempted to address whether the foliar application of zinc oxide nanoparticles (ZnONPs; 0, 10, 25, and 50 mgL^-1) can affect miRs, gene expression, and wheat grain quality. The seedlings were sprayed with ZnONPs (0, 10, 25, and 50 mgL^-1) or bulk counterpart (BZnO) five times at 72 h intervals. The application of ZnONPs at 10 mgL^-1 increased the number of spikelets and seed weight, while the nano-supplement at 50 mgL^-1 was accompanied by severe restriction on developing spikes and grains. ZnONPs, in a dose-dependent manner, transcriptionally influenced miR156 and miR171. The expression of miR171 showed a similar trend to that of miR156. The ZnONPs at optimum concentration upregulated the NAM transcription factor and sucrose transporter (SUT) at transcriptional levels. However, the transcription of both NAM and SUT genes displayed a downward trend in response to the toxic dose of ZnONPs (50 mgL^-1). Utilization of ZnONPs increased proline and total soluble phenolic content. Monitoring the accumulation of carbohydrates, including fructan, glucose, fructose, and sucrose, revealed that ZnONPs at 10 mgL^-1 modified the source/sink communication and nutrient remobilization. The molecular and physiological data revealed that the expression of miR156 and miR171 is tightly linked to seed grain development, remobilization of carbohydrates, and genes involved in nutrient transportation. This study establishes a novel strategy for obtaining higher yields in crops. This biological risk assessment investigation also displays the potential hazard of applying ZnONPs at the flowering developmental phase.

Genome-Wide Identification of DUF668 Gene Family and Expression Analysis under Drought and Salt Stresses in Sweet Potato [Ipomoea batatas (L.) Lam]

The domain of unknown function 668 (DUF668) is a gene family that plays a vital role in responses to adversity coercion stresses in plant. However, the function of the DUF668 gene family is not fully understood in sweet potato. In this study, bioinformatics methods were used to analyze the number, physicochemical properties, evolution, structure, and promoter cis-acting elements of the IbDUF668 family genes, and RNA-seq and qRT-PCR were performed to detect gene expression and their regulation under hormonal and abiotic stress. A total of 14 IbDUF668 proteins were identified in sweet potato, distributed on nine chromosomes. By phylogenetic analysis, IbDUF668 proteins can be divided into two subfamilies. Transcriptome expression profiling revealed that many genes from DUF668 in sweet potato showed specificity and differential expression under cold, heat, drought, salt and hormones (ABA, GA3 and IAA). Four genes (IbDUF668-6, 7, 11 and 13) of sweet potato were significantly upregulated by qRT-PCR under ABA, drought and NaCl stress. Results suggest that the DUF668 gene family is involved in drought and salt tolerance in sweet potato, and it will further provide the basic information of DUF668 gene mechanisms in plants.

Genome-Wide Analysis of SQUAMOSA-Promoter-Binding Protein-like Family in Flowering Pleioblastus pygmaeus

SQUAMOSA Promoter-Binding Protein-Like (SPL) family is well-known for playing an important role in plant growth and development, specifically in the reproductive process. Bamboo plants have special reproductive characteristics with a prolonged vegetative phase and uncertain flowering time. However, the underlying functions of SPL genes in reproductive growth are undisclosed in bamboo plants. In the study, a total of 28 SPLs were screened from an ornamental dwarf bamboo species, Pleioblastus pygmaeus. Phylogenetic analysis indicates that 183 SPLs from eight plant species can be classified into nine subfamilies, and the 28 PpSPLs are distributed among eight subfamilies. Homologous analysis shows that as many as 32 pairs of homologous genes were found between P. pygmaeus and rice, and 83 pairs were found between P. pygmaeus and Moso bamboo, whose Ka/Ks values are all <1. MiRNA target prediction reveals that 13 out of the 28 PpSPLs have recognition sites complementary to miRNA156. To screen the SPLs involved in the reproductive growth of bamboo plants, the mRNA abundance of the 28 PpSPLs was profiled in the different tissues of flowering P. pygmaeus and non-flowering plants by RNA-Seq. Moreover, the relative expression level of eight PpSPLs is significantly higher in flowering P. pygmaeus than that in non-flowering plants, which was also validated by RT-qPCR. Combined with phylogenetic analysis and homologous analysis, the eight significant, differentially expressed PpSPLs were identified to be associated with the reproductive process and flower organ development. Among them, there are four potential miRNA156-targeting PpSPLs involved in the flowering process. Of significant interest in the study is the identification of 28 SPLs and the exploration of four key flowering-related SPLs from P. pygmaeus, which provides a theoretic basis for revealing the underlying functions of SPLs in the reproductive growth of bamboo plants.

Analysis of Homologous Regions of Small RNAs MIR397 and MIR408 Reveals the Conservation of Microsynteny among Rice Crop-Wild Relatives

MIRNAs are small non-coding RNAs that play important roles in a wide range of biological processes in plant growth and development. MIR397 (involved in drought, low temperature, and nitrogen and copper (Cu) starvation) and MIR408 (differentially expressed in response to environmental stresses such as copper, light, mechanical stress, dehydration, cold, reactive oxygen species, and drought) belong to conserved MIRNA families that either negatively or positively regulate their target genes. In the present study, we identified the homologs of MIR397 and MIR408 in Oryza sativa and its six wild progenitors, three non-Oryza species, and one dicot species. We analyzed the 100 kb segments harboring MIRNA homologs from 11 genomes to obtain a comprehensive view of their community evolution around these loci in the farthest (distant) relatives of rice. Our study showed that mature MIR397 and MIR408 were highly conserved among all Oryza species. Comparative genomics analyses also revealed that the microsynteny of the 100 kb region surrounding MIRNAs was only conserved in Oryza spp.; disrupted in Sorghum, maize, and wheat; and completely lost in Arabidopsis. There were deletions, rearrangements, and translocations within the 100 kb segments in Oryza spp., but the overall microsynteny of the region was maintained. The phylogenetic analyses of the precursor regions of all MIRNAs under study revealed a bimodal clade of common origin. This comparative analysis of miRNA involved in abiotic stress tolerance in plants provides a powerful tool for future Oryza research. Crop wild relatives (CWRs) offer multiple traits with potential to decrease the amount of yield loss owing to biotic and abiotic stresses. Using a comparative genomics approach, the exploration of CWRs as a source of tolerance to these stresses by understanding their evolution can be further used to leverage their yield potential.

Genome-wide analysis of serine carboxypeptidase-like protein (SCPL) family and functional validation of Gh_SCPL42 unchromosome conferring cotton Verticillium der Verticillium wilt stress in Gossypium hirsutum

BACKGROUND

Serine carboxypeptidase-like protein (SCPL) plays an important role in response to stress in plant. However, our knowledge of the function of the SCPL gene family is limited.

RESULTS

In this study, a comprehensive and systematic analysis of SCPL gene family was conducted to explore the phylogeny and evolution of the SCPL gene in Gossypium hirsutum. The phenotype and molecular mechanism of silencing of the Gh_SCPL42 under Verticillium wilt stress was also studied. Our results showed that 96 SCPL genes were observed in genome of G. hirsutum, which distributed on 25 chromosomes and most of them were located in the nucleus. The phylogenetic tree analysis showed that members of SCPL gene family can be divided into three subgroups in G. hirsutum, which are relatively conservative in evolution. SCPL gene has a wide range of tissue expression types in G. hirsutum. Promoter analysis showed that the most cis-acting elements related to MeJA and ABA were contained. Through RNA-seq combined with genotyping, it was found that 11 GhSCPL genes not only had significant expression changes during Verticillium wilt stress but also had differential SNPs in the upstream, downstream, exonic or intronic regions. The expression of these 11 genes in the resistant (Zhongzhimian 2) and susceptible (Junmian 1) materials was further analyzed by qRT-PCR, it was found that 6 genes showed significant expression differences in the two materials. Among them, Gh_SCPL42 has the most obvious expression change. Furthermore, virus-induced gene silencing (VIGS) showed necrosis and yellowing of leaves and significantly higher disease severity index (DSI) and disease severity rate (DSR) values in VIGS plants than in control silenced Gh_SCPL42 plants. Moreover, the expression levels of genes related to the SA and JA pathways were significantly downregulated. These results show that Gh_SCPL42 might improve resistance to Verticillium wilt through the SA and JA pathways in G. hirsutum.

CONCLUSION

In conclusion, our findings indicated that Gh_SCPL42 gene plays an important role in resistance to Verticillium wilt in cotton. It was provided an important theoretical basis for further research on the function of SCPL gene family and the molecular mechanism of resistance to Verticillium wilt in cotton.

Genome-wide identification and characterization of the SPL gene family and its expression in the various developmental stages and stress conditions in foxtail millet (Setaria italica)

Background

Among the major transcription factors, SPL plays a crucial role in plant growth, development, and stress response. Foxtail millet (Setaria italica), as a C4 crop, is rich in nutrients and is beneficial to human health. However, research on the foxtail millet SPL (SQUAMOSA PROMOTER BINDING-LIKE) gene family is limited. 

Results

In this study, a total of 18 SPL genes were identified for the comprehensive analysis of the whole genome of foxtail millet. These SiSPL genes were divided into seven subfamilies (I, II, III, V, VI, VII, and VIII) according to the classification of the Arabidopsis thaliana SPL gene family. Structural analysis of the SiSPL genes showed that the number of introns in subfamilies I and II were much larger than others, and the promoter regions of SiSPL genes were rich in different cis-acting elements. Among the 18 SiSPL genes, nine genes had putative binding sites with foxtail millet miR156. No tandem duplication events were found between the SiSPL genes, but four pairs of segmental duplications were detected. The SiSPL genes expression were detected in different tissues, which was generally highly expressed in seeds development process, especially SiSPL6 and SiSPL16, which deserve further study. The results of the expression levels of SiSPL genes under eight types of abiotic stresses showed that many stress responsive genes, especially SiSPL9, SiSPL10, and SiSPL16, were highly expressed under multiple stresses, which deserves further attention.

Conclusions

In this research, 18 SPL genes were identified in foxtail millet, and their phylogenetic relationships, gene structural features, duplication events, gene expression and potential roles in foxtail millet development were studied. The findings provide a new perspective for the mining of the excellent SiSPL gene and the molecular breeding of foxtail millet.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08633-2.

Genome-Wide Identification, Characterization, and Expression Profiling Analysis of SPL Gene Family during the Inflorescence Development in Trifolium repens

Trifolium repens is the most widely cultivated perennial legume forage in temperate region around the world. It has rich nutritional value and good palatability, seasonal complementarity with grasses, and can improve the feed intake and digestibility of livestock. However, flowering time and inflorescence development directly affects the quality and yield of T. repens, as well as seed production. The Squa promoter binding protein-like (SPL) gene family is a plant specific transcription factor family, which has been proved to play a critical role in regulating plant formation time and development of flowers. In this study, a total of 37 TrSPL genes were identified from the whole genome of T. repens and were divided into nine clades based on phylogenetic tree. Seventeen TrSPL genes have potential target sites for miR156. The conserved motif of squamosa promoter binding protein (SBP) contains two zinc finger structures and one NLS structure. Gene structure analysis showed that all TrSPL genes contained SBP domain, while ankyrin repeat region was just distributed in part of genes. 37 TrSPL genes were relatively dispersedly distributed on 16 chromosomes, and 5 pairs of segmental repeat genes were found, which indicated that segmental duplication was the main way of gene expansion. Furthermore, the gene expression profiling showed that TrSPL11, TrSPL13, TrSPL22, and TrSPL26 were highly expressed only in the early stage of inflorescence development, while TrSPL1 and TrSPL6 are highly expressed only in the mature inflorescence. Significantly, the expression of TrSPL4 and TrSPL12 increased gradually with the development of inflorescences. The results of this study will provide valuable clues for candidate gene selection and elucidating the molecular mechanism of T. repens flowering regulation.

Genome-Wide Identification and Evolution Analysis of the Gibberellin Oxidase Gene Family in Six Gramineae Crops

The plant hormones gibberellins (GAs) regulate plant growth and development and are closely related to the yield of cash crops. The GA oxidases (GAoxs), including the GA2ox, GA3ox, and GA20ox subfamilies, play pivotal roles in GAs' biosynthesis and metabolism, but their classification and evolutionary pattern in Gramineae crops remain unclear. We thus conducted a comparative genomic study of GAox genes in six Gramineae representative crops, namely, Setaria italica (Si), Zea mays (Zm), Sorghum bicolor (Sb), Hordeum vulgare (Hv), Brachypodium distachyon (Bd), and Oryza sativa (Os). A total of 105 GAox genes were identified in these six crop genomes, belonging to the C19-GA2ox, C20-GA2ox, GA3ox, and GA20ox subfamilies. Based on orthogroup (OG) analysis, GAox genes were divided into nine OGs and the number of GAox genes in each of the OGs was similar among all tested crops, which indicated that GAox genes may have completed their family differentiations before the species differentiations of the tested species. The motif composition of GAox proteins showed that motifs 1, 2, 4, and 5, forming the 2OG-FeII_Oxy domain, were conserved in all identified GAox protein sequences, while motifs 11, 14, and 15 existed specifically in the GA20ox, C19-GA2ox, and C20-GA2ox protein sequences. Subsequently, the results of gene duplication events suggested that GAox genes mainly expanded in the form of WGD/SD and underwent purification selection and that maize had more GAox genes than other species due to its recent duplication events. The cis-acting elements analysis indicated that GAox genes may respond to growth and development, stress, hormones, and light signals. Moreover, the expression profiles of rice and maize showed that GAox genes were predominantly expressed in the panicles of the above two plants and the expression of several GAox genes was significantly induced by salt or cold stresses. In conclusion, our results provided further insight into GAox genes' evolutionary differences among six representative Gramineae and highlighted GAox genes that may play a role in abiotic stress.

Using iRNA-seq analysis to predict gene expression regulatory level and activity in Zea mays tissues

Plants regulate gene expression at the transcriptional and post-transcriptional levels to produce a variety of functionally diverse cells and tissues that ensure normal growth, development, and environmental response. Although distinct gene expression patterns have been characterized between different plant tissues, the specific role of transcriptional regulation of tissue-specific expression is not well-characterized in plants. RNA-seq, while widely used to assay for changes in transcript abundance, does not discriminate between differential expression caused by mRNA degradation and active transcription. Recently, the presence of intron sequences in RNA-seq analysis of libraries constructed with total RNA has been found to coincide with genes undergoing active transcription. We have adapted the intron RNA-sequencing analysis to determine genome-wide transcriptional activity in 2 different maize (Zea mays) tissues: husk and V2-inner stem tissue. A total of 5,341 genes were predicted to be transcriptionally differentially expressed between the 2 tissues, including many genes expected to have biological activity relevant to the functional and developmental identity of each tissue. Correlations with transcriptional enhancer and transcription factor activity support the validity of intron RNA-sequencing predictions of transcriptional regulation. A subset of transcription factors was further analyzed using gene regulatory network analysis to determine the possible impact of their activation. The predicted regulatory patterns between these genes were used to model a potential gene regulatory network of transcription factors and regulatory targets.

Osa-miR535 targets SQUAMOSA promoter binding protein-like 4 to regulate blast disease resistance in rice

Many rice microRNAs have been identified as fine-tuning factors in the regulation of agronomic traits and immunity. Among them, Osa-miR535 targets SQUAMOSA promoter binding protein-like 14 (OsSPL14) to positively regulate tillers but negatively regulate yield and immunity. Here, we uncovered that Osa-miR535 targets another SPL gene, OsSPL4, to suppress rice immunity against Magnaporthe oryzae. Overexpression of Osa-miR535 significantly decreased the accumulation of the fusion protein SPL4_TBS -YFP that contains the target site of Osa-miR535 in OsSPL4. Consistently, Osa-miR535 mediated the cleavage of OsSPL4 mRNA between the 10th and 11th base pair of the predicted binding site at the 3' untranslated region. Transgenic rice lines overexpressing OsSPL4 (OXSPL4) displayed enhanced blast disease resistance accompanied by enhanced immune responses, including increased expression of defense-relative genes and up-accumulated H₂ O₂ . By contrast, the knockout mutant osspl4 exhibited susceptibility. Moreover, OsSPL4 binds to the promoter of GH3.2, an indole-3-acetic acid-amido synthetase, and promotes its expression. Together, these data indicate that Os-miR535 targets OsSPL4 and OsSPL4-GH3.2, which may parallel the OsSPL14-WRKY45 module in rice blast disease resistance.

OsSPLs Regulate Male Fertility in Response to Different Temperatures by Flavonoid Biosynthesis and Tapetum PCD in PTGMS Rice

Photoperiod and thermo-sensitive genic male sterile (PTGMS) rice is an important resource for two line hybrid rice production. The SQUAMOSA–promoter binding, such as the (SPL) gene family, encode the plant specific transcription factors that regulate development and defense responses in plants. However, the reports about SPLs participating in male fertility regulation are limited. Here, we identified 19 OsSPL family members and investigated their involvement in the fertility regulation of the PTGMS rice lines, PA2364S and PA2864S, with different fertility transition temperatures. The results demonstrated that OsSPL2, OsSPL4, OsSPL16 and OsSPL17 affect male fertility in response to temperature changes through the MiR156-SPL module. WGCNA (weighted gene co-expression network analysis) revealed that CHI and APX1 were co-expressed with OsSPL17. Targeted metabolite and flavonoid biosynthetic gene expression analysis revealed that OsSPL17 regulates the expression of flavonoid biosynthesis genes CHI, and the up regulation of flavanones (eriodictvol and naringenin) and flavones (apigenin and luteolin) content contributed to plant fertility. Meanwhile, OsSPL17 negatively regulates APX1 to affect APX (ascorbate peroxidase) activity, thereby regulating ROS (reactive oxygen species) content in the tapetum, controlling the PCD (programmed cell death) process and regulating male fertility in rice. Overall, this report highlights the potential role of OsSPL for the regulation of male fertility in rice and provides a new insight for the further understanding of fertility molecular mechanisms in PTGMS rice.

Identification, evolutionary analysis and functional diversification of RAV gene family in cotton (G. hirsutum L.)

Genome wide analysis, expression pattern analysis, and functional characterization of RAV genes highlight their roles in roots, stem development and hormonal response. RAV (Related to ABI3 and VP1) gene family members have been involved in tissues/organs growth and hormone signaling in various plant species. Here, we identified 247 RAVs from 12 different species with 33 RAV genes from G. hirsutum. Phylogenetic analysis classified RAV genes into four distinct groups. Analysis of gene structure showed that most GhRAVs lack introns. Motif distribution pattern and protein sequence logos indicated that GhRAV genes were highly conserved during the process of evolution. Promotor cis-acting elements revealed that promotor regions of GhRAV genes encode numerous elements related to plant growth, abiotic stresses and phytohormones. Chromosomal location information showed uneven distribution of 33 GhRAV genes on different chromosomes. Collinearity analysis identified 628 and 52 orthologous/ paralogous gene pairs in G. hirsutum and G. barbadense, respectively. Ka/Ks values indicated that GhRAV and GbRAV genes underwent strong purifying selection pressure. Selecton model and codon model selection revealed that GhRAV amino acids were under purifying selection and adaptive evolution exists among GhRAV proteins. Three dimensional structure of GhRAVs indicated the presence of numerous alpha helix and beta-barrels. Expression level revealed that some GhRAV genes exhibited high expression in roots (GhRAV3, GhRAV4, GhRAV11, GhRAV18, GhRAV20 and GhRAV30) and stem (GhRAV3 and GhRAV18), indicating their potential role in roots and stem development. GhRAV genes can be regulated by phytohormonal stresses (BL, JA and IAA). Our study provides a reference for future studies related to the functional analysis of GhRAVs in cotton.

The Elite Alleles of OsSPL4 Regulate Grain Size and Increase Grain Yield in Rice

Grain weight and grain number, the two important yield traits, are mainly determined by grain size and panicle architecture in rice. Herein, we report the identification and functional analysis of OsSPL4 in panicle and grain development of rice. Using CRISPR/Cas9 system, two elite alleles of OsSPL4 were obtained, which exhibited an increasing number of grains per panicle and grain size, resulting in increase of rice yield. Cytological analysis showed that OsSPL4 could regulate spikelet development by promoting cell division. The results of RNA-seq and qRT-PCR validations also demonstrated that several MADS-box and cell-cycle genes were up-regulated in the mutation lines. Co-expression network revealed that many yield-related genes were involved in the regulation network of OsSPL4. In addition, OsSPL4 could be cleaved by the osa-miR156 in vivo, and the OsmiR156-OsSPL4 module might regulate the grain size in rice. Further analysis indicated that the large-grain allele of OsSPL4 in indica rice might introgress from aus varieties under artificial selection. Taken together, our findings suggested that OsSPL4 could be as a key regulator of grain size by acting on cell division control and provided a strategy for panicle architecture and grain size modification for yield improvement in rice.

Primary Bitter Taste of Citrus is Linked to a Functional Allele of the 1,2-Rhamnosyltransferase Gene Originating from Citrus grandis

1,2-Rhamnosyltransferase (1,2RhaT) catalyzes the final step of production of flavanone neohesperidoside (FNH) that is responsible for the primary bitter taste of citrus fruits. In this study, species-specific flavonoid profiles were determined in 87 Citrus accessions by identifying eight main flavanone glycosides (FGs). Accumulation of FNHs was completely correlated to the presence of the 1,2RhaT gene in 87 citrus accessions analyzed using a novel 1,2RhaT-specific DNA marker. Pummelo (Citrus grandis) was identified as the genetic origin for a function allele of 1,2RhaT that underpinned FNH-bitterness in modern citrus cultivars. In addition, genes encoding six MYB and five bHLH transcription factors were shown to coexpress with 1,2RhaT and other flavonoid pathway genes related to FNH accumulation, indicating that these transcription factors may affect the fruit taste of citrus. This study provides a better understanding of bitterness formation in Citrus varieties and a genetic marker for the early selection of nonbitterness lines in citrus breeding programs.

Multi-locus genome-wide association studies for five yield-related traits in rice

BACKGROUND

Improving the overall production of rice with high quality is a major target of breeders. Mining potential yield-related loci have been geared towards developing efficient rice breeding strategies. In this study, one single-locus genome-wide association studies (SL-GWAS) method (MLM) in conjunction with five multi-locus genome-wide association studies (ML-GWAS) approaches (mrMLM, FASTmrMLM, pLARmEB, pKWmEB, and ISIS EM-BLASSO) were conducted in a panel consisting of 529 rice core varieties with 607,201 SNPs.

RESULTS

A total of 152, 106, 12, 111, and 64 SNPs were detected by the MLM model associated with the five yield-related traits, namely grain length (GL), grain width (GW), grain thickness (GT), thousand-grain weight (TGW), and yield per plant (YPP), respectively. Furthermore, 74 significant quantitative trait nucleotides (QTNs) were presented across at least two ML-GWAS methods to be associated with the above five traits successively. Finally, 20 common QTNs were simultaneously discovered by both SL-GWAS and ML-GWAS methods. Based on genome annotation, gene expression analysis, and previous studies, two candidate key genes (LOC_Os09g02830 and LOC_Os07g31450) were characterized to affect GW and TGW, separately.

CONCLUSIONS

These outcomes will provide an indication for breeding high-yielding rice varieties in the immediate future.

Rice SPL10 positively regulates trichome development through expression of HL6 and auxin-related genes

Trichomes function in plant defenses against biotic and abiotic stresses; examination of glabrous lines, which lack trichomes, has revealed key aspects of trichome development and function. Tests of allelism in 51 glabrous rice (Oryza sativa) accessions collected worldwide identified OsSPL10 and OsWOX3B as regulators of trichome development in rice. Here, we report that OsSPL10 acts as a transcriptional regulator controlling trichome development. Haplotype and transient expression analyses revealed that variation in the approximately 700-bp OsSPL10 promoter region is the primary cause of the glabrous phenotype in the indica cultivar WD-17993. Disruption of OsSPL10 by genome editing decreased leaf trichome density and length in the NIL-HL6 background. Plants with genotype OsSPL10^WD-17993 /HL6 generated by crossing WD-17993 with NIL-HL6 also had fewer trichomes in the glumes. HAIRY LEAF6 (HL6) encodes another transcription factor that regulates trichome initiation and elongation, and OsSPL10 directly binds to the HL6 promoter to regulate its expression. Moreover, the transcript levels of auxin-related genes, such as OsYUCCA5 and OsPIN-FORMED1b, were altered in OsSPL10 overexpression and RNAi transgenic lines. Feeding tests using locusts (Locusta migratoria) demonstrated that non-glandular trichomes affect feeding by this herbivore. Our findings provide a molecular framework for trichome development and an ecological perspective on trichome functions.

Genome-wide identification of the DUF668 gene family in cotton and expression profiling analysis of GhDUF668 in Gossypium hirsutum under adverse stress

Background

Domain of unknown function 668 (DUF668) may play a crucial role in the plant growth and developmental response to adverse stress. However, our knowledge of the function of the DUF668 gene family is limited.

Results

Our study was conducted based on the DUF668 gene family identified from cotton genome sequencing. Phylogenetic analysis showed that the DUF668 family genes can be classified into four subgroups in cotton. We identified 32 DUF668 genes, which are distributed on 17 chromosomes and most of them located in the nucleus of Gossypium hirsutum. Gene structure and motif analyses revealed that the members of the DUF668 gene family can be clustered in G. hirsutum into two broad groups, which are relatively evolutionarily conserved. Transcriptome data analysis showed that the GhDUF668 genes are differentially expressed in different tissues under various stresses (cold, heat, drought, salt, and Verticillium dahliae), and expression is generally increased in roots and stems. Promoter and expression analyses indicated that Gh_DUF668–05, Gh_DUF668–08, Gh_DUF668–11, Gh_DUF668–23 and Gh_DUF668–28 in G. hirsutum might have evolved resistance to adverse stress. Additionally, qRT-PCR revealed that these 5 genes in four cotton lines, KK1543 (drought resistant), Xinluzao 26 (drought sensitive), Zhongzhimian 2 (disease resistant) and Simian 3 (susceptible), under drought and Verticillium wilt stress were all significantly induced. Roots had the highest expression of these 5 genes before and after the treatment. Among them, the expression levels of Gh_DUF668–08 and Gh_DUF668–23 increased sharply at 6 h and reached a maximum at 12 h under biotic and abiotic stress, which showed that they might be involved in the process of adverse stress resistance in cotton.

Conclusion

The significant changes in GhDUF668 expression in the roots after adverse stress indicate that GhDUF668 is likely to increase plant resistance to stress. This study provides an important theoretical basis for further research on the function of the DUF668 gene family and the molecular mechanism of adverse stress resistance in cotton.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07716-w.

Comparative analysis of the SPL gene family in five Rosaceae species: Fragaria vesca, Malus domestica, Prunus persica, Rubus occidentalis, and Pyrus pyrifolia

SQUAMOSA promoter-binding protein-like (SPL) transcription factors are very important for the plant growth and development. Here 15 RoSPLs were identified in Rubus occidentalis. The conserved domains and motifs, phylogenetic relationships, posttranscriptional regulation, and physiological function of the 92 SPL family genes in Fragaria vesca, Malus domestica, Prunus persica, R. occidentalis, and Pyrus pyrifolia were analyzed. Sequence alignment and phylogenetic analysis showed the SPL proteins had sequence conservation, some FvSPLs could be lost or developed, and there was a closer relationship between M. domestica and P. pyrifolia, F. vesca and R. occidentalis, respectively. Genes with similar motifs clustering together in the same group had their functional redundancy. Based on the function of SPLs in Arabidopsis thaliana, these SPLs could be involved in vegetative transition from juvenile to adult, morphological change in the reproductive phase, anthocyanin biosynthesis, and defense stress. Forty-eight SPLs had complementary sequences of miR156, of which nine PrpSPLs in P. persica and eight RoSPLs in R. occidentalis as the potential targets of miR156 were reported for the first time, suggesting the conservative regulatory effects of miR156 and indicating the roles of miR156-SPL modules in plant growth, development, and defense response. It provides a basic understanding of SPLs in Rosaceae plants.

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.

Phylogenetic Analysis of the SQUAMOSA Promoter-Binding Protein-Like Genes in Four Ipomoea Species and Expression Profiling of the IbSPLs During Storage Root Development in Sweet Potato (Ipomoea batatas)

As a major plant-specific transcription factor family, SPL genes play a crucial role in plant growth, development, and stress tolerance. The SPL transcription factor family has been widely studied in various plant species; however, systematic studies on SPL genes in the genus Ipomoea are lacking. Here, we identified a total of 29, 27, 26, and 23 SPLs in Ipomoea batatas, Ipomoea trifida, Ipomoea triloba, and Ipomoea nil, respectively. Based on the phylogenetic analysis of SPL proteins from model plants, the Ipomoea SPLs were classified into eight clades, which included conserved gene structures, domain organizations and motif compositions. Moreover, segmental duplication, which is derived from the Ipomoea lineage-specific whole-genome triplication event, was speculated to have a predominant role in Ipomoea SPL expansion. Particularly, tandem duplication was primarily responsible for the expansion of SPL subclades IV-b and IV-c. Furthermore, 25 interspecific orthologous groups were identified in Ipomoea, rice, Arabidopsis, and tomato. These findings support the expansion of SPLs in Ipomoea genus, with most of the SPLs being evolutionarily conserved. Of the 105 Ipomoea SPLs, 69 were predicted to be the targets of miR156, with seven IbSPLs being further verified as targets using degradome-seq data. Using transcriptomic data from aboveground and underground sweet potato tissues, IbSPLs showed diverse expression patterns, including seven highly expressed IbSPLs in the underground tissues. Furthermore, the expression of 11 IbSPLs was validated using qRT-PCR, and two (IbSPL17/IbSPL28) showed significantly increased expression during root development. Additionally, the qRT-PCR analysis revealed that six IbSPLs were strongly induced in the roots under phytohormone treatments, particularly zeatin and abscisic acid. Finally, the transcriptomic data of storage roots from 88 sweet potato accessions were used for weighted gene co-expression network analysis, which revealed four IbSPLs (IbSPL16/IbSPL17/IbSPL21/IbSPL28) clusters with genes involved in "regulation of root morphogenesis," "cell division," "cytoskeleton organization," and "plant-type cell wall organization or biogenesis," indicating their potential role in storage root development. This study not only provides novel insights into the evolutionary and functional divergence of the SPLs in the genus Ipomoea but also lays a foundation for further elucidation of the potential functional roles of IbSPLs on storage root development.

Characterization of Squamosa-Promoter Binding Protein-Box Family Genes Reveals the Critical Role of MsSPL20 in Alfalfa Flowering Time Regulation

SQUAMOSA Promoter-binding protein-Like (SPL) genes affect a broad range of plant biological processes and show potential application in crop improvement by genetic modification. As the most widely planted forage crop in the world, biomass and abiotic stresses tolerance are important breeding targets for alfalfa (Medicago sativa L.). Nevertheless, the systematic analysis of SPL genes in alfalfa genome remains lacking. In the present study, we characterized 22 putative non-redundant SPL genes in alfalfa genome and uncovered the abundant structural variation among MsSPL genes. The phylogenetic analysis of plant SPL proteins separated them into 10 clades and clade J was an alfalfa-specific clade, suggesting SPL genes in alfalfa might have experienced gene duplication and functional differentiation within the genome. Meanwhile, 11 MsSPL genes with perfect matches to miRNA response elements (MREs) could be degraded by miR156, and the cleavage sites were gene specific. In addition, we investigated the temporal and spatial expression patterns of MsSPL genes and their expression patterns in response to multiple treatments, characterizing candidate SPL genes in alfalfa development and abiotic stress tolerant regulation. More importantly, overexpression of the alfalfa-specific SPL gene (MsSPL20) showed stable delayed flowering time, as well as increased biomass. Further studies indicated that MsSPL20 delayed flowering time by regulating the expression of genes involved in floret development, including HD3A, FTIP1, TEM1, and HST1. Together, our findings provide valuable information for future research and utilization of SPL genes in alfalfa and elucidate a possibly alfalfa-specific flowering time regulation, thereby supplying candidate genes for alfalfa molecular-assisted breeding.

Genome-wide identification, phylogeny and expression analysis of the SPL gene family in wheat

BACKGROUND

Members of the plant-specific SPL gene family (squamosa promoter-binding protein -like) contain the SBP conserved domain and are involved in the regulation of plant growth and development, including the development of plant flowers and plant epidermal hair, the plant stress response, and the synthesis of secondary metabolites. This family has been identified in various plants. However, there is no systematic analysis of the SPL gene family at the genome-wide level of wheat.

RESULTS

In this study, 56 putative TaSPL genes were identified using the comparative genomics method; we renamed them TaSPL001 - TaSPL056 on their chromosomal distribution. According to the un-rooted neighbor joining phylogenetic tree, gene structure and motif analyses, the 56 TaSPL genes were divided into 8 subgroups. A total of 81 TaSPL gene pairs were designated as arising from duplication events and 64 interacting protein branches were identified as involve in the protein interaction network. The expression patterns of 21 randomly selected TaSPL genes in different tissues (roots, stems, leaves and inflorescence) and under 4 treatments (abscisic acid, gibberellin, drought and salt) were detected by quantitative real-time polymerase chain reaction (qRT-PCR).

CONCLUSIONS

The wheat genome contains 56 TaSPL genes and those in same subfamily share similar gene structure and motifs. TaSPL gene expansion occurred through segmental duplication events. Combining the results of transcriptional and qRT-PCR analyses, most of these TaSPL genes were found to regulate inflorescence and spike development. Additionally, we found that 13 TaSPLs were upregulated by abscisic acid, indicating that TaSPL genes play a positive role in the abscisic acid-mediated pathway of the seedling stage. This study provides comprehensive information on the SPL gene family of wheat and lays a solid foundation for elucidating the biological functions of TaSPLs and improvement of wheat yield.

Characterization and Functional Divergence of a Novel DUF668 Gene Family in Rice Based on Comprehensive Expression Patterns

The domain of unknown function (DUF) superfamily encodes proteins of unknown functions in plants. Among them, DUF668 family members in plants possess a 29 amino-acid conserved domain, and this family has not been described previously. Here, we report this plant-specific novel DUF668 gene family containing 12 OsDUF668 genes in rice (Oryza sativa) and 91 DUF668s for the other seven plant species. In our study, DUF668 genes were present in both dicot and monocot plants, indicating that DUF668 is a conserved gene family that originated by predating the dicot–monocot divergence. Based on the gene structure and motif composition, the DUF668 family consists of two distinct clades, I and II in the phylogenetic tree. Remarkably, OsDUF668 genes clustered on the chromosomes merely show close phylogenetic relationships, suggesting that gene duplications or collinearity seldom happened. Cis-elements prediction display that over 80% of DUF668s contain phytohormone and light responsiveness factors. Further comprehensive experimental analyses of the OsDUF668 family are implemented in 22 different tissues, five hormone treatments, seven environmental factor stresses, and two pathogen-defense related stresses. The OsDUF668 genes express ubiquitously in analyzed rice tissues, and seven genes show tissue-specific high expression profiles. All OsDUF668s respond to drought, and some of Avr9/Cf-9 rapidly elicited genes resist to salt, wound, and rice blast with rapidly altered expression patterns. These findings imply that OsDUF668 is essential for drought-enduring and plant defense. Together, our results bring the important role of the DUF668 gene family in rice development and fitness to the fore.