The Defects of Epigenetic Reprogramming in Dox-Dependent Porcine-iPSCs.

Porcine-induced pluripotent stem cells (piPSCs) are of great significance to animal breeding and human medicine; however, an important problem is that the maintenance of piPSCs mainly depends on exogenous expression of pluripotent transcription factors (TFs), and germline transmission-competent piPSCs have not yet been successfully established. In this study, we explore the defect of epigenetic reprogramming during piPSCs formation, including chromatin accessibility, DNA methylation, and imprinted gene expression, with high-throughput sequencing (ATAC-seq, WGBS, RNA-seq, and Re-seq) methods. We found the somatic features were successfully silenced by connecting closed chromatin loci with downregulated genes, while DNA methylation has limited effects on somatic silence. However, the incomplete chromatin remodeling and DNA demethylation in pluripotency genes hinder pluripotent activation, resulting in the low expression of endogenous pluripotency genes. In addition, the expression of potential imprinted genes was abnormal, and many allelic-biased expressed genes in porcine embryonic fibroblasts (PEFs) were erased, accompanied by establishment of new allelic-biased expressed genes in piPSCs. This study reveals the aberrant epigenetic reprogramming during dox-dependent piPSCs formation, which lays the foundation for research of porcine-iPSC reprogramming and genome imprinting.

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

Abstract 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.

The influence of FKBP5 genotype on expression of FKBP5 and other glucocorticoid-regulated genes, dependent on trauma exposure.

The FK506 binding protein 51 (FKBP5), an intrinsic regulator of the glucocorticoid receptor, has been associated with pathological behaviors particularly in the context of childhood trauma (CT), via a putatively regulatory polymorphism, rs1360780. However, trans- and cis-acting effects of this locus and its interaction with CT are incompletely understood. To study its effects on the expression of glucocorticoid-regulated genes including FKBP5, we used lymphoblastoid cell lines (LCLs) derived from 16 CT-exposed patients with greater than two substance dependence/suicidal behavior diagnoses (casesCT+) and 13 non-CT-exposed controls (controlsCT-). This study in LCLs measures long-term trait-like differences attributable to genotype or lasting epigenetic modification. Through analysis of differential allelic expression (DAE) using an FKBP5 3'-UTR reporter single nucleotide polymorphism (SNP), rs3800373, that is in strong linkage disequilibrium with rs1360780, we confirmed that the rs1360780 risk allele (A) (or conceivably that of a linked SNP) leads to higher FKBP5 expression in controlsCT-. Intriguingly, casesCT+ did not show DAE, perhaps because of a genotype-predicted difference in FKBP5 DNA methylation restricted to casesCT+. Furthermore, through correlation analyses on FKBP5 expression at baseline and after induction by dexamethasone, we observed that casesCT+ had lower induction of FKBP5 expression, indicating that overall they may have strong ultra-short negative-feedback. Only casesCT+ showed an effect of rs1360780 genotype on expression of FKBP5 and other glucocorticoid-regulated genes. Together, these results confirm that the rs1360780 locus alters FKBP5 expression and further that in trans-fashion this locus affects the expression of other glucocorticoid-regulated genes after a glucocorticoid challenge. The CT exposure appears to be essential for trans-effects of rs1360780 on glucocorticoid-regulated genes.

Association of breast cancer risk in BRCA1 and BRCA2 mutation carriers with genetic variants showing differential allelic expression: identification of a modifier of breast cancer risk at locus 11q22.3.

PURPOSE: Cis-acting regulatory SNPs resulting in differential allelic expression (DAE) may, in part, explain the underlying phenotypic variation associated with many complex diseases. To investigate whether common variants associated with DAE were involved in breast cancer susceptibility among BRCA1 and BRCA2 mutation carriers, a list of 175 genes was developed based of their involvement in cancer-related pathways. METHODS: Using data from a genome-wide map of SNPs associated with allelic expression, we assessed the association of ~320 SNPs located in the vicinity of these genes with breast and ovarian cancer risks in 15,252 BRCA1 and 8211 BRCA2 mutation carriers ascertained from 54 studies participating in the Consortium of Investigators of Modifiers of BRCA1/2. RESULTS: We identified a region on 11q22.3 that is significantly associated with breast cancer risk in BRCA1 mutation carriers (most significant SNP rs228595 p = 7 × 10-6). This association was absent in BRCA2 carriers (p = 0.57). The 11q22.3 region notably encompasses genes such as ACAT1, NPAT, and ATM. Expression quantitative trait loci associations were observed in both normal breast and tumors across this region, namely for ACAT1, ATM, and other genes. In silico analysis revealed some overlap between top risk-associated SNPs and relevant biological features in mammary cell data, which suggests potential functional significance. CONCLUSION: We identified 11q22.3 as a new modifier locus in BRCA1 carriers. Replication in larger studies using estrogen receptor (ER)-negative or triple-negative (i.e., ER-, progesterone receptor-, and HER2-negative) cases could therefore be helpful to confirm the association of this locus with breast cancer risk.

