BackgroundAt least 90% of human genes are alternatively spliced. Alternative splicing has an important function regulating gene expression and miss-splicing can contribute to risk for human diseases, including Alzheimer’s disease (AD).MethodsWe developed a splicing decision model as a molecular mechanism to identify functional exon skipping events and genetic variation affecting alternative splicing on a genome-wide scale by integrating genomics, transcriptomics, and neuroimaging data in a systems biology approach. In this study, we analyzed RNA-Seq data of hippocampus brain tissue from Alzheimer’s disease (AD; n = 24) and cognitively normal elderly controls (CN; n = 50) and identified three exon skipping events in two genes (RELN and NOS1) as significantly associated with AD (corrected p-value < 0.05 and fold change > 1.5). Next, we identified single-nucleotide polymorphisms (SNPs) affecting exon skipping events using the splicing decision model and then performed an association analysis of SNPs potentially affecting three exon skipping events with a global cortical measure of amyloid-β deposition measured by [18F] Florbetapir position emission tomography (PET) scan as an AD-related quantitative phenotype. A whole-brain voxel-based analysis was also performed.ResultsTwo exons in RELN and one exon in NOS1 showed significantly lower expression levels in the AD participants compared to CN participants, suggesting that the exons tend to be skipped more in AD. We also showed the loss of the core protein structure due to the skipped exons using the protein 3D structure analysis. The targeted SNP-based association analysis identified one intronic SNP (rs362771) adjacent to the skipped exon 24 in RELN as significantly associated with cortical amyloid-β levels (corrected p-value < 0.05). This SNP is within the splicing regulatory element, i.e., intronic splicing enhancer. The minor allele of rs362771 conferred decreases in cortical amyloid-β levels in the right temporal and bilateral parietal lobes.ConclusionsOur results suggest that exon skipping events and splicing-affecting SNPs in the human hippocampus may contribute to AD pathogenesis. Integration of multiple omics and neuroimaging data provides insights into possible mechanisms underlying AD pathophysiology through exon skipping and may help identify novel therapeutic targets.Electronic supplementary materialThe online version of this article (10.1186/s12920-018-0453-8) contains supplementary material, which is available to authorized users.

Background: Elucidating molecular mechanisms that are altered during HIV-1 infection may provide a better understanding of the HIV-1 life cycle and how it interacts with infected T-cells. One such mechanism is alternative splicing (AS), which has been studied for HIV-1 itself, but no systematic analysis has yet been performed on infected T-cells. We hypothesized that AS patterns in infected T-cells may illuminate the molecular mechanisms underlying HIV-1 infection and identify candidate molecular markers for specifically targeting infected T-cells. Methods: We downloaded previously published raw RNA-seq data obtained from HIV-1 infected and non-infected T-cells. We estimated percent spliced in (PSI) levels for each AS exon, then identified differential AS events in the infected cells (FDR < 0.05, PSI difference > 0.1). We performed functional gene set enrichment analysis on the genes with differentially expressed AS exons to identify their functional roles. In addition, we used RT-PCR to validate differential alternative splicing events in cyclin T1 (CCNT1) as a case study. Results: We identified 427 candidate genes with differentially expressed AS exons in infected T-cells, including 20 genes related to cell surface, 35 to kinases, and 121 to immune-related genes. In addition, protein-protein interaction analysis identified six essential subnetworks related to the viral life cycle, including Transcriptional regulation by TP53, Class I MHC mediated antigen, G2/M transition, and late phase of HIV life cycle. CCNT1 exon 7 was more frequently skipped in infected T-cells, leading to loss of the key Cyclin_N motif and affecting HIV-1 transcriptional elongation. Conclusions: Our findings may provide new insight into systemic host AS regulation under HIV-1 infection and may provide useful initial candidates for the discovery of new markers for specifically targeting infected T-cells.

