Heat shock transcription factor (Hsf) gene family in common bean (Phaseolus vulgaris): genome-wide identification, phylogeny, evolutionary expansion and expression analyses at the sprout stage under abiotic stress

Background Common bean (Phaseolus vulgaris) is an essential crop with high economic value. The growth of this plant is sensitive to environmental stress. Heat shock factor (Hsf) is a family of antiretroviral transcription factors that regulate plant defense system against biotic and abiotic stress. To date, few studies have identified and bio-analyzed Hsfs in common bean. Results In this study, 30 Hsf transcription factors (PvHsf1–30) were identified from the PFAM database. The PvHsf1–30 belonged to 14 subfamilies with similar motifs, gene structure and cis-acting elements. The Hsf members in Arabidopsis, rice (Oryza sativa), maize (Zea mays) and common bean were classified into 14 subfamilies. Collinearity analysis showed that PvHsfs played a role in the regulation of responses to abiotic stress. The expression of PvHsfs varied across different tissues. Moreover, quantitative real-time PCR (qRT-PCR) revealed that most PvHsfs were differentially expressed under cold, heat, salt and heavy metal stress, indicating that PvHsfs might play different functions depending on the type of abiotic stress. Conclusions In this study, we identified 30 Hsf transcription factors and determined their location, motifs, gene structure, cis-elements, collinearity and expression patterns. It was found that PvHsfs regulates responses to abiotic stress in common bean. Thus, this study provides a basis for further analysis of the function of PvHsfs in the regulation of abiotic stress in common bean.

Identification and expression analysis of the sucrose synthase gene family in pomegranate (Punica granatum L.)

Background Sucrose synthase (SUS, EC 2.4.1.13) is one of the major enzymes of sucrose metabolism in higher plants. It has been associated with C allocation, biomass accumulation, and sink strength. The SUS gene families have been broadly explored and characterized in a number of plants. The pomegranate (Punica granatum) genome is known, however, it lacks a comprehensive study on its SUS genes family. Methods PgSUS genes were identified from the pomegranate genome using a genome-wide search method. The PgSUS gene family was comprehensively analyzed by physicochemical properties, evolutionary relationship, gene structure, conserved motifs and domains, protein structure, syntenic relationships, and cis-acting elements using bioinformatics methods. The expression pattern of the PgSUS gene in different organs and fruit development stages were assayed with RNA-seq obtained from the NCBI SRA database as well as real-time quantitative polymerase chain reaction (qPCR). Results Five pomegranate SUS genes, located on four different chromosomes, were divided into three subgroupsaccording to the classification of other seven species. The PgSUS family was found to be highly conserved during evolution after studying the gene structure, motifs, and domain analysis. Furthermore, the predicted PgSUS proteins showed similar secondary and tertiary structures. Syntenic analysis demonstrated that four PgSUS genes showed syntenic relationships with four species, with the exception of PgSUS2. Predictive promoter analysis indicated that PgSUS genes may be responsive to light, hormone signaling, and stress stimulation. RNA-seq analysis revealed that PgSUS1/3/4 were highly expressed in sink organs, including the root, flower, and fruit, and particularly in the outer seed coats. qPCR analysis showed also that PgSUS1, PgSUS3, and PgSUS4 were remarkably expressed during fruit seed coat development. Our results provide a systematic overview of the PgSUS gene family in pomegranate, developing the framework for further research and use of functional PgSUS genes.

Genome-Wide Identification and Characterization of the RCI2 Gene Family in Allotetraploid Brassica napus Compared with Its Diploid Progenitors

