Differential allelic expression in the human genome: a robust approach to identify genetic and epigenetic cis-acting mechanisms regulating gene expression

Ornaments for efficient allele-specific expression estimation with bias correction

Allele-specific expression plays a crucial role in unraveling various biological mechanisms, including genomic imprinting and gene expression controlled by cis-regulatory variants. However, existing methods for quantification from RNA-sequencing (RNA-seq) reads do not adequately and efficiently remove various allele-specific read mapping biases, such as reference bias arising from reads containing the alternative allele that do not map to the reference transcriptome or ambiguous mapping bias caused by reads containing the reference allele that map differently from reads containing the alternative allele. We present Ornaments, a computational tool for rapid and accurate estimation of allele-specific transcript expression at unphased heterozygous loci from RNA-seq reads while correcting for allele-specific read mapping biases. Ornaments removes reference bias by mapping reads to a personalized transcriptome and ambiguous mapping bias by probabilistically assigning reads to multiple transcripts and variant loci they map to. Ornaments is a lightweight extension of kallisto, a popular tool for fast RNA-seq quantification, that improves the efficiency and accuracy of WASP, a popular tool for bias correction in allele-specific read mapping. In experiments with simulated and human lymphoblastoid cell-line RNA-seq reads with the genomes of the 1000 Genomes Project, we demonstrate that Ornaments improves the accuracy of WASP and kallisto, is nearly as efficient as kallisto, and is an order of magnitude faster than WASP per sample, with the additional cost of constructing a personalized index for multiple samples. Additionally, we show that Ornaments finds imprinted transcripts with higher sensitivity than WASP, which detects imprinted signals only at gene level.

Contribution of miR-124 rs531564 polymorphism to the occurrence of congenital Zika syndrome

Zika virus (ZIKV) cause Congenital Zika Syndrome (CZS) in individuals exposed during pregnancy. Studies have shown that ZIKV infection positively regulates the miR-124 expression in neural cells, which leads to a decrease of TFRC, a gene targeted of this miRNA. Both miR-124 and TFRC exhibit a pivotal role in nervous system development. Therefore, in this study we aimed to investigate whether genetic variants that affect the expression of these genes could act together with ZIKV to increase the risk of individuals developing CZS. TFRC rs406271 and MIR-124-1 rs531564 polymorphisms were genotyped, using TaqMan® Genotyping Assays, in a sample of children who were exposed to ZIKV during pregnancy, of whom 40 were born with CZS and 48 without congenital anomalies. We identified that individuals with CZS presented a higher frequency of CG genotype of rs531564 polymorphism in MIR-124-1 (p=0.048), which is associated with increased expression of miR-124. Since ZIKV also upregulates the expression of this miRNA, the presence of CG genotype in individuals exposed to the virus could lead to a scenario of overexpression of miR-124 in the brain. Since teratogenesis is a multifactorial event, this genetic finding could partly explain why such individuals are more susceptible to CZS, considering both the downregulation of important neurodevelopment genes, as well as deregulation of the neurogenesis process. Thus, we provide preliminary evidence about a possible genetic risk factor to CZS and highlight the importance of analyzing functional polymorphisms related to epigenetic modulators of neurodevelopment genes in the context of ZIKV teratogenesis.

Machine Learning Application Identifies Germline Markers of Hypertension in Patients With Ovarian Cancer Treated With Carboplatin, Taxane, and Bevacizumab

Pharmacogenomics studies how genes influence a person's response to treatment. When complex phenotypes are influenced by multiple genetic variations with little effect, a single piece of genetic information is often insufficient to explain this variability. The application of machine learning (ML) in pharmacogenomics holds great potential - namely, it can be used to unravel complicated genetic relationships that could explain response to therapy. In this study, ML techniques were used to investigate the relationship between genetic variations affecting more than 60 candidate genes and carboplatin-induced, taxane-induced, and bevacizumab-induced toxicities in 171 patients with ovarian cancer enrolled in the MITO-16A/MaNGO-OV2A trial. Single-nucleotide variation (SNV, formerly SNP) profiles were examined using ML to find and prioritize those associated with drug-induced toxicities, specifically hypertension, hematological toxicity, nonhematological toxicity, and proteinuria. The Boruta algorithm was used in cross-validation to determine the significance of SNVs in predicting toxicities. Important SNVs were then used to train eXtreme gradient boosting models. During cross-validation, the models achieved reliable performance with a Matthews correlation coefficient ranging from 0.375 to 0.410. A total of 43 SNVs critical for predicting toxicity were identified. For each toxicity, key SNVs were used to create a polygenic toxicity risk score that effectively divided individuals into high-risk and low-risk categories. In particular, compared with low-risk individuals, high-risk patients were 28-fold more likely to develop hypertension. The proposed method provided insightful data to improve precision medicine for patients with ovarian cancer, which may be useful for reducing toxicities and improving toxicity management.

Altering gene expression using antisense oligonucleotide therapy for hearing loss

Hearing loss affects more than 430 million people, worldwide, and is the third most common chronic physical condition in the United States and Europe (GBD Hearing Loss Collaborators, 2021; NIOSH, 2021; WHO, 2021). The loss of hearing significantly impacts motor and cognitive development, communication, education, employment, and overall quality of life. The inner ear houses the sensory organs for both hearing and balance and provides an accessible target for therapeutic delivery. Antisense oligonucleotides (ASOs) use various mechanisms to manipulate gene expression and can be tailor-made to treat disorders with defined genetic targets. In this review, we discuss the preclinical advancements within the field of the highly promising ASO-based therapies for hereditary hearing loss disorders. Particular focus is on ASO mechanisms of action, preclinical studies on ASO treatments of hearing loss, timing of therapeutic intervention, and delivery routes to the inner ear.

IRIDA Phenotype in TMPRSS6 Monoallelic-Affected Patients: Toward a Better Understanding of the Pathophysiology

Iron-refractory iron deficiency anemia (IRIDA) is an autosomal recessive inherited form of iron deficiency anemia characterized by discrepantly high hepcidin levels relative to body iron status. However, patients with monoallelic exonic TMPRSS6 variants have also been reported to express the IRIDA phenotype. The pathogenesis of an IRIDA phenotype in these patients is unknown and causes diagnostic uncertainty. Therefore, we retrospectively summarized the data of 16 patients (4 men, 12 women) who expressed the IRIDA phenotype in the presence of only a monoallelic TMPRSS6 variant. Eight unaffected relatives with identical exonic TMPRSS6 variants were used as controls. Haplotype analysis was performed to assess the (intra)genetic differences between patients and relatives. The expression and severity of the IRIDA phenotype were highly variable. Compared with their relatives, patients showed lower Hb, MCV, and TSAT/hepcidin ratios and inherited a different wild-type allele. We conclude that IRIDA in monoallelic TMPRSS6-affected patients is a phenotypically and genotypically heterogeneous disease that is more common in female patients. We hypothesize that allelic imbalance, polygenetic inheritance, or modulating environmental factors and their complex interplay are possible causes. This explorative study is the first step toward improved insights into the pathophysiology and improved diagnostic accuracy for patients presenting with IRIDA and a monoallelic exonic TMPRSS6 variant.

Muscle allele-specific expression QTLs may affect meat quality traits in Bos indicus

Single nucleotide polymorphisms (SNPs) located in transcript sequences showing allele-specific expression (ASE SNPs) were previously identified in the Longissimus thoracis muscle of a Nelore (Bos indicus) population consisting of 190 steers. Given that the allele-specific expression pattern may result from cis-regulatory SNPs, called allele-specific expression quantitative trait loci (aseQTLs), in this study, we searched for aseQTLs in a window of 1 Mb upstream and downstream from each ASE SNP. After this initial analysis, aiming to investigate variants with a potential regulatory role, we further screened our aseQTL data for sequence similarity with transcription factor binding sites and microRNA (miRNA) binding sites. These aseQTLs were overlapped with methylation data from reduced representation bisulfite sequencing (RRBS) obtained from 12 animals of the same population. We identified 1134 aseQTLs associated with 126 different ASE SNPs. For 215 aseQTLs, one allele potentially affected the affinity of a muscle-expressed transcription factor to its binding site. 162 aseQTLs were predicted to affect 149 miRNA binding sites, from which 114 miRNAs were expressed in muscle. Also, 16 aseQTLs were methylated in our population. Integration of aseQTL with GWAS data revealed enrichment for traits such as meat tenderness, ribeye area, and intramuscular fat . To our knowledge, this is the first report of aseQTLs identification in bovine muscle. Our findings indicate that various cis-regulatory and epigenetic mechanisms can affect multiple variants to modulate the allelic expression. Some of the potential regulatory variants described here were associated with the expression pattern of genes related to interesting phenotypes for livestock. Thus, these variants might be useful for the comprehension of the genetic control of these phenotypes.

