Identification of common genetic variation that modulates alternative splicing

Towards understanding sex differences in autism spectrum disorders

Autism Spectrum Disorder (ASD) is a neurodevelopmental condition characterized by social deficits, repetitive behaviours and lack of empathy. Its significant genetic heritability and potential comorbidities often lead to diagnostic and therapeutic challenges. This review addresses the biological basis of ASD, focusing on the sex differences in gene expression and hormonal influences. ASD is more commonly diagnosed in males at a ratio of 4:1, indicating a potential oversight in female-specific ASD research and a risk of underdiagnosis in females. We consider how ASD manifests differently across sexes by exploring differential gene expression in female and male brains and consider how variations in steroid hormones influence ASD characteristics. Synaptic function, including excitation/inhibition ratio imbalance, is influenced by gene mutations and this is explored as a key factor in the cognitive and behavioural manifestations of ASD. We also discuss the role of micro RNAs (miRNAs) and highlight a novel mutation in miRNA-873, which affects a suite of key synaptic genes, neurexin, neuroligin, SHANK and post-synaptic density proteins, implicated in the pathology of ASD. Our review suggests that genetic predisposition, sex differences in brain gene expression, and hormonal factors significantly contribute to the presentation, identification and severity of ASD, necessitating sex-specific considerations in diagnosis and treatments. These findings advocate for personalized interventions to improve the outcomes for individuals with ASD.

Circadian clock gene polymorphisms implicated in human pathologies

Circadian rhythms, ~24 h cycles of physiological and behavioral processes, can be synchronized by external signals (e.g., light) and persist even in their absence. Consequently, dysregulation of circadian rhythms adversely affects the well-being of the organism. This timekeeping system is generated and sustained by a genetically encoded endogenous mechanism composed of interlocking transcriptional/translational feedback loops that generate rhythmic expression of core clock genes. Genome-wide association studies (GWAS) and forward genetic studies show that SNPs in clock genes influence gene regulation and correlate with the risk of developing various conditions. We discuss genetic variations in core clock genes that are associated with various phenotypes, their implications for human health, and stress the need for thorough studies in this domain of circadian regulation.

A Novel Missense SNP in the Fatty Acid-Binding Protein 4 (FABP4) Gene is Associated with Growth Traits in Karakul and Awassi Sheep

The fatty acid-binding protein 4 (FABP4) plays a crucial role in the transportation and metabolism of fatty acids. It binds to long-chain fatty acids and facilitates their transport within cells. FABP4 is involved in lipid metabolism, insulin sensitivity, inflammation, and energy homeostasis. This study was conducted to assess the association between the FABP4 gene and growth traits in Karakul and Awassi sheep. A PCR-single strand conformation polymorphism (SSCP) protocol was utilized to assess the polymorphism of FABP4 PCR products with growth traits measured at birth, three, six, nine, and twelve-month intervals. One non-synonymous SNP was identified in the second exon, in which threonine was converted to aspartate in the 61st position in FABP4 (p.61Thr > Asp). This novel SNP showed significant associations with all growth traits measured at all age intervals. The results showed that lambs with the TT genotype exhibited higher growth traits than those with the GT and GG genotypes, respectively. The conducted prediction showed a clearly deleterious effect of p.61Thr > Asp on FABP4 structure, which was accompanied by reduced fatty acid binding efficiency. Owing to the predicted damaging effects caused by p.61Thr > Asp on FABP, lower levels of lipid transport and its consequent increased weight gain and other growth trait indices are expected. Therefore, a putative mechanism through which lambs with these genotypes exhibit higher growth traits is proposed. The FABP4 gene is suggested as a promising marker to improve growth traits in Karakul and Awassi sheep. However, more research is required to validate this mechanism.

Bioinformatics pipeline for the systematic mining genomic and proteomic variation linked to rare diseases: The example of monogenic diabetes

Monogenic diabetes is characterized as a group of diseases caused by rare variants in single genes. Like for other rare diseases, multiple genes have been linked to monogenic diabetes with different measures of pathogenicity, but the information on the genes and variants is not unified among different resources, making it challenging to process them informatically. We have developed an automated pipeline for collecting and harmonizing data on genetic variants linked to monogenic diabetes. Furthermore, we have translated variant genetic sequences into protein sequences accounting for all protein isoforms and their variants. This allows researchers to consolidate information on variant genes and proteins linked to monogenic diabetes and facilitates their study using proteomics or structural biology. Our open and flexible implementation using Jupyter notebooks enables tailoring and modifying the pipeline and its application to other rare diseases.

A missense SNP in the proopiomelanocortin (POMC) gene is associated with growth traits in Awassi and Karakul sheep

This study was conducted to assess the association between proopiomelanocortin (POMC) gene and growth traits in Awassi and Karakul sheep. PCR-single strand conformation polymorphism (SSCP) method was utilized to assess the polymorphism of POMC PCR amplicons with body weight and length, wither and rump height, chest and abdominal circumference measured at birth, 3, 6, 9, and 12 months intervals. Only one missense SNP (rs424417456:C > A) was detected in exon-2, in which glycine was converted to cysteine in the 65th position in POMC (p.65Gly > Cys). rs424417456 SNP showed significant associations with all growth traits in the third, sixth, ninth, and twelfth months. At the age of 3 months onward, lambs with CC genotype showed higher body weight, body length, wither and rump heights, and chest and abdominal circumferences than lambs with CA and AA genotypes, respectively. Prediction analyses indicated a deleterious effect of p.65Gly > Cys on POMC structure, function, and stability. Owing to the strong correlation between rs424417456:CC and better growth characteristics, this genotype is proposed as a promising marker to enhance growth traits in Awassi and Karakul sheep. The predicted damaging effects caused by rs424417456:CA and rs424417456:AA genotypes may entail a putative mechanism through which lambs with these genotypes exhibit lower growth traits.

Polymorphic variations and mRNA expression of the genes encoding interleukins as well as enzymes of oxidative and nitrative stresses as a potential risk of nephrolithiasis development

Urolithiasis is one of the most common urological diseases worldwide with an unclear aetiology. However, a growing body of evidence suggests the potential role of molecular disturbances of the inflammation as well as oxidative and nitrative stresses, in the pathogenesis of urolithiasis. Therefore, we aimed to detect the potential association between six selected single-nucleotide polymorphisms (SNPs) and the development of nephrolithiasis. Moreover, we verified the association of urolithiasis development and mRNA expression of IL-6, IL-8, SOD2, and NOS2 in peripheral blood mononuclear cells (PBMCs). Total genomic DNA and mRNA were isolated from the peripheral blood of 112 patients with urolithiasis and 114 healthy subjects. Using Taq-Man® probes, we genotyped the following SNPs: rs1800797 and rs2069845 in IL-6, rs2227307 in IL-8, rs4880 in SOD2, rs2297518 and rs2779249 in NOS2. In turn, the evaluation of mRNA expression was performed using real-time PCR and 2-ΔCt methods. We found that the C/T genotype of the c.47 T>C-SOD2 SNP increased the frequency of urolithiasis occurrence whereas the T/T homozygote of the same polymorphism decreased the risk of urolithiasis development in the Polish population. Moreover, our study confirmed that patients with urolithiasis were characterised by decreased IL-6, IL-8, and SOD2 mRNA expression levels compared to the controls. In conclusion, our results suggest that polymorphic variants and changes in mRNA expression of IL-6, IL8, SOD2, and NOS2 may be involved in the pathophysiology of urolithiasis.

Polymorphisms in corticotrophin-releasing hormone-proopiomalanocortin (CRH-POMC) system genes: Neuroimmune contributions to psoriasis disease

BACKGROUND

Skin is a target organ and source of the corticotropin-releasing hormone-proopiomelanocortin (CRH-POMC) system, operating as a coordinator and executor of responses to stress. Environmental stress exacerbates and triggers inflammatory skin diseases through modifying the cellular components of the immune system supporting the importance of CRH-POMC system in the pathogenesis of psoriasis. The aim of this study was to analyse the association of CRH-POMC polymorphisms with psoriasis and evaluate transcript expression of lesional psoriatic and normal skin in RNA-seq data.

METHODS

Samples of 104 patients with psoriasis and 174 healthy controls were genotyped for 42 single nucleotide polymorphisms (SNPs) of CRH-POMC using Applied Biosystems SNPlex™ method. The transcript quantification was performed using Salmon software v1.3.0.

RESULTS

