Increase of functional diversity by alternative splicing

Genome-wide identification, abiotic stress, and expression analysis of PYL family in Tartary buckwheat (Fagopyrum tataricum (L.) Gaertn.) during grain development

BACKGROUND

Abscisic acid (ABA) is a plant hormone that plays an important role in plant resistance to drought, salinity, cold, and pathogens. It is also important for regulating plant growth and development. Pyrabactin resistance/pyr1-like/regulatory components of the ABA receptor (PYL/RCAR) are ABA receptor proteins in plants and the core of ABA signal transduction pathways in plant regulatory factors. At present, there are no reports on the PYL family of Tartary buckwheat.

RESULTS

In this study, 19 paralogous form PYL genes in buckwheat were identified at the whole-genome level and named FtPYL1-FtPYL19 according to their positions on chromosomes. We further analyzed the gene structure, conserved motifs, cis-acting elements, gene duplication, phylogenetic relationships, and expression patterns under different stress treatments and during grain development of the 19 paralogous form PYL genes in Tartary buckwheat. The FtPYL gene exhibits a single exonic gene structure for about 68.4% of the duplicated forms from the total paralogous forms. The remaining subfamilies, such as I and II, contain three exons and two exons (e.g., FtPYL19), respectively. Nineteen FtPYL genes were evenly distributed across the eight chromosomes, with at least one FtPYL gene on each chromosome. In the FtPYL gene family, there was one tandem repeat event and five gene duplication events. We investigated the gene expression levels of FtPYL gene under four abiotic stresses and different stages of grain development. Under drought stress (PEG6000), the relative expression levels of FtPYL14 and FtPYL15 increased by fourfold. Under high temperature stress (38℃), the relative expression level of FtPYL16 dropped to 0.12, and that of FtPYL17 fell to 0.22. At different stages of grain development, the gene expression level of FtPY15 is extremely high at 19 D. The relative expression level of FtPYL7 in roots and stems reaches up to approximately 450, and the relative expression level of FtPYL10 in 13 D also reaches up to 248. In this study, the PYL gene family of Tartary buckwheat was identified and analyzed based on the whole genome, and 19 paralogous form FtPYL genes of Tartary buckwheat were bioinformatically analyzed. The expression patterns of 19 paralogous form FtPYL genes in Tartary buckwheat cultivars under different stress treatments and during grain development were analyzed. It was found that the FtPYL gene played an important role in grain development.

High Overexpression of SiAAP9 Leads to Growth Inhibition and Protein Ectopic Localization in Transgenic Arabidopsis

Amino acid permeases (AAPs) transporters are crucial for the long-distance transport of amino acids in plants, from source to sink. While Arabidopsis and rice have been extensively studied, research on foxtail millet is limited. This study identified two transcripts of SiAAP9, both of which were induced by NO3- and showed similar expression patterns. The overexpression of SiAAP9L and SiAAP9S in Arabidopsis inhibited plant growth and seed size, although SiAAP9 was found to transport more amino acids into seeds. Furthermore, SiAAP9-OX transgenic Arabidopsis showed increased tolerance to high concentrations of glutamate (Glu) and histidine (His). The high overexpression level of SiAAP9 suggested its protein was not only located on the plasma membrane but potentially on other organelles, as well. Interestingly, sequence deletion reduced SiAAP9's sensitivity to Brefeldin A (BFA), and SiAAP9 had ectopic localization on the endoplasmic reticulum (ER). Protoplast amino acid uptake experiments indicated that SiAAP9 enhanced Glu transport into foxtail millet cells. Overall, the two transcripts of SiAAP9 have similar functions, but SiAAP9L shows a higher colocalization with BFA compartments compared to SiAAP9S. Our research identifies a potential candidate gene for enhancing the nutritional quality of foxtail millet through breeding.

Widespread variation in molecular interactions and regulatory properties among transcription factor isoforms

Most human Transcription factors (TFs) genes encode multiple protein isoforms differing in DNA binding domains, effector domains, or other protein regions. The global extent to which this results in functional differences between isoforms remains unknown. Here, we systematically compared 693 isoforms of 246 TF genes, assessing DNA binding, protein binding, transcriptional activation, subcellular localization, and condensate formation. Relative to reference isoforms, two-thirds of alternative TF isoforms exhibit differences in one or more molecular activities, which often could not be predicted from sequence. We observed two primary categories of alternative TF isoforms: "rewirers" and "negative regulators", both of which were associated with differentiation and cancer. Our results support a model wherein the relative expression levels of, and interactions involving, TF isoforms add an understudied layer of complexity to gene regulatory networks, demonstrating the importance of isoform-aware characterization of TF functions and providing a rich resource for further studies.

Exonic splicing code and coordination of divalent metals in proteins

Exonic sequences contain both protein-coding and RNA splicing information but the interplay of the protein and splicing code is complex and poorly understood. Here, we have studied traditional and auxiliary splicing codes of human exons that encode residues coordinating two essential divalent metals at the opposite ends of the Irving-Williams series, a universal order of relative stabilities of metal-organic complexes. We show that exons encoding Zn2+-coordinating amino acids are supported much less by the auxiliary splicing motifs than exons coordinating Ca2+. The handicap of the former is compensated by stronger splice sites and uridine-richer polypyrimidine tracts, except for position -3 relative to 3' splice junctions. However, both Ca2+ and Zn2+ exons exhibit close-to-constitutive splicing in multiple tissues, consistent with their critical importance for metalloprotein function and a relatively small fraction of expendable, alternatively spliced exons. These results indicate that constraints imposed by metal coordination spheres on RNA splicing have been efficiently overcome by the plasticity of exon-intron architecture to ensure adequate metalloprotein expression.

Genome-wide identification, evolution, and role of SPL gene family in beet (Beta vulgaris L.) under cold stress

BACKGROUND

SPL transcription factors play vital roles in regulating plant growth, development, and abiotic stress responses. Sugar beet (Beta vulgaris L.), one of the world's main sugar-producing crops, is a major source of edible and industrial sugars for humans. Although the SPL gene family has been extensively identified in other species, no reports on the SPL gene family in sugar beet are available.

RESULTS

Eight BvSPL genes were identified at the whole-genome level and were renamed based on their positions on the chromosome. The gene structure, SBP domain sequences, and phylogenetic relationship with Arabidopsis were analyzed for the sugar beet SPL gene family. The eight BvSPL genes were divided into six groups (II, IV, V, VI, VII, and VIII). Of the BvSPL genes, no tandem duplication events were found, but one pair of segmental duplications was present. Multiple cis-regulatory elements related to growth and development were identified in the 2000-bp region upstream of the BvSPL gene start codon (ATG). Using quantitative real-time polymerase chain reaction (qRT-PCR), the expression profiles of the eight BvSPL genes were examined under eight types of abiotic stress and during the maturation stage. BvSPL transcription factors played a vital role in abiotic stress, with BvSPL3 and BvSPL6 being particularly noteworthy.

CONCLUSION

Eight sugar beet SPL genes were identified at the whole-genome level. Phylogenetic trees, gene structures, gene duplication events, and expression profiles were investigated. The qRT-PCR analysis indicated that BvSPLs play a substantial role in the growth and development of sugar beet, potentially participating in the regulation of root expansion and sugar accumulation.

Identification and characterization of mitogen-activated protein kinase kinase 4 (MKK4) from the mud crab Scylla paramamosain in response to Vibrio alginolyticus and White Spot Syndrome Virus (WSSV)

