Comparative analysis of intronless genes in teleost fish genomes: insights into their evolution and molecular function

A novel liver-specific immunoglobulin heavy chain-like gene in a cartilaginous fish

We identified a novel immunoglobulin (Ig) heavy chain-like gene (tsIgH) expressed in the liver of the banded houndshark Triakis scyllium by preliminary transcriptomic analysis. The tsIgH gene showed less than 30% of amino acid identities to Ig genes of the shark. The gene encodes one variable domain (VH) and three conserved domains (CH1-CH3) with a predicted signal peptide. Interestingly, this protein has only one cysteine residue in a linker region between VH and CH1 other than those required for the formation of the immunoglobulin domain. Genome sequencing revealed that each of the domains was encoded by a corresponding single exon, and the exon-intron structures of the homologues are conserved in the other cartilaginous fishes. By RT-qPCR analysis, the transcript of the tsIgH gene was observed only in the liver, while that of the IgM was mainly detected in the epigonal organ, liver, and spleen. The novel Ig-heavy chain-like gene in cartilaginous fish may provide new clues to the evolution of immunoglobulin genes.

Genome wide identification and characterization of the amino acid transporter (AAT) genes regulating seed protein content in chickpea (Cicer arietinum L.)

Amino acid transporters (AATs), besides, being a crucial component for nutrient partitioning system are also vital for growth and development of the plants and stress resilience. In order to understand the role of AAT genes in seed quality proteins, a comprehensive analysis of AAT gene family was carried out in chickpea leading to identification of 109 AAT genes, representing 10 subfamilies with random distribution across the chickpea genome. Several important stress responsive cis-regulatory elements like Myb, ABRE, ERE were detected in the promoter region of these CaAAT genes. Most of the genes belonging to the same sub-families shared the intron-exon distribution pattern owing to their conserved nature. Random distribution of these CaAAT genes was observed on plasma membrane, vacuolar membrane, Endoplasmic reticulum and Golgi membranes, which may be associated to distinct biochemical pathways. In total 92 out 109 CaAAT genes arise as result of duplication, among which segmental duplication was more prominent over tandem duplication. As expected, the phylogenetic tree was divided into 2 major clades, and further sub-divided into different sub-families. Among the 109 CaAAT genes, 25 were found to be interacting with 25 miRNAs, many miRNAs like miR156, miR159 and miR164 were interacting only with single AAT genes. Tissues specific expression pattern of many CaAAT genes was observed like CaAAP7 and CaAVT18 in nodules, CaAAP17, CaAVT5 and CaCAT9 in vegetative tissues while CaCAT10 and CaAAP23 in seed related tissues as per the expression analysis. Mature seed transcriptome data revealed that genotypes having high protein content (ICC 8397, ICC 13461) showed low CaAATs expression as compared to the genotypes having low protein content (FG 212, BG 3054). Amino acid profiling of these genotypes revealed a significant difference in amount of essential and non-essential amino acids, probably due to differential expression of CaAATs. Thus, the present study provides insights into the biological role of AAT genes in chickpea, which will facilitate their functional characterization and role in various developmental stages, stress responses and involvement in nutritional quality enhancement.

Evolutionary Perspective and Expression Analysis of Intronless Genes Highlight the Conservation of Their Regulatory Role

The structure of eukaryotic genes is generally a combination of exons interrupted by intragenic non-coding DNA regions (introns) removed by RNA splicing to generate the mature mRNA. A fraction of genes, however, comprise a single coding exon with introns in their untranslated regions or are intronless genes (IGs), lacking introns entirely. The latter code for essential proteins involved in development, growth, and cell proliferation and their expression has been proposed to be highly specialized for neuro-specific functions and linked to cancer, neuropathies, and developmental disorders. The abundant presence of introns in eukaryotic genomes is pivotal for the precise control of gene expression. Notwithstanding, IGs exempting splicing events entail a higher transcriptional fidelity, making them even more valuable for regulatory roles. This work aimed to infer the functional role and evolutionary history of IGs centered on the mouse genome. IGs consist of a subgroup of genes with one exon including coding genes, non-coding genes, and pseudogenes, which conform approximately 6% of a total of 21,527 genes. To understand their prevalence, biological relevance, and evolution, we identified and studied 1,116 IG functional proteins validating their differential expression in transcriptomic data of embryonic mouse telencephalon. Our results showed that overall expression levels of IGs are lower than those of MEGs. However, strongly up-regulated IGs include transcription factors (TFs) such as the class 3 of POU (HMG Box), Neurog1, Olig1, and BHLHe22, BHLHe23, among other essential genes including the β-cluster of protocadherins. Most striking was the finding that IG-encoded BHLH TFs fit the criteria to be classified as microproteins. Finally, predicted protein orthologs in other six genomes confirmed high conservation of IGs associated with regulating neural processes and with chromatin organization and epigenetic regulation in Vertebrata. Moreover, this study highlights that IGs are essential modulators of regulatory processes, such as the Wnt signaling pathway and biological processes as pivotal as sensory organ developing at a transcriptional and post-translational level. Overall, our results suggest that IG proteins have specialized, prevalent, and unique biological roles and that functional divergence between IGs and MEGs is likely to be the result of specific evolutionary constraints.

