IFN-γ activates an immune-like regulatory network in the cardiac vascular endothelium

The regulatory mechanisms underlying the response to pro-inflammatory cytokines during myocarditis are poorly understood. Here, we use iPSC-derived cardiovascular progenitor cells (CVPCs) to model the response to interferon gamma (IFN-γ) during myocarditis. We generate RNA-seq and ATAC-seq for four CVPCs that were treated with IFN-γ and compare them with paired untreated controls. Transcriptional differences after treatment show that IFN-γ initiates an innate immune cell-like response in the vascular cardiac endothelium. IFN-γ treatment also shifts the CVPC transcriptome towards the adult coronary artery and aorta profiles and expands the relative endothelial cell population in all four CVPC lines. Analysis of the accessible chromatin shows that IFN-γ is a potent chromatin remodeler and establishes an IRF-STAT immune-cell like regulatory network. Our findings reveal insights into the endothelial-specific protective mechanisms during myocarditis.

Multi-omic QTL mapping in early developmental tissues reveals phenotypic and temporal complexity of regulatory variants underlying GWAS loci

Most GWAS loci are presumed to affect gene regulation, however, only ∼43% colocalize with expression quantitative trait loci (eQTLs). To address this colocalization gap, we identify eQTLs, chromatin accessibility QTLs (caQTLs), and histone acetylation QTLs (haQTLs) using molecular samples from three early developmental (EDev) tissues. Through colocalization, we annotate 586 GWAS loci for 17 traits by QTL complexity, QTL phenotype, and QTL temporal specificity. We show that GWAS loci are highly enriched for colocalization with complex QTL modules that affect multiple elements (genes and/or peaks). We also demonstrate that caQTLs and haQTLs capture regulatory variations not associated with eQTLs and explain ∼49% of the functionally annotated GWAS loci. Additionally, we show that EDev-unique QTLs are strongly depleted for colocalizing with GWAS loci. By conducting one of the largest multi-omic QTL studies to date, we demonstrate that many GWAS loci exhibit phenotypic complexity and therefore, are missed by traditional eQTL analyses.

Complex regulatory networks influence pluripotent cell state transitions in human iPSCs

Stem cells exist in vitro in a spectrum of interconvertible pluripotent states. Analyzing hundreds of hiPSCs derived from different individuals, we show the proportions of these pluripotent states vary considerably across lines. We discover 13 gene network modules (GNMs) and 13 regulatory network modules (RNMs), which are highly correlated with each other suggesting that the coordinated co-accessibility of regulatory elements in the RNMs likely underlie the coordinated expression of genes in the GNMs. Epigenetic analyses reveal that regulatory networks underlying self-renewal and pluripotency are more complex than previously realized. Genetic analyses identify thousands of regulatory variants that overlapped predicted transcription factor binding sites and are associated with chromatin accessibility in the hiPSCs. We show that the master regulator of pluripotency, the NANOG-OCT4 Complex, and its associated network are significantly enriched for regulatory variants with large effects, suggesting that they play a role in the varying cellular proportions of pluripotency states between hiPSCs. Our work bins tens of thousands of regulatory elements in hiPSCs into discrete regulatory networks, shows that pluripotency and self-renewal processes have a surprising level of regulatory complexity, and suggests that genetic factors may contribute to cell state transitions in human iPSC lines.

Single-nucleotide variants within heart enhancers increase binding affinity and disrupt heart development

Transcriptional enhancers direct precise gene expression patterns during development and harbor the majority of variants associated with phenotypic diversity, evolutionary adaptations, and disease. Pinpointing which enhancer variants contribute to changes in gene expression and phenotypes is a major challenge. Here, we find that suboptimal or low-affinity binding sites are necessary for precise gene expression during heart development. Single-nucleotide variants (SNVs) can optimize the affinity of ETS binding sites, causing gain-of-function (GOF) gene expression, cell migration defects, and phenotypes as severe as extra beating hearts in the marine chordate Ciona robusta. In human induced pluripotent stem cell (iPSC)-derived cardiomyocytes, a SNV within a human GATA4 enhancer increases ETS binding affinity and causes GOF enhancer activity. The prevalence of suboptimal-affinity sites within enhancers creates a vulnerability whereby affinity-optimizing SNVs can lead to GOF gene expression, changes in cellular identity, and organismal-level phenotypes that could contribute to the evolution of novel traits or diseases.

eQTL mapping in fetal-like pancreatic progenitor cells reveals early developmental insights into diabetes risk

The impact of genetic regulatory variation active in early pancreatic development on adult pancreatic disease and traits is not well understood. Here, we generate a panel of 107 fetal-like iPSC-derived pancreatic progenitor cells (iPSC-PPCs) from whole genome-sequenced individuals and identify 4065 genes and 4016 isoforms whose expression and/or alternative splicing are affected by regulatory variation. We integrate eQTLs identified in adult islets and whole pancreas samples, which reveal 1805 eQTL associations that are unique to the fetal-like iPSC-PPCs and 1043 eQTLs that exhibit regulatory plasticity across the fetal-like and adult pancreas tissues. Colocalization with GWAS risk loci for pancreatic diseases and traits show that some putative causal regulatory variants are active only in the fetal-like iPSC-PPCs and likely influence disease by modulating expression of disease-associated genes in early development, while others with regulatory plasticity likely exert their effects in both the fetal and adult pancreas by modulating expression of different disease genes in the two developmental stages.

Analysis of regulatory network modules in hundreds of human stem cell lines reveals complex epigenetic and genetic factors contribute to pluripotency state differences between subpopulations