Mitogen-activated protein kinase kinase 4 (MKK4), serves as a critical component of the mitogen-activated protein kinase signaling pathway, facilitating the direct phosphorylation and activation of the c-Jun N-terminal kinase (JNK) and p38 families of MAP kinases in response to environmental stresses. In the current research, we identified two MKK4 subtypes, namely SpMKK4-1 and SpMKK4-2, from Scylla paramamosain, followed by the analysis of their molecular characteristics and tissue distributions. The expression of SpMKK4s was induced upon WSSV and Vibrio alginolyticus challenges, and the bacteria clearance capacity and antimicrobial peptide (AMP) genes' expression upon bacterial infection were significantly decreased after knocking down SpMKK4s. Additionally, the overexpression of both SpMKK4s remarkably activated NF-κB reporter plasmid in HEK293T cells, suggesting the activation of the NF-κB signaling pathway. These results indicated the participation of SpMKK4s in the innate immunity of crabs, which shed light on a better understanding of the mechanisms through which MKK4s regulate innate immunity.

Identification and verification of plasma protein biomarkers that accurately identify an ectopic pregnancy

BACKGROUND

Differentiating between a normal intrauterine pregnancy (IUP) and abnormal conditions including early pregnancy loss (EPL) or ectopic pregnancy (EP) is a major clinical challenge in early pregnancy. Currently, serial β-human chorionic gonadotropin (β-hCG) and progesterone are the most commonly used plasma biomarkers for evaluating pregnancy prognosis when ultrasound is inconclusive. However, neither biomarker can predict an EP with sufficient and reproducible accuracy. Hence, identification of new plasma biomarkers that can accurately diagnose EP would have great clinical value.

METHODS

Plasma was collected from a discovery cohort of 48 consenting women having an IUP, EPL, or EP. Samples were analyzed by liquid chromatography-tandem mass spectrometry (LC-MS/MS) followed by a label-free proteomics analysis to identify significant changes between pregnancy outcomes. A panel of 14 candidate biomarkers were then verified in an independent cohort of 74 women using absolute quantitation by targeted parallel reaction monitoring mass spectrometry (PRM-MS) which provided the capacity to distinguish between closely related protein isoforms. Logistic regression and Lasso feature selection were used to evaluate the performance of individual biomarkers and panels of multiple biomarkers to predict EP.

RESULTS

A total of 1391 proteins were identified in an unbiased plasma proteome discovery. A number of significant changes (FDR ≤ 5%) were identified when comparing EP vs. non-EP (IUP + EPL). Next, 14 candidate biomarkers (ADAM12, CGA, CGB, ISM2, NOTUM, PAEP, PAPPA, PSG1, PSG2, PSG3, PSG9, PSG11, PSG6/9, and PSG8/1) were verified as being significantly different between EP and non-EP in an independent cohort (FDR ≤ 5%). Using logistic regression models, a risk score for EP was calculated for each subject, and four multiple biomarker logistic models were identified that performed similarly and had higher AUCs than models with single predictors.

CONCLUSIONS

Overall, four multivariable logistic models were identified that had significantly better prediction of having EP than those logistic models with single biomarkers. Model 4 (NOTUM, PAEP, PAPPA, ADAM12) had the highest AUC (0.987) and accuracy (96%). However, because the models are statistically similar, all markers in the four models and other highly correlated markers should be considered in further validation studies.

Resurrecting the alternative splicing landscape of archaic hominins using machine learning

Alternative splicing contributes to adaptation and divergence in many species. However, it has not been possible to directly compare splicing between modern and archaic hominins. Here, we unmask the recent evolution of this previously unobservable regulatory mechanism by applying SpliceAI, a machine-learning algorithm that identifies splice-altering variants (SAVs), to high-coverage genomes from three Neanderthals and a Denisovan. We discover 5,950 putative archaic SAVs, of which 2,186 are archaic-specific and 3,607 also occur in modern humans via introgression (244) or shared ancestry (3,520). Archaic-specific SAVs are enriched in genes that contribute to traits potentially relevant to hominin phenotypic divergence, such as the epidermis, respiration and spinal rigidity. Compared to shared SAVs, archaic-specific SAVs occur in sites under weaker selection and are more common in genes with tissue-specific expression. Further underscoring the importance of negative selection on SAVs, Neanderthal lineages with low effective population sizes are enriched for SAVs compared to Denisovan and shared SAVs. Finally, we find that nearly all introgressed SAVs in humans were shared across the three Neanderthals, suggesting that older SAVs were more tolerated in human genomes. Our results reveal the splicing landscape of archaic hominins and identify potential contributions of splicing to phenotypic differences among hominins.

Identification of the NAC Transcription Factor Family during Early Seed Development in Akebia trifoliata (Thunb.) Koidz

This study aimed to gain an understanding of the possible function of NACs by examining their physicochemical properties, structure, chromosomal location, and expression. Being a family of plant-specific transcription factors, NAC (petunia no apical meristem and Arabidopsis thaliana ATAF1, ATAF2, and CUC2) is involved in plant growth and development. None of the NAC genes has been reported in Akebia trifoliata (Thunb.) Koidz (A. trifoliata). In this study, we identified 101 NAC proteins (AktNACs) in the A. trifoliata genome by bioinformatic analysis. One hundred one AktNACs were classified into the following twelve categories based on the phylogenetic analysis of NAC protein: NAC-a, NAC-b, NAC-c, NAC-d, NAC-e, NAC-f, NAC-g, NAC-h, NAC-i, NAC-j, NAC-k, and NAC-l. The accuracy of the clustering results was demonstrated based on the gene structure and conserved motif analysis of AktNACs. In addition, we identified 44 pairs of duplication genes, confirming the importance of purifying selection in the evolution of AktNACs. The morphology and microstructure of early A. trifoliata seed development showed that it mainly underwent rapid cell division, seed enlargement, embryo formation and endosperm development. We constructed AktNACs co-expression network and metabolite correlation network based on transcriptomic and metabolomic data of A. trifoliata seeds. The results of the co-expression network showed that 25 AtNAC genes were co-expressed with 233 transcription factors. Metabolite correlation analysis showed that 23 AktNACs were highly correlated with 28 upregulated metabolites. Additionally, 25 AktNACs and 235 transcription factors formed co-expression networks with 141 metabolites, based on correlation analysis involving AktNACs, transcription factors, and metabolites. Notably, AktNAC095 participates in the synthesis of 35 distinct metabolites. Eight of these metabolites, strongly correlated with AktNAC095, were upregulated during early seed development. These studies may provide insight into the evolution, possible function, and expression of AktNACs genes.

Identification of sex-specific splicing via comparative transcriptome analysis of gonads from sea cucumber Apostichopus japonicus

Alternative splicing (AS) is an essential post-transcriptional regulation mechanism for sex differentiation and gonadal development, which has rarely been reported in marine invertebrates. Sea cucumber (Apostichopus japonicus) is an economically important marine benthic echinoderm with a potential XX/XY sex determination mechanism, whose molecular mechanism in the gonadal differentiation has not been clearly understood. In this study, we analyzed available RNA-seq datasets of male and female gonads to explore if AS mechanism exerts an essential function in sex differentiation and gonadal development of A. japonicus. In our results, a total of 20,338 AS events from 7219 alternatively spliced genes, and 189 sexually differential alternative splicing (DAS) events from 156 genes were identified in gonadal transcriptome of sea cucumber. Gene Ontology analysis indicated that these DAS genes were significantly enriched in spermatogenesis-related GO terms. Maximal Clique Centrality (MCC) was then applied for protein-protein interaction (PPI) analysis to search for protein interactions and hub DAS gene. Among all DAS genes, we identified 10 DAS genes closely related to spermatogenesis and (or) sperm motility and a hub gene dnah1. Thus, this study revealed that alternative isoforms were generated from certain genes in female and male gonads through alternative splicing, which may provide direct evidence that alternative splicing mechanisms participate in female and male gonads. These results suggested a novel perspective for explaining the molecular mechanisms underlying gonadal differentiation between male and female sea cucumbers.

