The Evolutionary Relationship between Alternative Splicing and Gene Duplication

A serpin gene from a parasitoid wasp disrupts host immunity and exhibits adaptive alternative splicing

Alternative splicing (AS) is a major source of protein diversity in eukaryotes, but less is known about its evolution compared to gene duplication (GD). How AS and GD interact is also largely understudied. By constructing the evolutionary trajectory of the serpin gene PpSerpin-1 (Pteromalus puparum serpin 1) in parasitoids and other insects, we found that both AS and GD jointly contribute to serpin protein diversity. These two processes are negatively correlated and show divergent features in both protein and regulatory sequences. Parasitoid wasps exhibit higher numbers of serpin protein/domains than nonparasitoids, resulting from more GD but less AS in parasitoids. The potential roles of AS and GD in the evolution of parasitoid host-effector genes are discussed. Furthermore, we find that PpSerpin-1 shows an exon expansion of AS compared to other parasitoids, and that several isoforms are involved in the wasp immune response, have been recruited to both wasp venom and larval saliva, and suppress host immunity. Overall, our study provides an example of how a parasitoid serpin gene adapts to parasitism through AS, and sheds light on the differential features of AS and GD in the evolution of insect serpins and their associations with the parasitic life strategy.

Novel Gene Rearrangements in Mitochondrial Genomes of four families of Praying Mantises (Insecta, Mantodea) and Phylogenetic Relationships of Mantodea

The mitochondrial genome (mitogenome) plays an important role in phylogenetic studies of many species. The mitogenomes of many praying mantis groups have been well-studied, but mitogenomes of special mimic praying mantises, especially Acanthopoidea and Galinthiadoidea species, are still sorely lacking in the NCBI database. The present study analyzes five mitogenomes from four species of Acanthopoidea (Angela sp., Callibia diana, Coptopteryx sp., Raptrix fusca) and one of Galinthiadoidea (Galinthias amoena) that were sequenced by the primer-walking method. Among Angela sp. and Coptopteryx sp., three gene rearrangements were detected in ND3-A-R-N-S-E-F and COX1-L2-COX2 gene regions, two of which were novel. In addition, individual tandem repeats were found in control regions of four mitogenomes (Angela sp., C. diana, Coptopteryx sp., G. amoena). For those, plausible explanations were derived from the tandem duplication-random loss (TDRL) model and the slipped-strand mispairing model. One potential motif was found in Acanthopidae that was seen as a synapomorphy. Several conserved block sequences (CBSs) were detected within Acanthopoidea that paved the way for the design of specific primers. Via BI and ML analysis, based on four datasets (PCG12, PCG12R, PCG123, PCG123R), the merged phylogenetic tree within Mantodea was reconstructed. This showed that the monophyly of Acanthopoidea was supported and that the PCG12R dataset was the most suitable for reconstructing the phylogenetic tree within Mantodea.

Virus usurps alternative splicing to clear the decks for infection

Since invasion, there will be a tug-of-war between host and virus to scramble cellular resources, for either restraining or facilitating infection. Alternative splicing (AS) is a conserved and critical mechanism of processing pre-mRNA into mRNAs to increase protein diversity in eukaryotes. Notably, this kind of post-transcriptional regulatory mechanism has gained appreciation since it is widely involved in virus infection. Here, we highlight the important roles of AS in regulating viral protein expression and how virus in turn hijacks AS to antagonize host immune response. This review will widen the understandings of host-virus interactions, be meaningful to innovatively elucidate viral pathogenesis, and provide novel targets for developing antiviral drugs in the future.

The Role of non-muscle actin paralogs in cell cycle progression and proliferation

Uncontrolled cell proliferation leads to several pathologies, including cancer. Thus, this process must be tightly regulated. The cell cycle accounts for cell proliferation, and its progression is coordinated with changes in cell shape, for which cytoskeleton reorganization is responsible. Rearrangement of the cytoskeleton allows for its participation in the precise division of genetic material and cytokinesis. One of the main cytoskeletal components is filamentous actin-based structures. Mammalian cells have at least six actin paralogs, four of which are muscle-specific, while two, named β- and γ-actin, are abundantly present in all types of cells. This review summarizes the findings that establish the role of non-muscle actin paralogs in regulating cell cycle progression and proliferation. We discuss studies showing that the level of a given non-muscle actin paralog in a cell influences the cell's ability to progress through the cell cycle and, thus, proliferation. Moreover, we elaborate on the non-muscle actins' role in regulating gene transcription, interactions of actin paralogs with proteins involved in controlling cell proliferation, and the contribution of non-muscle actins to different structures in a dividing cell. The data cited in this review show that non-muscle actins regulate the cell cycle and proliferation through varying mechanisms. We point to the need for further studies addressing these mechanisms.

Genome-wide analysis of bromodomain gene family in Arabidopsis and rice

The bromodomain-containing proteins (BRD-proteins) belongs to family of 'epigenetic mark readers', integral to epigenetic regulation. The BRD-members contain a conserved 'bromodomain' (BRD/BRD-fold: interacts with acetylated-lysine in histones), and several additional domains, making them structurally/functionally diverse. Like animals, plants also contain multiple Brd-homologs, however the extent of their diversity and impact of molecular events (genomic duplications, alternative splicing, AS) therein, is relatively less explored. The present genome-wide analysis of Brd-gene families of Arabidopsis thaliana and Oryza sativa showed extensive diversity in structure of genes/proteins, regulatory elements, expression pattern, domains/motifs, and the bromodomain (w.r.t. length, sequence, location) among the Brd-members. Orthology analysis identified thirteen ortholog groups (OGs), three paralog groups (PGs) and four singleton members (STs). While more than 40% Brd-genes were affected by genomic duplication events in both plants, AS-events affected 60% A. thaliana and 41% O. sativa genes. These molecular events affected various regions (promoters, untranslated regions, exons) of different Brd-members with potential impact on expression and/or structure-function characteristics. RNA-Seq data analysis indicated differences in tissue-specificity and stress response of Brd-members. Analysis by RT-qPCR revealed differential abundance and salt stress response of duplicate A. thaliana and O. sativa Brd-genes. Further analysis of AtBrd gene, AtBrdPG1b showed salinity-induced modulation of splicing pattern. Bromodomain (BRD)-region based phylogenetic analysis placed the A. thaliana and O. sativa homologs into clusters/sub-clusters, mostly consistent with ortholog/paralog groups. The bromodomain-region displayed several conserved signatures in key BRD-fold elements (α-helices, loops), along with variations (1-20 sites) and indels among the BRD-duplicates. Homology modeling and superposition identified structural variations in BRD-folds of divergent and duplicate BRD-members, which might affect their interaction with the chromatin histones, and associated functions. The study also showed contribution of various duplication events in Brd-gene family expansion among diverse plants, including several monocot and dicot plant species.

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

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

