Search

Venous thromboembolism (VTE) is a frequent complication in patients with cancer. Homozygous carriers of the fibrinogen gamma gene (FGG) rs2066865 have a moderately increased risk of VTE, but the effect of the FGG variant in cancer is unknown. We aimed to investigate the effect of the FGG variant and active cancer on the risk of VTE. Cases with incident VTE (n=640) and a randomly selected age-weighted sub-cohort (n=3,734) were derived from a population-based cohort (the Tromsø study). Cox-regression was used to estimate hazard ratios (HR) with 95% confidence intervals (CI) for VTE according to categories of cancer and FGG In those without cancer, homozygosity at the FGG variant was associated with a 70% (HR 1.7, 95% CI: 1.2-2.3) increased risk of VTE compared to non-carriers. Cancer patients homozygous for the FGG variant had a two-fold (HR 2.0, 95% CI: 1.1-3.6) higher risk of VTE than cancer patients without the variant. Moreover, the six-months cumulative incidence of VTE among cancer patients was 6.4% (95% CI: 3.5-11.6) in homozygous carriers of FGG and 3.1% (95% CI: 2.3-4.7) in those without risk alleles. A synergistic effect was observed between rs2066865 and active cancer on the risk of VTE (synergy index: 1.81, 95% CI: 1.02-3.21, attributable proportion: 0.43, 95% CI: 0.11-0.74). In conclusion, homozygosity at the FGG variant and active cancer yielded a synergistic effect on the risk of VTE.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Structural variants (SVs) and short tandem repeats (STRs) comprise a broad group of diverse DNA variants which vastly differ in their sizes and distributions across the genome. Here, we identify genomic features of SV classes and STRs that are associated with gene expression and complex traits, including their locations relative to eGenes, likelihood of being associated with multiple eGenes, associated eGene types (e.g., coding, noncoding, level of evolutionary constraint), effect sizes, linkage disequilibrium with tagging single nucleotide variants used in GWAS, and likelihood of being associated with GWAS traits. We identify a set of high-impact SVs/STRs associated with the expression of three or more eGenes via chromatin loops and show that they are highly enriched for being associated with GWAS traits. Our study provides insights into the genomic properties of structural variant classes and short tandem repeats that are associated with gene expression and human traits.

Structural variants (SVs) and short tandem repeats (STRs) are important sources of genetic diversity but are not routinely analyzed in genetic studies because they are difficult to accurately identify and genotype. Because SVs and STRs range in size and type, it is necessary to apply multiple algorithms that incorporate different types of evidence from sequencing data and employ complex filtering strategies to discover a comprehensive set of high-quality and reproducible variants. Here we assemble a set of 719 deep whole genome sequencing (WGS) samples (mean 42×) from 477 distinct individuals which we use to discover and genotype a wide spectrum of SV and STR variants using five algorithms. We use 177 unique pairs of genetic replicates to identify factors that affect variant call reproducibility and develop a systematic filtering strategy to create of one of the most complete and well characterized maps of SVs and STRs to date.

The Genotype-Tissue Expression (GTEx) resource has provided insights into the regulatory impact of genetic variation on gene expression across human tissues; however, thus far has not considered how variation acts at the resolution of the different cell types. Here, using gene expression signatures obtained from mouse cell types, we deconvolute bulk RNA-seq samples from 28 GTEx tissues to quantify cellular composition, which reveals striking heterogeneity across these samples. Conducting eQTL analyses for GTEx liver and skin samples using cell composition estimates as interaction terms, we identify thousands of genetic associations that are cell-type-associated. The skin cell-type associated eQTLs colocalize with skin diseases, indicating that variants which influence gene expression in distinct skin cell types play important roles in traits and disease. Our study provides a framework to estimate the cellular composition of GTEx tissues enabling the functional characterization of human genetic variation that impacts gene expression in cell-type-specific manners.

