Origin, evolution, and diversification of inositol 1,4,5-trisphosphate 3-kinases in plants and animals

BACKGROUND

In Eukaryotes, inositol polyphosphates (InsPs) represent a large family of secondary messengers and play crucial roes in various cellular processes. InsPs are synthesized through a series of pohophorylation reactions catalyzed by various InsP kinases in a sequential manner. Inositol 1,4,5-trisphosphate 3-kinase (IP3 3-kinase/IP3K), one member of InsP kinase, plays important regulation roles in InsPs metabolism by specifically phosphorylating inositol 1,4,5-trisphosphate (IP3) to inositol 1,3,4,5-tetrakisphosphate (IP4) in animal cells. IP3Ks were widespread in fungi, plants and animals. However, its evolutionary history and patterns have not been examined systematically.

RESULTS

A total of 104 and 31 IP3K orthologues were identified across 57 plant genomes and 13 animal genomes, respectively. Phylogenetic analyses indicate that IP3K originated in the common ancestor before the divergence of fungi, plants and animals. In most plants and animals, IP3K maintained low-copy numbers suggesting functional conservation during plant and animal evolution. In Brassicaceae and vertebrate, IP3K underwent one and two duplication events, respectively, resulting in multiple gene copies. Whole-genome duplication (WGD) was the main mechanism for IP3K duplications, and the IP3K duplicates have experienced functional divergence. Finally, a hypothetical evolutionary model for the IP3K proteins is proposed based on phylogenetic theory.

CONCLUSION

Our study reveals the evolutionary history of IP3K proteins and guides the future functions of animal, plant, and fungal IP3K proteins.

Functional identification of specialized diterpene synthases from Chamaecyparis obtusa and C. obtusa var. formosana to illustrate the putative evolution of diterpene synthases in Cupressaceae

Chamaecyparis obtusa and C. obtusa var. formosana of the Cupressaceae family are well known for their fragrance and excellent physical properties. To investigate the biosynthesis of unique diterpenoid compounds, diterpene synthase genes for specialized metabolite synthesis were cloned from C. obtusa and C. obtusa var. formosana. Using an Escherichia coli co-expression system, eight diterpene synthases (diTPSs) were characterized. CoCPS and CovfCPS are class II monofunctional (+)-copalyl diphosphate synthases [(+)-CPSs]. Class I monofunctional CoLS and CovfLS convert (+)-copalyl diphosphate [(+)-CPP] to levopimaradiene, CoBRS, CovfBRS1, and CovfBRS3 convert (+)-CPP to (-)-beyerene, and CovfSDS converts (+)-CPP to (-)-sandaracopimaradiene. These enzymes are all monofunctional diterpene syntheses in Cupressaceae family of gymnosperm, and differ from those in Pinaceae. The discovery of the enzyme responsible for the biosynthesis of tetracyclic diterpene (-)-beyerene was characterized for the first time. Diterpene synthases with different catalytic functions exist in closely related species within the Cupressaceae family, indicating that this group of monofunctional diterpene synthases is particularly prone to the evolution of new functions and development of species-specific specialized diterpenoid constituents.

Genome-wide identification and expression analysis of the glycosyl hydrolase family 1 genes in Medicago sativa revealed their potential roles in response to multiple abiotic stresses

Glycoside hydrolase family 1 (GH1) β-glucosidases (BGLUs), are encoded by a large number of genes, which participate in the development and stress response of plants, particularly under biotic and abiotic stresses through the activation of phytohormones. However, there are few studies systematically analyzing stress or hormone-responsive BGLU genes in alfalfa. In this study, a total of 179 BGLU genes of the glycoside hydrolase family 1 were identified in the genome of alfalfa, and then were classified into five distinct clusters. Sequence alignments revealed several conserved and unique motifs among these MsBGLU proteins. Many cis-acting elements related to abiotic stresses and phytohormones were identified in the promoter of some MsBGLUs. Moreover, RNA-seq and RT-qPCR analyses showed that these MsBGLU genes exhibited distinct expression patterns in response to different abiotic stress and hormonal treatments. In summary, this study suggests that MsBGLU genes play crucial roles in response to various abiotic stresses and hormonal responses, and provides candidate genes for stress tolerance breeding in alfalfa.

Genome-wide identification of CaARR-Bs transcription factor gene family in pepper and their expression patterns under salinity stress

In plants, ARRs-B transcription factors play a crucial role in regulating cytokinin signal transduction, abiotic stress resistance, and plant development. A number of adverse environmental conditions have caused severe losses for the pepper (Capsicum annuum L.)-a significant and economically important vegetable. Among the transcription factors of the type B-ARRs family, multiple members have different functions. In pepper, only a few members of the ARRs-B family have been reported and characterized. The current study aimed to characterize ARRs-B transcription factors in C. annuum, including phylogenetic relationships, gene structures, protein motif arrangement, and RT-qPCR expression analyses and their role in salinity stress. In total, ten genes encode CaARRs-B transcription factors (CaARR1 to CaARR10) from the largest subfamily of type-B ARRs were identified in C. annum. The genome-wide analyses of the CaARRs-B family in C. annuum were performed based on the reported ARRs-B genes in Arabidopsis. An analysis of homologous alignments of candidate genes, including their phylogenetic relationships, gene structures, conserved domains, and qPCR expression profiles, was conducted. In comparison with other plant ARRs-B proteins, CaARRs-B proteins showed gene conservation and potentially specialized functions. In addition, tissue-specific expression profiles showed that CaARRs-B genes were differentially expressed, suggesting functionally divergent. CaARRs-B proteins had a typical conserved domain, including AAR-like (pfam: PF00072) and Myb DNA binding (pfam: PF00249) domains. Ten of the CaARRs-B genes were asymmetrically mapped on seven chromosomes in Pepper. Additionally, the phylogenetic tree of CaARRs-B genes from C. annuum and other plant species revealed that CaARRs-B genes were classified into four clusters, which may have evolved conservatively. Further, using quantitative real-time qRT-PCR, the study assessed the expression patterns of CaARRs-B genes in Capsicum annuum seedlings subjected to salt stress. The study used quantitative real-time qRT-PCR to examine CaARRs-B gene expression in Capsicum annuum seedlings under salt stress. Roots exhibited elevated expression of CaARR2 and CaARR9, while leaves showed decreased expression for CaARR3, CaARR4, CaARR7, and CaARR8. Notably, no amplification was observed for CaARR10. This research sheds light on the roles of CaARRs-B genes in pepper's response to salinity stress. These findings enrich our comprehension of the functional implications of CaARRs-B genes in pepper, especially in responding to salinity stress, laying a solid groundwork for subsequent in-depth studies and applications in the growth and development of Capsicum annuum.

Molecular analysis of the 14-3-3 genes in Panax ginseng and their responses to heat stress

Background

Panax Ginseng is a perennial and semi-shady herb with tremendous medicinal value. Due to its unique botanical characteristics, ginseng is vulnerable to various abiotic factors during its growth and development, especially in high temperatures. Proteins encoded by 14-3-3 genes form a highly conserved protein family that widely exists in eukaryotes. The 14-3-3 family regulates the vital movement of cells and plays an essential role in the response of plants to abiotic stresses, including high temperatures. Currently, there is no relevant research on the 14-3-3 genes of ginseng.

Methods

The identification of the ginseng 14-3-3 gene family was mainly based on ginseng genomic data and Hidden Markov Models (HMM). We used bioinformatics-related databases and tools to analyze the gene structure, physicochemical properties, cis-acting elements, gene ontology (GO), phylogenetic tree, interacting proteins, and transcription factor regulatory networks. We analyzed the transcriptome data of different ginseng tissues to clarify the expression pattern of the 14-3-3 gene family in ginseng. The expression level and modes of 14-3-3 genes under heat stress were analyzed by quantitative real-time PCR (qRT-PCR) technology to determine the genes in the 14-3-3 gene family responding to high-temperature stress.

Results

In this study, 42 14-3-3 genes were identified from the ginseng genome and renamed PgGF14-1 to PgGF14-42. Gene structure and evolutionary relationship research divided PgGF14s into epsilon (ε) and non-epsilon (non-ε) groups, mainly located in four evolutionary branches. The gene structure and motif remained highly consistent within a subgroup. The physicochemical properties and structure of the predicted PgGF14 proteins conformed to the essential characteristics of 14-3-3 proteins. RNA-seq results indicated that the detected PgGF14s existed in different organs and tissues but differed in abundance; their expression was higher in roots, stems, leaves, and fruits but lower in seeds. The analysis of GO, cis-acting elements, interacting proteins, and regulatory networks of transcription factors indicated that PgGF14s might participate in physiological processes, such as response to stress, signal transduction, material synthesis-metabolism, and cell development. The qRT-PCR results indicated PgGF14s had multiple expression patterns under high-temperature stress with different change trends in several treatment times, and 38 of them had an apparent response to high-temperature stress. Furthermore, PgGF14-5 was significantly upregulated, and PgGF14-4 was significantly downregulated in all treatment times. This research lays a foundation for further study on the function of 14-3-3 genes and provides theoretical guidance for investigating abiotic stresses in ginseng.

Aldolase of Mytilus galloprovincialis, Lamarck, 1819: Gene structure, tissue specificity of expression level and activity