A positive feedback loop between two C-type lectins originated from gene duplication and relish promotes the expression of antimicrobial peptides in Procambarus clarkii

Gene duplication (GD) leads to the expansion of gene families that contributes organisms adapting to stress or environment and dealing with the infection of various pathogens. C-type lectins (CTLs) in crustaceans undergo gene expansion and participate in various immune responses. However, the functions of different CTL produced by GD are not fully characterized. In the present study, two CTL genes (designated as PcLec-EPS and PcLec-QPS, respectively) were identified from Procambarus clarkii. PcLec-EPS and PcLec-QPS originate from GD and the main difference between them is exon 3. PcLec-EPS and PcLec-QPS respectively contains EPS and QPS motif in their carbohydrate recognition domain. The mRNA levels of PcLec-EPS and PcLec-QPS in hemocytes, gills, intestine and lymph underwent time-dependent enhancement after D-Mannose and D-Galactose challenge. Recombinant PcLec-EPS and PcLec-QPS could bind to carbohydrates and microbes, and agglutinate bacteria. The results of experiments on recombinant protein injection and RNA interference indicate that PcLec-EPS and PcLec-QPS can respectively strong recognize and bind D-Mannose and D-Galactose, activate the Relish transcriptional factor, and further upregulate the expression of different antimicrobial peptides (AMPs). In addition, these two CTLs and Relish could positively regulate the expression of each other, suggesting that there is a positive feedback loop between two CTLs and Relish that regulates the expression of AMPs. It may contribute to the expansion of the immune response for host quickly and efficiently eliminating pathogenic microorganisms. This study provides new knowledge for clear understanding the significance and function of different CTL generated by GD in immune defenses in crustacean.

Alteration of synonymous codon usage bias accompanies polyploidization in wheat

The diploidization of polyploid genomes is accompanied by genomic variation, including synonymous nucleotide substitutions that may lead to synonymous codon usage bias (SCUB). SCUB can mirror the evolutionary specialization of plants, but its effect on the formation of polyploidies is not well documented. We explored this issue here with hexaploid wheat and its progenitors. Synonymous codons (SCs) ending in either cytosine (NNC) or guanidine (NNG) were more frequent than those ending in either adenosine (NNA) or thymine (NNT), and the preference for NNC/G codons followed the increase in genome ploidy. The ratios between NNC/G and NNA/T codons gradually decreased in genes with more introns, and the difference in these ratios between wheat and its progenitors diminished with increasing ploidy. SCUB frequencies were heterogeneous among exons, and the bias preferred to NNA/T in more internal exons, especially for genes with more exons; while the preference did not appear to associate with ploidy. The SCUB alteration of the progenitors was different during the formation of hexaploid wheat, so that SCUB was the homogeneous among A, B and D subgenomes. DNA methylation-mediated conversion from cytosine to thymine weakened following the increase of genome ploidy, coinciding with the stronger bias for NNC/G SCs in the genome as a function of ploidy, suggesting that SCUB contribute to the epigenetic variation in hexaploid wheat. The patterns in SCUB mirrored the formation of hexaploid wheat, which provides new insight into genome shock-induced genetic variation during polyploidization. SCs representing non-neutral synonymous mutations can be used for genetic dissection and improvement of agricultural traits of wheat and other polyploidies.

Characterization of aquaporin-1ab (Aqp1ab) mRNA in mud loach (Misgurnus mizolepis) exposed to heavy metal and immunostimulant stimuli

Aquaporins (AQPs) are key proteins that regulate fluid homeostasis in cells via modulating osmotic water transport. In the present study, we identified three variants of Aqp1ab transcript (mmAQP1ab x1, mmAQP1ab x2, and mmAQP1ab x3) in mud loaches (Misgurnus mizolepis), and their expression patterns were examined in response to heavy metal and immunostimulant exposure. Mud loach Aqp1ab gene has a somewhat different organizational structure (i.e. five exons interrupted by four introns) compared to most other teleostean Aqp1ab orthologues, which have four exons. The 5'-flanking regulatory region of Aqp gene showed diverse transcription factor binding motifs, particularly those associated with stress/immune responses. Developmental expression patterns indicated that Aqp1ab mRNA was maternally inherited, presumably important for fine-tuning gene expression during embryonic and early larval developments. Expression of mud loach Aqp1ab mRNA was significantly and differentially modulated in several tissues (intestine, kidneys, spleen, and liver) in response to various heavy metal treatments. In addition, Aqp1ab gene expression was highly induced in response to immune challenge (LPS and polyI:C injections). Collectively, our results suggested that AQPs are multifunctional effectors playing diverse roles in cellular pathways relevant to immune and/or stress adaptation responses, in addition to their involvement in osmoregulation.

Viral Encoded miRNAs in Tumorigenesis: Theranostic Opportunities in Precision Oncology

About 15% of all human cancers have a viral etiology. Although progress has been made, understanding the viral oncogenesis and associated molecular mechanisms remain complex. The discovery of cellular miRNAs has led to major breakthroughs. Interestingly, viruses have also been discovered to encode their own miRNAs. These viral, small, non-coding miRNAs are also known as viral-miRNAs (v-miRNAs). Although the function of v-miRNAs largely remains to be elucidated, their role in tumorigenesis cannot be ignored. V-miRNAs have also been shown to exploit the cellular machinery to benefit viral replication and survival. Although the discovery of Hepatitis C virus (HCV), and its viral miRNAs, is a work in progress, the existence of HPV-, EBV-, HBV-, MCPyV- and KSHV-encoded miRNA has been documented. V-miRNAs have been shown to target host factors to advance tumorigenesis, evade and suppress the immune system, and deregulate both the cell cycle and the apoptotic machinery. Although the exact mechanisms of v-miRNAs-induced tumorigenesis are still unclear, v-miRNAs are active role-players in tumorigenesis, viral latency and cell transformation. Furthermore, v-miRNAs can function as posttranscriptional gene regulators of both viral and host genes. Thus, it has been proposed that v-miRNAs may serve as diagnostic biomarkers and therapeutic targets for cancers with a viral etiology. Although significant challenges exist in their clinical application, emerging reports demonstrate their potent role in precision medicine. This review will focus on the roles of HPV-, HCV-, EBV-, HBV-, MCPyV-, and KSHV-produced v-miRNAs in tumorigenesis, as effectors in immune evasion, as diagnostic biomarkers and as novel anti-cancer therapeutic targets. Finally, it will discuss the challenges and opportunities associated with v-miRNAs theranostics in precision oncology.