Signals of selection beyond bottlenecks between exotic populations of the bull-headed dung beetle, Onthophagus taurus

Colonization of new environments can lead to population bottlenecks and rapid phenotypic evolution that could be due to neutral and selective processes. Exotic populations of the bull-headed dung beetle (Onthophagus taurus) have differentiated in opposite directions from native beetles in male horn-to-body size allometry and female fecundity. Here we test for genetic and transcriptional differences among two exotic and one native O. taurus populations after three generations in common garden conditions. We sequenced RNA from 24 individuals for each of the three populations including both sexes, and spanning four developmental stages for the two exotic, differentiated populations. Identifying 270,400 high-quality single nucleotide polymorphisms, we revealed a strong signal of genetic differentiation between the three populations, and evidence of recent bottlenecks within and an excess of outlier loci between exotic populations. Differences in gene expression between populations were greatest in prepupae and early adult life stages, stages during which differences in male horn development and female fecundity manifest. Finally, genes differentially expressed between exotic populations also had greater genetic differentiation and performed functions related to chitin biosynthesis and nutrient sensing, possibly underlying allometry and fecundity trait divergences. Our results suggest that beyond bottlenecks, recent introductions have led to genetic and transcriptional differences in genes correlated with observed phenotypic differences.

Induced pluripotent stem cell reprogramming-associated methylation at the GABRA2 promoter and chr4p12 GABAA subunit gene expression in the context of alcohol use disorder

Twin studies indicate that there is a significant genetic contribution to the risk of developing alcohol use disorder (AUD). With the exception of coding variants in ADH1B and ALDH2, little is known about the molecular effects of AUD-associated loci. We previously reported that the AUD-associated synonymous polymorphism rs279858 within the GABAA α2 receptor subunit gene, GABRA2, was associated with gene expression of the chr4p12 GABAA subunit gene cluster in induced pluripotent stem cell (iPSC)-derived neural cultures. Based on this and other studies that showed changes in GABRA2 DNA methylation associated with schizophrenia and aging, we examined methylation in GABRA2. Specifically, using 69 iPSC lines and neural cultures derived from 47 of them, we examined whether GABRA2 rs279858 genotype predicted methylation levels and whether methylation was related to GABAA receptor subunit gene expression. We found that the GABRA2 CpG island undergoes random stochastic methylation during reprogramming and that methylation is associated with decreased GABRA2 gene expression, an effect that extends to the GABRB1 gene over 600 kb distal to GABRA2. Further, we identified additive effects of GABRA2 CpG methylation and GABRA2 rs279858 genotype on expression of the GABRB1 subunit gene in iPSC-derived neural cultures.

Transcriptional differentiation of UV-B protectant genes in maize landraces spanning an elevational gradient in Chiapas, Mexico

Globally, farmers cultivate and maintain crop landraces (i.e., traditional varieties). Landraces contain unique diversity shaped in part by natural and human-mediated selection and are an indispensable resource for farmers. Since environmental conditions change with elevation, crop landraces grown along elevational gradients have provided ideal locations to explore patterns of local adaptation. To further probe traits underlying this differentiation, transcriptome signatures can help provide a foundation for understanding the ways in which functional genetic diversity may be shaped by environment. In this study, we returned to an elevational gradient in Chiapas, Mexico, to assess transcriptional differentiation of genes underlying UV-B protection in locally adapted maize landraces from multiple elevations. We collected and planted landraces from three elevational zones (lowland, approximately 600 m; midland, approximately 1,550 m; highland approximately 2,100 m) in a common garden at 1,531 m. Using RNA-seq data derived from leaf tissue, we performed differential expression analysis between maize from these distinct elevations. Highland and lowland landraces displayed differential expression in phenylpropanoid and flavonoid biosynthesis genes involved in the production of UV-B protectants and did so at a rate greater than expected based on observed background transcriptional differentiation across the genome. These findings provide evidence for the differentiation of suites of genes involved in complex ecologically relevant pathways. Thus, while neutral evolutionary processes may have played a role in the observed patterns of differentiation, UV-B may have also acted as a selective pressure to differentiate maize landraces in the region. Studies of the distribution of functional crop genetic diversity across variable landscapes can aid us in understanding the response of diversity to abiotic/biotic change and, ultimately, may facilitate its conservation and utilization.

Association of VEGF Gene Family Variants with Central Macular Thickness and Visual Acuity after Aflibercept Short-Term Treatment in Diabetic Patients: A Pilot Study

INTRODUCTION

Diabetic retinopathy (DR) is one of the major vision-threatening causes worldwide. Searching for an individualized therapeutic strategy to prevent its progress is challenging.

OBJECTIVE

This work aimed to investigate the association of angiogenesis-inducer vascular endothelial growth factor (VEGF) gene family and related receptor variants (rs833069, rs12366035, rs7664413, rs7993418, and rs2305948) with susceptibility of DR and the response to 1 dose of aflibercept treatment in type 2 diabetes mellitus (T2DM).

METHODS

Consecutive eligible patients with T2DM (n = 125) and 110 unrelated controls were enrolled in this preliminary prospective case-controlled study. Genotyping was identified using TaqMan real-time PCR. Adjusted odds ratio (OR) with 95% confidence interval (CI) was applied to assess the strength of the association with the clinical/ophthalmological characteristics and early response to intravitreal aflibercept treatment in terms of improved visual acuity (BCVA) and central macular thickness (CMT).

RESULTS

