iPSCORE: A Resource of 222 iPSC Lines Enabling Functional Characterization of Genetic Variation across a Variety of Cell Types

Protocol for generating human assembloids to investigate thalamocortical and corticothalamic synaptic transmission and plasticity

Human induced pluripotent stem cells (hiPSCs) can be used to generate assembloids that recreate thalamocortical circuitry displaying short-term and long-term synaptic plasticity. Here, we describe a protocol for differentiating hiPSCs into thalamic and cortical organoids and then fusing them to generate thalamocortical assembloids. We detail the steps for using whole-cell patch-clamp electrophysiology to investigate the properties of synaptic transmission and synaptic plasticity in this model system. For complete details on the use and execution of this protocol, please refer to Patton et al.1.

Multiomic QTL mapping reveals phenotypic complexity of GWAS loci and prioritizes putative causal variants

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 map eQTLs, chromatin accessibility QTLs (caQTLs), and histone acetylation QTLs (haQTLs) using molecular samples from three early developmental-like tissues. Through colocalization, we annotate 10.4% (n = 540) of GWAS loci in 15 traits by QTL phenotype, temporal specificity, and complexity. We show that integration of chromatin QTLs results in a 2.3-fold higher annotation rate of GWAS loci because they capture distal GWAS loci missed by eQTLs, and that 5.4% (n = 13) of GWAS colocalizing eQTLs are early developmental specific. Finally, we utilize the iPSCORE multiomic QTLs to prioritize putative causal variants overlapping transcription factor motifs to elucidate the potential genetic underpinnings of 296 GWAS-QTL colocalizations.

Vinculin haploinsufficiency impairs integrin-mediated costamere remodeling on stiffer microenvironments

Vinculin (VCL) is a key adapter protein located in force-bearing costamere complexes, which mechanically couples the sarcomere to the ECM. Heterozygous vinculin frameshift genetic variants can contribute to cardiomyopathy when external stress is applied, but the mechanosensitive pathways underpinning VCL haploinsufficiency remain elusive. Here, we show that in response to extracellular matrix stiffening, heterozygous loss of VCL disrupts force-mediated costamere protein recruitment, thereby impairing cardiomyocyte contractility and sarcomere organization. Analyses of human pluripotent stem cell-derived cardiomyocytes (hPSC-CMs) harboring either VCL c.659dupA or VCL c.74del7 heterozygous VCL frameshift variants revealed that these VCL mutant hPSC-CMs exhibited heightened contractile strain energy, morphological maladaptation, and sarcomere disarray on stiffened matrix. Mechanosensitive recruitment of costameric talin 2, paxillin, focal adhesion kinase, and α-actinin was significantly reduced in vinculin variant cardiomyocytes. Despite poorly formed costamere complexes and sarcomeres, elevated expression of integrin β1 and cortical actin on stiff substrates may rescue force transmission on stiff substrates, an effect that is recapitulated in WT CMs by ligating integrin receptors and blocking mechanosensation. Together, these data support that heterozygous loss of VCL contributes to adverse cardiomyocyte remodeling by impairing adhesion-mediated force transmission from the costamere to the cytoskeleton. (191 words).

Cross cell-type systems genetics reveals the influence of eQTL at multiple points in the developmental trajectory of mouse neural progenitor cells

Genetic variation leads to phenotypic variability in pluripotent stem cells that presents challenges for regenerative medicine. Although recent studies have investigated the impact of genetic variation on pluripotency maintenance and differentiation capacity, less is known about how genetic variants affecting the pluripotent state influence gene regulation in later stages of development. Here, we characterized expression of more than 12,000 genes in 127 donor-matched Diversity Outbred (DO) mouse embryonic stem cell (mESC) and neural progenitor cell (mNPC) lines. Quantitative trait locus (QTL) mapping identified 2,947 expression QTL (eQTL) unique to DO mNPCs and 1,113 eQTL observed in both mNPCs and mESCs with highly concordant allele effects. We mapped three eQTL hotspots on Chromosomes (Chrs) 1, 10, and 11 that were unique to mNPCs. Target genes of the Chr 1 hotspot were overrepresented for those involved in mRNA processing, DNA repair, chromatin organization, protein degradation, and cell cycle. Mediation analysis of the Chr 1 hotspot identified Rnf152 as the best candidate mediator expressed in mNPCs, while cross-cell type mediation using mESC gene expression along with partial correlation analysis strongly implicated genetic variant(s) affecting Pign expression in the mESC state as regulating the mNPC Chr 1 eQTL hotspot. Together these findings highlight that many local eQTL confer similar effects on gene expression in multiple cell states; distant eQTL in DO mNPCs are numerous and largely unique to that cell state, with many co-localizing to mNPC-specific hotspots; and mediation analysis across cell types suggests that expression of Pign early in development (mESCs) shapes the transcriptome of the more specialized mNPC state.

xQTLatlas: a comprehensive resource for human cellular-resolution multi-omics genetic regulatory landscape

