Search

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

We examined genetic and epigenetic changes that occur during disease progression from indolent to aggressive forms of chronic lymphocytic leukemia (CLL) using serial samples from 27 patients. Analysis of DNA mutations grouped the leukemia cases into three categories: evolving (26%), expanding (26%) and static (47%). Thus, approximately three-quarters of the CLL cases had little to no genetic subclonal evolution. However, we identified significant recurrent DNA methylation changes during progression at 4752 CpGs enriched for regions near Polycomb 2 repressive complex (PRC2) targets. Progression-associated CpGs near the PRC2 targets undergo methylation changes in the same direction during disease progression as during normal development from naive to memory B cells. Our study shows that CLL progression does not typically occur via subclonal evolution, but that certain CpG sites undergo recurrent methylation changes. Our results suggest CLL progression may involve developmental processes shared in common with the generation of normal memory B cells.

Introduction: The increasing number of targeted therapies, together with a deeper understanding of cancer genetics and drug response, have prompted major healthcare centers to implement personalized treatment approaches relying on high-throughput tumor DNA sequencing. However, the optimal way to implement this transformative methodology is not yet clear. Current assays may miss important clinical information such as the mutation allelic fraction, the presence of sub-clones or chromosomal rearrangements, or the distinction between inherited variants and somatic mutations. Here, we present the evaluation of ultra-deep targeted sequencing (UDT-Seq) to generate and interpret the molecular profile of 38 breast cancer patients from two academic medical centers.Methods: We sequenced 47 genes in matched germline and tumor DNA samples from 38 breast cancer patients. The selected genes, or the pathways they belong to, can be targeted by drugs or are important in familial cancer risk or drug metabolism.Results: Relying on the added value of sequencing matched tumor and germline DNA and using a dedicated analysis, UDT-Seq has a high sensitivity to identify mutations in tumors with low malignant cell content. Applying UDT-Seq to matched tumor and germline specimens from the 38 patients resulted in a proposal for at least one targeted therapy for 22 patients, the identification of tumor sub-clones in 3 patients, the suggestion of potential adverse drug effects in 3 patients and a recommendation for genetic counseling for 2 patients.Conclusion: Overall our study highlights the additional benefits of a sequencing strategy, which includes germline DNA and is optimized for heterogeneous tumor tissues. © 2013 Harismendy et al.; licensee BioMed Central Ltd.

A major laminin-binding protein (LBP), distinct from previously described LBPs, has been isolated from chick and rat skeletal muscle (Mr 56,000 and 66,000, respectively). The purified LBPs from the two species were shown to be related antigenically and to have similar NH2-terminal amino acid sequences and total amino acid compositions. Protein blots using laminin and laminin fragments provided evidence that this LBP interacts with the major heparin-binding domain, E3, of laminin. Studies on the association of this LBP with muscle membrane fractions and reconstituted lipid vesicles indicate that this protein can interact with lipid bilayers and has properties of a peripheral, not an integral membrane protein. These properties are consistent with its amino acid sequence, determined from cDNAs (Clegg et al., 1988). Examination by light and electron microscopy of the LBP antigen distribution in skeletal muscle indicated that the protein is localized primarily extracellularly, near the extracellular matrix and myotube plasmalemma. While a form of this LBP has been identified in heart muscle, it is present at low or undetectable levels in other tissues examined by immunocytochemistry indicating that it is probably a muscle-specific protein. As this protein is localized extracellularly and can bind to both membranes and laminin, it may mediate myotube interactions with the extracellular matrix.

Payment of residential property rates continues to be a contentious issue among the leaseholders in South Africa. The study investigates the payment practices on residential properties in one municipality. The purpose of the study was to identify the causes of defaults in payment of rates and to proffer solutions to minimise defaults and to recover debts. Systems and contingency theories underpinned the study. A mixed method approach involving both qualitative and quantitative data was applied. The population comprised of five local municipalities, 185 268 residential leaseholders and 437 municipal workers. Data were collected using questionnaires and semi-structured interviews. Survey data were analysed using Statistical Package for the Social Science while data collected from the interviews were transcribed, coded and discussed. The study concluded that defaults in payment were, inter alia, caused by failure of municipalities to impress upon residents the importance of paying rates and to promote alternative payment options.Keywords: Debt Management, Communication, Municipal Financial AdministrationJEL Classifications: D83, H63

Genome analysis of individuals affected by retinitis pigmentosa (RP) identified two rare nucleotide substitutions at the same genomic location on chromosome 11 (g.61392563 [GRCh38]), 69 base pairs upstream of the start codon of the ciliopathy gene TMEM216 (c.-69G>A, c.-69G>T [GenBank: NM_001173991.3]), in individuals of South Asian and African ancestry, respectively. Genotypes included 71 homozygotes and 3 mixed heterozygotes in trans with a predicted loss-of-function allele. Haplotype analysis showed single-nucleotide variants (SNVs) common across families, suggesting ancestral alleles within the two distinct ethnic populations. Clinical phenotype analysis of 62 available individuals from 49 families indicated a similar clinical presentation with night blindness in the first decade and progressive peripheral field loss thereafter. No evident systemic ciliopathy features were noted. Functional characterization of these variants by luciferase reporter gene assay showed reduced promotor activity. Nanopore sequencing confirmed the lower transcription of the TMEM216 c.-69G>T allele in blood-derived RNA from a heterozygous carrier, and reduced expression was further recapitulated by qPCR, using both leukocytes-derived RNA of c.-69G>T homozygotes and total RNA from genome-edited hTERT-RPE1 cells carrying homozygous TMEM216 c.-69G>A. In conclusion, these variants explain a significant proportion of unsolved cases, specifically in individuals of African ancestry, suggesting that reduced TMEM216 expression might lead to abnormal ciliogenesis and photoreceptor degeneration., Competing Interests: Declaration of interests All the authors declare no competing interests., (Copyright © 2024 The Authors. Published by Elsevier Inc. All rights reserved.)

