Genome-wide analysis, evolutionary expansion, and expression of early auxin-responsive SAUR gene family in rice (Oryza sativa)

Genome-wide identification of the SAUR gene family and screening for SmSAURs involved in root development in Salvia miltiorrhiza

KEY MESSAGE

SmSAUR4, SmSAUR18, SmSAUR28, SmSAUR37, and SmSAUR38 were probably involved in the auxin-mediated root development in Salvia miltiorrhiza. Salvia miltiorrhiza is a widely utilized medicinal plant in China. Its roots and rhizomes are the main medicinal portions and are closely related to the quality of this herb. Previous studies have revealed that auxin plays pivotal roles in S. miltiorrhiza root development. Whether small auxin-up RNA genes (SAURs), which are crucial early auxin response genes, are involved in auxin-mediated root development in S. miltiorrhiza is worthy of investigation. In this study, 55 SmSAUR genes in S. miltiorrhiza were identified, and their physical and chemical properties, gene structure, cis-acting elements, and evolutionary relationships were analyzed. The expression levels of SmSAUR genes in different organs of S. miltiorrhiza were detected using RNA-seq combined with qRT‒PCR. The root development of S. miltiorrhiza seedlings was altered by the application of indole-3-acetic acid (IAA), and Pearson correlation coefficient analysis was conducted to screen SmSAURs that potentially participate in this physiological process. The diameter of primary lateral roots was positively correlated with SmSAUR4. The secondary lateral root number was positively correlated with SmSAUR18 and negatively correlated with SmSAUR4. The root length showed a positive correlation with SmSAUR28 and SmSAUR37 and a negative correlation with SmSAUR38. The fresh root biomass exhibited a positive correlation with SmSAUR38 and a negative correlation with SmSAUR28. The aforementioned SmSAURs were likely involved in auxin-mediated root development in S. miltiorrhiza. Our study provides a comprehensive overview of SmSAURs and provides the groundwork for elucidating the molecular mechanism underlying root morphogenesis in this species.

A Small Auxin-Up RNA Gene, IbSAUR36, Regulates Adventitious Root Development in Transgenic Sweet Potato

Small auxin-upregulated RNAs (SAURs), as the largest family of early auxin-responsive genes, play important roles in plant growth and development processes, such as auxin signaling and transport, hypocotyl development, and tolerance to environmental stresses. However, the functions of few SAUR genes are known in the root development of sweet potatoes. In this study, an IbSAUR36 gene was cloned and functionally analyzed. The IbSAUR36 protein was localized to the nucleus and plasma membrane. The transcriptional level of this gene was significantly higher in the pencil root and leaf.This gene was strongly induced by indole-3-acetic acid (IAA), but it was downregulated under methyl-jasmonate(MeJA) treatment. The promoter of IbSAUR36 contained the core cis-elements for phytohormone responsiveness. Promoter β-glucuronidase (GUS) analysis in Arabidopsis showed that IbSAUR36 is highly expressed in the young tissues of plants, such as young leaves, roots, and buds. IbSAUR36-overexpressing sweet potato roots were obtained by an efficient Agrobacterium rhizogenes-mediated root transgenic system. We demonstrated that overexpression of IbSAUR36 promoted the accumulation of IAA, upregulated the genes encoding IAA synthesis and its signaling pathways, and downregulated the genes encoding lignin synthesis and JA signaling pathways. Taken together, these results show that IbSAUR36 plays an important role in adventitious root (AR) development by regulating IAA signaling, lignin synthesis, and JA signaling pathways in transgenic sweet potatoes.

Genome-wide identification and expression-pattern analysis of sulfate transporter (SULTR) gene family in cotton under multiple abiotic stresses and fiber development

Sulfate transporter (SULTR) proteins are in charge of the transport and absorption on sulfate substances, and have been reported to play vital roles in the biological processes of plant growth and stress response. However, there were few reports of genome-wide identification and expression-pattern analysis of SULTRs in Hibiscus mutabilis. Gossypium genus is a ideal model for studying the allopolyploidy, therefore two diploid species (G. raimondii and G. arboreum) and two tetraploid species (G. hirsutum and G. barbadense) were chosen in this study to perform bioinformatic analyses, identifying 18, 18, 35, and 35 SULTR members, respectively. All the 106 cotton SULTR genes were utilized to construct the phylogenetic tree together with 11 Arabidopsis thaliana, 13 Oryza sativa, and 8 Zea mays ones, which was divided into Group1-Group4. The clustering analyses of gene structures and 10 conserved motifs among the cotton SULTR genes showed the consistent evolutionary relationship with the phylogenetic tree, and the results of gene-duplication identification among the four representative Gossypium species indicated that genome-wide or segment duplication might make main contributions to the expansion of SULTR gene family in cotton. Having conducted the cis-regulatory element analysis in promoter region, we noticed that the existing salicylic acid (SA), jasmonic acid (JA), and abscisic acid (ABA) elements could have influences with expression levels of cotton SULTR genes. The expression patterns of GhSULTR genes were also investigated on the 7 different tissues or organs and the developing ovules and fibers, most of which were highly expressed in root, stem, sepal, receptacel, ovule at 10 DPA, and fiber at 20 and 25 DPA. In addition, more active regulatory were observed in GhSULTR genes responding to multiple abiotic stresses, and 12 highly expressed genes showed the similar expression patterns in the quantitative Real-time PCR experiments under cold, heat, salt, and drought treatments. These findings broaden our insight into the evolutionary relationships and expression patterns of the SULTR gene family in cotton, and provide the valuable information for further screening the vital candidate genes on trait improvement.

Genome-Wide Analysis of the SAUR Gene Family and Its Expression Profiles in Response to Salt Stress in Santalum album

The SAUR (small auxin-up RNA) family constitutes a category of genes that promptly respond to the hormone auxin and play a pivotal role in diverse biological processes encompassing plant growth and the response to abiotic stress. Santalum album L., a semi-parasitic evergreen tree, is renowned for its economically valuable essential oils, positioning it among the most prized tree species. In this study, a meticulous identification and comprehensive analysis of 43 SAUR genes was conducted within S. album. Based on phylogenetic relationships, the SaSAUR genes were systematically categorized into five groups. A collinearity analysis revealed intriguing insights, disclosing 14 segmental duplications and 9 tandem duplications within the SaSAUR genes, emphasizing the pivotal role of duplication in the expansion of this gene family. Noteworthy variations in the expression levels of SaSAUR genes were observed by delving into the SaSAUR transcriptome data from various tissues, including leaves, roots, and heartwood, as well as under salt-stress conditions. Notably, SaSAUR08 and SaSAUR13 were significantly upregulated in heartwood compared with roots and leaves, while SaSAUR18 was markedly more expressed in roots compared with heartwood and leaves. Furthermore, SaSAUR27 and SaSAUR28 were found to respond closely to salt stress, hinting at their potential involvement in the salt-stress response mechanism. This research offers a comprehensive investigation of SAUR genes in S. album and establishes a foundation for future exploration of the SAUR gene family, particularly its relation to growth and salt-stress responses.

GRAS gene family in rye (Secale cereale L.): genome-wide identification, phylogeny, evolutionary expansion and expression analyses

BACKGROUND

The GRAS transcription factor family plays a crucial role in various biological processes in different plants, such as tissue development, fruit maturation, and environmental stress. However, the GRAS family in rye has not been systematically analyzed yet.

RESULTS

In this study, 67 GRAS genes in S. cereale were identified and named based on the chromosomal location. The gene structures, conserved motifs, cis-acting elements, gene replications, and expression patterns were further analyzed. These 67 ScGRAS members are divided into 13 subfamilies. All members include the LHR I, VHIID, LHR II, PFYRE, and SAW domains, and some nonpolar hydrophobic amino acid residues may undergo cross-substitution in the VHIID region. Interested, tandem duplications may have a more important contribution, which distinguishes them from other monocotyledonous plants. To further investigate the evolutionary relationship of the GRAS family, we constructed six comparative genomic maps of homologous genes between rye and different representative monocotyledonous and dicotyledonous plants. The response characteristics of 19 ScGRAS members from different subfamilies to different tissues, grains at filling stages, and different abiotic stresses of rye were systematically analyzed. Paclobutrazol, a triazole-based plant growth regulator, controls plant tissue and grain development by inhibiting gibberellic acid (GA) biosynthesis through the regulation of DELLA proteins. Exogenous spraying of paclobutrazol significantly reduced the plant height but was beneficial for increasing the weight of 1000 grains of rye. Treatment with paclobutrazol, significantly reduced gibberellin levels in grain in the filling period, caused significant alteration in the expression of the DELLA subfamily gene members. Furthermore, our findings with respect to genes, ScGRAS46 and ScGRAS60, suggest that these two family members could be further used for functional characterization studies in basic research and in breeding programmes for crop improvement.

CONCLUSIONS

We identified 67 ScGRAS genes in rye and further analysed the evolution and expression patterns of the encoded proteins. This study will be helpful for further analysing the functional characteristics of ScGRAS genes.

Growth and transcriptional response of wheat and rice to the tertiary amine BMVE

Introduction