We found that both VEGFB rs12366035 and VEGFC rs7664413 conferred higher risk for DR progression under allelic (OR [95% CI]: 1.71 [1.07-2.74]), homozygote comparison (3.55 [1.32-9.57]), and recessive (3.77 [1.43-9.93]) models for the former and under allelic (2.09 [1.25-3.490, homozygote comparison (2.76 [1.02-7.45]), and recessive (2.62 [0.98-6.98] models for the latter. In contrast, VEGFR1 rs7993418 conferred protection against DR under heterozygote comparison and dominant models. The rs12366035*T/T genotype showed the worst pretreatment BCVA score (0.35 ± 0.24) compared to other corresponding genotypes (0.66 ± 0.26 in C/T and 0.54 ± 0.25 in C/C carriers) (p = 0.008). Meanwhile, patients with rs7993418*G/G of VEGFR1 exhibited a significant reduction in CMT after aflibercept injection (12.26 ± 35.43 µ in G/G vs. 3.57 ± 8.74 µ in A/A) (p = 0.037).

CONCLUSIONS

Polymorphisms of the studied VEGF/receptors could be considered as genetic risk factors of DM/DR development and could play an important role in aflibercept early response for DR patients in the study population.

Correlation Between SNPs at the 3'UTR of the FGF2 Gene and Their Interaction with Environmental Factors in Han Chinese Diabetic Peripheral Neuropathy Patients

FGF2 is a neurotrophic factor that can act as a key regulatory molecule of neuroprotection, neurogenesis, and angiogenesis in various injuries. To explore the genetic background of the FGF2 gene on DPN development, this study analyzed the correlation between SNPs in the 3'UTR of the FGF2 gene and their interaction with environmental factors in DPN patients of Han Chinese nationality. Sanger sequencing was used to analyze the FGF2 genotypes at the rs1048201, rs3804158, rs41348645, rs6854081, rs3747676, rs7683093, rs1476215, and rs1476217 loci in 150 DPN patients, 150 NDPN patients, and 150 healthy control patients. Plasma FGF2 levels were measured in all subjects by using ELISAs. Subjects carrying the T allele at the rs1048201 locus in the FGF2 gene had a significantly lower risk of developing DPN compared with subjects carrying the C allele (OR = 0.43, 95% CI = 0.33-0.56, p < 0.01). Subjects with the G genotype at the rs6854081 locus had an exceptionally higher risk of developing DPN than subjects with the T allele (OR = 1.66, 95% CI = 1.39-1.89, p < 0.01). Individuals harboring the G allele at the rs7683093 locus had a markedly higher risk of DPN than patients with the C allele (OR = 1.63, 95% CI = 1.36-1.87, p < 0.01). Finally, individuals having the A genotype at the rs1476215 locus had a significantly higher risk of DPN than individuals carrying the T allele (OR = 1.82, 95% CI = 1.53-2.02, p < 0.01). There was an interaction between age and alcohol consumption and the SNP rs7683093. SNPs at rs1048201, rs6854081, rs7683093, and rs1476215 in the FGF2 3'UTR were strongly associated with plasma levels of FGF2 (p < 0.05). SNPs at the rs1048201, rs6854081, rs7683093, and rs1476215 loci in the FGF2 gene were significantly associated with the risk of DPN. A possible mechanism is that these SNPs affect the expression level of FGF2 by interrupting the binding of microRNAs to target sites in the 3'UTR.

Synonymous GATA2 mutations result in selective loss of mutated RNA and are common in patients with GATA2 deficiency

Deficiency of the transcription factor GATA2 is a highly penetrant genetic disorder predisposing to myelodysplastic syndromes (MDS) and immunodeficiency. It has been recognized as the most common cause underlying primary MDS in children. Triggered by the discovery of a recurrent synonymous GATA2 variant, we systematically investigated 911 patients with phenotype of pediatric MDS or cellular deficiencies for the presence of synonymous alterations in GATA2. In total, we identified nine individuals with five heterozygous synonymous mutations: c.351C>G, p.T117T (N = 4); c.649C>T, p.L217L; c.981G>A, p.G327G; c.1023C>T, p.A341A; and c.1416G>A, p.P472P (N = 2). They accounted for 8.2% (9/110) of cases with GATA2 deficiency in our cohort and resulted in selective loss of mutant RNA. While for the hotspot mutation (c.351C>G) a splicing error leading to RNA and protein reduction was identified, severe, likely late stage RNA loss without splicing disruption was found for other mutations. Finally, the synonymous mutations did not alter protein function or stability. In summary, synonymous GATA2 substitutions are a new common cause of GATA2 deficiency. These findings have broad implications for genetic counseling and pathogenic variant discovery in Mendelian disorders.

Statin Use Associates With Risk of Type 2 Diabetes via Epigenetic Patterns at ABCG1

Statin is the medication most widely prescribed to reduce plasma cholesterol levels. Yet, how the medication contributes to diabetes risk and impaired glucose metabolism is not clear. This study aims to examine the epigenetic mechanisms of ABCG1 through which statin use associates with risk of type 2 diabetes. We determined the association between the statin use, DNA methylation at ABCG1 and type 2 diabetes/glycemic traits in the Framingham Heart Study Offspring (FHS, n = 2741), with validation in the Women’s Health Initiative Study (WHI, n = 2020). The causal effect of statin use on the risk of type 2 diabetes was examined using a two-step Mendelian randomization approach. Next, based on transcriptome analysis, we determined the links between the medication-associated epigenetic status of ABCG1 and biological pathways on the pathogenesis of type 2 diabetes. Our results showed that DNA methylation levels at cg06500161 of ABCG1 were positively associated with the use of statin, type 2 diabetes and related traits (fasting glucose and insulin) in FHS and WHI. Two-step Mendelian randomization suggested a causal effect of statin use on type 2 diabetes and related traits through epigenetic mechanisms, specifically, DNA methylation at cg06500161. Our results highlighted that gene expression of ABCG1, ABCA1 and ACSL3, involved in both cholesterol metabolism and glycemic pathways, was inversely associated with statin use, CpG methylation, and diabetic signatures. We concluded that DNA methylation site cg06500161 at ABCG1 is a mediator of the association between statins and risk of type 2 diabetes.

Allele-Specific Expression of CD4+ T Cells in Response to Marek's Disease Virus Infection

Marek's disease (MD) is a T cell lymphoma disease induced by Marek's disease virus (MDV), a highly oncogenic α herpesvirus primarily affecting chickens. MD is a chronic infectious disease that threatens the poultry industry. However, the mechanisms of genetic resistance for MD are complex and not completely understood. In this study, to identify high-confidence candidate genes of MD genetic resistance, high throughput sequencing (RNA-seq) was used to obtain transcriptomic data of CD4+ T cells isolated from MDV-infected and non-infected groups of two reciprocal crosses of individuals mating by two highly inbred chicken lines (63 MD-resistant and 72 MD-susceptible). After RNA-seq analysis with two biological replicates in each group, we identified 61 and 123 single nucleotide polymorphisms (SNPs) (false discovery rate (FDR) < 0.05) annotated in 39 and 132 genes in intercrosses 63 × 72 and 72 × 63, respectively, which exhibited allele-specific expression (ASE) in response to MDV infection. Similarly, we identified 62 and 79 SNPs annotated in 66 and 96 genes in infected and non-infected groups, respectively. We identified 534 and 1543 differentially expressed genes (DEGs) (FDR < 0.05) related to MDV infection in intercrosses 63 × 72 and 72 × 63, respectively. We also identified 328 and 20 DEGs in infected and non-infected groups, respectively. The qRT-PCR using seven DEGs further verified our results of RNA-seq analysis. The qRT-PCR of 11 important ASE genes was performed for gene functional validation in CD4+ T cells and tumors. Combining the analyses, six genes (MCL1, SLC43A2, PDE3B, ADAM33, BLB1, and DMB2), especially MCL1, were highlighted as the candidate genes with the potential to be involved in MDV infection. Gene-set enrichment analysis revealed that many ASE genes are linked to T cell activation, T cell receptor (TCR), B cell receptor (BCR), ERK/MAPK, and PI3K/AKT-mTOR signaling pathways, which play potentially important roles in MDV infection. Our approach underlines the importance of comprehensive functional studies for gaining valuable biological insight into the genetic factors behind MD and other complex traits, and our findings provide additional insights into the mechanisms of MD and disease resistance breeding in poultry.

Sex-specific effects of a microsatellite polymorphism on human growth hormone receptor gene expression

Growth hormone (GH) binds to its specific receptor (GHR) at the surface of target cells activating multiple signaling pathways implicated in growth and metabolism. Dysregulation of GHRs leads to pathophysiological states that most commonly affect stature. We previously showed the association of a polymorphic (n = 15-37) GT microsatellite in the human GHR gene promoter with short stature in a sex-specific manner. In the present study we evaluated the functional relevance of this polymorphism in regulating GHR expression. Using luciferase reporter assays, we found that the GT repeat had a significant cis regulatory effect in response to HIF1α and a potential repressor role following C/EBPβ stimulation. Using a digital PCR application to measure allelic imbalance (AI), we showed a high prevalence of AI (∼76%) at the GHR locus in lymphoblastoid cell lines (LCLs), with a significantly higher degree of imbalance in LCLs derived from males. Examination of expression of GHR as well as other members of the GH-IGF1 axis in the LCLs revealed significant associations of GHR, IGF1 and BCL2 expression with GT genotype in a sex-specific manner. Our results suggest that this GT microsatellite exerts both cis and trans effects in a sex-specific context, revealing a new mechanism by which GHR gene expression is regulated.

A SNP upstream of the cyclic GMP-AMP synthase (cGAS) gene protects from relapse and extra-pulmonary TB and relates to BCG vaccination status in an Indian cohort

Tuberculosis (TB) caused by Mycobacterium tuberculosis (M.tb) is a major health care threat worldwide causing over a million deaths annually. Host-pathogen interaction is complex, and a strong genetic contribution to disease susceptibility has been proposed. We have investigated single-nucleotide polymorphisms (SNPs) within cGAS/STING in Indian TB patients and healthy cohorts from India and Germany by Lightcycler®480 genotyping technique. The cGAS/STING pathway is an essential defense pathway within the cytosol after M.tb is internalized and mycobacterial DNA is released inducing the production of type I IFNs. We found that the rs311686 SNP upstream of cGAS provides protection from getting TB overall and is differently distributed in pulmonary TB patients compared with extra-pulmonary and particularly relapse cases. This SNP furthermore differs in distribution when comparing individuals with respect to BCG vaccination status. Taken together, our results show that the presence of the rs311686 SNP influences the course of TB significantly. However, structural conformation changes were found only for the cGAS rs610913 SNP. These findings underscore the importance of M.tb DNA recognition for TB pathogenesis and may eventually help in risk stratification of individuals. This may ultimately help in prevention of disease and aid in developing new vaccination and treatment strategies.

Quantification of allelic differential expression using a simple Fluorescence primer PCR-RFLP-based method

Allelic differential expression (ADE) is common in diploid organisms, and is often the key reason for specific phenotype variations. Thus, ADE detection is important for identification of major genes and causal mutations. To date, sensitive and simple methods to detect ADE are still lacking. In this study, we have developed an accurate, simple, and sensitive method, named fluorescence primer PCR-RFLP quantitative method (fPCR-RFLP), for ADE analysis. This method involves two rounds of PCR amplification using a pair of primers, one of which is double-labeled with an overhang 6-FAM. The two alleles are then separated by RFLP and quantified by fluorescence density. fPCR-RFLP could precisely distinguish ADE cross a range of 1- to 32-fold differences. Using this method, we verified PLAG1 and KIT, two candidate genes related to growth rate and immune response traits of pigs, to be ADE both at different developmental stages and in different tissues. Our data demonstrates that fPCR-RFLP is an accurate and sensitive method for detecting ADE on both DNA and RNA level. Therefore, this powerful tool provides a way to analyze mutations that cause ADE.

Investigation of allele-specific expression of genes involved in adipogenesis and lipid metabolism suggests complex regulatory mechanisms of PPARGC1A expression in porcine fat tissues

Background

The expression of genes involved in regulating adipogenesis and lipid metabolism may affect economically important fatness traits in pigs. Allele-specific expression (ASE) reflects imbalance between allelic transcript levels and can be used to identify underlying cis-regulatory elements. ASE has not yet been intensively studied in pigs. The aim of this investigation was to analyze the differential allelic expression of four genes, PPARA, PPARG, SREBF1, and PPARGC1A, which are involved in the regulation of fat deposition in porcine subcutaneous and visceral fat and longissimus dorsi muscle.

Results

Quantification of allelic proportions by pyrosequencing revealed that both alleles of PPARG and SREBF1 are expressed at similar levels. PPARGC1A showed the greatest ASE imbalance in fat deposits in Polish Large White (PLW), Polish Landrace and Pietrain pigs; and PPARA in PLW pigs. Significant deviations of mean PPARGC1A allelic transcript ratio between cDNA and genomic DNA were detected in all tissues, with the most pronounced difference (p < 0.001) in visceral fat of PLW pigs. To search for potential cis-regulatory elements affecting ASE in the PPARGC1A gene we analyzed the effects of four SNPs (rs337351686, rs340650517, rs336405906 and rs345224049) in the promoter region, but none were associated with ASE in the breeds studied. DNA methylation analysis revealed significant CpG methylation differences between samples showing balanced (allelic transcript ratio ≈1) and imbalanced allelic expression for CpG site at the genomic position in chromosome 8 (SSC8): 18527678 in visceral fat (p = 0.017) and two CpG sites (SSC8:18525215, p = 0.030; SSC8:18525237, p = 0.031) in subcutaneous fat.

Conclusions

Our analysis of differential allelic expression suggests that PPARGC1A is subjected to cis-regulation in porcine fat tissues. Further studies are necessary to identify other regulatory elements localized outside the PPARGC1A proximal promoter region.

Electronic supplementary material

The online version of this article (10.1186/s12863-018-0696-6) contains supplementary material, which is available to authorized users.

Genetic Variation in Genes Underlying Diverse Dementias May Explain a Small Proportion of Cases in the Alzheimer's Disease Sequencing Project

BACKGROUND/AIMS

The Alzheimer's Disease Sequencing Project (ADSP) aims to identify novel genes influencing Alzheimer's disease (AD). Variants within genes known to cause dementias other than AD have previously been associated with AD risk. We describe evidence of co-segregation and associations between variants in dementia genes and clinically diagnosed AD within the ADSP.

METHODS

We summarize the properties of known pathogenic variants within dementia genes, describe the co-segregation of variants annotated as "pathogenic" in ClinVar and new candidates observed in ADSP families, and test for associations between rare variants in dementia genes in the ADSP case-control study. The participants were clinically evaluated for AD, and they represent European, Caribbean Hispanic, and isolate Dutch populations.

RESULTS/CONCLUSIONS

Pathogenic variants in dementia genes were predominantly rare and conserved coding changes. Pathogenic variants within ARSA, CSF1R, and GRN were observed, and candidate variants in GRN and CHMP2B were nominated in ADSP families. An independent case-control study provided evidence of an association between variants in TREM2, APOE, ARSA, CSF1R, PSEN1, and MAPT and risk of AD. Variants in genes which cause dementing disorders may influence the clinical diagnosis of AD in a small proportion of cases within the ADSP.

Dyslexia risk variant rs600753 is linked with dyslexia-specific differential allelic expression of DYX1C1

An increasing number of genetic variants involved in dyslexia development were discovered during the last years, yet little is known about the molecular functional mechanisms of these SNPs. In this study we investigated whether dyslexia candidate SNPs have a direct, disease-specific effect on local expression levels of the assumed target gene by using a differential allelic expression assay. In total, 12 SNPs previously associated with dyslexia and related phenotypes were suitable for analysis. Transcripts corresponding to four SNPs were sufficiently expressed in 28 cell lines originating from controls and a family affected by dyslexia. We observed a significant effect of rs600753 on expression levels of DYX1C1 in forward and reverse sequencing approaches. The expression level of the rs600753 risk allele was increased in the respective seven cell lines from members of the dyslexia family which might be due to a disturbed transcription factor binding sites. When considering our results in the context of neuroanatomical dyslexia-specific findings, we speculate that this mechanism may be part of the pathomechanisms underlying the dyslexia-specific brain phenotype. Our results suggest that allele-specific DYX1C1 expression levels depend on genetic variants of rs600753 and contribute to dyslexia. However, these results are preliminary and need replication.

Intrinsic MYH7 expression regulation contributes to tissue level allelic imbalance in hypertrophic cardiomyopathy

HCM, the most common inherited cardiac disease, is mainly caused by mutations in sarcomeric genes. More than a third of the patients are heterozygous for mutations in the MYH7 gene encoding for the β-myosin heavy chain. In HCM-patients, expression of the mutant and the wildtype allele can be unequal, thus leading to fractions of mutant and wildtype mRNA and protein which deviate from 1:1. This so-called allelic imbalance was detected in whole tissue samples but also in individual cells. There is evidence that the severity of HCM not only depends on the functional effect of the mutation itself, but also on the fraction of mutant protein in the myocardial tissue. Allelic imbalance has been shown to occur in a broad range of genes. Therefore, we aimed to examine whether the MYH7-alleles are intrinsically expressed imbalanced or whether the allelic imbalance is solely associated with the disease. We compared the expression of MYH7-alleles in non-HCM donors and in HCM-patients with different MYH7-missense mutations. In the HCM-patients, we identified imbalanced as well as equal expression of both alleles. Also at the protein level, allelic imbalance was determined. Most interestingly, we also discovered allelic imbalance and balance in non-HCM donors. Our findings therefore strongly indicate that apart from mutation-specific mechanisms, also non-HCM associated allelic-mRNA expression regulation may account for the allelic imbalance of the MYH7 gene in HCM-patients. Since the relative amount of mutant mRNA and protein or the extent of allelic imbalance has been associated with the severity of HCM, individual analysis of the MYH7-allelic expression may provide valuable information for the prognosis of each patient.

Bayesian Inference of Allele-Specific Gene Expression Indicates Abundant Cis-Regulatory Variation in Natural Flycatcher Populations

Polymorphism in cis-regulatory sequences can lead to different levels of expression for the two alleles of a gene, providing a starting point for the evolution of gene expression. Little is known about the genome-wide abundance of genetic variation in gene regulation in natural populations but analysis of allele-specific expression (ASE) provides a means for investigating such variation. We performed RNA-seq of multiple tissues from population samples of two closely related flycatcher species and developed a Bayesian algorithm that maximizes data usage by borrowing information from the whole data set and combines several SNPs per transcript to detect ASE. Of 2,576 transcripts analyzed in collared flycatcher, ASE was detected in 185 (7.2%) and a similar frequency was seen in the pied flycatcher. Transcripts with statistically significant ASE commonly showed the major allele in >90% of the reads, reflecting that power was highest when expression was heavily biased toward one of the alleles. This would suggest that the observed frequencies of ASE likely are underestimates. The proportion of ASE transcripts varied among tissues, being lowest in testis and highest in muscle. Individuals often showed ASE of particular transcripts in more than one tissue (73.4%), consistent with a genetic basis for regulation of gene expression. The results suggest that genetic variation in regulatory sequences commonly affects gene expression in natural populations and that it provides a seedbed for phenotypic evolution via divergence in gene expression.

Allele-specific expression at the androgen receptor alpha gene in a hybrid unisexual fish, the Amazon molly (Poecilia formosa)

The all-female Amazon molly (Poecilia formosa) is the result of a hybridization of the Atlantic molly (P. mexicana) and the sailfin molly (P. latipinna) approximately 120,000 years ago. As a gynogenetic species, P. formosa needs to copulate with heterospecific males including males from one of its bisexual ancestral species. However, the sperm only triggers embryogenesis of the diploid eggs. The genetic information of the sperm donor typically will not contribute to the next generation of P. formosa. Hence, P. formosa possesses generally one allele from each of its ancestral species at any genetic locus. This raises the question whether both ancestral alleles are equally expressed in P. formosa. Allele-specific expression (ASE) has been previously assessed in various organisms, e.g., human and fish, and ASE was found to be important in the context of phenotypic variability and disease. In this study, we utilized Real-Time PCR techniques to estimate ASE of the androgen receptor alpha (arα) gene in several distinct tissues of Amazon mollies. We found an allelic bias favoring the maternal ancestor (P. mexicana) allele in ovarian tissue. This allelic bias was not observed in the gill or the brain tissue. Sequencing of the promoter regions of both alleles revealed an association between an Indel in a known CpG island and differential expression. Future studies may reveal whether our observed cis-regulatory divergence is caused by an ovary-specific trans-regulatory element, preferentially activating the allele of the maternal ancestor.

Wilson disease: At the crossroads between genetics and epigenetics-A review of the evidence

Environmental factors, including diet, exercise, stress, and toxins, profoundly impact disease phenotypes. This review examines how Wilson disease (WD), an autosomal recessive genetic disorder, is influenced by genetic and environmental inputs. WD is caused by mutations in the copper-transporter gene ATP7B, leading to the accumulation of copper in the liver and brain, resulting in hepatic, neurological, and psychiatric symptoms. These symptoms range in severity and can first appear anytime between early childhood and old age. Over 300 disease-causing mutations in ATP7B have been identified, but attempts to link genotype to the phenotypic presentation have yielded little insight, prompting investigators to identify alternative mechanisms, such as epigenetics, to explain the highly varied clinical presentation. Further, WD is accompanied by structural and functional abnormalities in mitochondria, potentially altering the production of metabolites that are required for epigenetic regulation of gene expression. Notably, environmental exposure affects the regulation of gene expression and mitochondrial function. We present the "multi-hit" hypothesis of WD progression, which posits that the initial hit is an environmental factor that affects fetal gene expression and epigenetic mechanisms and subsequent "hits" are environmental exposures that occur in the offspring after birth. These environmental hits and subsequent changes in epigenetic regulation may impact copper accumulation and ultimately WD phenotype. Lifestyle changes, including diet, increased physical activity, stress reduction, and toxin avoidance, might influence the presentation and course of WD, and therefore may serve as potential adjunctive or replacement therapies.

Microsatellite and Single Nucleotide Polymorphisms in the Insulin-Like Growth Factor 1 Promoter with Insulin Sensitivity and Insulin Secretion

Background

To investigate associations of the CA microsatellite and rs35767, rs5742612, and rs2288377 polymorphisms and the single nucleotide polymorphism (SNP) haplotypes with and without the CA microsatellite in the IGF1 promoter with insulin sensitivity and secretion.

Material/Methods

The CA microsatellite and SNPs were genotyped in 389 type 2 diabetes mellitus (T2DM) patients. A 75 g oral glucose tolerance test (OGTT) was given to all the participants. Associations of the genotypes and haplotypes with insulin sensitivity, insulin secretion, glucose tolerance, and insulin-like growth factor 1 (IGF1) were analyzed by ANCOVA (general linear model) and multiple linear regression, after controlling for gender, age, and BMI.

Results

The CA microsatellite, rs35767 polymorphisms, and SNP haplotypes with or without CA showed no significant association with metabolic parameters. The C allele of rs5742612 was found to be associated with decreased insulin sensitivity (HOMA-S index, β=−0.131, P=0.008; fasting insulin level, β=0.022, P=0.006) and increased insulin secretion (HOMA-B index, β=0.099, P=0.008; insulin AUC, β=0.112, P=0.012). The linear regression model also indicated that the A allele of rs2288377 was associated with decreased insulin sensitivity (HOMA-S index, β=−0.159, P=0.001; fasting insulin, β=0.143, P=0.001) and increased insulin secretion (HOMA-B index, β=0.114, P=0.017; insulin AUC, β=0.042, P=0.002).

Conclusions

The CA microsatellite and rs35767 have no genotype-related difference in insulin sensitivity or secretion. The rs5742612 and rs2288377 polymorphisms are significantly associated with insulin biology, with the TT genotype exhibiting higher insulin sensitivity and lower insulin secretion compared with carriers of the C allele and A allele, respectively, mostly attributed to the direct functional roles of the two loci.

Population and allelic variation of A-to-I RNA editing in human transcriptomes

Background

A-to-I RNA editing is an important step in RNA processing in which specific adenosines in some RNA molecules are post-transcriptionally modified to inosines. RNA editing has emerged as a widespread mechanism for generating transcriptome diversity. However, there remain significant knowledge gaps about the variation and function of RNA editing.

Results

In order to determine the influence of genetic variation on A-to-I RNA editing, we integrate genomic and transcriptomic data from 445 human lymphoblastoid cell lines by combining an RNA editing QTL (edQTL) analysis with an allele-specific RNA editing (ASED) analysis. We identify 1054 RNA editing events associated with cis genetic polymorphisms. Additionally, we find that a subset of these polymorphisms is linked to genome-wide association study signals of complex traits or diseases. Finally, compared to random cis polymorphisms, polymorphisms associated with RNA editing variation are located closer spatially to their respective editing sites and have a more pronounced impact on RNA secondary structure.

Conclusions

Our study reveals widespread cis variation in RNA editing among genetically distinct individuals and sheds light on possible phenotypic consequences of such variation on complex traits and diseases.

Electronic supplementary material

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

Gene expression allelic imbalance in ovine brown adipose tissue impacts energy homeostasis

Heritable trait variation within a population of organisms is largely governed by DNA variations that impact gene transcription and protein function. Identifying genetic variants that affect complex functional traits is a primary aim of population genetics studies, especially in the context of human disease and agricultural production traits. The identification of alleles directly altering mRNA expression and thereby biological function is challenging due to difficulty in isolating direct effects of cis-acting genetic variations from indirect trans-acting genetic effects. Allele specific gene expression or allelic imbalance in gene expression (AI) occurring at heterozygous loci provides an opportunity to identify genes directly impacted by cis-acting genetic variants as indirect trans-acting effects equally impact the expression of both alleles. However, the identification of genes showing AI in the context of the expression of all genes remains a challenge due to a variety of technical and statistical issues. The current study focuses on the discovery of genes showing AI using single nucleotide polymorphisms as allelic reporters. By developing a computational and statistical process that addressed multiple analytical challenges, we ranked 5,809 genes for evidence of AI using RNA-Seq data derived from brown adipose tissue samples from a cohort of late gestation fetal lambs and then identified a conservative subgroup of 1,293 genes. Thus, AI was extensive, representing approximately 25% of the tested genes. Genes associated with AI were enriched for multiple Gene Ontology (GO) terms relating to lipid metabolism, mitochondrial function and the extracellular matrix. These functions suggest that cis-acting genetic variations causing AI in the population are preferentially impacting genes involved in energy homeostasis and tissue remodelling. These functions may contribute to production traits likely to be under genetic selection in the population.

Functional Validation of a Common Nonsynonymous Coding Variant in ZC3HC1 Associated With Protection From Coronary Artery Disease

Background—

Although virtually all coronary artery disease associated single-nucleotide polymorphisms identified by genome-wide association studies (GWAS) are in noncoding regions of the genome, a common polymorphism in ZC3HC1 (rs11556924), resulting in an arginine (Arg) to histidine (His) substitution in its encoded protein, NIPA (Nuclear Interacting Partner of Anaplastic Lyphoma Kinase) is linked to a protection from coronary artery disease. NIPA plays a role in cell cycle progression, but the functional consequences of this polymorphism have not been established.

Methods and Results—

Here we demonstrate that total ZC3HC1 expression in whole blood is similar across genotypes, despite expression being slightly biased toward the risk allele in heterozygotes. At the protein level, the protective His363 NIPA variant exhibits increased phosphorylation of a critical serine residue (Ser354) and higher protein expression as compared with the Arg363 variant. Binding experiments indicate that neither SKP1 (S-phase kinase-associated protein 1) nor CCNB1 binding were affected by the polymorphism. Despite similar nuclear distribution, NIPA His363 exhibits greater nuclear mobility. NIPA suppression results in a modest reduction of proliferation in vascular smooth muscle cells, but given low proliferative capacity, a significant effect of the variant was not noted. By contrast, we demonstrate that the protective variant reduces cell proliferation in HeLa cells.

Conclusions—

These findings extend the genetic association between rs11556924 and coronary artery disease risk by characterizing its effects on the encoded protein, NIPA. The resulting amino acid change Arg363His is associated with increased expression and nuclear mobility, as well as lower rates of cell growth in HeLa cells, further supporting a role for cell proliferation in atherosclerosis and its clinical consequences.

Iron refractory iron deficiency anemia: a heterogeneous disease that is not always iron refractory

TMPRSS6 variants that affect protein function result in impaired matriptase‐2 function and consequently uninhibited hepcidin production, leading to iron refractory iron deficiency anemia (IRIDA). This disease is characterized by microcytic, hypochromic anemia and serum hepcidin values that are inappropriately high for body iron levels. Much is still unknown about its pathophysiology, genotype–phenotype correlation, and optimal clinical management. We describe 14 different TMPRSS6 variants, of which 9 are novel, in 21 phenotypically affected IRIDA patients from 20 families living in the Netherlands; 16 out of 21 patients were female. In 7 out of 21 cases DNA sequencing and multiplex ligation dependent probe amplification demonstrated only heterozygous TMPRSS6 variants. The age at presentation, disease severity, and response to iron supplementation were highly variable, even for patients and relatives with similar TMPRSS6 genotypes. Mono‐allelic IRIDA patients had a milder phenotype with respect to hemoglobin and MCV and presented significantly later in life with anemia than bi‐allelic patients. Transferrin saturation (TSAT)/hepcidin ratios were lower in IRIDA probands than in healthy relatives. Most patients required parenteral iron. Genotype alone was not predictive for the response to oral iron. We conclude that IRIDA is a genotypically and phenotypically heterogeneous disease. The high proportion of female patients and the discrepancy between phenotypes of probands and relatives with the same genotype, suggest a complex interplay between genetic and acquired factors in the pathogenesis of IRIDA. In the absence of inflammation, the TSAT/hepcidin ratio is a promising diagnostic tool, even after iron supplementation has been given. Am. J. Hematol. 91:E482–E490, 2016. © 2016 Wiley Periodicals, Inc.

Latitudinal clines in gene expression and cis-regulatory element variation in Drosophila melanogaster

BACKGROUND

Organisms can rapidly adapt to their environment when colonizing a new habitat, and this could occur by changing protein sequences or by altering patterns of gene expression. The importance of gene expression in driving local adaptation is increasingly being appreciated, and cis-regulatory elements (CREs), which control and modify the expression of the nearby genes, are predicted to play an important role. Here we investigate genetic variation in gene expression in immune-challenged Drosophila melanogaster from temperate and tropical or sub-tropical populations in Australia and United States.

RESULTS

We find parallel latitudinal changes in gene expression, with genes involved in immunity, insecticide resistance, reproduction, and the response to the environment being especially likely to differ between latitudes. By measuring allele-specific gene expression (ASE), we show that cis-regulatory variation also shows parallel latitudinal differences between the two continents and contributes to the latitudinal differences in gene expression.

CONCLUSIONS

Both Australia and United States were relatively recently colonized by D. melanogaster, and it was recently shown that introductions of both African and European flies occurred, with African genotypes contributing disproportionately to tropical populations. Therefore, both the demographic history of the populations and local adaptation may be causing the patterns that we see.

Discovering Single Nucleotide Polymorphisms Regulating Human Gene Expression Using Allele Specific Expression from RNA-seq Data

The study of the genetics of gene expression is of considerable importance to understanding the nature of common, complex diseases. The most widely applied approach to identifying relationships between genetic variation and gene expression is the expression quantitative trait loci (eQTL) approach. Here, we increased the computational power of eQTL with an alternative and complementary approach based on analyzing allele specific expression (ASE). We designed a novel analytical method to identify cis-acting regulatory variants based on genome sequencing and measurements of ASE from RNA-sequencing (RNA-seq) data. We evaluated the power and resolution of our method using simulated data. We then applied the method to map regulatory variants affecting gene expression in lymphoblastoid cell lines (LCLs) from 77 unrelated northern and western European individuals (CEU), which were part of the HapMap project. A total of 2309 SNPs were identified as being associated with ASE patterns. The SNPs associated with ASE were enriched within promoter regions and were significantly more likely to signal strong evidence for a regulatory role. Finally, among the candidate regulatory SNPs, we identified 108 SNPs that were previously associated with human immune diseases. With further improvements in quantifying ASE from RNA-seq, the application of our method to other datasets is expected to accelerate our understanding of the biological basis of common diseases.

Erasure and reestablishment of random allelic expression imbalance after epigenetic reprogramming

Clonal level random allelic expression imbalance and random monoallelic expression provides cellular heterogeneity within tissues by modulating allelic dosage. Although such expression patterns have been observed in multiple cell types, little is known about when in development these stochastic allelic choices are made. We examine allelic expression patterns in human neural progenitor cells before and after epigenetic reprogramming to induced pluripotency, observing that loci previously characterized by random allelic expression imbalance (0.63% of expressed genes) are generally reset to a biallelic state in induced pluripotent stem cells (iPSCs). We subsequently neuralized the iPSCs and profiled isolated clonal neural stem cells, observing that significant random allelic expression imbalance is reestablished at 0.65% of expressed genes, including novel loci not found to show allelic expression imbalance in the original parental neural progenitor cells. Allelic expression imbalance was associated with altered DNA methylation across promoter regulatory regions, with clones characterized by skewed allelic expression being hypermethylated compared to their biallelic sister clones. Our results suggest that random allelic expression imbalance is established during lineage commitment and is associated with increased DNA methylation at the gene promoter.

Prevalent Accumulation of Non-Optimal Codons through Somatic Mutations in Human Cancers

Cancer is characterized by uncontrolled cell growth, and the cause of different cancers is generally attributed to checkpoint dysregulation of cell proliferation and apoptosis. Recent studies have shown that non-optimal codons were preferentially adopted by genes to generate cell cycle-dependent oscillations in protein levels. This raises the intriguing question of how dynamic changes of codon usage modulate the cancer genome to cope with a non-controlled proliferative cell cycle. In this study, we comprehensively analyzed the somatic mutations of codons in human cancers, and found that non-optimal codons tended to be accumulated through both synonymous and non-synonymous mutations compared with other types of genomic substitution. We further demonstrated that non-optimal codons were prevalently accumulated across different types of cancers, amino acids, and chromosomes, and genes with accumulation of non-optimal codons tended to be involved in protein interaction/signaling networks and encoded important enzymes in metabolic networks that played roles in cancer-related pathways. This study provides insights into the dynamics of codons in the cancer genome and demonstrates that accumulation of non-optimal codons may be an adaptive strategy for cancerous cells to win the competition with normal cells. This deeper interpretation of the patterns and the functional characterization of somatic mutations of codons will help to broaden the current understanding of the molecular basis of cancers.

Effects of IL1B single nucleotide polymorphisms on depressive and anxiety symptoms are determined by severity and type of life stress

Interleukin-1β is one of the main mediators in the cross-talk between the immune system and the central nervous system. Higher interleukin-1β levels are found in mood spectrum disorders, and the stress-induced expression rate of the interleukin-1β gene (IL1B) is altered by polymorphisms in the region. Therefore we examined the effects of rs16944 and rs1143643 single nucleotide polymorphisms (SNPs) within the IL1B gene on depressive and anxiety symptoms, as measured by the Brief Symptom Inventory, in a Hungarian population sample of 1053 persons. Distal and proximal environmental stress factors were also included in our analysis, namely childhood adversity and recent negative life-events. We found that rs16944 minor (A) allele specifically interacted with childhood adversity increasing depressive and anxiety symptoms, while rs1143643's minor (A) allele showed protective effect against depressive symptoms after recent life stress. The genetic main effects of the two SNPs were not significant in the main analysis, but the interaction effects remained significant after correction for multiple testing. In addition, the effect of rs16944 A allele was reversed in a subsample with low-exposure to life stress, suggesting a protective effect against depressive symptoms, in the post hoc analysis. In summary, both of the two IL1B SNPs showed specific environmental stressor-dependent effects on mood disorder symptoms. We also demonstrated that the presence of exposure to childhood adversity changed the direction of the rs16944 effect on depression phenotype. Therefore our results suggest that it is advisable to include environmental factors in genetic association studies when examining the effect of the IL1B gene.

Genome-wide association study of behavioral, physiological and gene expression traits in outbred CFW mice

Although mice are the most widely used mammalian model organism, genetic studies have suffered from limited mapping resolution due to extensive linkage disequilibrium (LD) that is characteristic of crosses among inbred strains. Carworth Farms White (CFW) mice are a commercially available outbred mouse population that exhibit rapid LD decay in comparison to other available mouse populations. We performed a genome-wide association study (GWAS) of behavioral, physiological and gene expression phenotypes using 1,200 male CFW mice. We used genotyping by sequencing (GBS) to obtain genotypes at 92,734 SNPs. We also measured gene expression using RNA sequencing in three brain regions. Our study identified numerous behavioral, physiological and expression quantitative trait loci (QTLs). We integrated the behavioral QTL and eQTL results to implicate specific genes, including Azi2 in sensitivity to methamphetamine and Zmynd11 in anxiety-like behavior. The combination of CFW mice, GBS and RNA sequencing constitutes a powerful approach to GWAS in mice.

Analysis of the porcine APOA2 gene expression in liver, polymorphism identification and association with fatty acid composition traits

APOA2 is a protein implicated in triglyceride, fatty acid and glucose metabolism. In pigs, the APOA2 gene is located on pig chromosome 4 (SSC4) in a QTL region affecting fatty acid composition, fatness and growth traits. In this study, we evaluated APOA2 as a candidate gene for meat quality traits in an Iberian × Landrace backcross population. The APOA2:c.131T>A polymorphism, located in exon 3 of APOA2 and determining a missense mutation, was associated with the percentage of hexadecenoic acid [C16:1(n-9)], linoleic acid [C18:2(n-6)], α-linolenic acid [C18:3(n-3)], dihomo-gamma-linolenic acid [C20:3(n-6)] and polyunsaturated fatty acids (PUFAs) in backfat. Furthermore, this SNP was associated with the global mRNA expression levels of APOA2 in liver and was used as a marker to determine allelic expression imbalance by pyrosequencing. We determined an overexpression of the T allele in heterozygous samples with a mean ratio of 2.8 (T/A), observing a high variability in the allelic expression among individuals. This result suggests that complex regulatory mechanisms, beyond a single polymorphism (e.g. epigenetic effects or multiple cis-acting polymorphisms), may be regulating APOA2 gene expression.

DNA methylation of loci within ABCG1 and PHOSPHO1 in blood DNA is associated with future type 2 diabetes risk

Identification of subjects with a high risk of developing type 2 diabetes (T2D) is fundamental for prevention of the disease. Consequently, it is essential to search for new biomarkers that can improve the prediction of T2D. The aim of this study was to examine whether 5 DNA methylation loci in blood DNA (ABCG1, PHOSPHO1, SOCS3, SREBF1, and TXNIP), recently reported to be associated with T2D, might predict future T2D in subjects from the Botnia prospective study. We also tested if these CpG sites exhibit altered DNA methylation in human pancreatic islets, liver, adipose tissue, and skeletal muscle from diabetic vs. non-diabetic subjects. DNA methylation at the ABCG1 locus cg06500161 in blood DNA was associated with an increased risk for future T2D (OR = 1.09, 95% CI = 1.02–1.16, P-value = 0.007, Q-value = 0.018), while DNA methylation at the PHOSPHO1 locus cg02650017 in blood DNA was associated with a decreased risk for future T2D (OR = 0.85, 95% CI = 0.75–0.95, P-value = 0.006, Q-value = 0.018) after adjustment for age, gender, fasting glucose, and family relation. Furthermore, the level of DNA methylation at the ABCG1 locus cg06500161 in blood DNA correlated positively with BMI, HbA1c, fasting insulin, and triglyceride levels, and was increased in adipose tissue and blood from the diabetic twin among monozygotic twin pairs discordant for T2D. DNA methylation at the PHOSPHO1 locus cg02650017 in blood correlated positively with HDL levels, and was decreased in skeletal muscle from diabetic vs. non-diabetic monozygotic twins. DNA methylation of cg18181703 (SOCS3), cg11024682 (SREBF1), and cg19693031 (TXNIP) was not associated with future T2D risk in subjects from the Botnia prospective study.

The interface between genetics and psychology: lessons from developmental dyslexia

Developmental dyslexia runs in families, and twin studies have confirmed that there is a substantial genetic contribution to poor reading. The way in which discoveries in molecular genetics are reported can be misleading, encouraging us to think that there are specific genes that might be used to screen for disorder. However, dyslexia is not a classic Mendelian disorder that is caused by a mutation in a single gene. Rather, like many other common disorders, it appears to involve combined effects of many genes and environmental factors, each of which has a small influence, possibly supplemented by rare variants that have larger effects but apply to only a minority of cases. Furthermore, to see clearer relationships between genotype and phenotype, we may need to move beyond the clinical category of dyslexia to look at underlying cognitive deficits that may be implicated in other neurodevelopmental disorders.

GeneiASE: Detection of condition-dependent and static allele-specific expression from RNA-seq data without haplotype information

Allele-specific expression (ASE) is the imbalance in transcription between maternal and paternal alleles at a locus and can be probed in single individuals using massively parallel DNA sequencing technology. Assessing ASE within a single sample provides a static picture of the ASE, but the magnitude of ASE for a given transcript may vary between different biological conditions in an individual. Such condition-dependent ASE could indicate a genetic variation with a functional role in the phenotypic difference. We investigated ASE through RNA-sequencing of primary white blood cells from eight human individuals before and after the controlled induction of an inflammatory response, and detected condition-dependent and static ASE at 211 and 13021 variants, respectively. We developed a method, GeneiASE, to detect genes exhibiting static or condition-dependent ASE in single individuals. GeneiASE performed consistently over a range of read depths and ASE effect sizes, and did not require phasing of variants to estimate haplotypes. We observed condition-dependent ASE related to the inflammatory response in 19 genes, and static ASE in 1389 genes. Allele-specific expression was confirmed by validation of variants through real-time quantitative RT-PCR, with RNA-seq and RT-PCR ASE effect-size correlations r = 0.67 and r = 0.94 for static and condition-dependent ASE, respectively.

GABRA2 Alcohol Dependence Risk Allele is Associated with Reduced Expression of Chromosome 4p12 GABAA Subunit Genes in Human Neural Cultures

BACKGROUND

Genetic variation in a region of chromosome 4p12 that includes the GABAA subunit gene GABRA2 has been reproducibly associated with alcohol dependence (AD). However, the molecular mechanisms underlying the association are unknown. This study examined correlates of in vitro gene expression of the AD-associated GABRA2 rs279858*C-allele in human neural cells using an induced pluripotent stem cell (iPSC) model system.

METHODS

We examined mRNA expression of chromosome 4p12 GABAA subunit genes (GABRG1, GABRA2, GABRA4, and GABRB1) in 36 human neural cell lines differentiated from iPSCs using quantitative polymerase chain reaction and next-generation RNA sequencing. mRNA expression in adult human brain was examined using the BrainCloud and BRAINEAC data sets.

RESULTS

We found significantly lower levels of GABRA2 mRNA in neural cell cultures derived from rs279858*C-allele carriers. Levels of GABRA2 RNA were correlated with those of the other 3 chromosome 4p12 GABAA genes, but not other neural genes. Cluster analysis based on the relative RNA levels of the 4 chromosome 4p12 GABAA genes identified 2 distinct clusters of cell lines, a low-expression cluster associated with rs279858*C-allele carriers and a high-expression cluster enriched for the rs279858*T/T genotype. In contrast, there was no association of genotype with chromosome 4p12 GABAA gene expression in postmortem adult cortex in either the BrainCloud or BRAINEAC data sets.

CONCLUSIONS

AD-associated variation in GABRA2 is associated with differential expression of the entire cluster of GABAA subunit genes on chromosome 4p12 in human iPSC-derived neural cell cultures. The absence of a parallel effect in postmortem human adult brain samples suggests that AD-associated genotype effects on GABAA expression, although not present in mature cortex, could have effects on regulation of the chromosome 4p12 GABAA cluster during neural development.

Analyzing allele specific RNA expression using mixture models

Background

Measuring allele-specific RNA expression provides valuable insights into cis-acting genetic and epigenetic regulation of gene expression. Widespread adoption of high-throughput sequencing technologies for studying RNA expression (RNA-Seq) permits measurement of allelic RNA expression imbalance (AEI) at heterozygous single nucleotide polymorphisms (SNPs) across the entire transcriptome, and this approach has become especially popular with the emergence of large databases, such as GTEx. However, the existing binomial-type methods used to model allelic expression from RNA-seq assume a strong negative correlation between reference and variant allele reads, which may not be reasonable biologically.

Results

Here we propose a new strategy for AEI analysis using RNA-seq data. Under the null hypothesis of no AEI, a group of SNPs (possibly across multiple genes) is considered comparable if their respective total sums of the allelic reads are of similar magnitude. Within each group of “comparable” SNPs, we identify SNPs with AEI signal by fitting a mixture of folded Skellam distributions to the absolute values of read differences. By applying this methodology to RNA-Seq data from human autopsy brain tissues, we identified numerous instances of moderate to strong imbalanced allelic RNA expression at heterozygous SNPs. Findings with SLC1A3 mRNA exhibiting known expression differences are discussed as examples.

Conclusion

The folded Skellam mixture model searches for SNPs with significant difference between reference and variant allele reads (adjusted for different library sizes), using information from a group of “comparable” SNPs across multiple genes. This model is particularly suitable for performing AEI analysis on genes with few heterozygous SNPs available from RNA-seq, and it can fit over-dispersed read counts without specifying the direction of the correlation between reference and variant alleles.

Electronic supplementary material

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

Allelome.PRO, a pipeline to define allele-specific genomic features from high-throughput sequencing data

Detecting allelic biases from high-throughput sequencing data requires an approach that maximises sensitivity while minimizing false positives. Here, we present Allelome.PRO, an automated user-friendly bioinformatics pipeline, which uses high-throughput sequencing data from reciprocal crosses of two genetically distinct mouse strains to detect allele-specific expression and chromatin modifications. Allelome.PRO extends approaches used in previous studies that exclusively analyzed imprinted expression to give a complete picture of the ‘allelome’ by automatically categorising the allelic expression of all genes in a given cell type into imprinted, strain-biased, biallelic or non-informative. Allelome.PRO offers increased sensitivity to analyze lowly expressed transcripts, together with a robust false discovery rate empirically calculated from variation in the sequencing data. We used RNA-seq data from mouse embryonic fibroblasts from F1 reciprocal crosses to determine a biologically relevant allelic ratio cutoff, and define for the first time an entire allelome. Furthermore, we show that Allelome.PRO detects differential enrichment of H3K4me3 over promoters from ChIP-seq data validating the RNA-seq results. This approach can be easily extended to analyze histone marks of active enhancers, or transcription factor binding sites and therefore provides a powerful tool to identify candidate cis regulatory elements genome wide.

Genetic variants in arachidonic acid pathway genes associated with NSAID-exacerbated respiratory disease

AIM

NSAIDs are the most frequent cause of hypersensitivity drug reactions. We have examined the association between NSAID-exacerbated respiratory disease (NERD) and genetic variants in arachidonic acid metabolism genes.

PATIENTS & METHODS

We included 250 NERD patients, 260 NSAID-tolerant asthmatic (NTA) subjects and 315 healthy controls.

RESULTS

Significant associations with NERD were identified for: ALOX15 rs3892408 C/C homozygous genotype (NERD vs NTA; p = 0.0001, pc = 0.0011; NERD vs controls; p = 0.0001, pc = 0.0011), PTGS-1 rs5789 A/A homozygous genotype (NERD vs NTA; p = 0.0001, pc = 0.0011; NERD vs controls; p = 0.0001, pc = 0.0011), PTGS-1 rs10306135 A/A homozygous genotype (NERD vs NTA; p = 0.0009, pc = 0.0091; NERD vs controls; p = 0.0064, pc = 0.045). Differences in ALOX5 copy number variations were also found (NERD vs NTA; p = 0.010; NERD vs controls; p = 0.0001).

CONCLUSION

These results improve our understanding of the underlying mechanisms of NERD and may help develop a predictive test for this pathology. Original submitted 3 November 2014; Revision submitted 2 April 2015.

Genomic imprinting: A missing piece of the Multiple Sclerosis puzzle?

Evidence for parent-of-origin effects in complex diseases such as Multiple Sclerosis (MS) strongly suggests a role for epigenetic mechanisms in their pathogenesis. In this review, we describe the importance of accounting for parent-of-origin when identifying new risk variants for complex diseases and discuss how genomic imprinting, one of the best-characterized epigenetic mechanisms causing parent-of-origin effects, may impact etiology of complex diseases. While the role of imprinted genes in growth and development is well established, the contribution and molecular mechanisms underlying the impact of genomic imprinting in immune functions and inflammatory diseases are still largely unknown. Here we discuss emerging roles of imprinted genes in the regulation of inflammatory responses with a particular focus on the Dlk1 cluster that has been implicated in etiology of experimental MS-like disease and Type 1 Diabetes. Moreover, we speculate on the potential wider impact of imprinting via the action of imprinted microRNAs, which are abundantly present in the Dlk1 locus and predicted to fine-tune important immune functions. Finally, we reflect on how unrelated imprinted genes or imprinted genes together with non-imprinted genes can interact in so-called imprinted gene networks (IGN) and suggest that IGNs could partly explain observed parent-of-origin effects in complex diseases. Unveiling the mechanisms of parent-of-origin effects is therefore likely to teach us not only about the etiology of complex diseases but also about the unknown roles of this fascinating phenomenon underlying uneven genetic contribution from our parents. This article is part of a Directed Issue entitled: Epigenetics dynamics in development and disease.

Polymorphisms in endothelin system genes, arsenic levels and obesity risk

BACKGROUND/OBJECTIVES

Obesity has been linked to morbidity and mortality through increased risk for many chronic diseases. Endothelin (EDN) system has been related to endothelial function but it can be involved in lipid metabolism regulation: Receptor type A (EDNRA) activates lipolysis in adipocytes, the two endothelin receptors mediate arsenic-stimulated adipocyte dysfunction, and endothelin system can regulate adiposity by modulating adiponectin activity in different situations and, therefore, influence obesity development. The aim of the present study was to analyze if single nucleotide polymorphisms (SNPs) in the EDN system could be associated with human obesity.

SUBJECTS/METHODS

We analyzed two samples of general-population-based studies from two different regions of Spain: the VALCAR Study, 468 subjects from the area of Valencia, and the Hortega Study, 1502 subjects from the area of Valladolid. Eighteen SNPs throughout five genes were analyzed using SNPlex.

RESULTS

We found associations for two polymorphisms of the EDNRB gene which codifies for EDN receptor type B. Genotypes AG and AA of the rs5351 were associated with a lower risk for obesity in the VALCAR sample (p=0.048, OR=0.63) and in the Hortega sample (p=0.001, OR=0.62). Moreover, in the rs3759475 polymorphism, genotypes CT and TT were also associated with lower risk for obesity in the Hortega sample (p=0.0037, OR=0.66) and in the VALCAR sample we found the same tendency (p=0.12, OR=0.70). Furthermore, upon studying the pooled population, we found a stronger association with obesity (p=0.0001, OR=0.61 and p=0.0008, OR=0.66 for rs5351 and rs3759475, respectively). Regarding plasma arsenic levels, we have found a positive association for the two SNPs studied with obesity risk in individuals with higher arsenic levels in plasma: rs5351 (p=0.0054, OR=0.51) and rs3759475 (p=0.009, OR=0.53).

CONCLUSIONS

Our results support the hypothesis that polymorphisms of the EDNRB gene may influence the susceptibility to obesity and can interact with plasma arsenic levels.

Differential allelic expression of SOS1 and hyperexpression of the activating SOS1 c.755C variant in a Noonan syndrome family

Noonan syndrome (NS) is a genetic condition characterized by congenital heart defects, short stature and characteristic facial features. We here present the case of a girl with moderate learning disabilities, delayed language development, craniofacial features and skin anomalies reminiscent of NS. After a mutation screening of the known NS genes PTPN11, SOS1, RAF1, KRAS, GRB2, BRAF and SHOC2 we found the heterozygous c.755T>C variant in SOS1 causing the p.I252T amino-acid substitution, which was considered possibly pathogenetic by bioinformatic predictions. The same variant was present in the proband's mother, displaying some NS features, and maternal grandfather showing no NS traits, but also by a healthy subject in 1000 genomes project database without phenotype informations. The functional analysis revealed that SOS1 c.755C activated the RAS-ERK intracellular pathway, whereas no effects on RAC-JNK cascade have been detected. After a comparison between the sequence of SOS1 cDNA from peripheral blood and SOS1 genomic DNA, we showed for the first time a differential allelic expression of the SOS1 gene in healthy individuals, thus occurring as a physiologic condition. Interestingly, we found that the mutated allele C was 50% more expressed than the wild-type allele T in all familial carriers. The comparable amount of SOS1 mRNA between mutated individuals and the controls indicates that the variant does not affect SOS1 expression. The present study provides a first evidence of allelic imbalance of SOS1 and pinpoints this condition as a possible mechanism underlying a different penetrance of some SOS1-mutated alleles in unrelated carriers.

The clinical presentation of Marfan syndrome is modulated by expression of wild-type FBN1 allele

Marfan syndrome is an autosomal dominant disorder mainly caused by mutations within FBN1 gene. The disease displays large variability in age of onset or severity and very poor phenotype/genotype correlations have been demonstrated. We investigated the hypothesis that phenotype severity could be related to the variable expression level of fibrillin-1 (FBN1) synthesized from the wild-type (WT) allele. Quantitative reverse-transcription and polymerase chain reaction was used to evaluate FBN1 levels in skin fibroblasts from 80 Marfan patients with premature termination codons and in skin fibroblasts from 80 controls. Results in controls showed a 3.9-fold variation in FBN1 mRNA synthesis level between subjects. A similar 4.4-fold variation was found in the Marfan population, but the mean level of FBN1 mRNA was a half of the control population. Differential allelic expression analysis in Marfan fibroblasts showed that over 90% of FBN1 mRNA was transcribed from the wild allele and the mutated allele was not detected. In the control population, independently of the expression level of FBN1, we observed steady-state equilibrium between the two allelic-mRNAs suggesting that FBN1 expression mainly depends on trans-acting regulators. Finally, we show that a low level of residual WT FBN1 mRNA accounts for a high risk of ectopia lentis and pectus abnormality and tends to increase the risk of aortic dilatation.