Brassica napus and its diploid progenitors (B. rapa and B. oleracea) are suitable for studying the problems associated with polyploidization. As an important anti-stress protein, RCI2 proteins widely exist in various tissues of plants, and are crucial to plant growth, development, and stress response. In this study, the RCI2 gene family was comprehensively identified and analyzed, and 9, 9, and 24 RCI2 genes were identified in B. rapa, B. oleracea, and B. napus, respectively. Phylogenetic analysis showed that all of the identified RCI2 genes were divided into two groups, and further divided into three subgroups. Ka/Ks analysis showed that most of the identified RCI2 genes underwent a purifying selection after the duplication events. Moreover, gene structure analysis showed that the structure of RCI2 genes is largely conserved during polyploidization. The promoters of the RCI2 genes in B. napus contained more cis-acting elements, which were mainly involved in plant development and growth, plant hormone response, and stress responses. Thus, B. napus might have potential advantages in some biological aspects. In addition, the changes of RCI2 genes during polyploidization were also discussed from the aspects of gene number, gene structure, gene relative location, and gene expression, which can provide reference for future polyploidization analysis.

Genomic profiling and expression analysis of the diacylglycerol kinase gene family in heterologous hexaploid wheat

The inositol phospholipid signaling system mediates plant growth, development, and responses to adverse conditions. Diacylglycerol kinase (DGK) is one of the key enzymes in the phosphoinositide-cycle (PI-cycle), which catalyzes the phosphorylation of diacylglycerol (DAG) to form phosphatidic acid (PA). To date, comprehensive genomic and functional analyses of DGKs have not been reported in wheat. In this study, 24 DGK gene family members from the wheat genome (TaDGKs) were identified and analyzed. Each putative protein was found to consist of a DGK catalytic domain and an accessory domain. The analyses of phylogenetic and gene structure analyses revealed that each TaDGK gene could be grouped into clusters I, II, or III. In each phylogenetic subgroup, the TaDGKs demonstrated high conservation of functional domains, for example, of gene structure and amino acid sequences. Four coding sequences were then cloned from Chinese Spring wheat. Expression analysis of these four genes revealed that each had a unique spatial and developmental expression pattern, indicating their functional diversification across wheat growth and development processes. Additionally, TaDGKs were also prominently up-regulated under salt and drought stresses, suggesting their possible roles in dealing with adverse environmental conditions. Further cis-regulatory elements analysis elucidated transcriptional regulation and potential biological functions. These results provide valuable information for understanding the putative functions of DGKs in wheat and support deeper functional analysis of this pivotal gene family. The 24 TaDGKs identified and analyzed in this study provide a strong foundation for further exploration of the biological function and regulatory mechanisms of TaDGKs in response to environmental stimuli.

Genome-wide Identification and Expression Analysis of the Cucumber PP2C Genes Family

Background: Type 2C protein phosphatase (PP2Cs) is a negative regulator of ABA signaling pathway, which play important roles in stress signal transduction in plants. However, cucumber (Cucumis sativus L.), as an important economic vegetable, has little research on its PP2C genes family. Results: This study conducted a genome-wide investigation of CsPP2C gene family. Through bioinformatics analysis, 56 CsPP2C genes were identified in cucumber. Based on phylogenetic analysis, the PP2C genes of cucumber and Arabidopsis were divided into 13 groups. Gene structure and conserved motif analysis showed that CsPP2C genes in the same group had similar gene structure and conserved domains. Collinearity analysis showed that segmental duplication events played a key role in the expansion of cucumber PP2C genes family. In addition, the expression of CsPP2Cs under different abiotic treatments was analyzed by qRT-PCR. The results showed that CsPP2C family genes showed different expression patterns under ABA, drought, salt and cold treatment, and a significantly responsive gene CsPP2Cs was obtained (CsPP2C3). By predicting the cis-elements in the promoter, we found that all CsPP2C members contained ABA response elements (ABRE) and drought response elements (MYC). Additionally, the expression patterns of CsPP2C genes were specific in different tissues. Conclusions: The results of this study provide a reference for the genome-wide identification of PP2C gene family in other species, and provide a basis for future studies on the function of PP2C gene in cucumber.

Genome-Wide Analysis of the Cytochrome P450 Gene Family Involved in Salt Tolerance in Gossypium hirsutum