Seed vigor is largely a product of sound seed development, maturation processes, genetics, and storage conditions. It is a crucial factor impacting plant growth and crop yield and is negatively affected by unfavorable environmental conditions, which can include drought and heat as well as cold wet conditions. The latter leads to slow germination and increased seedling susceptibility to pathogens. Prior research has shown that a class of plant growth regulators called substituted tertiary amines (STAs) can enhance seed germination, seedling growth, and crop productivity. However, inconsistent benefits have limited STA adoption on a commercial scale.

Methods

We developed a novel seed treatment protocol to evaluate the efficacy of 2-(N-methyl benzyl aminoethyl)-3-methyl butanoate (BMVE), which has shown promise as a crop seed treatment in field trials. Transcriptomic analysis of rice seedlings 24 h after BMVE treatment was done to identify the molecular basis for the improved seedling growth. The impact of BMVE on seed development was also evaluated by spraying rice panicles shortly after flower fertilization and subsequently monitoring the impact on seed traits.

Results

BMVE treatment of seeds 24 h after imbibition consistently improved wheat and rice seedling shoot and root growth in lab conditions. Treated wheat seedlings grown to maturity in a greenhouse also resulted in higher biomass than controls, though only under drought conditions. Treated seedlings had increased levels of transcripts involved in reactive oxygen species scavenging and auxin and gibberellic acid signaling. Conversely, several genes associated with increased reactive oxygen species/ROS load, abiotic stress responses, and germination hindering processes were reduced. BMVE spray increased both fresh and mature seed weights relative to the control for plants exposed to 96 h of heat stress. BMVE treatment during seed development also benefited germination and seedling growth in the next generation, under both ambient and heat stress conditions.

Discussion

The optimized experimental conditions we developed provide convincing evidence that BMVE does indeed have efficacy in plant growth enhancement. The results advance our understanding of how STAs work at the molecular level and provide insights for their practical application to improve crop growth.

Transcriptomics reveals a core transcriptional network of K-type cytoplasmic male sterility microspore abortion in wheat (Triticum aestivum L.)

BACKGROUND

Cytoplasmic male sterility (CMS) plays a crucial role in hybrid production. K-type CMS, a cytoplasmic male sterile line of wheat with the cytoplasms of Aegilops kotschyi, is widely used due to its excellent characteristics of agronomic performance, easy maintenance and easy restoration. However, the mechanism of its pollen abortion is not yet clear.

RESULTS

In this study, wheat K-type CMS MS(KOTS)-90-110 (MS line) and it's fertile near-isogenic line MR (KOTS)-90-110 (MR line) were investigated. Cytological analysis indicated that the anthers of MS line microspore nucleus failed to divide normally into two sperm nucleus and lacked starch in mature pollen grains, and the key abortive period was the uninucleate stage to dinuclear stage. Then, we compared the transcriptome of MS line and MR line anthers at these two stages. 11,360 and 5182 differentially expressed genes (DEGs) were identified between the MS and MR lines in the early uninucleate and binucleate stages, respectively. Based on GO enrichment and KEGG pathways analysis, it was evident that significant transcriptomic differences were "plant hormone signal transduction", "MAPK signaling pathway" and "spliceosome". We identified 17 and 10 DEGs associated with the IAA and ABA signal transduction pathways, respectively. DEGs related to IAA signal transduction pathway were downregulated in the early uninucleate stage of MS line. The expression level of DEGs related to ABA pathway was significantly upregulated in MS line at the binucleate stage compared to MR line. The determination of plant hormone content and qRT-PCR further confirmed that hormone imbalance in MS lines. Meanwhile, 1 and 2 DEGs involved in ABA and Ethylene metabolism were also identified in the MAPK cascade pathway, respectively; the significant up regulation of spliceosome related genes in MS line may be another important factor leading to pollen abortion.

CONCLUSIONS

We proposed a transcriptome-mediated pollen abortion network for K-type CMS in wheat. The main idea is hormone imbalance may be the primary factor, MAPK cascade pathway and alternative splicing (AS) may also play important regulatory roles in this process. These findings provided intriguing insights for the molecular mechanism of microspore abortion in K-type CMS, and also give useful clues to identify the crucial genes of CMS in wheat.

Genetic dissection of monosaccharides contents in rice whole grain using genome-wide association study

The simplest form of carbohydrates are monosaccharides which are the building blocks for the synthesis of polymers or complex carbohydrates. Monosaccharide contents of 197 rice accessions were quantified by HPAEC-PAD in rice (Oryza sativa L.) whole grain (RWG). A genome-wide association study (GWAS) was carried out using 33,812 single nucleotide polymorphisms (SNPs) to identify corresponding genomic regions influencing neutral monosaccharides contents. In total, 49 GWAS signals contained in 17 genomic regions (quantitative trait loci [QTLs]) on seven chromosomes of rice were determined to be associated with monosaccharides contents of whole grain. The QTLs were found for fucose (1), mannose (1), xylose (2), arabinose (2), galactose (4), and rhamnose (7) contents, all of which are novel. Based on co-location of annotated rice genes in the vicinity of GWAS signals, the constituents of the whole grain were associated with the following candidate genes: arabinose content with α-N-arabinofuranosidase, pectinesterase inhibitor, and glucosamine-fructose-6-phosphate aminotransferase 1; xylose content with ZOS1-10 (a C2H2 zinc finger transcription factor [TF]); mannose content with aldose 1-epimerase-like protein and a MYB family TF; galactose content with a GT8 family member (galacturonosyltransferase-like 3), a GRAS family TF, and a GH16 family member (xyloglucan endotransglucosylase/hydrolase xyloglucan 23); fucose content with gibberellin 20 oxidase and a lysine-rich arabinogalactan protein 19, and finally rhamnose content with myo-inositol-1-phosphate synthase, UDP-arabinopyranose mutase, and COBRA-like protein precursor. The results of this study should improve our understanding of the genetic basis of the factors that might be involved in the biosynthesis, regulation, and turnover of monosaccharides in RWG, aiming to enhance the nutritional value of rice grain and impact the related industries.

The Small Auxin-Up RNA SAUR10 Is Involved in the Promotion of Seedling Growth in Rice

Small auxin-up-regulated RNAs (SAURs) are genes rapidly activated in response to auxin hormones, significantly affecting plant growth and development. However, there is limited information available about the specific functions of SAURs in rice due to the presence of extensive redundant genes. In this study, we found that OsSAUR10 contains a conserved downstream element in its 3' untranslated region that causes its transcripts to be unstable, ultimately leading to the immediate degradation of the mRNA in rice. In our investigation, we discovered that OsSAUR10 is located in the plasma membrane, and its expression is regulated in a tissue-specific, developmental, and hormone-dependent manner. Additionally, we created ossaur10 mutants using the CRISPR/Cas9 method, which resulted in various developmental defects such as dwarfism, narrow internodes, reduced tillers, and lower yield. Moreover, histological observation comparing wild-type and two ossaur10 mutants revealed that OsSAUR10 was responsible for cell elongation. However, overexpression of OsSAUR10 resulted in similar phenotypes to the wild-type. Our research also indicated that OsSAUR10 plays a role in regulating the expression of two groups of genes involved in auxin biosynthesis (OsYUCCAs) and auxin polar transport (OsPINs) in rice. Thus, our findings suggest that OsSAUR10 acts as a positive plant growth regulator by contributing to auxin biosynthesis and polar transport.

Genome-Wide Identification and Expression Analysis of the Cucumber FKBP Gene Family in Response to Abiotic and Biotic Stresses

The FKBP (FK506-binding protein) gene family is an important member of the PPlase protease family and plays a vital role during the processes of plant growth and development. However, no studies of the FKBP gene family have been reported in cucumber. In this study, 19 FKBP genes were identified in cucumber, which were located on chromosomes 1, 3, 4, 6, and 7. Phylogenetic analysis divided the cucumber FKBP genes into three subgroups. The FKBP genes in the same subgroup exhibited similar structures and conserved motifs. The cis-acting elements analysis revealed that the promoters of cucumber FKBP genes contained hormone-, stress-, and development-related cis-acting elements. Synteny analysis of the FKBP genes among cucumber, Arabidopsis, and rice showed that 12 kinds of syntenic relationships were detected between cucumber and Arabidopsis FKBP genes, and 3 kinds of syntenic relationships were observed between cucumber and rice FKBP genes. The tissue-specific expression analysis showed that some FKBP genes were expressed in all tissues, while others were only highly expressed in part of the 10 types of tissues. The expression profile analysis of cucumber FKBP genes under 13 types of stresses showed that the CsaV3_1G007080 gene was differentially expressed under abiotic stresses (high temperature, NaCl, silicon, and photoperiod) and biotic stresses (downy mildew, green mottle mosaic virus, Fusarium wilt, phytophthora capsica, angular leaf spot, and root-knot nematode), which indicated that the CsaV3_1G007080 gene plays an important role in the growth and development of cucumber. The interaction protein analysis showed that most of the proteins in the FKBP gene family interacted with each other. The results of this study will lay the foundation for further research on the molecular biological functions of the cucumber FKBP gene family.

Identification and expression analysis of the small auxin-up RNA (SAUR) gene family in Lycium ruthenicum

