Genome-wide characterization and expression analysis of fatty acid desaturase gene family in Camelina sativa

Transcriptome Identification and Analysis of Fatty Acid Desaturase Gene Expression at Different Temperatures in Tausonia pullulans 6A7

Tausonia pullulans 6A7 is a low-temperature yeast strain that can produce lipases. Yeast, which is made up of chassis cells, is an important part of synthetic biology, and the use of the lipase-producing properties of T. pullulans 6A7 for the production of fatty acids provides a new pathway for targeted synthesis in yeast cell factories. In this study, we performed RNA-seq on lipase-producing T. pullulans 6A7 at different temperatures (15 °C, 20 °C, 20 °C without corn oil, and 25 °C). Therefore, a total of 8455 differentially expressed genes were screened, and 16 of them were FAD candidate genes. A Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis of group A (15 °C) vs. group D (25 °C) showed that the pathways of fatty acid biosynthesis (map00061) and the biosynthesis of unsaturated fatty acids (map01040) were significantly enriched. In the proposed temporal analysis of differentially expressed genes among the four temperature modulations, we found differentially expressed genes in nine clusters that had the same expression trends; these genes may be jointly involved in multiple biological processes in T. pullulans 6A7. In addition, we found 16 FAD candidate genes involved in fatty acid biosynthesis, and the expression of these genes had similar expression in the transcriptome trends with the different temperature treatments. These findings will help in future in-depth studies of the function and molecular mechanisms of these important FAD genes involved in fatty acid metabolism in yeast, and they could also be conducive to the establishment of a cellular factory for targeted fatty acid production by using yeast.

Key FAD2, FAD3, and SAD Genes Involved in the Fatty Acid Synthesis in Flax Identified Based on Genomic and Transcriptomic Data

FAD (fatty acid desaturase) and SAD (stearoyl-ACP desaturase) genes play key roles in the synthesis of fatty acids (FA) and determination of oil composition in flax (Linum usitatissimum L.). We searched for FAD and SAD genes in the most widely used flax genome of the variety CDC Bethune and three available long-read assembled flax genomes-YY5, 3896, and Atlant. We identified fifteen FAD2, six FAD3, and four SAD genes. Of all the identified genes, 24 were present in duplicated pairs. In most cases, two genes from a pair differed by a significant number of gene-specific SNPs (single nucleotide polymorphisms) or even InDels (insertions/deletions), except for FAD2a-1 and FAD2a-2, where only seven SNPs distinguished these genes. Errors were detected in the FAD2a-1, FAD2a-2, FAD3c-1, and FAD3d-2 sequences in the CDC Bethune genome assembly but not in the long-read genome assemblies. Expression analysis of the available transcriptomic data for different flax organs/tissues revealed that FAD2a-1, FAD2a-2, FAD3a, FAD3b, SAD3-1, and SAD3-2 were specifically expressed in embryos/seeds/capsules and could play a crucial role in the synthesis of FA in flax seeds. In contrast, FAD2b-1, FAD2b-2, SAD2-1, and SAD2-2 were highly expressed in all analyzed organs/tissues and could be involved in FA synthesis in whole flax plants. FAD2c-2, FAD2d-1, FAD3c-1, FAD3c-2, FAD3d-1, FAD3d-2, SAD3-1, and SAD3-2 showed differential expression under stress conditions-Fusarium oxysporum infection and drought. The obtained results are essential for research on molecular mechanisms of fatty acid synthesis, FAD and SAD editing, and marker-assisted and genomic selection for breeding flax varieties with a determined fatty acid composition of oil.

Genome-Wide Detection of SPX Family and Profiling of CoSPX-MFS3 in Regulating Low-Phosphate Stress in Tea-Oil Camellia