Many cellular models aimed at elucidating cancer biology do not recapitulate pathobiology including tumor heterogeneity, an inherent feature of cancer that underlies treatment resistance. Here we introduce a cancer modeling paradigm using genetically engineered human pluripotent stem cells (hiPSCs) that captures authentic cancer pathobiology. Orthotopic engraftment of the neural progenitor cells derived from hiPSCs that have been genome-edited to contain tumor-associated genetic driver mutations revealed by The Cancer Genome Atlas project for glioblastoma (GBM) results in formation of high-grade gliomas. Similar to patient-derived GBM, these models harbor inter-tumor heterogeneity resembling different GBM molecular subtypes, intra-tumor heterogeneity, and extrachromosomal DNA amplification. Re-engraftment of these primary tumor neurospheres generates secondary tumors with features characteristic of patient samples and present mutation-dependent patterns of tumor evolution. These cancer avatar models provide a platform for comprehensive longitudinal assessment of human tumor development as governed by molecular subtype mutations and lineage-restricted differentiation.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Induced pluripotent stem cell derived cardiovascular progenitor cells (iPSC-CVPCs) provide an unprecedented platform for examining the molecular underpinnings of cardiac development and disease etiology, but also have great potential to play pivotal roles in the future of regenerative medicine and pharmacogenomic studies. Biobanks like iPSCORE ( Stacey et al., 2013 ; Panopoulos et al., 2017 ), which contain iPSCs generated from hundreds of genetically and ethnically diverse individuals, are an invaluable resource for conducting these studies. Here, we present an optimized, cost-effective and highly standardized protocol for large-scale derivation of human iPSC-CVPCs using small molecules and purification using metabolic selection. We have successfully applied this protocol to derive iPSC-CVPCs from 154 different iPSCORE iPSC lines obtaining large quantities of highly pure cardiac cells. An important component of our protocol is Cell confluency estimates (ccEstimate), an automated methodology for estimating the time when an iPSC monolayer will reach 80% confluency, which is optimal for initiating iPSC-CVPC derivation, and enables the protocol to be readily used across iPSC lines with different growth rates. Moreover, we showed that cellular heterogeneity across iPSC-CVPCs is due to varying proportions of two distinct cardiac cell types: cardiomyocytes (CMs) and epicardium-derived cells (EPDCs), both of which have been shown to have a critical function in heart regeneration. This protocol eliminates the need of iPSC line-to-line optimization and can be easily adapted and scaled to high-throughput studies or to generate large quantities of cells suitable for regenerative medicine applications.

A substantial fraction of SNPs associated with human traits and diseases through genome-wide association studies (GWAS) are likely regulatory variants as they tend to be located within enhancers and associated with differential gene expression. Thus, as a key step in implementing personalized medicine, it is important to identify regulatory variants in the human genome, and characterize their underlying molecular mechanisms. However, identifying and elucidating the functions of regulatory variants is currently challenging as these variants show similar associations with many other neutral variants due to linkage disequilibrium, can be quite far from the gene(s) they regulate, and often have cell type-specific effects. In order to overcome these challenges and interrogate the function of these regulatory variants, it could be possible to examine and integrate epigenetic information in a cell type dependent manner. Here, I present three studies which focus on the functionality of the epigenome – specifically chromatin looping and co-accessibility – in the context of gene regulation, genetics, and disease. I present a tool for computationally working with chromatin loop data, and utilize this tool to show that genetic variation is not associated with large changes in chromatin looping, but rather small modulation in contact propensity which are associated with large changes in gene expression. I then examine chromatin co-accessibility, and show that genetic variants may be able to mediate long range effects on genes and accessible sites hundreds of megabases away – across entire chromosomes – which are associated with cell type relevant genes and diseases.

While genetic variation at chromatin loops is relevant for human disease, the relationships between contact propensity (the probability that loci at loops physically interact), genetics, and gene regulation are unclear. We quantitatively interrogate these relationships by comparing Hi-C and molecular phenotype data across cell types and haplotypes. While chromatin loops consistently form across different cell types, they have subtle quantitative differences in contact frequency that are associated with larger changes in gene expression and H3K27ac. For the vast majority of loci with quantitative differences in contact frequency across haplotypes, the changes in magnitude are smaller than those across cell types; however, the proportional relationships between contact propensity, gene expression, and H3K27ac are consistent. These findings suggest that subtle changes in contact propensity have a biologically meaningful role in gene regulation and could be a mechanism by which regulatory genetic variants in loop anchors mediate effects on expression.

The MHC region is highly associated with autoimmune and infectious diseases. Here we conduct an in-depth interrogation of associations between genetic variation, gene expression and disease. We create a comprehensive map of regulatory variation in the MHC region using WGS from 419 individuals to call eight-digit HLA types and RNA-seq data from matched iPSCs. Building on this regulatory map, we explored GWAS signals for 4083 traits, detecting colocalization for 180 disease loci with eQTLs. We show that eQTL analyses taking HLA type haplotypes into account have substantially greater power compared with only using single variants. We examined the association between the 8.1 ancestral haplotype and delayed colonization in Cystic Fibrosis, postulating that downregulation of RNF5 expression is the likely causal mechanism. Our study provides insights into the genetic architecture of the MHC region and pinpoints disease associations that are due to differential expression of HLA genes and non-HLA genes.