The plant hormone auxin regulates numerous aspects of plant growth and development, and small auxin-up RNA (SAUR) is the largest family of early auxin response genes in higher plants. SAUR has been implicated in the regulation of multiple biological processes. However, no comprehensive analysis of SAUR genes has been reported in Lycium ruthenicum. L. ruthenicum is a thorny shrub with very pronounced salt and drought tolerance, and studies have shown that stem thorns are related to drought tolerance in L. ruthenicum. In this study, the identification, phylogenetic analysis, and conserved motif prediction of SAUR genes were extensively explored. Furthermore, the tissue expression patterns of selected SAUR genes were assayed with quantitative real-time polymerase chain reaction (RT-qPCR). A total of 33 putative LrSAURs were identified and divided into three clusters in a phylogenetic tree of L. ruthenicum. MEME analysis identified 10 motifs in L. ruthenicum, and the results suggested that motif 1 and motif 3 were widely distributed. Analyzing the transcriptome data of stem thorns at four developmental stages indicated that LrSAURs were differentially expressed in L. ruthenicum, and could be divided into six expression patterns. The RT-qPCR analysis of 21 genes showed that LrSAUR2, LrSAUR8, LrSAUR9, LrSAUR11, LrSAUR12, and LrSAUR19 were mainly expressed in stems and stem thorns, and may be related to stem thorn development.

Overexpression of a novel small auxin-up RNA gene, OsSAUR11, enhances rice deep rootedness

BACKGROUND

Deep rooting is an important factor affecting rice drought resistance. However, few genes have been identified to control this trait in rice. Previously, we identified several candidate genes by QTL mapping of the ratio of deep rooting and gene expression analysis in rice.

RESULTS

In the present work, we cloned one of these candidate genes, OsSAUR11, which encodes a small auxin-up RNA (SAUR) protein. Overexpression of OsSAUR11 significantly enhanced the ratio of deep rooting of transgenic rice, but knockout of this gene did not significantly affect deep rooting. The expression of OsSAUR11 in rice root was induced by auxin and drought, and OsSAUR11-GFP was localized both in the plasma membrane and cell nucleus. Through an electrophoretic mobility shift assay and gene expression analysis in transgenic rice, we found that the transcription factor OsbZIP62 can bind to the promoter of OsSAUR11 and promote its expression. A luciferase complementary test showed that OsSAUR11 interacts with the protein phosphatase OsPP36. Additionally, expression of several auxin synthesis and transport genes (e.g., OsYUC5 and OsPIN2) were down-regulated in OsSAUR11-overexpressing rice plants.

CONCLUSIONS

This study revealed a novel gene OsSAUR11 positively regulates deep rooting in rice, which provides an empirical basis for future improvement of rice root architecture and drought resistance.

Genome-Wide Identification and Functional Analysis of the Roles of SAUR Gene Family Members in the Promotion of Cucumber Root Expansion

Auxin serves as an essential regulator of the expression of many different genes in plants, thereby regulating growth and development. The specific functional roles of members of the SAUR (small auxin-up RNA) auxin early response gene family in the development of cucumber plants, however, remain to be fully clarified. Here, 62 SAUR family genes were identified, followed by their classification into 7 groups that included several functionally associated cis-regulatory elements. Phylogenetic tree and chromosomal location-based analyses revealed a high degree of homology between two cucumber gene clusters and other plants in the Cucurbitaceae family. These findings, together with the results of an RNA-seq analysis, revealed high levels of CsSAUR31 expression within the root and male flower tissues. Plants overexpressing CsSAUR31 exhibited longer roots and hypocotyls. Together, these results provide a basis for further efforts to explore the roles that SAUR genes play in cucumber plants, while also expanding the pool of available genetic resources to guide research focused on plant growth and development.

Small Auxin Up RNA (SAUR) gene family identification and functional genes exploration during the floral organ and fruit developmental stages in pineapple (Ananas comosus L.) and its response to salinity and drought stresses

In plants, sexual reproduction relies on the proper development of floral organs that facilitate the successful development of fruits and seeds. Auxin responsive small auxin-up RNA (SAUR) genes play essential roles in floral organ formation and fruit development. However, little is known about the role of SAUR genes in pineapple floral organ formation and fruit development as well as stress responses. In this study, based on genome information and transcriptome datasets, 52 AcoSAUR genes were identified and grouped into 12 groups. The gene structure analysis revealed that most AcoSAUR genes did not have introns, although auxin-acting elements were abundant in the promoter region of AcoSAUR members. The expression analysis across the multiple flower and fruit development stages revealed differential expression of AcoSAUR genes, indicating a tissue and stage-specific function of AcoSAURs. Correlation analysis and pairwise comparisons between gene expression and tissue specificity identified stamen-, petal-, ovule-, and fruit-specific AcoSAURs involved in pineapple floral organs (AcoSAUR4/5/15/17/19) and fruit development (AcoSAUR6/11/36/50). RT-qPCR analysis revealed that AcoSAUR12/24/50 played positive roles in response to the salinity and drought treatment. This work provides an abundant genomic resource for functional analysis of AcoSAUR genes during the pineapple floral organs and fruit development stages. It also highlights the role of auxin signaling involved in pineapple reproductive organ growth.

Auxin participates in regulating the leaf curl development of Wucai (Brassica campestris L.)

Wucai (Brassica campestris L. ssp. chinensis var. rosularis Tsen) belongs to the Brassica genus of the Cruciferae family, and its leaf curl is a typical feature that distinguishes Wucai from other nonheading cabbage subspecies. Our previous research found that plant hormones were involved in the development of the leaf curl in Wucai. However, the molecular mechanisms and the hormones regulating the formation of leaf curl in Wucai have not yet been reported. This study aimed to understand the molecular functions related to hormone metabolism during the formation of leaf curl in Wucai. A total of 386 differentially expressed genes (DEGs) were identified by transcriptome sequencing of two different morphological parts of the same leaf of Wucai germplasm W7-2, and 50 DEGs were found to be related to plant hormones, which were mainly involved in the auxin signal transduction pathway. Then, we measured the content of endogenous hormones in two different forms of the same leaf of Wucai germplasm W7-2. A total of 17 hormones with differential content were identified, including auxin, cytokinins, jasmonic acids, salicylic acids, and abscisic acid. And we found that treatment with auxin transport inhibitor N-1-naphthylphthalamic acid can affect the leaf curl phenotype of Wucai and pak choi (Brassica rapa L. subsp. Chinensis). These results indicated that plant hormones, especially auxin, are involved in developing the leaf curl of Wucai. Our findings provide a potentially valuable reference for future research on the development of leaf curls.

Insight into the bZIP gene family in Lagenaria siceraria: Genome and transcriptome analysis to understand gene diversification in Cucurbitaceae and the roles of LsbZIP gene expression and function under cold stress

The basic leucine zipper (bZIP) as a well-known transcription factor family, figures prominently in diverse biological and developmental processes and response to abiotic/biotic stresses. However, no knowledge of the bZIP family is available for the important edible Cucurbitaceae crop bottle gourd. Herein, we identified 65 putative LsbZIP genes and characterized their gene structure, phylogenetic and orthologous relationships, gene expression profiles in different tissues and cultivars, and responsive genes under cold stress. The phylogenetic tree of 16 released Cucurbitaceae plant genomes revealed the evolutionary convergence and divergence of bZIP family. Based on the specific domains, LsbZIP family were classified into 12 clades (A-K, S) with similar motifs and exon-intron distribution. 65 LsbZIP genes have undergone 19 segmental and two tandem duplication events with purifying selection. The expression profiling of LsbZIP genes showed tissue-specific but no cultivar-specific pattern. The cold stress-responsive candidate LsbZIP genes were analyzed and validated by RNA-Seq and RT-PCR, providing new insights of transcriptional regulation of bZIP family genes in bottle gourd and their potential functions in cold-tolerant variety breeding.

Genome-wide analysis of the SAUR gene family and function exploration of DlSAUR32 during early longan somatic embryogenesis

The early auxin responsive small auxin up-regulated RNA (SAUR) family is an important gene family in the auxin signal transduction pathway. This study focused on the regulatory mechanism of DlSAUR genes during early somatic embryogenesis (SE) and its response to hormone treatment and abiotic stress. Mining of the available Dimocarpus longan Lour. (D. longan) genome sequence yielded 68 putative SAUR genes. Transcript profiles based on RNA-seq data showed that most of the 24 detected DlSAUR genes were highly expressed in the globular embryos (GE) (10) and most of them responded to heat stress and 2,4-D treatment. The results of qRT-PCR showed that most of DlSAUR genes were up-regulated under auxin inhibitor N-1-naphthylphthalamic acid (NPA) and auxin indole-3-acetic acid (IAA) treatments. Moreover, NPA could promote longan SE. The assay for ATAC-seq data analysis showed that chromatin accessibility of 19 of the 24 DlSAUR genes were open during early SE, and most DlSAUR genes differentially expressed during early SE were not associated with H3K4me1 signal enrichment. The DlSAUR32 was selected for subcellular localization and RNA-seq analysis, which encode a cell nuclear-localized protein. Dual-luciferase assays and transient transformation showed that the transcription factors (TFs) DlWRKY75-1 and DlWRKY75-2 might bind to the DlSAUR32 promoters to inhibition gene transcription. Transient overexpression of DlWRKY75-1 and DlWRKY75-2 decreased IAA content in N. benthamiana leaves. Thus, the regulatory network composed of DlSAUR32 and its related TFs may regulate the early longan SE and be involved in the auxin response regulatory pathway of longan.

