Alternative splicing modulated by genetic variants demonstrates accelerated evolution regulated by highly conserved proteins

Mapping genetic variants for nonsense-mediated mRNA decay regulation across human tissues

BACKGROUND

Nonsense-mediated mRNA decay (NMD) was originally conceived as an mRNA surveillance mechanism to prevent the production of potentially deleterious truncated proteins. Research also shows NMD is an important post-transcriptional gene regulation mechanism selectively targeting many non-aberrant mRNAs. However, how natural genetic variants affect NMD and modulate gene expression remains elusive.

RESULTS

Here we elucidate NMD regulation of individual genes across human tissues through genetical genomics. Genetic variants corresponding to NMD regulation are identified based on GTEx data through unique and robust transcript expression modeling. We identify genetic variants that influence the percentage of NMD-targeted transcripts (pNMD-QTLs), as well as genetic variants regulating the decay efficiency of NMD-targeted transcripts (dNMD-QTLs). Many such variants are missed in traditional expression quantitative trait locus (eQTL) mapping. NMD-QTLs show strong tissue specificity especially in the brain. They are more likely to overlap with disease single-nucleotide polymorphisms (SNPs). Compared to eQTLs, NMD-QTLs are more likely to be located within gene bodies and exons, especially the penultimate exons from the 3' end. Furthermore, NMD-QTLs are more likely to be found in the binding sites of miRNAs and RNA binding proteins.

CONCLUSIONS

We reveal the genome-wide landscape of genetic variants associated with NMD regulation across human tissues. Our analysis results indicate important roles of NMD in the brain. The preferential genomic positions of NMD-QTLs suggest key attributes for NMD regulation. Furthermore, the overlap with disease-associated SNPs and post-transcriptional regulatory elements implicates regulatory roles of NMD-QTLs in disease manifestation and their interactions with other post-transcriptional regulators.

Effect of OR51E1 Single Nucleotide Polymorphisms on Glioma Susceptibility in the Chinese Han Population

BACKGROUND

Glioma is one of the common primary intracranial tumors, which is heterogeneous among individuals with a low cure rate. Our study aimed to investigate the association between single nucleotide polymorphisms (SNPs) of the OR51E1 gene and glioma susceptibility in the Chinese Han population.

METHODS

A total of six SNPs on OR51E1 in 1,026 subjects (526 cases and 500 controls) were genotyped by MassARRAY iPLEX GOLD assay. The association between these SNPs and glioma susceptibility was analyzed using logistic regression, and odds ratios (ORs) and 95% confidence intervals (CIs) were also calculated. The multifactor dimensionality reduction (MDR) method was applied to detect "SNP-SNP" interactions.

RESULTS

In the overall sample, polymorphisms rs10768148, rs7102992, and rs10500608 were identified to be associated with glioma risk. In the stratified analysis based on gender, only polymorphism rs10768148 was observed to be associated with the risk of glioma. In the age-stratified analysis, rs7102992, rs74052483, and rs10500609 contributed to the risk of glioma in subjects aged > 40 years. And polymorphisms rs10768148 and rs7102992 were associated with the risk of glioma in subjects aged ≤ 40 years and subjects with astrocytoma. In addition, a strong synergistic relationship between rs74052483 and rs10768148, and a strong redundant relationship between rs7102992 and rs10768148 were identified in the study.

CONCLUSIONS

This study demonstrated the association of OR51E1 polymorphisms with glioma susceptibility, providing a basis for assessing glioma risk-associated variants in the Chinese Han population.

GALNT14-mediated O-glycosylation on PHB2 serine-161 enhances cell growth, migration and drug resistance by activating IGF1R cascade in hepatoma cells

The single nucleotide polymorphism (SNP) rs9679162 located on GALNT14 gene predicts therapeutic outcomes in patients with intermediate and advanced hepatocellular carcinoma (HCC), but the molecular mechanism remains unclear. Here, the associations between SNP genotypes, GALNT14 expression, and downstream molecular events were determined. A higher GALNT14 cancerous/noncancerous ratio was associated with the rs9679162-GG genotype, leading to an unfavorable postoperative prognosis. A novel exon-6-skipped GALNT14 mRNA variant was identified in patients carrying the rs9679162-TT genotype, which was associated with lower GALNT14 expression and favorable prognosis. Cell-based experiments showed that elevated levels of GALNT14 promoted HCC growth, migration, and resistance to anticancer drugs. Using a comparative lectin-capture glycoproteomic approach, PHB2 was identified as a substrate for GALNT14-mediated O-glycosylation. Site-directed mutagenesis experiments revealed that serine-161 (Ser161) was the O-glycosylation site. Further analysis showed that O-glycosylation of PHB2-Ser161 was required for the GALNT14-mediated growth-promoting phenotype. O-glycosylation of PHB2 was positively correlated with GALNT14 expression in HCC, resulting in increased interaction between PHB2 and IGFBP6, which in turn led to the activation of IGF1R-mediated signaling. In conclusion, the GALNT14-rs9679162 genotype was associated with differential expression levels of GALNT14 and the generation of a novel exon-6-skipped GALNT14 mRNA variant, which was associated with a favorable prognosis in HCC. The GALNT14/PHB2/IGF1R cascade modulated the growth, migration, and anticancer drug resistance of HCC cells, thereby opening the possibility of identifying new therapeutic targets against HCC.