Genome-wide analysis of European sea bass provides insights into the evolution and functions of single-exon genes

Several studies have attempted to understand the origin and evolution of single-exon genes (SEGs) in eukaryotic organisms, including fishes, but few have examined the functional and evolutionary relationships between SEGs and multiple-exon gene (MEG) paralogs, in particular the conservation of promoter regions. Given that SEGs originate via the reverse transcription of mRNA from a "parental" MEGs, such comparisons may enable identifying evolutionarily-related SEG/MEG paralogs, which might fulfill equivalent physiological functions. Here, the relationship of SEG proportion with MEG count, gene density, intron count, and chromosome size was assessed for the genome of the European sea bass, Dicentrarchus labrax. Then, SEGs with an MEG parent were identified, and promoter sequences of SEG/MEG paralogs were compared, to identify highly conserved functional motifs. The results revealed a total count of 1,585 (8.3% of total genes) SEGs in the European sea bass genome, which was correlated with MEG count but not with gene density. The significant correlation of SEG content with the number of MEGs suggests that SEGs were continuously and independently generated over evolutionary time following species divergence through retrotranscription events, followed by tandem duplications. Functional annotation showed that the majority of SEGs are functional, as is evident from their expression in RNA-seq data used to support homology-based genome annotation. Differences in 5'UTR and 3'UTR lengths between SEG/MEG paralogs observed in this study may contribute to gene expression divergence between them and therefore lead to the emergence of new SEG functions. The comparison of nonsynonymous to synonymous changes (Ka/Ks) between SEG/MEG parents showed that 74 of them are under positive selection (Ka/Ks > 1; p = .0447). An additional fifteen SEGs with an MEG parent have a common promoter, which implies that they are under the influence of common regulatory networks.

Identification and characterization of MADS box gene family in pigeonpea for their role during floral transition

MADS box genes are class of transcription factors involved in various physiological and developmental processes in plants. To understand their role in floral transition-related pathways, a genome-wide identification was done in Cajanus cajan, identifying 102 members which were classified into two different groups based on their gene structure. The status of all these genes was further analyzed in three wild species i.e. C. scarabaeoides, C. platycarpus and C. cajanifolius which revealed absence of 31-34 MADS box genes in them hinting towards their role in domestication and evolution. We could locate only a single copy of both FLOWERING LOCUS C (FLC) and SHORT VEGETATIVE PHASE (SVP) genes, while three paralogs of SUPPRESSOR OF ACTIVATION OF CONSTANS 1 (SOC1) were found in C. cajan genome. One of those SOC1 paralogs i.e. CcMADS1.5 was found to be missing in all three wild relatives, also forming separate clade in phylogeny. This SOC1 gene was also lacking the characteristic MADS box domain in it. Expression profiling of major MADS box genes involved in flowering was done in different tissues viz shoot apical meristem, vegetative leaf, reproductive meristem, and reproductive bud. Gene-based time tree of FLC and SOC1 gene dictates their divergence from Arabidopsis before 71 and 23 million year ago (mya), respectively. This study provides valuable insights into the functional characteristics, expression pattern, and evolution of MADS box proteins in grain legumes with emphasis on C. cajan, which may help in further characterizing these genes.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s13205-020-02605-7.

RIGD: A Database for Intronless Genes in the Rosaceae