In the present study, the structure of the fructose-1,6-bisphosphataldolase (FBA) gene in Mytilus galloprovincialis (Lamarck, 1819) was analyzed and its tissue specificity of expression level and activity was determined. A 1092 base pairs (bps) complete coding sequence of the FBA gene was assembled from M. galloprovincialis transcriptome. Only one gene encoding FBA (MgFBA) was identified in the M. galloprovincialis genome. The length of MgFBA was 363 amino acids with a molecular mass of 39.7 kDa. According to the amino acid residues, the detected MgFBA gene is a type I aldolase. The FBA gene in M. galloprovincialis had 7 exons; the maximum intron length was about 2.5 kbps. Intraspecific nucleotide diversity (15 mutations) between MgFBAs from the Mediterranean mussels and the Black Sea mussels (present study) was detected. All mutations were synonymous. Tissue specificity in FBA expression level and activity was established. No direct correlation between these functions was found. The highest level of FBA gene expression is found in muscle tissue. According to the phylogenetic analyses, FBA gene in invertebrates could be considered the ancestral gene of muscle type aldolase, which may explain the character of tissue-specific expression.

Evolution and phylogenetic diversity of the aquaporin gene family in arachnids

Water flux across cells predominantly occurs through the pore formed by the aquaporin channels. Since water balance is one of the most important challenges to terrestrial animals, aquaporin evolution and diversity is known to play roles in animal terrestrialisation. Arachnids (Arthropoda: Chelicerata: Arachnida) are the second most diverse group and represent the pioneer land colonists in animals; however, there remains no thorough investigation on aquaporin evolution and diversity in this evolutionarily important lineage. Here we reported a phylogenetic study of aquaporin evolution and diversity using genomic data from 116 arachnid species covering almost all (15/16) extant orders. A previously unrecognised subfamily related to aquaporin-4 (i.e. Aqp4-like subfamily) via phylogenetic analysis was identified, suggesting certain underestimate of the arachnid aquaporin diversity in earlier studies probably due to limited taxonomic sampling. Further analysis indicates that this subfamily emerged deep within the life tree of arthropods. Gene tree of another Aqp4-like subfamily (PripL) shows an unexpected basal split between acariform mites (Acariformes) and other arachnids. A closer inspection demonstrated that the PripL evolved quickly and has been under differential selection pressure in acariform mites. Evidence is provided that the evolutionarily ancient Glp subfamily (i.e. aquaglyceroporin) is significantly expanded in terrestrial arachnids compared with their marine relatives. Finally, in spite of the phylogenetic diversity, there exists conservation of some exons in size, functional domain, and intron-insertion phase: an 81-bp and a 218-bp exon, respectively, in apq4-like and glp genes across Eumetazoa lineages including arachnids and human beings. Both exons encode the carboxyl-terminal NPA motif, implying the coding and splicing pressure during hundreds of million years of animal evolution. Hypotheses were tested to explore the possible link between these findings and arachnid terrestrialisation.

Genome modeling: From chromatin fibers to genes

The intricacies of the 3D hierarchical organization of the genome have been approached by many creative modeling studies. The specific model/simulation technique combination defines and restricts the system and phenomena that can be investigated. We present the latest modeling developments and studies of the genome, involving models ranging from nucleosome systems and small polynucleosome arrays to chromatin fibers in the kb-range, chromosomes, and whole genomes, while emphasizing gene folding from first principles. Clever combinations allow the exploration of many interesting phenomena involved in gene regulation, such as nucleosome structure and dynamics, nucleosome-nucleosome stacking, polynucleosome array folding, protein regulation of chromatin architecture, mechanisms of gene folding, loop formation, compartmentalization, and structural transitions at the chromosome and genome levels. Gene-level modeling with full details on nucleosome positions, epigenetic factors, and protein binding, in particular, can in principle be scaled up to model chromosomes and cells to study fundamental biological regulation.

Comprehensive genome-wide identification and expression profiling of ADF gene family in Citrus sinensis, induced by endophytic colonization of Beauveria bassiana

Endophytic entomopathogenic species are known to systematically colonize host plants and form symbiotic associations that benefit the plants they live with. The actin-depolymerizing factors (ADFs) are a group of gene family that regulate growth, development, and defense-related functions in plants. Systematic studies of ADF family at the genome-wide level and their expression in response to endophytic colonization are essential to understand its functions but are currently lacking in this field. 14ADF genes were identified and characterized in the Citrus sinensis genome. The ADF genes of C. sinensis were classified into five groups according to the phylogenetic analysis of plant ADFs. Additionally, the cis-acting analysis revealed that these genes play essential role in plant growth/development, phytohormone, and biotic and abiotic responses; and the expression analysis showed that the symbiotic interactions generate a significant expression regulation level of ADF genes in leaves, stems and roots, compared to controls; thus enhancing seedlings' growth. Additionally, the 3D structures of the ADF domain were highly conserved during evolution. These results will be helpful for further functional validation of ADFs candidate genes and provide important insights into the vegetative growth, development and stress tolerance of C. sinensis in responses to endophytic colonization by B. bassiana.

The plant specific SHORT INTERNODES/STYLISH (SHI/STY) proteins: Structure and functions

Plant specific SHORT INTERNODES/STYLISH (SHI/STY) protein is a transcription factor involved in the formation and development of early lateral organs in plants. However, research on the SHI/STY protein family is not focused enough. In this article, we review recent studies on SHI/STY genes and explore the evolution and structure of SHI/STY. The biological functions of SHI/STYs are discussed in detail in this review, and the application of each biological function to modern agriculture is discussed. All SHI/STY proteins contain typical conserved RING-like zinc finger domain and IGGH domain. SHI/STYs are involved in the formation and development of lateral root, stem extension, leaf morphogenesis, and root nodule development. They are also involved in the regulation of pistil and stamen development and flowering time. At the same time, the regulation of some GA, JA, and auxin signals also involves these family proteins. For each aspect, unanswered or poorly understood questions were identified to help define future research areas. This review will provide a basis for further functional study of this gene family.

Genome-wide identification of the rubber tree superoxide dismutase (SOD) gene family and analysis of its expression under abiotic stress

Background

The rubber tree (Hevea brasiliensis) is the only species capable of producing high-quality natural rubber for commercial use, and is often subjected to various abiotic stresses in non-traditional rubber plantation areas. Superoxide dismutase (SOD) is a vital metalloenzyme translated by a SOD gene family member and acts as a first-line of protection in plant cells by catalysing the disproportionation of reactive oxygen species (ROS) to produce H₂O₂ and O₂. However, the SOD gene family is not reported in rubber trees.

Methods

Here, we used hidden markov model (HMM) and BLASTP methods to identify SOD genes in the H. brasiliensis genome. Phylogenetic tree, conserved motifs, gene structures, cis elements, and gene ontology annotation (GO) analyses were performed using MEGA 6.0, MEME, TBtools, PlantCARE, and eggNOG database, respectively. HbSOD gene expression profiles were analysed using quantitative reverse transcription polymerase chain reaction (qRT-PCR).

Results

We identified nine HbSOD genes in the rubber tree genome, including five HbCSDs, two HbFSDs, and two HbMSDs. Phylogenetic relationship analysis classified the SOD proteins from the rubber tree and other related species into three subfamilies. The results of gene structure and conserved motif analysis illustrated that most HbSOD genes have similar exon-intron numbers and conserved motifs in the same evolutionary branch. Five hormone-related, four stress-related, and light-responsive elements were detected in the HbSODs' promoters. HbSODs were expressed in different tissues, gradually increased with leaf development, and were abundantly expressed in mature leaves. HbCSD2 and HbCSD4 was significantly upregulated under low and high temperatures, and salt stress, except for HbCSD2, by heat. Furthermore, most HbSOD genes were significantly upregulated by drought, except HbMSD2. These findings imply that these genes may play vital roles in rubber tree stress resistance. Our results provide a basis for further studies on the functions of HbSOD genes in rubber trees and stress response mechanisms.

False gene and chromosome losses in genome assemblies caused by GC content variation and repeats

BACKGROUND

Many short-read genome assemblies have been found to be incomplete and contain mis-assemblies. The Vertebrate Genomes Project has been producing new reference genome assemblies with an emphasis on being as complete and error-free as possible, which requires utilizing long reads, long-range scaffolding data, new assembly algorithms, and manual curation. A more thorough evaluation of the recent references relative to prior assemblies can provide a detailed overview of the types and magnitude of improvements.

RESULTS