Camellia oleifera a member of the family Theaceae, is a phosphorus (P) tolerator native to southern China. The SPX gene family critically regulates plant growth and development and maintains phosphate (Pi) homeostasis. However, the involvement of SPX genes in Pi signaling in Tea-Oil Camellia remains unknown. In this work, 20 SPX genes were identified and categorized into four subgroups. Conserved domains, motifs, gene structure, chromosomal location and gene duplication events were also investigated in the SPX gene family. Defense and stress responsiveness cis-elements were identified in the SPX gene promoters, which participated in low-Pi stress responses. Based on transcriptome data and qRT-PCR results, nine CoSPX genes had similar expression patterns and eight genes (except CoPHO1H3) were up-regulated at 30 days after exposure to low-Pi stress. CoSPX-MFS3 was selected as a key candidate gene by WGCNA analysis. CoSPX-MFS3 was a tonoplast protein. Overexpression of CoSPX-MFS3 in Arabidopsis promoted the accumulation of total P content and decreased the anthocyanin content. Overexpression of CoSPX-MFS3 could enhance low-Pi tolerance by increased biomass and organic acid contents in transgenic Arabidopsis lines. Furthermore, the expression patterns of seven phosphate starvation genes were higher in transgenic Arabidopsis than those in the wild type. These results highlight novel physiological roles of the SPX family genes in C. oleifera under low-Pi stress, and lays the foundation for a deeper knowledge of the response mechanism of C. oleifera to low-Pi stress.

Comparative transcriptomic and evolutionary analysis of FAD-like genes of Brassica species revealed their role in fatty acid biosynthesis and stress tolerance

BACKGROUND

Fatty acid desaturases (FADs) are involved in regulating plant fatty acid composition by adding double bonds to growing hydrocarbon chain. Apart from regulating fatty acid composition FADs are of great importance, and are involved in stress responsiveness, plant development, and defense mechanisms. FADs have been extensively studied in crop plants, and are broadly classed into soluble and non-soluble fatty acids. However, FADs have not yet been characterized in Brassica carinata and its progenitors.

RESULTS

Here we have performed comparative genome-wide identification of FADs and have identified 131 soluble and 28 non-soluble FADs in allotetraploid B. carinata and its diploid parents. Most soluble FAD proteins are predicted to be resided in endomembrane system, whereas FAB proteins were found to be localized in chloroplast. Phylogenetic analysis classed the soluble and non-soluble FAD proteins into seven and four clusters, respectively. Positive type of selection seemed to be dominant in both FADs suggesting the impact of evolution on these gene families. Upstream regions of both FADs were enriched in stress related cis-regulatory elements and among them ABRE type of elements were in abundance. Comparative transcriptomic data analysis output highlighted that FADs expression reduced gradually in mature seed and embryonic tissues. Moreover, under heat stress during seed and embryo development seven genes remained up-regulated regardless of external stress. Three FADs were only induced under elevated temperature whereas five genes were upregulated under Xanthomonas campestris stress suggesting their involvement in abiotic and biotic stress response.

CONCLUSIONS

The current study provides insights into the evolution of FADs and their role in B. carinata under stress conditions. Moreover, the functional characterization of stress-related genes would exploit their utilization in future breeding programs of B. carinata and its progenitors.

Genome-Wide Characteristics of GH9B Family Members in Melon and Their Expression Profiles under Exogenous Hormone and Far-Red Light Treatment during the Grafting Healing Process

β-1,4-glucanase can not only promote the wound healing of grafted seedlings but can also have a positive effect on a plant's cell wall construction. As a critical gene of β-1,4-glucanase, GH9B is involved in cell wall remodeling and intercellular adhesion and plays a vital role in grafting healing. However, the GH9B family members have not yet been characterized for melons. In this study, 18 CmGH9Bs were identified from the melon genome, and these CmGH9Bs were located on 15 chromosomes. Our phylogenetic analysis of these CmGH9B genes and GH9B genes from other species divided them into three clusters. The gene structure and conserved functional domains of CmGH9Bs in different populations differed significantly. However, CmGH9Bs responded to cis elements such as low temperature, exogenous hormones, drought, and injury induction. The expression profiles of CmGH9Bs were different. During the graft healing process of the melon scion grafted onto the squash rootstock, both exogenous naphthyl acetic acid (NAA) and far-red light treatment significantly induced the upregulated expression of CmGH9B14 related to the graft healing process. The results provided a technical possibility for managing the graft healing of melon grafted onto squash by regulating CmGH9B14 expression.

Comprehensive Analysis of Calcium Sensor Families, CBL and CIPK, in Aeluropus littoralis and Their Expression Profile in Response to Salinity