Stem cells exist in vitro in a spectrum of interconvertible pluripotent states. Analyzing hundreds of hiPSCs derived from different individuals, we show the proportions of these pluripotent states vary considerably across lines. We discovered 13 gene network modules (GNMs) and 13 regulatory network modules (RNMs), which were highly correlated with each other suggesting that the coordinated co-accessibility of regulatory elements in the RNMs likely underlied the coordinated expression of genes in the GNMs. Epigenetic analyses revealed that regulatory networks underlying self-renewal and pluripotency have a surprising level of complexity. Genetic analyses identified thousands of regulatory variants that overlapped predicted transcription factor binding sites and were associated with chromatin accessibility in the hiPSCs. We show that the master regulator of pluripotency, the NANOG-OCT4 Complex, and its associated network were significantly enriched for regulatory variants with large effects, suggesting that they may play a role in the varying cellular proportions of pluripotency states between hiPSCs. Our work captures the coordinated activity of tens of thousands of regulatory elements in hiPSCs and bins these elements into discrete functionally characterized regulatory networks, shows that regulatory elements in pluripotency networks harbor variants with large effects, and provides a rich resource for future pluripotent stem cell research.

The human pangenome reference anticipates equitable and fundamental genomic insights

For the past few years, researchers in the Human Pangenome Reference Consortium (HPRC) have been working to catalog almost all human genomic diversity. Frazer and Schork preview an article recently published in Nature, "A draft human pangenome reference,"¹ which represents the initial release of 47 fully phased diploid assemblies of genomes of individuals with diverse ancestries.

Fine mapping spatiotemporal mechanisms of genetic variants underlying cardiac traits and disease

The causal variants and genes underlying thousands of cardiac GWAS signals have yet to be identified. Here, we leverage spatiotemporal information on 966 RNA-seq cardiac samples and perform an expression quantitative trait locus (eQTL) analysis detecting eQTLs considering both eGenes and eIsoforms. We identify 2,578 eQTLs associated with a specific developmental stage-, tissue- and/or cell type. Colocalization between eQTL and GWAS signals of five cardiac traits identified variants with high posterior probabilities for being causal in 210 GWAS loci. Pulse pressure GWAS loci are enriched for colocalization with fetal- and smooth muscle- eQTLs; pulse rate with adult- and cardiac muscle- eQTLs; and atrial fibrillation with cardiac muscle- eQTLs. Fine mapping identifies 79 credible sets with five or fewer SNPs, of which 15 were associated with spatiotemporal eQTLs. Our study shows that many cardiac GWAS variants impact traits and disease in a developmental stage-, tissue- and/or cell type-specific fashion.

Type 1 diabetes risk genes mediate pancreatic beta cell survival in response to proinflammatory cytokines

We combined functional genomics and human genetics to investigate processes that affect type 1 diabetes (T1D) risk by mediating beta cell survival in response to proinflammatory cytokines. We mapped 38,931 cytokine-responsive candidate cis-regulatory elements (cCREs) in beta cells using ATAC-seq and snATAC-seq and linked them to target genes using co-accessibility and HiChIP. Using a genome-wide CRISPR screen in EndoC-βH1 cells, we identified 867 genes affecting cytokine-induced survival, and genes promoting survival and up-regulated in cytokines were enriched at T1D risk loci. Using SNP-SELEX, we identified 2,229 variants in cytokine-responsive cCREs altering transcription factor (TF) binding, and variants altering binding of TFs regulating stress, inflammation, and apoptosis were enriched for T1D risk. At the 16p13 locus, a fine-mapped T1D variant altering TF binding in a cytokine-induced cCRE interacted with SOCS1, which promoted survival in cytokine exposure. Our findings reveal processes and genes acting in beta cells during inflammation that modulate T1D risk.

Systems genomics in age-related macular degeneration

Genomic studies in age-related macular degeneration (AMD) have identified genetic variants that account for the majority of AMD risk. An important next step is to understand the functional consequences and downstream effects of the identified AMD-associated genetic variants. Instrumental for this next step are 'omics' technologies, which enable high-throughput characterization and quantification of biological molecules, and subsequent integration of genomics with these omics datasets, a field referred to as systems genomics. Single cell sequencing studies of the retina and choroid demonstrated that the majority of candidate AMD genes identified through genomic studies are expressed in non-neuronal cells, such as the retinal pigment epithelium (RPE), glia, myeloid and choroidal cells, highlighting that many different retinal and choroidal cell types contribute to the pathogenesis of AMD. Expression quantitative trait locus (eQTL) studies in retinal tissue have identified putative causal genes by demonstrating a genetic overlap between gene regulation and AMD risk. Linking genetic data to complement measurements in the systemic circulation has aided in understanding the effect of AMD-associated genetic variants in the complement system, and supports that protein QTL (pQTL) studies in plasma or serum samples may aid in understanding the effect of genetic variants and pinpointing causal genes in AMD. A recent epigenomic study fine-mapped AMD causal variants by determing regulatory regions in RPE cells differentiated from induced pluripotent stem cells (iPSC-RPE). Another approach that is being employed to pinpoint causal AMD genes is to produce synthetic DNA assemblons representing risk and protective haplotypes, which are then delivered to cellular or animal model systems. Pinpointing causal genes and understanding disease mechanisms is crucial for the next step towards clinical translation. Clinical trials targeting proteins encoded by the AMD-associated genomic loci C3, CFB, CFI, CFH, and ARMS2/HTRA1 are currently ongoing, and a phase III clinical trial for C3 inhibition recently showed a modest reduction of lesion growth in geographic atrophy. The EYERISK consortium recently developed a genetic test for AMD that allows genotyping of common and rare variants in AMD-associated genes. Polygenic risk scores (PRS) were applied to quantify AMD genetic risk, and may aid in predicting AMD progression. In conclusion, genomic studies represent a turning point in our exploration of AMD. The results of those studies now serve as a driving force for several clinical trials. Expanding to omics and systems genomics will further decipher function and causality from the associations that have been reported, and will enable the development of therapies that will lessen the burden of AMD.