Novel Mitochondrial Gene Rearrangement and Intergenic Regions Exist in the Mitochondrial Genomes from Four Newly Established Families of Praying Mantises (Insecta: Mantodea)

Simple Summary

Mantodea is regarded as an excellent material to study the gene rearrangements and large non-coding regions (LNCRs) in mitochondrial genomes. Meanwhile, as a result of the convergent evolution and parallelism, the gene rearrangements and LNCRs are specific to some taxonomic groups within Mantodea, which play an important role in phylogenetic relationship research. Nine mitochondrial genomes (mitogenomes) from four newly established families of praying mantises are obtained and annotated. Eight types of gene rearrangements, including four novel types of gene rearrangements in Mantodea, are detected, which can be explained by the tandem replication-random loss (TDRL) model. Moreover, one conserved motif between trnI-trnQ is detected in Toxoderidae. This study shed light on the formation mechanisms of these gene rearrangements and LNCRs in four newly established families of praying mantises.

Abstract

Long non-coding regions (NCRs) and gene rearrangements are commonly seen in mitochondrial genomes of Mantodea and are primarily focused on three regions: CR-I-Q-M-ND2, COX2-K-D-ATP8, and ND3-A-R-N-S-E-F-ND5. In this study, eight complete and one nearly complete mitochondrial genomes of praying mantises were acquired for the purpose of discussing mitochondrial gene rearrangements and phylogenetic relationships within Mantodea, primarily in the newly established families Haaniidae and Gonypetidae. Except for Heterochaeta sp. JZ-2017, novel mitochondrial gene arrangements were detected in Cheddikulama straminea, Sinomiopteryx graham, Pseudovates chlorophaea, Spilomantis occipitalis. Of note is the fact that one type of novel arrangement was detected for the first time in the Cyt b-S2-ND1 region. This could be reliably explained by the tandem replication-random loss (TDRL) model. The long NCR between trnT and trnP was generally found in Iridopteryginae and was similar to the ND4L or ND6 gene. Combined with gene rearrangements and intergenic regions, the monophyly of Haaniidae was supported, whereas the paraphyly of Gonypetidae was recovered. Furthermore, several synapomorphies unique to some clades were detected that conserved block sequences between trnI and trnQ and gaps between trnT and trnP in Toxoderidae and Iridopteryginae, respectively.

Co-Expression Analysis Reveals Differential Expression of Homologous Genes Associated with Specific Terpenoid Biosynthesis in Rehmannia glutinosa

Terpenoids are naturally occurring compounds involved in respiration, photosynthesis, membrane fluidity, and pathogen interactions and are classified according to the structure of their carbon skeleton. Although most terpenoids possess pharmacological activity, knowledge about terpenoid metabolism in medicinal plants is insufficient. Rehmannia glutinosa (R. glutinosa) is a traditional herb that is widely used in East Asia and has been reported to contain various terpenoids. In this study, we performed a comprehensive transcriptome analysis of terpenoid metabolism in R. glutinosa using two RNA sequencing platforms: Illumina and PacBio. The results show that the sterol, saponin, iridoid, and carotenoid pathways are active in R. glutinosa. Sterol and saponin biosynthesis were mevalonate pathway dependent, whereas iridoid and carotenoid biosynthesis were methylerythritol 4-phosphate pathway dependent. In addition, we found that the homologous genes of key enzymes involved in terpenoid metabolism were expressed differentially and that the differential expression of these genes was associated with specific terpenoid biosynthesis. The different expression of homologous genes encoding acetyl-CoA acetyltransferase, 3-hydroxy-3-methylglutaryl-CoA reductase, mevalonate kinase, mevalonate diphosphate decarboxylase, farnesyl pyrophosphate synthase, squalene synthase, and squalene epoxidase was associated with sterol and saponin biosynthesis. Homologous genes encoding 1-deoxy-D-xylulose 5-phosphate synthase were also differentially expressed and were associated with carotenoid and iridoid biosynthesis. These results suggest that the biosynthesis of specific terpenoids can be regulated by the homologous of key enzymes involved in plant terpenoid metabolism.

ATTED-II v11: A Plant Gene Coexpression Database Using a Sample Balancing Technique by Subagging of Principal Components