Plants have acquired sets of highly regulated and complex signaling pathways to respond to unfavorable environmental conditions during evolution. Calcium signaling, as a vital mechanism, enables plants to respond to external stimuli, including abiotic and biotic stresses, and coordinate the basic processes of growth and development. In the present study, two calcium sensor families, CBL and CIPK, were investigated in a halophyte plant, Aeluropus littoralis, with a comprehensive analysis. Here, six AlCBL genes, and twenty AlCIPK genes were studied. The analysis of the gene structure and conserved motifs, as well as physicochemical properties, showed that these genes are highly conserved during evolution. The expression levels of AlCBL genes and AlCIPK genes were evaluated under salt stress in leaf and root tissue. Based on the real-time RT-PCR results, the AlCIPK gene family had a higher variation in mRNA abundance than the AlCBL gene family. AlCIPK genes were found to have a higher abundance in leaves than in roots. The results suggest that the correlation between AlCBL genes and AlCIPK is tissue-specific, and different correlations can be expected in leaves and roots. Based on these correlations, AlCIPK3.1-AlCBL4.1 and AlCIPK1.2-AlCBL4.4 can be co-expressed in the root tissue, while AlCBL10 has the potential to be co-expressed with AlCIPK5, AlCIPK26, and AlCIPK12.3 in the leaf tissue. Our findings reveal valuable information on the structure and function of calcium sensor families in A. littoralis, a halophyte plant, that can be used in future research on the biological function of CBLs and CIPKs on salt stress resistance.

Genome-wide analysis of the laccase (LAC) gene family in Aeluropus littoralis: A focus on identification, evolution and expression patterns in response to abiotic stresses and ABA treatment

Laccases are plant enzymes with essential functions during growth and development. These monophenoloxidases are involved in lignin polymerization, and their expression respond to environmental stress. However, studies of laccases in some plants and fungi have highlighted that many structural and functional aspects of these genes are still unknown. Here, the laccase gene family in Aeluropus littoralis (AlLAC) is described based on sequence structure and expression patterns under abiotic stresses and ABA treatment. Fifteen non-redundant AlLACs were identified from the A. littoralis genome, which showed differences in physicochemical characteristics and gene structure. Based on phylogenetic analysis, AlLACs and their orthologues were classified into five groups. A close evolutionary relationship was observed between LAC gene family members in rice and A. littoralis. According to the interaction network, AlLACs interact more with proteins involved in biological processes such as iron incorporation into the metallo-sulfur cluster, lignin catabolism, regulation of the symbiotic process and plant-type primary cell wall biogenesis. Gene expression analysis of selected AlLACs using real-time RT (reverse transcription)-PCR revealed that AlLACs are induced in response to abiotic stresses such as cold, salt, and osmotic stress, as well as ABA treatment. Moreover, AlLACs showed differential expression patterns in shoot and root tissues. Our findings indicate that AlLACs are preferentially involved in the late response of A. littoralis to abiotic stress.

Genome-Wide Characterization of the Sulfate Transporter Gene Family in Oilseed Crops: Camelina sativa and Brassica napus

Sulfate transporters (SULTRs) are responsible for the uptake of sulfate (SO₄^2-) ions in the rhizosphere by roots and their distribution to plant organs. In this study, SULTR family members in the genomes of two oilseed crops (Camelina sativa and Brassica napus) were identified and characterized based on their sequence structures, duplication events, phylogenetic relationships, phosphorylation sites, and expression levels. In total, 36 and 45 putative SULTR genes were recognized in the genomes of C. sativa and B. napus, respectively. SULTR proteins were predicted to be basophilic proteins with low hydrophilicity in both studied species. According to the observed phylogenetic relationships, we divided the SULTRs into five groups, out of which the SULTR 3 group showed the highest variation. Additionally, several duplication events were observed between the SULTRs. The first duplication event occurred approximately five million years ago between three SULTR 3.1 genes in C. sativa. Furthermore, two subunits were identified in the 3D structures of the SULTRs, which demonstrated that the active binding sites differed between C. sativa and B. napus. According to the available RNA-seq data, the SULTRs showed diverse expression levels in tissues and diverse responses to stimuli. SULTR 3 was expressed in all tissues. SULTR 3.1 was more upregulated in response to abiotic stresses in C. sativa, while SULTR 3.3 and SULTR 2.1 were upregulated in B. napus. Furthermore, SULTR 3 and SULTR 4.1 were upregulated in response to biotic stresses in B. napus. Additionally, the qPCR data showed that the SULTRs in C. sativa were involved in the plant's response to salinity. Based on the distribution of cis-regulatory elements in the promoter region, we speculated that SULTRs might be controlled by phytohormones, such as ABA and MeJA. Therefore, it seems likely that SULTR genes in C. sativa have been more heavily influenced by evolutionary processes and have acquired further diversity. The results reveal new insights of the structures and functions of SULTRs in oilseed crops. However, further analyses, related to functional studies, are needed to uncover the role of SULTRs in the plants' development and growth processes, as well as in their response to stimuli.