Here we evaluate new vertebrate genome references relative to the previous assemblies for the same species and, in two cases, the same individuals, including a mammal (platypus), two birds (zebra finch, Anna's hummingbird), and a fish (climbing perch). We find that up to 11% of genomic sequence is entirely missing in the previous assemblies. In the Vertebrate Genomes Project zebra finch assembly, we identify eight new GC- and repeat-rich micro-chromosomes with high gene density. The impact of missing sequences is biased towards GC-rich 5'-proximal promoters and 5' exon regions of protein-coding genes and long non-coding RNAs. Between 26 and 60% of genes include structural or sequence errors that could lead to misunderstanding of their function when using the previous genome assemblies.

CONCLUSIONS

Our findings reveal novel regulatory landscapes and protein coding sequences that have been greatly underestimated in previous assemblies and are now present in the Vertebrate Genomes Project reference genomes.

The anaplerotic pyruvate carboxylase from white shrimp Litopenaeus vannamei: Gene structure, molecular characterization, protein modelling and expression during hypoxia

Hypoxic zones are spreading worldwide in marine environments affecting many organisms. Shrimp and other marine crustaceans can withstand environmental hypoxia using several strategies, including the regulation of energy producing metabolic pathways. Pyruvate carboxylase (PC) catalyzes the first reaction of gluconeogenesis to produce oxaloacetate from pyruvate. In mammals, PC also participates in lipogenesis, insulin secretion and other processes, but this enzyme has been scarcely studied in marine invertebrates. In this work, we characterized the gene encoding PC in the white shrimp Litopenaeus vannamei, modelled the protein structure and evaluated its gene expression in hepatopancreas during hypoxia, as well as glucose and lactate concentrations. The PC gene codes for a mitochondrial protein and has 21 coding exons and 4 non-coding exons that generate three transcript variants with differences only in the 5'-UTR. Total PC expression is more abundant in hepatopancreas compared to gills or muscle, indicating tissue-specific expression. Under hypoxic conditions of 1.53 mg/L dissolved oxygen, PC expression is maintained in hepatopancreas, indicating its key role even in energy-limited conditions. Finally, both glucose and lactate concentrations were maintained under hypoxia for 24-48 h in hepatopancreas.

CFVisual: an interactive desktop platform for drawing gene structure and protein architecture

BACKGROUND

When researchers perform gene family analysis, they often analyze the structural characteristics of the gene, such as the distribution of introns and exons. At the same time, characteristic structural analysis of amino acid sequence is also essential, for example, motif and domain features. Researchers often integrate these analyses into one image to dig out more information, but the tools responsible for this integration are lacking.

RESULTS

Here, we developed a tool (CFVisual) for drawing gene structure and protein architecture. CFVisual can draw the phylogenetic tree, gene structure, and protein architecture in one picture, and has rich interactive capabilities, which can meet the work needs of researchers. Furthermore, it also supports arbitrary stitching of the above analysis images. It has become a useful helper in gene family analysis. The CFVisual package was implemented in Python and is freely available from https://github.com/ChenHuilong1223/CFVisual/ .

CONCLUSION

CFVisual has been used by some researchers and cited by some articles. In the future, CFVisual will continue to serve as a good helper for researchers in the study of gene structure and protein architecture.

Genome-wide identification and expression profile of GhGRF gene family in Gossypium hirsutum L

Background

Cotton is the primary source of renewable natural fiber in the textile industry and an important biodiesel crop. Growth regulating factors (GRFs) are involved in regulating plant growth and development.

Methods

Using genome-wide analysis, we identified 35 GRF genes in Gossypium hirsutum.

Results

Chromosomal location information revealed an uneven distribution of GhGRF genes, with maximum genes on chromosomes A02, A05, and A12 from the At sub-genome and their corresponding D05 and D12 from the Dt sub-genome. In the phylogenetic tree, 35 GRF genes were divided into five groups, including G1, G2, G3, G4, and G5. The majority of GhGRF genes have two to three introns and three to four exons, and their deduced proteins contained conserved QLQ and WRC domains in the N-terminal end of GRFs in Arabidopsis and rice. Sequence logos revealed that GRF genes were highly conserved during the long-term evolutionary process. The CDS of the GhGRF gene can complement MiRNA396a. Moreover, most GhGRF genes transcripts developed high levels of ovules and fibers. Analyses of promoter cis-elements and expression patterns indicated that GhGRF genes play an essential role in regulating plant growth and development by coordinating the internal and external environment and multiple hormone signaling pathways. Our analysis indicated that GhGRFs are ideal target genes with significant potential for improving the molecular structure of cotton.

NPM2 in malignant peritoneal mesothelioma: from basic tumor biology to clinical medicine

BACKGROUND

This review systematically summarizes gene biology features and protein structure of nucleoplasmin2 (NPM2) and the relationship between NPM2 and malignant peritoneal mesothelioma (MPM), in order to explore the molecular pathological mechanism of MPM and explore new therapeutic targets.

METHODS

NCBI PubMed database was used for the literature search. NCBI Gene and Protein databases, Ensembl Genome Browser, UniProt, and RCSB PDB database were used for gene and protein review. Three online tools (Consurf, DoGSiteScorer, and ZdockServer), the GEPIA database, and the Cancer Genome Atlas were used to analyze bioinformatics characteristics for NPM2 protein.

RESULTS

The main structural domains of NPM2 protein include the N-terminal core region, acidic region, and motif and disordered region. The N-terminal core region, involved in histone binding, is the most conserved domain in the nucleoplasmin (NPM) family. NPM2 with a large acidic tract in its C-terminal tail (NPM2-A2) is able to bind histones and form large complexes. Bioinformatics results indicated that NPM2 expression was correlated with the pathology of multiple tumors. Among mesothelioma patients, 5-year survival of patients with low-NPM2-expression was significantly higher than that of the high-NPM2-expression patients. NPM2 can facilitate the formation of histone deacetylation. NPM2 may promote histone deacetylation and inhibit the related-gene transcription, thus leading to abnormal proliferation, invasion, and metastasis of MPM.

CONCLUSION

NPM2 may play a key role in the development and progression of MPM.

Identification of ARF family in blueberry and its potential involvement of fruit development and pH stress response

Background

Auxin responsive factor (ARF) family is one of core components in auxin signalling pathway, which governs diverse developmental processes and stress responses. Blueberry is an economically important berry-bearing crop and prefers to acidic soil. However, the understandings of ARF family has not yet been reported in blueberry.

Results

In the present study, 60 ARF genes (VcARF) were identified in blueberry, and they showed diverse gene structures and motif compositions among the groups and similar within each group in the phylogenetic tree. Noticeably, 9 digenic, 5 trigenic and 6 tetragenic VcARF pairs exhibited more than 95% identity to each other. Computational analysis indicated that 23 VcARFs harbored the miRNA responsive element (MRE) of miR160 or miR167 like other plant ARF genes. Interestingly, the MRE of miR156d/h-3p was observed in the 5’UTR of 3 VcARFs, suggesting a potentially novel post-transcriptional control. Furthermore, the transcript accumulations of VcARFs were investigated during fruit development, and three categories of transcript profiles were observed, implying different functional roles. Meanwhile, the expressions of VcARFs to different pH conditions (pH4.5 and pH6.5) were surveyed in pH-sensitive and tolerant blueberry species, and a number of VcARFs showed different transcript accumulations. More importantly, distinct transcriptional response to pH stress (pH6.5) were observed for several VcARFs (such as VcARF6s and VcARF19-3/19–4) between pH-sensitive and tolerant species, suggesting their potential roles in adaption to pH stress.

Conclusions

Sixty VcARF genes were identified and characterized, and their transcript profiles were surveyed during fruit development and in response to pH stress. These findings will contribute to future research for eliciting the functional roles of VcARFs and regulatory mechanisms, especially fruit development and adaption to pH stress.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08556-y.

Molecular identification and expression analysis of five sucrose synthase genes in Sorghum Bicolor

In higher plants, sucrose synthase (Susy, EC 2.4.1.13) as an enzyme with a core function, involved in the synthesis and breakdown of sugars, and plays an important role in growth and metabolism. Although, the different genes encoding Susy isozyme proteins have been cloned and functionally verified in several plant species, to date detailed information about the Susy genes is lacking in Sorghum. Here, we demonstrated the identification of five novel Susy genes from the sorghum genome database. Sequence, structure and phylogenetic analyses of these five SbSusy genes revealed evolutionary conservation through Susy gene family members across Sorghum and other crop plants. The expression of sorghum Susy genes was investigated via transcriptome database in various developmental stages and different tissues. Further qRT-PCR was performed to reveal the induction of SbSusy genes under salt, drought and sugar induction. The results indicated that all Susy genes were differentially expressed in various tissues and highly associated with sucrose metabolism. This study shows a theoretical reference of Susy genes in Sorghum, which provides new insights for the knowledge of the evolution relationships, and basic information to help clarify the molecular mechanism of Susy synthase genes in Sorghum.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-022-01166-8.

Identifying genes with conserved splicing structure and orthologous isoforms in human, mouse and dog

BACKGROUND

In eukaryote transcriptomes, a significant amount of transcript diversity comes from genes' capacity to generate different transcripts through alternative splicing. Identifying orthologous alternative transcripts across multiple species is of particular interest for genome annotators. However, there is no formal definition of transcript orthology based on the splicing structure conservation. Likewise there is no public dataset benchmark providing groups of orthologous transcripts sharing a conserved splicing structure.

RESULTS

We introduced a formal definition of splicing structure orthology and we predicted transcript orthologs in human, mouse and dog. Applying a selective strategy, we analyzed 2,167 genes and their 18,109 known transcripts and identified a set of 253 gene orthologs that shared a conserved splicing structure in all three species. We predicted 6,861 transcript CDSs (coding sequence), mainly for dog, an emergent model species. Each predicted transcript was an ortholog of a known transcript: both share the same CDS splicing structure. Evidence for the existence of the predicted CDSs was found in external data.

CONCLUSIONS

We generated a dataset of 253 gene triplets, structurally conserved and sharing all their CDSs in human, mouse and dog, which correspond to 879 triplets of spliced CDS orthologs. We have released the dataset both as an SQL database and as tabulated files. The data consists of the 879 CDS orthology groups with their detailed splicing structures, and the predicted CDSs, associated with their experimental evidence. The 6,861 predicted CDSs are provided in GTF files. Our data may contribute to compare highly conserved genes across three species, for comparative transcriptomics at the isoform level, or for benchmarking splice aligners and methods focusing on the identification of splicing orthologs. The data is available at https://data-access.cesgo.org/index.php/s/V97GXxOS66NqTkZ .

The impacts of allopolyploidization on Methyl-CpG-Binding Domain (MBD) gene family in Brassica napus

BACKGROUND

Polyploidization promotes species formation and is widespread in angiosperms. Genome changes dramatically bring opportunities and challenges to plants after polyploidy. Methyl-CpG-Binding Domain (MBD) proteins can recognize and bind to methylation sites and they play an important role in the physiological process related to methylation in animals and plants. However, research on the influence of the allopolyploidization process on the MBD gene family is still lacking, so it is necessary to conduct a comprehensive analysis.

RESULTS

In this study, twenty-two, ten and eleven MBD genes were identified in the genome of allotetraploid B. napus and its diploid ancestors, B. rapa and B. oleracea, respectively. Based on the clades of the MBD gene in Arabidopsis, rice and maize, we divided the new phylogenetic tree into 8 clades. Among them, the true MBD genes in Brassica existed in only 5 clades. Clade IV and Clade VI were unique in term of MBD genes in dicotyledons. Ka/Ks calculations showed that MBD genes underwent purifying selection in Brassica and may retain genes through sequence or functional differentiation early in evolution. In the process of allopolyploidization, the number of MBD gene introns increased, and the protein motifs changed. The MBD proteins had their own special motifs in each clade, and the MBD domains were only conserved in their clades. At the same time, the MBD genes were expressed in flower, leaf, silique, and stem tissues, and the expression levels of the different genes were significantly different, while the tissue specificity was not obvious. The allopolyploidization process may increase the number of cis-acting elements and activate the transposable elements. During allopolyploidization, the expression pattern of the MBD gene changes, which may be regulated by cis-acting elements and transposable elements. The number imbalance of cis-acting elements and transposable elements in A_n and C_n subgenomes may also lead to biased A_n subgenome expression of the MBD gene in B. napus.

CONCLUSIONS

In this study, by evaluating the number, structure, phylogeny and expression of the MBD gene in B. napus and its diploid ancestors, we increased the understanding of MBD genes in allopolyploids and provided a reference for future analysis of allopolyploidization.

De Novo Assembly of the Nearly Complete Fathead Minnow Reference Genome Reveals a Repetitive but Compact Genome

The fathead minnow is a widely used model organism in environmental toxicology. The lack of a high-quality fathead minnow reference genome, however, has severely hampered its uses in toxicogenomics. We present the de novo assembly and annotation of the fathead minnow genome using long PacBio reads, Bionano and Hi-C scaffolding data, and large RNA-sequencing data sets from different tissues and life stages. The new annotated fathead minnow reference genome has a scaffold N50 of 12.0 Mbp and a complete benchmarking universal single-copy orthologs score of 95.1%. The completeness of annotation for the new reference genome is comparable to that of the zebrafish GRCz11 reference genome. The fathead minnow genome, revealed to be highly repetitive and sharing extensive syntenic regions with the zebrafish genome, has a much more compact gene structure than the zebrafish genome. Particularly, comparative genomic analysis with zebrafish, mouse, and human showed that fathead minnow homologous genes are relatively conserved in exon regions but had strikingly shorter intron regions. The new fathead minnow reference genome and annotation data, publicly available from the National Center for Biotechnology Information and the University of California Santa Cruz genome browser, provides an essential resource for aquatic toxicogenomic studies in ecotoxicology and public health. Environ Toxicol Chem 2022;41:448-461. Published 2021. This article is a U.S. Government work and is in the public domain in the USA.

Comparative Analysis of Complete Mitogenomes of Two Oxyurichthys Gobies and Their Phylogenetic Implication

Oxyurichthys is a genus of goby that is widespread in the tropical Indo-West Pacific region. Oxyurichthys species are usually found in estuarine and coastal marine habitats. In Southeast Asia, they are commercial fishes and often collected by trawling to serve the market's demand. The mitogenome serves as a good marker for investigating the systematics and evolution of fishes, but the mitogenome of Oxyurichthys species remains unknown. In this study, mitogenomes of two Oxyurichthys gobies, O. ophthalmonema and O. microlepis, were characterized and compared. The sizes of the mitogenomes were 16,504 bp and 16,506 bp for O. ophthalmonema and O. microlepis, respectively. Mitogenomes of these two species were similar in gene content and structure. Both included 37 genes and a control region. The two Oxyurichthys mitogenomes shared similar gene features and base composition with other recorded gobies. Typical conserved blocks (CSB-1, CSB-2, CSB-3 and CSB-D) were found in the control region of both species. Phylogenetic analyses based on concatenation of 13 protein-coding genes and 2 rRNAs revealed that the two Oxyurichthys species clustered together and were sister to species of the genera Sicydium, Sicyopterus and Stiphodon. The findings of the present study support previous evolutionary studies of gobies using other molecular markers.

Characterization of the complete mitochondrial genome of Diploastrea heliopora and phylogeny of the scleractinia species which have group I introns in their COI genes

Mitochondrial genome DNA is a powerful marker for resolving phylogenetic relationships among scleractinian corals. Here, we decode the complete mitochondrial genome of Diploastrea heliopora (Lamarck, 1816) for the first time. The general features are 18 363 bp in length, and conventionally, with 13 protein coding genes, two ribosomal RNAs, and two transfer RNAs. Gene arrangement and distribution are similar to other scleractinian corals. Moreover, the COI gene of D. heliopora is broken up into two parts by a complex group I intron. This intron is 1076 bases in length and contains helical structures (P1-P10, except P2) and four conserved regions (P, Q, R, and S). The mitochondrial genome of D. heliopora has asymmetric base composition (13.03% C, 20.29% G, 25.91% A, and 40.77% for T). Based on concatenated protein coding genes, ML and BI trees show similar phylogenetic relationship: D. heliopora clustered closely with Sclerophyllia maxima and Echinophyllia aspera into the robust branch. The data and conclusion in this study are reference for further phylogenetic studies of corals.

The Evolution of Hemocyanin Genes in Caenogastropoda: Gene Duplications and Intron Accumulation in Highly Diverse Gastropods

Hemocyanin is the oxygen transport protein of most molluscs and represents an important physiological factor that has to be well-adapted to their environments because of the strong influences of abiotic factors on its oxygen affinity. Multiple independent gene duplications and intron gains have been reported for hemocyanin genes of Tectipleura (Heterobranchia) and the caenogastropod species Pomacea canaliculata, which contrast with the uniform gene architectures of hemocyanins in Vetigastropoda. The goal of this study was to analyze hemocyanin gene evolution within the diverse group of Caenogastropoda in more detail. Our findings reveal multiple gene duplications and intron gains and imply that these represent general features of Apogastropoda hemocyanins. Whereas hemocyanin exon–intron structures are identical within different Tectipleura lineages, they differ strongly within Caenogastropoda among phylogenetic groups as well as between paralogous hemocyanin genes of the same species. Thus, intron accumulation took place more gradually within Caenogastropoda but finally led to a similar consequence, namely, a multitude of introns. Since both phenomena occurred independently within Heterobranchia and Caenogastropoda, the results support the hypothesis that introns may contribute to adaptive radiation by offering new opportunities for genetic variability (multiple paralogs that may evolve differently) and regulation (multiple introns). Our study indicates that adaptation of hemocyanin genes may be one of several factors that contributed to the evolution of the large diversity of Apogastropoda. While questions remain, this hypothesis is presented as a starting point for the further study of hemocyanin genes and possible correlations between hemocyanin diversity and adaptive radiation.

Genome-wide identification and analysis of chemokine receptor superfamily in miiuy croaker, Miichthys miiuy

The chemokine receptor (ChemR) superfamily, which is divided into 4 subfamilies (CXCR, CCR, XCR, and CX3CR), is the main receptors of chemokines in innate immune responses. In the current study, we have identified 27 ChemRs in miiuy croaker: 13 CCR genes, 11 CXCR genes, and 3 XCR genes. Multiple characteristics of these genes, including phylogeny, gene structures, conserved motifs, chromosome locations, evolutionary mechanism, and expression levels upon the bacterial challenge were analyzed. Gene structure and location analysis showed that all ChemR genes contain fewer introns (≤4) and they are unevenly distributed on the 12 chromosomes. And the XCR subfamily of miiuy croaker don't have the DRY motif of ChemR. Phylogenetic and synteny analysis showed that these genes experienced tandem and segmental duplication event in several species, and tandem duplication might be the main expansion way in miiuy croaker. The major ChemRs of each orthologous group in vertebrates were selected for molecular evolution analysis, the results of which indicated that compared with vertebrates, ChemRs of teleost fishes may have a relatively high evolutionary dynamic. In addition, a total of 21 positively selected codons were detected in vertebrate ChemRs under Model 8. RNA-Seq analysis and qRT-PCR verification demonstrated that CXCR3.2, CXCR5, and XCR1 genes were up-regulated significantly upon the Vibrio harveyi infection. These results provide valuable information for investigating the evolutionary relationships of chemokine receptor superfamily in miiuy croaker and laid the basis for further functional analysis.

Genome-wide analysis of HSP70 gene superfamily in Pyropia yezoensis (Bangiales, Rhodophyta): identification, characterization and expression profiles in response to dehydration stress

BACKGROUND

Heat shock proteins (HSPs) perform a fundamental role in protecting plants against abiotic stresses. Individual family members have been analyzed in previous studies, but there has not yet been a comprehensive analysis of the HSP70 gene family in Pyropia yezoensis.

RESULTS

We investigated 15 putative HSP70 genes in Py. yezoensis. These genes were classified into two sub-families, denoted as DnaK and Hsp110. In each sub-family, there was relative conservation of the gene structure and motif. Synteny-based analysis indicated that seven and three PyyHSP70 genes were orthologous to HSP70 genes in Pyropia haitanensis and Porphyra umbilicalis, respectively. Most PyyHSP70s showed up-regulated expression under different degrees of dehydration stress. PyyHSP70-1 and PyyHSP70-3 were expressed in higher degrees compared with other PyyHSP70s in dehydration treatments, and then expression degrees somewhat decreased in rehydration treatment. Subcellular localization showed PyyHSP70-1-GFP and PyyHSP70-3-GFP were in the cytoplasm and nucleus/cytoplasm, respectively. Similar expression patterns of paired orthologs in Py. yezoensis and Py. haitanensis suggest important roles for HSP70s in intertidal environmental adaptation during evolution.

CONCLUSIONS

These findings provide insight into the evolution and modification of the PyyHSP70 gene family and will help to determine the functions of the HSP70 genes in Py. yezoensis growth and development.

Genome-wide analysis and expression patterns of lipid phospholipid phospholipase gene family in Brassica napus L

BACKGROUND

Lipid phosphate phosphatases (LPP) are critical for regulating the production and degradation of phosphatidic acid (PA), an essential signaling molecule under stress conditions. Thus far, the LPP family genes have not been reported in rapeseed (Brassica napus L.).

RESULTS

In this study, a genome-wide analysis was carried out to identify LPP family genes in rapeseed that respond to different stress conditions. Eleven BnLPPs genes were identified in the rapeseed genome. Based on phylogenetic and synteny analysis, BnLPPs were classified into four groups (Group I-Group IV). Gene structure and conserved motif analysis showed that similar intron/exon and motifs patterns occur in the same group. By evaluating cis-elements in the promoters, we recognized six hormone- and seven stress-responsive elements. Further, six putative miRNAs were identified targeting three BnLPP genes. Gene ontology analysis disclosed that BnLPP genes were closely associated with phosphatase/hydrolase activity, membrane parts, phosphorus metabolic process, and dephosphorylation. The qRT-PCR based expression profiles of BnLPP genes varied in different tissues/organs. Likewise, several gene expression were significantly up-regulated under NaCl, PEG, cold, ABA, GA, IAA, and KT treatments.

CONCLUSIONS

This is the first report to describe the comprehensive genome-wide analysis of the rapeseed LPP gene family. We identified different phytohormones and abiotic stress-associated genes that could help in enlightening the plant tolerance against phytohormones and abiotic stresses. The findings unlocked new gaps for the functional verification of the BnLPP gene family during stresses, leading to rapeseed improvement.

Exon-1 skipping and intron-1 retaining by alternative splicing of the c-KIT gene encodes a novel splice variant in the skin of Merino sheep (Ovis aries)

c-KIT, a type III receptor protein tyrosine kinase, plays an essential role in melanocyte development, migration, and survival. Mutations within the c-KIT gene are previously shown to cause the white coat color phenotypes in pigs, mice, goats, and humans. However, up so far, the splicing isoform(s), genomic architecture of c-KIT have not been characterized well in merino sheep. Reverse transcriptase (RT)-PCR analysis with molecular prediction identified two basic splice variants: Transcript Variant-1, 2 for 12 bp insertion coding sequences (CDS) corresponding to the four amino acids 'GNSK', respectively. Using 5' RACE, here we report for the first time a novel c-KIT 'Transcript Variant-3' from the skin of merino sheep by comparative genome analyses at exon(1)-intron(1)-exon(2) boundaries. In contrast, a single product of 795 bp was characterized by 3' RACE. We also demonstrated that the c-KIT gene expression at the transcript level is not mediated via an intron-9 splicing event. Overall, beyond what was observed in other mammals, our data provide novel insights into the molecular structure of the c-KIT gene in sheep.

Identification and Expression Analysis of miR160 and Their Target Genes in Cucumber

miR160 plays a crucial role in various biological processes by regulating their target gene auxin response factor (ARF) in plants. However, little is known about miR160 and ARF in cucumber fruit expansion. Here, 4 Csa-MIR160 family members and 17 CsARFs were identified through a genome-wide search. Csa-miR160 showed a closer relationship with those in melon. Phylogenetic analysis revealed that CsARFs were divided into four classes and most of CsARFs presented a closer evolutionary relationship with those from tomato. Putative cis-elements analysis predicted that Csa-MIR160 and CsARFs were involved in light, phytohormone and stress response, which proved that they might take part in light, phytohormone and stress signaling pathway by the miR160-ARF module. In addition, CsARF5, CsARF11, CsARF13 and CsARF14 were predicted as the target genes of Csa-miR160. qRT-PCR revealed that Csa-miR160 and their target gene CsARFs were differentially expressed in differential cucumber tissues and developmental stages. Csa-miR160d was only expressed in the expanded cucumber fruit. CsARF5, CsARF11 and CsARF13 exhibited the lower expression in the expanded fruit than those in the ovary, while, CsARF14 showed the reverse trend. Our results suggested that Csa-miR160d might play a crucial role in cucumber fruit expansion by negatively targeting CsARF5, CsARF11 and CsARF13. This is the first genome-wide analysis of miR160 in cucumber. These findings provide useful information and resources for further studying the role of miR160 and ARF in cucumber fruit expansion.

Genome-wide identification and characterization of basic helix-loop-helix genes in nine molluscs

The basic helix-loop-helix (bHLH) transcription factors form a large superfamily that plays an important role in numerous physiological processes, including development and response to environmental stresses. In this study, the distribution of bHLH genes in nine molluscs was systematically investigated (including five bivalves, three gastropods and one cephalopod). Finally, 53-85 bHLH genes were identified from each genome and classified into corresponding families by using phylogenetic analysis. The results of gene structure and conserved motif analysis illustrated the hereditary conservation of bHLH transcription factors during evolution but showed low similarity in group C. Through transcription profile analysis of C. gigas and T. granosa, we found a important role of bHLH genes in responding to multiple external challenges and development; meanwhile, they also exhibited tissue-specific expression. Interestingly, we were also surprised to find different bHLH genes from the same group generally possess similar patterns expression that tends to simultaneously present high or lower expression of multiple challenges and different tissues in this study. In summary, this study lays the foundation for further investigation of the biological functions and evolution of molluscan bHLH genes.

The evolution of hemocyanin genes in Tectipleura: a multitude of conserved introns in highly diverse gastropods

Background

Hemocyanin is the oxygen transporter of most molluscs. Since the oxygen affinity of hemocyanin is strongly temperature-dependent, this essential protein needs to be well-adapted to the environment. In Tectipleura, a very diverse group of gastropods with > 27,000 species living in all kinds of habitats, several hemocyanin genes have already been analyzed. Multiple independent duplications of this gene have been identified and may represent potential adaptations to different environments and lifestyles. The aim of this study is to further explore the evolution of these genes by analyzing their exon–intron architectures.

Results

We have reconstructed the gene architectures of ten hemocyanin genes from four Tectipleura species: Aplysia californica, Lymnaea stagnalis, Cornu aspersum and Helix pomatia. Their hemocyanin genes each contain 53 introns, significantly more than in the hemocyanin genes of Cephalopoda (9–11), Vetigastropoda (15) and Caenogastropoda (28–33). The gene structures of Tectipleura hemocyanins are identical in terms of intron number and location, with the exception of one out of two hemocyanin genes of L. stagnalis that comprises one additional intron. We found that gene structures that differ between molluscan lineages most probably evolved more recently through independent intron gains.

Conclusions

The strict conservation of the large number of introns in Tectipleura hemocyanin genes over 200 million years suggests the influence of a selective pressure on this gene structure. While we could not identify conserved sequence motifs within these introns, it may be simply the great number of introns that offers increased possibilities of gene regulation relative to hemocyanin genes with less introns and thus may have facilitated habitat shifts and speciation events. This hypothesis is supported by the relatively high number of introns within the hemocyanin genes of Pomacea canaliculata that has evolved independently of the Tectipleura. Pomacea canaliculata belongs to the Caenogastropoda, the sister group of Heterobranchia (that encompass Tectipleura) which is also very diverse and comprises species living in different habitats. Our findings provide a hint to some of the molecular mechanisms that may have supported the spectacular radiation of one of Metazoa’s most species rich groups.

2020 update on human coronaviruses: One health, one world

Since human coronavirus (HCoVs) was first described in the 1960s, seven strains of respiratory human coronaviruses have emerged and caused human infections. After the emergence of severe acute respiratory syndrome coronavirus (SARS-CoV) and Middle East respiratory syndrome coronavirus (MERS-CoV), a pneumonia outbreak of coronavirus disease 2019 (COVID-19) caused by a novel coronavirus (SARS-CoV-2) has represented a pandemic threat to global public health in the 21st century. Without effectively prophylactic and therapeutic strategies including vaccines and antiviral drugs, these three coronaviruses have caused severe respiratory syndrome and high case-fatality rates around the world. In this review, we detail the emergence event, origin and reservoirs of all HCoVs, compare the differences with regard to structure and receptor usage, and summarize therapeutic strategies for COVID-19 that cause severe pneumonia and global pandemic.

Genome-wide identification and analysis of the thiolase family in insects

Thiolases are important enzymes involved in lipid metabolism in both prokaryotes and eukaryotes, and are essential for a range of metabolic pathways, while, little is known for this important family in insects. To shed light on the evolutionary models and functional diversities of the thiolase family, 137 thiolase genes were identified in 20 representative insect genomes. They were mainly classified into five classes, namely cytosolic thiolase (CT-thiolase), T1-thiolase, T2-thiolase, trifunctional enzyme thiolase (TFE-thiolase), and sterol carrier protein 2 thiolase (SCP2-thiolase). The intron number and exon/intron structures of the thiolase genes reserve large diversification. Subcellular localization prediction indicated that all the thiolase proteins were mitochondrial, cytosolic, or peroxisomal enzymes. Four highly conserved sequence fingerprints were found in the insect thiolase proteins, including CxS-, NEAF-, GHP-, and CxGGGxG-motifs. Homology modeling indicated that insect thiolases share similar 3D structures with mammals, fishes, and microorganisms. In Bombyx mori, microarray data and reverse transcription-polymerase chain reaction (RT-PCR) analysis suggested that some thiolases might be involved in steroid metabolism, juvenile hormone (JH), and sex pheromone biosynthesis pathways. In general, sequence and structural characteristics were relatively conserved among insects, bacteria and vertebrates, while different classes of thiolases might have differentiation in specific functions and physiological processes. These results will provide an important foundation for future functional validation of insect thiolases.

Molecular evolution of melatonin receptor genes (mtnr) in vertebrates and its shedding light on mtnr1c

Melatonin receptors (MTNRs) play important roles in regulation of circadian rhythms and seasonal reproduction. However, their origin and evolution in vertebrates have not been investigated. Here, we performed a comprehensive examination by comparative genome mining of MTNRs in vertebrates. We successfully extracted 164 putative encoding sequences for MTNRs (including 57 mtnr1a, 59 mtnr1b and 48 mtnr1c) from 45 high-quality representative genomes. Interestingly, the putative expansions of mtnr1a and mtnr1b in zebrafish were also identified in other Cyprinifomes, but not in other orders of teleost. Using phylogenetic interference, we observed this expansion to be clustered into a primitive position of the Actinopterygii, which may be resulted from teleost-specific genome duplication. The C-terminal extension of MTNR1C, predicted to be proteoglycan 4 (PRG4), originated after the speciation of Monotremata or Marsupialia. Our present genomics survey provides novel insights into the evolution of MTNRs in vertebrates and updates our understanding of these proteins.

Genome-wide identification, genomic organization and expression profiles of SlARR-B gene family in tomato

The type-B authentic response regulators (ARR-Bs) function as positive regulators of cytokinin signal transduction and play important roles in abiotic stress resistance and plant development. However, little of ARR-B family is known in tomato. In this study, we performed a comprehensive analysis of ARR-B family factors in tomato. In total, 12 genes encoding ARR-B transcription factors (named as SlARR-B1-SlARR-B12) were identified from tomato. We analyzed the structures, chromosome locations, phylogeny, protein motifs, and expression profiles of these SlARR-B genes. Gene structure analysis showed that 5-12 exons and 4-11 introns existed in the SlARR-B genes. These SlARR-B genes were asymmetrically distributed on eight chromosomes in tomato. Phylogenetic tree of SlARR-B genes from tomato and other plant species revealed that SlARR-B genes were classified into 6 subfamilies. SlARR-B proteins had typical conserved domains, including Motif 1 and Motif 2. The investigation of the expression profiles of SlARR-B genes in all the examined tissues demonstrated that these genes were differentially expressed, including roots, stems, leaves, flowers, and fruits at developmental stages. Notably, the expression of SlARR-B11 and SlARR-B12 exhibited high expression levels in flowers. Each gene was induced by at least one of different phytohormones (SA, IAA, ABA, IBA, 6-BA, JA, GA, and ETH) and four abiotic stress treatments (heat, drought, salt, and cold). This study sets a good foundation for further characterization of the SlARR-B transcription factors in plant development and abiotic stress responses of tomato.

Genome-wide characterization of the abscisic acid-, stress- and ripening-induced (ASR) gene family in wheat (Triticum aestivum L.)

BACKGROUND

Abscisic acid-, stress-, and ripening-induced (ASR) genes are a class of plant specific transcription factors (TFs), which play important roles in plant development, growth and abiotic stress responses. The wheat ASRs have not been described in genome-wide yet.

METHODS

We predicted the transmembrane regions and subcellular localization using the TMHMM server, and Plant-mPLoc server and CELLO v2.5, respectively. Then the phylogeny tree was built by MEGA7. The exon-intron structures, conserved motifs and TFs binding sites were analyzed by GSDS, MEME program and PlantRegMap, respectively.

RESULTS

In wheat, 33ASR genes were identified through a genome-wide survey and classified into six groups. Phylogenetic analyses revealed that the TaASR proteins in the same group tightly clustered together, compared with those from other species. Duplication analysis indicated that the TaASR gene family has expanded mainly through tandem and segmental duplication events. Similar gene structures and conserved protein motifs of TaASRs in wheat were identified in the same groups. ASR genes contained various TF binding cites associated with the stress responses in the promoter region. Gene expression was generally associated with the expected group-specific expression pattern in five tissues, including grain, leaf, root, spike and stem, indicating the broad conservation of ASR genes function during wheat evolution. The qRT-PCR analysis revealed that several ASRs were up-regulated in response to NaCl and PEG stress.

CONCLUSION

We identified ASR genes in wheat and found that gene duplication events are the main driving force for ASR gene evolution in wheat. The expression of wheat ASR genes was modulated in responses to multiple abiotic stresses, including drought/osmotic and salt stress. The results provided important information for further identifications of the functions of wheat ASR genes and candidate genes for high abiotic stress tolerant wheat breeding.

Molecular and immunological characterizations of interleukin-11 in large yellow croaker (Larimichthys crocea)

Interleukin (IL)-11 is a multifunctional cytokine that exerts a series of important immunomodulatory effects and exists in many tissues and cells. A 1106-bp nucleotide sequence representing the complete cDNA of IL-11 was obtained from large yellow croaker (Larimichthys crocea), containing an open reading frame (ORF) of 603 bp encoding for 200 amino acids (aa). The predicted LcIL-11 protein included a 12aa signal peptide and a conserved IL-11 domain. The polypeptide sequence identities between LcIL-11 and its counterparts in mammals and other fish are from 84% to 92% with known fish IL-11a and 22%-27% with fish IL-11b. LcIL-11 mRNA existed in most tissues with the most predominant expression in the gill. After immune challenge, the expression levels of LcIL-11 were induced largely in vivo and in vitro, with the peak-value of 32 times as much as the control in the liver at 24 h after Vibrio parahaemolyticus injection (p < 0.05) and the greatest value of 13.9 times as much as the control in LCK cells at 12 h after poly I:C stimulation (p < 0.05). Furthermore, the overexpression vector pcDNA3.1-LcIL-11 was constructed and transfected to LCK cells. Our results showed that the transcriptional expression levels of tumor necrosis factor (TNF)-α and myxovirus resistant protein (Mx) significantly up-regulated in LCK cells after LcIL-11 overexpression (p < 0.05). However, no significant changes of IL-1β, janus kinase (JAK)2 and signal transducers and activators of transcription (STAT)5 was detected. Our finding indicated that LcIL-11 might enhance TNF-α and antiviral protein Mx expression in large yellow croaker.

Zmat2 in mammals: conservation and diversification among genes and Pseudogenes

BACKGROUND

Recent advances in genetics and genomics present unique opportunities for enhancing our understanding of mammalian biology and evolution through detailed multi-species comparative analysis of gene organization and expression. Yet, of the more than 20,000 protein coding genes found in mammalian genomes, fewer than 10% have been examined in any detail. Here we elucidate the power of data available in publicly-accessible genomic and genetic resources by querying them to evaluate Zmat2, a minimally studied gene whose human ortholog has been implicated in spliceosome function and in keratinocyte differentiation.

RESULTS

We find extensive conservation in coding regions and overall structure of Zmat2 in 18 mammals representing 13 orders and spanning ~ 165 million years of evolutionary development, and in their encoded proteins. We identify a tandem duplication in the Zmat2 gene and locus in opossum, but not in other monotremes, marsupials, or other mammals, indicating that this event occurred subsequent to the divergence of these species from one another. We also define a collection of Zmat2 pseudogenes in half of the mammals studied, and suggest based on phylogenetic analysis that they each arose independently in the recent evolutionary past.

CONCLUSIONS

Mammalian Zmat2 genes and ZMAT2 proteins illustrate conservation of structure and sequence, along with the development and diversification of pseudogenes in a large fraction of species. Collectively, these observations also illustrate how the focused identification and interpretation of data found in public genomic and gene expression resources can be leveraged to reveal new insights of potentially high biological significance.

Characterization of two kcnk3 genes in Nile tilapia (Oreochromis niloticus): Molecular cloning, tissue distribution, and transcriptional changes in various salinity of seawater

As one important member of the two-pore-domain potassium channel (K_2P) family, potassium channel subfamily K member 3 (KCNK3) has been reported for thermogenesis regulation, energy homeostasis, membrane potential conduction, and pulmonary hypertension in mammals. However, its roles in fishes are far less examined and published. In the present study, we identified two kcnk3 genes (kcnk3a and kcnk3b) in an euryhaline fish, Nile tilapia (Oreochromis niloticus), by molecular cloning, genomic survey and laboratory experiments to investigate their potential roles for osmoregulation. We obtained full-length coding sequences of the kcnk3a and kcnk3b genes (1209 and 1173 bp), which encode 402 and 390 amino acids, respectively. Subsequent multiple sequence alignments, putative 3D-structure model prediction, genomic survey and phylogenetic analysis confirmed that two kcnk3 paralogs are widely presented in fish genomes. Interestingly, a DNA fragment inversion of a kcnk3a cluster was found in Cypriniforme in comparison with other fishes. Quantitative real-time PCRs demonstrated that both the tilapia kcnk3 genes were detected in all the examined tissues with a similar distribution pattern, and the highest transcriptions were observed in the heart. Meanwhile, both kcnk3 genes in the gill were proved to have a similar transcriptional change pattern in response to various salinity of seawater, implying that they might be involved in osmoregulation. Furthermore, three predicted transcription factors (arid3a, arid3b, and arid5a) of both kcnk3 genes also showed a similar pattern as their target genes in response to the various salinity, suggesting their potential positive regulatory roles. In summary, we for the first time characterized the two kcnk3 genes in Nile tilapia, and demonstrated their potential involvement in osmoregulation for this economically important fish.

Genome-wide identification and transcriptional expression analysis of superoxide dismutase (SOD) family in wheat (Triticum aestivum)

Superoxide dismutases (SODs) are a family of key antioxidant enzymes that play a crucial role in plant growth and development. Previously, this gene family has been investigated in Arabidopsis and rice. In the present study, a genome-wide analysis of the SOD gene family in wheat were performed. Twenty-six SOD genes were identified from the whole genome of wheat, including 17 Cu/Zn-SODs, six Fe-SODs, and three Mn-SODs. The chromosomal location mapping analysis indicated that these three types of SOD genes were only distributed on 2, 4, and 7 chromosomes, respectively. Phylogenetic analyses of wheat SODs and several other species revealed that these SOD proteins can be assigned to two major categories. SOD1 mainly comprises of Cu/Zn-SODs, and SOD2 mainly comprises of Fe-SODs and Mn-SODs. Gene structure and motif analyses indicated that most of the SOD genes showed a relatively conserved exon/intron arrangement and motif composition. Analyses of transcriptional data indicated that most of the wheat SOD genes were expressed in almost all of the examined tissues and had important functions in abiotic stress resistance. Finally, quantitative real-time polymerase chain reaction (qRT-PCR) analysis was used to reveal the regulating roles of wheat SOD gene family in response to NaCl, mannitol, and polyethylene glycol stresses. qRT-PCR showed that eight randomly selected genes with relatively high expression levels responded to all three stresses based on released transcriptome data. However, their degree of response and response patterns were different. Interestingly, among these genes, TaSOD1.7, TaSOD1.9, TaSOD2.1, and TaSOD2.3 feature research value owing to their remarkable expression-fold change in leaves or roots under different stresses. Overall, our results provide a basis of further functional research on the SOD gene family in wheat and facilitate their potential use for applications in the genetic improvement on wheat in drought and salt stress environments.

Tick mitochondrial genomes: structural characteristics and phylogenetic implications

Ticks are obligate blood-sucking arachnid ectoparasites from the order Acarina, and many are notorious as vectors of a wide variety of zoonotic pathogens. However, the systematics of ticks in several genera is still controversial. The mitochondrial genome (mt-genome) has been widely used in arthropod phylogeny, molecular evolution and population genetics. With the development of sequencing technologies, an increasing number of tick mt-genomes have been sequenced and annotated. To date, 63 complete tick mt-genomes are available in the NCBI database, and these genomes have become an increasingly important genetic resource and source of molecular markers in phylogenetic studies of ticks in recent years. The present review summarizes all available complete mt-genomes of ticks in the NCBI database and analyses their characteristics, including structure, base composition and gene arrangement. Furthermore, a phylogenetic tree was constructed using mitochondrial protein-coding genes (PCGs) and ribosomal RNA (rRNA) genes from ticks. The results will provide important clues for deciphering new tick mt-genomes and establish a foundation for subsequent taxonomic research.

Genome-wide characterization and expression analysis of aquaporins in salt cress (Eutrema salsugineum)

Aquaporins (AQPs) serve as water channel proteins and belong to major intrinsic proteins (MIPs) family, functioning in rapidly and selectively transporting water and other small solutes across biological membranes. Importantly, AQPs have been shown to play a critical role in abiotic stress response pathways of plants. As a species closely related to Arabidopsis thaliana, Eutrema salsugineum has been proposed as a model for studying salt resistance in plants. Here we surveyed 35 full-length AQP genes in E. salsugineum, which could be grouped into four subfamilies including 12 plasma membrane intrinsic proteins (PIPs), 11 tonoplast intrinsic proteins (TIPs), nine NOD-like intrinsic proteins (NIPs), and three small basic intrinsic proteins (SIPs) by phylogenetic analysis. EsAQPs were comprised of 237-323 amino acids, with a theoretical molecular weight (MW) of 24.31-31.80 kDa and an isoelectric point (pI) value of 4.73-10.49. Functional prediction based on the NPA motif, aromatic/arginine (ar/R) selectivity filter, Froger's position and specificity-determining position suggested quite differences in substrate specificities of EsAQPs. EsAQPs exhibited global expressions in all organs as shown by gene expression profiles and should be play important roles in response to salt, cold and drought stresses. This study provides comprehensive bioinformation on AQPs in E. salsugineum, which would be helpful for gene function analysis for further studies.

Genome-wide identification, structure characterization, and expression pattern profiling of aquaporin gene family in cucumber

BACKGROUND

Aquaporin (AQP) proteins comprise a group of membrane intrinsic proteins (MIPs) that are responsible for transporting water and other small molecules, which is crucial for plant survival under stress conditions including salt stress. Despite the vital role of AQPs, little is known about them in cucumber (Cucumis sativus L.).

RESULTS

In this study, we identified 39 aquaporin-encoding genes in cucumber that were separated by phylogenetic analysis into five sub-families (PIP, TIP, NIP, SIP, and XIP). Their substrate specificity was then assessed based on key amino acid residues such as the aromatic/Arginine (ar/R) selectivity filter, Froger's positions, and specificity-determining positions. The putative cis-regulatory motifs available in the promoter region of each AQP gene were analyzed and results revealed that their promoter regions contain many abiotic related cis-regulatory elements. Furthermore, analysis of previously released RNA-seq data revealed tissue- and treatment-specific expression patterns of cucumber AQP genes (CsAQPs). Three aquaporins (CsTIP1;1, CsPIP2;4, and CsPIP1;2) were the most transcript abundance genes, with CsTIP1;1 showing the highest expression levels among all aquaporins. Subcellular localization analysis in Nicotiana benthamiana epidermal cells revealed the diverse and broad array of sub-cellular localizations of CsAQPs. We then performed RNA-seq to identify the expression pattern of CsAQPs under salt stress and found a general decreased expression level of root CsAQPs. Moreover, qRT-PCR revealed rapid changes in the expression levels of CsAQPs in response to diverse abiotic stresses including salt, polyethylene glycol (PEG)-6000, heat, and chilling stresses. Additionally, transient expression of AQPs in N. benthamiana increased leaf water loss rate, suggesting their potential roles in the regulation of plant water status under stress conditions.

CONCLUSIONS

Our results indicated that CsAQPs play important roles in response to salt stress. The genome-wide identification and primary function characterization of cucumber aquaporins provides insight to elucidate the complexity of the AQP gene family and their biological functions in cucumber.

Tracing the origin and evolution history of methylation-related genes in plants

BACKGROUND

DNA methylation is a crucial epigenetic modification, which is involved in many biological processes, including gene expression regulation, embryonic development, cell differentiation and genomic imprinting etc. And it also involves many key regulatory genes in eukaryotes. By tracing the evolutionary history of methylation-related genes, we can understand the origin and expansion time of these genes, which helps to understand the evolutionary history of plants, and we can also understand the changes of DNA methylation patterns in different species. However, most studies on the evolution of methylation-related genes failed to be carried out for the whole DNA methylation pathway.

RESULTS

In this study, we conducted a comprehensive identification of 33 methylation-related genes in 77 species, and investigated gene origin and evolution throughout the plant kingdom. We found that the origin of genes responsible for methylation maintenance and demethylation evolved early, while most de novo methylation-related genes appeared late. The methylation-related genes were expanded by whole genome duplication and tandem replication, but were also accompanied by a large number of gene absence events in different species. The gene length and intron length varied a lot in different species, but exon structure and functional domains were relatively conserved. The phylogenetic relationships of methylation-related genes were traced to reveal the evolution history of DNA methylation in different species. The expression patterns of methylation-related genes have changed during the evolution of species, and the expression patterns of these genes in different species can be clustered into four categories.

CONCLUSIONS

The study describes a global characterization of DNA methylation-related genes in the plant kingdom. The similarities and differences in origin time, gene structure and phylogenetic relationship of these genes lead us to understand the evolutionary conservation and dynamics of DNA methylation in plants.

Genomic, molecular evolution, and expression analysis of NOX genes in soybean (Glycine max)

Reactive oxygen species (ROS) are versatile signaling molecules in sensing stresses and play critical roles in signaling and development. Plasma membrane NADPH oxidases (NOXs) are key producers of ROS, and play important roles in the regulation of plant-pathogen interactions. Here, we performed a comprehensive analysis of the NOX gene family in the soybean genome (Glycine max) and 17 NOX (GmNOX) genes were identified. Structural analysis revealed that the GmNOX proteins in soybean were as conserved as those in other plants. 8 duplicated gene pairs were formed by a Glycine-specific whole-genome duplication (WGD) event approximately 13 million years ago (Mya). The Ka/Ks ratios of GmNOX genes ranged from 0.04 to 0.28, suggesting that the GmNOX family had undergone purifying selection in soybean. Gene expression patterns showed different expression of these duplicate genes, suggesting that the GmNOXs were retained by substantial subfunctionalization during the soybean evolutionary processes. Subsequently, the expression of GmNOXs in response to drought and phytohormones were characterized via qPCR. Importantly, four GmNOXs showed strong expression in nodules, pointing to their probable involvement in nodulation. Thus, our results shed light on the evolutionary history of this family in soybean and contribute to the functional characterization of GmNOX genes in soybean.

Genome-wide identification and expression profile analysis of the Hsp20 gene family in Barley (Hordeum vulgare L.)

In plants, heat shock proteins (Hsps) play important roles in response to diverse stresses. Hsp20 is the major family of Hsps, but their role remains poorly understood in barley (Hordeum vulgare L.). To reveal the mechanisms of barley Hsp20s (HvHsp20s) response to stress conditions, we performed a comprehensive genome-wide analysis of the HvHsp20 gene family using bioinformatics-based methods. In total, 38 putative HvHsp20s were identified in barley and grouped into four subfamilies (C, CP, PX, and MT) based on predicted subcellular localization and their phylogenetic relationships. A sequence analysis indicated that most HvHsp20 genes have no intron or one with a relatively short length. In addition, the same group of HvHsp20 proteins in the phylogenetic tree shared similar gene structure and motifs, indicating that they were highly conserved and might have similar function. Based on RNA-seq data analysis, we showed that the transcript levels of HvHsp20 genes could be induced largely by abiotic and biotic stresses such as heat, salt, and powdery mildew. Three HvHsp20 genes, HORVU7Hr1G036540, HORVU7Hr1G036470, and HORVU3Hr1G007500, were up-regulated under biotic and abiotic stresses, suggesting their potential roles in mediating the response of barley plants to environment stresses. These results provide valuable information for further understanding the complex mechanisms of HvHsp20 gene family in barley.

Genome-wide identification and expression analysis of sucrose synthase genes in allotetraploid Brassica juncea

Sucrose plays pivotal role in energy metabolism and regulating gene expression of several physiological processes in higher plants. Here, fourteen sucrose synthase (SUS) genes have been identified in the allotetraploid genome of Indian mustard, Brassica juncea. The identified SUS genes in B. juncea (BjSUS) were derived from the two-progenitor species, B. rapa and B. nigra. Intron-exon analysis indicated loss or gain of 1-3 introns in diversification of SUS gene family. Phylogenetic analysis revealed discrete evolutionary paths for the BjSUS genes, originating from three ancestor groups, SUS I, SUS II and SUS III. Gene expression study revealed significant variability in expression of the BjSUS paralogs across the different tissues. BjSUS genes showed transcriptional activation in response to defense hormones and a late response to wounding. Tissue and temporal specificity of expression revealed importance of specific SUS paralogs at different developmental stages and under different stress conditions. The study highlighted differential involvement of SUS paralogs in sucrose metabolism across the tissues and stress-responses, in a major oilseed crop B. juncea.

Molecular evolutionary trends and feeding ecology diversification in the Hemiptera, anchored by the milkweed bug genome

BACKGROUND

The Hemiptera (aphids, cicadas, and true bugs) are a key insect order, with high diversity for feeding ecology and excellent experimental tractability for molecular genetics. Building upon recent sequencing of hemipteran pests such as phloem-feeding aphids and blood-feeding bed bugs, we present the genome sequence and comparative analyses centered on the milkweed bug Oncopeltus fasciatus, a seed feeder of the family Lygaeidae.

RESULTS

The 926-Mb Oncopeltus genome is well represented by the current assembly and official gene set. We use our genomic and RNA-seq data not only to characterize the protein-coding gene repertoire and perform isoform-specific RNAi, but also to elucidate patterns of molecular evolution and physiology. We find ongoing, lineage-specific expansion and diversification of repressive C2H2 zinc finger proteins. The discovery of intron gain and turnover specific to the Hemiptera also prompted the evaluation of lineage and genome size as predictors of gene structure evolution. Furthermore, we identify enzymatic gains and losses that correlate with feeding biology, particularly for reductions associated with derived, fluid nutrition feeding.

CONCLUSIONS

With the milkweed bug, we now have a critical mass of sequenced species for a hemimetabolous insect order and close outgroup to the Holometabola, substantially improving the diversity of insect genomics. We thereby define commonalities among the Hemiptera and delve into how hemipteran genomes reflect distinct feeding ecologies. Given Oncopeltus's strength as an experimental model, these new sequence resources bolster the foundation for molecular research and highlight technical considerations for the analysis of medium-sized invertebrate genomes.

Molecular evolutionary trends and feeding ecology diversification in the Hemiptera, anchored by the milkweed bug genome

BACKGROUND

The Hemiptera (aphids, cicadas, and true bugs) are a key insect order, with high diversity for feeding ecology and excellent experimental tractability for molecular genetics. Building upon recent sequencing of hemipteran pests such as phloem-feeding aphids and blood-feeding bed bugs, we present the genome sequence and comparative analyses centered on the milkweed bug Oncopeltus fasciatus, a seed feeder of the family Lygaeidae.

RESULTS

The 926-Mb Oncopeltus genome is well represented by the current assembly and official gene set. We use our genomic and RNA-seq data not only to characterize the protein-coding gene repertoire and perform isoform-specific RNAi, but also to elucidate patterns of molecular evolution and physiology. We find ongoing, lineage-specific expansion and diversification of repressive C2H2 zinc finger proteins. The discovery of intron gain and turnover specific to the Hemiptera also prompted the evaluation of lineage and genome size as predictors of gene structure evolution. Furthermore, we identify enzymatic gains and losses that correlate with feeding biology, particularly for reductions associated with derived, fluid nutrition feeding.

CONCLUSIONS

With the milkweed bug, we now have a critical mass of sequenced species for a hemimetabolous insect order and close outgroup to the Holometabola, substantially improving the diversity of insect genomics. We thereby define commonalities among the Hemiptera and delve into how hemipteran genomes reflect distinct feeding ecologies. Given Oncopeltus's strength as an experimental model, these new sequence resources bolster the foundation for molecular research and highlight technical considerations for the analysis of medium-sized invertebrate genomes.

Conserved Critical Evolutionary Gene Structures in Orthologs

Unravelling gene structure requires the identification and understanding of the constraints that are often associated with the evolutionary history and functional domains of genes. We speculated in this manuscript with the possibility of the existence in orthologs of an emergent highly conserved gene structure that might explain their coordinated evolution during speciation events and their parental function. Here, we will address the following issues: (1) is there any conserved hypothetical structure along ortholog gene sequences? (2) If any, are such conserved structures maintained and conserved during speciation events? The data presented show evidences supporting this hypothesis. We have found that, (1) most orthologs studied share highly conserved compositional structures not observed previously. (2) While the percent identity of nucleotide sequences of orthologs correlates with the percent identity of composon sequences, the number of emergent compositional structures conserved during speciation does not correlate with the percent identity. (3) A broad range of species conserves the emergent compositional stretches. We will also discuss the concept of critical gene structure.