ATTED-II (https://atted.jp) is a gene coexpression database for nine plant species based on publicly available RNAseq and microarray data. One of the challenges in constructing condition-independent coexpression data based on publicly available gene expression data is managing the inherent sampling bias. Here, we report ATTED-II version 11, wherein we adopted a coexpression calculation methodology to balance the samples using principal component analysis and ensemble calculation. This approach has two advantages. First, omitting principal components with low contribution rates reduces the main contributors of noise. Second, balancing large differences in contribution rates enables considering various sample conditions entirely. In addition, based on RNAseq- and microarray-based coexpression data, we provide species-representative, integrated coexpression information to enhance the efficiency of interspecies comparison of the coexpression data. These coexpression data are provided as a standardized z-score to facilitate integrated analysis with different data sources. We believe that with these improvements, ATTED-II is more valuable and powerful for supporting interspecies comparative studies and integrated analyses using heterogeneous data.

Correlated Expression of the Opsin Retrogene LWS-R and its Host Gene in Two Poeciliid Fishes

The important role of retrogenes in genome evolution and species differentiation is becoming increasingly accepted. One synapomorphy among cyprinodontoid fish is a retrotransposed version of a long-wavelength sensitive (LWS) opsin gene, LWS-R, within an intron of the gephyrin (GPHN) gene. These two genes display opposing orientations. It had been speculated that LWS-R hijacks the cis-regulatory elements of GPHN for transcription, but whether their expression is correlated had remained unclear. Here, in silico predictions identified putative promoters upstream of the translation start site of LWS-R, indicating that its transcription is driven by its own promoter rather than by the GPHN promoter. However, consistent expression ratios of LWS-R:GPHN in the eyeball and brain of fishes indicate that the respective gene transcriptions are correlated. Co-expression is potentially modulated by histone exchange during GPHN transcription. Two isoforms were detected in this study, i.e., intron-free and intron-retaining. Intron-free LWS-R was only expressed in the eyeball of fishes, whereas intron-retaining LWS-R occurred in both eyeball and brain. Expression of vision-associated LWS-R beyond the eyeball supports that it is co-expressed with more ubiquitous GPHN.

Expression Analysis of the TdDRF1 Gene in Field-Grown Durum Wheat under Full and Reduced Irrigation

Some of the key genes and regulatory mechanisms controlling drought response in durum wheat have been identified. One of the major challenges for breeders is how to use this knowledge for the achievement of drought stress tolerance. In the present study, we report the expression profiles of the TdDRF1 gene, at consecutive plant growth stages, from different durum wheat genotypes evaluated in two different field environments. The expression of a possible target gene (Wdnh13) of the TdDRF1 gene was also investigated and analogies with the transcript profiles were found. The results of the qRT-PCR highlighted differences in molecular patterns, thus suggesting a genotype dependency of the TdDRF1 gene expression in response to the stress induced. Furthermore, a statistical association between the expression of TdDRF1 transcripts and agronomic traits was also performed and significant differences were found among genotypes, suggesting a relationship. One of the genotypes was found to combine molecular and agronomic characteristics.

The role of alternative splicing in adaptation and evolution

Regulation of gene expression plays a central role in adaptive divergence and evolution. Although the role of gene regulation in microevolutionary processes is gaining wide acceptance, most studies have only investigated the evolution of transcript levels, ignoring the potentially significant role of transcript structures. We argue that variation in alternative splicing plays an important and widely unexplored role in adaptation (e.g., by increasing transcriptome and/or proteome diversity, or buffering potentially deleterious genetic variation). New studies increasingly highlight the potential for independent evolution in alternative splicing and transcript level, providing alternative paths for selection to act upon. We propose that alternative splicing and transcript levels can provide contrasting, nonredundant mechanisms of equal importance for adaptive diversification of gene function and regulation.

PIC-Me: paralogs and isoforms classifier based on machine-learning approaches

Background

Paralogs formed through gene duplication and isoforms formed through alternative splicing have been important processes for increasing protein diversity and maintaining cellular homeostasis. Despite their recognized importance and the advent of large-scale genomic and transcriptomic analyses, paradoxically, accurate annotations of all gene loci to allow the identification of paralogs and isoforms remain surprisingly incomplete. In particular, the global analysis of the transcriptome of a non-model organism for which there is no reference genome is especially challenging.

Results

To reliably discriminate between the paralogs and isoforms in RNA-seq data, we redefined the pre-existing sequence features (sequence similarity, inverse count of consecutive identical or non-identical blocks, and match-mismatch fraction) previously derived from full-length cDNAs and EST sequences and described newly discovered genomic and transcriptomic features (twilight zone of protein sequence alignment and expression level difference). In addition, the effectiveness and relevance of the proposed features were verified with two widely used support vector machine (SVM) and random forest (RF) models. From nine RNA-seq datasets, all AUC (area under the curve) scores of ROC (receiver operating characteristic) curves were over 0.9 in the RF model and significantly higher than those in the SVM model.

Conclusions

In this study, using an RF model with five proposed RNA-seq features, we implemented our method called Paralogs and Isoforms Classifier based on Machine-learning approaches (PIC-Me) and showed that it outperformed an existing method. Finally, we envision that our tool will be a valuable computational resource for the genomics community to help with gene annotation and will aid in comparative transcriptomics and evolutionary genomics studies, especially those on non-model organisms.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-021-04229-x.

Alternative splicing and gene expression play contrasting roles in the parallel phenotypic evolution of a salmonid fish

Understanding the contribution of different molecular processes to evolution and development is crucial for identifying the mechanisms of adaptation. Here, we used RNA-sequencing data to test the importance of alternative splicing and differential gene expression in a case of parallel adaptive evolution, the replicated postglacial divergence of the salmonid fish Arctic charr (Salvelinus alpinus) into sympatric benthic and pelagic ecotypes across multiple independent lakes. We found that genes differentially spliced between ecotypes were mostly not differentially expressed (<6% overlap) and were involved in different biological processes. Differentially spliced genes were primarily enriched for muscle development and functioning, while differentially expressed genes were involved in metabolism, immunity and growth. Furthermore, alternative splicing and gene expression were mostly controlled by independent cis-regulatory quantitative trait loci (<3.4% overlap). Cis-regulatory regions were associated with the parallel divergence in splicing (16.5% of intron clusters) and expression (6.7%-10.1% of differentially expressed genes), indicating shared regulatory variation across ecotype pairs. Contrary to theoretical expectation, we found that differentially spliced genes tended to be highly central in regulatory networks ("hub genes") and were annotated to significantly more gene ontology terms compared to nondifferentially spliced genes, consistent with a higher level of pleiotropy. Together, our results suggest that the concerted regulation of alternative splicing and differential gene expression through different regulatory regions leads to the divergence of complementary processes important for local adaptation. This provides novel insights into the importance of contrasting but putatively complementary molecular processes in rapid parallel adaptive evolution.

Isolation and Characterization of APETALA3 Orthologs and Promoters from the Distylous Fagopyrum esculentum

Common buckwheat (Fagopyrum esculentum) produces distylous flowers with undifferentiated petaloid tepals, which makes it obviously different from flowers of model species. In model species Arabidopsis, APETALA3 (AP3) is expressed in petal and stamen and specifies petal and stamen identities during flower development. Combining with our previous studies, we found that small-scale gene duplication (GD) event and alternative splicing (AS) of common buckwheat AP3 orthologs resulted in FaesAP3_1, FaesAP3_2 and FaesAP3_2a. FaesAP3_2 and FaesAP3_2a were mainly expressed in the stamen of thrum and pin flower. Promoters functional analysis suggested that intense GUS staining was observed in the whole stamen in pFaesAP3_2::GUS transgenic Arabidopsis, while intense GUS staining was observed only in the filament of stamen in pFaesAP3_1::GUS transgenic Arabidopsis. These suggested that FaesAP3_1 and FaesAP3_2 had overlapping functions in specifying stamen filament identity and work together to determine normal stamen development. Additionally, FaesAP3_2 and FaesAP3_2a owned the similar ability to rescue stamen development of Arabidopsis ap3-3 mutant, although AS resulted in a frameshift mutation and consequent omission of the complete PI-derived motif and euAP3 motif of FaesAP3_2a. These suggested that the MIK region of AP3-like proteins was crucial for determining stamen identity, while the function of AP3-like proteins in specifying petal identity was gradually obtained after AP3 Orthologs acquiring a novel C-terminal euAP3 motif during the evolution of core eudicots. Our results also provide a clue to understanding the early evolution of the functional specificity of euAP3-type proteins involving in floral organ development in core eudicots, and also suggested that FaesAP3_2 holds the potential application for biotechnical engineering to develop a sterile male line of F. esculentum.

Transposable Elements in the Genome of Human Parasite Schistosoma mansoni: A Review

Transposable elements (TEs) are DNA sequences able to transpose within the host genome and, consequently, influence the dynamics of evolution in the species. Among the possible effects, TEs insertions may alter the expression and coding patterns of genes, leading to genomic innovations. Gene-duplication events, resulting from DNA segmental duplication induced by TEs transposition, constitute another important mechanism that contributes to the plasticity of genomes. This review aims to cover the current knowledge regarding TEs in the genome of the parasite Schistosoma mansoni, an agent of schistosomiasis-a neglected tropical disease affecting at least 250 million people worldwide. In this context, the literature concerning TEs description and TEs impact on the genomic architecture for S. mansoni was revisited, displaying evidence of TEs influence on schistosome speciation-mediated by bursts of transposition-and in gene-duplication events related to schistosome-host coevolution processes, as well several instances of TEs contribution into the coding sequences of genes. These findings indicate the relevant role of TEs in the evolution of the S. mansoni genome.

The divergence of alternative splicing between ohnologs in teleost fishes

BACKGROUND

Gene duplication and alternative splicing (AS) are two distinct mechanisms generating new materials for genetic innovations. The evolutionary link between gene duplication and AS is still controversial, due to utilizing duplicates from inconsistent ages of duplication events in earlier studies. With the aid of RNA-seq data, we explored evolutionary scenario of AS divergence between duplicates with ohnologs that resulted from the teleost genome duplication event in zebrafish, medaka, and stickleback.

RESULTS

Ohnologs in zebrafish have fewer AS forms compared to their singleton orthologs, supporting the function-sharing model of AS divergence between duplicates. Ohnologs in stickleback have more AS forms compared to their singleton orthologs, which supports the accelerated model of AS divergence between duplicates. The evolution of AS in ohnologs in medaka supports a combined scenario of the function-sharing and the accelerated model of AS divergence between duplicates. We also found a small number of ohnolog pairs in each of the three teleosts showed significantly asymmetric AS divergence. For example, the well-known ovary-factor gene cyp19a1a has no AS form but its ohnolog cyp19a1b has multiple AS forms in medaka, suggesting that functional divergence between duplicates might have result from AS divergence.

CONCLUSIONS

We found that a combined scenario of function-sharing and accelerated models for AS evolution in ohnologs in teleosts and rule out the independent model that assumes a lack of correlation between gene duplication and AS. Our study thus provided insights into the link between gene duplication and AS in general and ohnolog divergence in teleosts from AS perspective in particular.

Cancer type-specific alterations in actin genes: Worth a closer look?

Actins form a strongly conserved family of proteins that are central to the functioning of the actin cytoskeleton partaking in natural processes such as cell division, adhesion, contraction and migration. These processes, however, also occur during the various phases of cancer progression. Yet, surprisingly, alterations in the six human actin genes in cancer studies have received little attention and the focus was mostly on deregulated expression levels of actins and even more so of actin-binding or regulatory proteins. Starting from the early mutation work in the 1980s, we propose based on reviewing literature and data from patient cancer genomes that alterations in actin genes are different in distinct cancer subtypes, suggesting some specificity. These actin gene alterations include (missense) mutations, gene fusions and copy number alterations (deletions and amplifications) and we illustrate their occurrence for a limited number of examples including actin mutations in lymphoid cancers and nonmelanoma skin cancer and actin gene copy number alterations for breast, prostate and liver cancers. A challenge in the future will be to further sort out the specificity per actin gene, alteration type and cancer subtype. Even more challenging is (experimentally) distinguishing between cause and consequence: which alterations are passengers and which are involved in tumor progression of particular cancer subtypes?

Sense-oriented AluYRa1 elements provide a lineage-specific transcription environment for polyadenylation

Transposable elements cause alternative splicing (AS) in different ways, contributing to transcript diversification. Alternative polyadenylation (APA), one of the AS events, is related to the generation of mRNA isoforms in 70% of human genes. In this study, we tried to investigate AluYRa1s located at the terminal region of cynomolgus monkey genes, utilizing both computational analysis and molecular experimentation. We found that ten genes had AluYRa1 at their 3' end, and nine of these AluYRa1s were sense-oriented. Furthermore, in seven genes, AluYRa1s were expected to have a similar consensus sequence for polyadenylation cleavage. Additional computational analysis using the annotation files from the UCSC database showed that AluYRa1 was more involved in polyadenylation than in open reading frame exon splicing. To examine the extent of AluYRa1 involvement in polyadenylation, RNA-seq data from 30 normal cynomolgus monkeys were analyzed using TAPAS, a recently devised software that detects all the promising polyadenylation sites including APA sites. We observed that approximately 74% of possible polyadenylation sites in the analyzed genes were provided by sense-oriented AluYRa1. In conclusion, AluYRa1 is an Old-World monkey-specific TE, and its sense-oriented insertion at the 3'UTR region tends to provide a favorable environment for polyadenylation, diversifying gene transcripts.

Aberrant Splicing Events and Epigenetics in Viral Oncogenomics: Current Therapeutic Strategies

Global cancer incidence and mortality are on the rise. Although cancer is fundamentally a non-communicable disease, a large number of cancers are known to have a viral aetiology. A high burden of infectious agents (Human immunodeficiency virus (HIV), human papillomavirus (HPV), hepatitis B virus (HBV)) in certain Sub-Saharan African countries drives the rates of certain cancers. About one-third of all cancers in Africa are attributed to infection. Seven viruses have been identified with carcinogenic characteristics, namely the HPV, HBV, Hepatitis C virus (HCV), Epstein–Barr virus (EBV), Human T cell leukaemia virus 1 (HTLV-1), Kaposi’s Sarcoma Herpesvirus (KSHV), and HIV-1. The cellular splicing machinery is compromised upon infection, and the virus generates splicing variants that promote cell proliferation, suppress signalling pathways, inhibition of tumour suppressors, alter gene expression through epigenetic modification, and mechanisms to evade an immune response, promoting carcinogenesis. A number of these splice variants are specific to virally-induced cancers. Elucidating mechanisms underlying how the virus utilises these splice variants to maintain its latent and lytic phase will provide insights into novel targets for drug discovery. This review will focus on the splicing genomics, epigenetic modifications induced by and current therapeutic strategies against HPV, HBV, HCV, EBV, HTLV-1, KSHV and HIV-1.

Splicing Genomics Events in Cervical Cancer: Insights for Phenotypic Stratification and Biomarker Potency

Gynaecological cancers are attributed to the second most diagnosed cancers in women after breast cancer. On a global scale, cervical cancer is the fourth most common cancer and the most common cancer in developing countries with rapidly increasing mortality rates. Human papillomavirus (HPV) infection is a major contributor to the disease. HPV infections cause prominent cellular changes including alternative splicing to drive malignant transformation. A fundamental characteristic attributed to cancer is the dysregulation of cellular transcription. Alternative splicing is regulated by several splicing factors and molecular changes in these factors lead to cancer mechanisms such as tumour development and progression and drug resistance. The serine/arginine-rich (SR) proteins and heterogeneous ribonucleoproteins (hnRNPs) have prominent roles in modulating alternative splicing. Evidence shows molecular alteration and expression levels in these splicing factors in cervical cancer. Furthermore, aberrant splicing events in cancer-related genes lead to chemo- and radioresistance. Identifying clinically relevant modifications in alternative splicing events and splicing variants, in cervical cancer, as potential biomarkers for their role in cancer progression and therapy resistance is scrutinised. This review will focus on the molecular mechanisms underlying the aberrant splicing events in cervical cancer that may serve as potential biomarkers for diagnosis, prognosis, and novel drug targets.

Revising transcriptome assemblies with phylogenetic information

A common transcriptome assembly error is to mistake different transcripts of the same gene as transcripts from multiple closely related genes. This error is difficult to identify during assembly, but in a phylogenetic analysis such errors can be diagnosed from gene phylogenies where they appear as clades of tips from the same species with improbably short branch lengths. treeinform is a method that uses phylogenetic information across species to refine transcriptome assemblies within species. It identifies transcripts of the same gene that were incorrectly assigned to multiple genes and reassign them as transcripts of the same gene. The treeinform method is implemented in Agalma, available at https://bitbucket.org/caseywdunn/agalma, and the general approach is relevant in a variety of other contexts.

Changes in Alternative Splicing in Response to Domestication and Polyploidization in Wheat

Alternative splicing (AS) occurs extensively in eukaryotes as an important mechanism for regulating transcriptome complexity and proteome diversity, but variation in the AS landscape in response to domestication and polyploidization in crops is unclear. Hexaploid wheat (AABBDD, Triticum aestivum) has undergone two separate allopolyploidization events, providing an ideal model for studying AS changes during domestication and polyploidization events. In this study, we performed high-throughput transcriptome sequencing of roots and leaves from wheat species with varied ploidies, including wild diploids (A^bA^b, Triticum boeoticum) and tetraploids (AABB, Triticum dicoccoides), domesticated diploids (A^mA^m, Triticum monococcum) and tetraploids (AABB, Triticum dicoccum), hexaploid wheat (AABBDD, T aestivum), as well as newly synthesized hexaploids together with their parents. Approximately 22.1% of genes exhibited AS, with the major AS type being intron retention. The number of AS events decreased after domestication in both diploids and tetraploids. Moreover, the frequency of AS occurrence tended to decrease after polyploidization, consistent with the functional sharing model that proposes AS and duplicated genes are complementary in regulating transcriptome plasticity in polyploid crops. In addition, the subgenomes exhibited biased AS responses to polyploidization, and ∼87.1% of homeologs showed AS partitioning in hexaploid wheat. Interestingly, substitution of the D-subgenome modified 42.8% of AS patterns of the A- and B-subgenomes, indicating subgenome interplay reprograms AS profiles at a genome-wide level, although the causal-consequence relationship requires further study. Conclusively, our study shows that AS variation occurs extensively after polyploidization and domestication in wheat species.

Transcriptome-wide analysis and modelling of prognostic alternative splicing signatures in invasive breast cancer: a prospective clinical study

Aberrant alternative splicing (AS) has been highly involved in the tumorigenesis and progression of most cancers. The potential role of AS in invasive breast cancer (IBC) remains largely unknown. In this study, RNA sequencing of IBC samples from The Cancer Genome Atlas was acquired. AS events were screened by conducting univariate and multivariate Cox analysis and least absolute shrinkage and selection operator regression. In total, 2146 survival-related AS events were identified from 1551 parental genes, of which 93 were related to prognosis, and a prognostic marker model containing 14 AS events was constructed. We also constructed the regulatory network of splicing factors (SFs) and AS events, and identified DDX39B as the node SF gene, and verified the accuracy of the network through experiments. Next, we performed quantitative real-time reverse transcription polymerase chain reaction (qRT-PCR) in triple negative breast cancer patients with different responses to neoadjuvant chemotherapy, and found that the exon-specific expression of EPHX2, C6orf141, and HERC4 was associated with the different status of patients that received neoadjuvant chemotherapy. In conclusion, this study found that DDX39B, EPHX2 (exo7), and HERC4 (exo23) can be used as potential targets for the treatment of breast cancer, which provides a new idea for the treatment of breast cancer.

Salmo salar glucocorticoid receptors analyses of alternative splicing variants under stress conditions

Cortisol is the main corticosteroid in teleosts, exerting multiple functions by activating glucocorticoid receptors (GR). Most teleost species have two GR genes, gr-1 and gr-2. Some teleost also presents two splice variants for gr-1; gr-1a and gr-1b. In this study, we report for first time the presence of 2 homeologous genes for gr-1 and gr-2, located on chromosomes 4q-13q (gr-1) and 5p-9q (gr-2) of the Salmo salar genome. Furthermore, our results describe gr-1 splice variants derived from chromosome 4 and 13, sharing typical teleost GR elements such as the 9 amino acid insertion in the DNA binding domain (DBD) and variations in the length of the ligand binding domain (LBD). Three splice variants were predicted for the gr-2 homeologous gene in chromosome 5, with differences of a 5 amino acid insertion in the DBD. We also identified an uncommon truncated gr-2 gene in chromosome 9 in salmon, which lacked the DBD and LBD domains. Finally, by designing specific primers for each predicted splice variant, we validated and evaluated the expression of their transcripts in S. salar subjected to stress caused by stocking density. Differences were observed in the expression of all identified mRNAs, revealing that gr-1 and gr-2 splice variants were upregulated in head kidney and gills of post-stressed fish. In conclusion, our findings suggest that from specific salmonid genomic duplication (125 MYA), two gene copies of each GR receptor were generated in S. salar. The identified splice variants could contribute to the variability of GR receptor complex modulation expression during stressful events, leading to variations in physiological responses in fish.

Gene Fusions Derived by Transcriptional Readthrough are Driven by Segmental Duplication in Human

Gene fusion occurs when two or more individual genes with independent open reading frames becoming juxtaposed under the same open reading frame creating a new fused gene. A small number of gene fusions described in detail have been associated with novel functions, for example, the hominid-specific PIPSL gene, TNFSF12, and the TWE-PRIL gene family. We use Sequence Similarity Networks and species level comparisons of great ape genomes to identify 45 new genes that have emerged by transcriptional readthrough, that is, transcription-derived gene fusion. For 35 of these putative gene fusions, we have been able to assess available RNAseq data to determine whether there are reads that map to each breakpoint. A total of 29 of the putative gene fusions had annotated transcripts (9/29 of which are human-specific). We carried out RT-qPCR in a range of human tissues (placenta, lung, liver, brain, and testes) and found that 23 of the putative gene fusion events were expressed in at least one tissue. Examining the available ribosome foot-printing data, we find evidence for translation of three of the fused genes in human. Finally, we find enrichment for transcription-derived gene fusions in regions of known segmental duplication in human. Together, our results implicate chromosomal structural variation brought about by segmental duplication with the emergence of novel transcripts and translated protein products.

Single-molecule real-time sequencing identifies massive full-length cDNAs and alternative-splicing events that facilitate comparative and functional genomics study in the hexaploid crop sweet potato

Background

Sweet potato (Ipomoea batatas (L.) Lam.) is one of the most important crops in many developing countries and provides a candidate source of bioenergy. However, neither a complete reference genome nor large-scale full-length cDNA sequences for this outcrossing hexaploid crop are available, which in turn impedes progress in research studies in I. batatas functional genomics and molecular breeding.

Methods

In this study, we sequenced full-length transcriptomes in I. batatas and its diploid ancestor I. trifida by single-molecule real-time sequencing and Illumina second-generation sequencing technologies. With the generated datasets, we conducted comprehensive intraspecific and interspecific sequence analyses and experimental characterization.

Results

A total of 53,861/51,184 high-quality long-read transcripts were obtained, which covered about 10,439/10,452 loci in the I. batatas/I. trifida genome. These datasets enabled us to predict open reading frames successfully in 96.83%/96.82% of transcripts and identify 34,963/33,637 full-length cDNA sequences, 1,401/1,457 transcription factors, 25,315/27,090 simple sequence repeats, 1,656/1,389 long non-coding RNAs, and 5,251/8,901 alternative splicing events. Approximately, 32.34%/38.54% of transcripts and 46.22%/51.18% multi-exon transcripts underwent alternative splicing in I. batatas/I. trifida. Moreover, we validated one alternative splicing event in each of 10 genes and identified tuberous-root-specific expressed isoforms from a starch-branching enzyme, an alpha-glucan phosphorylase, a neutral invertase, and several ABC transporters. Overall, the collection and analysis of large-scale long-read transcripts generated in this study will serve as a valuable resource for the I. batatas research community, which may accelerate the progress in its structural, functional, and comparative genomics studies.

Cooperative evolution of two different TEs results in lineage-specific novel transcripts in the BLOC1S2 gene

BACKGROUND

The BLOC1S2 gene encodes the multifunctional protein BLOS2, a shared subunit of two lysosomal trafficking complexes: i) biogenesis of lysosome-related organelles complex-1 and i) BLOC-1-related complex. In our previous study, we identified an intriguing unreported transcript of the BLOC1S2 gene that has a novel exon derived from two transposable elements (TEs), MIR and AluSp. To investigate the evolutionary footprint and molecular mechanism of action of this transcript, we performed PCR and RT-PCR experiments and sequencing analyses using genomic DNA and RNA samples from humans and various non-human primates.