Re-evaluating the impact of alternative RNA splicing on proteomic diversity

Alternative splicing (AS) constitutes a mechanism by which protein-coding genes and long non-coding RNA (lncRNA) genes produce more than a single mature transcript. From plants to humans, AS is a powerful process that increases transcriptome complexity. Importantly, splice variants produced from AS can potentially encode for distinct protein isoforms which can lose or gain specific domains and, hence, differ in their functional properties. Advances in proteomics have shown that the proteome is indeed diverse due to the presence of numerous protein isoforms. For the past decades, with the help of advanced high-throughput technologies, numerous alternatively spliced transcripts have been identified. However, the low detection rate of protein isoforms in proteomic studies raised debatable questions on whether AS contributes to proteomic diversity and on how many AS events are really functional. We propose here to assess and discuss the impact of AS on proteomic complexity in the light of the technological progress, updated genome annotation, and current scientific knowledge.

Isoform function prediction by Gene Ontology embedding

MOTIVATION

High-resolution annotation of gene functions is a central task in functional genomics. Multiple proteoforms translated from alternatively spliced isoforms from a single gene are actual function performers and greatly increase the functional diversity. The specific functions of different isoforms can decipher the molecular basis of various complex diseases at a finer granularity. Multi-instance learning (MIL)-based solutions have been developed to distribute gene(bag)-level Gene Ontology (GO) annotations to isoforms(instances), but they simply presume that a particular annotation of the gene is responsible by only one isoform, neglect the hierarchical structures and semantics of massive GO terms (labels), or can only handle dozens of terms.

RESULTS

We propose an efficacy approach IsofunGO to differentiate massive functions of isoforms by GO embedding. Particularly, IsofunGO first introduces an attributed hierarchical network to model massive GO terms, and a GO network embedding strategy to learn compact representations of GO terms and project GO annotations of genes into compressed ones, this strategy not only explores and preserves hierarchy between GO terms but also greatly reduces the prediction load. Next, it develops an attention-based MIL network to fuse genomics and transcriptomics data of isoforms and predict isoform functions by referring to compressed annotations. Extensive experiments on benchmark datasets demonstrate the efficacy of IsofunGO. Both the GO embedding and attention mechanism can boost the performance and interpretability.

AVAILABILITYAND IMPLEMENTATION

The code of IsofunGO is available at http://www.sdu-idea.cn/codes.php?name=IsofunGO.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

SRRM2 organizes splicing condensates to regulate alternative splicing

SRRM2 is a nuclear-speckle marker containing multiple disordered domains, whose dysfunction is associated with several human diseases. Using mainly EGFP-SRRM2 knock-in HEK293T cells, we show that SRRM2 forms biomolecular condensates satisfying most hallmarks of liquid-liquid phase separation, including spherical shape, dynamic rearrangement, coalescence and concentration dependence supported by in vitro experiments. Live-cell imaging shows that SRRM2 organizes nuclear speckles along the cell cycle. As bona-fide splicing factor present in spliceosome structures, SRRM2 deficiency induces skipping of cassette exons with short introns and weak splice sites, tending to change large protein domains. In THP-1 myeloid-like cells, SRRM2 depletion compromises cell viability, upregulates differentiation markers, and sensitizes cells to anti-leukemia drugs. SRRM2 induces a FES splice isoform that attenuates innate inflammatory responses, and MUC1 isoforms that undergo shedding with oncogenic properties. We conclude that SRRM2 acts as a scaffold to organize nuclear speckles, regulating alternative splicing in innate immunity and cell homeostasis.

Genome-wide identification and characterization of the SPL gene family and its expression in the various developmental stages and stress conditions in foxtail millet (Setaria italica)

BACKGROUND

Among the major transcription factors, SPL plays a crucial role in plant growth, development, and stress response. Foxtail millet (Setaria italica), as a C4 crop, is rich in nutrients and is beneficial to human health. However, research on the foxtail millet SPL (SQUAMOSA PROMOTER BINDING-LIKE) gene family is limited.  RESULTS: In this study, a total of 18 SPL genes were identified for the comprehensive analysis of the whole genome of foxtail millet. These SiSPL genes were divided into seven subfamilies (I, II, III, V, VI, VII, and VIII) according to the classification of the Arabidopsis thaliana SPL gene family. Structural analysis of the SiSPL genes showed that the number of introns in subfamilies I and II were much larger than others, and the promoter regions of SiSPL genes were rich in different cis-acting elements. Among the 18 SiSPL genes, nine genes had putative binding sites with foxtail millet miR156. No tandem duplication events were found between the SiSPL genes, but four pairs of segmental duplications were detected. The SiSPL genes expression were detected in different tissues, which was generally highly expressed in seeds development process, especially SiSPL6 and SiSPL16, which deserve further study. The results of the expression levels of SiSPL genes under eight types of abiotic stresses showed that many stress responsive genes, especially SiSPL9, SiSPL10, and SiSPL16, were highly expressed under multiple stresses, which deserves further attention.

CONCLUSIONS

In this research, 18 SPL genes were identified in foxtail millet, and their phylogenetic relationships, gene structural features, duplication events, gene expression and potential roles in foxtail millet development were studied. The findings provide a new perspective for the mining of the excellent SiSPL gene and the molecular breeding of foxtail millet.

Genome-wide identification of associations between enhancer and alternative splicing in human and mouse

Background

Alternative splicing (AS) increases the diversity of transcriptome and could fine-tune the function of genes, so that understanding the regulation of AS is vital. AS could be regulated by many different cis-regulatory elements, such as enhancer. Enhancer has been experimentally proved to regulate AS in some genes. However, there is a lack of genome-wide studies on the association between enhancer and AS (enhancer-AS association). To bridge the gap, here we developed an integrative analysis on a genome-wide scale to identify enhancer-AS associations in human and mouse.

Result

