Wibowo Arindrarto

Disease variants alter transcription factor levels and methylation of their binding sites

Most disease-associated genetic variants are noncoding, making it challenging to design experiments to understand their functional consequences. Identification of expression quantitative trait loci (eQTLs) has been a powerful approach to infer the downstream effects of disease-associated variants, but most of these variants remain unexplained. The analysis of DNA methylation, a key component of the epigenome, offers highly complementary data on the regulatory potential of genomic regions. Here we show that disease-associated variants have widespread effects on DNA methylation in trans that likely reflect differential occupancy of trans binding sites by cis-regulated transcription factors. Using multiple omics data sets from 3,841 Dutch individuals, we identified 1,907 established trait-associated SNPs that affect the methylation levels of 10,141 different CpG sites in trans (false discovery rate (FDR) < 0.05). These included SNPs that affect both the expression of a nearby transcription factor (such as NFKB1, CTCF and NKX2-3) and methylation of its respective binding site across the genome. Trans methylation QTLs effectively expose the downstream effects of disease-associated variants.

Identification of context-dependent expression quantitative trait loci in whole blood

Genetic risk factors often localize to noncoding regions of the genome with unknown effects on disease etiology. Expression quantitative trait loci (eQTLs) help to explain the regulatory mechanisms underlying these genetic associations. Knowledge of the context that determines the nature and strength of eQTLs may help identify cell types relevant to pathophysiology and the regulatory networks underlying disease. Here we generated peripheral blood RNA–seq data from 2,116 unrelated individuals and systematically identified context-dependent eQTLs using a hypothesis-free strategy that does not require previous knowledge of the identity of the modifiers. Of the 23,060 significant cis-regulated genes (false discovery rate (FDR) ≤ 0.05), 2,743 (12%) showed context-dependent eQTL effects. The majority of these effects were influenced by cell type composition. A set of 145 cis-eQTLs depended on type I interferon signaling. Others were modulated by specific transcription factors binding to the eQTL SNPs.

KeyGenes, a Tool to Probe Tissue Differentiation Using a Human Fetal Transcriptional Atlas

Summary Differentiated derivatives of human pluripotent stem cells in culture are generally phenotypically immature compared to their adult counterparts. Their identity is often difficult to determine with certainty because little is known about their human fetal equivalents in vivo. Cellular identity and signaling pathways directing differentiation are usually determined by extrapolating information from either human adult tissue or model organisms, assuming conservation with humans. To resolve this, we generated a collection of human fetal transcriptional profiles at different developmental stages. Moreover, we developed an algorithm, KeyGenes, which uses this dataset to quantify the extent to which next-generation sequencing or microarray data resemble specific cell or tissue types in the human fetus. Using KeyGenes combined with the human fetal atlas, we identified multiple cell and tissue samples unambiguously on a limited set of features. We thus provide a flexible and expandable platform to monitor and evaluate the efficiency of differentiation in vitro.

Uncoupling protein 2 gene polymorphisms are associated with obesity

BackgroundUncoupling protein 2 (UCP2) gene polymorphisms have been reported as genetic risk factors for obesity and type 2 diabetes mellitus (T2DM). We examined the association of commonly observed UCP2 G(−866)A (rs659366) and Ala55Val (C > T) (rs660339) single nucleotide polymorphisms (SNPs) with obesity, high fasting plasma glucose, and serum lipids in a Balinese population.MethodsA total of 603 participants (278 urban and 325 rural subjects) were recruited from Bali Island, Indonesia. Fasting plasma glucose (FPG), triglyceride (TG), high density lipoprotein cholesterol (HDL-C), low density lipoprotein cholesterol (LDL-C) and total cholesterol (TC) were measured. Obesity was determined based on WHO classifications for adult Asians. Participants were genotyped for G(−866)A and Ala55Val polymorphisms of the UCP2 gene.ResultsObesity prevalence was higher in urban subjects (51%) as compared to rural subjects (23%). The genotype, minor allele (MAF), and heterozygosity frequencies were similar between urban and rural subjects for both SNPs. All genotype frequencies were in Hardy-Weinberg equilibrium. A combined analysis of genotypes and environment revealed that the urban subjects carrying the A/A genotype of the G(−866)A SNP have higher BMI than the rural subjects with the same genotype. Since the two SNPs showed strong linkage disequilibrium (D’ = 0.946, r2 = 0.657), a haplotype analysis was performed. We found that the AT haplotype was associated with high BMI only when the urban environment was taken into account.ConclusionsWe have demonstrated the importance of environmental settings in studying the influence of the common UCP2 gene polymorphisms in the development of obesity in a Balinese population.

Transcriptional signatures of somatic neoblasts and germline cells in Macrostomum lignano

The regeneration-capable flatworm Macrostomum lignano is a powerful model organism to study the biology of stem cells in vivo. As a flatworm amenable to transgenesis, it complements the historically used planarian flatworm models, such as Schmidtea mediterranea. However, information on the transcriptome and markers of stem cells in M. lignano is limited. We generated a de novo transcriptome assembly and performed the first comprehensive characterization of gene expression in the proliferating cells of M. lignano, represented by somatic stem cells, called neoblasts, and germline cells. Knockdown of a selected set of neoblast genes, including Mlig-ddx39, Mlig-rrm1, Mlig-rpa3, Mlig-cdk1, and Mlig-h2a, confirmed their crucial role for the functionality of somatic neoblasts during homeostasis and regeneration. The generated M. lignano transcriptome assembly and gene expression signatures of somatic neoblasts and germline cells will be a valuable resource for future molecular studies in M. lignano. DOI: http://dx.doi.org/10.7554/eLife.20607.001