scAllele: A versatile tool for the detection and analysis of variants in scRNA-seq

Single-cell RNA sequencing (scRNA-seq) data contain rich information at the gene, transcript, and nucleotide levels. Most analyses of scRNA-seq have focused on gene expression profiles, and it remains challenging to extract nucleotide variants and isoform-specific information. Here, we present scAllele, an integrative approach that detects single-nucleotide variants, insertions, deletions, and their allelic linkage with splicing patterns in scRNA-seq. We demonstrate that scAllele achieves better performance in identifying nucleotide variants than other commonly used tools. In addition, the read-specific variant calls by scAllele enables allele-specific splicing analysis, a unique feature not afforded by other methods. Applied to a lung cancer scRNA-seq dataset, scAllele identified variants with strong allelic linkage to alternative splicing, some of which are cancer specific and enriched in cancer-relevant pathways. scAllele represents a versatile tool to uncover multilayer information and previously unidentified biological insights from scRNA-seq data.

Alternatively Spliced Isoforms of the P2X7 Receptor: Structure, Function and Disease Associations

The P2X7 receptor (P2X7R) is an ATP-gated membrane ion channel that is expressed by multiple cell types. Following activation by extracellular ATP, the P2X7R mediates a broad range of cellular responses including cytokine and chemokine release, cell survival and differentiation, the activation of transcription factors, and apoptosis. The P2X7R is made up of three P2X7 subunits that contain specific domains essential for the receptor’s varied functions. Alternative splicing produces P2X7 isoforms that exclude one or more of these domains and assemble in combinations that alter P2X7R function. The modification of the structure and function of the P2X7R may adversely affect cellular responses to carcinogens and pathogens, and alternatively spliced (AS) P2X7 isoforms have been associated with several cancers. This review summarizes recent advances in understanding the structure and function of AS P2X7 isoforms and their associations with cancer and potential role in modulating the inflammatory response.

Comprehensive characterization genetic regulation and chromatin landscape of enhancer-associated long non-coding RNAs and their implication in human cancer