This study demonstrated the associations between melanocortin 1 receptor (MC1R) polymorphisms rs2228479, rs3212369, dopachrome tautomerase (DCT) polymorphisms rs7987802, rs2031526, rs9524501 and psoriasis in the Tatar population. Very strong association was evident for the SNP rs7987802 in the DCT gene (pc  = 5.95е-006) in psoriasis patients. Additionally, the haplotype analysis provided AT DCT (rs7992630 and rs7987802) and AGA MC1R (rs3212358, 2228479 and 885479) haplotypes significantly associated (pc  ˂ 0.05) with psoriasis in the Tatar population, supporting the involvement of DCT and MC1R to the psoriasis susceptibility. Moreover, MC1R-203 and DCT-201 expression levels were decreased in psoriasis lesional skin compared with healthy control skin.

CONCLUSIONS

This study is the first to identify genetic variants of the MC1R and DCT genes significantly associated with psoriasis in Tatar population. Our results support potential roles of CRH-POMC system genes and DCT in the pathogenesis of psoriasis.

Identification of disease-related aberrantly spliced transcripts in myeloma and strategies to target these alterations by RNA-based therapeutics

Novel drug discoveries have shifted the treatment paradigms of most hematological malignancies, including multiple myeloma (MM). However, this plasma cell malignancy remains incurable, and novel therapies are therefore urgently needed. Whole-genome transcriptome analyses in a large cohort of MM patients demonstrated that alterations in pre-mRNA splicing (AS) are frequent in MM. This manuscript describes approaches to identify disease-specific alterations in MM and proposes RNA-based therapeutic strategies to eradicate such alterations. As a "proof of concept", we examined the causes of aberrant HMMR (Hyaluronan-mediated motility receptor) splicing in MM. We identified clusters of single nucleotide variations (SNVs) in the HMMR transcript where the altered splicing took place. Using bioinformatics tools, we predicted SNVs and splicing factors that potentially contribute to aberrant HMMR splicing. Based on bioinformatic analyses and validation studies, we provided the rationale for RNA-based therapeutic strategies to selectively inhibit altered HMMR splicing in MM. Since splicing is a hallmark of many cancers, strategies described herein for target identification and the design of RNA-based therapeutics that inhibit gene splicing can be applied not only to other genes in MM but also more broadly to other hematological malignancies and solid tumors as well.

Variation of genes encoding nitric oxide synthases and antioxidant enzymes as potential risks of multiple sclerosis development: a preliminary study

Multiple sclerosis (MS) is a neurodegenerative disease characterized by a variable clinical course and diverse pathophysiology, including nitrative and oxidative stresses as well as inflammation. We aimed to detect the potential association between five selected single-nucleotide polymorphisms (SNPs) in genes encoding nitric oxide synthetases as well as antioxidant enzymes and the development of MS in a Polish population. Genomic DNA was isolated from peripheral blood collected from 142 MS patients and 140 controls. Using Taq-Man® probes, we genotyped the following SNPs: rs1879417 in NOS1, and rs2297518 in NOS2 as well as rs4880 in SOD2, rs7943316 in CAT, rs713041 in GPX4. In the case of rs2297518, the C/C genotype and C allele SNP were associated with an enhanced occurrence of MS, while the C/T, T/T genotypes, and T allele of the same polymorphism reduced this risk. Moreover, the C/C homozygote and C allele of the rs4880 SNP reduced MS risk, while the T allele increased the risk. In addition, the A/T heterozygote of rs7943316 polymorphism was associated with an increased risk of MS occurrence. We also detected that the C/C genotype and C allele of rs713041 decreased the risk of MS, whereas the T/T genotype and T allele increased this risk. In conclusion, the results of our study suggest some links between polymorphic variability in the nitrative/oxidative stress-related genes and the risk of MS development in the Polish population.

Single-Nucleotide Polymorphisms in Oxidative Stress-Related Genes and the Risk of a Stroke in a Polish Population-A Preliminary Study

The present preliminary case-control study was undertaken to detect the potential association of six single nucleotide polymorphisms (SNPs) in oxidative stress-related genes: SOD2 (c.47T > C; rs4880), CAT (c.-89A > T; rs7943316), GPX4 (c.660T > A; rs713041), NOS1 (g.117803515C > T; rs1879417) and NOS2 (c.1823C > T; rs2297518 and c.-227G > C; rs10459953) and the occurrence of a stroke. The SNPs were determined using the TaqMan^® Allelic Discrimination Assay in 107 patients with strokes and 107 age- and sex-matched individuals who had not experienced cerebrovascular accidents. The T alleles of the rs4880 were positively correlated with a stroke (bootstrap OR 1.31; 1.07-1.59 95% CI). In the case of the rs713041, an association with the T allele was found (bootstrap OR 1.36; 1.12-1.67). In addition, the occurrence of a stroke was associated with the presence of the C allele of the rs1879417 (bootstrap OR 1.32; 1.09-1.61). We also found that the C/C genotype and C allele of the rs2297518 increased the risk of a stroke (bootstrap ORs 7.00; 4.34-11.29 and 4.96; 3.88-6.34, respectively). Moreover, the C allele of the rs10459953 was associated with an increased occurrence of this disease (bootstrap OR 1.31; 1.08-1.60). These results indicated that genetics variants in the SOD2, GPX4, NOS1 and NOS2 might be associated with susceptibility to strokes in the Polish population.

Identification and analysis of splicing quantitative trait loci across multiple tissues in the human genome

Alternative splicing (AS) is a fundamental step in eukaryotic mRNA biogenesis. Here, we develop an efficient and reproducible pipeline for the discovery of genetic variants that affect AS (splicing QTLs, sQTLs). We use it to analyze the GTEx dataset, generating a comprehensive catalog of sQTLs in the human genome. Downstream analysis of this catalog provides insight into the mechanisms underlying splicing regulation. We report that a core set of sQTLs is shared across multiple tissues. sQTLs often target the global splicing pattern of genes, rather than individual splicing events. Many also affect the expression of the same or other genes, uncovering regulatory loci that act through different mechanisms. sQTLs tend to be located in post-transcriptionally spliced introns, which would function as hotspots for splicing regulation. While many variants affect splicing patterns by altering the sequence of splice sites, many more modify the binding sites of RNA-binding proteins. Genetic variants affecting splicing can have a stronger phenotypic impact than those affecting gene expression.

Expression Changes Confirm Genomic Variants Predicted to Result in Allele-Specific, Alternative mRNA Splicing