Most eukaryotic genes are interrupted by one or more introns, and only prokaryotic genomes are composed of mainly single-exon genes without introns. Due to the absence of introns, intronless genes in eukaryotes have become important materials for comparative genomics and evolutionary biology. There is currently no cohesive database that collects intronless genes in plants into a single database, although many databases on exons and introns exist. In this study, we constructed the Rosaceae Intronless Genes Database (RIGD), a user-friendly web interface to explore and collect information on intronless genes from different plants. Six Rosaceae species, Pyrus bretschneideri, Pyrus communis, Malus domestica, Prunus persica, Prunus mume, and Fragaria vesca, are included in the current release of the RIGD. Sequence data and gene annotation were collected from different databases and integrated. The main purpose of this study is to provide gene sequence data. In addition, attribute analysis, functional annotations, subcellular localization prediction, and GO analysis are reported. The RIGD allows users to browse, search, and download data with ease. Blast and comparative analyses are also provided through this online database, which is available at http://www.rigdb.cn/.

Polymorphism of MHC class IIB in an acheilognathid species, Rhodeus sinensis shaped by historical selection and recombination

BACKGROUND

Rhodeus sinensis is a bitterling species occurring throughout the numerous freshwater systems on the East Asia. Here, we analyzed the diversity of the MHC class IIB (DAB) genes from this species, which may offer meaningful insights into evolutionary processes in this species as well as other bitterlings.

RESULTS

Using cDNA and gDNA samples from 50 individuals, we discovered classical 140 allelic sequences that could be allocated into either DAB1 (Rhsi-DAB1) or DAB3 (Rhsi-DAB3). DAB sequences completely lacking the intron, but identical or similar to Rhsi-DAB1, were also discovered from our gDNA samples, and this intron loss likely originated from the retrotransposition events of processed mDNA. The β1 domain was the most polymorphic in both Rhsi-DAB1 and -DAB3. Putative peptide biding residues (PBRs) in Rhsi-DAB1, but not in Rhsi-DAB3, exhibited a significant dN/dS, presumably indicating that different selection pressures have acted on those two DABs. Recombination between different alleles seemed to have contributed to the increase of diversity in Rhsi-DABs. Upon phylogenetic analysis, Rhsi-DAB1 and -DAB3 formed independent clusters. Several alleles from other species of Cypriniformes were embedded in the clade of Rhsi-DAB1, whereas Rhsi-DAB3 clustered with alleles from the wider range of taxa (Cyprinodontiformes), indicating that these two Rhsi-DABs have taken different historical paths.

CONCLUSIONS

A great deal of MHC class IIB allelic diversity was found in R. sinensis, and gene duplication, selection and recombination may have contributed to this diversity. Based on our data, it is presumed that such historical processes have commonly or differently acted on the polymorphism of Rhsi-DAB1 and -DAB3.

Genome-wide identification and characterization of toll-like receptor genes in spotted sea bass (Lateolabrax maculatus) and their involvement in the host immune response to Vibrio harveyi infection

Toll-like receptor (TLR) genes are the earliest reported pathogen recognition receptors (PRRs) and have been extensively studied. These genes play pivotal roles in the innate immune defense against pathogen invasion. In this study, a total of 16 tlr genes were identified and characterized in spotted sea bass (Lateolabrax maculatus). The tlr genes of spotted sea bass were classified into five subfamilies (tlr1-subfamily, tlr3-subfamily, tlr5-subfamily, tlr7-subfamily, and tlr11-subfamily) according to the phylogenetic analysis, and their annotations were confirmed by a syntenic analysis. The protein domain analysis indicated that most tlr genes had the following three major TLR protein domains: a leucine-rich repeat (LRR) domain, a transmembrane region (TM) and a Toll/interleukin-1 receptor (TIR) domain. The tlr genes in spotted sea bass were distributed in 11 of 24 chromosomes. The mRNA expression levels of 16 tlr genes in response to Vibrio harveyi infection were quantified in the head kidney. Most genes were downregulated following V. harveyi infection, while only 5 tlr genes, including tlr1-1, tlr1-2, tlr2-2, tlr5, and tlr7, were significantly upregulated. Collectively, these results help elucidate the crucial roles of tlr genes in the immune response of spotted sea bass and may supply valuable genomic resources for future studies investigating fish disease management.

Improved ontology for eukaryotic single-exon coding sequences in biological databases