Long non-coding RNAs (lncRNAs) that emanate from enhancer regions (defined as enhancer-associated lncRNAs, or elncRNAs) are emerging as critical regulators in disease progression. However, their biological characteristics and clinical relevance have not been fully portrayed. Here, based on the traditional expression quantitative loci (eQTL) and our optimized residual eQTL method, we comprehensively described the genetic effect on elncRNA expression in more than 300 lymphoblastoid cell lines. Meanwhile, a chromatin atlas of elncRNAs relative to the genetic regulation state was depicted. By applying the maximum likelihood estimate method, we successfully identified causal elncRNAs for protein-coding gene expression reprogramming and showed their associated single nucleotide polymorphisms (SNPs) favor binding of transcription factors. Further epigenome analysis revealed two immune-associated elncRNAs AL662844.4 and LINC01215 possess high levels of H3K27ac and H3K4me1 in human cancer. Besides, pan-cancer analysis of 3D genome, transcriptome, and regulatome data showed they potentially regulate tumor-immune cell interaction through affecting MHC class I genes and CD47, respectively. Moreover, our study showed there exist associations between elncRNA and patient survival. Finally, we made a user-friendly web interface available for exploring the regulatory relationship of SNP-elncRNA-protein-coding gene triplets (http://bio-bigdata.hrbmu.edu.cn/elncVarReg). Our study provides critical mechanistic insights for elncRNA function and illustrates their implications in human cancer.

Integrated Genome and Transcriptome Sequencing to Solve a Neuromuscular Puzzle: Miyoshi Muscular Dystrophy and Early Onset Primary Dystonia in Siblings of the Same Family

BACKGROUND

Neuromuscular disorders (NMD), many of which are hereditary, affect muscular function. Due to advances in high-throughput sequencing technologies, the diagnosis of hereditary NMDs has dramatically improved in recent years.

METHODS AND RESULTS

In this study, we report an family with two siblings exhibiting two different NMD, Miyoshi muscular dystrophy (MMD) and early onset primary dystonia (EOPD). Whole exome sequencing (WES) identified a novel monoallelic frameshift deletion mutation (dysferlin: c.4404delC/p.I1469Sfs∗17) in the Dysferlin gene in the index patient who suffered from MMD. This deletion was inherited from his unaffected father and was carried by his younger sister with EOPD. However, immunostaining staining revealed an absence of dysferlin expression in the proband's muscle tissue and thus suggested the presence of the second underlying mutant allele in dysferlin. Using integrated RNA sequencing (RNA-seq) and whole genome sequencing (WGS) of muscle tissue, a novel deep intronic mutation in dysferlin (dysferlin: c.5341-415A > G) was discovered in the index patient. This mutation caused aberrant mRNA splicing and inclusion of an additional pseudoexon (PE) which we termed PE48.1. This PE was inherited from his unaffected mother. PE48.1 inclusion altered the Dysferlin sequence, causing premature termination of translation.

CONCLUSION

Using integrated genome and transcriptome sequencing, we discovered hereditary MMD and EOPD affecting two siblings of same family. Our results added further weight to the combined use of RNA-seq and WGS as an important method for detection of deep intronic gene mutations, and suggest that integrated sequencing assays are an effective strategy for the diagnosis of hereditary NMDs.

Next Generation Sequencing Based Multiplex Long-Range PCR for Routine Genotyping of Autoinflammatory Disorders

Background

During the last decade, remarkable progress with massive sequencing has been made in the identification of disease-associated genes for AIDs using next-generation sequencing technologies (NGS). An international group of experts described the ideal genetic screening method which should give information about SNVs, InDels, Copy Number Variations (CNVs), GC rich regions. We aimed to develop and validate a molecular diagnostic method in conjunction with the NGS platform as an inexpensive, extended and uniform coverage and fast screening tool which consists of nine genes known to be associated with various AIDs.

Methods

For the validation of basic and expanded panels, long-range multiplex models were setup on healthy samples without any known variations for MEFV, MVK, TNFRSF1A, NLRP3, PSTPIP1, IL1RN, NOD2, NLRP12 and LPIN2 genes. Patients with AIDs who had already known causative variants in these genes were sequenced for analytical validation. As a last step, multiplex models were validated on patients with pre-diagnosis of AIDs. All sequencing steps were performed on the Illumina NGS platform. Validity steps included the selection of related candidate genes, primer design, development of screening methods, validation and verification of the product. The GDPE (Gentera) bioinformatics pipeline was followed.

Results

Although there was no nonsynonymous variation in 21 healthy samples, 107 synonymous variant alleles and some intronic and UTR variants were detected. In 10 patients who underwent analytical validation, besides the 11 known nonsynonymous variant alleles, 11 additional nonsynonymous variant alleles and a total of 81 synonymous variants were found. In the clinical validation phase, 46 patients sequenced with multiplex panels, genetic and clinical findings were combined for diagnosis.

Conclusion

In this study, we describe the development and validation of an NGS-based multiplex array enabling the "long-amplicon" approach for targeted sequencing of nine genes associated with common AIDs. This screening tool is less expensive and more comprehensive compared to other methods and more informative than traditional sequencing. The proposed panel offers advantages to WES or hybridization probe equivalents in terms of CNV analysis, high sensitivity and uniformity, GC-rich region sequencing, InDel detection and intron covering.

Systematic evaluation of the effect of polyadenylation signal variants on the expression of disease-associated genes

Single nucleotide variants (SNVs) within polyadenylation signals (PASs), a specific six-nucleotide sequence required for mRNA maturation, can impair RNA-level gene expression and cause human diseases. However, there is a lack of genome-wide investigation and systematic confirmation tools for identifying PAS variants. Here, we present a computational strategy to integrate the most reliable resources for discovering distinct genomic features of PAS variants and also develop a credible and convenient experimental tool to validate the effect of PAS variants on expression of disease-associated genes. This approach will greatly accelerate the deciphering of PAS variation-related human diseases.

Allele-specific alternative splicing and its functional genetic variants in human tissues

Alternative splicing is an RNA processing mechanism that affects most genes in human, contributing to disease mechanisms and phenotypic diversity. The regulation of splicing involves an intricate network of cis-regulatory elements and trans-acting factors. Due to their high sequence specificity, cis-regulation of splicing can be altered by genetic variants, significantly affecting splicing outcomes. Recently, multiple methods have been applied to understanding the regulatory effects of genetic variants on splicing. However, it is still challenging to go beyond apparent association to pinpoint functional variants. To fill in this gap, we utilized large-scale data sets of the Genotype-Tissue Expression (GTEx) project to study genetically modulated alternative splicing (GMAS) via identification of allele-specific splicing events. We demonstrate that GMAS events are shared across tissues and individuals more often than expected by chance, consistent with their genetically driven nature. Moreover, although the allelic bias of GMAS exons varies across samples, the degree of variation is similar across tissues versus individuals. Thus, genetic background drives the GMAS pattern to a similar degree as tissue-specific splicing mechanisms. Leveraging the genetically driven nature of GMAS, we developed a new method to predict functional splicing-altering variants, built upon a genotype-phenotype concordance model across samples. Complemented by experimental validations, this method predicted >1000 functional variants, many of which may alter RNA-protein interactions. Lastly, 72% of GMAS-associated SNPs were in linkage disequilibrium with GWAS-reported SNPs, and such association was enriched in tissues of relevance for specific traits/diseases. Our study enables a comprehensive view of genetically driven splicing variations in human tissues.

Identification of Anthocyanins-Related Glutathione S-Transferase (GST) Genes in the Genome of Cultivated Strawberry (Fragaria × ananassa)

Anthocyanins are responsible for the red color of strawberry, they are a subclass of flavonoids synthesized in cytosol and transferred to vacuole to form the visible color. Previous studies in model and ornamental plants indicated members of the glutathione S-transferase (GST) gene family were involved in vacuolar accumulation of anthocyanins. In the present study, a total of 130 FaGST genes were identified in the genome of cultivated strawberry (Fragaria × ananassa), which were unevenly distributed across the 28 chromosomes from the four subgenomes. Evolutionary analysis revealed the expansion of FaGST family was under stable selection and mainly drove by WGD/segmental duplication event. Classification and phylogenetic analysis indicated that all the FaGST genes were clarified into seven subclasses, among which FaGST1, FaGST37, and FaGST97 belonging to Phi class were closely related to FvRAP, an anthocyanin-related GST of wildwood strawberry, and this clade was clustered with other known anthocyanin-related GSTs. RNAseq-based expression analysis at different developmental stages of strawberry revealed that the expression of FaGST1, FaGST37, FaGST39, FaGST73, and FaGST97 was gradually increased during the fruit ripening, consistent with the anthocyanins accumulation. These expression patterns of those five FaGST genes were also significantly correlated with those of other anthocyanin biosynthetic genes such as FaCHI, FaCHS, and FaANS, as well as anthocyanin regulatory gene FaMYB10. These results indicated FaGST1, FaGST37, FaGST39, FaGST73, and FaGST97 may function in vacuolar anthocyanin accumulation in cultivated strawberry.

A Comprehensive Analysis of 2013 Dystrophinopathies in China: A Report From National Rare Disease Center

Duchenne muscular dystrophy (DMD) and Becker muscular dystrophy (BMD) are X-linked recessive neuromuscular disorders caused by mutations in DMD. A high-quality database of DMD/BMD is essential not only for clinical practice but also for fundamental research. Here, we aimed to build the largest Chinese national dystrophinopathy database using the National Rare Diseases Registry System of China. Peking Union Medical College Hospital (PUMCH) was the National Rare Diseases Center of China. This research involved 2013 patients with dystrophinopathies, whose diagnoses were confirmed; they were registered and followed up at PUMCH from March 2011 to December 2018. Family history, clinical signs, and treatment data were reported for patients with DMD and BMD at different rates. All six serum biochemical indexes could accurately distinguish between DMD and BMD patients. Copy number variations were the most frequent mutation type (79.2% in DMD and 84.3% in BMD), of which large deletions accounted for 88.4 and 88.6%, large duplications accounted for 11.6 and 11.4% in DMD and BMD, respectively. An exon deletion hotspot, located in exons 45-54, was observed in DMD, and intron 44 was the most frequent deletion starting point (26.5%). Duplication and single nucleotide variations appeared to be uniformly distributed among all exons. Eleven patients were identified to have ultrarare mutation types. Eleven other patients suffered from two separate mutations simultaneously, some of which may have taken place via dependent mechanisms. Thus, we have established the largest hospital-based Chinese dystrophinopathy database via the National Rare Diseases Registry System. This study provides valuable information for further diagnostic and therapeutic studies of dystrophinopathy.

Comprehensive fundamental somatic variant calling and quality management strategies for human cancer genomes

Next-generation sequencing (NGS) technology has revolutionised human cancer research, particularly via detection of genomic variants with its ultra-high-throughput sequencing and increasing affordability. However, the inundation of rich cancer genomics data has resulted in significant challenges in its exploration and translation into biological insights. One of the difficulties in cancer genome sequencing is software selection. Currently, multiple tools are widely used to process NGS data in four stages: raw sequence data pre-processing and quality control (QC), sequence alignment, variant calling and annotation and visualisation. However, the differences between these NGS tools, including their installation, merits, drawbacks and application, have not been fully appreciated. Therefore, a systematic review of the functionality and performance of NGS tools is required to provide cancer researchers with guidance on software and strategy selection. Another challenge is the multidimensional QC of sequencing data because QC can not only report varied sequence data characteristics but also reveal deviations in diverse features and is essential for a meaningful and successful study. However, monitoring of QC metrics in specific steps including alignment and variant calling is neglected in certain pipelines such as the 'Best Practices Workflows' in GATK. In this review, we investigated the most widely used software for the fundamental analysis and QC of cancer genome sequencing data and provided instructions for selecting the most appropriate software and pipelines to ensure precise and efficient conclusions. We further discussed the prospects and new research directions for cancer genomics.

Genetic variability in postoperative nausea and vomiting: A systematic review

BACKGROUND

Postoperative nausea and vomiting (PONV) is the most frequent side effect following anaesthesia. Predisposition to developing PONV is multifactorial with patient risk factors and anaesthetic techniques both being contributory. However, there is also a genetic susceptibility to PONV, and several studies have aimed to identify polymorphisms contributing to a genetic PONV risk.

OBJECTIVE

We summarised previous published studies investigating genetic contribution to PONV risk.

DESIGN

Systematic review without meta-analysis.

DATA SOURCE

We searched MEDLINE until June 2019.

ELIGIBILITY CRITERIA

Articles were chosen for review when PONV and polymorphisms were included. Exclusion criteria were reviews/meta-analysis/comments, articles not in the English language, nonappropriate content (e.g. PONV not as primary aim of the study, study investigated opioid-induced nausea) or if articles were pharmacogenetic studies addressing treatment of PONV.

RESULTS

A total of 59 studies were screened and 14 articles were reviewed including one genome-wide association study (GWAS). Seven studies were performed in East Asians, and seven in Caucasians. Seventeen polymorphisms have been positively associated with PONV in at least one study. Allele frequency of the investigated polymorphisms differs widely between the ethnicities. Furthermore, the anaesthesia regimen and the postoperative time point at which the association with PONV was reported were quite different. Only two polymorphisms, the CHRM3 rs2165870 and the KCNB2 rs349358 (both first associated with PONV in a GWAS), have been significantly associated with PONV incidence in Caucasians in independent studies.

CONCLUSION

There is a genetic susceptibility to the development of PONV. Two single nucleotide polymorphisms (SNPs), the CHRM3 rs2165870 and the KCNB2 rs349358 SNP, seem to have a major influence on PONV incidence, at least in Caucasians. Both SNPs were primarily identified in a GWAS and this association may lead to a better understanding of the disease aetiology. Further high-quality studies are needed to reveal more insights in genetic PONV susceptibility, particularly so in non-Caucasian ethnicities.

Global Analysis of Alternative Splicing Difference in Peripheral Immune Organs between Tongcheng Pigs and Large White Pigs Artificially Infected with PRRSV In Vivo

Alternative splicing (AS) plays a significant role in regulating gene expression at the transcriptional level in eukaryotes. Flexibility and diversity of transcriptome and proteome can be significantly increased through alternative splicing of genes. In the present study, transcriptome data of peripheral immune organs including spleen and inguinal lymph nodes (ILN) were used to identify AS difference between PRRSV-resistant Tongcheng (TC) pigs and PRRSV-susceptible Large White (LW) pigs artificially infected with porcine reproductive and respiratory syndrome virus (PRRSV) in vivo. The results showed that PRRSV infection induced global alternative splicing events (ASEs) with different modes. Among them, 373 genes and 595 genes in the spleen and ILN of TC pigs, while 458 genes and 560 genes in the spleen and ILN of LW pigs had significantly differential ASEs. Alternative splicing was subject to tissue-specific and lineage-specific regulation in response to PRRSV infection. Enriched GO terms and pathways showed that genes with differential ASEs played important roles in transcriptional regulation, immune response, metabolism, and apoptosis. Furthermore, a splicing factor associated with apoptosis, SRSF4, was significantly upregulated in LW pigs. Functional analysis on apoptosis associated genes was validated by RT-PCR and DNA sequencing. These findings revealed different response to PRRSV between PRRSV-resistant TC pigs and PRRSV-susceptible LW pigs at the level of alternative splicing, suggesting the potential relationship between AS and disease resistance to PRRSV.

A practical view of fine-mapping and gene prioritization in the post-genome-wide association era

Over the past 15 years, genome-wide association studies (GWASs) have enabled the systematic identification of genetic loci associated with traits and diseases. However, due to resolution issues and methodological limitations, the true causal variants and genes associated with traits remain difficult to identify. In this post-GWAS era, many biological and computational fine-mapping approaches now aim to solve these issues. Here, we review fine-mapping and gene prioritization approaches that, when combined, will improve the understanding of the underlying mechanisms of complex traits and diseases. Fine-mapping of genetic variants has become increasingly sophisticated: initially, variants were simply overlapped with functional elements, but now the impact of variants on regulatory activity and direct variant-gene 3D interactions can be identified. Moreover, gene manipulation by CRISPR/Cas9, the identification of expression quantitative trait loci and the use of co-expression networks have all increased our understanding of the genes and pathways affected by GWAS loci. However, despite this progress, limitations including the lack of cell-type- and disease-specific data and the ever-increasing complexity of polygenic models of traits pose serious challenges. Indeed, the combination of fine-mapping and gene prioritization by statistical, functional and population-based strategies will be necessary to truly understand how GWAS loci contribute to complex traits and diseases.

Analysis of genetically driven alternative splicing identifies FBXO38 as a novel COPD susceptibility gene

While many disease-associated single nucleotide polymorphisms (SNPs) are associated with gene expression (expression quantitative trait loci, eQTLs), a large proportion of complex disease genome-wide association study (GWAS) variants are of unknown function. Some of these SNPs may contribute to disease by regulating gene splicing. Here, we investigate whether SNPs that are associated with alternative splicing (splice QTL or sQTL) can identify novel functions for existing GWAS variants or suggest new associated variants in chronic obstructive pulmonary disease (COPD). RNA sequencing was performed on whole blood from 376 subjects from the COPDGene Study. Using linear models, we identified 561,060 unique sQTL SNPs associated with 30,333 splice sites corresponding to 6,419 unique genes. Similarly, 708,928 unique eQTL SNPs involving 15,913 genes were detected at 10% FDR. While there is overlap between sQTLs and eQTLs, 55.3% of sQTLs are not eQTLs. Co-localization analysis revealed that 7 out of 21 loci associated with COPD (p<1x10-6) in a published GWAS have at least one shared causal variant between the GWAS and sQTL studies. Among the genes identified to have splice sites associated with top GWAS SNPs was FBXO38, in which a novel exon was discovered to be protective against COPD. Importantly, the sQTL in this locus was validated by qPCR in both blood and lung tissue, demonstrating that splice variants relevant to lung tissue can be identified in blood. Other identified genes included CDK11A and SULT1A2. Overall, these data indicate that analysis of alternative splicing can provide novel insights into disease mechanisms. In particular, we demonstrated that SNPs in a known COPD GWAS locus on chromosome 5q32 influence alternative splicing in the gene FBXO38.

A transcriptome-wide association study of high-grade serous epithelial ovarian cancer identifies new susceptibility genes and splice variants

We sought to identify susceptibility genes for high-grade serous ovarian cancer (HGSOC) by performing a transcriptome-wide association study of gene expression and splice junction usage in HGSOC-relevant tissue types (N = 2,169) and the largest genome-wide association study available for HGSOC (N = 13,037 cases and 40,941 controls). We identified 25 transcriptome-wide association study significant genes, 7 at the junction level only, including LRRC46 at 19q21.32, (P = 1 × 10-9), CHMP4C at 8q21 (P = 2 × 10-11) and a PRC1 junction at 15q26 (P = 7 × 10-9). In vitro assays for CHMP4C showed that the associated variant induces allele-specific exon inclusion (P = 0.0024). Functional screens in HGSOC cell lines found evidence of essentiality for three of the new genes we identified: HAUS6, KANSL1 and PRC1, with the latter comparable to MYC. Our study implicates at least one target gene for 6 out of 13 distinct genome-wide association study regions, identifying 23 new candidate susceptibility genes for HGSOC.

Production of an aberrant splice variant of CCL5 is not caused by genetic mutation in the mammary glands of mastitis‑infected Holstein cows

Genetic mutations, including single nucleotide polymorphisms (SNPs), result in aberrant alternatively splicing of gene and involves in susceptibility of inflammatory diseases, including bovine mastitis. C‑C motif chemokine ligand 5 (CCL5) is an immune‑associated gene, but its alternative splicing (AS) mechanism of gene expression has not yet been understood. The present study identified the splice variant of CCL5 and the compared differential expression of various transcripts between healthy and mastitic mammary gland tissue from cows. A novel transcript lacking exon 2 with a deletion of 112 bp (referred to as CCL5‑AS) was identified in the mammary gland. The expression of CCL5‑AS was lower compared with CCL5‑reference in the healthy and mastitic mammary tissues. A total of two novel SNPs (g.1647 C>T and g.1804 G>A) were identified in exon 2 of CCL5. Using the splicing mini‑gene reporter assay in bovine mammary epithelial cells (MAC‑T) and 293T cells, it was confirmed that the production of CCL5‑AS was not caused by the two SNPs. The present findings suggested that alternative splicing is one of the mechanisms of CCL5 expression regulation and is involved in mastitis infection, but that genetic mutation was not responsible for the generation of the abnormal transcript of CCL5 in cows.

Allele-specific binding of RNA-binding proteins reveals functional genetic variants in the RNA

Allele-specific protein-RNA binding is an essential aspect that may reveal functional genetic variants (GVs) mediating post-transcriptional regulation. Recently, genome-wide detection of in vivo binding of RNA-binding proteins is greatly facilitated by the enhanced crosslinking and immunoprecipitation (eCLIP) method. We developed a new computational approach, called BEAPR, to identify allele-specific binding (ASB) events in eCLIP-Seq data. BEAPR takes into account crosslinking-induced sequence propensity and variations between replicated experiments. Using simulated and actual data, we show that BEAPR largely outperforms often-used count analysis methods. Importantly, BEAPR overcomes the inherent overdispersion problem of these methods. Complemented by experimental validations, we demonstrate that the application of BEAPR to ENCODE eCLIP-Seq data of 154 proteins helps to predict functional GVs that alter splicing or mRNA abundance. Moreover, many GVs with ASB patterns have known disease relevance. Overall, BEAPR is an effective method that helps to address the outstanding challenge of functional interpretation of GVs.

Non-coding and coding genomic variants distinguish prostate cancer, castration-resistant prostate cancer, familial prostate cancer, and metastatic castration-resistant prostate cancer from each other

A considerable number of deposited variants has provided new possibilities for knowledge discovery in different types of prostate cancer. Here, we analyzed variants located on 3'UTR, 5'UTR, CDs, Intergenic, and Intronic regions in castration-resistant prostate cancer (8496 variants), familial prostate cancer (3241 variants), metastatic castration-resistant prostate cancer (3693 variants), and prostate cancer (16599 variants). Chromosome regions 10p15-p14 and 2p13 were highly enriched (P < 0.00001) for variants located in 3'UTR, 5'UTR, CDs, intergenic, and intronic regions in castration-resistant prostate cancer. In contrast, 10p15-p14, 10q23.3, 12q13.11, 13q12.3, 1q25, and 8p22 regions were enriched (P < 0.001) in familial prostate cancer. In metastatic castration-resistant prostate cancer, 10p15-p14, 10q23.3, 11q22-q23, 14q21.1, and 14q32.13 were highly variant regions (P < 0.001). Chromosome 2 and chromosome 1 hosted many enriched variant regions. AKR1C3, BRCA1, BRCA2, CHGA, CYP19A1, HOXB13, KLK3, and PTEN contained the highest number of 3'UTR, 5'UTR, CDs, Intergenic, and Intronic variants. Network analysis showed that these genes are upstream of important functions including prostate gland development, tumor recurrence, prostate cancer-specific survival, tumor progression, cancer mortality, long-term survival, cancer recurrence, angiogenesis, and AR. Interestingly, all of EGFR, JAK2, NR3C1, PDZD2, and SEMA3C genes had single nucleotide polymorphisms (SNP) in castration-resistant prostate cancer, consistent with high selection pressure on these genes during drug treatment and consequent resistance. High occurrence of variants in 3'UTRs suggests the importance of regulatory variants in different types of prostate cancer; an area that has been neglected compared with coding variants. This study provides a comprehensive overview of genomic regions contributing to different types of prostate cancer.

Convergent evolution on the hypoxia-inducible factor (HIF) pathway genes EGLN1 and EPAS1 in high-altitude ducks

During periods of reduced O₂ supply, the most profound changes in gene expression are mediated by hypoxia-inducible factor (HIF) transcription factors that play a key role in cellular responses to low-O₂ tension. Using target-enrichment sequencing, we tested whether variation in 26 genes in the HIF signaling pathway was associated with high altitude and therefore corresponding O₂ availability in three duck species that colonized the Andes from ancestral low-altitude habitats in South America. We found strong support for convergent evolution in the case of two of the three duck species with the same genes (EGLN1, EPAS1), and even the same exons (exon 12, EPAS1), exhibiting extreme outliers with a high probability of directional selection in the high-altitude populations. These results mirror patterns of adaptation seen in human populations, which showed mutations in EPAS1, and transcriptional regulation differences in EGLN1, causing changes in downstream target transactivation, associated with a blunted hypoxic response.

Rescue of recurrent deep intronic mutation underlying cell type-dependent quantitative NEMO deficiency

X-linked dominant incontinentia pigmenti (IP) and X-linked recessive anhidrotic ectodermal dysplasia with immunodeficiency (EDA-ID) are caused by loss-of-function and hypomorphic IKBKG (also known as NEMO) mutations, respectively. We describe a European mother with mild IP and a Japanese mother without IP, whose 3 boys with EDA-ID died from ID. We identify the same private variant in an intron of IKBKG, IVS4+866 C>T, which was inherited from and occurred de novo in the European mother and Japanese mother, respectively. This mutation creates a new splicing donor site, giving rise to a 44-nucleotide pseudoexon (PE) generating a frameshift. Its leakiness accounts for NF-κB activation being impaired but not abolished in the boys' cells. However, aberrant splicing rates differ between cell types, with WT NEMO mRNA and protein levels ranging from barely detectable in leukocytes to residual amounts in induced pluripotent stem cell-derived (iPSC-derived) macrophages, and higher levels in fibroblasts and iPSC-derived neuronal precursor cells. Finally, SRSF6 binds to the PE, facilitating its inclusion. Moreover, SRSF6 knockdown or CLK inhibition restores WT NEMO expression and function in mutant cells. A recurrent deep intronic splicing mutation in IKBKG underlies a purely quantitative NEMO defect in males that is most severe in leukocytes and can be rescued by the inhibition of SRSF6 or CLK.

Architecture of polymorphisms in the human genome reveals functionally important and positively selected variants in immune response and drug transporter genes

Background

Genetic polymorphisms can contribute to phenotypic differences amongst individuals, including disease risk and drug response. Characterization of genetic polymorphisms that modulate gene expression and/or protein function may facilitate the identification of the causal variants. Here, we present the architecture of genetic polymorphisms in the human genome focusing on those predicted to be potentially functional/under natural selection and the pathways that they reside.

Results

In the human genome, polymorphisms that directly affect protein sequences and potentially affect function are the most constrained variants with the lowest single-nucleotide variant (SNV) density, least population differentiation and most significant enrichment of rare alleles. SNVs which potentially alter various regulatory sites, e.g. splicing regulatory elements, are also generally under negative selection.

Interestingly, genes that regulate the expression of transcription/splicing factors and histones are conserved as a higher proportion of these genes is non-polymorphic, contain ultra-conserved elements (UCEs) and/or has no non-synonymous SNVs (nsSNVs)/coding INDELs. On the other hand, major histocompatibility complex (MHC) genes are the most polymorphic with SNVs potentially affecting the binding of transcription/splicing factors and microRNAs (miRNA) exhibiting recent positive selection (RPS). The drug transporter genes carry the most number of potentially deleterious nsSNVs and exhibit signatures of RPS and/or population differentiation. These observations suggest that genes that interact with the environment are highly polymorphic and targeted by RPS.

Conclusions

In conclusion, selective constraints are observed in coding regions, master regulator genes, and potentially functional SNVs. In contrast, genes that modulate response to the environment are highly polymorphic and under positive selection.

Electronic supplementary material

The online version of this article (10.1186/s40246-018-0175-1) contains supplementary material, which is available to authorized users.

Splicing heterogeneity: separating signal from noise

Single-cell analyses have revealed a tremendous variety among cells in the abundance and chemical composition of RNA. Much of this heterogeneity is due to alternative splicing by the spliceosome. Little is known about how many of the resulting isoforms are biologically functional or just provide noise with little to no impact. The dynamic nature of the spliceosome provides numerous opportunities for regulation but is also the source of stochastic fluctuations. We discuss possible origins of splicing stochasticity, the experimental approaches for studying heterogeneity in isoforms, and the potential biological significance of noisy splicing in development and disease.

RNA editing in nascent RNA affects pre-mRNA splicing

In eukaryotes, nascent RNA transcripts undergo an intricate series of RNA processing steps to achieve mRNA maturation. RNA editing and alternative splicing are two major RNA processing steps that can introduce significant modifications to the final gene products. By tackling these processes in isolation, recent studies have enabled substantial progress in understanding their global RNA targets and regulatory pathways. However, the interplay between individual steps of RNA processing, an essential aspect of gene regulation, remains poorly understood. By sequencing the RNA of different subcellular fractions, we examined the timing of adenosine-to-inosine (A-to-I) RNA editing and its impact on alternative splicing. We observed that >95% A-to-I RNA editing events occurred in the chromatin-associated RNA prior to polyadenylation. We report about 500 editing sites in the 3' acceptor sequences that can alter splicing of the associated exons. These exons are highly conserved during evolution and reside in genes with important cellular function. Furthermore, we identified a second class of exons whose splicing is likely modulated by RNA secondary structures that are recognized by the RNA editing machinery. The genome-wide analyses, supported by experimental validations, revealed remarkable interplay between RNA editing and splicing and expanded the repertoire of functional RNA editing sites.

The Expanding Landscape of Alternative Splicing Variation in Human Populations

Alternative splicing is a tightly regulated biological process by which the number of gene products for any given gene can be greatly expanded. Genomic variants in splicing regulatory sequences can disrupt splicing and cause disease. Recent developments in sequencing technologies and computational biology have allowed researchers to investigate alternative splicing at an unprecedented scale and resolution. Population-scale transcriptome studies have revealed many naturally occurring genetic variants that modulate alternative splicing and consequently influence phenotypic variability and disease susceptibility in human populations. Innovations in experimental and computational tools such as massively parallel reporter assays and deep learning have enabled the rapid screening of genomic variants for their causal impacts on splicing. In this review, we describe technological advances that have greatly increased the speed and scale at which discoveries are made about the genetic variation of alternative splicing. We summarize major findings from population transcriptomic studies of alternative splicing and discuss the implications of these findings for human genetics and medicine.

Annotation-free quantification of RNA splicing using LeafCutter

The excision of introns from pre-mRNA is an essential step in mRNA processing. We developed LeafCutter to study sample and population variation in intron splicing. LeafCutter identifies variable splicing events from short-read RNA-seq data and finds events of high complexity. Our approach obviates the need for transcript annotations and circumvents the challenges in estimating relative isoform or exon usage in complex splicing events. LeafCutter can be used both to detect differential splicing between sample groups and to map splicing quantitative trait loci (sQTLs). Compared with contemporary methods, our approach identified 1.4-2.1 times more sQTLs, many of which helped us ascribe molecular effects to disease-associated variants. Transcriptome-wide associations between LeafCutter intron quantifications and 40 complex traits increased the number of associated disease genes at a 5% false discovery rate by an average of 2.1-fold compared with that detected through the use of gene expression levels alone. LeafCutter is fast, scalable, easy to use, and available online.

Whole gene sequencing identifies deep-intronic variants with potential functional impact in patients with hypertrophic cardiomyopathy

BACKGROUND

High throughput sequencing technologies have revolutionized the identification of mutations responsible for genetic diseases such as hypertrophic cardiomyopathy (HCM). However, approximately 50% of individuals with a clinical diagnosis of HCM have no causal mutation identified. This may be due to the presence of pathogenic mutations located deep within the introns, which are not detected by conventional sequencing analysis restricted to exons and exon-intron boundaries.

OBJECTIVE

The aim of this study was to develop a whole-gene sequencing strategy to prioritize deep intronic variants that may play a role in HCM pathogenesis.

METHODS AND RESULTS

The full genomic DNA sequence of 26 genes previously associated with HCM was analysed in 16 unrelated patients. We identified likely pathogenic deep intronic variants in VCL, PRKAG2 and TTN genes. These variants, which are predicted to act through disruption of either splicing or transcription factor binding sites, are 3-fold more frequent in our cohort of probands than in normal European populations. Moreover, we found a patient that is compound heterozygous for a splice site mutation in MYBPC3 and the deep intronic VCL variant. Analysis of family members revealed that carriers of the MYBPC3 mutation alone do not manifest the disease, while family members that are compound heterozygous are clinically affected.

CONCLUSION

This study provides a framework for scrutinizing variation along the complete intronic sequence of HCM-associated genes and prioritizing candidates for mechanistic and functional analysis. Our data suggest that deep intronic variation contributes to HCM phenotype.

The Evolutionary Relationship between Alternative Splicing and Gene Duplication

The protein diversity that exists today has resulted from various evolutionary processes. It is well known that gene duplication (GD) along with the accumulation of mutations are responsible, among other factors, for an increase in the number of different proteins. The gene structure in eukaryotes requires the removal of non-coding sequences, introns, to produce mature mRNAs. This process, known as cis-splicing, referred to here as splicing, is regulated by several factors which can lead to numerous splicing arrangements, commonly designated as alternative splicing (AS). AS, producing several transcripts isoforms form a single gene, also increases the protein diversity. However, the evolution and manner for increasing protein variation differs between AS and GD. An important question is how are patterns of AS affected after a GD event. Here, we review the current knowledge of AS and GD, focusing on their evolutionary relationship. These two processes are now considered the main contributors to the increasing protein diversity and therefore their relationship is a relevant, yet understudied, area of evolutionary study.

The roles of RNA processing in translating genotype to phenotype

A goal of human genetics studies is to determine the mechanisms by which genetic variation produces phenotypic differences that affect human health. Efforts in this respect have previously focused on genetic variants that affect mRNA levels by altering epigenetic and transcriptional regulation. Recent studies show that genetic variants that affect RNA processing are at least equally as common as, and are largely independent from, those variants that affect transcription. We highlight the impact of genetic variation on pre-mRNA splicing and polyadenylation, and on the stability, translation and structure of mRNAs as mechanisms that produce phenotypic traits. These results emphasize the importance of including RNA processing signals in analyses to identify functional variants.

The roles of RNA processing in translating genotype to phenotype

A goal of human genetics studies is to determine the mechanisms by which genetic variation produces phenotypic differences that affect human health. Efforts in this respect have previously focused on genetic variants that affect mRNA levels by altering epigenetic and transcriptional regulation. Recent studies show that genetic variants that affect RNA processing are at least equally as common as, and are largely independent from, those variants that affect transcription. We highlight the impact of genetic variation on pre-mRNA splicing and polyadenylation, and on the stability, translation and structure of mRNAs as mechanisms that produce phenotypic traits. These results emphasize the importance of including RNA processing signals in analyses to identify functional variants.