Splice isoform structure and abundance can be affected by either noncoding or masquerading coding variants that alter the structure or abundance of transcripts. When these variants are common in the population, these nonconstitutive transcripts are sufficiently frequent so as to resemble naturally occurring, alternative mRNA splicing. Prediction of the effects of such variants has been shown to be accurate using information theory-based methods. Single nucleotide polymorphisms (SNPs) predicted to significantly alter natural and/or cryptic splice site strength were shown to affect gene expression. Splicing changes for known SNP genotypes were confirmed in HapMap lymphoblastoid cell lines with gene expression microarrays and custom designed q-RT-PCR or TaqMan assays. The majority of these SNPs (15 of 22) as well as an independent set of 24 variants were then subjected to RNAseq analysis using the ValidSpliceMut web beacon (http://validsplicemut.cytognomix.com), which is based on data from the Cancer Genome Atlas and International Cancer Genome Consortium. SNPs from different genes analyzed with gene expression microarray and q-RT-PCR exhibited significant changes in affected splice site use. Thirteen SNPs directly affected exon inclusion and 10 altered cryptic site use. Homozygous SNP genotypes resulting in stronger splice sites exhibited higher levels of processed mRNA than alleles associated with weaker sites. Four SNPs exhibited variable expression among individuals with the same genotypes, masking statistically significant expression differences between alleles. Genome-wide information theory and expression analyses (RNAseq) in tumor exomes and genomes confirmed splicing effects for 7 of the HapMap SNP and 14 SNPs identified from tumor genomes. q-RT-PCR resolved rare splice isoforms with read abundance too low for statistical significance in ValidSpliceMut. Nevertheless, the web-beacon provides evidence of unanticipated splicing outcomes, for example, intron retention due to compromised recognition of constitutive splice sites. Thus, ValidSpliceMut and q-RT-PCR represent complementary resources for identification of allele-specific, alternative splicing.

Sequence analysis and variant identification at the APOC3 gene locus indicates association of rs5218 with BMI in a sample of Kuwaiti's

Background

APOC3 is important in lipid transport and metabolism with limited studies reporting genetic sequence variations in specific ethnic groups. The present study aimed to analyze the full APOC3 sequence among Kuwaiti Arabs and test the association of selected variants with lipid levels and BMI.

Methods

Variants were identified by Sanger sequencing the entire APOC3 gene in 100 Kuwaiti Arabs. Variants and their genotypes were fully characterized and used to construct haplotype blocks. Four variants (rs5128, rs2854117, rs2070668, KUAPOC3N3 g.5196 A > G) were selected for testing association with serum lipid levels and BMI in a cohort (n = 733).

Results

APOC3 sequence (4.3 kb) of a Kuwaiti Arab was deposited in Genbank (accession number KJ437193). Forty-two variants including 3 novels were identified including an “A” insertion at genomic positions 116,700,599–116,700,600 (promoter region) and two substitutions in intron 1 at genomic positions 116,700,819 and 116,701,159. Only three variants, (rs5128, rs2854117, and rs2070668) were analyzed for association of which rs5128 showed a trend for association with increased BMI, TG and VLDL levels that was further investigated using multivariate analysis. A significant association of rs5128 with BMI (p <  0.05) was observed following a dominant genetic model with increased risk by an OR of 4.022 (CI: 1.13–14.30).

Conclusion

The present study is the first to report sequence analysis of APOC3 in an Arab ethnic group. This study supports the inclusion of rs5128 as a marker for assessing genetic risk to dyslipidemia and obesity and the inclusion of the novel variant g.5196 A > G for population stratification of Arabs.

Genetic association of molecular traits: A help to identify causative variants in complex diseases

In the past 15 years, major progresses have been made in the understanding of the genetic basis of regulation of gene expression. These new insights have revolutionized our approach to resolve the genetic variation underlying complex diseases. Gene transcript levels were the first expression phenotypes that were studied. They are heritable and therefore amenable to genome-wide association studies. The genetic variants that modulate them are called expression quantitative trait loci. Their study has been extended to other molecular quantitative trait loci (molQTLs) that regulate gene expression at the various levels, from chromatin state to cellular responses. Altogether, these studies have generated a wealth of basic information on the genome-wide patterns of gene expression and their inter-individual variation. Most importantly, molQTLs have become an invaluable asset in the genetic study of complex diseases. Although the identification of the disease-causing variants on the basis of their overlap with molQTLs requires caution, molQTLs can help to prioritize the relevant candidate gene(s) in the disease-associated regions and bring a functional interpretation of the associated variants, therefore, bridging the gap between genotypes and clinical phenotypes.

Association of BCR/ABL transcript variants with different blood parameters and demographic features in Iraqi chronic myeloid leukemia patients

Background

Chronic myeloid leukemia (CML) is a myeloproliferative neoplasm characterized by the presence of BCR‐ABL fusion gene (GenBank accession NC_000022.11). In the vast majority of CML patients, the typical subtype of BCR‐ABL transcript are b3a2, b2a2 or both. The aim of this study was to determine the different subtypes of BCR‐ABL transcript and their impact on the demographic and hematological parameters in Iraqi patients with CML.

Methods

One hundred patients with chronic phase CML (11 newly diagnosed and 89 imatinib‐resistant) were enrolled in this study. Ribonucleic acid (RNA) was extracted from leukocytes, and complementary DNA was created using reverse transcriptase polymerase chain reaction technique. A multiplex polymerase chain reaction with four specific primers was used to determine the BCR‐ABL fusion subtypes in each patient.

Results

Male to female ratio was 1.38:1. Fifty‐nine patients expressed b3a2 transcript, whereas 39 of the remaining cases were positive for b2a2 variant. One case expressed b2a3 transcript, while the last case coexpressed the two subtypes of mRNA b3a2/b2a2. Male and female were significantly associated with b3a2 and b2a2 subtypes, respectively. The b3a2 subtype showed higher total leukocyte count than b2a2 subgroup, while b2a2 variant demonstrated significantly elevated platelet counts compared to those with b3a2 transcript. A significantly higher plateletcrit percentage (PCT%) was found in patients with b2a2 transcript whereas.

Conclusions

The testified Iraqi group expressed M‐BCR‐ABL type with preponderance of b3a2 over b2a2 subtype. There was a gender‐skewed distribution in BCR‐ABL transcript types with b3a2 transcript more prevalent in males. The type of BCR‐ABL transcript is reflected by different leukocyte and platelet counts at diagnosis, which might represent a distinct phenotype and disease biology.

Current perspectives in assessing humoral immunity after measles vaccination

INTRODUCTION

Repeated measles outbreaks in countries with relatively high vaccine coverage are mainly due to failure to vaccinate and importation; however, cases in immunized individuals exist raising questions about suboptimal measles vaccine-induced humoral immunity and/or waning immunity in a low measles-exposure environment.

AREAS COVERED

The plaque reduction neutralization measurement of functional measles-specific antibodies correlates with protection is the gold standard in measles serology, but it does not assess cellular-immune or other parameters that may be associated with durable and/or protective immunity after vaccination. Additional correlates of protection and long-term immunity and new determinants/signatures of vaccine responsiveness such as specific CD46 and IFI44L genetic variants associated with neutralizing antibody titers after measles vaccination are under investigation. Current and future systems biology studies, coupled with new technology/assays and analytical approaches, will lead to an increasingly sophisticated understanding of measles vaccine-induced humoral immunity and will identify 'signatures' of protective and durable immune responses.

EXPERT OPINION

This will translate into the development of highly predictive assays of measles vaccine efficacy, effectiveness, and durability for prospective identification of potential low/non-responders and susceptible individuals who require additional vaccine doses. Such new advances may drive insights into the development of new/improved vaccine formulations and delivery systems.

Natural variation in genes potentially involved in plant architecture and adaptation in switchgrass (Panicum virgatum L.)

Background

Advances in genomic technologies have expanded our ability to accurately and exhaustively detect natural genomic variants that can be applied in crop improvement and to increase our knowledge of plant evolution and adaptation. Switchgrass (Panicum virgatum L.), an allotetraploid (2n = 4× = 36) perennial C4 grass (Poaceae family) native to North America and a feedstock crop for cellulosic biofuel production, has a large potential for genetic improvement due to its high genotypic and phenotypic variation. In this study, we analyzed single nucleotide polymorphism (SNP) variation in 372 switchgrass genotypes belonging to 36 accessions for 12 genes putatively involved in biomass production to investigate signatures of selection that could have led to ecotype differentiation and to population adaptation to geographic zones.

Results

A total of 11,682 SNPs were mined from ~ 15 Gb of sequence data, out of which 251 SNPs were retained after filtering. Population structure analysis largely grouped upland accessions into one subpopulation and lowland accessions into two additional subpopulations. The most frequent SNPs were in homozygous state within accessions. Sixty percent of the exonic SNPs were non-synonymous and, of these, 45% led to non-conservative amino acid changes. The non-conservative SNPs were largely in linkage disequilibrium with one haplotype being predominantly present in upland accessions while the other haplotype was commonly present in lowland accessions. Tajima’s test of neutrality indicated that PHYB, a gene involved in photoperiod response, was under positive selection in the switchgrass population. PHYB carried a SNP leading to a non-conservative amino acid change in the PAS domain, a region that acts as a sensor for light and oxygen in signal transduction.

Conclusions

Several non-conservative SNPs in genes potentially involved in plant architecture and adaptation have been identified and led to population structure and genetic differentiation of ecotypes in switchgrass. We suggest here that PHYB is a key gene involved in switchgrass natural selection. Further analyses are needed to determine whether any of the non-conservative SNPs identified play a role in the differential adaptation of upland and lowland switchgrass.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1193-2) contains supplementary material, which is available to authorized users.

The development of functional mapping by three sex-related loci on the third whorl of different sex types of Carica papaya L

Carica papaya L. is an important economic crop worldwide and is used as a model plant for sex-determination research. To study the different flower sex types, we screened sex-related genes using alternative splicing sequences (AS-seqs) from a transcriptome database of the three flower sex types, i.e., males, females, and hermaphrodites, established at 28 days before flowering using 15 bacterial artificial chromosomes (BACs) of C. papaya L. After screening, the cDNA regions of the three sex-related loci, including short vegetative phase-like (CpSVPL), the chromatin assembly factor 1 subunit A-like (CpCAF1AL), and the somatic embryogenesis receptor kinase (CpSERK), which contained eight sex-related single-nucleotide polymorphisms (SNPs) from the different sex types of C. papaya L., were genotyped using high-resolution melting (HRM). The three loci were examined regarding the profiles of the third whorl, as described below. CpSVPL, which had one SNP associated with the three sex genotypes, was highly expressed in the male and female sterile flowers (abnormal hermaphrodite flowers) that lacked the fourth whorl structure. CpCAF1AL, which had three SNPs associated with the male genotype, was highly expressed in male and normal hermaphrodite flowers, and had no AS-seqs, whereas it exhibited low expression and an AS-seqs in intron 11 in abnormal hermaphrodite flowers. Conversely, carpellate flowers (abnormal hermaphrodite flowers) showed low expression of CpSVPL and AS-seqs in introns 5, 6, and 7 of CpSERK, which contained four SNPs associated with the female genotype. Specifically, the CpSERK and CpCAF1AL loci exhibited no AS-seq expression in the third whorl of the male and normal hermaphrodite flowers, respectively, and variance in the AS-seq expression of all other types of flowers. Functional mapping of the third whorl of normal hermaphrodites indicated no AS-seq expression in CpSERK, low CpSVPL expression, and, for CpCAF1AL, high expression and no AS-seq expression on XY^h-type chromosomes.