We collected enhancer datasets which include 28 human and 24 mouse tissues and cell lines, and RNA-seq datasets which are paired with the selected tissues. Combining with data integration and statistical analysis, we identified 3,242 human and 7,716 mouse genes which have significant enhancer-AS associations in at least one tissue. On average, for each gene, about 6% of enhancers in human (5% in mouse) are associated to AS change and for each enhancer, approximately one gene is identified to have enhancer-AS association in both human and mouse. We found that 52% of the human significant (34% in mouse) enhancer-AS associations are the co-existence of homologous genes and homologous enhancers. We further constructed a user-friendly platform, named Visualization of Enhancer-associated Alternative Splicing (VEnAS, http://venas.iis.sinica.edu.tw/), to provide genomic architecture, intuitive association plot, and contingency table of the significant enhancer-AS associations.

Conclusion

This study provides the first genome-wide identification of enhancer-AS associations in human and mouse. The results suggest that a notable portion of enhancers are playing roles in AS regulations. The analyzed results and the proposed platform VEnAS would provide a further understanding of enhancers on regulating alternative splicing.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08537-1.

Genome-wide identification, phylogenetic and expression pattern analysis of MADS-box family genes in foxtail millet (Setaria italica)

Foxtail millet (Setaria italica) is rich in nutrients and extremely beneficial to human health. We identified and comprehensively analyzed 89 MADS-box genes in the foxtail millet genome. According to the classification of MADS-box genes in Arabidopsis thaliana and rice, the SiMADS-box genes were divided into M-type (37) and MIKC-type (52). During evolution, the differentiation of MIKC-type MADS-box genes occurred before that of monocotyledons and dicotyledons. The SiMADS-box gene structure has undergone much differentiation, and the number of introns in the MIKC-type subfamily is much greater than that in the M-type subfamily. Analysis of gene duplication events revealed that MIKC-type MADS-box gene segmental duplication accounted for the vast majority of gene duplication events, and MIKC-type MADS-box genes played a major role in the amplification of SiMADS-box genes. Collinearity analysis showed highest collinearity between foxtail millet and maize MADS-box genes. Analysis of tissue-specific expression showed that SiMADS-box genes are highly expressed throughout the grain-filling process. Expression analysis of SiMADS-box genes under eight different abiotic stresses revealed many stress-tolerant genes, with induced expression of SiMADS33 and SiMADS78 under various stresses warranting further attention. Further, some SiMADS-box proteins may interact under external stress. This study provides insights for MADS-box gene mining and molecular breeding of foxtail millet in the future.

Comprehensive deletion landscape of CRISPR-Cas9 identifies minimal RNA-guided DNA-binding modules

Proteins evolve through the modular rearrangement of elements known as domains. Extant, multidomain proteins are hypothesized to be the result of domain accretion, but there has been limited experimental validation of this idea. Here, we introduce a technique for genetic minimization by iterative size-exclusion and recombination (MISER) for comprehensively making all possible deletions of a protein. Using MISER, we generate a deletion landscape for the CRISPR protein Cas9. We find that the catalytically-dead Streptococcus pyogenes Cas9 can tolerate large single deletions in the REC2, REC3, HNH, and RuvC domains, while still functioning in vitro and in vivo, and that these deletions can be stacked together to engineer minimal, DNA-binding effector proteins. In total, our results demonstrate that extant proteins retain significant modularity from the accretion process and, as genetic size is a major limitation for viral delivery systems, establish a general technique to improve genome editing and gene therapy-based therapeutics.

The genome, transcriptome, and proteome of the fish parasite Pomphorhynchus laevis (Acanthocephala)

Thorny-headed worms (Acanthocephala) are endoparasites exploiting Mandibulata (Arthropoda) and Gnathostomata (Vertebrata). Despite their world-wide occurrence and economic relevance as a pest, genome and transcriptome assemblies have not been published before. However, such data might hold clues for a sustainable control of acanthocephalans in animal production. For this reason, we present the first draft of an acanthocephalan nuclear genome, besides the mitochondrial one, using the fish parasite Pomphorhynchus laevis (Palaeacanthocephala) as a model. Additionally, we have assembled and annotated the transcriptome of this species and the proteins encoded. A hybrid assembly of long and short reads resulted in a near-complete P. laevis draft genome of ca. 260 Mb, comprising a large repetitive portion of ca. 63%. Numbers of transcripts and translated proteins (35,683) were within the range of other members of the Rotifera-Acanthocephala clade. Our data additionally demonstrate a significant reorganization of the acanthocephalan gene repertoire. Thus, more than 20% of the usually conserved metazoan genes were lacking in P. laevis. Ontology analysis of the retained genes revealed many connections to the incorporation of carotinoids. These are probably taken up via the surface together with lipids, thus accounting for the orange coloration of P. laevis. Furthermore, we found transcripts and protein sequences to be more derived in P. laevis than in rotifers from Monogononta and Bdelloidea. This was especially the case in genes involved in energy metabolism, which might reflect the acanthocephalan ability to use the scarce oxygen in the host intestine for respiration and simultaneously carry out fermentation. Increased plasticity of the gene repertoire through the integration of foreign DNA into the nuclear genome seems to be another underpinning factor of the evolutionary success of acanthocephalans. In any case, energy-related genes and their proteins may be considered as candidate targets for the acanthocephalan control.

Intrinsic Disorder in Tetratricopeptide Repeat Proteins

Among the realm of repeat containing proteins that commonly serve as “scaffolds” promoting protein-protein interactions, there is a family of proteins containing between 2 and 20 tetratricopeptide repeats (TPRs), which are functional motifs consisting of 34 amino acids. The most distinguishing feature of TPR domains is their ability to stack continuously one upon the other, with these stacked repeats being able to affect interaction with binding partners either sequentially or in combination. It is known that many repeat-containing proteins are characterized by high levels of intrinsic disorder, and that many protein tandem repeats can be intrinsically disordered. Furthermore, it seems that TPR-containing proteins share many characteristics with hybrid proteins containing ordered domains and intrinsically disordered protein regions. However, there has not been a systematic analysis of the intrinsic disorder status of TPR proteins. To fill this gap, we analyzed 166 human TPR proteins to determine the degree to which proteins containing TPR motifs are affected by intrinsic disorder. Our analysis revealed that these proteins are characterized by different levels of intrinsic disorder and contain functional disordered regions that are utilized for protein-protein interactions and often serve as targets of various posttranslational modifications.

Transcriptome analysis of lateral buds from Phyllostachys edulis rhizome during germination and early shoot stages

BACKGROUND

The vegetative growth is an important stage for plants when they conduct photosynthesis, accumulate and collect all resources needed and prepare for reproduction stage. Bamboo is one of the fastest growing plant species. The rapid growth of Phyllostachys edulis results from the expansion of intercalary meristem at the basal part of nodes, which are differentiated from the apical meristem of rhizome lateral buds. However, little is known about the major signaling pathways and players involved during this rapid development stage of bamboo. To study this question, we adopted the high-throughput sequencing technology and compared the transcriptomes of Moso bamboo rhizome buds in germination stage and late development stage.

RESULTS

We found that the development of Moso bamboo rhizome lateral buds was coordinated by multiple pathways, including meristem development, sugar metabolism and phytohormone signaling. Phytohormones have fundamental impacts on the plant development. We found the evidence of several major hormones participating in the development of Moso bamboo rhizome lateral bud. Furthermore, we showed direct evidence that Gibberellic Acids (GA) signaling participated in the Moso bamboo stem elongation.

CONCLUSION

Significant changes occur in various signaling pathways during the development of rhizome lateral buds. It is crucial to understand how these changes are translated to Phyllostachys edulis fast growth. These results expand our knowledge on the Moso bamboo internodes fast growth and provide research basis for further study.

Altered expression of OsNPF7.1 and OsNPF7.4 differentially regulates tillering and grain yield in rice

The rice nitrate and di/tripeptide transporter (NPF) gene family plays an indispensable role in nitrogen transport and plant growth. In this study, 18 alternatively spliced OsNPF genes with 36 different forms of mRNAs were identified, and two of these, namely OsNPF7.1 and OsNPF7.4, showed opposite expression patterns in axillary buds under different nitrogen concentrations. Our results indicate that the expression levels of OsNPF7.1 and OsNPF7.4 determine the axillary bud outgrowth, especially for the second bud, and subsequently influence the tiller number in rice. The overexpression of either of the variants of OsNPF7.1 or the knockout of OsNPF7.4 increased the seedling biomass as well as the tiller number, filled grain number, and grain yield in rice. However, the RNAi-mediated silencing of OsNPF7.1 or the overexpression of either of the variants of OsNPF7.4 had an opposite effect. The overexpression of OsNPF7.1 or OsNPF7.4 could improve the uptake of nitrate, but the OsNPF7.4-overexpressing plants had lower biomass. It is possible that excessive nitrate in the OsNPF7.4-overexpressing plants led to the accumulation of amino acids in the leaf sheath, which inhibited seedling biomass. In addition, the reduced reutilization of nitrate in the seedlings also limited the plant biomass. However, the moderate increase in nitrate and amino acid concentrations in the OsNPF7.1-overexpressing plants could promote seedling biomass and enhance grain yield. In conclusion, our data suggest that different members in the NPF family have different roles in rice, and this study suggests an alternative way to modify rice architecture and enhance grain yield by regulating the expression of OsNPF7.1 and OsNPF7.4.

Novel alternative splice variants of the human protein arginine methyltransferase 1 (PRMT1) gene, discovered using next-generation sequencing

Next-generation sequencing (NGS) technology is highly expected to help researchers disclose the complexity of alternative splicing and understand its association with carcinogenesis. Alternative splicing alterations are firmly associated with multiple malignancies, in terms of functional roles in malignant transformation, motility, and/or metastasis of cancer cells. One perfect example illustrating the connection between alternative splicing and cancer is the human protein arginine methyltransferase 1 (PRMT1) gene, previously cloned from members of our research group and involved in a variety of processes including transcription, DNA repair, and signal transduction. Two splice variants of PRMT1 (variants v.1 and v.2) are downregulated in breast cancer. In addition, PRMT1 v.2 promotes the survival and invasiveness of breast cancer cells, while it could serve as a biomarker of unfavorable prognosis in colon cancer patients. The aim of this study was the molecular cloning of novel alternative splice variants of PRMT1 with the use of 3' RACE coupled with NGS technology. Extensive bioinformatics and computational analysis revealed a significant number of 19 novel alternative splicing events between annotated exons of PRMT1 as well as one novel exon, resulting in the discovery of multiple PRMT1 transcripts. In order to validate the full sequence of the novel transcripts, RT-PCR was carried out with the use of variant-specific primers. As a result, 58 novel PRMT1 transcripts were identified, 34 of which are mRNAs encoding new protein isoforms, whereas the rest 24 transcripts are candidates for nonsense-mediated mRNA decay (NMD).

Developmental processes regulate craniofacial variation in disease and evolution

Variation in development mediates phenotypic differences observed in evolution and disease. Although the mechanisms underlying phenotypic variation are still largely unknown, recent research suggests that variation in developmental processes may play a key role. Developmental processes mediate genotype-phenotype relationships and consequently play an important role regulating phenotypes. In this review, we provide an example of how shared and interacting developmental processes may explain convergence of phenotypes in spliceosomopathies and ribosomopathies. These data also suggest a shared pathway to disease treatment. We then discuss three major mechanisms that contribute to variation in developmental processes: genetic background (gene-gene interactions), gene-environment interactions, and developmental stochasticity. Finally, we comment on evolutionary alterations to developmental processes, and the evolution of disease buffering mechanisms.

The Putative Smallest Introns in the Arabidopsis Genome

Most eukaryotic genes contain introns, which are noncoding sequences that are removed during premRNA processing. Introns are usually preserved across evolutionary time. However, the sizes of introns vary greatly. In Arabidopsis, some introns are longer than 10 kilo base pairs (bp) and others are predicted to be shorter than 10 bp. To identify the shortest intron in the genome, we analyzed the predicted introns in annotated version 10 of the Arabidopsis thaliana genome and found 103 predicted introns that are 30 bp or shorter, which make up only 0.08% of all introns in the genome. However, our own bioinformatics and experimental analyses found no evidence for the existence of these predicted introns. The predicted introns of 30–39 bp, 40–49 bp, and 50–59 bp in length are also rare and constitute only 0.07%, 0.2%, and 0.28% of all introns in the genome, respectively. An analysis of 30 predicted introns 31–59 bp long verified two in this range, both of which were 59 bp long. Thus, this study suggests that there is a limit to how small introns in A. thaliana can be, which is useful for the understanding of the evolution and processing of small introns in plants in general.

Systematic evaluation of isoform function in literature reports of alternative splicing

BACKGROUND

Although most genes in mammalian genomes have multiple isoforms, an ongoing debate is whether these isoforms are all functional as well as the extent to which they increase the functional repertoire of the genome. To ground this debate in data, it would be helpful to have a corpus of experimentally-verified cases of genes which have functionally distinct splice isoforms (FDSIs).

RESULTS

We established a curation framework for evaluating experimental evidence of FDSIs, and analyzed over 700 human and mouse genes, strongly biased towards genes that are prominent in the alternative splicing literature. Despite this bias, we found experimental evidence meeting the classical definition for functionally distinct isoforms for ~ 5% of the curated genes. If we relax our criteria for inclusion to include weaker forms of evidence, the fraction of genes with evidence of FDSIs remains low (~ 13%). We provide evidence that this picture will not change substantially with further curation and conclude there is a large gap between the presumed impact of splicing on gene function and the experimental evidence. Furthermore, many functionally distinct isoforms were not traceable to a specific isoform in Ensembl, a database that forms the basis for much computational research.

CONCLUSIONS

We conclude that the claim that alternative splicing vastly increases the functional repertoire of the genome is an extrapolation from a limited number of empirically supported cases. We also conclude that more work is needed to integrate experimental evidence and genome annotation databases. Our work should help shape research around the role of splicing on gene function from presuming large general effects to acknowledging the need for stronger experimental evidence.

Two Splicing Variants of OsNPF7.7 Regulate Shoot Branching and Nitrogen Utilization Efficiency in Rice

Rice includes 93 nitrate and peptide transporters family (NPF) members that facilitate the soil uptake and internal reallocation of nitrogen for growth and development. This study demonstrated that OsNPF7.7 had two splicing variants, and altered expression of each variant could regulate shoot branching and nitrogen utilization efficiency (NUtE) in rice. The expression of both variants was down-regulated in the buds by increased nitrogen level in the Japonica rice variety ZH11. The expression level of long-variant OsNPF7.7-1 was higher in panicles at reproductive stage, however, the expression level of short-variant OsNPF7.7-2 was higher in buds and leaves at vegetative stage compared to each other in ZH11. OsNPF7.7-1 was localized in the plasma membrane, whereas OsNPF7.7-2 was localized in the vacuole membrane. Furthermore, the results indicated that the expression level of each variant for OsNPF7.7 determined axillary bud outgrowth, and then influenced the rice tiller number. Overexpression of OsNPF7.7-1 could promote nitrate influx and concentration in root, whereas overexpression of OsNPF7.7-2 could improve ammonium influx and concentration in root. RNAi and osnpf7.7 lines of OsNPF7.7 showed an increased amount of amino acids in leaf sheaths, but a decreased amount in leaf blades, which affected nitrogen allocation and plant growth. The elevated expression of each variant for OsNPF7.7 in ZH11 enhanced NUtE using certain fertilization regimes under paddy field conditions. Moreover, overexpression of each variant for OsNPF7.7 in KY131 increased significantly the filled grain number per plant. Thus, increased each variant of OsNPF7.7 has the potential to improve grain yield and NUtE in rice.

Conserved molecular structure of the centromeric histone CENH3 in Secale and its phylogenetic relationships

It has been repeatedly demonstrated that the centromere-specific histone H3 (CENH3), a key component of the centromere, shows considerable variability between species within taxa. We determined the molecular structure and phylogenetic relationships of CENH3 in 11 Secale species and subspecies that possess distinct pollination systems and are adapted to a wide range of abiotic and biotic stresses. The rye (Secale cereale) genome encodes two paralogous CENH3 genes, which differ in intron-exon structure and are transcribed into two main forms of the protein, αCENH3 and βCENH3. These two forms differ in size and amino acid substitutions. In contrast to the reported differences in CENH3 structure between species within other taxa, the main forms of this protein in Secale species and subspecies have a nearly identical structure except some nonsynonymous substitutions. The CENH3 proteins are strictly controlled by genetic factors responsible for purifying selection. A comparison between Hordeum, Secale and Triticum species demonstrates that the structure of CENH3 in the subtribes Hordeinae and Triticinae evolved at different rates. The assumption that reticulate evolution served as a factor stabilizing the structure and evolutionary rate of CENH3 and that this factor was more powerful within Secale and Triticum than in Hordeum, is discussed.

Intron retention enhances gene regulatory complexity in vertebrates

Background

While intron retention (IR) is now widely accepted as an important mechanism of mammalian gene expression control, it remains the least studied form of alternative splicing. To delineate conserved features of IR, we performed an exhaustive phylogenetic analysis in a highly purified and functionally defined cell type comprising neutrophilic granulocytes from five vertebrate species spanning 430 million years of evolution.

Results

Our RNA-sequencing-based analysis suggests that IR increases gene regulatory complexity, which is indicated by a strong anti-correlation between the number of genes affected by IR and the number of protein-coding genes in the genome of individual species. Our results confirm that IR affects many orthologous or functionally related genes in granulocytes. Further analysis uncovers new and unanticipated conserved characteristics of intron-retaining transcripts. We find that intron-retaining genes are transcriptionally co-regulated from bidirectional promoters. Intron-retaining genes have significantly longer 3′ UTR sequences, with a corresponding increase in microRNA binding sites, some of which include highly conserved sequence motifs. This suggests that intron-retaining genes are highly regulated post-transcriptionally.

Conclusions

Our study provides unique insights concerning the role of IR as a robust and evolutionarily conserved mechanism of gene expression regulation. Our findings enhance our understanding of gene regulatory complexity by adding another contributor to evolutionary adaptation.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-017-1339-3) contains supplementary material, which is available to authorized users.