Efficient extraction of knowledge from biological data requires the development of structured vocabularies to unambiguously define biological terms. This paper proposes descriptions and definitions to disambiguate the term 'single-exon gene'. Eukaryotic Single-Exon Genes (SEGs) have been defined as genes that do not have introns in their protein coding sequences. They have been studied not only to determine their origin and evolution but also because their expression has been linked to several types of human cancer and neurological/developmental disorders and many exhibit tissue-specific transcription. Unfortunately, the term 'SEGs' is rife with ambiguity, leading to biological misinterpretations. In the classic definition, no distinction is made between SEGs that harbor introns in their untranslated regions (UTRs) versus those without. This distinction is important to make because the presence of introns in UTRs affects transcriptional regulation and post-transcriptional processing of the mRNA. In addition, recent whole-transcriptome shotgun sequencing has led to the discovery of many examples of single-exon mRNAs that arise from alternative splicing of multi-exon genes, these single-exon isoforms are being confused with SEGs despite their clearly different origin. The increasing expansion of RNA-seq datasets makes it imperative to distinguish the different SEG types before annotation errors become indelibly propagated in biological databases. This paper develops a structured vocabulary for their disambiguation, allowing a major reassessment of their evolutionary trajectories, regulation, RNA processing and transport, and provides the opportunity to improve the detection of gene associations with disorders including cancers, neurological and developmental diseases.

An integrated and comparative approach towards identification, characterization and functional annotation of candidate genes for drought tolerance in sorghum (Sorghum bicolor (L.) Moench)

BACKGROUND

Drought is the most disastrous abiotic stress that severely affects agricultural productivity worldwide. Understanding the biological basis of drought-regulated traits, requires identification and an in-depth characterization of genetic determinants using model organisms and high-throughput technologies. However, studies on drought tolerance have generally been limited to traditional candidate gene approach that targets only a single gene in a pathway that is related to a trait. In this study, we used sorghum, one of the model crops that is well adapted to arid regions, to mine genes and define determinants for drought tolerance using drought expression libraries and RNA-seq data.

RESULTS

We provide an integrated and comparative in silico candidate gene identification, characterization and annotation approach, with an emphasis on genes playing a prominent role in conferring drought tolerance in sorghum. A total of 470 non-redundant functionally annotated drought responsive genes (DRGs) were identified using experimental data from drought responses by employing pairwise sequence similarity searches, pathway and interpro-domain analysis, expression profiling and orthology relation. Comparison of the genomic locations between these genes and sorghum quantitative trait loci (QTLs) showed that 40% of these genes were co-localized with QTLs known for drought tolerance. The genome reannotation conducted using the Program to Assemble Spliced Alignment (PASA), resulted in 9.6% of existing single gene models being updated. In addition, 210 putative novel genes were identified using AUGUSTUS and PASA based analysis on expression dataset. Among these, 50% were single exonic, 69.5% represented drought responsive and 5.7% were complete gene structure models. Analysis of biochemical metabolism revealed 14 metabolic pathways that are related to drought tolerance and also had a strong biological network, among categories of genes involved. Identification of these pathways, signifies the interplay of biochemical reactions that make up the metabolic network, constituting fundamental interface for sorghum defence mechanism against drought stress.

CONCLUSIONS

This study suggests untapped natural variability in sorghum that could be used for developing drought tolerance. The data presented here, may be regarded as an initial reference point in functional and comparative genomics in the Gramineae family.

Identification and characteristic of three members of the C1q/TNF-related proteins (CTRPs) superfamily in Eudontomyzon morii

C1q is the target recognition protein of the classical complement pathway and a major connecting link between innate and acquired immunity. C1q and the multifunctional tumor necrosis factor (TNF) ligand family is of similar crystal structures, are designated the C1q/TNF-related proteins (CTRPs) superfamily. They are involved in processes as diverse as host defense, inflammation, apoptosis, autoimmunity, cell differentiation, organogenesis, hibernation and insulinresistant obesity. In this study, three members of the CTRPs superfamily were isolated and characterized in Yalu River lampreys (Eudontomyzon morii), and are respectively named LaC1qC, LaCTRP1, LaCTRP9. The full-length cDNAs of C1qC-like (LaC1qAL), CTRP1-like (LaCTRP1), and CTRP9-like (LaCTRP9) consist of 723, 762 and 825 bp of nucleotide sequence encoding polypeptides of 241, 254 and 275 amino acids, respectively. All-three proteins share three common domains: a signal peptide at the N terminus, a collagenous domain (characteristic Gly-X-Y repeats), and a C-terminal globular domain. In addition, the higher expression level of the three proteins in heart by RT-PCR and real-time PCR tissue profiling implied that they might involve in immune response or injury repair of the heart in lamprey.

IGDD: a database of intronless genes in dicots