Dominant functional role of the novel phosphorylation site S811 in the human renal NaCl cotransporter

The NaCl cotransporter (NCC) is essential for electrolyte homeostasis and control of blood pressure. The human SLC12A3 gene, which encodes NCC, gives rise to 3 isoforms, of which only the shortest isoform [NaCl cotransporter isoform 3 (NCC₃)] has been studied extensively. All NCC isoforms share key phosphorylation sites at T55 and T60 that are essential mediators of NCC function. Recently, a novel phosphorylation site at S811 was identified in isoforms 1 and 2 [NaCl cotransporter splice variant (NCC_SV)], which are only present in humans and higher primates. The aim of the current study, therefore, is to investigate the role of S811 phosphorylation in the regulation of NCC by a combination of biochemical and fluorescent microscopy analyses. We demonstrate that hypotonic low-chloride buffer increases S811 phosphorylation, whereas phosphorylation-deficient S811A mutant hinders phosphorylation at T55 and T60 in NCC_SV and NCC₃. NCC_SV S811A impairs NCC₃ activity in a dominant-negative fashion, although it does not affect plasma membrane abundance. This effect may be explained by the heterodimerization of NCC_SV with NCC₃. Taken together, our study highlights the dominant-negative effect of NCC_SV on T55 and T60 phosphorylation and NCC activity. Here, we reveal a new function of NCC_SV in humans that broadens the understanding on NCC regulation in blood pressure control.-Tutakhel, O. A. Z., Bianchi, F., Smits, D. A., Bindels, R. J. M., Hoenderop, J. G. J., van der Wijst, J. Dominant functional role of the novel phosphorylation site S811 in the human renal NaCl cotransporter.

Haplotype-specific MAPT exon 3 expression regulated by common intronic polymorphisms associated with Parkinsonian disorders

BACKGROUND

Genome wide association studies have identified microtubule associated protein tau (MAPT) H1 haplotype single nucleotide polymorphisms (SNPs) as leading common risk variants for Parkinson's disease, progressive supranuclear palsy and corticobasal degeneration. The MAPT risk variants fall within a large 1.8 Mb region of high linkage disequilibrium, making it difficult to discern the functionally important risk variants. Here, we leverage the strong haplotype-specific expression of MAPT exon 3 to investigate the functionality of SNPs that fall within this H1 haplotype region of linkage disequilibrium.

METHODS

In this study, we dissect the molecular mechanisms by which haplotype-specific SNPs confer allele-specific effects on the alternative splicing of MAPT exon 3. Firstly, we use haplotype-hybrid whole-locus genomic MAPT vectors studies to identify functional SNPs. Next, we characterise the RNA-protein interactions at two loci by mass spectrometry. Lastly, we knockdown candidate splice factors to determine their effect on MAPT exon 3 using a novel allele-specific qPCR assay.

RESULTS

Using whole-locus genomic DNA expression vectors to express MAPT haplotype variants, we demonstrate that rs17651213 regulates exon 3 inclusion in a haplotype-specific manner. We further investigated the functionality of this region using RNA-electrophoretic mobility shift assays to show differential RNA-protein complex formation at the H1 and H2 sequence variants of SNP rs17651213 and rs1800547 and subsequently identified candidate trans-acting splicing factors interacting with these functional SNPs sequences by RNA-protein pull-down experiment and mass spectrometry. Finally, gene knockdown of candidate splice factors identified by mass spectrometry demonstrate a role for hnRNP F and hnRNP Q in the haplotype-specific regulation of exon 3 inclusion.

CONCLUSIONS

We identified common splice factors hnRNP F and hnRNP Q regulating the haplotype-specific splicing of MAPT exon 3 through intronic variants rs1800547 and rs17651213. This work demonstrates an integrated approach to characterise the functionality of risk variants in large regions of linkage disequilibrium.

Domestication reduces alternative splicing expression variations in sorghum

Domestication is known to strongly reduce genomic diversity through population bottlenecks. The resulting loss of polymorphism has been thoroughly documented in numerous cultivated species. Here we investigate the impact of domestication on the diversity of alternative transcript expressions using RNAseq data obtained on cultivated and wild sorghum accessions (ten accessions for each pool). In that aim, we focus on genes expressing two isoforms in sorghum and estimate the ratio between expression levels of those isoforms in each accession. Noticeably, for a given gene, one isoform can either be overexpressed or underexpressed in some wild accessions, whereas in the cultivated accessions, the balance between the two isoforms of the same gene appears to be much more homogenous. Indeed, we observe in sorghum significantly more variation in isoform expression balance among wild accessions than among domesticated accessions. The possibility exists that the loss of nucleotide diversity due to domestication could affect regulatory elements, controlling transcription or degradation of these isoforms. Impact on the isoform expression balance is discussed. As far as we know, this is the first time that the impact of domestication on transcript isoform balance has been studied at the genomic scale. This could pave the way towards the identification of key domestication genes with finely tuned isoform expressions in domesticated accessions while being highly variable in their wild relatives.

Transcriptional signature of lymphoblastoid cell lines of BRCA1, BRCA2 and non-BRCA1/2 high risk breast cancer families

Approximately 25% of hereditary breast cancer cases are associated with a strong familial history which can be explained by mutations in BRCA1 or BRCA2 and other lower penetrance genes. The remaining high-risk families could be classified as BRCAX (non-BRCA1/2) families.

Gene expression involving alternative splicing represents a well-known mechanism regulating the expression of multiple transcripts, which could be involved in cancer development. Thus using RNA-seq methodology, the analysis of transcriptome was undertaken to potentially reveal transcripts implicated in breast cancer susceptibility and development.

RNA was extracted from immortalized lymphoblastoid cell lines of 117 women (affected and unaffected) coming from BRCA1, BRCA2 and BRCAX families. Anova analysis revealed a total of 95 transcripts corresponding to 85 different genes differentially expressed (Bonferroni corrected p-value <0.01) between those groups. Hierarchical clustering allowed distinctive subgrouping of BRCA1/2 subgroups from BRCAX individuals. We found 67 transcripts, which could discriminate BRCAX from BRCA1/BRCA2 individuals while 28 transcripts discriminate affected from unaffected BRCAX individuals.

To our knowledge, this represents the first study identifying transcripts differentially expressed in lymphoblastoid cell lines from major classes of mutation-related breast cancer subgroups, namely BRCA1, BRCA2 and BRCAX. Moreover, some transcripts could discriminate affected from unaffected BRCAX individuals, which could represent potential therapeutic targets for breast cancer treatment.

Genome-wide associations of CD46 and IFI44L genetic variants with neutralizing antibody response to measles vaccine

Population-based studies have revealed 2-10% measles vaccine failure rate even after two vaccine doses. While the mechanisms behind this remain unknown, we hypothesized that host genetic factors are likely to be involved. We performed a genome-wide association study of measles specific neutralizing antibody and IFNγ ELISPOT response in a combined sample of 2872 subjects. We identified two distinct chromosome 1 regions (previously associated with MMR-related febrile seizures), associated with vaccine-induced measles neutralizing antibody titers. The 1q32 region contained 20 significant SNPs in/around the measles virus receptor-encoding CD46 gene, including the intronic rs2724384 (p value = 2.64 × 10^-09) and rs2724374 (p value = 3.16 × 10^-09) SNPs. The 1q31.1 region contained nine significant SNPs in/around IFI44L, including the intronic rs1333973 (p value = 1.41 × 10^-10) and the missense rs273259 (His73Arg, p value = 2.87 × 10^-10) SNPs. Analysis of differential exon usage with mRNA-Seq data and RT-PCR suggests the involvement of rs2724374 minor G allele in the CD46 STP region exon B skipping, resulting in shorter CD46 isoforms. Our study reveals common CD46 and IFI44L SNPs associated with measles-specific humoral immunity, and highlights the importance of alternative splicing/virus cellular receptor isoform usage as a mechanism explaining inter-individual variation in immune response after live measles vaccine.

Using RNA-Seq to Discover Genetic Polymorphisms That Produce Hidden Splice Variants