Ultra-Sharp Nanowire Arrays Natively Permeate, Record, and Stimulate Intracellular Activity in Neuronal and Cardiac Networks

We report innovative scalable, vertical, ultra-sharp nanowire arrays that are individually addressable to enable long-term, native recordings of intracellular potentials. Stable amplitudes of intracellular potentials from 3D tissue-like networks of neurons and cardiomyocytes are obtained. Individual electrical addressability is necessary for high-fidelity intracellular electrophysiological recordings. This study paves the way toward predictive, high-throughput, and low-cost electrophysiological drug screening platforms.

SARS-CoV-2 susceptibility and COVID-19 disease severity are associated with genetic variants affecting gene expression in a variety of tissues

Variability in SARS-CoV-2 susceptibility and COVID-19 disease severity between individuals is partly due to genetic factors. Here, we identify 4 genomic loci with suggestive associations for SARS-CoV-2 susceptibility and 19 for COVID-19 disease severity. Four of these 23 loci likely have an ethnicity-specific component. Genome-wide association study (GWAS) signals in 11 loci colocalize with expression quantitative trait loci (eQTLs) associated with the expression of 20 genes in 62 tissues/cell types (range: 1:43 tissues/gene), including lung, brain, heart, muscle, and skin as well as the digestive system and immune system. We perform genetic fine mapping to compute 99% credible SNP sets, which identify 10 GWAS loci that have eight or fewer SNPs in the credible set, including three loci with one single likely causal SNP. Our study suggests that the diverse symptoms and disease severity of COVID-19 observed between individuals is associated with variants across the genome, affecting gene expression levels in a wide variety of tissue types.

Deciphering the genetic architecture and ethnographic distribution of IRD in three ethnic populations by whole genome sequence analysis

Patients with inherited retinal dystrophies (IRDs) were recruited from two understudied populations: Mexico and Pakistan as well as a third well-studied population of European Americans to define the genetic architecture of IRD by performing whole-genome sequencing (WGS). Whole-genome analysis was performed on 409 individuals from 108 unrelated pedigrees with IRDs. All patients underwent an ophthalmic evaluation to establish the retinal phenotype. Although the 108 pedigrees in this study had previously been examined for mutations in known IRD genes using a wide range of methodologies including targeted gene(s) or mutation(s) screening, linkage analysis and exome sequencing, the gene mutations responsible for IRD in these 108 pedigrees were not determined. WGS was performed on these pedigrees using Illumina X10 at a minimum of 30X depth. The sequence reads were mapped against hg19 followed by variant calling using GATK. The genome variants were annotated using SnpEff, PolyPhen2, and CADD score; the structural variants (SVs) were called using GenomeSTRiP and LUMPY. We identified potential causative sequence alterations in 61 pedigrees (57%), including 39 novel and 54 reported variants in IRD genes. For 57 of these pedigrees the observed genotype was consistent with the initial clinical diagnosis, the remaining 4 had the clinical diagnosis reclassified based on our findings. In seven pedigrees (12%) we observed atypical causal variants, i.e. unexpected genotype(s), including 4 pedigrees with causal variants in more than one IRD gene within all affected family members, one pedigree with intrafamilial genetic heterogeneity (different affected family members carrying causal variants in different IRD genes), one pedigree carrying a dominant causative variant present in pseudo-recessive form due to consanguinity and one pedigree with a de-novo variant in the affected family member. Combined atypical and large structural variants contributed to about 20% of cases. Among the novel mutations, 75% were detected in Mexican and 50% found in European American pedigrees and have not been reported in any other population while only 20% were detected in Pakistani pedigrees and were not previously reported. The remaining novel IRD causative variants were listed in gnomAD but were found to be very rare and population specific. Mutations in known IRD associated genes contributed to pathology in 63% Mexican, 60% Pakistani and 45% European American pedigrees analyzed. Overall, contribution of known IRD gene variants to disease pathology in these three populations was similar to that observed in other populations worldwide. This study revealed a spectrum of mutations contributing to IRD in three populations, identified a large proportion of novel potentially causative variants that are specific to the corresponding population or not reported in gnomAD and shed light on the genetic architecture of IRD in these diverse global populations.

The study was performed to identify the underlying cause of inherited retinal degeneration (IRD) in 409 individuals from 108 families. Primarily, these families were recruited from three different geographic regions: Mexico, Pakistan and European Americans from the United States. Blood samples were collected from all individuals for genome analysis. This analysis detected causative variants in 61 out of the 108 pedigrees. A total of 93 gene variants were found in the 61 families. Among these, 54 were previously reported as causative variants and the remaining 39 have not been reported in IRD pedigrees. Interestingly, 54% of these novel variants were not listed in gnomAD. In addition to these findings, complex causative genotypes were observed in 20% of pedigrees. Overall, causative variants were detected in 63% Mexican, 60% Pakistani and 45% European American pedigrees. This study revealed the distribution of IRD causative variants in pedigrees with diverse ethnic and geographic backgrounds.

