Genome-wide identification, phylogenetic, and expression analysis under abiotic stress conditions of LIM gene family in Medicago sativa L

Global identification of LIM genes in response to different heat stress regimes in Lactuca sativa

BACKGROUND

LIM (Lineage-11 (LIN-11), Insulin-1 (ISL-1), and Mechanotransduction-3 (MEC-3)) genes belong to a family that hold ubiquitous properties contributing to organ, seed, and pollen development as well as developmental and cellular responses to biotic and abiotic stresses. Lettuce (Lactuca sativa) is a highly consumed vegetable crop susceptible heat stress. High temperatures limit lettuce's overall yield, quality and marketability. Lettuce LIM genes have not been identified and their role in response to high temperatures is not known. Aiming to identify potential new targets for thermoresilience, we searched for LIM genes in lettuce and compared them with orthologous of several dicotyledons and monocotyledons plant species.

RESULTS

We identified fourteen lettuce LIM genes distributed into eight different subgroups using a genome-wide analysis strategy. Three belonging to DAR (DA means "large" in Chinese) class I, two DAR class II, one in the WLIM1, two in the WLIM2, one in the PLIM1, two in the PLIM2 class, one ßLIM and two δLIMs. No DAR-like were identified in any of the species analyzed including lettuce. Interestingly, unlike other gene families in lettuce which underwent large genome tandem duplications, LIM genes did not increase in number compared to other plant species. The response to heat stress induced a dynamic transcriptional response on LsLIM genes. All heat stress regimes, including night stress, day stress and day and night stress were largely responsible for changes in LIM transcriptional expression.

CONCLUSIONS

Our global analysis at the genome level provides a detailed identification of LIM genes in lettuce and other dicotyledonous and monocotyledonous plant species. Gene structure, physical and chemical properties as well as chromosomal location and Cis-regulatory element analysis together with our gene expression analysis under different temperature regimes identified LsWLIM1, LsWLIM2b, LsDAR3 and LsDAR5 as candidate genes that could be used by breeding programs aiming to produce lettuce varieties able to withstand high temperatures.

The functions of laccase gene GhLAC15 in fiber colouration and development in brown-colored cotton

The monotonicity of color type in naturally colored cottons (NCCs) has become the main limiting factor to their widespread use, simultaneously coexisting with poor fiber quality. The synchronous improvement of fiber quality and color become more urgent and crucial as the demand for sustainable development increases. The homologous gene of wild cotton Gossypium stocksii LAC15 in G. hirsutum, GhLAC15, was also dominantly expressed in the developing fibers of brown cotton XC20 from 5 DPA (day post anthesis) to 25 DPA, especially at the secondary cell wall thickening stage (20 DPA and 25 DPA). In XC20 plants with downregulated GhLAC15 (GhLAC15i), a remarkable reduction in proanthocyanidins (PAs) and lignin contents was observed. Some of the key genes in the phenylpropane and flavonoid biosynthesis pathway were down-regulated in GhLAC15i plants. Notably, the fiber length of GhLAC15i plants showed an obvious increase and the fiber color was lightened. Moreover, we found that the thickness of cotton fiber cell wall was decreased in GhLAC15i plants and the fiber surface became smoother compared to that of WT. Taken together, this study revealed that GhLAC15 played an important role in PAs and lignin biosynthesis in naturally colored cotton fibers. It might mediate fiber color and fiber quality by catalyzing PAs oxidation and lignin polymerization, ultimately regulating fiber colouration and development.

Identification and expression analysis of XIP gene family members in rice