RESULTS

The results showed that the MIR element had integrated into the genome of our common ancestor, specifically in the BLOC1S2 gene region, before the radiation of all primate lineages and that the AluSp element had integrated into the genome of our common ancestor, fortunately in the middle of the MIR sequences, after the divergence of Old World monkeys and New World monkeys. The combined MIR and AluSp sequences provide a 3' splice site (AG) and 5' splice site (GT), respectively, and generate the Old World monkey-specific transcripts. Moreover, branch point sequences for the intron removal process are provided by the MIR and AluSp combination.

CONCLUSIONS

We show for the first time that sequential integration into the same location and sequence divergence events of two different TEs generated lineage-specific transcripts through sequence collaboration during primate evolution.

Evolutionary Divergence of Duplicated Hsf Genes in Populus

Heat shock transcription factors (Hsfs), which function as the activator of heat shock proteins (Hsps), play multiple roles in response to environmental stress and the development of plants. The Hsf family had experienced gene expansion via whole-genome duplication from a single cell algae to higher plants. However, how the Hsf gene family went through evolutionary divergence after genome duplication is unknown. As a model wood species, Populus trichocarpa is widely distributed in North America with various ecological and climatic environments. In this study, we used P. trichocarpa as materials and identified the expression divergence of the PtHsf gene family in developmental processes, such as dormant bud formation and opening, catkins development, and in response to environments. Through the co-expression network, we further discovered the divergent co-expressed genes that related to the functional divergence of PtHsfs. Then, we studied the alternative splicing events, single nucleotide polymorphism distribution and tertiary structures of members of the PtHsf gene family. In addition to expression divergence, we uncovered the evolutionary divergence in the protein level which may be important to new function formations and for survival in changing environments. This study comprehensively analyzed the evolutionary divergence of a member of the PtHsf gene family after genome duplication, paving the way for further gene function analysis and genetic engineering.