Understanding how genetic variants influence molecular phenotypes in different cellular contexts is crucial for elucidating the molecular and cellular mechanisms behind complex traits, which in turn has spurred significant advances in research into molecular quantitative trait locus (xQTL) at the cellular level. With the rapid proliferation of data, there is a critical need for a comprehensive and accessible platform to integrate this information. To meet this need, we developed xQTLatlas (http://www.hitxqtl.org.cn/), a database that provides a multi-omics genetic regulatory landscape at cellular resolution. xQTLatlas compiles xQTL summary statistics from 151 cell types and 339 cell states across 55 human tissues. It organizes these data into 20 xQTL types, based on four distinct discovery strategies, and spans 13 molecular phenotypes. Each entry in xQTLatlas is meticulously annotated with comprehensive metadata, including the origin of the tissue, cell type, cell state and the QTL discovery strategies utilized. Additionally, xQTLatlas features multiscale data exploration tools and a suite of interactive visualizations, facilitating in-depth analysis of cell-level xQTL. xQTLatlas provides a valuable resource for deepening our understanding of the impact of functional variants on molecular phenotypes in different cellular environments, thereby facilitating extensive research efforts.

Efficient site-specific integration of large genes in mammalian cells via continuously evolved recombinases and prime editing

Methods for the targeted integration of genes in mammalian genomes suffer from low programmability, low efficiencies or low specificities. Here we show that phage-assisted continuous evolution enhances prime-editing-assisted site-specific integrase gene editing (PASSIGE), which couples the programmability of prime editing with the ability of recombinases to precisely integrate large DNA cargoes exceeding 10 kilobases. Evolved and engineered Bxb1 recombinase variants (evoBxb1 and eeBxb1) mediated up to 60% donor integration (3.2-fold that of wild-type Bxb1) in human cell lines with pre-installed recombinase landing sites. In single-transfection experiments at safe-harbour and therapeutically relevant sites, PASSIGE with eeBxb1 led to an average targeted-gene-integration efficiencies of 23% (4.2-fold that of wild-type Bxb1). Notably, integration efficiencies exceeded 30% at multiple sites in primary human fibroblasts. PASSIGE with evoBxb1 or eeBxb1 outperformed PASTE (for 'programmable addition via site-specific targeting elements', a method that uses prime editors fused to recombinases) on average by 9.1-fold and 16-fold, respectively. PASSIGE with continuously evolved recombinases is an unusually efficient method for the targeted integration of genes in mammalian cells.

Leveraging large-scale multi-omics evidences to identify therapeutic targets from genome-wide association studies

BACKGROUND

Therapeutic targets supported by genetic evidence from genome-wide association studies (GWAS) show higher probability of success in clinical trials. GWAS is a powerful approach to identify links between genetic variants and phenotypic variation; however, identifying the genes driving associations identified in GWAS remains challenging. Integration of molecular quantitative trait loci (molQTL) such as expression QTL (eQTL) using mendelian randomization (MR) and colocalization analyses can help with the identification of causal genes. Careful interpretation remains warranted because eQTL can affect the expression of multiple genes within the same locus.

METHODS

We used a combination of genomic features that include variant annotation, activity-by-contact maps, MR, and colocalization with molQTL to prioritize causal genes across 4,611 disease GWAS and meta-analyses from biobank studies, namely FinnGen, Estonian Biobank and UK Biobank.

RESULTS

Genes identified using this approach are enriched for gold standard causal genes and capture known biological links between disease genetics and biology. In addition, we find that eQTL colocalizing with GWAS are statistically enriched for corresponding disease-relevant tissues. We show that predicted directionality from MR is generally consistent with matched drug mechanism of actions (> 85% for approved drugs). Compared to the nearest gene mapping method, genes supported by multi-omics evidences displayed higher enrichment in approved therapeutic targets (risk ratio 1.75 vs. 2.58 for genes with the highest level of support). Finally, using this approach, we detected anassociation between the IL6 receptor signal transduction gene IL6ST and polymyalgia rheumatica, an indication for which sarilumab, a monoclonal antibody against IL-6, has been recently approved.

CONCLUSIONS

Combining variant annotation, activity-by-contact maps, and molQTL increases performance to identify causal genes, while informing on directionality which can be translated to successful target identification and drug development.

Extensive co-regulation of neighboring genes complicates the use of eQTLs in target gene prioritization

Identifying causal genes underlying genome-wide association studies (GWASs) is a fundamental problem in human genetics. Although colocalization with gene expression quantitative trait loci (eQTLs) is often used to prioritize GWAS target genes, systematic benchmarking has been limited due to unavailability of large ground truth datasets. Here, we re-analyzed plasma protein QTL data from 3,301 individuals of the INTERVAL cohort together with 131 eQTL Catalog datasets. Focusing on variants located within or close to the affected protein identified 793 proteins with at least one cis-pQTL where we could assume that the most likely causal gene was the gene coding for the protein. We then benchmarked the ability of cis-eQTLs to recover these causal genes by comparing three Bayesian colocalization methods (coloc.susie, coloc.abf, and CLPP) and five Mendelian randomization (MR) approaches (three varieties of inverse-variance weighted MR, MR-RAPS, and MRLocus). We found that assigning fine-mapped pQTLs to their closest protein coding genes outperformed all colocalization methods regarding both precision (71.9%) and recall (76.9%). Furthermore, the colocalization method with the highest recall (coloc.susie - 46.3%) also had the lowest precision (45.1%). Combining evidence from multiple conditionally distinct colocalizing QTLs with MR increased precision to 81%, but this was accompanied by a large reduction in recall to 7.1%. Furthermore, the choice of the MR method greatly affected performance, with the standard inverse-variance-weighted MR often producing many false positives. Our results highlight that linking GWAS variants to target genes remains challenging with eQTL evidence alone, and prioritizing novel targets requires triangulation of evidence from multiple sources.

Advancing Parkinson's disease treatment: cell replacement therapy with neurons derived from pluripotent stem cells

The motor symptoms of Parkinson's disease (PD) are caused by the progressive loss of dopamine neurons from the substantia nigra. There are currently no treatments that can slow or reverse the neurodegeneration. To restore the lost neurons, international groups have initiated clinical trials using human embryonic or induced pluripotent stem cells (PSCs) to derive dopamine neuron precursors that are used as transplants to replace the lost neurons. Proof-of-principle experiments in the 1980s and 1990s showed that grafts of fetal ventral mesencephalon, which contains the precursors of the substantial nigra, could, under rare circumstances, reverse symptoms of the disease. Improvements in PSC technology and genomics have inspired researchers to design clinical trials using PSC-derived dopamine neuron precursors as cell replacement therapy for PD. We focus here on 4 such first-in-human clinical trials that have begun in the US, Europe, and Japan. We provide an overview of the sources of PSCs and the methods used to generate cells for transplantation. We discuss pros and cons of strategies for allogeneic, immune-matched, and autologous approaches and novel methods for overcoming rejection by the immune system. We consider challenges for safety and efficacy of the cells for durable engraftment, focusing on the genomics-based quality control methods to assure that the cells will not become cancerous. Finally, since clinical trials like these have never been undertaken before, we comment on the value of cooperation among rivals to contribute to advancements that will finally provide relief for the millions suffering from the symptoms of PD.

DNA damage-associated protein co-expression network in cardiomyocytes informs on tolerance to genetic variation and disease

Cardiovascular disease (CVD) is associated with both genetic variants and environmental factors. One unifying consequence of the molecular risk factors in CVD is DNA damage, which must be repaired by DNA damage response proteins. However, the impact of DNA damage on global cardiomyocyte protein abundance, and its relationship to CVD risk remains unclear. We therefore treated induced pluripotent stem cell-derived cardiomyocytes with the DNA-damaging agent Doxorubicin (DOX) and a vehicle control, and identified 4,178 proteins that contribute to a network comprising 12 co-expressed modules and 403 hub proteins with high intramodular connectivity. Five modules correlate with DOX and represent distinct biological processes including RNA processing, chromatin regulation and metabolism. DOX-correlated hub proteins are depleted for proteins that vary in expression across individuals due to genetic variation but are enriched for proteins encoded by loss-of-function intolerant genes. While proteins associated with genetic risk for CVD, such as arrhythmia are enriched in specific DOX-correlated modules, DOX-correlated hub proteins are not enriched for known CVD risk proteins. Instead, they are enriched among proteins that physically interact with CVD risk proteins. Our data demonstrate that DNA damage in cardiomyocytes induces diverse effects on biological processes through protein co-expression modules that are relevant for CVD, and that the level of protein connectivity in DNA damage-associated modules influences the tolerance to genetic variation.

Genetic analyses of inflammatory polyneuropathy and chronic inflammatory demyelinating polyradiculoneuropathy identified candidate genes

Chronic inflammatory demyelinating polyneuropathy (CIDP) is a rare, immune-mediated disorder in which an aberrant immune response causes demyelination and axonal damage of the peripheral nerves. Genetic contribution to CIDP is unclear and no genome-wide association study (GWAS) has been reported so far. In this study, we aimed to identify CIDP-related risk loci, genes, and pathways. We first focused on CIDP, and 516 CIDP cases and 403,545 controls were included in the GWAS analysis. We also investigated genetic risk for inflammatory polyneuropathy (IP), in which we performed a GWAS study using FinnGen data and combined the results with GWAS from the UK Biobank using a fixed-effect meta-analysis. A total of 1,261 IP cases and 823,730 controls were included in the analysis. Stratified analyses by gender were performed. Mendelian randomization (MR), colocalization, and transcriptome-wide association study (TWAS) analyses were performed to identify associated genes. Gene-set analyses were conducted to identify associated pathways. We identified one genome-wide significant locus at 20q13.33 for CIDP risk among women, the top variant located at the intron region of gene CDH4. Sex-combined MR, colocalization, and TWAS analyses identified three candidate pathogenic genes for CIDP and five genes for IP. MAGMA gene-set analyses identified a total of 18 pathways related to IP or CIDP. Sex-stratified analyses identified three genes for IP among males and two genes for IP among females. Our study identified suggestive risk genes and pathways for CIDP and IP. Functional analyses should be conducted to further confirm these associations.

On-Chip Neural Induction Boosts Neural Stem Cell Commitment: Toward a Pipeline for iPSC-Based Therapies

The clinical translation of induced pluripotent stem cells (iPSCs) holds great potential for personalized therapeutics. However, one of the main obstacles is that the current workflow to generate iPSCs is expensive, time-consuming, and requires standardization. A simplified and cost-effective microfluidic approach is presented for reprogramming fibroblasts into iPSCs and their subsequent differentiation into neural stem cells (NSCs). This method exploits microphysiological technology, providing a 100-fold reduction in reagents for reprogramming and a ninefold reduction in number of input cells. The iPSCs generated from microfluidic reprogramming of fibroblasts show upregulation of pluripotency markers and downregulation of fibroblast markers, on par with those reprogrammed in standard well-conditions. The NSCs differentiated in microfluidic chips show upregulation of neuroectodermal markers (ZIC1, PAX6, SOX1), highlighting their propensity for nervous system development. Cells obtained on conventional well plates and microfluidic chips are compared for reprogramming and neural induction by bulk RNA sequencing. Pathway enrichment analysis of NSCs from chip showed neural stem cell development enrichment and boosted commitment to neural stem cell lineage in initial phases of neural induction, attributed to a confined environment in a microfluidic chip. This method provides a cost-effective pipeline to reprogram and differentiate iPSCs for therapeutics compliant with current good manufacturing practices.

Optimizing Nodal, Wnt and BMP signaling pathways for robust and efficient differentiation of human induced pluripotent stem cells to intermediate mesoderm cells

Several differentiation protocols have enabled the generation of intermediate mesoderm (IM)-derived cells from human pluripotent stem cells (hPSC). However, the substantial variability between existing protocols for generating IM cells compromises their efficiency, reproducibility, and overall success, potentially hindering the utility of urogenital system organoids. Here, we examined the role of high levels of Nodal signaling and BMP activity, as well as WNT signaling in the specification of IM cells derived from a UCSD167i-99-1 human induced pluripotent stem cells (hiPSC) line. We demonstrate that precise modulation of WNT and BMP signaling significantly enhances IM differentiation efficiency. Treatment of hPSC with 3 μM CHIR99021 induced TBXT+/MIXL1+ mesoderm progenitor (MP) cells after 48 h of differentiation. Further treatment with a combination of 3 μM CHIR99021 and 4 ng/mL BMP4 resulted in the generation of OSR1+/GATA3+/PAX2+ IM cells within a subsequent 48 h period. Molecular characterization of differentiated cells was confirmed through immunofluorescence staining and RT-qPCR. Hence, this study establishes a consistent and reproducible protocol for differentiating hiPSC into IM cells that faithfully recapitulates the molecular signatures of IM development. This protocol holds promise for improving the success of protocols designed to generate urogenital system organoids in vitro, with potential applications in regenerative medicine, drug discovery, and disease modeling.

Recruitment of CTCF to the SIRT1 promoter after Oxidative Stress mediates Cardioprotective Transcription

Because most DNA-binding transcription factors (dbTFs), including the architectural regulator CTCF, bind RNA and exhibit di-/multimerization, a central conundrum is whether these distinct properties are regulated post-transcriptionally to modulate transcriptional programs. Here, investigating stress-dependent activation of SIRT1, encoding an evolutionarily-conserved protein deacetylase, we show that induced phosphorylation of CTCF acts as a rheostat to permit CTCF occupancy of low-affinity promoter DNA sites to precisely the levels necessary. This CTCF recruitment to the SIRT1 promoter is eliciting a cardioprotective cardiomyocyte transcriptional activation program and provides resilience against the stress of the beating heart in vivo . Mice harboring a mutation in the conserved low-affinity CTCF promoter binding site exhibit an altered, cardiomyocyte-specific transcriptional program and a systolic heart failure phenotype. This transcriptional role for CTCF reveals that a covalent dbTF modification regulating signal-dependent transcription serves as a previously unsuspected component of the oxidative stress response.

Personalized transcriptome signatures in a cardiomyopathy stem cell biobank

BACKGROUND

There is growing evidence that pathogenic mutations do not fully explain hypertrophic (HCM) or dilated (DCM) cardiomyopathy phenotypes. We hypothesized that if a patient's genetic background was influencing cardiomyopathy this should be detectable as signatures in gene expression. We built a cardiomyopathy biobank resource for interrogating personalized genotype phenotype relationships in human cell lines.

METHODS

We recruited 308 diseased and control patients for our cardiomyopathy stem cell biobank. We successfully reprogrammed PBMCs (peripheral blood mononuclear cells) into induced pluripotent stem cells (iPSCs) for 300 donors. These iPSCs underwent whole genome sequencing and were differentiated into cardiomyocytes for RNA-seq. In addition to annotating pathogenic variants, mutation burden in a panel of cardiomyopathy genes was assessed for correlation with echocardiogram measurements. Line-specific co-expression networks were inferred to evaluate transcriptomic subtypes. Drug treatment targeted the sarcomere, either by activation with omecamtiv mecarbil or inhibition with mavacamten, to alter contractility.

RESULTS

We generated an iPSC biobank from 300 donors, which included 101 individuals with HCM and 88 with DCM. Whole genome sequencing of 299 iPSC lines identified 78 unique pathogenic or likely pathogenic mutations in the diseased lines. Notably, only DCM lines lacking a known pathogenic or likely pathogenic mutation replicated a finding in the literature for greater nonsynonymous SNV mutation burden in 102 cardiomyopathy genes to correlate with lower left ventricular ejection fraction in DCM. We analyzed RNA-sequencing data from iPSC-derived cardiomyocytes for 102 donors. Inferred personalized co-expression networks revealed two transcriptional subtypes of HCM. The first subtype exhibited concerted activation of the co-expression network, with the degree of activation reflective of the disease severity of the donor. In contrast, the second HCM subtype and the entire DCM cohort exhibited partial activation of the respective disease network, with the strength of specific gene by gene relationships dependent on the iPSC-derived cardiomyocyte line. ADCY5 was the largest hubnode in both the HCM and DCM networks and partially corrected in response to drug treatment.

CONCLUSIONS

We have a established a stem cell biobank for studying cardiomyopathy. Our analysis supports the hypothesis the genetic background influences pathologic gene expression programs and support a role for ADCY5 in cardiomyopathy.

Brain organoids engineered to give rise to glia and neural networks after 90 days in culture exhibit human-specific proteoforms

Human brain organoids are emerging as translationally relevant models for the study of human brain health and disease. However, it remains to be shown whether human-specific protein processing is conserved in human brain organoids. Herein, we demonstrate that cell fate and composition of unguided brain organoids are dictated by culture conditions during embryoid body formation, and that culture conditions at this stage can be optimized to result in the presence of glia-associated proteins and neural network activity as early as three-months in vitro. Under these optimized conditions, unguided brain organoids generated from induced pluripotent stem cells (iPSCs) derived from male-female siblings are similar in growth rate, size, and total protein content, and exhibit minimal batch-to-batch variability in cell composition and metabolism. A comparison of neuronal, microglial, and macroglial (astrocyte and oligodendrocyte) markers reveals that profiles in these brain organoids are more similar to autopsied human cortical and cerebellar profiles than to those in mouse cortical samples, providing the first demonstration that human-specific protein processing is largely conserved in unguided brain organoids. Thus, our organoid protocol provides four major cell types that appear to process proteins in a manner very similar to the human brain, and they do so in half the time required by other protocols. This unique copy of the human brain and basic characteristics lay the foundation for future studies aiming to investigate human brain-specific protein patterning (e.g., isoforms, splice variants) as well as modulate glial and neuronal processes in an in situ-like environment.

Generation of Skeletal Muscle Organoids from Human Pluripotent Stem Cells

Various protocols have been proven effective in the directed differentiation of mouse and human pluripotent stem cells into skeletal muscles and used to study myogenesis. Current 2D myogenic differentiation protocols can mimic muscle development and its alteration under pathological conditions such as muscular dystrophies. 3D skeletal muscle differentiation approaches can, in addition, model the interaction between the various cell types within the developing organoid. Our protocol ensures the differentiation of human embryonic/induced pluripotent stem cells (hESC/hiPSC) into skeletal muscle organoids (SMO) via cells with paraxial mesoderm and neuromesodermal progenitors' identity and further production of organized structures of the neural plate margin and the dermomyotome. Continuous culturing omits neural lineage differentiation and promotes fetal myogenesis, including the maturation of fibroadipogenic progenitors and PAX7-positive myogenic progenitors. The PAX7 progenitors resemble the late fetal stages of human development and, based on single-cell transcriptomic profiling, cluster close to adult satellite cells of primary muscles. To overcome the limited availability of muscle biopsies from patients with muscular dystrophy during disease progression, we propose to use the SMO system, which delivers a stable population of skeletal muscle progenitors from patient-specific iPSCs to investigate human myogenesis in healthy and diseased conditions. Key features • Development of skeletal muscle organoid differentiation from human pluripotent stem cells, which recapitulates myogenesis. • Analysis of early embryonic and fetal myogenesis. • Provision of skeletal muscle progenitors for in vitro and in vivo analysis for up to 14 weeks of organoid culture. • In vitro myogenesis from patient-specific iPSCs allows to overcome the bottleneck of muscle biopsies of patients with pathological conditions.

Integrating population genetics, stem cell biology and cellular genomics to study complex human diseases

Human pluripotent stem (hPS) cells can, in theory, be differentiated into any cell type, making them a powerful in vitro model for human biology. Recent technological advances have facilitated large-scale hPS cell studies that allow investigation of the genetic regulation of molecular phenotypes and their contribution to high-order phenotypes such as human disease. Integrating hPS cells with single-cell sequencing makes identifying context-dependent genetic effects during cell development or upon experimental manipulation possible. Here we discuss how the intersection of stem cell biology, population genetics and cellular genomics can help resolve the functional consequences of human genetic variation. We examine the critical challenges of integrating these fields and approaches to scaling them cost-effectively and practically. We highlight two areas of human biology that can particularly benefit from population-scale hPS cell studies, elucidating mechanisms underlying complex disease risk loci and evaluating relationships between common genetic variation and pharmacotherapeutic phenotypes.

Sex-biased gene expression during neural differentiation of human embryonic stem cells

Sex differences in the developing human brain are primarily attributed to hormonal influence. Recently however, genetic differences and their impact on the developing nervous system have attracted increased attention. To understand genetically driven sexual dimorphisms in neurodevelopment, we investigated genome-wide gene expression in an in vitro differentiation model of male and female human embryonic stem cell lines (hESC), independent of the effects of human sex hormones. Four male and four female-derived hESC lines were differentiated into a population of mixed neurons over 37 days. Differential gene expression and gene set enrichment analyses were conducted on bulk RNA sequencing data. While similar differentiation tendencies in all cell lines demonstrated the robustness and reproducibility of our differentiation protocol, we found sex-biased gene expression already in undifferentiated ESCs at day 0, but most profoundly after 37 days of differentiation. Male and female cell lines exhibited sex-biased expression of genes involved in neurodevelopment, suggesting that sex influences the differentiation trajectory. Interestingly, the highest contribution to sex differences was found to arise from the male transcriptome, involving both Y chromosome and autosomal genes. We propose 13 sex-biased candidate genes (10 upregulated in male cell lines and 3 in female lines) that are likely to affect neuronal development. Additionally, we confirmed gene dosage compensation of X/Y homologs escaping X chromosome inactivation through their Y homologs and identified a significant overexpression of the Y-linked demethylase UTY and KDM5D in male hESC during neuron development, confirming previous results in neural stem cells. Our results suggest that genetic sex differences affect neuronal differentiation trajectories, which could ultimately contribute to sex biases during human brain development.

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.

Strengthening cardiac therapy pipelines using human pluripotent stem cell-derived cardiomyocytes

Advances in hiPSC isolation and reprogramming and hPSC-CM differentiation have prompted their therapeutic application and utilization for evaluating potential cardiovascular safety liabilities. In this perspective, we showcase key efforts toward the large-scale production of hiPSC-CMs, implementation of hiPSC-CMs in industry settings, and recent clinical applications of this technology. The key observations are a need for traceable gender and ethnically diverse hiPSC lines, approaches to reduce cost of scale-up, accessible clinical trial datasets, and transparent guidelines surrounding the safety and efficacy of hiPSC-based therapies.

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.

Gene editing improves endoplasmic reticulum-mitochondrial contacts and unfolded protein response in Friedreich's ataxia iPSC-derived neurons

Friedreich ataxia (FRDA) is a multisystemic, autosomal recessive disorder caused by homozygous GAA expansion mutation in the first intron of frataxin (FXN) gene. FXN is a mitochondrial protein critical for iron-sulfur cluster biosynthesis and deficiency impairs mitochondrial electron transport chain functions and iron homeostasis within the organelle. Currently, there is no effective treatment for FRDA. We have previously demonstrated that single infusion of wild-type hematopoietic stem and progenitor cells (HSPCs) resulted in prevention of neurologic and cardiac complications of FRDA in YG8R mice, and rescue was mediated by FXN transfer from tissue engrafted, HSPC-derived microglia/macrophages to diseased neurons/myocytes. For a future clinical translation, we developed an autologous stem cell transplantation approach using CRISPR/Cas9 for the excision of the GAA repeats in FRDA patients' CD34+ HSPCs; this strategy leading to increased FXN expression and improved mitochondrial functions. The aim of the current study is to validate the efficiency and safety of our gene editing approach in a disease-relevant model. We generated a cohort of FRDA patient-derived iPSCs and isogenic lines that were gene edited with our CRISPR/Cas9 approach. iPSC derived FRDA neurons displayed characteristic apoptotic and mitochondrial phenotype of the disease, such as non-homogenous microtubule staining in neurites, increased caspase-3 expression, mitochondrial superoxide levels, mitochondrial fragmentation, and partial degradation of the cristae compared to healthy controls. These defects were fully prevented in the gene edited neurons. RNASeq analysis of FRDA and gene edited neurons demonstrated striking improvement in gene clusters associated with endoplasmic reticulum (ER) stress in the isogenic lines. Gene edited neurons demonstrated improved ER-calcium release, normalization of ER stress response gene, XBP-1, and significantly increased ER-mitochondrial contacts that are critical for functional homeostasis of both organelles, as compared to FRDA neurons. Ultrastructural analysis for these contact sites displayed severe ER structural damage in FRDA neurons, that was undetected in gene edited neurons. Taken together, these results represent a novel finding for disease pathogenesis showing dramatic ER structural damage in FRDA, validate the efficacy profile of our FXN gene editing approach in a disease relevant model, and support our approach as an effective strategy for therapeutic intervention for Friedreich's ataxia.

Anthracyclines induce cardiotoxicity through a shared gene expression response signature

TOP2 inhibitors (TOP2i) are effective drugs for breast cancer treatment. However, they can cause cardiotoxicity in some women. The most widely used TOP2i include anthracyclines (AC) Doxorubicin (DOX), Daunorubicin (DNR), Epirubicin (EPI), and the anthraquinone Mitoxantrone (MTX). It is unclear whether women would experience the same adverse effects from all drugs in this class, or if specific drugs would be preferable for certain individuals based on their cardiotoxicity risk profile. To investigate this, we studied the effects of treatment of DOX, DNR, EPI, MTX, and an unrelated monoclonal antibody Trastuzumab (TRZ) on iPSC-derived cardiomyocytes (iPSC-CMs) from six healthy females. All TOP2i induce cell death at concentrations observed in cancer patient serum, while TRZ does not. A sub-lethal dose of all TOP2i induces limited cellular stress but affects calcium handling, a function critical for cardiomyocyte contraction. TOP2i induce thousands of gene expression changes over time, giving rise to four distinct gene expression response signatures, denoted as TOP2i early-acute, early-sustained, and late response genes, and non-response genes. There is no drug- or AC-specific signature. TOP2i early response genes are enriched in chromatin regulators, which mediate AC sensitivity across breast cancer patients. However, there is increased transcriptional variability between individuals following AC treatments. To investigate potential genetic effects on response variability, we first identified a reported set of expression quantitative trait loci (eQTLs) uncovered following DOX treatment in iPSC-CMs. Indeed, DOX response eQTLs are enriched in genes that respond to all TOP2i. Next, we identified 38 genes in loci associated with AC toxicity by GWAS or TWAS. Two thirds of the genes that respond to at least one TOP2i, respond to all ACs with the same direction of effect. Our data demonstrate that TOP2i induce thousands of shared gene expression changes in cardiomyocytes, including genes near SNPs associated with inter-individual variation in response to DOX treatment and AC-induced cardiotoxicity.

High-dimensional phenotyping to define the genetic basis of cellular morphology

The morphology of cells is dynamic and mediated by genetic and environmental factors. Characterizing how genetic variation impacts cell morphology can provide an important link between disease association and cellular function. Here, we combine genomic sequencing and high-content imaging approaches on iPSCs from 297 unique donors to investigate the relationship between genetic variants and cellular morphology to map what we term cell morphological quantitative trait loci (cmQTLs). We identify novel associations between rare protein altering variants in WASF2, TSPAN15, and PRLR with several morphological traits related to cell shape, nucleic granularity, and mitochondrial distribution. Knockdown of these genes by CRISPRi confirms their role in cell morphology. Analysis of common variants yields one significant association and nominate over 300 variants with suggestive evidence (P < 10-6) of association with one or more morphology traits. We then use these data to make predictions about sample size requirements for increasing discovery in cellular genetic studies. We conclude that, similar to molecular phenotypes, morphological profiling can yield insight about the function of genes and variants.

A patterned human primitive heart organoid model generated by pluripotent stem cell self-organization

Pluripotent stem cell-derived organoids can recapitulate significant features of organ development in vitro. We hypothesized that creating human heart organoids by mimicking aspects of in utero gestation (e.g., addition of metabolic and hormonal factors) would lead to higher physiological and anatomical relevance. We find that heart organoids produced using this self-organization-driven developmental induction strategy are remarkably similar transcriptionally and morphologically to age-matched human embryonic hearts. We also show that they recapitulate several aspects of cardiac development, including large atrial and ventricular chambers, proepicardial organ formation, and retinoic acid-mediated anterior-posterior patterning, mimicking the developmental processes found in the post-heart tube stage primitive heart. Moreover, we provide proof-of-concept demonstration of the value of this system for disease modeling by exploring the effects of ondansetron, a drug administered to pregnant women and associated with congenital heart defects. These findings constitute a significant technical advance for synthetic heart development and provide a powerful tool for cardiac disease modeling.

Human skeletal muscle organoids model fetal myogenesis and sustain uncommitted PAX7 myogenic progenitors

In vitro culture systems that structurally model human myogenesis and promote PAX7+ myogenic progenitor maturation have not been established. Here we report that human skeletal muscle organoids can be differentiated from induced pluripotent stem cell lines to contain paraxial mesoderm and neuromesodermal progenitors and develop into organized structures reassembling neural plate border and dermomyotome. Culture conditions instigate neural lineage arrest and promote fetal hypaxial myogenesis toward limb axial anatomical identity, with generation of sustainable uncommitted PAX7 myogenic progenitors and fibroadipogenic (PDGFRa+) progenitor populations equivalent to those from the second trimester of human gestation. Single-cell comparison to human fetal and adult myogenic progenitor /satellite cells reveals distinct molecular signatures for non-dividing myogenic progenitors in activated (CD44High/CD98+/MYOD1+) and dormant (PAX7High/FBN1High/SPRY1High) states. Our approach provides a robust 3D in vitro developmental system for investigating muscle tissue morphogenesis and homeostasis.

Genome-wide association study analysis of disease severity in Acne reveals novel biological insights

Acne vulgaris is a common skin disease that affects >85% of teenage young adults among which >8% develop severe lesions that leaves permanent scars. Genetic heritability studies of acne in twin cohorts have estimated that the heritability for acne is 80%. Previous genome-wide association studies (GWAS) have identified 50 genetic loci associated with increased risk of developing acne when compared to healthy individuals. However only a few studies have investigated genetic association with disease severity. GWAS of disease progression may provide a more effective approach to unveil potential disease modifying therapeutic targets. Here, we performed a multi-ethnic GWAS analysis to capture disease severity in acne patients by using individuals with normal acne as a control. Our cohort consists of a total of 2,956 participants, including 290 severe acne cases and 930 normal acne controls from FinnGen, and 522 cases and 1,214 controls from BioVU. We also performed mendelian randomization (MR), colocalization analyses and transcriptome-wide association study (TWAS) to identify putative causal genes. Lastly, we performed gene-set enrichment analysis using MAGMA to implicate biological pathways that drive disease severity in Acne. We identified two new loci associated with acne severity at the genome-wide significance level, six novel associated genes by MR, colocalization and TWAS analyses, including genes CDC7, SLC7A1, ADAM23, TTLL10, CDK20 and DNAJA4 , and 5 novel pathways by MAGMA analyses. Our study suggests that the etiologies of acne susceptibility and severity have limited overlap, with only 26% of known acne risk loci presenting nominal association with acne severity and none of the novel severity associated genes reported as associated with acne risk in previous GWAS.

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.

Reliable multiplex generation of pooled induced pluripotent stem cells

Reprogramming somatic cells into pluripotent stem cells (iPSCs) enables the study of systems in vitro. To increase the throughput of reprogramming, we present induction of pluripotency from pooled cells (iPPC)-an efficient, scalable, and reliable reprogramming procedure. Using our deconvolution algorithm that employs pooled sequencing of single-nucleotide polymorphisms (SNPs), we accurately estimated individual donor proportions of the pooled iPSCs. With iPPC, we concurrently reprogrammed over one hundred donor lymphoblastoid cell lines (LCLs) into iPSCs and found strong correlations of individual donors' reprogramming ability across multiple experiments. Individual donors' reprogramming ability remains consistent across both same-day replicates and multiple experimental runs, and the expression of certain immunoglobulin precursor genes may impact reprogramming ability. The pooled iPSCs were also able to differentiate into cerebral organoids. Our procedure enables a multiplex framework of using pooled libraries of donor iPSCs for downstream research and investigation of in vitro phenotypes.

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.

eQTL Catalogue 2023: New datasets, X chromosome QTLs, and improved detection and visualisation of transcript-level QTLs

The eQTL Catalogue is an open database of uniformly processed human molecular quantitative trait loci (QTLs). We are continuously updating the resource to further increase its utility for interpreting genetic associations with complex traits. Over the past two years, we have increased the number of uniformly processed studies from 21 to 31 and added X chromosome QTLs for 19 compatible studies. We have also implemented Leafcutter to directly identify splice-junction usage QTLs in all RNA sequencing datasets. Finally, to improve the interpretability of transcript-level QTLs, we have developed static QTL coverage plots that visualise the association between the genotype and average RNA sequencing read coverage in the region for all 1.7 million fine mapped associations. To illustrate the utility of these updates to the eQTL Catalogue, we performed colocalisation analysis between vitamin D levels in the UK Biobank and all molecular QTLs in the eQTL Catalogue. Although most GWAS loci colocalised both with eQTLs and transcript-level QTLs, we found that visual inspection could sometimes be used to distinguish primary splicing QTLs from those that appear to be secondary consequences of large-effect gene expression QTLs. While these visually confirmed primary splicing QTLs explain just 6/53 of the colocalising signals, they are significantly less pleiotropic than eQTLs and identify a prioritised causal gene in 4/6 cases.

Discovery of target genes and pathways at GWAS loci by pooled single-cell CRISPR screens

Most variants associated with complex traits and diseases identified by genome-wide association studies (GWAS) map to noncoding regions of the genome with unknown effects. Using ancestrally diverse, biobank-scale GWAS data, massively parallel CRISPR screens, and single-cell transcriptomic and proteomic sequencing, we discovered 124 cis-target genes of 91 noncoding blood trait GWAS loci. Using precise variant insertion through base editing, we connected specific variants with gene expression changes. We also identified trans-effect networks of noncoding loci when cis target genes encoded transcription factors or microRNAs. Networks were themselves enriched for GWAS variants and demonstrated polygenic contributions to complex traits. This platform enables massively parallel characterization of the target genes and mechanisms of human noncoding variants in both cis and trans.

Genetic dissection of the pluripotent proteome through multi-omics data integration

Genetic background drives phenotypic variability in pluripotent stem cells (PSCs). Most studies to date have used transcript abundance as the primary molecular readout of cell state in PSCs. We performed a comprehensive proteogenomics analysis of 190 genetically diverse mouse embryonic stem cell (mESC) lines. The quantitative proteome is highly variable across lines, and we identified pluripotency-associated pathways that were differentially activated in the proteomics data that were not evident in transcriptome data from the same lines. Integration of protein abundance to transcript levels and chromatin accessibility revealed broad co-variation across molecular layers as well as shared and unique drivers of quantitative variation in pluripotency-associated pathways. Quantitative trait locus (QTL) mapping localized the drivers of these multi-omic signatures to genomic hotspots. This study reveals post-transcriptional mechanisms and genetic interactions that underlie quantitative variability in the pluripotent proteome and provides a regulatory map for mESCs that can provide a basis for future mechanistic studies.

Systematic visualisation of molecular QTLs reveals variant mechanisms at GWAS loci

Splicing quantitative trait loci (QTLs) have been implicated as a common mechanism underlying complex trait associations. However, utilising splicing QTLs in target discovery and prioritisation has been challenging due to extensive data normalisation which often renders the direction of the genetic effect as well as its magnitude difficult to interpret. This is further complicated by the fact that strong expression QTLs often manifest as weak splicing QTLs and vice versa, making it difficult to uniquely identify the underlying molecular mechanism at each locus. We find that these ambiguities can be mitigated by visualising the association between the genotype and average RNA sequencing read coverage in the region. Here, we generate these QTL coverage plots for 1.7 million molecular QTL associations in the eQTL Catalogue identified with five quantification methods. We illustrate the utility of these QTL coverage plots by performing colocalisation between vitamin D levels in the UK Biobank and all molecular QTLs in the eQTL Catalogue. We find that while visually confirmed splicing QTLs explain just 6/53 of the colocalising signals, they are significantly less pleiotropic than eQTLs and identify a prioritised causal gene in 4/6 cases. All our association summary statistics and QTL coverage plots are freely available at https://www.ebi.ac.uk/eqtl/.

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.

A review on qualifications and cost effectiveness of induced pluripotent stem cells (IPSCs)-induced cardiomyocytes in drug screening tests

Human induced pluripotent stem cells (hIPSCs) have initiated a higher degree of successes in disease modelling, preclinical evaluation of drug therapy and pharmaco-toxicological testing. Since the discovery of iPSCs in 2006, many advanced techniques have been introduced to differentiate iPSCs to cardiomyocytes, which have been progressively improved. The disease models from iPSC-induced cardiomyocytes (iPSC-CM) have been successfully helping to study a variety of cardiac diseases such as long QT syndrome, drug-induced long QT, different cardiomyopathies related to mutations in mitochondria or desmosomal proteins and other rare genetic diseases. IPSC-CMs have also been used to screen the role of chemicals in cardiovascular drug discovery and individualisation of drug dosages. In this review, the quality of current procedures for characterisation and maturation of iPSC-CM lines will be discussed. Also, we will focus on time efficiency and cost of standard differentiation methods after reprogramming.

Human In Vitro Models of Epilepsy Using Embryonic and Induced Pluripotent Stem Cells

The challenges in making animal models of complex human epilepsy phenotypes with varied aetiology highlights the need to develop alternative disease models that can address the limitations of animal models by effectively recapitulating human pathophysiology. The advances in stem cell technology provide an opportunity to use human iPSCs to make disease-in-a-dish models. The focus of this review is to report the current information and progress in the generation of epileptic patient-specific iPSCs lines, isogenic control cell lines, and neuronal models. These in vitro models can be used to study the underlying pathological mechanisms of epilepsies, anti-seizure medication resistance, and can also be used for drug testing and drug screening with their isogenic control cell lines.

TERT promoter C228T mutation in neural progenitors confers growth advantage following telomere shortening in vivo

BACKGROUND

Heterozygous TERT (telomerase reverse transcriptase) promoter mutations (TPMs) facilitate TERT expression and are the most frequent mutation in glioblastoma (GBM). A recent analysis revealed this mutation is one of the earliest events in gliomagenesis. However, no appropriate human models have been engineered to study the role of this mutation in the initiation of these tumors.

METHOD

We established GBM models by introducing the heterozygous TPM in human induced pluripotent stem cells (hiPSCs) using a two-step targeting approach in the context of GBM genetic alterations, CDKN2A/B and PTEN deletion, and EGFRvIII overexpression. The impact of the mutation was evaluated through the in vivo passage and in vitro experiment and analysis.

RESULTS

Orthotopic injection of neuronal precursor cells (NPCs) derived from hiPSCs with the TPM into immunodeficient mice did not enhance tumorigenesis compared to TERT promoter wild type NPCs at initial in vivo passage presumably due to relatively long telomeres. However, the mutation recruited GA-Binding Protein and engendered low-level TERT expression resulting in enhanced tumorigenesis and maintenance of short telomeres upon secondary passage as observed in human GBM. These results provide the first insights regarding increased tumorigenesis upon introducing a TPM compared to isogenic controls without TPMs.

CONCLUSION

Our novel GBM models presented the growth advantage of heterozygous TPMs for the first time in the context of GBM driver mutations relative to isogenic controls, thereby allowing for the identification and validation of TERT promoter-specific vulnerabilities in a genetically accurate background.

The consequences of recurrent genetic and epigenetic variants in human pluripotent stem cells

It is well established that human pluripotent stem cells (hPSCs) can acquire genetic and epigenetic changes during culture in vitro. Given the increasing use of hPSCs in research and therapy and the vast expansion in the number of hPSC lines available for researchers, the International Society for Stem Cell Research has recognized the need to reassess quality control standards for ensuring the genetic integrity of hPSCs. Here, we summarize current knowledge of the nature of recurrent genetic and epigenetic variants in hPSC culture, the methods for their detection, and what is known concerning their effects on cell behavior in vitro or in vivo. We argue that the potential consequences of low-level contamination of cell therapy products with cells bearing oncogenic variants are essentially unknown at present. We highlight the key challenges facing the field with particular reference to safety assessment of hPSC-derived cellular therapeutics.

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.

A non-coding GWAS variant impacts anthracycline-induced cardiotoxic phenotypes in human iPSC-derived cardiomyocytes

Anthracyclines, widely used to treat breast cancer, have the potential for cardiotoxicity. We have previously identified and validated a germline single nucleotide polymorphism, rs28714259, associated with an increased risk of anthracycline-induced heart failure. We now provide insights into the mechanism by which rs28714259 might confer increased risk of cardiac damage. Using hiPSC-derived cardiomyocyte cell lines with either intrinsic polymorphism or CRISPR-Cas9-mediated deletion of rs28714259 locus, we demonstrate that glucocorticoid receptor signaling activated by dexamethasone pretreatment prior to doxorubicin exposure preserves cardiomyocyte viability and contractility in cardiomyocytes containing the major allele. Homozygous loss of the rs28714259 major allele diminishes dexamethasone's protective effect. We further demonstrate that the risk allele of rs28714259 disrupts glucocorticoid receptor and rs28714259 binding affinity. Finally, we highlight the activation of genes and pathways involved in cardiac hypertrophy signaling that are blocked by the risk allele, suggesting a decreased adaptive survival response to doxorubicin-related stress.

Modeling the Effects of Maternal Diabetes on the Developing Human Heart Using Pluripotent Stem Cell-Derived Heart Organoids

Congenital heart defects (CHD) constitute the most common type of birth defect in humans. Maternal diabetes during the first trimester of pregnancy (pregestational diabetes, or PGD) is one of the most prominent factors contributing to CHD, and is present in a significant population of female patients with diabetes in reproductive age. PGD is challenging to manage clinically due to the extreme sensitivity of the developing embryo to glucose oscillations, and constitutes a critical health problem for the mother and the fetus. The prevalence of PGD-induced CHD is increasing due to the ongoing diabetes epidemic. While studies using animal models and cells in culture have demonstrated that PGD alters critical cellular and developmental processes, the mechanisms remain obscure, and it is unclear to what extent these models recapitulate PGD-induced CHD in humans. Clinical practice precludes direct studies in developing human embryos, further highlighting the need for physiologically relevant models. To bypass many of these technical and ethical limitations, we describe here a human pluripotent stem cell (hPSC)-based method to generate developmentally relevant self-organizing human heart organoids. By using glucose and insulin to mimic the diabetic environment that the embryo faces in PGD, this system allows modeling critical features of PGD in a human system with relevant physiology, structure, and cell types. The protocol starts with the generation of hPSC-derived embryoid bodies in a 96-well plate, followed by a small molecule-based three-step Wnt activation/inhibition/activation strategy. Organoids are then differentiated under healthy (normal insulin and glucose) and diabetic conditions (high insulin and glucose) over time, allowing for the study of the effects of pregestational diabetes on the developing human heart. We also provide an immunofluorescence protocol for comparing, characterizing, and analyzing the differences between the healthy and diabetic organoids, and comment on additional steps for preparing the organoids for analysis by other techniques after differentiation. © 2022 The Authors. Current Protocols published by Wiley Periodicals LLC. Basic Protocol 1: Generation of hPSC-derived embryoid bodies Basic Protocol 2: Differentiation of EBs into heart organoids under healthy and diabetes-like conditions Basic Protocol 3: Immunofluorescence and organoid preparation for other assays.

Functional genomics and the future of iPSCs in disease modeling

Induced pluripotent stem cells (iPSCs) are valuable in disease modeling because of their potential to expand and differentiate into virtually any cell type and recapitulate key aspects of human biology. Functional genomics are genome-wide studies that aim to discover genotype-phenotype relationships, thereby revealing the impact of human genetic diversity on normal and pathophysiology. In this review, we make the case that human iPSCs (hiPSCs) are a powerful tool for functional genomics, since they provide an in vitro platform for the study of population genetics. We describe cutting-edge tools and strategies now available to researchers, including multi-omics technologies, advances in hiPSC culture techniques, and innovations in drug development. Functional genomics approaches based on hiPSCs hold great promise for advancing drug discovery, disease etiology, and the impact of genetic variation on human biology.

'Channeling' therapeutic discovery for epileptic encephalopathy through iPSC technologies

Induced pluripotent stem cell (iPSC) and gene editing technologies have revolutionized the field of in vitro disease modeling, granting us access to disease-pertinent human cells of the central nervous system. These technologies are particularly well suited for the study of diseases with strong monogenic etiologies. Epilepsy is one of the most common neurological disorders in children, with approximately half of all genetic cases caused by mutations in ion channel genes. These channelopathy-associated epilepsies are clinically diverse, mechanistically complex, and hard to treat. Here, we review the genetic links to epilepsy, the opportunities and challenges of iPSC-based approaches for developing in vitro models of channelopathy-associated disorders, the available tools for effective phenotyping of iPSC-derived neurons, and discuss the potential therapeutic approaches for these devastating diseases.

The Promoting Equity in Stem Cell Genomics Survey

Aim: This study aimed to determine knowledge and attitudes toward induced pluripotent stem cell technology and biobanking. Methods: A survey instrument was developed to determine individuals' knowledge of and attitudes toward these technologies. Results: Results from 276 ethnically diverse participants who took the online survey demonstrated significant associations (p ≤ 0. 05) in knowledge by ethnicity and race regarding properties of stem cells, different types of stem cells and previous sample donation behavior. Significantly more Whites 39% (n = 53) compared with Blacks or African-Americans 19.2% (n = 14) had previous knowledge of induced pluripotent stem cells (χ² = 8.544; p = 0.003) Conclusion: Overall, White race was associated with greater knowledge about stem cells and biobanks and greater willingness to donate samples for future research.

In heart failure reactivation of RNA-binding proteins is associated with the expression of 1,523 fetal-specific isoforms

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.

Heart failure is a chronic condition in which the heart does not pump enough blood. It has been shown that in heart failure, the adult heart reverts to a fetal-like metabolic state and oxygen consumption. Additionally, genes and isoforms that are expressed in the heart only during fetal development (i.e. not in the healthy adult heart) are turned on in heart failure. However, the underlying molecular mechanisms and the extent to which the switch to a fetal gene program occurs remains unclear. In this study, we initially characterized the differences between the fetal and adult heart transcriptomes (entire set of expressed genes and isoforms). We found that RNA binding proteins (RBPs), a family of genes that regulate multiple aspects of a transcript’s maturation, including transcription, splicing and post-transcriptional modifications, play a central role in the differences between fetal and adult heart tissues. We observed that many RBPs that are only expressed in the heart during fetal development become reactivated in heart failure, resulting in the expression of 1,523 fetal-specific isoforms. These findings suggest that reactivation of fetal-specific RBPs in heart failure drives a transcriptome-wide switch to expression of fetal-specific isoforms; and hence that RBPs could potentially serve as novel therapeutic targets.