Genomics Links Inflammation With Neurocognitive Impairment in Children Living With Human Immunodeficiency Virus Type-1

BACKGROUND

We identified host single-nucleotide variants (SNVs) associated with neurocognitive impairment (NCI) in perinatally HIV-infected (PHIV) children.

METHODS

Whole-exome sequencing (WES) was performed on 217 PHIV with cognitive score for age (CSA) < 70 and 247 CSA ≥ 70 (discovery cohort [DC]). SNVs identified in DC were evaluated in 2 validation cohorts (VC). Logistic regression was used to estimate adjusted odds ratios (ORs) for NCI. A human microglia NLRP3 inflammasome assay characterized the role of identified genes.

RESULTS

Twenty-nine SNVs in 24 genes reaching P ≤ .002 and OR ≥ 1.5 comparing CSA < 70 to CSA ≥ 70 were identified in the DC, of which 3 SNVs were identified in VCs for further study. Combining the 3 cohorts, SNV in CCRL2 (rs3204849) was associated with decreased odds of NCI (P < .0001); RETREG1/FAM134B (rs61733811) and YWHAH (rs73884247) were associated with increased risk of NCI (P < .0001 and P < .001, respectively). Knockdown of CCRL2 led to decreased microglial release of IL-1β following exposure to ssRNA40 while knockdown of RETREG1 and YWHAH resulted in increased IL-1β release.

CONCLUSIONS

Using WES and 2 VCs, and gene silencing of microglia we identified 3 genetic variants associated with NCI and inflammation in HIV-infected children.

Enhancer release and retargeting activates disease-susceptibility genes

The functional engagement between an enhancer and its target promoter ensures precise gene transcription1. Understanding the basis of promoter choice by enhancers has important implications for health and disease. Here we report that functional loss of a preferred promoter can release its partner enhancer to loop to and activate an alternative promoter (or alternative promoters) in the neighbourhood. We refer to this target-switching process as 'enhancer release and retargeting'. Genetic deletion, motif perturbation or mutation, and dCas9-mediated CTCF tethering reveal that promoter choice by an enhancer can be determined by the binding of CTCF at promoters, in a cohesin-dependent manner-consistent with a model of 'enhancer scanning' inside the contact domain. Promoter-associated CTCF shows a lower affinity than that at chromatin domain boundaries and often lacks a preferred motif orientation or a partnering CTCF at the cognate enhancer, suggesting properties distinct from boundary CTCF. Analyses of cancer mutations, data from the GTEx project and risk loci from genome-wide association studies, together with a focused CRISPR interference screen, reveal that enhancer release and retargeting represents an overlooked mechanism that underlies the activation of disease-susceptibility genes, as exemplified by a risk locus for Parkinson's disease (NUCKS1-RAB7L1) and three loci associated with cancer (CLPTM1L-TERT, ZCCHC7-PAX5 and PVT1-MYC).

Insights into genetic factors contributing to variability in SARS-CoV-2 susceptibility and COVID-19 disease severity

Variability in SARS-CoV-2 susceptibility and COVID-19 disease severity between individuals is partly due to genetic factors. Here, we applied colocalization to compare summary statistics for 16 GWASs from the COVID-19 Host Genetics Initiative to investigate similarities and differences in their genetic signals. We identified 9 loci associated with susceptibility (one with two independent GWAS signals; one with an ethnicity-specific signal), 14 associated with severity (one with two independent GWAS signals; two with ethnicity-specific signals) and one harboring two discrepant GWAS signals (one for susceptibility; one for severity). Utilizing colocalization we also identified 45 GTEx tissues that had eQTL(s) for 18 genes strongly associated with GWAS signals in eleven loci (1-4 genes per locus). Some of these genes showed tissue-specific altered expression and others showed altered expression in up to 41 different tissue types. Our study provides insights into the complex molecular mechanisms underlying inherited predispositions to COVID-19-disease phenotypes.

Identification of rare and common regulatory variants in pluripotent cells using population-scale transcriptomics

Induced pluripotent stem cells (iPSCs) are an established cellular system to study the impact of genetic variants in derived cell types and developmental contexts. However, in their pluripotent state, the disease impact of genetic variants is less known. Here, we integrate data from 1,367 human iPSC lines to comprehensively map common and rare regulatory variants in human pluripotent cells. Using this population-scale resource, we report hundreds of novel colocalization events for human traits specific to iPSCs, and find increased power to identify rare regulatory variants compared with somatic tissues. Finally, we demonstrate how iPSCs enable the identification of causal genes for rare diseases.

Systematic analysis of binding of transcription factors to noncoding variants

Many sequence variants have been linked to complex human traits and diseases1, but deciphering their biological functions remains challenging, as most of them reside in noncoding DNA. Here we have systematically assessed the binding of 270 human transcription factors to 95,886 noncoding variants in the human genome using an ultra-high-throughput multiplex protein-DNA binding assay, termed single-nucleotide polymorphism evaluation by systematic evolution of ligands by exponential enrichment (SNP-SELEX). The resulting 828 million measurements of transcription factor-DNA interactions enable estimation of the relative affinity of these transcription factors to each variant in vitro and evaluation of the current methods to predict the effects of noncoding variants on transcription factor binding. We show that the position weight matrices of most transcription factors lack sufficient predictive power, whereas the support vector machine combined with the gapped k-mer representation show much improved performance, when assessed on results from independent SNP-SELEX experiments involving a new set of 61,020 sequence variants. We report highly predictive models for 94 human transcription factors and demonstrate their utility in genome-wide association studies and understanding of the molecular pathways involved in diverse human traits and diseases.