Genome-wide identification and expression analysis of the SAUR gene family in foxtail millet (Setaria italica L.)

BACKGROUND

Auxin performs important functions in plant growth and development processes, as well as abiotic stress. Small auxin-up RNA (SAUR) is the largest gene family of auxin-responsive factors. However, the knowledge of the SAUR gene family in foxtail millet is largely obscure.

RESULTS

In the current study, 72 SiSAUR genes were identified and renamed according to their chromosomal distribution in the foxtail millet genome. These SiSAUR genes were unevenly distributed on nine chromosomes and were classified into three groups through phylogenetic tree analysis. Most of the SiSAUR members from the same group showed similar gene structure and motif composition characteristics. Analysis of cis-acting elements showed that many hormone and stress response elements were identified in the promoter region of SiSAURs. Gene replication analysis revealed that many SiSAUR genes were derived from gene duplication events. We also found that the expression of 10 SiSAURs was induced by abiotic stress and exogenous hormones, which indicated that SiSAUR genes may participated in complex physiological processes.

CONCLUSIONS

Overall, these results will be valuable for further studies on the biological role of SAUR genes in foxtail development and response to stress conditions and may shed light on the improvement of the genetic breeding of foxtail millet.

Genome-Wide Identification, Characterization, and Expression Analysis of GRAS Gene Family in Ginger (Zingiber officinale Roscoe)

GRAS family proteins are one of the most abundant transcription factors in plants; they play crucial roles in plant development, metabolism, and biotic- and abiotic-stress responses. The GRAS family has been identified and functionally characterized in some plant species. However, this family in ginger (Zingiber officinale Roscoe), a medicinal crop and non-prescription drug, remains unknown to date. In the present study, 66 GRAS genes were identified by searching the complete genome sequence of ginger. The GRAS family is divided into nine subfamilies based on the phylogenetic analyses. The GRAS genes are distributed unevenly across 11 chromosomes. By analyzing the gene structure and motif distribution of GRAS members in ginger, we found that the GRAS genes have more than one cis-acting element. Chromosomal location and duplication analysis indicated that whole-genome duplication, tandem duplication, and segmental duplication may be responsible for the expansion of the GRAS family in ginger. The expression levels of GRAS family genes are different in ginger roots and stems, indicating that these genes may have an impact on ginger development. In addition, the GRAS genes in ginger showed extensive expression patterns under different abiotic stresses, suggesting that they may play important roles in the stress response. Our study provides a comprehensive analysis of GRAS members in ginger for the first time, which will help to better explore the function of GRAS genes in the regulation of tissue development and response to stress in ginger.

Genome-Wide Analysis of SAUR Gene Family Identifies a Candidate Associated with Fruit Size in Loquat (Eriobotrya japonica Lindl.)

Fruit size is an important fruit quality trait that influences the production and commodity values of loquats (Eriobotrya japonica Lindl.). The Small Auxin Upregulated RNA (SAUR) gene family has proven to play a vital role in the fruit development of many plant species. However, it has not been comprehensively studied in a genome-wide manner in loquats, and its role in regulating fruit size remains unknown. In this study, we identified 95 EjSAUR genes in the loquat genome. Tandem duplication and segmental duplication contributed to the expansion of this gene family in loquats. Phylogenetic analysis grouped the SAURs from Arabidopsis, rice, and loquat into nine clusters. By analyzing the transcriptome profiles in different tissues and at different fruit developmental stages and comparing two sister lines with contrasting fruit sizes, as well as by functional predictions, a candidate gene (EjSAUR22) highly expressed in expanding fruits was selected for further functional investigation. A combination of Indoleacetic acid (IAA) treatment and virus-induced gene silencing revealed that EjSAUR22 was not only responsive to auxin, but also played a role in regulating cell size and fruit expansion. The findings from our study provide a solid foundation for understanding the molecular mechanisms controlling fruit size in loquats, and also provide potential targets for manipulation of fruit size to accelerate loquat breeding.

Tissue-specific transcriptome responses to Fusarium head blight and Fusarium root rot

Fusarium head blight (FHB) and Fusarium root rot (FRR) are important diseases of small-grain cereals caused by Fusarium species. While host response to FHB has been subject to extensive study, very little is known about response to FRR and the transcriptome responses of FHB and FRR have not been thoroughly compared. Brachypodium distachyon (Bd) is an effective model for investigating host responses to both FHB and FRR. In this study the transcriptome response of Bd to F. graminearum (Fg) infection of heads and roots was investigated. An RNA-seq analysis was performed on both Bd FHB and FRR during the early infection. Additionally, an RNA-seq analysis was performed on in vitro samples of Fg for comparison with Fg gene expression in planta. Differential gene expression and gene-list enrichment analyses were used to compare FHB and FRR transcriptome responses in both Bd and Fg. Differential expression of selected genes was confirmed using RT-qPCR. Most genes associated with receptor signalling, cell-wall modification, oxidative stress metabolism, and cytokinin and auxin biosynthesis and signalling genes were generally upregulated in FHB or were downregulated in FRR. In contrast, Bd genes involved in jasmonic acid and ethylene biosynthesis and signalling, and antimicrobial production were similarly differentially expressed in both tissues in response to infection. A transcriptome analysis of predicted Fg effectors with the same infected material revealed elevated expression of core tissue-independent genes including cell-wall degradation enzymes and the gene cluster for DON production but also several tissue-dependent genes including those for aurofusarin production and cutin degradation. This evidence suggests that Fg modulates its transcriptome to different tissues of the same host.

SAUR63 stimulates cell growth at the plasma membrane

In plants, regulated cell expansion determines organ size and shape. Several members of the family of redundantly acting Small Auxin Up RNA (SAUR) proteins can stimulate plasma membrane (PM) H⁺-ATPase proton pumping activity by inhibiting PM-associated PP2C.D phosphatases, thereby increasing the PM electrochemical potential, acidifying the apoplast, and stimulating cell expansion. Similarly, Arabidopsis thaliana SAUR63 was able to increase growth of various organs, antagonize PP2C.D5 phosphatase, and increase H⁺-ATPase activity. Using a gain-of-function approach to bypass genetic redundancy, we dissected structural requirements for SAUR63 growth-promoting activity. The divergent N-terminal domain of SAUR63 has a predicted basic amphipathic α-helix and was able to drive partial PM association. Deletion of the N-terminal domain decreased PM association of a SAUR63 fusion protein, as well as decreasing protein level and eliminating growth-promoting activity. Conversely, forced PM association restored ability to promote H⁺-ATPase activity and cell expansion, indicating that SAUR63 is active when PM-associated. Lipid binding assays and perturbations of PM lipid composition indicate that the N-terminal domain can interact with PM anionic lipids. Mutations in the conserved SAUR domain also reduced PM association in root cells. Thus, both the N-terminal domain and the SAUR domain may cooperatively mediate the SAUR63 PM association required to promote growth.

Plant organs reach their final shape and size after substantial cell expansion. Proton pumps at the plasma membrane promote cell expansion by acidifying the cell wall to loosen it, and by increasing electrochemical potential across the plasma membrane for solute uptake that maintains intracellular turgor. Plasma-membrane-associated proteins tightly regulate proton pump activity, in order for organs to grow to an appropriate extent. We have studied requirements for activity of one such regulatory protein in the model plant Arabidopsis called SAUR63. This protein is made rapidly in response to plant growth hormones, and it increases proton pump activity to promote organ growth. These activities depend on its binding to anionic lipids in the plasma membrane, and forced plasma membrane association of SAUR63 can increase growth. Many proteins in the same family are found within Arabidopsis and in all land plants, and likely differ in their affinity for the plasma membrane or in other properties. Further studies of other family members may show how such proteins regulate growth under diverse physiological contexts.

A Long Noncoding RNA Derived from lncRNA-mRNA Networks Modulates Seed Vigor

The discovery of long noncoding RNAs (lncRNAs) has filled a great gap in our understanding of posttranscriptional gene regulation in a variety of biological processes related to plant stress responses. However, systematic analyses of the lncRNAs expressed in rice seeds that germinate under cold stress have been elusive. In this study, we performed strand-specific whole transcriptome sequencing in germinated rice seeds under cold stress and normal temperature. A total of 6258 putative lncRNAs were identified and expressed in a stage-specific manner compared to mRNA. By investigating the targets of differentially expressed (DE) lncRNAs of LT-I (phase I of low temperature)/NT-I (phase I of normal temperature), it was shown that the auxin-activated signaling pathway was significantly enriched, and twenty-three protein-coding genes with most of the members of the SAUR family located in chromosome 9 were identified as the candidate target genes that may interact with five lncRNAs. A seed vigor-related lncRNA, SVR, which interplays with the members of the SAUR gene family in cis was eventually identified. The CRISPR/Cas 9 engineered mutations in SVR cause delay of germination. The findings provided new insights into the connection between lncRNAs and the auxin-activated signaling pathway in the regulation of rice seed vigor.

Genome-Wide Identification of TaSAUR Gene Family Members in Hexaploid Wheat and Functional Characterization of TaSAUR66-5B in Improving Nitrogen Use Efficiency