Plant cytochrome P450 (P450) participates in a wide range of biosynthetic reactions and targets a variety of biological molecules. These reactions lead to various fatty acid conjugates, plant hormones, secondary metabolites, lignin, and various defensive compounds. In our previous research, transcriptome analysis was performed on the salt-tolerant upland cotton “Tongyan No. 1.” Many differentially expressed genes (DEGs) belong to the P450 family, and their domains occur widely in plants. In this current research, P450 genes were identified in Gossypium hirsutum with the aid of bioinformatics methods for investigating phylogenetic relations, gene structure, cis-elements, chromosomal localization, and collinearity within a genome. qRT-PCR was conducted to analyze P450 gene expression patterns under salt stress. The molecular weights of the 156 P450 genes were in the range of 5,949.6–245,576.3 Da, and the length of the encoded amino acids for all the identified P450 genes ranged from 51 to 2,144. P450 proteins are divided into four different subfamilies based on phylogenetic relationship, gene structure, and chromosomal localization of gene replication. The length of P450 genes in upland cotton differs greatly, ranging from 1,500 to 13,000 bp. The number of exons in the P450 family genes ranged from 1 to 9, while the number of introns ranged from 0 to 8, and there were similar trends within clusters. A total of 31 cis-acting elements were identified by analyzing 1,500 bp promoter sequences. Differences were found in cis-acting elements among genes. The consistency between qRT-PCR and previous transcriptome analysis of salt tolerance DEGs indicated that they were likely to be involved in the salt tolerance of cotton seedlings. Our results provide valuable information on the evolutionary relationships of genes and functional characteristics of the gene family, which is beneficial for further study of the cotton P450 gene family.

Genome-wide investigation of the AP2/ERF gene family in ginger: evolution and expression profiling during development and abiotic stresses

Background AP2/ERF transcription factors (TFs) constitute one of the largest TF families in plants, which play crucial roles in plant metabolism, growth, and development as well as biotic and abiotic stresses responses. Although the AP2/ERF family has been thoroughly identified in many plant species and several AP2/ERF TFs have been functionally characterized, little is known about this family in ginger (Zingiber officinale Roscoe), an important affinal drug and diet vegetable. Recent completion of the ginger genome sequencing provides an opportunity to investigate the expression profiles of AP2/ERF genes in ginger on a genome-wide basis. Results A total of 163 AP2/ERF genes were obtained in the Z.officinale genome and renamed according to the chromosomal distribution of the ZoAP2/ERF genes. Phylogenetic analysis divided them into three subfamilies, of which 35 belonged to the AP2 subfamily, 120 to ERF, three to RAV, and five to Sololist, respectively, which is in accordance with the number of conserved domains and gene structure analysis. A total of 10 motifs were detected in ZoAP2/ERF genes, and some of the unique motifs were found to be important for the function of ZoAP2/ERF genes. The chromosomal localization, gene structure, and conserved protein motif analyses, as well as the characterization of gene duplication events provided deep insight into the evolutionary features of these ZoAP2/ERF genes. The expression profiles derived from the RNA-seq data and quantitative reserve transcription (qRT-PCR) analysis of ZoAP2/ERFs during development and responses to abiotic stresses were investigated in ginger. Conclusion A comprehensive analysis of the AP2/ERF gene expression patterns in various tissues by RNA-seq and qRT-PCR showed that they played an important role in the growth and development of ginger, and genes that might regulate rhizome and flower development were preliminary identified. In additionally, the ZoAP2/ERF family genes that responded to abiotic stresses were also identified. This study is the first time to identify the ZoAP2/ERF family, which contributes to research on evolutionary characteristics and better understanding the molecular basis for development and abiotic stress response, as well as further functional characterization of ZoAP2/ERF genes with an aim of ginger crop improvement.

Stress Resistance Heat Shock Protein 70 (HSP70) Analysis in Sorghum (Sorghum bicolor L.) at Genome-Wide Level