Variants in the IL17 pathway genes are associated with atopic asthma and atopy makers in a South American population

BACKGROUND

Asthma is a complex disorder with multiple phenotypes which can influence its severity and response to treatment. The TH17 lymphocytes producing IL-17A and IL17-F cytokines, may have a role on asthma inflammation. The aim of our study was to evaluate the association between genetic variants in IL17 pathway genes with asthma and atopy markers.

MATERIALS AND METHODS

Genotyping was performed using a commercial panel in 1245 participants of SCAALA cohort. The study included 91 SNVs in IL-17 pathway genes. Logistic regressions for asthma and atopy markers were performed using PLINK 1.9. In silico analyses were performed using rSNPbase, RegulomeDB, and Gtex portal for in silico gene expression.

RESULTS AND DISCUSSION

The T allele of rs1974226 in IL17A was positively associated with asthma (OR: 1.37; 95% CI 1.02-1.82). Also, the T allele of rs279548 was positively associated with asthma (OR: 1.30; 95% CI 1.02-1.64), atopy (OR: 1.62; 95% CI 1.05-2.50) and increased expression of the IL17RC in lung and whole blood tissues. The others genetic variants in the IL17 pathways genes were associated with both protection and risk for asthma development as well as with IgE levels.

CONCLUSION

The genetic variants in IL-17-related genes are associated with the atopic asthma phenotype and IgE production.