Deep characterization of a common D4Z4 variant identifies biallelic DUX4 expression as a modifier for disease penetrance in FSHD2

Facioscapulohumeral muscular dystrophy is caused by incomplete repression of the transcription factor DUX4 in skeletal muscle as a consequence of D4Z4 macrosatellite repeat contraction in chromosome 4q35 (FSHD1) or variants in genes encoding D4Z4 chromatin repressors (FSHD2). A clinical hallmark of FSHD is variability in onset and progression suggesting the presence of disease modifiers. A well-known cis modifier is the polymorphic DUX4 polyadenylation signal (PAS) that defines FSHD permissive alleles: D4Z4 chromatin relaxation on non-permissive alleles which lack the DUX4-PAS cannot cause disease in the absence of stable DUX4 mRNA. We have explored the nature and relevance of a common variant of the major FSHD haplotype 4A161, which is defined by 1.6 kb size difference of the most distal D4Z4 repeat unit. While the short variant (4A161S) has been extensively studied, we demonstrate that the long variant (4A161L) is relatively common in the European population, is capable of expressing DUX4, but that DUX4 mRNA processing differs from 4A161S. While we do not find evidence for a difference in disease severity between FSHD carriers of an 4A161S or 4A161L allele, our study does uncover biallelic DUX4 expression in FSHD2 patients. Compared to control individuals, we observed an increased frequency of FSHD2 patients homozygous for disease permissive alleles, and who are thus capable of biallelic DUX4 expression, while SMCHD1 variant carriers with only one permissive allele were significantly more often asymptomatic. This suggests that biallelic DUX4 expression lowers the threshold for disease presentation and is a modifier for disease severity in FSHD2.

Genes associated with venous thromboembolism in colorectal cancer patients

Essentials The underlying pathophysiological mechanisms behind cancer‐associated thrombosis are unknown. We compared expression profiles in tumor cells from patients with and without thrombosis. Tumors from patients with thrombosis showed significant differential gene expression profiles. Patients with thrombosis had a proinflammatory status and increased fibrin levels in the tumor.

Molecular signatures of age-associated chronic degeneration of shoulder muscles

Chronic muscle diseases are highly prevalent in the elderly causing severe mobility limitations, pain and frailty. The intrinsic molecular mechanisms are poorly understood due to multifactorial causes, slow progression with age and variations between individuals. Understanding the underlying molecular mechanisms could lead to new treatment options which are currently limited. Shoulder complaints are highly common in the elderly, and therefore, muscles of the shoulders rotator cuff could be considered as a model for chronic age-associated muscle degeneration. Diseased shoulder muscles were characterized by muscle atrophy and fatty infiltration compared with unaffected shoulder muscles. We confirmed fatty infiltration using histochemical analysis. Additionally, fibrosis and loss of contractile myosin expression were found in diseased muscles. Most cellular features, including proliferation rate, apoptosis and cell senescence, remained unchanged and genome-wide molecular signatures were predominantly similar between diseased and intact muscles. However, we found down-regulation of a small subset of muscle function genes, and up-regulation of extracellular region genes. Myogenesis was defected in muscle cell culture from diseased muscles but was restored by elevating MyoD levels. We suggest that impaired muscle functionality in a specific environment of thickened extra-cellular matrix is crucial for the development of chronic age-associated muscle degeneration.

Unraveling the polygenic architecture of complex traits using blood eQTL meta-analysis

Summary While many disease-associated variants have been identified through genome-wide association studies, their downstream molecular consequences remain unclear. To identify these effects, we performed cis- and trans-expression quantitative trait locus (eQTL) analysis in blood from 31,684 individuals through the eQTLGen Consortium. We observed that cis-eQTLs can be detected for 88% of the studied genes, but that they have a different genetic architecture compared to disease-associated variants, limiting our ability to use cis-eQTLs to pinpoint causal genes within susceptibility loci. In contrast, trans-eQTLs (detected for 37% of 10,317 studied trait-associated variants) were more informative. Multiple unlinked variants, associated to the same complex trait, often converged on trans-genes that are known to play central roles in disease etiology. We observed the same when ascertaining the effect of polygenic scores calculated for 1,263 genome-wide association study (GWAS) traits. Expression levels of 13% of the studied genes correlated with polygenic scores, and many resulting genes are known to drive these traits.

Parental haplotype-specific single-cell transcriptomics reveal incomplete epigenetic reprogramming in human female germ cells

In contrast to mouse, human female germ cells develop asynchronously. Germ cells transition to meiosis, erase genomic imprints, and reactivate the X chromosome. It is unknown if these events all appear asynchronously, and how they relate to each other. Here we combine exome sequencing of human fetal and maternal tissues with single-cell RNA-sequencing of five donors. We reconstruct full parental haplotypes and quantify changes in parental allele-specific expression, genome-wide. First we distinguish primordial germ cells (PGC), pre-meiotic, and meiotic transcriptional stages. Next we demonstrate that germ cells from various stages monoallelically express imprinted genes and confirm this by methylation patterns. Finally, we show that roughly 30% of the PGCs are still reactivating their inactive X chromosome and that this is related to transcriptional stage rather than fetal age. Altogether, we uncover the complexity and cell-to-cell heterogeneity of transcriptional and epigenetic remodeling in female human germ cells.In mammalian female germ cells, parent-specific epigenetic marks are erased and the X chromosome reactivated before entry into meiosis. Here, by combining parental haplotype reconstruction with single-cell transcriptomics of human female embryonic germ cells, the authors demonstrate that epigenetic reprogramming occurs in a heterogeneous fashion and during a broad time window up to week 14.