Excessive input of nitrogen fertilizer not only causes a great waste of resources but brings about a series of ecological and environmental problems. Although Small Auxin Up-regulated RNAs (SAURs) participate in diverse biological processes, the function of SAURs in the nitrogen starvation response has not been well-studied. Here, we identified 308 TaSAURs in wheat and divided them into 10 subfamilies. The promoter regions of most TaSAURs contain hormone responsive elements, and their expression levels change under the treatment of different hormones, such as IAA, MeJA, and ABA. Interestingly, overexpression of one of the TaSAUR family members, a nitrogen starvation responsive gene, TaSAUR66-5B, can promote the growth of Arabidopsis and wheat roots. In addition, overexpression of TaSAUR66-5B in Arabidopsis up-regulates the expression levels of auxin biosynthesis related genes, suggesting that overexpression TaSAUR66-5B may promote root growth by increasing the biosynthesis of auxin. Furthermore, overexpression of TaSAUR66-5B in wheat can increase the biomass and grain yields of transgenic plants, as well as the nitrogen concentration and accumulation of both shoots and grains, especially under low nitrogen conditions. This study provides important genomic information of the TaSAUR gene family and lays a foundation for elucidating the functions of TaSAURs in improving nitrogen utilization efficiency in wheat.

Genome–Wide Identification of the GRAS Family Genes in Melilotus albus and Expression Analysis under Various Tissues and Abiotic Stresses

The GRAS gene family is a plant–specific family of transcription factors, which play an important role in many metabolic pathways, such as plant growth and development and stress response. However, there is no report on the comprehensive study of the GRAS gene family of Melilotus albus. Here, we identified 55 MaGRAS genes, which were classified into 8 subfamilies by phylogenetic analysis, and unevenly distributed on 8 chromosomes. The structural analysis indicated that 87% of MaGRAS genes have no intron, which is highly conservative in different species. MaGRAS proteins of the same subfamily have similar protein motifs, which are the source of functional differences of different genomes. Transcriptome and qRT–PCR data were combined to determine the expression of 12 MaGRAS genes in 6 tissues, including flower, seed, leaf, stem, root and nodule, which indicated the possible roles in plant growth and development. Five and seven MaGRAS genes were upregulated under ABA, drought, and salt stress treatments in the roots and shoots, respectively, indicating that they play vital roles in the response to ABA and abiotic stresses in M. albus. Furthermore, in yeast heterologous expression, MaGRAS12, MaGRAS34 and MaGRAS33 can enhance the drought or salt tolerance of yeast cells. Taken together, these results provide basic information for understanding the underlying molecular mechanisms of GRAS proteins and valuable information for further studies on the growth, development and stress responses of GRAS proteins in M. albus.

SAUR8, a small auxin-up RNA Gene in poplar, confers drought tolerance to transgenic Arabidopsis plants

SAUR (small auxin-up RNA) is an early auxin-responsive gene. In this study, a novel SAUR gene PtSAUR8 was cloned from poplar (Populus trichocarpa), and subcellular location analysis showed that it is targeted to the nuclear membrane. In addition, PtSAUR8 overexpression in Arabidopsis improved the plant resistance to drought stress. Meanwhile, the loss of function mutant saur53 showed more drought sensitivity compared to the WT. PtSAUR8 conferred drought tolerance in transgenic Arabidopsis, as determined through phenotypic and stress-associated physiological indicator analyses, namely, root length, germination rate, relative water content, proline content, CAT content, POD content, malondialdehyde content, hydrogen peroxide content, and relative conductivity. In addition, after the 1 μM abscisic acid (ABA) treatment, the PtSAUR8-OE lines promoted stomata closure. Quantitative fluorescence analysis of related genes induced by drought mutant stress further confirmed that overexpression of PtSAUR8 can improve drought resistance in transgenic Arabidopsis lines. Therefore, PtSAUR8 may play a role in plant drought resistance through ABA-mediated pathways; thus, providing new research materials for molecular breeding of poplar resistance.

Comparative Genomic Analysis of SAUR Gene Family, Cloning and Functional Characterization of Two Genes (PbrSAUR13 and PbrSAUR52) in Pyrus bretschneideri

The SAUR (small auxin-up RNA) gene family is the biggest family of early auxin response genes in higher plants and has been associated with the control of a variety of biological processes. Although SAUR genes had been identified in several genomes, no systematic analysis of the SAUR gene family has been reported in Chinese white pear. In this study, comparative and systematic genomic analysis has been performed in the SAUR gene family and identified a total of 116 genes from the Chinese white pear. A phylogeny analysis revealed that the SAUR family could be classified into four groups. Further analysis of gene structure (introns/exons) and conserved motifs showed that they are diverse functions and SAUR-specific domains. The most frequent mechanisms are whole-genome duplication (WGD) and dispersed duplication (DSD), both of which may be important in the growth of the SAUR gene family in Chinese white pear. Moreover, cis-acting elements of the PbrSAUR genes were found in promoter regions associated with the auxin-responsive elements that existed in most of the upstream sequences. Remarkably, the qRT-PCR and transcriptomic data indicated that PbrSAUR13 and PbrSAUR52 were significantly expressed in fruit ripening. Subsequently, subcellular localization experiments revealed that PbrSAUR13 and PbrSAUR52 were localized in the nucleus. Moreover, PbrSAUR13 and PbrSAUR52 were screened for functional verification, and Dangshan pear and frandi strawberry were transiently transformed. Finally, the effects of these two genes on stone cells and lignin were analyzed by phloroglucinol staining, Fourier infrared spectroscopy, and qRT-PCR. It was found that PbrSAUR13 promoted the synthesis and accumulation of stone cells and lignin, PbrSAUR52 inhibited the synthesis and accumulation of stone cells and lignin. In conclusion, these results indicate that PbrSAUR13 and PbrSAUR52 are predominantly responsible for lignin inhibit synthesis, which provides a basic mechanism for further study of PbrSAUR gene functions.

Genome-wide characterization and expression analysis of SAUR gene family in Melon (Cucumis melo L.)

We identified 66 melon SAUR genes by bioinformatic analyses. CmSAUR19, 38, 58, 62 genes are specifically expressed in different stages of fruit growth, suggesting their participation in regulating fruit development. Auxin plays a crucial role in plant growth by regulating the multiple auxin response genes. However, in melon (Cucumis melo L.), the functions of the auxin early response gene family SAUR (Small auxin up RNA) genes in fruit development are still poorly understood. Through genome-wide characterization of CmSAUR family in melon, we identified a total of 66 CmSAUR genes. The open reading frames of the CmSAUR genes ranged from 234 to 525 bp, containing only one exon and lacking introns. Chromosomal position and phylogenetic tree analyses found that the two gene clusters in the melon chromosome are highly homologous in the Cucurbitaceae plants. Among the four conserved motifs in CmSAUR proteins, motif 1, motif 2, and motif 3 located in most of the family protein sequences, and motif 4 showed a close correlation with the two gene clusters. The CmSAUR28 and CmSAUR58 genes have auxin response elements located in the promoters, suggesting they may be involved in the auxin signaling pathway to regulate fruit development. Through transcriptomic profiling in the four developmental stages of fruit and different lateral organs, we selected 16 differentially-expressed SAUR genes for performing further expression analyses. qRT-PCR results showed that five SAUR genes are specifically expressed in flower organs and ovaries. CmSAUR19 and CmSAUR58 were significantly accumulated in the early developmental stage of the fruit. CmSAUR38 and CmAUR62 showed high expression in the climacteric and post-climacteric stages, suggesting their specific role in controlling fruit ripening. This work provides a foundation for further exploring the function of the SAUR gene in fruit development.

Genome-wide survey and expression analysis of GRAS transcription factor family in sweetpotato provides insights into their potential roles in stress response

BACKGROUND

The plant-specific GRAS transcription factors play pivotal roles in various adverse environmental conditions. Numerous GRAS genes have been explored and characterized in different plants, however, comprehensive survey on GRASs in sweetpotato is lagging.

RESULTS

In this study, 72 putative sweetpotato IbGRAS genes with uneven distribution were isolated on 15 chromosomes and classified into 12 subfamilies supported by gene structures and motif compositions. Moreover, both tandem duplication and segmental duplication events played critical roles in the expansion of sweetpotato GRAS genes, and the collinearity between IbGRAS genes and the related orthologs from nine other plants further depicted evolutionary insights into GRAS gene family. RNA-seq analysis under salt stress and qRT-PCR detection of 12 selected IbGRAS genes demonstrated their significant and varying inductions under multiple abiotic stresses (salt, drought, heat and cold) and hormone treatments (ABA, ACC and JA). Consistently, the promoter regions of IbGRAS genes harbored a series of stress- and hormone-associated cis-acting elements. Among them, IbGRAS71, the potential candidate for breeding tolerant plants, was characterized as having transactivation activity in yeasts, while IbGRAS-2/-4/-9 did not. Moreover, a complex interaction relationship between IbGRASs was observed through the interaction network analysis and yeast two-hybrid assays.

CONCLUSIONS

Our results laid a foundation for further functional identifications of IbGRAS genes, and multiple members may serve as potential regulators for molecular breeding of tolerant sweetpotato.

CBL-Interacting Protein Kinase OsCIPK18 Regulates the Response of Ammonium Toxicity in Rice Roots