Insights into the Sulfate Transporter Gene Family and Its Expression Patterns in Durum Wheat Seedlings under Salinity

Sulfate transporters (SULTRs) are an essential plant transporter class responsible for the absorption and distribution of sulfur, an essential plant growth element. SULTRs are also involved in processes related to growth and development and in response to environmental stimuli. In the present study, 22 TdSULTR family members were identified and characterized in the genome of Triticum turgidum L. ssp. durum (Desf.) using available bioinformatics tools. The expression levels of candidate TdSULTR genes were investigated under salt treatments of 150 and 250 mM NaCl after several different exposure times. TdSULTRs showed diversity in terms of physiochemical properties, gene structure, and pocket sites. TdSULTRs and their orthologues were classified into the known five main plant groups of highly diverse subfamilies. In addition, it was noted that segmental duplication events could lengthen TdSULTR family members under evolutionary processes. Based on pocket site analysis, the amino acids leucine (L), valine (V), and serine (S) were most often detected in TdSULTR protein binding sites. Moreover, it was predicted that TdSULTRs have a high potential to be targeted by phosphorylation modifications. According to promoter site analysis, the plant bioregulators ABA and MeJA were predicted to affect TdSULTR expression patterns. Real-time PCR analysis revealed TdSULTR genes are differentially expressed at 150 mM NaCl but show similar expression in response to 250 mM NaCl. TdSULTR reached a maximum level of expression 72 h after the 250 mM salt treatment. Overall, we conclude that TdSULTR genes are involved in the response to salinity in durum wheat. However, additional studies of functionality are needed to determine their precise function and linked-interaction pathways.

Longer Duration of Active Oil Biosynthesis during Seed Development Is Crucial for High Oil Yield-Lessons from Genome-Wide In Silico Mining and RNA-Seq Validation in Sesame

Sesame, one of the ancient oil crops, is an important oilseed due to its nutritionally rich seeds with high protein content. Genomic scale information for sesame has become available in the public databases in recent years. The genes and their families involved in oil biosynthesis in sesame are less studied than in other oilseed crops. Therefore, we retrieved a total of 69 genes and their translated amino acid sequences, associated with gene families linked to the oil biosynthetic pathway. Genome-wide in silico mining helped identify key regulatory genes for oil biosynthesis, though the findings require functional validation. Comparing sequences of the SiSAD (stearoyl-acyl carrier protein (ACP)-desaturase) coding genes with known SADs helped identify two SiSAD family members that may be palmitoyl-ACP-specific. Based on homology with lysophosphatidic acid acyltransferase (LPAAT) sequences, an uncharacterized gene has been identified as SiLPAAT1. Identified key regulatory genes associated with high oil content were also validated using publicly available transcriptome datasets of genotypes contrasting for oil content at different developmental stages. Our study provides evidence that a longer duration of active oil biosynthesis is crucial for high oil accumulation during seed development. This underscores the importance of early onset of oil biosynthesis in developing seeds. Up-regulating, identified key regulatory genes of oil biosynthesis during early onset of seed development, should help increase oil yields.

Genome-wide identification of XTH genes in Liriodendron chinense and functional characterization of LcXTH21