Genome-wide identification and analysis of the ALTERNATIVE OXIDASE gene family in diploid and hexaploid wheat

A comprehensive understanding of wheat responses to environmental stress will contribute to the long-term goal of feeding the planet. ALERNATIVE OXIDASE (AOX) genes encode proteins involved in a bypass of the electron transport chain and are also known to be involved in stress tolerance in multiple species. Here, we report the identification and characterization of the AOX gene family in diploid and hexaploid wheat. Four genes each were found in the diploid ancestors Triticum urartu, and Aegilops tauschii, and three in Aegilops speltoides. In hexaploid wheat (Triticum aestivum), 20 genes were identified, some with multiple splice variants, corresponding to a total of 24 proteins for those with observed transcription and translation. These proteins were classified as AOX1a, AOX1c, AOX1e or AOX1d via phylogenetic analysis. Proteins lacking most or all signature AOX motifs were assigned to putative regulatory roles. Analysis of protein-targeting sequences suggests mixed localization to the mitochondria and other organelles. In comparison to the most studied AOX from Trypanosoma brucei, there were amino acid substitutions at critical functional domains indicating possible role divergence in wheat or grasses in general. In hexaploid wheat, AOX genes were expressed at specific developmental stages as well as in response to both biotic and abiotic stresses such as fungal pathogens, heat and drought. These AOX expression patterns suggest a highly regulated and diverse transcription and expression system. The insights gained provide a framework for the continued and expanded study of AOX genes in wheat for stress tolerance through breeding new varieties, as well as resistance to AOX-targeted herbicides, all of which can ultimately be used synergistically to improve crop yield.