Xylanase inhibitor proteins (XIP) are widely distributed in the plant kingdom, and also exist in rice. However, a systematic bioinformatics analysis of this gene family in rice (OsXIP) has not been conducted to date. In this study, we identified 32 members of the OsXIP gene family and analyzed their physicochemical properties, chromosomal localization, gene structure, protein structure, expression profiles, and interaction networks. Our results indicated that OsXIP genes exhibit an uneven distribution across eight rice chromosomes. These genes generally feature a low number of introns or are intronless, all family members, except for OsXIP20, contain two highly conserved motifs, namely Motif 8 and Motif 9. In addition, it is worth noting that the promoter regions of OsXIP gene family members feature a widespread presence of abscisic acid response elements (ABRE) and gibberellin response elements (GARE-motif and TATC-box). Quantitative Real-time PCR (qRT-PCR) analysis unveiled that the expression of OsXIP genes exhibited higher levels in leaves and roots, with considerable variation in the expression of each gene in these tissues both prior to and following treatments with abscisic acid (ABA) and gibberellin (GA3). Protein interaction studies and microRNA (miRNA) target prediction showed that OsXIP engages with key elements within the hormone-responsive and drought signaling pathways. The qRT-PCR suggested osa-miR2927 as a potential key regulator in the rice responding to drought stress, functioning as tissue-specific and temporally regulation. This study provides a theoretical foundation for further analysis of the functions within the OsXIP gene family.

An integrated multi-omics approach reveals polymethoxylated flavonoid biosynthesis in Citrus reticulata cv. Chachiensis

Citrus reticulata cv. Chachiensis (CRC) is an important medicinal plant, its dried mature peels named "Guangchenpi", has been used as a traditional Chinese medicine to treat cough, indigestion, and lung diseases for several hundred years. However, the biosynthesis of the crucial natural products polymethoxylated flavonoids (PMFs) in CRC remains unclear. Here, we report a chromosome-scale genome assembly of CRC with the size of 314.96 Mb and a contig N50 of 16.22 Mb. Using multi-omics resources, we discover a putative caffeic acid O-methyltransferase (CcOMT1) that can transfer a methyl group to the 3-hydroxyl of natsudaidain to form 3,5,6,7,8,3',4'-heptamethoxyflavone (HPMF). Based on transient overexpression and virus-induced gene silencing experiments, we propose that CcOMT1 is a candidate enzyme in HPMF biosynthesis. In addition, a potential gene regulatory network associated with PMF biosynthesis is identified. This study provides insights into PMF biosynthesis and may assist future research on mining genes for the biosynthesis of plant-based medicines.

Oat AsDA1-2D enhances heat stress tolerance and negatively regulates seed-storage globulin

The importance of oats has increased because of their high nutritional value and health benefits in the human diet. High-temperature stress during the reproductive growth period has a detrimental effect on grain morphology by changing the structure and concentration of several seed-storage proteins. DA1, a conserved ubiquitin-proteasome pathway component, plays an important role in regulating grain size by controlling cell proliferation in maternal integuments during the grain-filling stage. However, there have been no reports or studies on oat DA1 genes. In this study, we identified three DA1-like genes (AsDA1-2D, AsDA1-5A, and AsDA1-1D) using genome-wide analysis. Among these, AsDA1-2D was found to be responsible for high-temperature stress tolerance via a yeast thermotolerance assay. The physical interaction of AsDA1-2D with oat-storage-globulin (AsGL-4D) and a protease inhibitor (AsPI-4D) was observed using yeast two-hybrid screening. A subcellular localization assay revealed that AsDA1-2D and its interacting proteins are localized in the cytosol and plasma membrane. An in vitro pull-down assay showed that AsDA1-2D forms a complex with both AsPI-4D and AsGL-4D. An in vitro cell-free degradation assay showed that AsGL-4D was degraded by AsDA1-2D under high-temperature conditions and that AsPI-4D inhibited the function of AsDA1-2D. These results suggest that AsDA1-2D acts as a cysteine protease and negatively regulates oat-grain-storage-globulin under heat stress.

Cloning and functional analysis of GhDFR1, a key gene of flavonoid synthesis pathway in naturally colored cotton

BACKGROUND

The naturally colored brown cotton fiber is the most widely used environmentally friendly textile material, which primarily contains proanthocyanidins and their derivatives. Many structural genes in the flavonoid synthesis pathway are known to improve the genetic resources of naturally colored cotton. Among them, DFR is a crucial late enzyme to synthesis both anthocyanins and proanthocyanidins in the plant flavonoid pathway.