Ammonium ( NH 4 + ) is one of the major nitrogen sources for plants. However, excessive ammonium can cause serious harm to the growth and development of plants, i.e., ammonium toxicity. The primary regulatory mechanisms behind ammonium toxicity are still poorly characterized. In this study, we showed that OsCIPK18, a CBL-interacting protein kinase, plays an important role in response to ammonium toxicity by comparative analysis of the physiological and whole transcriptome of the T-DNA insertion mutant (cipk18) and the wild-type (WT). Root biomass and length of cipk18 are less inhibited by excess NH 4 + compared with WT, indicating increased resistance to ammonium toxicity. Transcriptome analysis reveals that OsCIPK18 affects the NH 4 + uptake by regulating the expression of OsAMT1;2 and other NH 4 + transporters, but does not affect ammonium assimilation. Differentially expressed genes induced by excess NH 4 + in WT and cipk18 were associated with functions, such as ion transport, metabolism, cell wall formation, and phytohormones signaling, suggesting a fundamental role for OsCIPK18 in ammonium toxicity. We further identified a transcriptional regulatory network downstream of OsCIPK18 under NH 4 + stress that is centered on several core transcription factors. Moreover, OsCIPK18 might function as a transmitter in the auxin and abscisic acid (ABA) signaling pathways affected by excess ammonium. These data allowed us to define an OsCIPK18-regulated/dependent transcriptomic network for the response of ammonium toxicity and provide new insights into the mechanisms underlying ammonium toxicity.

PpSAUR43, an Auxin-Responsive Gene, Is Involved in the Post-Ripening and Softening of Peaches

Auxin’s role in the post-ripening of peaches is widely recognized as important. However, little is known about the processes by which auxin regulates fruit post-ripening. As one of the early auxin-responsive genes, it is critical to understand the role of small auxin-up RNA (SAUR) genes in fruit post-ripening and softening. Herein, we identified 72 PpSAUR auxin-responsive factors in the peach genome and divided them into eight subfamilies based on phylogenetic analysis. Subsequently, the members related to peach post-ripening in the PpSAUR gene family were screened, and we targeted PpSAUR43. The expression of PpSAUR43 was decreased with fruit post-ripening in melting flesh (MF) fruit and was high in non-melting flesh (NMF) fruit. The overexpression of PpSAUR43 showed a slower rate of firmness decline, reduced ethylene production, and a delayed fruit post-ripening process. The MADS-box gene family plays an important regulatory role in fruit ripening. In this study, we showed with yeast two-hybrid (Y2H) and bimolecular fluorescence complementation (BIFC) experiments that PpSAUR43 can interact with the MADS-box transcription factor PpCMB1(PpMADS2), which indicates that PpSAUR43 may inhibit fruit ripening by suppressing the function of the PpCMB1 protein. Together, these results indicate that PpSAUR43 acts as a negative regulator involved in the peach post-ripening process.

Transcriptomic identification of long noncoding RNAs and their hormone-associated nearby coding genes involved in the differential development of caryopses localized on different branches in rice

Long noncoding RNAs (lncRNAs) play important regulatory roles in caryopsis development and grain size in rice. However, whether there exist differences in lncRNA expression between caryopses located on primary branches (CPB) and caryopses located on secondary branches (CSB) that contribute to their differential development remains elusive. Here, we performed transcriptome-wide analysis to identify 2,273 lncRNAs expressed in CPB and CSB at 0, 5, 12, and 20 days after flowering (DAF). Although these lncRNAs were widely distributed, the majority were located in intergenic regions of the 12 rice chromosomes. Based on gene expression cluster analysis, lncRNAs expressed in CPB and CSB were clustered into two subtypes in a position-independent manner: one includes 0- and 5-DAF CPB and CSB, and 12-DAF CSB; the second includes 12-DAF CPB and 20-DAF CPB and CSB. Furthermore, according to the expression value of each lncRNA, K-means cluster analysis revealed 135 early-stage, 116 middle-stage, and 114 late-stage expression-delayed lncRNAs in CSB. Then, we analyzed the expression values of the expression-delayed lncRNAs and nearby coding genes (100 kb upstream and downstream of the lncRNAs), and found 631 lncRNA-mRNA pairs, including 258 lncRNAs and 571 nearby coding genes, some of which are related to hormone-regulated grain development. These results suggested that expression-delayed lncRNAs in CSB may regulate the development of CPB and CSB, providing insight into the mechanism underlying the developmental differences between CPB and CSB, and the differences in grain yield.

Chinese Cherry (Cerasus pseudocerasus Lindl.) ARF7 Participates in Root Development and Responds to Drought and Low Phosphorus

In this paper, an auxin-responsive transcription factor, CpARF7, was isolated from the roots of Chinese cherry (Cerasus pseudocerasus Lindl. Cv. “Manao Hong”). CpARF7 is highly homologous to AtARF7 or AtARF19 in Arabidopsis, and PavARF1 or PavARF14 in sweet cherry. However, in the phenotype of transgenic tomatoes, the root morphology changed, the main root elongated, and the lateral root increased. Both drought treatment and low-phosphorus conditions can elongate the roots of transgenic tomatoes. In addition, the drought resistance and low-phosphorus tolerance of the transgenic lines are improved, and the POD, SOD, and CAT activities under drought and low-phosphorus environments are increased. There is an effect on the tomato somatotropin suppressor gene, SlIAAs, in which SlIAA1/14/19/29 are up-regulated and SlIAA2/11/12/16 are down-regulated. These results indicate that CpARF7 plays an essential regulatory role in root formation and abiotic stress response, and deepens the understanding of auxin-responsive genes in root growth and abiotic stress.

The SAUR gene family in coffee: genome-wide identification and gene expression analysis during somatic embryogenesis

BACKGROUND

Small auxin-up RNA (SAUR) genes form a wide family supposedly involved in different physiological and developmental processes in plants such as leaf senescence, auxin signaling and transport, hypocotyl development and tolerance to abiotic stresses. The transcription of SAUR genes is quickly induced by auxins, a group of phytohormones of major importance on embryo development. To better understand the distribution and expression profile of such still not explored family in Coffea sp., especially during the development of somatic embryogenesis (SE), SAUR members were characterized in silico using the available Coffea canephora genome data and analyzed for gene expression by RT-qPCR in C. arabica embryogenic samples.

METHODS AND RESULTS

Over C. canephora genome 31 CcSAURs were distributed by 11 chromosomes. Out of these 31 gene members, 5 SAURs were selected for gene expression analysis in C. arabica embryogenic materials. CaSAUR12 and CaSAUR18 were the members highly expressed through almost all plant materials. The other genes had more expression in at least one of the developing embryo stages or plantlets. The CaSAUR12 was the only member to exhibit an increased expression in both non-embryogenic calli and the developing embryo stages.

CONCLUSION

The identification of SAUR family on C. canephora genome followed by the analysis of gene expression profile across coffee somatic embryogenesis process on C. arabica represents a further additional step towards a better comprehension of molecular components acting on SE. Along with new research about this gene family such knowledge may support studies about clonal propagation methods via somatic embryogenesis to help the scientific community towards improvements into coffee crop.

Genetic Divergence of Lineage-Specific Tandemly Duplicated Gene Clusters in Four Diploid Potato Genotypes

Potato (Solanum tuberosum L.) is the most important non-grain food crop. Tandem duplication significantly contributes to genome evolution. The objectives of this study were to (i) identify tandemly duplicated genes and compare their genomic distributions across potato genotypes, (ii) investigate the bias in functional specificities, (iii) explore the relationships among coding sequence, promoter and expression divergences associated with tandemly duplicated genes, (iv) examine the role of tandem duplication in generating and expanding lineage-specific gene families, (v) investigate the evolutionary forces affecting tandemly duplicated genes, and (vi) assess the similarities and differences with respect to above mentioned aspects between cultivated genotypes and their wild-relative. In this study, we used well-annotated and chromosome-scale de novo genome assemblies of multiple potato genotypes. Our results showed that tandemly duplicated genes are abundant and dispersed through the genome. We found that several functional specificities, such as disease resistance, stress-tolerance, and biosynthetic pathways of tandemly duplicated genes were differentially enriched across multiple potato genomes. Our results indicated the existence of a significant correlation among expression, promoter, and protein divergences in tandemly duplicated genes. We found about one fourth of tandemly duplicated gene clusters as lineage-specific among multiple potato genomes, and these tended to localize toward centromeres and revealed distinct selection signatures and expression patterns. Furthermore, our results showed that a majority of duplicated genes were retained through sub-functionalization followed by genetic redundancy, while only a small fraction of duplicated genes was retained though neo-functionalization. The lineage-specific expansion of gene families by tandem duplication coupled with functional bias might have significantly contributed to potato's genotypic diversity, and, thus, to adaption to environmental stimuli.

Overexpression of auxin response gene MdIAA24 enhanced cadmium tolerance in apple (Malus domestica)