RNA-seq is a powerful and popular technology for studying posttranscriptional regulation of gene expression, such as alternative splicing. The first step in analyzing RNA-seq data is to map the sequenced reads back to the genome. However, commonly used RNA-seq aligners use the consensus splice site dinucleotide motifs to map reads across splice junctions. This can be deceiving due to genomic variants that create novel splice site dinucleotides, leaving the personal splice junction reads un-mapped to the reference genome. We developed and evaluated a method called RNA Personal Genome Alignment Analyzer (rPGA) to identify "hidden" splicing variations in personal transcriptomes, by mapping personal RNA-seq data to personal genomes. Our work demonstrates that the personal genome approach to RNA-seq read alignment enables the discovery of a large but previously unknown catalog of splicing variations in human populations.

Sequence variation between 462 human individuals fine-tunes functional sites of RNA processing

Recent advances in the cost-efficiency of sequencing technologies enabled the combined DNA- and RNA-sequencing of human individuals at the population-scale, making genome-wide investigations of the inter-individual genetic impact on gene expression viable. Employing mRNA-sequencing data from the Geuvadis Project and genome sequencing data from the 1000 Genomes Project we show that the computational analysis of DNA sequences around splice sites and poly-A signals is able to explain several observations in the phenotype data. In contrast to widespread assessments of statistically significant associations between DNA polymorphisms and quantitative traits, we developed a computational tool to pinpoint the molecular mechanisms by which genetic markers drive variation in RNA-processing, cataloguing and classifying alleles that change the affinity of core RNA elements to their recognizing factors. The in silico models we employ further suggest RNA editing can moonlight as a splicing-modulator, albeit less frequently than genomic sequence diversity. Beyond existing annotations, we demonstrate that the ultra-high resolution of RNA-Seq combined from 462 individuals also provides evidence for thousands of bona fide novel elements of RNA processing-alternative splice sites, introns, and cleavage sites-which are often rare and lowly expressed but in other characteristics similar to their annotated counterparts.

ENPP1 and ESR1 genotypes influence temporomandibular disorders development and surgical treatment response in dentofacial deformities

Dentofacial deformities are dys-morpho-functional disorders involving the temporomandibular joints (TMJ). Many authors have reported a TMJ improvement in dysfunctional subjects with malocclusion after orthodontic or combined orthodontic and surgical treatment particularly for the relief of pain. In particular, few studies have highlighted the demographic and clinical predictors of response to surgical treatment. To date, no genetic factor has yet been identified as a predictor of response to surgical treatment. The aim of this cohort study is therefore to identify single-nucleotide polymorphisms associated with postoperative temporomandibular disorders (TMD) or with TMJ symptoms after orthognathic surgery. Here, we found the AA genotype of SNP rs1643821 (ESR1 gene) as a risk factor for dysfunctional worsening after orthognathic surgery. In addition, we have identified TT genotype of SNP rs858339 (ENPP1 gene) as a protective factor against TMD in a population of patients with dentofacial deformities. Conversely, the heterozygous genotype AT was identified as a risk factor of TMD with respect to the rest of our population. All these elements are particularly important to bring new screening strategies and tailor future treatment.

PERSPECTIVE

This study allows us to identify sub-populations at high risk of developing postoperative temporomandibular disorders after orthognathic surgery procedures. Many other genes of interest could be potential factors influencing the dysfunctional response to orthognathic surgery, particularly genes of the Opera cohort.

Genetic polymorphisms in circadian negative feedback regulation genes predict overall survival and response to chemotherapy in gastric cancer patients

Circadian negative feedback loop (CNFL) genes play important roles in cancer development and progression. To evaluate the effects of single nucleotide polymorphisms (SNPs) in CNFL genes on the survival of GC patients, 13 functional SNPs from 5 CNFL genes were genotyped in a cohort of 1030 resected GC patients (704 in the training set, 326 in the validation set) to explore the association of SNPs with overall survival (OS). Among the 13 SNPs, three SNPs (rs1056560 in CRY1, rs3027178 in PER1 and rs228729 in PER3) were significantly associated with OS of GC in the training set, and verified in the validation set and pooled analysis. Furthermore, a dose-dependent cumulative effect of these SNPs on GC survival was observed, and survival tree analysis showed higher order interactions between these SNPs. In addition, protective effect conferred by adjuvant chemotherapy (ACT) on GC was observed in patients with variant alleles (TG/GG) of rs1056560, but not in those with homozygous wild (TT) genotype. Functional assay suggested rs1056560 genotypes significantly affect CRY1 expression in cancer cells. Our study presents that SNPs in the CNFL genes may be associated with GC prognosis, and provides the guidance in selecting potential GC patients most likely responsive to ACT.

Strictly co-isogenic C57BL/6J-Prnp-/- mice: A rigorous resource for prion science

Although its involvement in prion replication and neurotoxicity during transmissible spongiform encephalopathies is undisputed, the physiological role of the cellular prion protein (PrP(C)) remains enigmatic. A plethora of functions have been ascribed to PrP(C) based on phenotypes of Prnp(-/-) mice. However, all currently available Prnp(-/-) lines were generated in embryonic stem cells from the 129 strain of the laboratory mouse and mostly crossed to non-129 strains. Therefore, Prnp-linked loci polymorphic between 129 and the backcrossing strain resulted in systematic genetic confounders and led to erroneous conclusions. We used TALEN-mediated genome editing in fertilized mouse oocytes to create the Zurich-3 (ZH3) Prnp-ablated allele on a pure C57BL/6J genetic background. Genomic, transcriptional, and phenotypic characterization of Prnp(ZH3/ZH3) mice failed to identify phenotypes previously described in non-co-isogenic Prnp(-/-) mice. However, aged Prnp(ZH3/ZH3) mice developed a chronic demyelinating peripheral neuropathy, confirming the crucial involvement of PrP(C) in peripheral myelin maintenance. This new line represents a rigorous genetic resource for studying the role of PrP(C) in physiology and disease.

A polymorphism in the base excision repair gene PARP2 is associated with differential prognosis by chemotherapy among postmenopausal breast cancer patients

BACKGROUND

Personalized therapy considering clinical and genetic patient characteristics will further improve breast cancer survival. Two widely used treatments, chemotherapy and radiotherapy, can induce oxidative DNA damage and, if not repaired, cell death. Since base excision repair (BER) activity is specific for oxidative DNA damage, we hypothesized that germline genetic variation in this pathway will affect breast cancer-specific survival depending on treatment.

METHODS

We assessed in 1,408 postmenopausal breast cancer patients from the German MARIE study whether cancer specific survival after adjuvant chemotherapy, anthracycline chemotherapy, and radiotherapy is modulated by 127 Single Nucleotide Polymorphisms (SNPs) in 21 BER genes. For SNPs with interaction terms showing p<0.1 (likelihood ratio test) using multivariable Cox proportional hazard analyses, replication in 6,392 patients from nine studies of the Breast Cancer Association Consortium (BCAC) was performed.

RESULTS

rs878156 in PARP2 showed a differential effect by chemotherapy (p=0.093) and was replicated in BCAC studies (p=0.009; combined analysis p=0.002). Compared to non-carriers, carriers of the variant G allele (minor allele frequency=0.07) showed better survival after chemotherapy (combined allelic hazard ratio (HR)=0.75, 95% 0.53-1.07) and poorer survival when not treated with chemotherapy (HR=1.42, 95% 1.08-1.85). A similar effect modification by rs878156 was observed for anthracycline-based chemotherapy in both MARIE and BCAC, with improved survival in carriers (combined allelic HR=0.73, 95% CI 0.40-1.32). None of the SNPs showed significant differential effects by radiotherapy.

CONCLUSIONS

Our data suggest for the first time that a SNP in PARP2, rs878156, may together with other genetic variants modulate cancer specific survival in breast cancer patients depending on chemotherapy. These germline SNPs could contribute towards the design of predictive tests for breast cancer patients.

Discover hidden splicing variations by mapping personal transcriptomes to personal genomes

RNA-seq has become a popular technology for studying genetic variation of pre-mRNA alternative splicing. Commonly used RNA-seq aligners rely on the consensus splice site dinucleotide motifs to map reads across splice junctions. Consequently, genomic variants that create novel splice site dinucleotides may produce splice junction RNA-seq reads that cannot be mapped to the reference genome. We developed and evaluated an approach to identify ‘hidden’ splicing variations in personal transcriptomes, by mapping personal RNA-seq data to personal genomes. Computational analysis and experimental validation indicate that this approach identifies personal specific splice junctions at a low false positive rate. Applying this approach to an RNA-seq data set of 75 individuals, we identified 506 personal specific splice junctions, among which 437 were novel splice junctions not documented in current human transcript annotations. 94 splice junctions had splice site SNPs associated with GWAS signals of human traits and diseases. These involve genes whose splicing variations have been implicated in diseases (such as OAS1), as well as novel associations between alternative splicing and diseases (such as ICA1). Collectively, our work demonstrates that the personal genome approach to RNA-seq read alignment enables the discovery of a large but previously unknown catalog of splicing variations in human populations.