Background

Intronless genes are a significant characteristic of prokaryotes. Systematic identification and annotation are primary and crucial steps for determining the functions of intronless genes and understanding their occurrence in eukaryotes.

Description

In this paper, we describe the construction of the Intronless Genes Database in Dicots (IGDD; available at http://bio.njfu.edu.cn/igdd/), which contains data for five well-annotated plants including Arabidopsis thaliana, Carica papaya, Populus trichocarpa, Salix suchowensis and Vitis vinifera. Using highly visual settings, IGDD displays the structural and functional annotations, the homolog groups, the syntenic relationships, the expression patterns, and the statistical characteristics of intronless genes. In addition, useful tools such as an advanced search and local BLAST are available through a user-friendly and intuitive web interface.

Conclusion

In conclusion, the IGDD provides a comprehensive and up-to-date platform for researchers to assist the exploration of intronless genes in dicot plants.

Molecular cloning, structure and expressional profiles of two novel single-exon genes (PoCCR6A and PoCCR6B) in the Japanese flounder (Paralichthys olivaceus)

CCR6 is an important binding receptor of CCL20 and beta-defensins, and has multiple functions in the innate and acquired immune responses. In this study, we cloned the PoCCR6A and PoCCR6B genes of the Japanese flounder and studied the gene structure and expression patterns of these two genes in bacterial infection. The full-length PoCCR6A cDNA is 1415 bp and the open reading frame (ORF) is 1113 bp, encoding a 370-amino-acid peptide. The full-length PoCCR6B cDNA is 2193 bp and the ORF is 1029 bp, encoding a 363-amino-acid peptide. The structures of PoCCR6A and PoCCR6B indicate that they are single-exon genes. The predicted proteins encoded by PoCCR6A and PoCCR6B have the typical G-protein-coupled receptor (GPCR) family signature of seven transmembrane domains and several conserved structural features. A tissue distribution analysis showed that PoCCR6A is predominately expressed in the intestine, gill, and blood, and PoCCR6B in the gill, spleen, and liver. The expression patterns of the two chemokine receptors were analyzed during bacterial infection. In spleen and kidney, the expression of PoCCR6A was significantly upregulated at 24 h after infection, whereas the expression of PoCCR6B was steady at these time points. While in intestine, both of them were upregulated at 6 h-12 h after infection, and in gill the expression levels of them were upregulated at 24 h. The patterns of expression suggested that PoCCR6A and PoCCR6B play an important role in the immune response of the Japanese flounder, especially in the mucosal tissues.

PIGD: a database for intronless genes in the Poaceae

Background

Intronless genes are a feature of prokaryotes; however, they are widespread and unequally distributed among eukaryotes and represent an important resource to study the evolution of gene architecture. Although many databases on exons and introns exist, there is currently no cohesive database that collects intronless genes in plants into a single database.

Description

In this study, we present the Poaceae Intronless Genes Database (PIGD), a user-friendly web interface to explore information on intronless genes from different plants. Five Poaceae species, Sorghum bicolor, Zea mays, Setaria italica, Panicum virgatum and Brachypodium distachyon, are included in the current release of PIGD. Gene annotations and sequence data were collected and integrated from different databases. The primary focus of this study was to provide gene descriptions and gene product records. In addition, functional annotations, subcellular localization prediction and taxonomic distribution are reported. PIGD allows users to readily browse, search and download data. BLAST and comparative analyses are also provided through this online database, which is available at http://pigd.ahau.edu.cn/.

Conclusion

PIGD provides a solid platform for the collection, integration and analysis of intronless genes in the Poaceae. As such, this database will be useful for subsequent bio-computational analysis in comparative genomics and evolutionary studies.

Involvement of claudins in zebrafish brain ventricle morphogenesis

Zebrafish brain ventricle morphogenesis involves an initial circulation-independent opening followed by a blood flow- and circulation-dependent expansion process. Zebrafish claudin-5a is required for the establishment of a neuroepithelial-ventricular barrier, which maintains the hydrostatic pressure within the ventricular cavity, thereby contributing to brain ventricle opening and expansion. In mammalia, several claudin family members, including claudin-3 and claudin-5, are expressed within microvessel endothelial cells of the blood-brain barrier. Whether zebrafish brain ventricle morphogenesis provides a model for studying these claudins during early embryonic development was unknown. This review focuses on the expression and function of these zebrafish claudins during brain ventricle morphogenesis.