BackgroundHepatitis B virus (HBV), hepatitis C virus (HCV), and alcohol consumption are predominant causes of hepatocellular carcinoma (HCC). However, the molecular mechanisms underlying how differently these causes are implicated in HCC development are not fully understood. Therefore, we investigated differential alternative splicing (AS) regulation among HCC patients with these risk factors.MethodsWe conducted a genome-wide survey of AS events associated with HCCs among HBV (n = 95), HCV (n = 47), or alcohol (n = 76) using RNA-sequencing data obtained from The Cancer Genome Atlas.ResultsIn three group comparisons of HBV vs. HCV, HBV vs. alcohol, and HCV vs. alcohol for RNA seq (ΔPSI> 0.05, FDR < 0.05), 133, 93, and 29 differential AS events (143 genes) were identified, respectively. Of 143 AS genes, eight and one gene were alternatively spliced specific to HBV and HCV, respectively. Through functional analysis over the canonical pathways and gene ontologies, we identified significantly enriched pathways in 143 AS genes including immune system, mRNA splicing-major pathway, and nonsense-mediated decay, which may be important to carcinogenesis in HCC risk factors. Among eight genes with HBV-specific splicing events, HLA-A, HLA-C, and IP6K2 exhibited more differential expression of AS events (ΔPSI> 0.1). Intron retention of HLA-A was observed more frequently in HBV-associated HCC than HCV- or alcohol-associated HCC, and intron retention of HLA-C showed vice versa. Exon 3 (based on ENST00000432678) of IP6K2 was less skipped in HBV-associated in HCC compared to HCV- or alcohol-associated HCC.ConclusionAS may play an important role in regulating transcription differences implicated in HBV-, HCV-, and alcohol-related HCC development.

Gene duplication is widely accepted as a key evolutionary process, leading to new genes and novel protein functions. By providing the raw genetic material necessary for functional expansion, the mechanisms that involve the retention and functional diversification of duplicate genes are one of the central topics in evolutionary and comparative genomics. One proposed source of retention and functional diversification is protein subcellular relocalization (PSR). PSR postulates that changes in the subcellular location of eukaryotic duplicate proteins can positively modify function and therefore be beneficial to the organism. As such, PSR would promote retention of those relocalized duplicates and result in significantly lower death rates compared with death rates of nonrelocalized duplicate pairs. We surveyed both relocalized and nonrelocalized duplicate proteins from the available genomes and proteomes of 59 eukaryotic species and compared their relative death rates over a Ks range between 0 and 1. Using the Cox proportional hazard model, we observed that the death rates of relocalized duplicate pairs were significantly lower than the death rates of the duplicates without relocalization in most eukaryotic species examined in this study. These observations suggest that PSR significantly increases retention of duplicate genes and that it plays an important, but currently underappreciated, role in the evolution of eukaryotic genomes.

Neuroendocrine tumors (NETs) of the small intestine undergo large chromosomal and methylation changes. The objective of this study was to identify methylation differences in NETs and consider how the differentially methylated genes may impact patient survival. Genome-wide methylation and chromosomal copy number variation (CNV) of NETs from the small intestine and appendix were measured. Tumors were divided into three molecular subtypes according to CNV results: chromosome 18 loss (18LOH), Multiple CNV, and No CNV. Comparison of 18LOH tumors with MultiCNV and NoCNV tumors identified 901 differentially methylated genes. Genes from the G-protein coupled receptor (GPCR) pathways are statistically overrepresented in the differentially methylated genes. One of the highlighted genes from the GPCR pathway is somatostatin (SST), a clinical target for NETs. Patient survival based on low versus high methylation in all samples identified four significant genes (p < 0.05) OR2S2, SMILR, RNU6-653P, and AC010543.1. Within the 18LOH molecular subtype tumors, survival differences were identified in high versus low methylation of 24 genes. The most significant is TRHR (p < 0.01), a GPCR with multiple FDA-approved drugs. By separating NETs into different molecular subtypes based on chromosomal changes, we find that multiple GPCRs and their ligands appear to be regulated through methylation and correlated with survival. These results suggest opportunities for better treatment strategies for NETs based on molecular features.