Revealing Instability: Genetic Variation Underlies Variability in mESC Pluripotency

mESCs can self-renew indefinitely in vitro. However, depending on culture conditions some strains are more unstable than others. In this issue of Cell Stem Cell, Skelly et al. (2020) and Ortmann et al. (2020) shed light into the role genetic variation plays in control of ground state pluripotency.

Detection and validation of novel mutations in MERTK in a simplex case of retinal degeneration using WGS and hiPSC-RPEs model

Inherited retinal degenerations (IRDs) are a group of genetically heterogeneous conditions with a broad phenotypic heterogeneity. Here, we report detection and validation of the underlying cause of progressive retinal degeneration in a nuclear family of European descent with a single affected individual. Whole genome sequencing of the proband and her unaffected sibling identified a novel intron 8 donor splice site variant (c.1296 + 1G>A) and a novel 731 base pair deletion encompassing exon 9 (Chr2:g.112751488_112752218 del) resulting in c.1297_1451del; p.K433_G484fsTer3 in the Mer tyrosine kinase protooncogene (MERTK), which is highly expressed in the retinal pigment epithelium (RPE). The proband carried both variants in the heterozygous state, which segregated with disease in the pedigree. These MERTK variants are predicted to result in the defective splicing of exon 8 and loss of exon 9 respectively. To evaluate the impact of these novel variants, peripheral blood mononuclear cells of the proband and her parents were reprogrammed to humaninduced pluripotent stem cell (hiPSC) lines, which were subsequently differentiated to hiPSC-RPE. Analysis of the proband's hiPSC-RPE revealed the absence of both MERTK transcript and its respective protein as well as abnormal phagocytosis when compared with the parental hiPSC-RPE. In summary, whole genome sequencing identified novel compound heterozygous variants in MERTK as the underlying cause of progressive retinal degeneration in a simplex case. Further, analysis using an hiPSC-RPE model established the functional impact of novel MERTK mutations and revealed the potential mechanism underlying pathology in the proband.

In vitro Differentiation of Human iPSC-derived Cardiovascular Progenitor Cells (iPSC-CVPCs)

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.

Discovery and quality analysis of a comprehensive set of structural variants and short tandem repeats

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.

Properties of structural variants and short tandem repeats associated with gene expression and complex traits

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.

Genomic integrity of human induced pluripotent stem cells across nine studies in the NHLBI NextGen program

Human induced pluripotent stem cell (hiPSC) lines have previously been generated through the NHLBI sponsored NextGen program at nine individual study sites. Here, we examined the structural integrity of 506 hiPSC lines as determined by copy number variations (CNVs). We observed that 149 hiPSC lines acquired 258 CNVs relative to donor DNA. We identified six recurrent regions of CNVs on chromosomes 1, 2, 3, 16 and 20 that overlapped with cancer associated genes. Furthermore, the genes mapping to regions of acquired CNVs show an enrichment in cancer related biological processes (IL6 production) and signaling cascades (JNK cascade & NFκB cascade). The genomic region of instability on chr20 (chr20q11.2) includes transcriptomic signatures for cancer associated genes such as ID1, BCL2L1, TPX2, PDRG1 and HCK. Of these HCK shows statistically significant differential expression between carrier and non-carrier hiPSC lines. Overall, while a low level of genomic instability was observed in the NextGen generated hiPSC lines, the observation of structural instability in regions with known cancer associated genes substantiates the importance of systematic evaluation of genetic variations in hiPSCs before using them as disease/research models.

Cellular deconvolution of GTEx tissues powers discovery of disease and cell-type associated regulatory variants

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.

Longitudinal assessment of tumor development using cancer avatars derived from genetically engineered pluripotent stem cells

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.

In vitro Differentiation of Human iPSC-derived Retinal Pigment Epithelium Cells (iPSC-RPE)

Induced Pluripotent Stem Cells (iPSCs) serve as an excellent model system for studying the molecular underpinnings of tissue development. Human iPSC-derived retinal pigment epithelium (iPSC-RPE) cells have fetal-like molecular profiles. Hence, biobanks like iPSCORE, which contain iPSCs generated from hundreds of individuals, are an invaluable resource for examining how common genetic variants exert their effects during RPE development resulting in individuals having different propensities to develop Age-related Macular Degeneration (AMD) as adults. Here, we present an optimized, cost-effective and highly reproducible protocol for derivation of human iPSC-RPE cells using small molecules under serum-free condition and for their quality control using flow cytometry and immunofluorescence. While most previous protocols have required laborious manual selection to enrich for iPSC-RPE cells, our protocol uses whole culture passaging and yields a large number of iPSC-RPE cells with high purity (88-98.1% ZO-1 and MiTF double positive cells). The simplicity and robustness of this protocol would enable its adaption for high-throughput applications involving the generation of iPSC-RPE samples from hundreds of individuals.

Systematic genetic analysis of the MHC region reveals mechanistic underpinnings of HLA type associations with disease

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.

A Network of microRNAs Acts to Promote Cell Cycle Exit and Differentiation of Human Pancreatic Endocrine Cells

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.

Association of Human iPSC Gene Signatures and X Chromosome Dosage with Two Distinct Cardiac Differentiation Trajectories

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.

•

Cellular heterogeneity across iPSC-CVPCs due to varying fractions of CMs and EPDCs