Liriodendron chinense is a relic tree species of the family Magnoliaceae with multiple uses in timber production, landscape decoration, and afforestation. L. chinense often experiences drought stress in arid areas. However, the molecular basis underlying the drought response of L. chinense remains unclear. Many studies have reported that the xyloglucan endotransglucosylase/hydrolase (XTH) family plays an important role in drought stress resistance. Hereby, to explore the drought resistance mechanism of L. chinense, we identify XTH genes on a genome-wide scale in L. chinense. A total of 27 XTH genes were identified in L. chinense, and these genes were classified into three subfamilies. Drought treatment and RT-qPCR analysis revealed that six LcXTH genes significantly responded to drought stress, especially LcXTH21. Hence, we cloned the LcXTH21 gene and overexpressed it in tobacco via gene transfer to analyze its function. The roots of transgenic plants were more developed than those of wild-type plants under different polyethylene glycol (PEG) concentration, and further RT-qPCR analysis showed that LcXTH21 highly expressed in root compared to aboveground organs, indicating that LcXTH21 may play a role in drought resistance through promoting root development. The results of this study provide new insights into the roles of LcXTH genes in the drought stress response. Our findings will also aid future studies of the molecular mechanisms by which LcXTH genes contribute to the drought response.

Genome-wide identification and expression profiling analysis of sucrose synthase (SUS) and sucrose phosphate synthase (SPS) genes family in Actinidia chinensis and A. eriantha

Sucrose synthase (SUS) is a common sugar-base transfer enzyme in plants, and sucrose phosphate synthase (SPS) is one of the major enzymes in higher plants that regulates sucrose synthesis. However, information of the SPS and SUS gene families in Actinidia, as well as their evolutionary and functional properties, is limited. According to the SPS and SUS proteins conserved domain of Arabidopsis thaliana, we found 6 SPS genes and 6 SUS genes from A. chinensis (cultivar: 'Hongyang'), and 3 SPS genes and 6 SUS genes from A. eriantha (cultivar: 'White'). The novel CDC50 conserved domains were discovered on AcSUS2, and all members of the gene family contain similar distinctive conserved domains. The majority of SUS and SPS proteins were hydrophilic, lipid-soluble enzymes that were expected to be found in the cytoplasm. The tertiary structure of SPS and SUS protein indicated that there were many tertiary structures in SPS, and there were windmill-type and spider-type tertiary structures in SUS. The phylogenetic tree was created using the neighbor-joining method, and members of the SPS and SUS gene families are grouped into three subgroups. Genes with comparable intron counts, conserved motifs, and phosphorylation sites were clustered together first. SPS and SUS were formed through replication among their own family members. AcSPS1, AcSPS2, AcSPS4, AcSPS5, AcSUS5, AcSUS6, AeSPS3, AeSUS3 and AeSUS4 were the important genes in regulating the synthesis and accumulation of sucrose for Actinidia during the fruit growth stages.

In Silico Analysis of Fatty Acid Desaturases Structures in Camelina sativa, and Functional Evaluation of Csafad7 and Csafad8 on Seed Oil Formation and Seed Morphology

Fatty acid desaturases add a second bond into a single bond of carbon atoms in fatty acid chains, resulting in an unsaturated bond between the two carbons. They are classified into soluble and membrane-bound desaturases, according to their structure, subcellular location, and function. The orthologous genes in Camelina sativa were identified and analyzed, and a total of 62 desaturase genes were identified. It was revealed that they had the common fatty acid desaturase domain, which has evolved separately, and the proteins of the same family also originated from the same ancestry. A mix of conserved, gained, or lost intron structure was obvious. Besides, conserved histidine motifs were found in each family, and transmembrane domains were exclusively revealed in the membrane-bound desaturases. The expression profile analysis of C. sativa desaturases revealed an increase in young leaves, seeds, and flowers. C. sativa ω3-fatty acid desaturases CsaFAD7 and CsaDAF8 were cloned and the subcellular localization analysis showed their location in the chloroplast. They were transferred into Arabidopsis thaliana to obtain transgenic lines. It was revealed that the ω3-fatty acid desaturase could increase the C18:3 level at the expense of C18:2, but decreases in oil content and seed weight, and wrinkled phenotypes were observed in transgenic CsaFAD7 lines, while no significant change was observed in transgenic CsaFAD8 lines in comparison to the wild-type. These findings gave insights into the characteristics of desaturase genes, which could provide an excellent basis for further investigation for C. sativa improvement, and overexpression of ω3-fatty acid desaturases in seeds could be useful in genetic engineering strategies, which are aimed at modifying the fatty acid composition of seed oil.

The GASA Gene Family in Cacao (Theobroma cacao, Malvaceae): Genome Wide Identification and Expression Analysis