Analysis of Genes with Alternatively Spliced Transcripts in the Leaf, Root, Panicle and Seed of Rice Using a Long Oligomer Microarray and RNA-Seq

Pre-mRNA splicing further increases protein diversity acquired through evolution. The underlying driving forces for this phenomenon are unknown, especially in terms of gene expression. A rice alternatively spliced transcript detection microarray (ASDM) and RNA sequencing (RNA-Seq) were applied to differentiate the transcriptome of 4 representative organs of Oryza sativa L. cv. Ilmi: leaves, roots, 1-cm-stage panicles and young seeds at 21 days after pollination. Comparison of data obtained by microarray and RNA-Seq showed a bell-shaped distribution and a co-lineation for highly expressed genes. Transcripts were classified according to the degree of organ enrichment using a coefficient value (CV, the ratio of the standard deviation to the mean values): highly variable (CVI), variable (CVII), and constitutive (CVIII) groups. A higher index of the portion of loci with alternatively spliced transcripts in a group (IAST) value was observed for the constitutive group. Genes of the highly variable group showed the characteristics of the examined organs, and alternatively spliced transcripts tended to exhibit the same organ specificity or less organ preferences, with avoidance of 'organ distinctness'. In addition, within a locus, a tendency of higher expression was found for transcripts with a longer coding sequence (CDS), and a spliced intron was the most commonly found type of alternative splicing for an extended CDS. Thus, pre-mRNA splicing might have evolved to retain maximum functionality in terms of organ preference and multiplicity.