METHODS

The protein sequences of GhDFRs were analyzed using bioinformatic tools. The expression levels of GhDFRs in various tissues and organs of upland cotton Zongxu1 (ZX1), were analyzed by quantitative real-time PCR, and the expression pattern of GhDFR1 during fiber development of white cotton and brown cotton was analyzed further. The function of GhDFR1 in NCC ZX1 was preliminarily analyzed by virus induced gene silencing (VIGS) technology.

RESULTS

Bioinformatic analysis revealed that GhDFRs sequences in upland cotton genome were extremely conserved. Furthermore, evolutionary tree analysis revealed that the functions of GhDFR1 and GhDFR2, and GhDFR3 and GhDFR4, presented different and shared some similarities. Our study showed GhDFR1 and GhDFR2 were specifically expressed in fibers, while GhDFR3 and GhDFR4 were specifically expressed in petals. GhDFR1 was exclusively expressed in brown cotton fiber at various stages of development and progressively increased with the growth of fiber, but the trend of expression in white cotton was quite the opposite. We silenced GhDFR1 expression in brown cotton fiber using VIGS technology, and observed the VIGS-interference plants. After reducing the expression level of GhDFR1, the period for significant GhDFR1 expression in the developing fibers changed, reducing the content of anthocyanins, and lightening the color of mature cotton fibers.

CONCLUSION

GhDFR1 was preferentially expressed in brown cotton during fiber development. The timing of GhDFR1 expression for flavonoid synthesis altered, resulting in anthocyanin contents reduced and the fiber color of the GhDFR1i lines lightened. These findings showed the role of GhDFR1 in fiber coloration of NCC and provided a new candidate for NCC genetic improvement.

Genome-wide identification, phylogenetic, and expression analysis under abiotic stress conditions of Whirly (WHY) gene family in Medicago sativa L

The WHY family is a group of plant-specific transcription factors, that can bind to single-stranded DNA molecules and play a variety of functions in plant nuclei and organelles, participating in the regulation of plant leaf senescence. It has been identified and analyzed in many species, however, the systematic identification and analysis of the WHY genes family have not yet been reported in alfalfa (Medicago sativa L.). Therefore, to explore the function of alfalfa the WHY genes, and 10 MsWHY genes were identified and further characterized their evolutionary relationship and expression patterns by analyzing the recently published genome of alfalfa. Comprehensive analysis of the chromosome location, physicochemical properties of the protein, evolutionary relationship, conserved motifs, and responses to abiotic stresses of the WHY gene family in alfalfa using bioinformatics methods. The results showed that 10 MsWHY genes were distributed on 10 chromosomes, and collinearity analysis showed that many MsWHYs might be derived from segmental duplications, and these genes are under purifying selection. Based on phylogenetic analyses, the WHY gene family of alfalfa can be divided into four subfamilies: I-IV subfamily, and approximately all the WHY genes within the same subfamily share similar gene structures. The 10 MsWHY gene family members contained 10 motifs, of which motif 2 and motif 4 are the conserved motifs shared by these genes. Furthermore, the analysis of cis-regulatory elements indicated that regulatory elements related to transcription, cell cycle, development, hormone, and stress response are abundant in the promoter sequence of the MsWHY genes. Real-time quantitative PCR demonstrated that MsWHYs gene expression is induced by drought, salt, and methyl jasmonate. The present study serves as a basic foundation for future functional studies on the alfalfa WHY family.

ggmotif: An R Package for the extraction and visualization of motifs from MEME software