Cadmium (Cd), a phytotoxic heavy metal accumulated in plants and fruits, has significant adverse effects on plant growth and development as well as human health. In particular, Cd pollution has become a serious agricultural issue in recent years. Apple is one of the most popular fruits consumed at the global scale. Improving apple Cd resistance via reductions in Cd absorption can benefit apple tree growth and ensure fruit safety. In this study, we determined that, under the 200 μM Cd treatment, 35S::MdIAA24 apple plants exhibited more biomass and less Cd accumulation in the tested tissues compared to wild type (WT). Furthermore, the 35S::MdIAA24 apple plants demonstrated more favorable photosynthesis characteristics, less reactive oxygen species (ROS) and a greater amount of active antioxidant enzymes under the Cd condition than WT. The expression levels of the Cd uptake genes were observed to be lower in the 35S::MdIAA24 apple plants compared with those of the WT under the Cd treatment. The results highlight the ability of the overexpression of MdIAA24 to enhance apple Cd resistance by improving antioxidant capacity and reducing Cd absorption.

Comparative transcriptome analysis of root types in salt tolerant and sensitive rice varieties in response to salinity stress

Salinity tolerance in rice is a very important trait, especially in areas that are affected by soil salinity, such as tsunami-devastated areas and coastal regions in rice-producing countries. The roots are the key organs that first detect and respond to salinity stress; thus, it is important to have an understanding of how roots contribute to salinity tolerance in agricultural crops. After salinity treatment of the salt tolerant (Mulai) and sensitive (IR29) rice varieties, it appeared that among the three types of roots, the L-type lateral roots (LLR) were the most sensitive to salinity stress in Mulai and the most tolerant in IR29. The nodal roots (NR) and the S-type lateral roots (SLR) were all negatively affected by salinity treatment in both rice varieties. In order to elucidate the molecular mechanism of the difference in stress response among rice root types, the RNA-seq transcriptome profiles of NR, LLR, and SLR were analyzed in Mulai and IR29. Between the two rice varieties, more transporters were found to participate in the regulation of salt tolerance in Mulai roots, such as those involved in ion and sugar transport. In IR29, many of the genes detected were associated with transcription regulation, including stress-inducible genes such as NAC, WRKY and MYB. Among the different root types, gene expression in LLR and SLR were significantly regulated in both rice varieties. Taken together, the genes identified in this study may be utilized in the varietal improvement of rice with very specific root traits that can enhance tolerance to salinity stress.

Genome-Wide Identification of Hsp90 Gene Family in Perennial Ryegrass and Expression Analysis under Various Abiotic Stresses

The heat shock protein 90 (Hsp90) is a protein produced in plants in response to stress. This study identified and analyzed Hsp90 gene family members in the perennial ryegrass genome. From the results, eight Hsp90 proteins were obtained and their MW, pI and number of amino acid bases varied. The amino acid bases ranged from 526 to 862. The CDS also ranged from 20 (LpHsp0-4) to 1 (LpHsp90-5). The least number of CDS regions was 1 (LpHsp90-5) with 528 kb amino acids, while the highest was 20 (LpHsp90-4) with 862 kb amino acids, which showed diversity among the protein sequences. The phylogenetic tree revealed that Hsp90 genes in Lolium perenne, Arabidopsis thaliana, Oryza sativa and Brachypodium distachyon could be divided into two groups with five paralogous gene pairs and three orthologous gene pairs. The expression analysis after perennial ryegrass was subjected to heat, salt, chromium (Cr), cadmium (Cd), polyethylene glycol (PEG) and abscisic acid (ABA) revealed that LpHsp90 genes were generally highly expressed under heat stress, but only two LpHsp90 proteins were expressed under Cr stresses. Additionally, the expression of the LpHsp90 proteins differed at each time point in all treatments. This study provides the basis for an understanding of the functions of LpHsp90 proteins in abiotic stress studies and in plant breeding.

Insight Between the Epigenetics and Transcription Responding of Cotton Hypocotyl Cellular Elongation Under Salt-Alkaline Stress

Gossypium barbadense is a cultivated cotton not only known for producing superior fiber but also for its salt and alkaline resistance. Here, we used Whole Genome Bisulfite Sequencing (WGBS) technology to map the cytosine methylation of the whole genome of the G. barbadense hypocotyl at single base resolution. The methylation sequencing results showed that the mapping rates of the three samples were 75.32, 77.54, and 77.94%, respectively. In addition, the Bisulfite Sequence (BS) conversion rate was 99.78%. Approximately 71.03, 53.87, and 6.26% of the cytosine were methylated at CG, CHG, and CHH sequence contexts, respectively. A comprehensive analysis of DNA methylation and transcriptome data showed that the methylation level of the promoter region was a positive correlation in the CHH context. Saline-alkaline stress was related to the methylation changes of many genes, transcription factors (TFs) and transposable elements (TEs), respectively. We explored the regulatory mechanism of DNA methylation in response to salt and alkaline stress during cotton hypocotyl elongation. Our data shed light into the relationship of methylation regulation at the germination stage of G. barbadense hypocotyl cell elongation and salt-alkali treatment. The results of this research help understand the early growth regulation mechanism of G. barbadense in response to abiotic stress.

Genome-wide transcriptome signatures of ant-farmed Squamellaria epiphytes reveal key functions in a unique symbiosis

Farming of fungi by ants, termites, or beetles has led to ecologically successful societies fueled by industrial-scale food production. Another type of obligate insect agriculture in Fiji involves the symbiosis between the ant Philidris nagasau and epiphytes in the genus Squamellaria (Rubiaceae) that the ants fertilize, defend, harvest, and depend on for nesting. All farmed Squamellaria form tubers (domatia) with preformed entrance holes and complex cavity networks occupied by P. nagasau. The inner surface of the domatia consists of smooth-surfaced walls where the ants nest and rear their brood, and warty-surfaced walls where they fertilize their crop by defecation. Here, we use RNA sequencing to identify gene expression patterns associated with the smooth versus warty wall types. Since wall differentiation occurred in the most recent common ancestor of all farmed species of Squamellaria, our study also identifies genetic pathways co-opted following the emergence of agriculture. Warty-surfaced walls show many upregulated genes linked to auxin transport, root development, and nitrogen transport consistent with their root-like function; their defense-related genes are also upregulated, probably to protect these permeable areas from pathogen entry. In smooth-surfaced walls, genes functioning in suberin and wax biosynthesis are upregulated, contributing to the formation of an impermeable ant-nesting area in the domatium. This study throws light on a number of functional characteristics of plant farming by ants and illustrates the power of genomic studies of symbiosis.

The Melastoma dodecandrum genome and the evolution of Myrtales

Melastomataceae have abundant morphological diversity with high economic and ornamental merit in Myrtales. The phylogenetic position of Myrtales is still contested. Here, we report the first chromosome-level genome assembly of Melastoma dodecandrum in Melastomataceae. The assembled genome size was 299.81 Mb with a contig N50 value of 3.00 Mb. Genome evolution analysis indicated that M. dodecandrum, Eucalyptus grandis and Punica granatum were clustered into a clade of Myrtales and formed a sister group with the ancestor of fabids and malvids. We found that M. dodecandrum experienced four whole-genome polyploidization events: the ancient event was shared with most eudicots, one event was shared with Myrtales, and the other two events were unique to M. dodecandrum. Moreover, we identified MADS-box genes and found that the AP1-like genes expanded, and AP3-like genes might have undergone subfunctionalization. We found that the SUAR63-like genes and AG-like genes showed different expression patterns in stamens, which may be associated with heteranthery. In addition, we found that LAZY1-like genes were involved in the negative regulation of stem branching development, which may be related to its creeping features. Our study sheds new light on the evolution of Melastomataceae and Myrtales, which provides a comprehensive genetic resource for future research.

Genome-Wide Identification of GRAS Gene Family and Their Responses to Abiotic Stress in Medicago sativa

Alfalfa (Medicago sativa) is a high-quality legume forage crop worldwide, and alfalfa production is often threatened by abiotic environmental stresses. GRAS proteins are important transcription factors that play a vital role in plant development, as well as in response to environmental stress. In this study, the availability of alfalfa genome "Zhongmu No.1" allowed us to identify 51 GRAS family members, i.e., MsGRAS. MsGRAS proteins could be classified into nine subgroups with distinct conserved domains, and tandem and segmental duplications were observed as an expansion strategy of this gene family. In RNA-Seq analysis, 14 MsGRAS genes were not expressed in the leaf or root, 6 GRAS genes in 3 differentially expressed gene clusters were involved in the salinity stress response in the leaf. Moreover, qRT-PCR results confirmed that MsGRAS51 expression was induced under drought stress and hormone treatments (ABA, GA and IAA) but down-regulated in salinity stress. Collectively, our genome-wide characterization, evolutionary, and expression analysis suggested that the MsGRAS proteins might play crucial roles in response to abiotic stresses and hormonal cues in alfalfa. For the breeding of alfalfa, it provided important information on stress resistance and functional studies on MsGRAS and hormone signaling.

Genome-wide identification, expression analysis, and functional study of the GRAS transcription factor family and its response to abiotic stress in sorghum [Sorghum bicolor (L.) Moench]

Background

GRAS, an important family of transcription factors, have played pivotal roles in regulating numerous intriguing biological processes in plant development and abiotic stress responses. Since the sequencing of the sorghum genome, a plethora of genetic studies were mainly focused on the genomic information. The indepth identification or genome-wide analysis of GRAS family genes, especially in Sorghum bicolor, have rarely been studied.

Results