The gibberellic acid-stimulated Arabidopsis (GASA/GAST) gene family is widely distributed in plants and involved in various physiological and biological processes. These genes also provide resistance to abiotic and biotic stresses, including antimicrobial, antiviral, and antifungal. We are interested in characterizing the GASA gene family and determining its role in various physiological and biological process in Theobroma cacao. Here, we report 17 tcGASA genes distributed on six chromosomes in T. cacao. The gene structure, promoter region, protein structure and biochemical properties, expression, and phylogenetics of all tcGASAs were analyzed. Phylogenetic analyses divided tcGASA proteins into five groups. Among 17 tcGASA genes, nine segmentally duplicating genes were identified which formed four pairs and cluster together in phylogenetic tree. Differential expression analyses revealed that most of the tcGASA genes showed elevated expression in the seeds (cacao food), implying their role in seed development. The differential expression of tcGASAs was recorded between the tolerant and susceptible cultivars of cacao, which indicating their possible role as fungal resistant. Our findings provide new insight into the function, evolution, and regulatory system of the GASA family genes in T.cacao and may suggest new target genes for development of fungi-resistant cacao varieties in breeding programs.

Genome-Wide Analysis of Potassium Channel Genes in Rice: Expression of the OsAKT and OsKAT Genes under Salt Stress

Potassium (K+), as a vital element, is involved in regulating important cellular processes such as enzyme activity, cell turgor, and nutrient movement in plant cells, which affects plant growth and production. Potassium channels are involved in the transport and release of potassium in plant cells. In the current study, three OsKAT genes and two OsAKT genes, along with 11 nonredundant putative potassium channel genes in the rice genome, were characterized based on their physiochemical properties, protein structure, evolution, duplication, in silico gene expression, and protein-protein interactions. In addition, the expression patterns of OsAKTs and OsKATs were studied in root and shoot tissues under salt stress using real-time PCR in three rice cultivars. K+ channel genes were found to have diverse functions and structures, and OsKATs showed high genetic divergence from other K+ channel genes. Furthermore, the Ka/Ks ratios of duplicated gene pairs from the K+ channel gene family in rice suggested that these genes underwent purifying selection. Among the studied K+ channel proteins, OsKAT1 and OsAKT1 were identified as proteins with high potential N-glycosylation and phosphorylation sites, and LEU, VAL, SER, PRO, HIS, GLY, LYS, TYR, CYC, and ARG amino acids were predicted as the binding residues in the ligand-binding sites of K+ channel proteins. Regarding the coexpression network and KEGG ontology results, several metabolic pathways, including sugar metabolism, purine metabolism, carbon metabolism, glycerophospholipid metabolism, monoterpenoid biosynthesis, and folate biosynthesis, were recognized in the coexpression network of K+ channel proteins. Based on the available RNA-seq data, the K+ channel genes showed differential expression levels in rice tissues in response to biotic and abiotic stresses. In addition, the real-time PCR results revealed that OsAKTs and OsKATs are induced by salt stress in root and shoot tissues of rice cultivars, and OsKAT1 was identified as a key gene involved in the rice response to salt stress. In the present study, we found that the repression of OsAKTs, OsKAT2, and OsKAT2 in roots was related to salinity tolerance in rice. Our findings provide valuable insights for further structural and functional assays of K+ channel genes in rice.

New Insights Into Structure and Function of TIFY Genes in Zea mays and Solanum lycopersicum: A Genome-Wide Comprehensive Analysis

The TIFY gene family, a key plant-specific transcription factor (TF) family, is involved in diverse biological processes including plant defense and growth regulation. Despite TIFY proteins being reported in some plant species, a genome-wide comparative and comprehensive analysis of TIFY genes in plant species can reveal more details. In the current study, the members of the TIFY gene family were significantly increased by the identification of 18 and six new members using maize and tomato reference genomes, respectively. Thus, a genome-wide comparative analysis of the TIFY gene family between 48 tomato (Solanum lycopersicum, a dicot plant) genes and 26 maize (Zea mays, a monocot plant) genes was performed in terms of sequence structure, phylogenetics, expression, regulatory systems, and protein interaction. The identified TIFYs were clustered into four subfamilies, namely, TIFY-S, JAZ, ZML, and PPD. The PPD subfamily was only detected in tomato. Within the context of the biological process, TIFY family genes in both studied plant species are predicted to be involved in various important processes, such as reproduction, metabolic processes, responses to stresses, and cell signaling. The Ka/Ks ratios of the duplicated paralogous gene pairs indicate that all of the duplicated pairs in the TIFY gene family of tomato have been influenced by an intense purifying selection, whereas in the maize genome, there are three duplicated blocks containing Ka/Ks > 1, which are implicated in evolution with positive selection. The amino acid residues present in the active site pocket of TIFY proteins partially differ in each subfamily, although the Mg or Ca ions exist heterogeneously in the centers of the active sites of all the predicted TIFY protein models. Based on the expression profiles of TIFY genes in both plant species, JAZ subfamily proteins are more associated with the response to abiotic and biotic stresses than other subfamilies. In conclusion, globally scrutinizing and comparing the maize and tomato TIFY genes showed that TIFY genes play a critical role in cell reproduction, plant growth, and responses to stress conditions, and the conserved regulatory mechanisms may control their expression.