Association of breast cancer risk with genetic variants showing differential allelic expression: Identification of a novel breast cancer susceptibility locus at 4q21.

There are significant inter-individual differences in the levels of gene expression. Through modulation of gene expression, cis-acting variants represent an important source of phenotypic variation. Consequently, cis-regulatory SNPs associated with differential allelic expression are functional candidates for further investigation as disease-causing variants. To investigate whether common variants associated with differential allelic expression were involved in breast cancer susceptibility, a list of genes was established on the basis of their involvement in cancer related pathways and/or mechanisms. Thereafter, using data from a genome-wide map of allelic expression associated SNPs, 313 genetic variants were selected and their association with breast cancer risk was then evaluated in 46,451 breast cancer cases and 42,599 controls of European ancestry ascertained from 41 studies participating in the Breast Cancer Association Consortium. The associations were evaluated with overall breast cancer risk and with estrogen receptor negative and positive disease. One novel breast cancer susceptibility locus on 4q21 (rs11099601) was identified (OR = 1.05, P = 5.6x10-6). rs11099601 lies in a 135 kb linkage disequilibrium block containing several genes, including, HELQ, encoding the protein HEL308 a DNA dependant ATPase and DNA Helicase involved in DNA repair, MRPS18C encoding the Mitochondrial Ribosomal Protein S18C and FAM175A (ABRAXAS), encoding a BRCA1 BRCT domain-interacting protein involved in DNA damage response and double-strand break (DSB) repair. Expression QTL analysis in breast cancer tissue showed rs11099601 to be associated with HELQ (P = 8.28x10-14), MRPS18C (P = 1.94x10-27) and FAM175A (P = 3.83x10-3), explaining about 20%, 14% and 1%, respectively of the variance inexpression of these genes in breast carcinomas.

The abundance of cis-acting loci leading to differential allele expression in F1 mice and their relationship to loci harboring genes affecting complex traits.

BACKGROUND: Genome-wide surveys have detected cis-acting quantitative trait loci altering levels of RNA transcripts (RNA-eQTLs) by associating SNV alleles to transcript levels. However, the sensitivity and specificity of detection of cis- expression quantitative trait loci (eQTLs) by genetic approaches, reliant as it is on measurements of transcript levels in recombinant inbred strains or offspring from arranged crosses, is unknown, as is their relationship to QTL's for complex phenotypes. RESULTS: We used transcriptome-wide differential allele expression (DAE) to detect cis-eQTLs in forebrain and kidney from reciprocal crosses between three mouse inbred strains, 129S1/SvlmJ, DBA/2J, and CAST/EiJ and C57BL/6 J. Two of these crosses were previously characterized for cis-eQTLs and QTLs for various complex phenotypes by genetic analysis of recombinant inbred (RI) strains. 5.4 %, 1.9 % and 1.5 % of genes assayed in forebrain of B6/129SF1, B6/DBAF1, and B6/CASTF1 mice, respectively, showed differential allelic expression, indicative of cis-acting alleles at these genes. Moreover, the majority of DAE QTLs were observed to be tissue-specific with only a small fraction showing cis-effects in both tissues. Comparing DAE QTLs in F1 mice to cis-eQTLs previously mapped in RI strains we observed that many of the cis-eQTLs were not confirmed by DAE. Additionally several novel DAE-QTLs not identified as cis-eQTLs were identified suggesting that there are differences in sensitivity and specificity for QTL detection between the two methodologies. Strain specific DAE QTLs in B6/DBAF1 mice were located in excess at candidate genes for alcohol use disorders, seizures, and angiogenesis previously implicated by genetic linkage in C57BL/6J × DBA/2JF2 mice or BXD RI strains. CONCLUSIONS: Via a survey for differential allele expression in F1 mice, a substantial proportion of genes were found to have alleles altering expression in cis-acting fashion. Comparing forebrain and kidney, many or most of these alleles were tissue-specific in action. The identification of strain specific DAE QTLs, can assist in assessment of candidate genes located within the large intervals associated with trait QTLs.