Mis-Spliced Lr34 Transcript Events in Winter Wheat

Lr34 in wheat is a non-race-specific gene that confers resistance against multiple fungal pathogens. The resistant allele Lr34 and the susceptible allele Lr34s can be distinguished by three polymorphisms that cause alternation of deduced amino acid sequences of Lr34 at the protein level. In seedlings of a cultivar carrying the resistant Lr34r allele, only a portion (35%) of its transcripts was correctly spliced and the majority (65%) of its transcripts were incorrectly spliced due to multiple mis-splicing events. Lr34 mis-splicing events were also observed at adult plant age when this gene exerts its function. All of the mis-spliced Lr34r cDNA transcripts observed in this study resulted in a premature stop codon due to a shift of the open reading frame; hence, the mis-spliced Lr34r cDNAs were deduced to encode incomplete proteins. Even if a cultivar has a functional Lr34 gene, its transcripts might not completely splice in a correct pattern. These findings suggested that the partial resistance conferred by a quantitative gene might be due to mis-splicing events in its transcripts; hence, the resistance of the gene could be increased by eliminating or mutating regulators that cause mis-splicing events in wheat.

Fluorescent Protein-Based Quantification of Alternative Splicing of a Target Cassette Exon in Mammalian Cells

Alternative splicing is an important mechanism of regulation of gene expression and expansion of proteome complexity. Recently we developed a new fluorescence reporter for quantitative analysis of alternative splicing of a target cassette exon in live cells (Gurskaya et al., 2012). It consists of a specially designed minigene encoding red and green fluorescent proteins (Katushka and TagGFP2) and a fragment of the target gene between them. Skipping or inclusion of the alternative exon induces a frameshift; ie, alternative exon length must not be a multiple of 3. Finally, red and green fluorescence intensities of cells expressing this reporter are used to estimate the percentage of alternative (exon-skipped) and normal (exon-retained) transcripts. Here, we provide a detailed description of design and application of the fluorescence reporter of a target alternative exon splicing in mammalian cell lines.

A proteogenomic approach to understand splice isoform functions through sequence and expression-based computational modeling

The products of multi-exon genes are a mixture of alternatively spliced isoforms, from which the translated proteins can have similar, different or even opposing functions. It is therefore essential to differentiate and annotate functions for individual isoforms. Computational approaches provide an efficient complement to expensive and time-consuming experimental studies. The input data of these methods range from DNA sequence, to RNA selection pressure, to expressed sequence tags, to full-length complementary DNA, to exon array, to RNA-seq expression, to proteomic data. Notably, RNA-seq technology generates quantitative profiling of transcript expression at the genome scale, with an unprecedented amount of expression data available for developing isoform function prediction methods. Integrative analysis of these data at different molecular levels enables a proteogenomic approach to systematically interrogate isoform functions. Here, we briefly review the state-of-the-art methods according to their input data sources, discuss their advantages and limitations and point out potential ways to improve prediction accuracies.

Network-Based Isoform Quantification with RNA-Seq Data for Cancer Transcriptome Analysis

High-throughput mRNA sequencing (RNA-Seq) is widely used for transcript quantification of gene isoforms. Since RNA-Seq data alone is often not sufficient to accurately identify the read origins from the isoforms for quantification, we propose to explore protein domain-domain interactions as prior knowledge for integrative analysis with RNA-Seq data. We introduce a Network-based method for RNA-Seq-based Transcript Quantification (Net-RSTQ) to integrate protein domain-domain interaction network with short read alignments for transcript abundance estimation. Based on our observation that the abundances of the neighboring isoforms by domain-domain interactions in the network are positively correlated, Net-RSTQ models the expression of the neighboring transcripts as Dirichlet priors on the likelihood of the observed read alignments against the transcripts in one gene. The transcript abundances of all the genes are then jointly estimated with alternating optimization of multiple EM problems. In simulation Net-RSTQ effectively improved isoform transcript quantifications when isoform co-expressions correlate with their interactions. qRT-PCR results on 25 multi-isoform genes in a stem cell line, an ovarian cancer cell line, and a breast cancer cell line also showed that Net-RSTQ estimated more consistent isoform proportions with RNA-Seq data. In the experiments on the RNA-Seq data in The Cancer Genome Atlas (TCGA), the transcript abundances estimated by Net-RSTQ are more informative for patient sample classification of ovarian cancer, breast cancer and lung cancer. All experimental results collectively support that Net-RSTQ is a promising approach for isoform quantification. Net-RSTQ toolbox is available at http://compbio.cs.umn.edu/Net-RSTQ/.

Domain atrophy creates rare cases of functional partial protein domains

BACKGROUND

Protein domains display a range of structural diversity, with numerous additions and deletions of secondary structural elements between related domains. We have observed a small number of cases of surprising large-scale deletions of core elements of structural domains. We propose a new concept called domain atrophy, where protein domains lose a significant number of core structural elements.

RESULTS

Here, we implement a new pipeline to systematically identify new cases of domain atrophy across all known protein sequences. The output of this pipeline was carefully checked by hand, which filtered out partial domain instances that were unlikely to represent true domain atrophy due to misannotations or un-annotated sequence fragments. We identify 75 cases of domain atrophy, of which eight cases are found in a three-dimensional protein structure and 67 cases have been inferred based on mapping to a known homologous structure. Domains with structural variations include ancient folds such as the TIM-barrel and Rossmann folds. Most of these domains are observed to show structural loss that does not affect their functional sites.