•

iPSC non-genetic factors (clone and passage) associated with cardiac cell fate

•

Expression levels of signature genes in iPSCs associated with cardiac lineage fate

•

iPSC donor sex plays a role in cardiac lineage fate

Identification of Common and Rare Genetic Variation Associated With Plasma Protein Levels Using Whole-Exome Sequencing and Mass Spectrometry

BACKGROUND

Identifying genetic variation associated with plasma protein levels, and the mechanisms by which they act, could provide insight into alterable processes involved in regulation of protein levels. Although protein levels can be affected by genetic variants, their estimation can also be biased by missense variants in coding exons causing technical artifacts. Integrating genome sequence genotype data with mass spectrometry-based protein level estimation could reduce bias, thereby improving detection of variation that affects RNA or protein metabolism.

METHODS

Here, we integrate the blood plasma protein levels of 664 proteins from 165 participants of the Tromsø Study, measured via tandem mass tag mass spectrometry, with whole-exome sequencing data to identify common and rare genetic variation associated with peptide and protein levels (protein quantitative trait loci [pQTLs]). We additionally use literature and database searches to prioritize putative functional variants for each pQTL.

RESULTS

We identify 109 independent associations (36 protein and 73 peptide) and use genotype data to exclude 49 (4 protein and 45 peptide) as technical artifacts. We describe 2 particular cases of rare variation: 1 associated with the complement pathway and 1 with platelet degranulation. We identify putative functional variants and show that pQTLs act through diverse molecular mechanisms that affect both RNA and protein metabolism.

CONCLUSIONS

We show that although the majority of pQTLs exert their effects by modulating RNA metabolism, many affect protein levels directly. Our work demonstrates the extent by which pQTL studies are affected by technical artifacts and highlights how prioritizing the functional variant in pQTL studies can lead to insights into the molecular steps by which a protein may be regulated.

Fibrinogen gamma gene rs2066865 and risk of cancer-related venous thromboembolism

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.

Allele-specific NKX2-5 binding underlies multiple genetic associations with human electrocardiographic traits

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.

Activation of hedgehog signaling associates with early disease progression in chronic lymphocytic leukemia

Targeted sequencing of 103 leukemia-associated genes in leukemia cells from 841 treatment-naive patients with chronic lymphocytic leukemia (CLL) identified 89 (11%) patients as having CLL cells with mutations in genes encoding proteins that putatively are involved in hedgehog (Hh) signaling. Consistent with this finding, there was a significant association between the presence of these mutations and the expression of GLI1 (χ2 test, P < .0001), reflecting activation of the Hh pathway. However, we discovered that 38% of cases without identified mutations also were GLI1+ Patients with GLI1+ CLL cells had a shorter median treatment-free survival than patients with CLL cells lacking expression of GLI1 independent of IGHV mutation status. We found that GANT61, a small molecule that can inhibit GLI1, was highly cytotoxic for GLI1+ CLL cells relative to that of CLL cells without GLI1. Collectively, this study shows that a large proportion of patients have CLL cells with activated Hh signaling, which is associated with early disease progression and enhanced sensitivity to inhibition of GLI1.

A large-scale exome array analysis of venous thromboembolism

Although recent Genome-Wide Association Studies have identified novel associations for common variants, there has been no comprehensive exome-wide search for low-frequency variants that affect the risk of venous thromboembolism (VTE). We conducted a meta-analysis of 11 studies comprising 8,332 cases and 16,087 controls of European ancestry and 382 cases and 1,476 controls of African American ancestry genotyped with the Illumina HumanExome BeadChip. We used the seqMeta package in R to conduct single variant and gene-based rare variant tests. In the single variant analysis, we limited our analysis to the 64,794 variants with at least 40 minor alleles across studies (minor allele frequency [MAF] ~0.08%). We confirmed associations with previously identified VTE loci, including ABO, F5, F11, and FGA. After adjusting for multiple testing, we observed no novel significant findings in single variant or gene-based analysis. Given our sample size, we had greater than 80% power to detect minimum odds ratios greater than 1.5 and 1.8 for a single variant with MAF of 0.01 and 0.005, respectively. Larger studies and sequence data may be needed to identify novel low-frequency and rare variants associated with VTE risk.

Human iPSC-Derived Retinal Pigment Epithelium: A Model System for Prioritizing and Functionally Characterizing Causal Variants at AMD Risk Loci

We evaluate whether human induced pluripotent stem cell-derived retinal pigment epithelium (iPSC-RPE) cells can be used to prioritize and functionally characterize causal variants at age-related macular degeneration (AMD) risk loci. We generated iPSC-RPE from six subjects and show that they have morphological and molecular characteristics similar to those of native RPE. We generated RNA-seq, ATAC-seq, and H3K27ac ChIP-seq data and observed high similarity in gene expression and enriched transcription factor motif profiles between iPSC-RPE and human fetal RPE. We performed fine mapping of AMD risk loci by integrating molecular data from the iPSC-RPE, adult retina, and adult RPE, which identified rs943080 as the probable causal variant at VEGFA. We show that rs943080 is associated with altered chromatin accessibility of a distal ATAC-seq peak, decreased overall gene expression of VEGFA, and allele-specific expression of a non-coding transcript. Our study thus provides a potential mechanism underlying the association of the VEGFA locus with AMD.

Effect of prothrombotic genotypes on the risk of venous thromboembolism in patients with and without ischemic stroke. The Tromsø Study