The Role of Transposable Elements in Emergence of Metazoa

Systems initially emerged for protecting genomes against insertions of transposable elements and represented by mechanisms of splicing regulation, RNA-interference, and epigenetic factors have played a key role in the evolution of animals. Many studies have shown inherited transpositions of mobile elements in embryogenesis and preservation of their activities in certain tissues of adult organisms. It was supposed that on the emergence of Metazoa the self-regulation mechanisms of transposons related with the gene networks controlling their activity could be involved in intercellular cell coordination in the cascade of successive divisions with differentiated gene expression for generation of tissues and organs. It was supposed that during evolution species-specific features of transposons in the genomes of eukaryotes could form the basis for creation of dynamically related complexes of systems for epigenetic regulation of gene expression. These complexes could be produced due to the influence of noncoding transposon-derived RNAs on DNA methylation, histone modifications, and processing of alternative splicing variants, whereas the mobile elements themselves could be directly involved in the regulation of gene expression in cis and in trans. Transposons are widely distributed in the genomes of eukaryotes; therefore, their activation can change the expression of specific genes. In turn, this can play an important role in cell differentiation during ontogenesis. It is supposed that transposons can form a species-specific pattern for control of gene expression, and that some variants of this pattern can be favorable for adaptation. The presented data indicate the possible influence of transposons in karyotype formation. It is supposed that transposon localization relative to one another and to protein-coding genes can influence the species-specific epigenetic regulation of ontogenesis.

Alternative splicing after gene duplication drives CEACAM1-paralog diversification in the horse

Background

The CEA gene family is one of the most rapidly evolving gene families in the human genome. The founder gene of the family is thought to be an ancestor of the inhibitory immune checkpoint molecule CEACAM1. Comprehensive analyses of mammalian genomes showed that the CEA gene family is subject to tremendous gene family expansion and contraction events in different mammalian species. While in some species (e.g. rabbits) less than three CEACAM1 related genes exist, were in others (certain microbat species) up to 100 CEACAM1 paralogs identified. We have recently reported that the horse has also an extended CEA gene family. Since mechanisms of gene family expansion and diversification are not well understood we aimed to analyze the equine CEA gene family in detail.

Results

We found that the equine CEA gene family contains 17 functional CEACAM1-related genes. Nine of them were secreted molecules and eight CEACAMs contain transmembrane and cytoplasmic domain exons, the latter being in the focus of the present report. Only one (CEACAM41) gene has exons coding for activating signaling motifs all other CEACAM1 paralogs contain cytoplasmic exons similar to that of the inhibitory receptor CEACAM1. However, cloning of cDNAs showed that only one CEACAM1 paralog contain functional immunoreceptor tyrosine-based inhibitory motifs in its cytoplasmic tail. Three receptors have acquired a stop codon in the transmembrane domain and two have lost their inhibitory motifs due to alternative splicing events. In addition, alternative splicing eliminated the transmembrane exon sequence of the putative activating receptor, rendering it to a secreted molecule. Transfection of eukaryotic cells with FLAG-tagged alternatively spliced CEACAMs indicates that they can be expressed in vivo. Thus detection of CEACAM41 mRNA in activated PBMC suggests that CEACAM41 is secreted by lymphoid cells upon activation.

Conclusions

The results of our study demonstrate that alternative splicing after gene duplication is a potent mechanism to accelerate functional diversification of the equine CEA gene family members. This potent mechanism has created novel CEACAM receptors with unique signaling capacities and secreted CEACAMs which potentially enables equine lymphoid cells to control distantly located immune cells.

Genome-wide characterization of differential transcript usage in Arabidopsis thaliana

Alternative splicing and the usage of alternate transcription start- or stop sites allows a single gene to produce multiple transcript isoforms. Most plant genes express certain isoforms at a significantly higher level than others, but under specific conditions this expression dominance can change, resulting in a different set of dominant isoforms. These events of differential transcript usage (DTU) have been observed for thousands of Arabidopsis thaliana, Zea mays and Vitis vinifera genes, and have been linked to development and stress response. However, neither the characteristics of these genes, nor the implications of DTU on their protein coding sequences or functions, are currently well understood. Here we present a dataset of isoform dominance and DTU for all genes in the AtRTD2 reference transcriptome based on a protocol that was benchmarked on simulated data and validated through comparison with a published reverse transciptase-polymerase chain reaction panel. We report DTU events for 8148 genes across 206 public RNA-Seq samples, and find that protein sequences are affected in 22% of the cases. The observed DTU events show high consistency across replicates, and reveal reproducible patterns in response to treatment and development. We also demonstrate that genes with different evolutionary ages, expression breadths and functions show large differences in the frequency at which they undergo DTU, and in the effect that these events have on their protein sequences. Finally, we showcase how the generated dataset can be used to explore DTU events for genes of interest or to find genes with specific DTU in samples of interest.