Genome-wide comparative analysis of Mg transporter gene family between Triticum turgidum and Camelina sativa

Magnesium (Mg) as a bimetal plays critical roles in biochemical processes, membrane stability, and enzyme activity. Mg transporters (MGTs) are involving in maintaining Mg homeostasis in cells. Although the MGT family members have been identified in different plant species, there is no comprehensive analysis of the other plants' MGT genes. In the current study, 62 and 41 non-redundant putative MGT proteins were recognized into the genome of Camelina sativa, and Triticum turgidum and they were compared based on physicochemical properties, protein structure, expression, and interaction. All identified MGTs were classified into three subgroups, NIPA, CorA, and MRS2/MGT, based on conserved-motifs distribution. The results showed that the secondary structure pattern in NIPA and MRS2 subfamily members in both studied plant species were highly similar. Furthermore, MGTs encompass the conserved structures and the critical sites mainly in the metal ion and Mg²⁺ binding centers as well as the catalytic sites were observed. The highest numbers of protein channels were predicted in CorA proteins in both C. sativa and T. turgidum with 24 and 17 channel numbers, respectively. The Ser, Pro, Gly, Lys, Tyr, and Arg amino acids were predicted as the binding residues in MGTs channel regions. The expression pattern of identified genes demonstrated that MGT genes have diverse tissue-specific expression and stress response expression patterns. Besides, 147 co-expressed genes with MGTs were clustered into the eight co-expression nodes involved in N-glycan biosynthesis, protein processing in the endoplasmic reticulum, carbon metabolism, biosynthesis of amino acids, and endocytosis. In the present study, all interpretations are based on in silico predictions, which can be used in further studies related to functional genomics of MGT genes.

The GASA Gene Family in Theobroma cacao: Genome wide Identification and Expression Analysis.. [DOI: 10.1101/2021.01.27.425041]

The gibberellic acid-stimulated Arabidopsis (GASA/GAST) gene family is widely distributed in plants. The role of the GASA gene family has been reported previously in various physiological and biological processes, such as cell division, root and seed development, stem growth, and fruit ripening. These genes also provide resistance to abiotic and biotic stresses including antimicrobial, antiviral, and antifungal. Here, we report 17 tcGASA genes in Theobroma cacao L. distributed on six chromosomes. The gene structure, promoter-region sequences, protein structure, and biochemical properties, expression, and phylogenetics of all tcGASAs were analyzed. Phylogenetic analyses divided tcGASA proteins into five groups. The nine segmentally duplicating genes form four pairs and cluster together in phylogenetic tree. Purifying selection pressure was recorded on tcGASA, including duplicated genes. Several stress/hormone-responsive cis-regulatory elements were also recognized in the promoter region of tcGASAs. Differential expression analyses revealed that most of the tcGASA genes showed elevated expression in the seeds (cacao food), implying their role in seed development. The black rod disease of genus Phytophthora caused up to 20–25% loss (700,000 metric tons) in world cacao production. The role of tcGASA genes in conferring fungal resistance was also explored based on RNAseq data against Phytophthora megakarya. The differential expression of tcGASA genes was recorded between the tolerant and susceptible cultivars of cacao plants, which were inoculated with the fungus for 24h and 72h. This differential expression indicating possible role of tcGASA genes to fungal resistant in cacao. Our findings provide new insight into the function, evolution, and regulatory system of the GASA family genes in T. cacao and provide new target genes for development of fungi-resistant cacao varieties in breeding programs.