Despite the importance of understanding how variability across induced pluripotent stem cell (iPSC) lines due to non-genetic factors (clone and passage) influences their differentiation outcome, large-scale studies capable of addressing this question have not yet been conducted. Here, we differentiated 191 iPSC lines to generate iPSC-derived cardiovascular progenitor cells (iPSC-CVPCs). We observed cellular heterogeneity across the iPSC-CVPC samples due to varying fractions of two cell types: cardiomyocytes (CMs) and epicardium-derived cells (EPDCs). Comparing the transcriptomes of CM-fated and EPDC-fated iPSCs, we discovered that 91 signature genes and X chromosome dosage differences are associated with these two distinct cardiac developmental trajectories. In an independent set of 39 iPSCs differentiated into CMs, we confirmed that sex and transcriptional differences affect cardiac-fate outcome. Our study provides novel insights into how iPSC transcriptional and X chromosome gene dosage differences influence their response to differentiation stimuli and, hence, cardiac cell fate.

Genetic variants affecting pancreatic islet enhancers are central to T2D risk, but the gene targets of islet enhancer activity are largely unknown. We generate a high-resolution map of islet chromatin loops using Hi-C assays in three islet samples and use loops to annotate target genes of islet enhancers defined using ATAC-seq and published ChIP-seq data. We identify candidate target genes for thousands of islet enhancers, and find that enhancer looping is correlated with islet-specific gene expression. We fine-map T2D risk variants affecting islet enhancers, and find that candidate target genes of these variants defined using chromatin looping and eQTL mapping are enriched in protein transport and secretion pathways. At IGF2BP2, a fine-mapped T2D variant reduces islet enhancer activity and IGF2BP2 expression, and conditional inactivation of IGF2BP2 in mouse islets impairs glucose-stimulated insulin secretion. Our findings provide a resource for studying islet enhancer function and identifying genes involved in T2D risk.

The cardiac transcription factor (TF) gene NKX2-5 has been associated with electrocardiographic (EKG) traits through genome-wide association studies (GWASs), but the extent to which differential binding of NKX2-5 at common regulatory variants contributes to these traits has not yet been studied. We analyzed transcriptomic and epigenomic data from induced pluripotent stem cell-derived cardiomyocytes from seven related individuals, and identified ~2,000 single-nucleotide variants associated with allele-specific effects (ASE-SNVs) on NKX2-5 binding. NKX2-5 ASE-SNVs were enriched for altered TF motifs, for heart-specific expression quantitative trait loci and for EKG GWAS signals. Using fine-mapping combined with epigenomic data from induced pluripotent stem cell-derived cardiomyocytes, we prioritized candidate causal variants for EKG traits, many of which were NKX2-5 ASE-SNVs. Experimentally characterizing two NKX2-5 ASE-SNVs (rs3807989 and rs590041) showed that they modulate the expression of target genes via differential protein binding in cardiac cells, indicating that they are functional variants underlying EKG GWAS signals. Our results show that differential NKX2-5 binding at numerous regulatory variants across the genome contributes to EKG phenotypes.

Pancreatic endocrine cell differentiation is orchestrated by the action of transcription factors that operate in a gene regulatory network to activate endocrine lineage genes and repress lineage-inappropriate genes. MicroRNAs (miRNAs) are important modulators of gene expression, yet their role in endocrine cell differentiation has not been systematically explored. Here we characterize miRNA-regulatory networks active in human endocrine cell differentiation by combining small RNA sequencing, miRNA over-expression, and network modeling approaches. Our analysis identified Let-7g, Let-7a, miR-200a, miR-127, and miR-375 as endocrine-enriched miRNAs that drive endocrine cell differentiation-associated gene expression changes. These miRNAs are predicted to target different transcription factors, which converge on genes involved in cell cycle regulation. When expressed in human embryonic stem cell-derived pancreatic progenitors, these miRNAs induce cell cycle exit and promote endocrine cell differentiation. Our study delineates the role of miRNAs in human endocrine cell differentiation and identifies miRNAs that could facilitate endocrine cell reprogramming.