A total of 81 SbGRAS genes were identified based on the S. bicolor genome. They were named SbGRAS01 to SbGRAS81 and grouped into 13 subfamilies (LISCL, DLT, OS19, SCL4/7, PAT1, SHR, SCL3, HAM-1, SCR, DELLA, HAM-2, LAS and OS4). SbGRAS genes are not evenly distributed on the chromosomes. According to the results of the gene and motif composition, SbGRAS members located in the same group contained analogous intron/exon and motif organizations. We found that the contribution of tandem repeats to the increase in sorghum GRAS members was slightly greater than that of fragment repeats. By quantitative (q) RT-PCR, the expression of 13 SbGRAS members in different plant tissues and in plants exposed to six abiotic stresses at the seedling stage were quantified. We further investigated the relationship between DELLA genes, GAs and grain development in S. bicolor. The paclobutrazol treatment significantly increased grain weight, and affected the expression levels of all DELLA subfamily genes. SbGRAS03 is the most sensitive to paclobutrazol treatment, but also has a high response to abiotic stresses.

Conclusions

Collectively, SbGRAs play an important role in plant development and response to abiotic stress. This systematic analysis lays the foundation for further study of the functional characteristics of GRAS genes of S. bicolor.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07848-z.

Genome-wide identification and expression analysis of Raffinose synthetase family in cotton

Background

The Raffinose synthetase (RAFS) genes superfamily is critical for the synthesis of raffinose, which accumulates in plant leaves under abiotic stress. However, it remains unclear whether RAFS contributes to resistance to abiotic stress in plants, specifically in the Gossypium species.

Results

In this study, we identified 74 RAFS genes from G. hirsutum, G. barbadense, G. arboreum and G. raimondii by using a series of bioinformatic methods. Phylogenetic analysis showed that the RAFS gene family in the four Gossypium species could be divided into four major clades; the relatively uniform distribution of the gene number in each species ranged from 12 to 25 based on species ploidy, most likely resulting from an ancient whole-genome polyploidization. Gene motif analysis showed that the RAFS gene structure was relatively conservative. Promoter analysis for cis-regulatory elements showed that some RAFS genes might be regulated by gibberellins and abscisic acid, which might influence their expression levels. Moreover, we further examined the functions of RAFS under cold, heat, salt and drought stress conditions, based on the expression profile and co-expression network of RAFS genes in Gossypium species. Transcriptome analysis suggested that RAFS genes in clade III are highly expressed in organs such as seed, root, cotyledon, ovule and fiber, and under abiotic stress in particular, indicating the involvement of genes belonging to clade III in resistance to abiotic stress. Gene co-expressed network analysis showed that GhRFS2A-GhRFS6A, GhRFS6D, GhRFS7D and GhRFS8A-GhRFS11A were key genes, with high expression levels under salt, drought, cold and heat stress.

Conclusion

The findings may provide insights into the evolutionary relationships and expression patterns of RAFS genes in Gossypium species and a theoretical basis for the identification of stress resistance materials in cotton.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-021-04276-4.

The GRAS gene family in watermelons: identification, characterization and expression analysis of different tissues and root-knot nematode infestations

The family of GRAS plant-specific transcription factor plays diverse roles in numerous biological processes. Despite the identification and characterization of GRAS genes family in dozens of plant species, until now, GRAS members in watermelon (Citrullus lanatus) have not been investigated comprehensively. In this study, using bioinformatic analysis, we identified 37 GRAS genes in the watermelon genome (ClGRAS). These genes are classified into 10 distinct subfamilies based on previous research, and unevenly distributed on 11 chromosomes. Furthermore, a complete analysis was conducted to characterize conserved motifs and gene structures, which revealed the members within same subfamily that have analogous conserved gene structure and motif composition. Additionally, the expression pattern of ClGRAS genes was characterized in fruit flesh and rind tissues during watermelon fruit development and under red light (RL) as well as root knot nematode infestation. Finally, for verification of the availability of public transcriptome data, we also evaluated the expression levels of randomly selected four ClGRAS genes under RL and nematode infection by using qRT-PCR. The qRT-PCR results indicated that several ClGRAS genes were differentially expressed, implying their vital role in RL induction of watermelon resistance against root-knot nematodes. The results obtained in this study could be useful in improving the quality of watermelon.

Genome-Wide Analysis of Tubulin Gene Family in Cassava and Expression of Family Member FtsZ2-1 during Various Stress

Filamentous temperature-sensitive protein Z (Tubulin/FtsZ) family is a group of conserved GTP-binding (guanine nucleotide-binding) proteins, which are closely related to plant tissue development and organ formation as the major component of the cytoskeleton. According to the published genome sequence information of cassava (Manihot esculenta Crantz), 23 tubulin genes (MeTubulins) were identified, which were divided into four main groups based on their type and phylogenetic characteristics. The same grouping generally has the same or similar motif composition and exon–intron structure. Collinear analysis showed that fragment repetition event is the main factor in amplification of cassava tubulin superfamily gene. The expression profiles of MeTubulin genes in various tissue were analyzed, and it was found that MeTubulins were mainly expressed in leaf, petiole, and stem, while FtsZ2-1 was highly expressed in storage root. The qRT-PCR results of the FtsZ2-1 gene under hormone and abiotic stresses showed that indole-3-acetic acid (IAA) and gibberellin A3 (GA3) stresses could significantly increase the expression of the FtsZ2-1 gene, thereby revealing the potential role of FtsZ2-1 in IAA and GA3 stress-induced responses.

A comprehensive transcriptome analysis of contrasting rice cultivars highlights the role of auxin and ABA responsive genes in heat stress response

Sensing a change in ambient temperature is key to survival among all living organisms. Temperature fluctuations due to climate change are a matter of grave concern since it adversely affects growth and eventually the yield of crop plants, including two of the major cereals, i.e., rice and wheat. Thus, to understand the response of rice seedlings to elevated temperatures, we performed microarray-based transcriptome analysis of two contrasting rice cultivars, Annapurna (heat tolerant) and IR64 (heat susceptible), by subjecting their seedlings to 37 °C and 42 °C, sequentially. The transcriptome analyses revealed a set of uniquely regulated genes and related pathways in red rice cultivar Annapurna, particularly associated with auxin and ABA as a part of heat stress response in rice. The changes in expression of few auxin and ABA associated genes, such as OsIAA13, OsIAA20, ILL8, OsbZIP12, OsPP2C51, OsDi19-1 and OsHOX24, among others, were validated under high-temperature conditions using RT-qPCR. In particular, the expression of auxin-inducible SAUR genes was enhanced considerably at both elevated temperatures. Further, using genes that expressed inversely under heat vs. cold temperature conditions, we built a regulatory network between transcription factors (TF) such as HSFs, NAC, WRKYs, bHLHs or bZIPs and their target gene pairs and determined regulatory coordination in their expression under varying temperature conditions. Our work thus provides useful insights into temperature-responsive genes, particularly under elevated temperature conditions, and could serve as a resource of candidate genes associated with thermotolerance or downstream components of temperature sensors in rice.

RNA-Sequencing Reveals Differentially Expressed Rice Genes Functionally Associated with Defense against BPH and WBPH in RILs Derived from a Cross between RP2068 and TN1

BACKGROUND

Rice is staple food for over two billion people. Planthoppers like BPH and WBPH occur together in most of rice growing regions across Asia and cause extensive yield loss by feeding and transmission of disease-causing viruses. Chemical control of the pest is expensive and ecologically disastrous; breeding resistant varieties is an acceptable option. But most of such efforts are focused on BPH with an assumption that these varieties will also be effective against WBPH. No critical studies are available to understand rice resistance, common or otherwise, against these two planthoppers.

RESULTS

Our studies aimed to understand the defense mechanisms in rice line RP2068 against BPH and WBPH through RNA sequencing analysis of a RIL line TR3RR derived from the cross TN1 (susceptible) and RP2068 (resistant) after infestation with BPH or WBPH. Results revealed higher number of differentially expressed genes (DEGs) in BPH infested plants than in WBPH infested plants when compared with the uninfested plants. These DEGs could be grouped into UPUP, DNDN, UPDN and DNUP groups based on whether the DEGs were up (UP) or down (DN) regulated against BPH and WBPH, respectively. Gene ontology analysis, specially of members of the last two groups, revealed differences in plant response to the two planthoppers. Abundance of miRNAs and detection of their target genes also indicated that separate sets of genes were suppressed or induced against BPH and WBPH. These results were validated through the analysis of expression of 27 genes through semi-quantitative and quantitative real-time RT-PCR using a set of five RILs that were genetically identical but with different reaction against the two planthoppers. Coupled with data obtained through pathway analysis involving these 27 genes, expression studies revealed common and differential response of rice RP2068 against BPH and WBPH. Trehalose biosynthesis, proline transport, methylation were key pathways commonly upregulated; glucosinolate biosynthesis, response to oxidative stress, proteolysis, cytokinesis pathways were commonly down regulated; photosynthesis, regulation of transcription, expression and transport of peptides and defense related pathways were exclusively upregulated against WBPH; MYB transcription factor mediated defense induction was exclusive to BPH.

CONCLUSION

Rice defense against the two sympatric planthoppers: BPH and WBPH has distinct features in RP2068. Hence, a conscious combination of resistance to these two pests is essential for effective field management.