Essentials Prothrombotic genotypes may agument the risk of venous thromboembolism (VTE) after ischemic stroke. We studied this effect in a case-cohort study using a genetic risk score. In stroke patients, a one-category increase in the genetic risk score was associated with a 50% higher relative risk of VTE. The risk of VTE in stroke patients increased with an increasing number of risk alleles. SUMMARY: Background Patients with ischemic stroke have a transiently increased risk of subsequent venous thromboembolism (VTE). Prothrombotic genotypes may augment VTE risk under conditions of high thrombosis risk related to stroke. Aims To investigate the effect of prothrombotic genotypes in patients with ischemic stroke on the risk of VTE in a population-based case-cohort study. Methods Cases with incident VTE (n = 664) and a randomly selected age-weighted subcohort (n = 1817) were sampled from three surveys of the Tromsø Study (1994-2008). Participants were genotyped for ABO (rs8176719), F5 (rs6025), F2 (rs1799963), FGG (rs2066865) and F11 (rs2036914) single-nucleotide polymorphisms (SNPs). Cox regression models were used to calculate hazard ratios (HRs) for incident VTE according to individual SNPs and categories of risk alleles (5-SNP score; 0-1, 2, 3-4 and ≥ 5) in participants with and without ischemic stroke. Results There were 192 patients with incident stroke, of whom 43 developed VTE during a median of 15.2 years of follow-up. The risk alleles of individual SNPs augmented the elevated VTE risk brought about by ischemic stroke. In stroke patients, a one-category increase in the genetic risk score was associated with a 50% higher relative risk of overall VTE (HR 1.5, 95% confidence interval [CI] 1.3-1.8) and an 80% higher relative risk of provoked VTE (HR 1.8, 95% CI 1.5-2.1). Stroke patients with ≥ 5 risk alleles had a 12-fold (HR 11.7, 95% CI 4.1-33.3) higher relative risk of VTE than stroke-free participants with 0-1 risk alleles. Conclusions Prothrombotic genotypes increased the risk of VTE in stroke patients, and the risk increased with an increasing number of risk alleles.

Subtle changes in chromatin loop contact propensity are associated with differential gene regulation and expression

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.

Discovery of novel plasma biomarkers for future incident venous thromboembolism by untargeted synchronous precursor selection mass spectrometry proteomics

Essentials Discovery of predictive biomarkers of venous thromboembolism (VTE) may aid risk stratification. A case-control study where plasma was sampled before the occurrence of VTE was established. We generated untargeted plasma proteomic profiles of 200 individuals by use of mass spectrometry. Assessment of the biomarker potential of 501 proteins yielded 46 biomarker candidates.

ABSTRACT

Background Prophylactic anticoagulant treatment may substantially reduce the incidence of venous thromboembolism (VTE) but entails considerable risk of severe bleeding. Identification of individuals at high risk of VTE through the use of predictive biomarkers is desirable in order to achieve a favorable benefit-to-harm ratio. Objective We aimed to identify predictive protein biomarker candidates of VTE. Methods We performed a case-control study of 200 individuals that participated in the Tromsø Study, a population-based cohort, where blood samples were collected before the VTE events occurred. Untargeted tandem mass tag-synchronous precursor selection-mass spectrometry (TMT-SPS-MS3)-based proteomic profiling was used to study the plasma proteomes of each individual. Results Of the 501 proteins detected in a sufficient number of samples to allow multivariate analysis, 46 proteins were associated with VTE case-control status with P-values below the 0.05 significance threshold. The strongest predictive biomarker candidates, assessed by statistical significance, were transthyretin, vitamin K-dependent protein Z and protein/nucleic acid deglycase DJ-1. Conclusions Our untargeted approach of plasma proteome profiling revealed novel predictive biomarker candidates of VTE and confirmed previously reported candidates, thereby providing conceptual support for the validity of the study. A larger nested case-control study will be conducted to validate our findings.

Cell-Surface Marker Signature for Enrichment of Ventricular Cardiomyocytes Derived from Human Embryonic Stem Cells

To facilitate understanding of human cardiomyocyte (CM) subtype specification, and the study of ventricular CM biology in particular, we developed a broadly applicable strategy for enrichment of ventricular cardiomyocytes (VCMs) derived from human embryonic stem cells (hESCs). A bacterial artificial chromosome transgenic H9 hESC line in which GFP expression was driven by the human ventricular-specific myosin light chain 2 (MYL2) promoter was generated, and screened to identify cell-surface markers specific for MYL2-GFP-expressing VCMs. A CD77⁺/CD200^- cell-surface signature facilitated isolation of >97% cardiac troponin I-positive cells from H9 hESC differentiation cultures, with 65% expressing MYL2-GFP. This study provides a tool for VCM enrichment when using some, but not all, human pluripotent stem cell lines. Tools generated in this study can be utilized toward understanding CM subtype specification, and enriching for VCMs for therapeutic applications.

Insights into the Mutational Burden of Human Induced Pluripotent Stem Cells from an Integrative Multi-Omics Approach

To understand the mutational burden of human induced pluripotent stem cells (iPSCs), we sequenced genomes of 18 fibroblast-derived iPSC lines and identified different classes of somatic mutations based on structure, origin, and frequency. Copy-number alterations affected 295 kb in each sample and strongly impacted gene expression. UV-damage mutations were present in ∼45% of the iPSCs and accounted for most of the observed heterogeneity in mutation rates across lines. Subclonal mutations (not present in all iPSCs within a line) composed 10% of point mutations and, compared with clonal variants, showed an enrichment in active promoters and increased association with altered gene expression. Our study shows that, by combining WGS, transcriptome, and epigenome data, we can understand the mutational burden of each iPSC line on an individual basis and suggests that this information could be used to prioritize iPSC lines for models of specific human diseases and/or transplantation therapy.