Heat shock proteins (HSP70) play an important role in many biological processes. However, as typical in Sorghum bicolor, the systematic identification of the HSP70 gene is very limited, and the role of the Hsp70 gene in the evolution of Sorghum bicolor has not been described systematically a lot. To overcome the gap, Insilco analysis of HSP70 gene family was conducted.The investigation was utilizing the bioinformatics method to analyze the HSP70 gene family and it has been identified that 30 HSP70 genes from the genome sequence of Sorghum bicolor. A comprehensive analysis of these 30 identified genes undertaking the analysis of gene structure, phylogeny, and physicochemical properties, subcellular localization, and promoter region analysis. The gene structure visualization analyses revealed that 22 genes contains both 5’ and 3’ UTRS and one 5’ and one 3’ gene and 6 genes without UTR. The highest number of introns was recorded as 12 and those genes have shown that without in any intron. In the promoter region analysis, ten protein motifs are identified and characterized and 2219 cis-acting elements are identified. Among those, the promoter enhancer elements share the highest number (1411) and light-responsive elements share the next value (335). The physicochemical properties analysis revealed that 23 families have an acidic nature while four families are basic and the rests are neutral. In general, the different analyses performed disclosed their structural organization, subcellular localization, physicochemical properties, cis-acting elements, phylogenetic, and understress conditions. This study provides further information for the functional characterization of HSP70 and helps to understand the mechanisms of abiotic stress tolerance under diverse stress conditions in Sorghum bicolor.

Pre-metazoan origin of neuropeptide signalling

Neuropeptides are a diverse class of signalling molecules in metazoans. They occur in all animals with a nervous system and also in neuron-less placozoans. However, their origin has remained unclear because no neuropeptide shows deep homology across lineages and none have been found in sponges. Here, we identify two neuropeptide precursors, phoenixin and nesfatin, with broad evolutionary conservation. By database searches, sequence alignments and gene-structure comparisons we show that both precursors are present in bilaterians, cnidarians, ctenophores and sponges. We also found phoenixin and a secreted nesfatin precursor homolog in the choanoflagellate Salpingoeca rosetta. Phoenixin in particular, is highly conserved, including its cleavage sites, suggesting that prohormone processing occurs also in choanoflagellates. In addition, based on phyletic patterns and negative pharmacological assays we question the originally proposed GPR-173 (SREB3) as a phoenixin receptor. Our findings indicate that signalling by secreted neuropeptide homologs has pre-metazoan origins and thus evolved before neurons.

Whokaryote: distinguishing eukaryotic and prokaryotic contigs in metagenomes based on gene structure

Metagenomics has become a prominent technology to study the functional potential of all organisms in a microbial community. Most studies focus on the bacterial content of these communities, while ignoring eukaryotic microbes. Indeed, many metagenomics analysis pipelines silently assume that all contigs in a metagenome are prokaryotic. However, because of marked differences in gene structure, prokaryotic gene prediction tools fail to accurately predict eukaryotic genes. Here, we developed a classifier that distinguishes eukaryotic from prokaryotic contigs based on foundational differences between these taxa in gene structure. We first developed a random forest classifier that uses intergenic distance, gene density and gene length as the most important features. We show that, with an estimated accuracy of 97%, this classifier with principled features grounded in biology can perform almost as well as the classifiers EukRep and Tiara, which use k-mer frequencies as features. By re-training our classifier with Tiara predictions as additional feature, weaknesses of both types of classifiers are compensated; the result is an enhanced classifier that outperforms all individual classifiers, with an F1-score of 1.00 on precision, recall and accuracy for both eukaryotes and prokaryotes, while still being fast. In a reanalysis of metagenome data from a disease-suppressive plant endosphere microbial community, we show how using Whokaryote to select contigs for eukaryotic gene prediction facilitates the discovery of several biosynthetic gene clusters that were missed in the original study. Our enhanced classifier, which we call ′Whokaryote′, is wrapped in an easily installable package and is freely available from https://git.wageningenur.nl/lotte.pronk/whokaryote.