CONCLUSION

Our analysis has significantly increased the known cases of domain atrophy. We discuss specific instances of domain atrophy and see that there has often been a compensatory mechanism that helps to maintain the stability of the partial domain. Our study indicates that although domain atrophy is an extremely rare phenomenon, protein domains under certain circumstances can tolerate extreme mutations giving rise to partial, but functional, domains.

Comprehensive analysis of alternative splicing in Digitalis purpurea by strand-specific RNA-Seq

Digitalis purpurea (D. purpurea) is one of the most important medicinal plants and is well known in the treatment of heart failure because of the cardiac glycosides that are its main active compounds. However, in the absence of strand specific sequencing information, the post-transcriptional mechanism of gene regulation in D. purpurea thus far remains unknown. In this study, a strand-specific RNA-Seq library was constructed and sequenced using Illumina HiSeq platforms to characterize the transcriptome of D. purpurea with a focus on alternative splicing (AS) events and the effect of AS on protein domains. De novo RNA-Seq assembly resulted in 48,475 genes. Based on the assembled transcripts, we reported a list of 3,265 AS genes, including 5,408 AS events in D. purpurea. Interestingly, both glycosyltransferases and monooxygenase, which were involved in the biosynthesis of cardiac glycosides, are regulated by AS. A total of 2,422 AS events occurred in coding regions, and 959 AS events were located in the regions of 882 unique protein domains, which could affect protein function. This D. purpurea transcriptome study substantially increased the expressed sequence resource and presented a better understanding of post-transcriptional regulation to further facilitate the medicinal applications of D. purpurea for human health.

Gene duplication followed by exon structure divergence substitutes for alternative splicing in zebrafish

In this study we report novel findings regarding the evolutionary relationship between gene duplication and alternative splicing, two processes that increase proteomic diversity. By studying teleost fish, we find that gene duplication followed by exon structure divergence between paralogs, but not gene duplication alone, leads to a significant reduction in alternative splicing, as measured by both the proportion of genes that undergo alternative splicing as well as mean number of transcripts per gene. Additionally, we show that this effect is independent of gene family size and gene function. Furthermore, we provide evidence that the reduction in alternative splicing may be due to the partitioning of ancestral splice forms among the duplicate genes - a form of subfunctionalization. Taken together these results indicate that exon structure evolution subsequent to gene duplication may be a common substitute for alternative splicing.

Evolution at protein ends: major contribution of alternative transcription initiation and termination to the transcriptome and proteome diversity in mammals

Alternative splicing (AS), alternative transcription initiation (ATI) and alternative transcription termination (ATT) create the extraordinary complexity of transcriptomes and make key contributions to the structural and functional diversity of mammalian proteomes. Analysis of mammalian genomic and transcriptomic data shows that contrary to the traditional view, the joint contribution of ATI and ATT to the transcriptome and proteome diversity is quantitatively greater than the contribution of AS. Although the mean numbers of protein-coding constitutive and alternative nucleotides in gene loci are nearly identical, their distribution along the transcripts is highly non-uniform. On average, coding exons in the variable 5' and 3' transcript ends that are created by ATI and ATT contain approximately four times more alternative nucleotides than core protein-coding regions that diversify exclusively via AS. Short upstream exons that encompass alternative 5'-untranslated regions and N-termini of proteins evolve under strong nucleotide-level selection whereas in 3'-terminal exons that encode protein C-termini, protein-level selection is significantly stronger. The groups of genes that are subject to ATI and ATT show major differences in biological roles, expression and selection patterns.

Correcting for differential transcript coverage reveals a strong relationship between alternative splicing and organism complexity

What at the genomic level underlies organism complexity? Although several genomic features have been associated with organism complexity, in the case of alternative splicing, which has long been proposed to explain the variation in complexity, no such link has been established. Here, we analyzed over 39 million expressed sequence tags available for 47 eukaryotic species with fully sequenced genomes to obtain a comparable index of alternative splicing estimates, which corrects for the distorting effect of a variable number of transcripts per species—an important obstacle for comparative studies of alternative splicing. We find that alternative splicing has steadily increased over the last 1,400 My of eukaryotic evolution and is strongly associated with organism complexity, assayed as the number of cell types. Importantly, this association is not explained as a by-product of covariance between alternative splicing with other variables previously linked to complexity including gene content, protein length, proteome disorder, and protein interactivity. In addition, we found no evidence to suggest that the relationship of alternative splicing to cell type number is explained by drift due to reduced N_e in more complex species. Taken together, our results firmly establish alternative splicing as a significant predictor of organism complexity and are, in principle, consistent with an important role of transcript diversification through alternative splicing as a means of determining a genome’s functional information capacity.

Analysis of the overall structure of the multi-domain amyloid precursor protein (APP)

The amyloid precursor protein (APP) and its processing by the α-, β- and γ-secretases is widely believed to play a central role during the development of Alzheimer´s disease. The three-dimensional structure of the entire protein, its physiologic function and the regulation of its proteolytic processing remain, however, largely unclear to date. To gain a deeper understanding of the structure of APP that underlies all of its functions, we first cloned and recombinantly expressed different constructs in E. coli. Using limited proteolysis followed by mass spectrometry and Edman degradation as well as analytical gel permeation chromatography coupled static light scattering, we experimentally analyzed the structural domain boundaries and determined that the large ectodomain of APP consists of exactly two rigidly folded domains - the E1-domain (Leu18-Ala190) and the E2-domain (Ser295-Asp500). Both, the acidic domain (AcD) connecting E1 and E2 as well as the juxtamembrane region (JMR) connecting E2 to the single transmembrane helix are highly flexible and extended. We identified in-between the E1-domain and the AcD an additional domain of conservation and partial flexibility that we termed extension domain (ED, Glu191-Glu227). Using Bio-layer interferometry, pull-down assays and analytical gel filtration experiments we demonstrated that the E1-domain does not tightly interact with the E2-domain, both in the presence and in the absence of heparin. APP hence forms an extended molecule that is flexibly tethered to the membrane. Its multi-domain architecture enables together with the many known functionalities the concomitant performance of several, independent functions, which might be regulated by cellular, compartment specific pH-changes.

Protein isoform-specific validation defines multiple chloride intracellular channel and tropomyosin isoforms as serological biomarkers of ovarian cancer

New serological biomarkers for early detection and clinical management of ovarian cancer are urgently needed, and many candidates have been reported. A major challenge frequently encountered when validating candidates in patients is establishing quantitative assays that distinguish between highly homologous proteins. The current study tested whether multiple members of two recently discovered ovarian cancer biomarker protein families, chloride intracellular channel (CLIC) proteins and tropomyosins (TPM), were detectable in ovarian cancer patient sera. A multiplexed, label-free multiple reaction monitoring (MRM) assay was established to target peptides specific to all detected CLIC and TPM family members, and their serum levels were quantitated for ovarian cancer patients and non-cancer controls. In addition to CLIC1 and TPM1, which were the proteins initially discovered in a xenograft mouse model, CLIC4, TPM2, TPM3, and TPM4 were present in ovarian cancer patient sera at significantly elevated levels compared with controls. Some of the additional biomarkers identified in this homolog-centric verification and validation approach may be superior to the previously identified biomarkers at discriminating between ovarian cancer and non-cancer patients. This demonstrates the importance of considering all potential protein homologs and using quantitative assays for cancer biomarker validation with well-defined isoform specificity.

BIOLOGICAL SIGNIFICANCE