Telomeres are the protective nucleoprotein structures at the end of linear eukaryotic chromosomes. Telomeres' repetitive nature and length have traditionally challenged the precise assessment of the composition and length of individual human telomeres. Here, we present Telo-seq to resolve bulk, chromosome arm-specific and allele-specific human telomere lengths using Oxford Nanopore Technologies' native long-read sequencing. Telo-seq resolves telomere shortening in five population doubling increments and reveals intrasample, chromosome arm-specific, allele-specific telomere length heterogeneity. Telo-seq can reliably discriminate between telomerase- and ALT-positive cancer cell lines. Thus, Telo-seq is a tool to study telomere biology during development, aging, and cancer at unprecedented resolution., (© 2024. The Author(s).)

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., (© 2024. The Author(s).)

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., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2023 The Authors. Published by Elsevier Inc. All rights reserved.)

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., (© 2023. The Author(s).)

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., (© 2023. The Author(s).)

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., Competing Interests: Declaration of competing interest AIdH is currently an employee of AbbVie; LDO, HHC, HC and LH are employees of Genentech; SK is an employee of Gemini Therapeutics. GSH is a shareholder, consultant, and co-founder of Perceive Biotherapeutics, Inc. and an inventor on patents and patent applications owned by the University of Iowa and the University of Utah. RFM, APV, TS, FG, JLH, JJW, SJT, RA, DS, CCWK, JDB, KAF, BHFW and MBG have no competing interests to declare., (Copyright © 2022 The Author(s). Published by Elsevier Ltd.. All rights reserved.)

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., Competing Interests: K.J.G. is a consultant of Genentech and holds stock in Neurocrine Biosciences. B.R. is a consultant of Arima Genomics and a co-founder of Epigenome Technologies. P.B. is an employee of Shoreline Bioscience. N.N. is an employee of Guardant Health. E.B. is an employee and shareholder of Aetion. K.K. is an employee of Cartography Bio. Y.Q. is an employee of Sana Biotechnology. M.D. is an employee and shareholder of Seer. J.C. is an employee and shareholder of Pfizer., (© 2022.)

Reactivation of fetal-specific genes and isoforms occurs during heart failure. However, the underlying molecular mechanisms and the extent to which the fetal program switch occurs remains unclear. Limitations hindering transcriptome-wide analyses of alternative splicing differences (i.e. isoform switching) in cardiovascular system (CVS) tissues between fetal, healthy adult and heart failure have included both cellular heterogeneity across bulk RNA-seq samples and limited availability of fetal tissue for research. To overcome these limitations, we have deconvoluted the cellular compositions of 996 RNA-seq samples representing heart failure, healthy adult (heart and arteria), and fetal-like (iPSC-derived cardiovascular progenitor cells) CVS tissues. Comparison of the expression profiles revealed that reactivation of fetal-specific RNA-binding proteins (RBPs), and the accompanied re-expression of 1,523 fetal-specific isoforms, contribute to the transcriptome differences between heart failure and healthy adult heart. Of note, isoforms for 20 different RBPs were among those that reverted in heart failure to the fetal-like expression pattern. We determined that, compared with adult-specific isoforms, fetal-specific isoforms encode proteins that tend to have more functions, are more likely to harbor RBP binding sites, have canonical sequences at their splice sites, and contain typical upstream polypyrimidine tracts. Our study suggests that compared with healthy adult, fetal cardiac tissue requires stricter transcriptional regulation, and that during heart failure reversion to this stricter transcriptional regulation occurs. Furthermore, we provide a resource of cardiac developmental stage-specific and heart failure-associated genes and isoforms, which are largely unexplored and can be exploited to investigate novel therapeutics for heart failure., Competing Interests: The authors have declared that no competing interests exist.

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.

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., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2021 The Author(s). Published by Elsevier Inc. All rights reserved.)

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., Competing Interests: The authors have declared that no competing interests exist.

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., (© The Author(s) 2020. Published by Oxford University Press for the Infectious Diseases Society of America. All rights reserved. For permissions, e-mail: journals.permissions@oup.com.)

The functional engagement between an enhancer and its target promoter ensures precise gene transcription 1 . 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)., (© 2021. The Author(s), under exclusive licence to Springer Nature Limited.)

Many sequence variants have been linked to complex human traits and diseases 1 , 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.

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 well 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 new 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.

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&#95;112752218 del) resulting in c.1297&#95;1451del; p.K433&#95;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., (© 2020 Wiley Periodicals LLC.)

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., Competing Interests: Competing interestsAuthors are co-founders and shareholders in ATGCures, Inc. Otherwise, authors declare no conflicts or competing interests., (Copyright © 2020 The Authors; exclusive licensee Bio-protocol LLC.)

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., (Copyright © 2020 Elsevier Inc. All rights reserved.)

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

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., (Copyright© 2020 Ferrata Storti Foundation.)

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., Competing Interests: Declarations of Competing Interest None., (Copyright © 2020 The Author(s). Published by Elsevier B.V. All rights reserved.)

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.

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

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.

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., Competing Interests: Competing interestsAuthors declare no competing interests., (Copyright © 2019 The Authors; exclusive licensee Bio-protocol LLC.)

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., (Copyright © 2019 The Authors. Published by Elsevier Inc. All rights reserved.)

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., Competing Interests: MD, JR, DJ, MD, MB, HM, OS, NN, AD, KF No competing interests declared, (© 2019, D'Antonio et al.)