Interpreting functional effects of coding variants: challenges in proteome-scale prediction, annotation and assessment

Accurate assessment of genetic variation in human DNA sequencing studies remains a nontrivial challenge in clinical genomics and genome informatics. Ascribing functional roles and/or clinical significances to single nucleotide variants identified from a next-generation sequencing study is an important step in genome interpretation. Experimental characterization of all the observed functional variants is yet impractical; thus, the prediction of functional and/or regulatory impacts of the various mutations using in silico approaches is an important step toward the identification of functionally significant or clinically actionable variants. The relationships between genotypes and the expressed phenotypes are multilayered and biologically complex; such relationships present numerous challenges and at the same time offer various opportunities for the design of in silico variant assessment strategies. Over the past decade, many bioinformatics algorithms have been developed to predict functional consequences of single nucleotide variants in the protein coding regions. In this review, we provide an overview of the bioinformatics resources for the prediction, annotation and visualization of coding single nucleotide variants. We discuss the currently available approaches and major challenges from the perspective of protein sequence, structure, function and interactions that require consideration when interpreting the impact of putatively functional variants. We also discuss the relevance of incorporating integrated workflows for predicting the biomedical impact of the functionally important variations encoded in a genome, exome or transcriptome. Finally, we propose a framework to classify variant assessment approaches and strategies for incorporation of variant assessment within electronic health records.

Association between the Advanced Glycosylation End Product-Specific Receptor Gene and Cardiovascular Death in Older Men

Advanced glycosylation end product-specific receptor (AGER) signaling has been implicated in atherosclerosis. The aim of this study was to evaluate whether a common genetic variation in the AGER gene is associated with cardiovascular (CV) death. We included 1304 older men who were genotyped for rs1035798:C>T, which is a single nucleotide polymorphism (SNP) mapped to the third intron of AGER. Cox proportional hazard analysis was used to estimate the association of rs1035798:C>T with CV death. In addition we analyzed total RNA extracted from carotid atherosclerosis biopsies of 18 patients that did or did not have recent symptoms of cerebral embolization by quantitative real-time reverse transcription PCR (qRT-PCR). The minor T-allele of rs1035798:C>T was found to be associated with CV death under dominant (HR = 1.43, 95% CI: 1.01–2.02, P = 0.04) and recessive (HR = 2.05, 95% CI: 1.11–3.81, P = 0.02) models of inheritance even after adjustment for traditional cardiovascular risk factors. No association was found between rs1035798:C>T and non-CV death. qRT-PCR results suggested that median relative expression of AGER isoform 1 and isoform 6 transcripts were approximately 6- (P = 0.01) and 2-fold (P = 0.02) greater, respectively, within carotid biopsies of symptomatic compared to asymptomatic patients. These data suggest that the minor (T) allele of rs1035798:C>T represents an independent susceptibility factor for CV death. The expression of AGER isoforms is different in atheroma from patients with recent symptoms. Further studies are needed to investigate if rs1035798:C>T influences the alternative splicing of AGER.

Characterization of a long-chain fatty acid-CoA ligase 1 gene and association between its SNPs and growth traits in the clam Meretrix meretrix

Long-chain fatty acid-CoA ligases (ACSLs) play crucial roles in fatty acid (FA) metabolism. They convert free long-chain FA into acyl-CoAs, which are key intermediates in both anabolic and catabolic pathways. A long-chain fatty acid-CoA ligase gene was cloned in the clam Meretrix meretrix (MmeACSL1), with a full-length cDNA of 1865 bp encoding 475 amino acids. Its expression was only detected in hepatopancreas by semi-quantitative reverse transcription PCR. Expression level of MmeACSL1 exhibited a significant increase in a starvation experiment (P<0.05). This indicates that MmeACSL1 plays an important role in normal metabolism of M. meretrix and may be involved in energy supply and storage. Two exon SNPs and six intron SNPs were developed in this gene by direct sequencing. A marker-trait association analysis showed that five of these SNPs were significantly associated with growth traits (P<0.05). Also, haplotypes comprised of the five SNPs were revealed to be significantly growth-related (P<0.05), which further corroborates the relationship of this gene to growth traits of M. meretrix.

An interleukin-6 receptor polymorphism is associated with opisthorchiasis-linked cholangiocarcinoma risk in Thailand

The cholangiocarcinoma (CCA) is a relatively rare cancer worldwide but it is highly prevalent in Thailand where the liver fluke, Opisthorchis viverrini is endemic. There are reports that interleukin 6 (IL-6) may play an important role in the pathogenesis of opisthorchiasis associated CCA. Functionally, IL-6 can act on target cells through its receptor, IL-6R, and IL-6R polymorphisms may affect the functional activity of IL-6 leading to susceptibility to cholangiocarcinogenesis. Therefore, we assessed the association of the 48892 A/C (Asp358Ala) polymorphism in exon 9 of the IL-6R gene in 79 CCA cases compared to 80 healthy controls using the PCR- RFLP technique. The results showed significant differences between CCA cases and controls in overall genotype (p=0.001) and allele frequencies (p=0.0002). Chi-square for trend test revealed a significant association between genotype and CCA susceptibility (p=0.0002). The odds ratios (ORs) for genotype were 0.283 (95% CI=0.131-0.605, AC vs. AA; p=0.0003) and 0.206 (95% CI=0.196-1.245, CC vs. AA; p=0.0416), the OR for alleles was 0.347 (95% CI=0.187-0.633, allele C vs. allele A; p=0.0002) and that for the carrier C variant was 0.272 (95% CI=0.130-0.564; p=0.0001). This study demonstrated a close association between an IL-6R polymorphism, specifically higher A allele, and cholangiocarcinoma.

Genetic mutations potentially cause two novel NCF1 splice variants up-regulated in the mammary gland, blood and neutrophil of cows infected by Escherichia coli

Neutrophil cytosolic factor 1 (NCF1) plays a crucial role in host defense against microbial pathogens. In this study, we investigated the potential alternative splicing patterns, expression and splice-relevant single nucleotide polymorphisms (SNPs) of the bovine NCF1 gene to increase insights into its potential role against bovine mastitis caused by Escherichia coli infection. Using RT-PCR and clone sequencing methods, we found two novel splice variants designed as NCF1-TV1 (retained intron 6) and NCF1-TV2 (retained part of intron 8), respectively, encoding two putative truncated proteins (239AA and 283AA). Two splice variants were drastically up-regulated in the mastitis-infected cows' mammary tissues, blood and neutrophils compared with these of healthy cows using real-time RT-PCR. Genomic sequencing analysis identified four novel SNPs g.10112 G>A, g.10766 T>C, SNPs g.12085 G>A and g.12430 T>C at the ends of intron 6 and intron 8 of NCF1. ESE motif predicted that the SNP (g.10766 T>C) might affect the binding with splicing-related factors and subsequently caused the production of aberrant splice variant NCF1-TV1, which is one of the potential reasons that the functional SNP was associated with increased milk somatic cell score in cow. Our results would help in better understanding the NCF1 gene function in the process against pathogen infection, and the effect of splicing-related SNP on the production of aberrant splice variant and careful functional characterization in dairy cattle.

Natural variation of gene models in Drosophila melanogaster

Background

Variation within splicing regulatory sequences often leads to differences in gene models among individuals within a species. Two alleles of the same gene may express transcripts with different exon/intron structures and consequently produce functionally different proteins. Matching genomic and transcriptomic data allows us to identify putative regulatory variants associated with changes in splicing patterns.

Results

Here we analyzed natural variation of splicing patterns in the transcriptomes of 81 natural strains of Drosophila melanogaster with known genotypes. We identified dozens of genotype-specific splicing patterns associated with putative cis-splicing quantitative trait loci (sQTL). The majority of changes can be explained by mutations in splice sites. Allelic-imbalance in splicing patterns confirmed that the majority are regulated mainly by cis-genetic effects. Remarkably, allele-specific splicing changes often lead to qualitative changes in gene models, yielding many isoforms not previously annotated. The observed alterations are typically outside protein-coding regions or affect only very short protein segments.

Conclusions

Overall, the sets of gene models appear to be flexible within D. melanogaster populations. The observed variation in splicing patterns are predicted to have limited effects on the encoded protein sequences. To our knowledge, this is the first sQTL mapping study in Drosophila.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1415-6) contains supplementary material, which is available to authorized users.