IFT88 mutations identified in individuals with non-syndromic recessive retinal degeneration result in abnormal ciliogenesis

Whole genome sequencing (WGS) was performed to identify the variants responsible for inherited retinal degeneration (IRD) in a Caucasian family. Segregation analysis of selected rare variants with pathogenic potential identified a set of compound heterozygous changes p.Arg266*:c.796C>T and p.Ala568Thr:c.1702G>A in the intraflagellar transport protein-88 (IFT88) gene segregating with IRD. Expression of IFT88 with the p.Arg266* and p.Ala568Thr mutations in mIMDC3 cells by transient transfection and in HeLa cells by introducing the mutations using CRISPR-cas9 system suggested that both mutations result in the formation of abnormal ciliary structures. The introduction of the IFT88 p.Arg266* variant in the homozygous state in HeLa cells by CRISPR-Cas9 genome-editing revealed that the mutant transcript undergoes nonsense-mediated decay leading to a significant depletion of IFT88 transcript. Additionally, abnormal ciliogenesis was observed in these cells. These observations suggest that the rare and unique combination of IFT88 alleles observed in this study provide insight into the physiological role of IFT88 in humans and the likely mechanism underlying retinal pathology in the pedigree with IRD.

Publisher Correction: PTEN regulates glioblastoma oncogenesis through chromatin-associated complexes of DAXX and histone H3.3

This corrects the article DOI: 10.1038/ncomms15223.

A mutation in IFT43 causes non-syndromic recessive retinal degeneration

The aim of this work is to identify the molecular cause of autosomal recessive early onset retinal degeneration in a consanguineous pedigree. Seventeen members of a four-generation Pakistani family were recruited and underwent a detailed ophthalmic examination. Exomes of four affected and two unaffected individuals were sequenced. Variants were filtered using exomeSuite to identify rare potentially pathogenic variants in genes expressed in the retina and/or brain and consistent with the pattern of inheritance. Effect of the variant observed in the gene Intraflagellar Transport Protein 43 (IFT43) was studied by heterologous expression in mIMCD3 and MDCK cells. Expression and sub-cellular localization of IFT43 in the retina and transiently transfected cells was examined by RT-PCR, western blot analysis, and immunohistochemistry. Affected members were diagnosed with early onset non-syndromic progressive retinal degeneration and the presence of bone spicules distributed throughout the retina at younger ages while the older affected members showed severe central choroidal atrophy. Whole-exome sequencing analysis identified a novel homozygous c.100 G > A change in IFT43 segregating with retinal degeneration and not present in ethnicity-matched controls. Immunostaining showed IFT43 localized in the photoreceptors, and to the tip of the cilia in transfected mIMCD3 and MDCK cells. The cilia in mIMCD3 and MDCK cells expressing mutant IFT43 were found to be significantly shorter (P < 0.001) than cells expressing wild-type IFT43. Our studies identified a novel homozygous mutation in the ciliary protein IFT43 as the underlying cause of recessive inherited retinal degeneration. This is the first report demonstrating the involvement of IFT43 in retinal degeneration.

Joint effects of prothrombotic genotypes and body height on the risk of venous thromboembolism: the Tromsø study

Essentials Body height and prothrombotic genotypes are associated with risk of venous thromboembolism (VTE). The joint effect of prothrombotic genotypes and tall stature on VTE risk is scarcely investigated. We investigated the joint effect of prothrombotic genotypes and tall stature on VTE risk. Prothrombotic genotypes did not yield excess risk of VTE in subjects with a tall stature.

SUMMARY

Background Studies have reported synergistic effects of prothrombotic single-nucleotide polymorphisms (SNPs) and obesity on the risk of venous thromboembolism (VTE). Tall stature is associated with an increased VTE risk, but the joint effect of prothrombotic genotypes and tall stature on the VTE risk is unknown. Aims To investigate the joint effects of prothrombotic genotypes and tall stature on the VTE risk. Methods Cases with incident VTE (n = 676) and a randomly selected age-weighted subcohort (n = 1842) were sampled from the Tromsø study (cohort follow-up: 1994-2012). DNA was genotyped for rs6025 (factor V Leiden), rs1799963 (FII), rs8176719 (ABO blood group), rs2066865 (fibrinogen-γ), and rs2036914 (FIX). Age-adjusted and sex-adjusted hazard ratios (HRs) of VTE were calculated by categories of risk alleles (de Haan 5-SNP score: 0-1, 2-3, and ≥ 4) and body height (< 40th, 40th-80th and > 80th percentiles). Results The VTE risk increased by increasing category of body height, and subjects with height ≥ 178 cm had a two-fold higher VTE risk (HR 2.03; 95% confidence interval [CI] 1.51-2.73) than those with height ≤ 165 cm. The VTE risk also increased across categories of risk alleles. However, the combination of a tall stature and risk alleles, either individual SNPs or risk score, did not result in an excess VTE risk. Subjects with four or more risk alleles and height ≥ 178 cm had a two-fold (HR 2.08; 95% CI 1.24-3.52) higher VTE risk than subjects ≤ 165 cm with no risk allele or one risk allele. Conclusions In contrast to obesity, the presence of prothrombotic genotypes did not result in an excess VTE risk in subjects with a tall stature.