Differential Allelic Expression of HTR1B in Suicide Victims: Genetic and Epigenetic Effect of the Cis-Acting Variants.

OBJECTIVES: In the present study, we tested the allelic imbalance of the C861G single nucleotide polymorphism (SNP) of HTR1B in the frontal cortex of suicide victims. METHODS: The study was conducted using 3 sets of samples. First, C861G allele-specific mRNA levels in the frontal cortex were compared between suicide (n = 13) and nonsuicide controls (n = 13) from the Stanley Medical Research postmortem brain collection. Second, we tested common variants in the HTR1B promoter for linkage disequilibrium (LD) with the C861G variant in an unrelated sample of suicide attempters (SA; n = 38) and non-SA (NSA; n = 42). Finally, we performed a family-based association study of the C861G and promoter variants in 162 nuclear families using suicidal behavior severity scores as phenotype. RESULTS: We observed no alterations in the C/G expression ratio in suicide victims compared to nonsuicide controls (p = 0.370). When comparing the LD between the C861G and cis-acting SNPs, we did not find any differences in SA and NSA. There was no association between preferential transmission of cis-acting SNPs and suicidal behavior severity scores in both maternal and paternal meiosis. CONCLUSIONS: We found several promoter variants in LD that may potentially influence the allelic imbalance in the C861G variant. However, no evidence of allelic imbalance nor parent-of-origin effects of the C861G variant was observed in suicidal behavior. Further research is required to assess this marker in larger cohorts.

Overall and allele-specific expression of the SMC1A gene in female Cornelia de Lange syndrome patients and healthy controls.

Cornelia de Lange syndrome (CdLS) is a rare multisystem disorder characterized by facial dysmorphisms, limb anomalies, and growth and cognitive deficits. Mutations in genes encoding subunits (SMC1A, SMC3, RAD21) or regulators (NIPBL, HDAC8) of the cohesin complex account for approximately 65% of clinically diagnosed CdLS cases. The SMC1A gene (Xp11.22), responsible for 5% of CdLS cases, partially escapes X chromosome inactivation in humans and the allele on the inactive X chromosome is variably expressed. In this study, we evaluated overall and allele-specific SMC1A expression. Real-time PCR analysis conducted on 17 controls showed that SMC1A expression in females is 50% higher than in males. Immunoblotting experiments confirmed a 44% higher protein level in healthy females than in males, and showed no significant differences in SMC1A protein levels between controls and patients. Pyrosequencing was used to assess the reciprocal level of allelic expression in six female carriers of different SMC1A mutations and 15 controls who were heterozygous at a polymorphic transcribed SMC1A locus. The two alleles were expressed at a 1:1 ratio in the control group and at a 2:1 ratio in favor of the wild type allele in the test group. Since a dominant negative effect is considered the pathogenic mechanism in SMC1A-defective female patients, the level of allelic preferential expression might be one of the factors contributing to the wide phenotypic variability observed in these patients. An extension of this study to a larger cohort containing mild to borderline cases could enhance our understanding of the clinical spectrum of SMC1A-linked CdLS.

Allelic differences in bovine kappa - CN gene which may regulate gene expression.

The most common kappa casein (κ-CN) variants, κ-CN A and κ-CN B, are synthesised differentially in the lactating mammary gland of heterozygous animals (κ-CN AB). In this study we evaluated several approaches for quantification of allele specific mRNA transcripts. The most consistent results were obtained using allele specific RT-PCR and capillary electrophoresis. On average, 13.4% more allele B specific than A specific transcripts were found. DNA sequencing of the proximal promoter region in several homozygous animals (κ-CN AA, BB, EE) did not reveal any allele specific polymorphisms. Using the EMSA and DNase I footprinting we confirmed functional binding sites for three transcription factors (AP-2, NF1 and MGF) within the κ-CN proximal promoter region. Sequence analysis of the 3'-UTR of the κ-CN gene revealed seven allele specific sites. Two of these allelic differences were close to previously identified 3'-end regulatory sequences. In addition, allele specific differences in length between mRNAs of both variants were found. The two later findings suggest a possible post translational control determining content differences of κ-CN in milk.