Integrative genomics and transcriptomics analysis of human embryonic and induced pluripotent stem cells

Background

Human genomic variations, including single nucleotide polymorphisms (SNPs) and copy number variations (CNVs), are associated with several phenotypic traits varying from mild features to hereditary diseases. Several genome-wide studies have reported genomic variants that correlate with gene expression levels in various tissue and cell types.

Results

We studied human embryonic stem cells (hESCs) and human induced pluripotent stem cells (hiPSCs) measuring the SNPs and CNVs with Affymetrix SNP 6 microarrays and expression values with Affymetrix Exon microarrays. We computed the linear relationships between SNPs and expression levels of exons, transcripts and genes, and the associations between gene CNVs and gene expression levels. Further, for a few of the resulted genes, the expression value was associated with both CNVs and SNPs. Our results revealed altogether 217 genes and 584 SNPs whose genomic alterations affect the transcriptome in the same cells. We analyzed the enriched pathways and gene ontologies within these groups of genes, and found out that the terms related to alternative splicing and development were enriched.

Conclusions

Our results revealed that in the human pluripotent stem cells, the expression values of several genes, transcripts and exons were affected due to the genomic variation.

Electronic supplementary material

The online version of this article (doi:10.1186/s13040-014-0032-2) contains supplementary material, which is available to authorized users.

SNPlice: variants that modulate Intron retention from RNA-sequencing data

RATIONALE

The growing recognition of the importance of splicing, together with rapidly accumulating RNA-sequencing data, demand robust high-throughput approaches, which efficiently analyze experimentally derived whole-transcriptome splice profiles.

RESULTS

We have developed a computational approach, called SNPlice, for identifying cis-acting, splice-modulating variants from RNA-seq datasets. SNPlice mines RNA-seq datasets to find reads that span single-nucleotide variant (SNV) loci and nearby splice junctions, assessing the co-occurrence of variants and molecules that remain unspliced at nearby exon-intron boundaries. Hence, SNPlice highlights variants preferentially occurring on intron-containing molecules, possibly resulting from altered splicing. To illustrate co-occurrence of variant nucleotide and exon-intron boundary, allele-specific sequencing was used. SNPlice results are generally consistent with splice-prediction tools, but also indicate splice-modulating elements missed by other algorithms. SNPlice can be applied to identify variants that correlate with unexpected splicing events, and to measure the splice-modulating potential of canonical splice-site SNVs.

AVAILABILITY AND IMPLEMENTATION

SNPlice is freely available for download from https://code.google.com/p/snplice/ as a self-contained binary package for 64-bit Linux computers and as python source-code.

CONTACT

pmudvari@gwu.edu or horvatha@gwu.edu

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Eight SNVs in NF-κB pathway genes and their different performances between subclinical mastitis and mixed Chinese Holstein cows

The nuclear factor-kappa B (NF-κB) pathway proteins are key players in controlling both innate and adaptive immunity. However, the information on NF-κB pathway genes is very limited in mastitis resistance and milk production of Chinese Holstein cows. In this study, we examine the association of the NF-κB pathway gene variants with milk quality traits and somatic cell score (SCS) in Chinese Holstein cows. Eight single nucleotide variants (SNVs) were identified within the bovine NF-κB pathway genes, using DNA pooled sequencing, PCR-RFLP, and forced PCR-RFLP methods. These SNVs include SNV1: g. 536 C>T (exon 10 of Rel), SNV2: g. 94 G>A (exon 20 of p100), SNV3: g. 43 T>C (intron 6 of p105), SNV4: g. 2397 T>G (intron 9 of p105), SNV5: g. 382 G>C (intron 1 of IκBδ), SNV6: g. 21 C>T (exon 5 of IκBζ), SNV7: g. 272 G>A (intron 6 of IκBζ), and SNV8: g. 18 C>T (intron 10 of IκBζ). The association analysis in mixed Chinese Holstein population showed that SNV1 was significantly or highly significantly associated (P<0.01 and P<0.05) with fat rate, protein rate and SCS. Furthermore, the SNV1-CC (wild genotype) determined serine showed the significantly lower SCS and higher milk production traits compared to TT and TC. SNV2 was significantly associated (P<0.05) with SCS; SNV3 was significantly associated (P<0.05) with fat rate; and SNV4 was significantly associated (P<0.05) with fat rate and SCS. In 199 subclinical mastitis Chinese Holstein cows, the statistical results absolutely differed from the mixed Chinese Holstein individuals. Splice-site prediction by SplicePort showed that single nucleotide difference at eight SNVs results in the acceptor score and donor score changing obviously that may lead to alternative splicing. In brief, SNV1, SNV2, SNV3 and SNV4 could be useful genetic markers for mastitis resistance selection and breeding in Chinese Holstein cows. Furthermore, whether these SNVs lead to alternative splicing need further research.

Functional polymorphisms of circadian negative feedback regulation genes are associated with clinical outcome in hepatocellular carcinoma patients receiving radical resection

Previous studies have demonstrated that circadian negative feedback loop genes play an important role in the development and progression of many cancers. However, the associations between single-nucleotide polymorphisms (SNPs) in these genes and the clinical outcomes of hepatocellular carcinoma (HCC) after surgical resection have not been studied so far. Thirteen functional SNPs in circadian genes were genotyped using the Sequenom iPLEX genotyping system in a cohort of 489 Chinese HCC patients who received radical resection. Multivariate Cox proportional hazards model and Kaplan-Meier curve were used for the prognosis analysis. Cumulative effect analysis and survival tree analysis were used for the multiple SNPs analysis. Four individual SNPs, including rs3027178 in PER1, rs228669 and rs2640908 in PER3 and rs3809236 in CRY1, were significantly associated with overall survival (OS) of HCC patients, and three SNPs, including rs3027178 in PER1, rs228729 in PER3 and rs3809236 in CRY1, were significantly associated with recurrence-free survival (RFS). Moreover, we observed a cumulative effect of significant SNPs on OS and RFS (P for trend < 0.001 for both). Survival tree analysis indicated that wild genotype of rs228729 in PER3 was the primary risk factor contributing to HCC patients' RFS. Our study suggests that the polymorphisms in circadian negative feedback loop genes may serve as independent prognostic biomarkers in predicting clinical outcomes for HCC patients who received radical resection. Further studies with different ethnicities are needed to validate our findings and generalize its clinical utility.

Integrative post-genome-wide association analysis of CDKN2A and TP53 SNPs and risk of esophageal adenocarcinoma

Incidence of esophageal adenocarcinoma (EA) in Western countries has increased markedly in recent decades. Although several risk factors have been identified for EA and its precursor, Barrett's esophagus (BE), including reflux, Caucasian race, male gender, obesity, and smoking, less is known about the role of inherited genetic variation. Frequent somatic mutations in the tumor suppressor genes CDKN2A and TP53 were recently reported in EA tumors, while somatic alterations at 9p (CDKN2A) and 17p (TP53) have been implicated as predictors of progression from BE to EA. Motivated by these findings, we used data from a genome-wide association study of 2515 EA cases and 3207 controls to analyze 37 germline single nucleotide polymorphisms at the CDKN2A and TP53 loci. Three CDKN2A polymorphisms were nominally associated (P < 0.05) with reduced risk of EA: rs2518720 C>T [intronic, odds ratio 0.90, P = 0.0121, q = 0.3059], rs3088440 G>A (3'UTR, odds ratio 0.84, P = 0.0186, q = 0.3059), and rs4074785 C>T (intronic, odds ratio 0.85, P = 0.0248, q = 0.3059). None of the TP53 single nucleotide polymorphisms reached nominal significance. Two of the CDKN2A variants identified were also associated with reduced risk of progression from BE to EA, when assessed in a prospective cohort of 408 BE patients: rs2518720 (hazard ratio 0.57, P = 0.0095, q = 0.0285) and rs3088440 (hazard ratio 0.34, P = 0.0368, q = 0.0552). In vitro functional studies of rs3088440, a single nucleotide polymorphism located in the seed sequence of a predicted miR-663b binding site, suggested a mechanism whereby the G>A substitution may attenuate miR-663b-mediated repression of the CDKN2A transcript. This study provides the first evidence that germline variation at the CDKN2A locus may influence EA susceptibility.

A new anti-lipopolysaccharide factor (ALF) gene with its SNP polymorphisms related to WSSV-resistance of Litopenaeus vannamei