MEME (Multiple Em for Motif Elicitation) is the most commonly used tool to identify motifs within deoxyribonucleic acid (DNA) or protein sequences. However, the results generated by the MEMEare saved using file formats .xml and .txt, which are difficult to read, visualize, or integrate with other widely used phylogenetic tree packages, such as ggtree. To overcome this problem, we developed the ggmotif R package, which provides two easy-to-use functions that can facilitate the extraction and visualization of motifs from the results files generated by the MEME. ggmotif can extract the information of the location of motif(s) on the corresponding sequence(s) from the .xml format file and visualize it. Additionally, the data extracted by ggmotif can be easily integrated with the phylogenetic data. On the other hand, ggmotif can obtain the sequence of each motif from the .txt format file and draw the sequence logo with the function ggseqlogo from the ggseqlogo R package. The ggmotif R package is freely available (including examples and vignettes) from GitHub at https://github.com/lixiang117423/ggmotif or from CRAN at https://CRAN.R-project.org/package=ggmotif.

Genome-Wide Identification and Expression Analysis of HSF Transcription Factors in Alfalfa (Medicago sativa) under Abiotic Stress

Alfalfa (Medicago sativa) is one of the most important legume forage species in the world. It is often affected by several abiotic stressors that result in reduced yields and poor growth. Therefore, it is crucial to study the resistance of M. sativa to abiotic stresses. Heat shock transcription factors (HSF) are key players in a number of transcriptional regulatory pathways. These pathways play an essential role in controlling how plants react to different abiotic stressors. Studies on the HSF gene family have been reported in many species but have not yet undergone a thorough analysis in M. sativa. Therefore, in order to identify a more comprehensive set of HSF genes, from the genomic data, we identified 16 members of the MsHSF gene, which were unevenly distributed over six chromosomes. We also looked at their gene architectures and protein motifs, and phylogenetic analysis allowed us to divide them into 3 groups with a total of 15 subgroups. Along with these aspects, we then examined the physicochemical properties, subcellular localization, synteny analysis, GO annotation and enrichment, and protein interaction networks of amino acids. Finally, the analysis of 16 MsHSF genes' expression levels across all tissues and under four abiotic stresses using publicly available RNA-Seq data revealed that these genes had significant tissue-specific expression. Moreover, the expression of most MsHSF genes increased dramatically under abiotic stress, further validating the critical function played by the MsHSF gene family in abiotic stress. These results provided basic information about MsHSF gene family and laid a foundation for further study on the biological role of MsHSF gene in response to stress in M. sativa.

Genome-wide identification of ABA receptor PYL/RCAR gene family and their response to cold stress in Medicago sativa L

Abscisic acid (ABA) is an important phytohormone involved in plant growth, plant development, and the protection of plants against abiotic stresses. PYL/RCAR (pyrabactin resistance/pyr1-like/regulatory components of ABA receptor) is the receptor protein of ABA and the core component of the ABA signal transduction network. The PYL gene family has been identified and analyzed in many species, however, there is no report about the research on the whole genome-wide identification of the alfalfa (Medicago sativa L.) PYL gene family. Therefore, to explore the function of alfalfa PYL genes, 39 MsPYL genes were identified by analyzing the recently published genome of alfalfa. Using bioinformatics methods, we systematically analyzed the chromosome location, protein physicochemical properties, evolutionary relationship, conserved motifs, and response to low-temperature stress of the MsPYL family of alfalfa. The results showed that 39 alfalfa MsPYL genes were distributed on 24 chromosomes, and the analysis of gene duplication events showed that fragment duplication was predominant duplication in alfalfa MsPYL family gene expansion. The phylogenetic tree of MsPYL protein of alfalfa and the phylogenetic tree of PYL genes of 3 species show that the MsPYL gene family can be divided into 3 subfamilies, and the structures of the same subfamilies are relatively similar. The 39 MsPYL gene family members of alfalfa contain 10 Motifs. Motif1, Motif2, Motif3, and Motif5 are the conserved motifs shared by these genes; cis-regulatory elements in promoter regions indicate that regulatory elements related to transcription, cell cycle, development, hormone, and stress response are abundantly present in the MsPYL promoter sequences; Real-time fluorescence quantitative PCR analysis showed that the expression of MsPYL genes can be induced by low-temperature treatment. This study provides a reference for further exploring the structural and functional characterization of the alfalfa PYL gene family.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-01066-3.