This manuscript addresses the importance of distinguishing between protein homologs and isoforms when identifying and validating cancer biomarkers in plasma or serum. Specifically, it describes the use of targeted in-depth LC-MS/MS analysis to determine the members of two protein families, chloride intracellular channel (CLIC) and tropomyosin (TPM) proteins that are detectable in sera of ovarian cancer patients. It then establishes a multiplexed isoform- and homology-specific MRM assay to quantify all observed gene products in these two protein families as well as many of the closely related tropomyosin isoforms. Using this assay, levels of all detected CLICs and TPMs were quantified in ovarian cancer patient and control subject sera. These results demonstrate that in addition to the previously known CLIC1, multiple tropomyosins and CLIC4 are promising new ovarian cancer biomarkers. Based on these initial validation studies, these new ovarian cancer biomarkers appear to be superior to most previously known ovarian cancer biomarkers.

Distinct types of disorder in the human proteome: functional implications for alternative splicing

Intrinsically disordered regions have been associated with various cellular processes and are implicated in several human diseases, but their exact roles remain unclear. We previously defined two classes of conserved disordered regions in budding yeast, referred to as "flexible" and "constrained" conserved disorder. In flexible disorder, the property of disorder has been positionally conserved during evolution, whereas in constrained disorder, both the amino acid sequence and the property of disorder have been conserved. Here, we show that flexible and constrained disorder are widespread in the human proteome, and are particularly common in proteins with regulatory functions. Both classes of disordered sequences are highly enriched in regions of proteins that undergo tissue-specific (TS) alternative splicing (AS), but not in regions of proteins that undergo general (i.e., not tissue-regulated) AS. Flexible disorder is more highly enriched in TS alternative exons, whereas constrained disorder is more highly enriched in exons that flank TS alternative exons. These latter regions are also significantly more enriched in potential phosphosites and other short linear motifs associated with cell signaling. We further show that cancer driver mutations are significantly enriched in regions of proteins associated with TS and general AS. Collectively, our results point to distinct roles for TS alternative exons and flanking exons in the dynamic regulation of protein interaction networks in response to signaling activity, and they further suggest that alternatively spliced regions of proteins are often functionally altered by mutations responsible for cancer.

The impact of splicing on protein domain architecture

Many proteins are composed of protein domains, functional units of common descent. Multidomain forms are common in all eukaryotes making up more than half of the proteome and the evolution of novel domain architecture has been accelerated in metazoans. It is also becoming increasingly clear that alternative splicing is prevalent among vertebrates. Given that protein domains are defined as structurally, functionally and evolutionarily distinct units, one may speculate that some alternative splicing events may lead to clean excisions of protein domains, thus generating a number of different domain architectures from one gene template. However, recent findings indicate that smaller alternative splicing events, in particular in disordered regions, might be more prominent than domain architectural changes. The problem of identifying protein isoforms is, however, still not resolved. Clearly, many splice forms identified through detection of mRNA sequences appear to produce 'nonfunctional' proteins, such as proteins with missing internal secondary structure elements. Here, we review the state of the art methods for identification of functional isoforms and present a summary of what is known, thus far, about alternative splicing with regard to protein domain architectures.

Molecular cloning and expression profile of sex-specific genes, Figla and Dmrt1, in the protogynous hermaphroditic fish, Halichoeres poecilopterus

The genes folliculogenesis specific basic helix-loop-helix (facor in the germline alpha, Figla) and doublesex and mab-3 related transcription factor 1 (Dmrt1) are female- and male-specific genes that play key roles in sex differentiation. To obtain a better understanding of the molecular mechanisms underlying female-to-male sex change, we cloned the cDNAs of these genes from an ovary and a testis of the protogynus wrasse, Halichoeres poecilopterus. This fish has two isoforms of Dmrt1, Dmrt1a and Dmrt1b, caused by an alternative splicing. The Dmrt1b has an insertion of three nucleotides (CAG) in the open reading frame. Figla and Dmrt1 displayed gonadal-specific expression and abundant in the ovaries and in the testes, respectively. In particular, levels of Figla expression in the ovaries were higher in the spawning season than in the non-spawning season. Once sex change began, Figla mRNA decreased and Dmrt1 mRNA increased with progression of oocyte degeneration and spermatogenesis. These expression levels were maintained until the completion of the sex change. Low Figla and high Dmrt1 were also observed in testes of primary males, which functioned as a gonochoristic male throughout its life span in this wrasse. The results of this study suggest that these genes may regulate the gonadal transition from ovary to testis by the same mechanism as that of formation and maintenance of the primary testis in H. poecilopterus.

Whence genes in pieces: reconstruction of the exon-intron gene structures of the last eukaryotic common ancestor and other ancestral eukaryotes

In eukaryotes, protein-coding sequences are interrupted by non-coding sequences known as introns. During mRNA maturation, introns are excised by the spliceosome and the coding regions, exons, are spliced to form the mature coding region. The intron densities widely differ between eukaryotic lineages, from 6 to 7 introns per kb of coding sequence in vertebrates, some invertebrates and green plants, to only a few introns across the entire genome in many unicellular eukaryotes. Evolutionary reconstructions using maximum likelihood methods suggest intron-rich ancestors for each major group of eukaryotes. For the last common ancestor of animals, the highest intron density of all extant and extinct eukaryotes was inferred, at 120-130% of the human intron density. Furthermore, an intron density within 53-74% of the human values was inferred for the last eukaryotic common ancestor. Accordingly, evolution of eukaryotic genes in all lines of descent involved primarily intron loss, with substantial gain only at the bases of several branches including plants and animals. These conclusions have substantial biological implications indicating that the common ancestor of all modern eukaryotes was a complex organism with a gene architecture resembling those in multicellular organisms. Alternative splicing most likely initially appeared as an inevitable result of splicing errors and only later was employed to generate structural and functional diversification of proteins.

Tissue-specific splicing of disordered segments that embed binding motifs rewires protein interaction networks

Alternative inclusion of exons increases the functional diversity of proteins. Among alternatively spliced exons, tissue-specific exons play a critical role in maintaining tissue identity. This raises the question of how tissue-specific protein-coding exons influence protein function. Here we investigate the structural, functional, interaction, and evolutionary properties of constitutive, tissue-specific, and other alternative exons in human. We find that tissue-specific protein segments often contain disordered regions, are enriched in posttranslational modification sites, and frequently embed conserved binding motifs. Furthermore, genes containing tissue-specific exons tend to occupy central positions in interaction networks and display distinct interaction partners in the respective tissues, and are enriched in signaling, development, and disease genes. Based on these findings, we propose that tissue-specific inclusion of disordered segments that contain binding motifs rewires interaction networks and signaling pathways. In this way, tissue-specific splicing may contribute to functional versatility of proteins and increases the diversity of interaction networks across tissues.

Identification of novel splice variants and exons of human endothelial cell-specific chemotaxic regulator (ECSCR) by bioinformatics analysis

Recent discovery of biological function of endothelial cell-specific chemotaxic regulator (ECSCR), previously known as endothelial cell-specific molecule 2 (ECSM2), in modulating endothelial cell migration, apoptosis, and angiogenesis, has made it an attractive molecule in vascular research. Thus, identification of splice variants of ECSCR could provide new strategies for better understanding its roles in health and disease. In this study, we performed a series of blast searches on the human EST database with known ECSCR cDNA sequence (Variant 1), and identified additional three splice variants (Variants 2-4). When examining the ECSCR gene in the human genome assemblies, we found a large unknown region between Exons 9 and 11. By PCR amplification and sequencing, we partially mapped Exon 10 within this previously unknown region of the ECSCR gene. Taken together, in addition to previously reported human ECSCR, we identified three novel full-length splice variants potentially encoding different protein isoforms. We further defined a total of twelve exons and nearly all exon-intron boundaries of the gene, of which only eight are annotated in current public databases. Our work provides new information on gene structure and alternative splicing of the human ECSCR, which may imply its functional complexity. This undoubtedly opens new opportunities for future investigation of the biological and pathological significance of these ECSCR splice variants.