Anti-lipopolysaccharide factors (ALFs) of crustacean play an important role against bacteria or virus infection. In this study, the cDNA sequence and genomic sequence of one new isoform of ALF designated as nLvALF1 were reported. The open reading frame (ORF) of nLvALF1 consisted of 369 bp encoding 123 amino acids and the genomic structure of nLvALF1 comprised four introns and three exons. The predicted pI of the deduced protein was 8.82 and the molecular weight (MW) was 13.72 KDa. The deduced amino acid sequence of nLvALF1 contained a typical functional domain of ALF: LPS-binding domain. Phylogenetic analysis showed that nLvALF1 had the closest relationship with FcALF1 from Fenneropenaeus chinensis. The nLvALF1 was specifically expressed in lymphoid organ (Oka) of shrimp. Its transcriptional level was significantly up-regulated after white spot syndrome virus (WSSV) challenge, suggesting that nLvALF1 might participate in defense against WSSV in Litopenaeus vannamei. In order to search potential genetic markers associated with WSSV-resistance, we scanned the polymorphisms of the genomic fragment with 397 bp where the LPS-binding domain encoding sequence located and 18 SNPs were found. The distribution frequency of these SNPs was analyzed in WSSV susceptible shrimp and resistant shrimp separately. Significant differences existed in allelic frequencies at loci g.1361-T > C, g.1370-T > C, g.1419-T > A between the WSSV-resistant group and the WSSV-susceptible/susceptible group (P < 0.05). The specific haplotype CT consisted of g.1415-C > A and g.1419-T > A was associated with susceptibility to WSSV (P < 0.05). These findings provide theoretical support for selection of WSSV-resistant varieties of L. vannamei.

A mutation in a splicing factor that causes retinitis pigmentosa has a transcriptome-wide effect on mRNA splicing

BACKGROUND

Substantial progress has been made in the identification of sequence elements that control mRNA splicing and the genetic variants in these elements that alter mRNA splicing (referred to as splicing quantitative trait loci - sQTLs). Genetic variants that affect mRNA splicing in trans are harder to identify because their effects can be more subtle and diffuse, and the variants are not co-located with their targets. We carried out a transcriptome-wide analysis of the effects of a mutation in a ubiquitous splicing factor that causes retinitis pigmentosa (RP) on mRNA splicing, using exon microarrays.

RESULTS

Exon microarray data was generated from whole blood samples obtained from four individuals with a mutation in the splicing factor PRPF8 and four sibling controls. Although the mutation has no known phenotype in blood, there was evidence of widespread differences in splicing between cases and controls (affecting approximately 20% of exons). Most probesets with significantly different inclusion (defined as the expression intensity of the exon divided by the expression of the corresponding transcript) between cases and controls had higher inclusion in cases and corresponded to exons that were shorter than average, AT rich, located towards the 5' end of the gene and flanked by long introns. Introns flanking affected probesets were particularly depleted for the shortest category of introns, associated with splicing via intron definition.

CONCLUSIONS

Our results show that a mutation in a splicing factor, with a phenotype that is restricted to retinal tissue, acts as a trans-sQTL cluster in whole blood samples. Characteristics of the affected exons suggest that they are spliced co-transcriptionally and via exon definition. However, due to the small sample size available for this study, further studies are required to confirm the widespread impact of this PRPF8 mutation on mRNA splicing outside the retina.

A candidate gene association study for growth performance in an improved giant freshwater prawn (Macrobrachium rosenbergii ) culture line

A candidate gene approach using type I single nucleotide polymorphism (SNP) markers can provide an effective method for detecting genes and gene regions that underlie phenotypic variation in adaptively significant traits. In the absence of available genomic data resources, transcriptomes were recently generated in Macrobrachium rosenbergii to identify candidate genes and markers potentially associated with growth. The characterisation of 47 candidate loci by ABI re-sequencing of four cultured and eight wild samples revealed 342 putative SNPs. Among these, 28 SNPs were selected in 23 growth-related candidate genes to genotype in 200 animals selected for improved growth performance in an experimental GFP culture line in Vietnam. The associations between SNP markers and individual growth performance were then examined. For additive and dominant effects, a total of three exonic SNPs in glycogen phosphorylase (additive), heat shock protein 90 (additive and dominant) and peroxidasin (additive), and a total of six intronic SNPs in ankyrin repeats-like protein (additive and dominant), rolling pebbles (dominant), transforming growth factor-β induced precursor (dominant), and UTP-glucose-1-phosphate uridylyltransferase 2 (dominant) genes showed significant associations with the estimated breeding values in the experimental animals (P =0.001-0.031). Individually, they explained 2.6-4.8 % of the genetic variance (R²=0.026-0.048). This is the first large set of SNP markers reported for M. rosenbergii and will be useful for confirmation of associations in other samples or culture lines as well as having applications in marker-assisted selection in future breeding programs.

Transcriptome sequencing of Chinese and Caucasian population identifies ethnic-associated differential transcript abundance of heterogeneous nuclear ribonucleoprotein K (hnRNPK)

Gene expression variations (GEV) among different ethnic groups have been a subject matter for extensive study. Relatively less known is the extent of alternative splicing variations (ASV) in the context of ethnicity. We conducted a transcriptome sequencing study of 20 lymphoblastoid cell lines obtained from Caucasian and Han Chinese, and identified known genes that exhibit differential isoform abundance between the two ethnic groups. Among them hnRNPK, a co-factor of p53 (TP53), could be further replicated in a 39-sample cohort with TaqMan assay. Although within-population novel splice variants are common, inter-population novel splice variants are rare. We further analyzed 5.63 billion sequencing reads retrieved from the NCBI Sequence Read Archive and identified potential ethnic-specific transcribed regions.

Key role of CRF in the skin stress response system

The discovery of corticotropin-releasing factor (CRF) or CRH defining the upper regulatory arm of the hypothalamic-pituitary-adrenal (HPA) axis, along with the identification of the corresponding receptors (CRFRs 1 and 2), represents a milestone in our understanding of central mechanisms regulating body and local homeostasis. We focused on the CRF-led signaling systems in the skin and offer a model for regulation of peripheral homeostasis based on the interaction of CRF and the structurally related urocortins with corresponding receptors and the resulting direct or indirect phenotypic effects that include regulation of epidermal barrier function, skin immune, pigmentary, adnexal, and dermal functions necessary to maintain local and systemic homeostasis. The regulatory modes of action include the classical CRF-led cutaneous equivalent of the central HPA axis, the expression and function of CRF and related peptides, and the stimulation of pro-opiomelanocortin peptides or cytokines. The key regulatory role is assigned to the CRFR-1α receptor, with other isoforms having modulatory effects. CRF can be released from sensory nerves and immune cells in response to emotional and environmental stressors. The expression sequence of peptides includes urocortin/CRF→pro-opiomelanocortin→ACTH, MSH, and β-endorphin. Expression of these peptides and of CRFR-1α is environmentally regulated, and their dysfunction can lead to skin and systemic diseases. Environmentally stressed skin can activate both the central and local HPA axis through either sensory nerves or humoral factors to turn on homeostatic responses counteracting cutaneous and systemic environmental damage. CRF and CRFR-1 may constitute novel targets through the use of specific agonists or antagonists, especially for therapy of skin diseases that worsen with stress, such as atopic dermatitis and psoriasis.

Association of osteopontin polymorphisms with nasopharyngeal carcinoma risk

No previous study has reported the association of osteopontin polymorphisms with nasopharyngeal carcinoma (NPC) risk. We aimed to investigate the association in a Chinese population. Four variants of osteopontin, rs11730582, rs1126772, rs9138, and rs4754 polymorphisms, were assessed in a case-control study which consists of 108 NPC patients and 210 health controls, by using polymerase chain reaction--restriction fragment length polymorphism method. Serum osteopontin levels were measured by enzyme-linked immunosorbent assay. The serum osteopontin levels were significantly higher in NPC patients than those in controls (P<0.01). Carriers of CC and CT genotype of rs11730582 presented lower serum osteopontin levels than those of TT genotype carriers (P<0.05). Genotypes CT and CT+CC of rs11730582 were associated with the risk of NPC (CT:OR = 0.57, 95%CI = 0.34-0.94; CC+CT:OR = 0.54, 95%CI = 0.34-0.87). Haplotype analysis revealed that haplotype T-A-A-C of rs11730582, rs1126772, rs9138, and rs4754 was associated with NPC risk (OR = 0.49, 95%CI = 0.27-0.86). Stratification analysis showed that genotypes CT and CT+CC of rs11730582 were associated with tumor stage and lymph node metastasis (P<0.05). No associations were found between rs1126772, rs9138, and rs4754 polymorphisms and NPC risk (P>0.05). The variant rs11730582 of osteopontin is associated with NPC risk. It potentially serves as a genetic marker of NPC predisposition.