Mechanosensitive super-enhancers regulate genes linked to atherosclerosis in endothelial cells

Vascular homeostasis and pathophysiology are tightly regulated by mechanical forces generated by hemodynamics. Vascular disorders such as atherosclerotic diseases largely occur at curvatures and bifurcations where disturbed blood flow activates endothelial cells while unidirectional flow at the straight part of vessels promotes endothelial health. Integrated analysis of the endothelial transcriptome, the 3D epigenome, and human genetics systematically identified the SNP-enriched cistrome in vascular endothelium subjected to well-defined atherosclerosis-prone disturbed flow or atherosclerosis-protective unidirectional flow. Our results characterized the endothelial typical- and super-enhancers and underscored the critical regulatory role of flow-sensitive endothelial super-enhancers. CRISPR interference and activation validated the function of a previously unrecognized unidirectional flow-induced super-enhancer that upregulates antioxidant genes NQO1, CYB5B, and WWP2, and a disturbed flow-induced super-enhancer in endothelium which drives prothrombotic genes EDN1 and HIVEP in vascular endothelium. Our results employing multiomics identify the cis-regulatory architecture of the flow-sensitive endothelial epigenome related to atherosclerosis and highlight the regulatory role of super-enhancers in mechanotransduction mechanisms.

Dietary macronutrient composition impacts gene regulation in adipose tissue

Diet is a key lifestyle component that influences metabolic health through several factors, including total energy intake and macronutrient composition. While the impact of caloric intake on gene expression and physiological phenomena in various tissues is well described, the influence of dietary macronutrient composition on these parameters is less well studied. Here, we use the Nutritional Geometry framework to investigate the role of macronutrient composition on metabolic function and gene regulation in adipose tissue. Using ten isocaloric diets that vary systematically in their proportion of energy from fat, protein, and carbohydrates, we find that gene expression and splicing are highly responsive to macronutrient composition, with distinct sets of genes regulated by different macronutrient interactions. Specifically, the expression of many genes associated with Bardet-Biedl syndrome is responsive to dietary fat content. Splicing and expression changes occur in largely separate gene sets, highlighting distinct mechanisms by which dietary composition influences the transcriptome and emphasizing the importance of considering splicing changes to more fully capture the gene regulation response to environmental changes such as diet. Our study provides insight into the gene regulation plasticity of adipose tissue in response to macronutrient composition, beyond the already well-characterized response to caloric intake.

Upregulation of SYNGAP1 expression in mice and human neurons by redirecting alternative splicing

The Ras GTPase-activating protein SYNGAP1 plays a central role in synaptic plasticity, and de novo SYNGAP1 mutations are among the most frequent causes of autism and intellectual disability. How SYNGAP1 is regulated during development and how to treat SYNGAP1-associated haploinsufficiency remain challenging questions. Here, we characterize an alternative 3' splice site (A3SS) of SYNGAP1 that induces nonsense-mediated mRNA decay (A3SS-NMD) in mouse and human neural development. We demonstrate that PTBP1/2 directly bind to and promote SYNGAP1 A3SS inclusion. Genetic deletion of the Syngap1 A3SS in mice upregulates Syngap1 protein and alleviates the long-term potentiation and membrane excitability deficits caused by a Syngap1 knockout allele. We further report a splice-switching oligonucleotide (SSO) that converts SYNGAP1 unproductive isoform to the functional form in human iPSC-derived neurons. This study describes the regulation and function of SYNGAP1 A3SS-NMD, the genetic rescue of heterozygous Syngap1 knockout mice, and the development of an SSO to potentially alleviate SYNGAP1-associated haploinsufficiency.

Single cell profiling at the maternal-fetal interface reveals a deficiency of PD-L1+ non-immune cells in human spontaneous preterm labor

The mechanisms that underlie the timing of labor in humans are largely unknown. In most pregnancies, labor is initiated at term (≥ 37 weeks gestation), but in a signifiicant number of women spontaneous labor occurs preterm and is associated with increased perinatal mortality and morbidity. The objective of this study was to characterize the cells at the maternal-fetal interface (MFI) in term and preterm pregnancies in both the laboring and non-laboring state in Black women, who have among the highest preterm birth rates in the U.S. Using mass cytometry to obtain high-dimensional single-cell resolution, we identified 31 cell populations at the MFI, including 25 immune cell types and six non-immune cell types. Among the immune cells, maternal PD1⁺ CD8 T cell subsets were less abundant in term laboring compared to term non-laboring women. Among the non-immune cells, PD-L1⁺ maternal (stromal) and fetal (extravillous trophoblast) cells were less abundant in preterm laboring compared to term laboring women. Consistent with these observations, the expression of CD274, the gene encoding PD-L1, was significantly depressed and less responsive to fetal signaling molecules in cultured mesenchymal stromal cells from the decidua of preterm compared to term women. Overall, these results suggest that the PD1/PD-L1 pathway at the MFI may perturb the delicate balance between immune tolerance and rejection and contribute to the onset of spontaneous preterm labor.

Male reproductive traits are differentially affected by dietary macronutrient balance but unrelated to adiposity

Dietary factors influence male reproductive function in both experimental and epidemiological studies. However, there are currently no specific dietary guidelines for male preconception health. Here, we use the Nutritional Geometry framework to examine the effects of dietary macronutrient balance on reproductive traits in C57BL/6 J male mice. Dietary effects are observed in a range of morphological, testicular and spermatozoa traits, although the relative influence of protein, fat, carbohydrate, and their interactions differ depending on the trait being examined. Interestingly, dietary fat has a positive influence on sperm motility and antioxidant capacity, differing to typical high fat diet studies where calorie content is not controlled for. Moreover, body adiposity is not significantly correlated with any of the reproductive traits measured in this study. These results demonstrate the importance of macronutrient balance and calorie intake on reproductive function and support the need to develop specific, targeted, preconception dietary guidelines for males.

A non-coding variant linked to metabolic obesity with normal weight affects actin remodelling in subcutaneous adipocytes

Recent large-scale genomic association studies found evidence for a genetic link between increased risk of type 2 diabetes and decreased risk for adiposity-related traits, reminiscent of metabolically obese normal weight (MONW) association signatures. However, the target genes and cellular mechanisms driving such MONW associations remain to be identified. Here, we systematically identify the cellular programmes of one of the top-scoring MONW risk loci, the 2q24.3 risk locus, in subcutaneous adipocytes. We identify a causal genetic variant, rs6712203, an intronic single-nucleotide polymorphism in the COBLL1 gene, which changes the conserved transcription factor motif of POU domain, class 2, transcription factor 2, and leads to differential COBLL1 gene expression by altering the enhancer activity at the locus in subcutaneous adipocytes. We then establish the cellular programme under the genetic control of the 2q24.3 MONW risk locus and the effector gene COBLL1, which is characterized by impaired actin cytoskeleton remodelling in differentiating subcutaneous adipocytes and subsequent failure of these cells to accumulate lipids and develop into metabolically active and insulin-sensitive adipocytes. Finally, we show that perturbations of the effector gene Cobll1 in a mouse model result in organismal phenotypes matching the MONW association signature, including decreased subcutaneous body fat mass and body weight along with impaired glucose tolerance. Taken together, our results provide a mechanistic link between the genetic risk for insulin resistance and low adiposity, providing a potential therapeutic hypothesis and a framework for future identification of causal relationships between genome associations and cellular programmes in other disorders.

African-specific alleles modify risk for asthma at the 17q12-q21 locus in African Americans

BACKGROUND

Asthma is the most common chronic disease in children, occurring at higher frequencies and with more severe disease in children with African ancestry.

METHODS

We tested for association with haplotypes at the most replicated and significant childhood-onset asthma locus at 17q12-q21 and asthma in European American and African American children. Following this, we used whole-genome sequencing data from 1060 African American and 100 European American individuals to identify novel variants on a high-risk African American-specific haplotype. We characterized these variants in silico using gene expression and ATAC-seq data from airway epithelial cells, functional annotations from ENCODE, and promoter capture (pc)Hi-C maps in airway epithelial cells. Candidate causal variants were then assessed for correlation with asthma-associated phenotypes in African American children and adults.

RESULTS

Our studies revealed nine novel African-specific common variants, enriched on a high-risk asthma haplotype, which regulated the expression of GSDMA in airway epithelial cells and were associated with features of severe asthma. Using ENCODE annotations, ATAC-seq, and pcHi-C, we narrowed the associations to two candidate causal variants that are associated with features of T2 low severe asthma.

CONCLUSIONS

Previously unknown genetic variation at the 17q12-21 childhood-onset asthma locus contributes to asthma severity in individuals with African ancestries. We suggest that many other population-specific variants that have not been discovered in GWAS contribute to the genetic risk for asthma and other common diseases.

Abstract P3072: Human Induced Pluripotent Stem Cell Derived Cardiomyocytes (iPSC-CMs) Respond To Inflammatory Stimuli Through Innate Immune Responses

Introduction: Cells respond to inflammatory stimuli in settings of sterile and pathologic cardiac inflammation. Nuclear factor Kappa B (NFKB) is a transcriptional mediator of innate immune responses. Lipopolysaccharide (LPS) promotes NFKB activation through toll-like receptors (TLR) to initiate cell autonomous processes that resist infection. High-mobility box group 1 (HMGB1) is a chromatin associated protein released upon nuclear rupture that similarly promotes cytokine release in response to injury. Few studies have interrogated multiple inflammatory stimuli and their impacts on iPSC-CM function. We examined innate immune activation using distinct TLR agonists and assessed functional responses in both 2D iPSC-CM monoculture and 3D collagen-based engineered heart tissues.

Methods and Results: We measured electric field potential duration (Nanion CardioExcyte 96) of iPSC-CM monolayers after exposure to LPS and HMGB1. We observed a four fold increase in downstroke velocity after exposure to HMGB1 (512.7±126.9 V/s; p=0.0002 n=12) and a 30% decrease in downstroke velocity after LPS treatment (-188.9±61.4 V/s; p=0.004 n=12). This response was associated with nuclear localization of phosphorylated NFKB. Collagen-based engineered heart tissues demonstrated contractile defects in response to 48h of exposure to TLR agonists. A 50% reduction in fractional shortening was observed following HMGB1 treatment (-0.0060±0.0009; p<0.0001 n=6) and an 80% reduction was associated with LPS exposure (-0.0088±0.0006; p<0.0001 n=6). To identify gene expression regulatory regions associated with this response, we examined transcriptional activity of three candidate enhancer regions implicated in cytokine production in iPSC-CMs. These regions demonstrated a 0.5- to 20-fold increase in enhancer activity in iPSC-CMs. One enhancer demonstrated a 22% increase in activity following 48h of LPS, consistent with a direct role for NFKB in innate immune activation in human iPSC-CMs; this enhancer regulates a gene with increased expression in dilated cardiomyopathy.

Conclusion: These data indicate that cardiomyocytes respond directly to inflammatory stimuli and demonstrate the utility of 2D and 3D human tissue culture models for cardiac inflammation.

Genome-wide association and multi-omics studies identify MGMT as a novel risk gene for Alzheimer's disease among women

INTRODUCTION

Variants in the tau gene (MAPT) region are associated with breast cancer in women and Alzheimer's disease (AD) among persons lacking apolipoprotein E ε4 (ε4-).

METHODS

To identify novel genes associated with tau-related pathology, we conducted two genome-wide association studies (GWAS) for AD, one among 10,340 ε4- women in the Alzheimer's Disease Genetics Consortium (ADGC) and another in 31 members (22 women) of a consanguineous Hutterite kindred.

RESULTS

We identified novel associations of AD with MGMT variants in the ADGC (rs12775171, odds ratio [OR] = 1.4, P = 4.9 × 10-8 ) and Hutterite (rs12256016 and rs2803456, OR = 2.0, P = 1.9 × 10-14 ) datasets. Multi-omics analyses showed that the most significant and largest number of associations among the single nucleotide polymorphisms (SNPs), DNA-methylated CpGs, MGMT expression, and AD-related neuropathological traits were observed among women. Furthermore, promoter capture Hi-C analyses revealed long-range interactions of the MGMT promoter with MGMT SNPs and CpG sites.

DISCUSSION

These findings suggest that epigenetically regulated MGMT expression is involved in AD pathogenesis, especially in women.

Common Genetic Variants Contribute to Risk of Transposition of the Great Arteries

RATIONALE

Dextro-transposition of the great arteries (D-TGA) is a severe congenital heart defect which affects approximately 1 in 4,000 live births. While there are several reports of D-TGA patients with rare variants in individual genes, the majority of D-TGA cases remain genetically elusive. Familial recurrence patterns and the observation that most cases with D-TGA are sporadic suggest a polygenic inheritance for the disorder, yet this remains unexplored.

OBJECTIVE

We sought to study the role of common single nucleotide polymorphisms (SNPs) in risk for D-TGA.

METHODS AND RESULTS

We conducted a genome-wide association study in an international set of 1,237 patients with D-TGA and identified a genome-wide significant susceptibility locus on chromosome 3p14.3, which was subsequently replicated in an independent case-control set (rs56219800, meta-analysis P=8.6x10-10, OR=0.69 per C allele). SNP-based heritability analysis showed that 25% of variance in susceptibility to D-TGA may be explained by common variants. A genome-wide polygenic risk score derived from the discovery set was significantly associated to D-TGA in the replication set (P=4x10-5). The genome-wide significant locus (3p14.3) co-localizes with a putative regulatory element that interacts with the promoter of WNT5A, which encodes the Wnt Family Member 5A protein known for its role in cardiac development in mice. We show that this element drives reporter gene activity in the developing heart of mice and zebrafish and is bound by the developmental transcription factor TBX20. We further demonstrate that TBX20 attenuates Wnt5a expression levels in the developing mouse heart.

CONCLUSIONS

This work provides support for a polygenic architecture in D-TGA and identifies a susceptibility locus on chromosome 3p14.3 near WNT5A. Genomic and functional data support a causal role of WNT5A at the locus.

Pluripotent stem cell-derived endometrial stromal fibroblasts in a cyclic, hormone-responsive, coculture model of human decidua

Various human diseases and pregnancy-related disorders reflect endometrial dysfunction. However, rodent models do not share fundamental biological processes with the human endometrium, such as spontaneous decidualization, and no existing human cell cultures recapitulate the cyclic interactions between endometrial stromal and epithelial compartments necessary for decidualization and implantation. Here we report a protocol differentiating human pluripotent stem cells into endometrial stromal fibroblasts (PSC-ESFs) that are highly pure and able to decidualize. Coculture of PSC-ESFs with placenta-derived endometrial epithelial cells generated organoids used to examine stromal-epithelial interactions. Cocultures exhibited specific endometrial markers in the appropriate compartments, organization with cell polarity, and hormone responsiveness of both cell types. Furthermore, cocultures recapitulate a central feature of the human decidua by cyclically responding to hormone withdrawal followed by hormone retreatment. This advance enables mechanistic studies of the cyclic responses that characterize the human endometrium.

Genetic Variation in Enhancers Modifies Cardiomyopathy Gene Expression and Progression

BACKGROUND

Inherited cardiomyopathy associates with a range of phenotypes, mediated by genetic and nongenetic factors. Noninherited cardiomyopathy also displays varying progression and outcomes. Expression of cardiomyopathy genes is under the regulatory control of promoters and enhancers, and human genetic variation in promoters and enhancers may contribute to this variability.

METHODS

We superimposed epigenomic profiling from hearts and cardiomyocytes, including promoter-capture chromatin conformation information, to identify enhancers for 2 cardiomyopathy genes, MYH7 and LMNA. Enhancer function was validated in human cardiomyocytes derived from induced pluripotent stem cells. We also conducted a genome-wide search to ascertain genomic variation in enhancers positioned to alter cardiac expression and correlated one of these variants to cardiomyopathy progression using biobank data.

RESULTS

Multiple enhancers were identified and validated for LMNA and MYH7, including a key enhancer that regulates the switch from MYH6 expression to MYH7 expression. Deletion of this enhancer resulted in a dose-dependent increase in MYH6 and faster contractile rate in engineered heart tissues. We searched for genomic variation in enhancer sequences across the genome, with a focus on nucleotide changes that create or interrupt transcription factor binding sites. The sequence variant, rs875908, disrupts a T-Box Transcription Factor 5 binding motif and maps to an enhancer region 2 kilobases from the transcriptional start site of MYH7. Gene editing to remove the enhancer that harbors this variant markedly reduced MYH7 expression in human cardiomyocytes. Using biobank-derived data, rs875908 associated with longitudinal echocardiographic features of cardiomyopathy.

CONCLUSIONS

Enhancers regulate cardiomyopathy gene expression, and genomic variation within these enhancer regions associates with cardiomyopathic progression over time. This integrated approach identified noncoding modifiers of cardiomyopathy and is applicable to other cardiac genes.

Altered transcriptional and chromatin responses to rhinovirus in bronchial epithelial cells from adults with asthma

There is a life-long relationship between rhinovirus (RV) infection and the development and clinical manifestations of asthma. In this study we demonstrate that cultured primary bronchial epithelial cells from adults with asthma (n = 9) show different transcriptional and chromatin responses to RV infection compared to those without asthma (n = 9). Both the number and magnitude of transcriptional and chromatin responses to RV were muted in cells from asthma cases compared to controls. Pathway analysis of the transcriptionally responsive genes revealed enrichments of apoptotic pathways in controls but inflammatory pathways in asthma cases. Using promoter capture Hi-C we tethered regions of RV-responsive chromatin to RV-responsive genes and showed enrichment of these regions and genes at asthma GWAS loci. Taken together, our studies indicate a delayed or prolonged inflammatory state in cells from asthma cases and highlight genes that may contribute to genetic risk for asthma.

Altered Enhancer and Promoter Usage Leads to Differential Gene Expression in the Normal and Failed Human Heart

BACKGROUND

The failing heart is characterized by changes in gene expression. However, the regulatory regions of the genome that drive these gene expression changes have not been well defined in human hearts.

METHODS

To define genome-wide enhancer and promoter use in heart failure, cap analysis of gene expression sequencing was applied to 3 healthy and 4 failed human hearts to identify promoter and enhancer regions used in left ventricles. Healthy hearts were derived from donors unused for transplantation and failed hearts were obtained as discarded tissue after transplantation.

RESULTS

Cap analysis of gene expression sequencing identified a combined potential for ≈23 000 promoters and ≈5000 enhancers active in human left ventricles. Of these, 17 000 promoters and 1800 enhancers had additional support for their regulatory function. Comparing promoter usage between healthy and failed hearts highlighted promoter shifts which altered aminoterminal protein sequences. Enhancer usage between healthy and failed hearts identified a majority of differentially used heart failure enhancers were intronic and primarily localized within the first intron, revealing this position as a common feature associated with tissue-specific gene expression changes in the heart.

CONCLUSIONS

This data set defines the dynamic genomic regulatory landscape underlying heart failure and serves as an important resource for understanding genetic contributions to cardiac dysfunction. Additionally, regulatory changes contributing to heart failure are attractive therapeutic targets for controlling ventricular remodeling and clinical progression.

Advancing human health in the decade ahead: pregnancy as a key window for discovery: A Burroughs Wellcome Fund Pregnancy Think Tank

Recent revolutionary advances at the intersection of medicine, omics, data sciences, computing, epidemiology, and related technologies inspire us to ponder their impact on health. Their potential impact is particularly germane to the biology of pregnancy and perinatal medicine, where limited improvement in health outcomes for women and children has remained a global challenge. We assembled a group of experts to establish a Pregnancy Think Tank to discuss a broad spectrum of major gestational disorders and adverse pregnancy outcomes that affect maternal-infant lifelong health and should serve as targets for leveraging the many recent advances. This report reflects avenues for future effects that hold great potential in 3 major areas: developmental genomics, including the application of methodologies designed to bridge genotypes, physiology, and diseases, addressing vexing questions in early human development; gestational physiology, from immune tolerance to growth and the timing of parturition; and personalized and population medicine, focusing on amalgamating health record data and deep phenotypes to create broad knowledge that can be integrated into healthcare systems and drive discovery to address pregnancy-related disease and promote general health. We propose a series of questions reflecting development, systems biology, diseases, clinical approaches and tools, and population health, and a call for scientific action. Clearly, transdisciplinary science must advance and accelerate to address adverse pregnancy outcomes. Disciplines not traditionally involved in the reproductive sciences, such as computer science, engineering, mathematics, and pharmacology, should be engaged at the study design phase to optimize the information gathered and to identify and further evaluate potentially actionable therapeutic targets. Information sources should include noninvasive personalized sensors and monitors, alongside instructive "liquid biopsies" for noninvasive pregnancy assessment. Future research should also address the diversity of human cohorts in terms of geography, racial and ethnic distributions, and social and health disparities. Modern technologies, for both data-gathering and data-analyzing, make this possible at a scale that was previously unachievable. Finally, the psychosocial and economic environment in which pregnancy takes place must be considered to promote the health and wellness of communities worldwide.

Abstract MP115: Transcriptional Patterning of the Ventricular Cardiac Conduction System

The heart beats 2 billion times over the average human lifespan and is the manifestation of a pattern of cardiomyocyte depolarization organized by a specialized network of cardiomyocytes, the cardiac conduction system (CCS). Patterning of the CCS into atrial node versus ventricular conduction system (VCS) components with distinct physiology is essential for the normal heartbeat and has been recognized for more than a century. However molecular basis of this regional patterning is not well understood. Using mouse genetics, we found that the ratio between T-box transcriptional activator, Tbx5 , and T-box transcriptional repressor, Tbx3 , determines the molecular and functional output of VCS myocytes. Adult VCS-specific removal of Tbx5 or overexpression of Tbx3 re-patterned the fast VCS into slow, nodal-like cells based on molecular, cellular and functional criteria. Specifically, the ventricular fast conduction gene expression network was lost whereas the slow conduction nodal gene expression network was retained. Action potentials (APs) of Tbx5 -deficient VCS myocytes adopted nodal-specific characteristics, including increased AP duration and cellular automaticity. Removal of Tbx5 in-vivo precipitated inappropriate depolarizations in the His-bundle that initiated lethal ventricular arrhythmias. A T-box rheostat mechanism for CCS patterning was confirmed by Tbx5/Tbx3 genetic interaction studies. TBX5 bound and directly activated cis -regulatory elements at fast conduction loci genome-wide, defining the identity of the adult VCS. Furthermore, TBX5 bound and activated cis-regulatory elements at Tbx5 itself, describing a multi-tiered T-box-dependent fast conduction gene regulatory network (GRN). The hierarchical GRN established a bi-stable network in mathematical modeling, with only high or low fast conduction gene expression states, suggesting a genomic mechanism for the establishment of VCS versus nodal states in-vivo . Thus, the CCS is patterned entirely as a slow, nodal ground state, with a T-box dependent, physiologically dominant, fast conduction network driven specifically in the VCS. Disruption of the fast VCS GRN allowed nodal physiology to emerge, providing a molecular mechanism for some lethal ventricular arrhythmias.

Transcriptional Patterning of the Ventricular Cardiac Conduction System

RATIONALE

The heartbeat is organized by the cardiac conduction system (CCS), a specialized network of cardiomyocytes. Patterning of the CCS into atrial node versus ventricular conduction system (VCS) components with distinct physiology is essential for the normal heartbeat. Distinct node versus VCS physiology has been recognized for more than a century, but the molecular basis of this regional patterning is not well understood.

OBJECTIVE

To study the genetic and genomic mechanisms underlying node versus VCS distinction and investigate rhythm consequences of failed VCS patterning.

METHODS AND RESULTS

Using mouse genetics, we found that the balance between T-box transcriptional activator, Tbx5, and T-box transcriptional repressor, Tbx3, determined the molecular and functional output of VCS myocytes. Adult VCS-specific removal of Tbx5 or overexpression of Tbx3 re-patterned the fast VCS into slow, nodal-like cells based on molecular and functional criteria. In these cases, gene expression profiling showed diminished expression of genes required for VCS-specific fast conduction but maintenance of expression of genes required for nodal slow conduction physiology. Action potentials of Tbx5-deficient VCS myocytes adopted nodal-specific characteristics, including increased action potential duration and cellular automaticity. Removal of Tbx5 in vivo precipitated inappropriate depolarizations in the atrioventricular (His)-bundle associated with lethal ventricular arrhythmias. TBX5 bound and directly activated cis-regulatory elements at fast conduction channel genes required for fast physiological characteristics of the VCS action potential, defining the identity of the adult VCS.

CONCLUSIONS

The CCS is patterned entirely as a slow, nodal ground state, with a T-box dependent, physiologically dominant, fast conduction network driven specifically in the VCS. Disruption of the fast VCS gene regulatory network allowed nodal physiology to emerge, providing a plausible molecular mechanism for some lethal ventricular arrhythmias.

Identification of atrial fibrillation associated genes and functional non-coding variants

Disease-associated genetic variants that lie in non-coding regions found by genome-wide association studies are thought to alter the functionality of transcription regulatory elements and target gene expression. To uncover causal genetic variants, variant regulatory elements and their target genes, here we cross-reference human transcriptomic, epigenomic and chromatin conformation datasets. Of 104 genetic variant regions associated with atrial fibrillation candidate target genes are prioritized. We optimize EMERGE enhancer prediction and use accessible chromatin profiles of human atrial cardiomyocytes to more accurately predict cardiac regulatory elements and identify hundreds of sub-threshold variants that co-localize with regulatory elements. Removal of mouse homologues of atrial fibrillation-associated regions in vivo uncovers a distal regulatory region involved in Gja1 (Cx43) expression. Our analyses provide a shortlist of genes likely affected by atrial fibrillation-associated variants and provide variant regulatory elements in each region that link genetic variation and target gene regulation, helping to focus future investigations.

Tbx20 Is Required in Mid-Gestation Cardiomyocytes and Plays a Central Role in Atrial Development

Supplemental Digital Content is available in the text.

Rationale:

Mutations in the transcription factor TBX20 (T-box 20) are associated with congenital heart disease. Germline ablation of Tbx20 results in abnormal heart development and embryonic lethality by embryonic day 9.5. Because Tbx20 is expressed in multiple cell lineages required for myocardial development, including pharyngeal endoderm, cardiogenic mesoderm, endocardium, and myocardium, the cell type–specific requirement for TBX20 in early myocardial development remains to be explored.

Objective:

Here, we investigated roles of TBX20 in midgestation cardiomyocytes for heart development.

Methods and Results:

Ablation of Tbx20 from developing cardiomyocytes using a doxycycline inducible cTnTCre transgene led to embryonic lethality. The circumference of developing ventricular and atrial chambers, and in particular that of prospective left atrium, was significantly reduced in Tbx20 conditional knockout mutants. Cell cycle analysis demonstrated reduced proliferation of Tbx20 mutant cardiomyocytes and their arrest at the G1-S phase transition. Genome-wide transcriptome analysis of mutant cardiomyocytes revealed differential expression of multiple genes critical for cell cycle regulation. Moreover, atrial and ventricular gene programs seemed to be aberrantly regulated. Putative direct TBX20 targets were identified using TBX20 ChIP-Seq (chromatin immunoprecipitation with high throughput sequencing) from embryonic heart and included key cell cycle genes and atrial and ventricular specific genes. Notably, TBX20 bound a conserved enhancer for a gene key to atrial development and identity, COUP-TFII/Nr2f2 (chicken ovalbumin upstream promoter transcription factor 2/nuclear receptor subfamily 2, group F, member 2). This enhancer interacted with the NR2F2 promoter in human cardiomyocytes and conferred atrial specific gene expression in a transgenic mouse in a TBX20-dependent manner.

Conclusions:

Myocardial TBX20 directly regulates a subset of genes required for fetal cardiomyocyte proliferation, including those required for the G1-S transition. TBX20 also directly downregulates progenitor-specific genes and, in addition to regulating genes that specify chamber versus nonchamber myocardium, directly activates genes required for establishment or maintenance of atrial and ventricular identity. TBX20 plays a previously unappreciated key role in atrial development through direct regulation of an evolutionarily conserved COUPT-FII enhancer.

A novel role of human lung endothelial cells in allergic airway disease by producing and responding to IL-33

Allergic asthma affects more than 300 million people worldwide and is characterized by airway hypersensitivity and eosinophilia. A prevalent feature of allergic asthma is increased serum IL-33 levels, and asthma GWAS have demonstrated that SNPs in the IL-33 locus are significantly associated with disease. While IL-33 and its downstream type 2 responses have been extensively studied in murine models of asthma, much less is known about the regulation of human IL33. Analysis of the mouse and human IL33 loci reveals little genomic conservation between the two species in non-coding regions. We generated a novel BAC transgenic mouse strain containing the human IL-33 locus with a fluorescent reporter to interrogate the regulation and expression of human IL-33. Surprisingly, the mice expressed human IL-33 primarily in endothelial cells, whereas murine IL-33 was expressed primarily in epithelial cells. These results mirror the expression profiles of IL-33 in primary human lung cells from the LungGENS database, thus demonstrating that our novel mouse model faithfully replicates human IL-33 expression. To understand how human IL-33 in the lung is regulated and expressed during inflammation, we examined BAC transgenic mice challenged with either house dust mite extract (HDM) or poly(I:C). In contrast to murine IL-33, expression of human IL-33 was reduced during allergic inflammation. We tested whether IL-33 is able to autoregulate itself via a negative feedback loop. Indeed, IL-33 administration downregulated the expression of human IL-33 in lung endothelial cells. Together, these data emphasize a distinct and novel role in humans for lung endothelial cells in allergic airway disease by producing and responding to IL-33.

DNA methylation in lung cells is associated with asthma endotypes and genetic risk

The epigenome provides a substrate through which environmental exposures can exert their effects on gene expression and disease risk, but the relative importance of epigenetic variation on human disease onset and progression is poorly characterized. Asthma is a heterogeneous disease of the airways, for which both onset and clinical course result from interactions between host genotype and environmental exposures, yet little is known about the molecular mechanisms for these interactions. We assessed genome-wide DNA methylation using the Infinium Human Methylation 450K Bead Chip and characterized the transcriptome by RNA sequencing in primary airway epithelial cells from 74 asthmatic and 41 nonasthmatic adults. Asthma status was based on doctor's diagnosis and current medication use. Genotyping was performed using various Illumina platforms. Our study revealed a regulatory locus on chromosome 17q12-21 associated with asthma risk and epigenetic signatures of specific asthma endotypes and molecular networks. Overall, these data support a central role for DNA methylation in lung cells, which promotes distinct molecular pathways of asthma pathogenesis and modulates the effects of genetic variation on disease risk and clinical heterogeneity.

Probing chromatin landscape reveals roles of endocardial TBX20 in septation

Mutations in the T-box transcription factor TBX20 are associated with multiple forms of congenital heart defects, including cardiac septal abnormalities, but our understanding of the contributions of endocardial TBX20 to heart development remains incomplete. Here, we investigated how TBX20 interacts with endocardial gene networks to drive the mesenchymal and myocardial movements that are essential for outflow tract and atrioventricular septation. Selective ablation of Tbx20 in murine endocardial lineages reduced the expression of extracellular matrix and cell migration genes that are critical for septation. Using the assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq), we identified accessible chromatin within endocardial lineages and intersected these data with TBX20 ChIP-seq and chromatin loop maps to determine that TBX20 binds a conserved long-range enhancer to regulate versican (Vcan) expression. We also observed reduced Vcan expression in Tbx20-deficient mice, supporting a direct role for TBX20 in Vcan regulation. Further, we show that the Vcan enhancer drove reporter gene expression in endocardial lineages in a TBX20-binding site-dependent manner. This work illuminates gene networks that interact with TBX20 to orchestrate cardiac septation and provides insight into the chromatin landscape of endocardial lineages during septation.

HIF1α Represses Cell Stress Pathways to Allow Proliferation of Hypoxic Fetal Cardiomyocytes

Transcriptional mediators of cell stress pathways, including HIF1α, ATF4, and p53, are key to normal development and play critical roles in disease, including ischemia and cancer. Despite their importance, mechanisms by which pathways mediated by these transcription factors interact with one another are not fully understood. In addressing the controversial role of HIF1α in cardiomyocytes (CMs) during heart development, we discovered a mid-gestational requirement for HIF1α for proliferation of hypoxic CMs, involving metabolic switching and a complex interplay among HIF1α, ATF4, and p53. Loss of HIF1α resulted in activation of ATF4 and p53, the latter inhibiting CM proliferation. Bioinformatic and biochemical analyses revealed unexpected mechanisms by which HIF1α intersects with ATF4 and p53 pathways. Our results highlight previously undescribed roles of HIF1α and interactions among major cell stress pathways that could be targeted to enhance proliferation of CMs in ischemia and may have relevance to other diseases, including cancer.

Targeted Germline Modifications in Rats Using CRISPR/Cas9 and Spermatogonial Stem Cells

Organisms with targeted genomic modifications are efficiently produced by gene editing in embryos using CRISPR/Cas9 RNA-guided DNA endonuclease. Here, to facilitate germline editing in rats, we used CRISPR/Cas9 to catalyze targeted genomic mutations in rat spermatogonial stem cell cultures. CRISPR/Cas9-modified spermatogonia regenerated spermatogenesis and displayed long-term sperm-forming potential following transplantation into rat testes. Targeted germline mutations in Epsti1 and Erbb3 were vertically transmitted from recipients to exclusively generate "pure," non-mosaic mutant progeny. Epsti1 mutant rats were produced with or without genetic selection of donor spermatogonia. Monoclonal enrichment of Erbb3 null germlines unmasked recessive spermatogenesis defects in culture that were buffered in recipients, yielding mutant progeny isogenic at targeted alleles. Thus, spermatogonial gene editing with CRISPR/Cas9 provided a platform for generating targeted germline mutations in rats and for studying spermatogenesis.

Early-life physical activity reverses metabolic and Foxo1 epigenetic misregulation induced by gestational sleep disturbance

Sleep disorders are highly prevalent during late pregnancy and can impose adverse effects, such as preeclampsia and diabetes. However, the consequences of sleep fragmentation (SF) on offspring metabolism and epigenomic signatures are unclear. We report that physical activity during early life, but not later, reversed the increased body weight, altered glucose and lipid homeostasis, and increased visceral adipose tissue in offspring of mice subjected to gestational SF (SFo). The reversibility of this phenotype may reflect epigenetic mechanisms induced by SF during gestation. Accordingly, we found that the metabolic master switch Foxo1 was epigenetically misregulated in SFo livers in a temporally regulated fashion. Temporal Foxo1 analysis and its gluconeogenetic targets revealed that the epigenetic abnormalities of Foxo1 precede the metabolic syndrome phenotype. Importantly, regular physical activity early, but not later in life, reversed Foxo1 epigenetic misregulation and altered the metabolic phenotype in gestationally SF-exposed offspring. Thus, we have identified a restricted postnatal period during which lifestyle interventions may reverse the Foxo1 epigenetically mediated risk for metabolic dysfunction later in the life, as induced by gestational sleep disorders.

Evidence of non-pancreatic beta cell-dependent roles of Tcf7l2 in the regulation of glucose metabolism in mice

Non-coding variation within TCF7L2 remains the strongest genetic determinant of type 2 diabetes risk in humans. A considerable effort has been placed in understanding the functional roles of TCF7L2 in pancreatic beta cells, despite evidence of TCF7L2 expression in various peripheral tissues important in glucose homeostasis. Here, we use a humanized mouse model overexpressing Tcf7l2, resulting in glucose intolerance, to infer the contribution of Tcf7l2 overexpression in beta cells and in other tissues to the metabolic phenotypes displayed by these mice. Restoring Tcf7l2 expression specifically in beta cells to endogenous levels, in face of its overexpression elsewhere, results in impaired insulin secretion, reduced beta cell number and islet area, corroborating data obtained in humans showing similar phenotypes as a result of manipulations leading to Tcf7l2 loss of function. Interestingly, the persistent overexpression of Tcf7l2 in non-pancreatic tissues results in a significant worsening in glucose tolerance in vivo, indicating that Tcf7l2 overexpression in beta cells does not account for the glucose intolerance in the Tcf7l2 overexpression mouse model. Collectively, these data posit that Tcf7l2 plays key roles in glucose metabolism through actions beyond pancreatic beta cells, and further points to functionally opposing cell-type specific effects for Tcf7l2 on the maintenance of balanced glucose metabolism, thereby urging a careful examination of its role in non-pancreatic tissues as well as its composite metabolic effects across distinct tissues. Uncovering these roles may lead to new therapeutic targets for type 2 diabetes.

A common genetic variant within SCN10A modulates cardiac SCN5A expression

Variants in SCN10A, which encodes a voltage-gated sodium channel, are associated with alterations of cardiac conduction parameters and the cardiac rhythm disorder Brugada syndrome; however, it is unclear how SCN10A variants promote dysfunctional cardiac conduction. Here we showed by high-resolution 4C-seq analysis of the Scn10a-Scn5a locus in murine heart tissue that a cardiac enhancer located in Scn10a, encompassing SCN10A functional variant rs6801957, interacts with the promoter of Scn5a, a sodium channel-encoding gene that is critical for cardiac conduction. We observed that SCN5A transcript levels were several orders of magnitude higher than SCN10A transcript levels in both adult human and mouse heart tissue. Analysis of BAC transgenic mouse strains harboring an engineered deletion of the enhancer within Scn10a revealed that the enhancer was essential for Scn5a expression in cardiac tissue. Furthermore, the common SCN10A variant rs6801957 modulated Scn5a expression in the heart. In humans, the SCN10A variant rs6801957, which correlated with slowed conduction, was associated with reduced SCN5A expression. These observations establish a genomic mechanism for how a common genetic variation at SCN10A influences cardiac physiology and predisposes to arrhythmia.

Beyond the ENCODE project: using genomics and epigenomics strategies to study enhancer evolution

The complex expression patterns observed for many genes are often regulated by distal transcription enhancers. Changes in the nucleotide sequences of enhancers may therefore lead to changes in gene expression, representing a central mechanism by which organisms evolve. With the development of the experimental technique of chromatin immunoprecipitation (ChIP), in which discrete regions of the genome bound by specific proteins can be identified, it is now possible to identify transcription factor binding events (putative cis-regulatory elements) in entire genomes. Comparing protein-DNA binding maps allows us, for the first time, to attempt to identify regulatory differences and infer global patterns of change in gene expression across species. Here, we review studies that used genome-wide ChIP to study the evolution of enhancers. The trend is one of high divergence of cis-regulatory elements between species, possibly compensated by extensive creation and loss of regulatory elements and rewiring of their target genes. We speculate on the meaning of the differences observed and discuss that although ChIP experiments identify the biochemical event of protein-DNA interaction, it cannot determine whether the event results in a biological function, and therefore more studies are required to establish the effect of divergence of binding events on species-specific gene expression.

Enhancers: five essential questions

It is estimated that the human genome contains hundreds of thousands of enhancers, so understanding these gene-regulatory elements is a crucial goal. Several fundamental questions need to be addressed about enhancers, such as how do we identify them all, how do they work, and how do they contribute to disease and evolution? Five prominent researchers in this field look at how much we know already and what needs to be done to answer these questions.

Regulatory variation in a TBX5 enhancer leads to isolated congenital heart disease

Recent studies have identified the genetic underpinnings of a growing number of diseases through targeted exome sequencing. However, this strategy ignores the large component of the genome that does not code for proteins, but is nonetheless biologically functional. To address the possible involvement of regulatory variation in congenital heart diseases (CHDs), we searched for regulatory mutations impacting the activity of TBX5, a dosage-dependent transcription factor with well-defined roles in the heart and limb development that has been associated with the Holt-Oram syndrome (heart-hand syndrome), a condition that affects 1/100 000 newborns. Using a combination of genomics, bioinformatics and mouse genetic engineering, we scanned ∼700 kb of the TBX5 locus in search of cis-regulatory elements. We uncovered three enhancers that collectively recapitulate the endogenous expression pattern of TBX5 in the developing heart. We re-sequenced these enhancer elements in a cohort of non-syndromic patients with isolated atrial and/or ventricular septal defects, the predominant cardiac defects of the Holt-Oram syndrome, and identified a patient with a homozygous mutation in an enhancer ∼90 kb downstream of TBX5. Notably, we demonstrate that this single-base-pair mutation abrogates the ability of the enhancer to drive expression within the heart in vivo using both mouse and zebrafish transgenic models. Given the population-wide frequency of this variant, we estimate that 1/100 000 individuals would be homozygous for this variant, highlighting that a significant number of CHD associated with TBX5 dysfunction might arise from non-coding mutations in TBX5 heart enhancers, effectively decoupling the heart and hand phenotypes of the Holt-Oram syndrome.

TBX5 drives Scn5a expression to regulate cardiac conduction system function

Cardiac conduction system (CCS) disease, which results in disrupted conduction and impaired cardiac rhythm, is common with significant morbidity and mortality. Current treatment options are limited, and rational efforts to develop cell-based and regenerative therapies require knowledge of the molecular networks that establish and maintain CCS function. Recent genome-wide association studies (GWAS) have identified numerous loci associated with adult human CCS function, including TBX5 and SCN5A. We hypothesized that TBX5, a critical developmental transcription factor, regulates transcriptional networks required for mature CCS function. We found that deletion of Tbx5 from the mature murine ventricular conduction system (VCS), including the AV bundle and bundle branches, resulted in severe VCS functional consequences, including loss of fast conduction, arrhythmias, and sudden death. Ventricular contractile function and the VCS fate map remained unchanged in VCS-specific Tbx5 knockouts. However, key mediators of fast conduction, including Nav1.5, which is encoded by Scn5a, and connexin 40 (Cx40), demonstrated Tbx5-dependent expression in the VCS. We identified a TBX5-responsive enhancer downstream of Scn5a sufficient to drive VCS expression in vivo, dependent on canonical T-box binding sites. Our results establish a direct molecular link between Tbx5 and Scn5a and elucidate a hierarchy between human GWAS loci that affects function of the mature VCS, establishing a paradigm for understanding the molecular pathology of CCS disease.

Glyoxalase 1 increases anxiety by reducing GABAA receptor agonist methylglyoxal

Glyoxalase 1 (Glo1) expression has previously been associated with anxiety in mice; however, its role in anxiety is controversial, and the underlying mechanism is unknown. Here, we demonstrate that GLO1 increases anxiety by reducing levels of methylglyoxal (MG), a GABAA receptor agonist. Mice overexpressing Glo1 on a Tg bacterial artificial chromosome displayed increased anxiety-like behavior and reduced brain MG concentrations. Treatment with low doses of MG reduced anxiety-like behavior, while higher doses caused locomotor depression, ataxia, and hypothermia, which are characteristic effects of GABAA receptor activation. Consistent with these data, we found that physiological concentrations of MG selectively activated GABAA receptors in primary neurons. These data indicate that GLO1 increases anxiety by reducing levels of MG, thereby decreasing GABAA receptor activation. More broadly, our findings potentially link metabolic state, neuronal inhibitory tone, and behavior. Finally, we demonstrated that pharmacological inhibition of GLO1 reduced anxiety, suggesting that GLO1 is a possible target for the treatment of anxiety disorders.

An in vivo cis-regulatory screen at the type 2 diabetes associated TCF7L2 locus identifies multiple tissue-specific enhancers

Genome-wide association studies (GWAS) have repeatedly shown an association between non-coding variants in the TCF7L2 locus and risk for type 2 diabetes (T2D), implicating a role for cis-regulatory variation within this locus in disease etiology. Supporting this hypothesis, we previously localized complex regulatory activity to the TCF7L2 T2D-associated interval using an in vivo bacterial artificial chromosome (BAC) enhancer-trapping reporter strategy. To follow-up on this broad initial survey of the TCF7L2 regulatory landscape, we performed a fine-mapping enhancer scan using in vivo mouse transgenic reporter assays. We functionally interrogated approximately 50% of the sequences within the T2D-associated interval, utilizing sequence conservation within this 92-kb interval to determine the regulatory potential of all evolutionary conserved sequences that exhibited conservation to the non-eutherian mammal opossum. Included in this study was a detailed functional interrogation of sequences spanning both protective and risk alleles of single nucleotide polymorphism (SNP) rs7903146, which has exhibited allele-specific enhancer function in pancreatic beta cells. Using these assays, we identified nine segments regulating various aspects of the TCF7L2 expression profile and that constitute nearly 70% of the sequences tested. These results highlight the regulatory complexity of this interval and support the notion that a TCF7L2 cis-regulatory disruption leads to T2D predisposition.

A fine-scale chimpanzee genetic map from population sequencing

To study the evolution of recombination rates in apes, we developed methodology to construct a fine-scale genetic map from high-throughput sequence data from 10 Western chimpanzees, Pan troglodytes verus. Compared to the human genetic map, broad-scale recombination rates tend to be conserved, but with exceptions, particularly in regions of chromosomal rearrangements and around the site of ancestral fusion in human chromosome 2. At fine scales, chimpanzee recombination is dominated by hotspots, which show no overlap with those of humans even though rates are similarly elevated around CpG islands and decreased within genes. The hotspot-specifying protein PRDM9 shows extensive variation among Western chimpanzees, and there is little evidence that any sequence motifs are enriched in hotspots. The contrasting locations of hotspots provide a natural experiment, which demonstrates the impact of recombination on base composition.

Dual transcriptional activator and repressor roles of TBX20 regulate adult cardiac structure and function

The ongoing requirement in adult heart for transcription factors with key roles in cardiac development is not well understood. We recently demonstrated that TBX20, a transcriptional regulator required for cardiac development, has key roles in the maintenance of functional and structural phenotypes in adult mouse heart. Conditional ablation of Tbx20 in adult cardiomyocytes leads to a rapid onset and progression of heart failure, with prominent conduction and contractility phenotypes that lead to death. Here we describe a more comprehensive molecular characterization of the functions of TBX20 in adult mouse heart. Coupling genome-wide chromatin immunoprecipitation and transcriptome analyses (RNA-Seq), we identified a subset of genes that change expression in Tbx20 adult cardiomyocyte-specific knockout hearts which are direct downstream targets of TBX20. This analysis revealed a dual role for TBX20 as both a transcriptional activator and a repressor, and that each of these functions regulates genes with very specialized and distinct molecular roles. We also show how TBX20 binds to its targets genome-wide in a context-dependent manner, using various cohorts of co-factors to either promote or repress distinct genetic programs within adult heart. Our integrative approach has uncovered several novel aspects of TBX20 and T-box protein function within adult heart. Sequencing data accession number (http://www.ncbi.nlm.nih.gov/geo): GSE30943.

Transcriptional enhancers in development and disease

Distal transcription enhancers are cis-regulatory elements that promote gene expression, enabling spatiotemporal control of genetic programs such as those required in metazoan developmental processes. Because of their importance, their disruption can lead to disease.

Tbx20 regulates a genetic program essential to adult mouse cardiomyocyte function

Human mutations in or variants of TBX20 are associated with congenital heart disease, cardiomyopathy, and arrhythmias. To investigate whether cardiac disease in patients with these conditions results from an embryonic or ongoing requirement for Tbx20 in myocardium, we ablated Tbx20 specifically in adult cardiomyocytes in mice. This ablation resulted in the onset of severe cardiomyopathy accompanied by arrhythmias, with death ensuing within 1 to 2 weeks of Tbx20 ablation. Accounting for this dramatic phenotype, we identified molecular signatures that posit Tbx20 as a central integrator of a genetic program that maintains cardiomyocyte function in the adult heart. Expression of a number of genes encoding critical transcription factors, ion channels, and cytoskeletal/myofibrillar proteins was downregulated consequent to loss of Tbx20. Genome-wide ChIP analysis of Tbx20-binding regions in the adult heart revealed that many of these genes were direct downstream targets of Tbx20 and uncovered a previously undescribed DNA-binding site for Tbx20. Bioinformatics and in vivo functional analyses revealed a cohort of transcription factors that, working with Tbx20, integrated multiple environmental signals to maintain ion channel gene expression in the adult heart. Our data provide insight into the mechanisms by which mutations in TBX20 cause adult heart disease in humans.

Modulation ofTcf7l2 expression alters behavior in mice

The comorbidity of type 2 diabetes (T2D) with several psychiatric diseases is well established. While environmental factors may partially account for these co-occurrences, common genetic susceptibilities could also be implicated in the confluence of these diseases. In support of shared genetic burdens, TCF7L2, the strongest genetic determinant for T2D risk in the human population, has been recently implicated in schizophrenia (SCZ) risk, suggesting that this may be one of many loci that pleiotropically influence both diseases. To investigate whether Tcf7l2 is involved in behavioral phenotypes in addition to its roles in glucose metabolism, we conducted several behavioral tests in mice with null alleles of Tcf7l2 or overexpressing Tcf7l2. We identified a role for Tcf7l2 in anxiety-like behavior and a dose-dependent effect of Tcf7l2 alleles on fear learning. None of the mutant mice showed differences in prepulse inhibition (PPI), which is a well-established endophenotype for SCZ. These results show that Tcf7l2 alters behavior in mice. Importantly, these differences are observed prior to the onset of detectable glucose metabolism abnormalities. Whether these differences are related to human anxiety-disorders or schizophrenia remains to be determined. These animal models have the potential to elucidate the molecular basis of psychiatric comorbidities in diabetes and should therefore be studied further.

Genome-wide identification of conserved regulatory function in diverged sequences

Plasticity of gene regulatory encryption can permit DNA sequence divergence without loss of function. Functional information is preserved through conservation of the composition of transcription factor binding sites (TFBS) in a regulatory element. We have developed a method that can accurately identify pairs of functional noncoding orthologs at evolutionarily diverged loci by searching for conserved TFBS arrangements. With an estimated 5% false-positive rate (FPR) in approximately 3000 human and zebrafish syntenic loci, we detected approximately 300 pairs of diverged elements that are likely to share common ancestry and have similar regulatory activity. By analyzing a pool of experimentally validated human enhancers, we demonstrated that 7/8 (88%) of their predicted functional orthologs retained in vivo regulatory control. Moreover, in 5/7 (71%) of assayed enhancer pairs, we observed concordant expression patterns. We argue that TFBS composition is often necessary to retain and sufficient to predict regulatory function in the absence of overt sequence conservation, revealing an entire class of functionally conserved, evolutionarily diverged regulatory elements that we term "covert."

Alterations in TCF7L2 expression define its role as a key regulator of glucose metabolism

Genome-wide association studies (GWAS) have consistently implicated noncoding variation within the TCF7L2 locus with type 2 diabetes (T2D) risk. While this locus represents the strongest genetic determinant for T2D risk in humans, it remains unclear how these noncoding variants affect disease etiology. To test the hypothesis that the T2D-associated interval harbors cis-regulatory elements controlling TCF7L2 expression, we conducted in vivo transgenic reporter assays to characterize the TCF7L2 regulatory landscape. We found that the 92-kb genomic interval associated with T2D harbors long-range enhancers regulating various aspects of the spatial-temporal expression patterns of TCF7L2, including expression in tissues involved in the control of glucose homeostasis. By selectively deleting this interval, we establish a critical role for these enhancers in robust TCF7L2 expression. To further determine whether variation in Tcf7l2 expression may lead to diabetes, we developed a Tcf7l2 copy-number allelic series in mice. We show that a null Tcf7l2 allele leads, in a dose-dependent manner, to lower glycemic profiles. Tcf7l2 null mice also display enhanced glucose tolerance coupled to significantly lowered insulin levels, suggesting that these mice are protected against T2D. Confirming these observations, transgenic mice harboring multiple Tcf7l2 copies and overexpressing this gene display reciprocal phenotypes, including glucose intolerance. These results directly demonstrate that Tcf7l2 plays a role in regulating glucose tolerance, suggesting that overexpression of this gene is associated with increased risk of T2D. These data highlight the role of enhancer elements as mediators of T2D risk in humans, strengthening the evidence that variation in cis-regulatory elements may be a paradigm for genetic predispositions to common disease.

The emerging genetic landscape underlying cardiac conduction system function

Proper function of an organized Cardiac Conduction System (CCS) is vital to the survival of metazoans ranging from fly to man. The routine use of non-invasive electrocardiogram measures in the diagnosis and monitoring of cardiovascular health has established a trove of reliable CCS functional data in both normal and diseased cardiac states. Recent combination of echocardiogram (ECG) data with genome-wide association studies has identified genomic regions implicated in ECG variability which impact CCS function. In this study, we review the substantial recent progress in this area, highlighting the identification of novel loci, confirming the importance of previously implicated loci in CCS function, and exploring potential links between genes with important roles in developmental processes and variation in function of the CCS.

Genome-wide maps of transcription regulatory elements

Expression of eukaryotic genes with complex spatial-temporal regulation during development requires finer regulation than that of genes with simpler expression patterns. Given the high degree of conservation of the developmental gene set across distantly related phylogenetic taxa, it is argued that evolutionary variation has occurred by tweaking regulation of expression of developmental genes, rather than by changes in genes themselves. Complex regulation is often achieved through the coordinated action of transcription regulatory elements spread across the genome up to tens of kilobases from the promoters of their target genes. Disruption of regulatory elements has been implicated in several diseases and studies showing associations between disease traits and nonprotein coding variation hint for a role of regulatory elements as cause of diseases. Therefore, the identification and mapping of regulatory elements in genome scale is crucial to understand how gene expression is regulated, how organisms evolve, and to identify sequence variation causing diseases. Previously developed experimental techniques have been adapted to identify regulatory elements in genome scale and high-throughput, allowing a global view of their biological roles. We review methods as chromatin immunoprecipitation, DNase I hypersensitivity, and computational approaches and how they have been employed to generate maps of histone modifications, open chromatin, nucleosome positioning, and transcription factor binding regions in whole mammalian genomes. Given the importance of non-promoter elements in gene regulation and the recent explosion in the number of studies devoted to them, we focus on these elements and discuss the insights on gene regulation being obtained by these studies.

An 8q24 gene desert variant associated with prostate cancer risk confers differential in vivo activity to a MYC enhancer

Genome-wide association studies (GWAS) routinely identify risk variants in noncoding DNA, as exemplified by reports of multiple single nucleotide polymorphisms (SNPs) associated with prostate cancer in five independent regions in a gene desert on 8q24. Two of these regions also have been associated with breast and colorectal cancer. These findings implicate functional variation within long-range cis-regulatory elements in disease etiology. We used an in vivo bacterial artificial chromosome (BAC) enhancer-trapping strategy in mice to scan a half-megabase of the 8q24 gene desert encompassing the prostate cancer-associated regions for long-range cis-regulatory elements. These BAC assays identified both prostate and mammary gland enhancer activities within the region. We demonstrate that the 8q24 cancer-associated variant rs6983267 lies within an in vivo prostate enhancer whose expression mimics that of the nearby MYC proto-oncogene. Additionally, we show that the cancer risk allele increases prostate enhancer activity in vivo relative to the non-risk allele. This allele-specific enhancer activity is detectable during early prostate development and throughout prostate maturation, raising the possibility that this SNP could assert its influence on prostate cancer risk before tumorigenesis occurs. Our study represents an efficient strategy to build experimentally on GWAS findings with an in vivo method for rapidly scanning large regions of noncoding DNA for functional cis-regulatory sequences harboring variation implicated in complex diseases.

Homotypic clusters of transcription factor binding sites are a key component of human promoters and enhancers

Clustering of multiple transcription factor binding sites (TFBSs) for the same transcription factor (TF) is a common feature of cis-regulatory modules in invertebrate animals, but the occurrence of such homotypic clusters of TFBSs (HCTs) in the human genome has remained largely unknown. To explore whether HCTs are also common in human and other vertebrates, we used known binding motifs for vertebrate TFs and a hidden Markov model-based approach to detect HCTs in the human, mouse, chicken, and fugu genomes, and examined their association with cis-regulatory modules. We found that evolutionarily conserved HCTs occupy nearly 2% of the human genome, with experimental evidence for individual TFs supporting their binding to predicted HCTs. More than half of the promoters of human genes contain HCTs, with a distribution around the transcription start site in agreement with the experimental data from the ENCODE project. In addition, almost half of the 487 experimentally validated developmental enhancers contain them as well--a number more than 25-fold larger than expected by chance. We also found evidence of negative selection acting on TFBSs within HCTs, as the conservation of TFBSs is stronger than the conservation of sequences separating them. The important role of HCTs as components of developmental enhancers is additionally supported by a strong correlation between HCTs and the binding of the enhancer-associated coactivator protein Ep300 (also known as p300). Experimental validation of HCT-containing elements in both zebrafish and mouse suggest that HCTs could be used to predict both the presence of enhancers and their tissue specificity, and are thus a feature that can be effectively used in deciphering the gene regulatory code. In conclusion, our results indicate that HCTs are a pervasive feature of human cis-regulatory modules and suggest that they play an important role in gene regulation in the human and other vertebrate genomes.

Genome-wide discovery of human heart enhancers

The various organogenic programs deployed during embryonic development rely on the precise expression of a multitude of genes in time and space. Identifying the cis-regulatory elements responsible for this tightly orchestrated regulation of gene expression is an essential step in understanding the genetic pathways involved in development. We describe a strategy to systematically identify tissue-specific cis-regulatory elements that share combinations of sequence motifs. Using heart development as an experimental framework, we employed a combination of Gibbs sampling and linear regression to build a classifier that identifies heart enhancers based on the presence and/or absence of various sequence features, including known and putative transcription factor (TF) binding specificities. In distinguishing heart enhancers from a large pool of random noncoding sequences, the performance of our classifier is vastly superior to four commonly used methods, with an accuracy reaching 92% in cross-validation. Furthermore, most of the binding specificities learned by our method resemble the specificities of TFs widely recognized as key players in heart development and differentiation, such as SRF, MEF2, ETS1, SMAD, and GATA. Using our classifier as a predictor, a genome-wide scan identified over 40,000 novel human heart enhancers. Although the classifier used no gene expression information, these novel enhancers are strongly associated with genes expressed in the heart. Finally, in vivo tests of our predictions in mouse and zebrafish achieved a validation rate of 62%, significantly higher than what is expected by chance. These results support the existence of underlying cis-regulatory codes dictating tissue-specific transcription in mammalian genomes and validate our enhancer classifier strategy as a method to uncover these regulatory codes.

Distinct genetic regulation of progression of diabetes and renal disease in the Goto-Kakizaki rat

Goto-Kakizaki (GK) rats develop early-onset type 2 diabetes (T2D) symptoms, with signs of diabetic nephropathy becoming apparent with aging. To determine whether T2D and renal disease share similar genetic architecture, we ran a quantitative trait locus (QTL) analysis in the F2 progeny of a GK x Brown Norway (BN) rat cross. Further, to determine whether genetic components change over time, we ran the QTL analysis on phenotypes collected longitudinally, at 3, 6, 9 and 12 mo, from the same animals. We confirmed three chromosomal regions that are linked to early diabetes phenotypes (chromosomes 1, 5, and 10) and a single region involved in the late progression of the disorder (chromosome 4). A single region was identified for the onset of the renal phenotype proteinuria (chromosome 5). This region overlaps the diabetic QTL, although it is not certain whether similar genes are involved in both phenotypes. A second QTL linked to the progression of the renal phenotype was found on chromosome 7. Linkage for triglyceride and cholesterol levels were also identified (chromosomes 7 and 8, respectively). These results demonstrate that, in general, different genetic components control diabetic and renal phenotypes in a diabetic nephropathy model. Furthermore, these results demonstrate that, over time, different genetic components are involved in progression of disease from those that were involved in disease onset. This observation would suggest that clinical studies collecting participants over a wide age distribution may be diluting genetic effects and reducing power to detect true effects.

Further delineation of the 15q13 microdeletion and duplication syndromes: a clinical spectrum varying from non-pathogenic to a severe outcome

BACKGROUND

Recurrent 15q13.3 microdeletions were recently identified with identical proximal (BP4) and distal (BP5) breakpoints and associated with mild to moderate mental retardation and epilepsy.

METHODS

To assess further the clinical implications of this novel 15q13.3 microdeletion syndrome, 18 new probands with a deletion were molecularly and clinically characterised. In addition, we evaluated the characteristics of a family with a more proximal deletion between BP3 and BP4. Finally, four patients with a duplication in the BP3-BP4-BP5 region were included in this study to ascertain the clinical significance of duplications in this region.

RESULTS

The 15q13.3 microdeletion in our series was associated with a highly variable intra- and inter-familial phenotype. At least 11 of the 18 deletions identified were inherited. Moreover, 7 of 10 siblings from four different families also had this deletion: one had a mild developmental delay, four had only learning problems during childhood, but functioned well in daily life as adults, whereas the other two had no learning problems at all. In contrast to previous findings, seizures were not a common feature in our series (only 2 of 17 living probands). Three patients with deletions had cardiac defects and deletion of the KLF13 gene, located in the critical region, may contribute to these abnormalities. The limited data from the single family with the more proximal BP3-BP4 deletion suggest this deletion may have little clinical significance. Patients with duplications of the BP3-BP4-BP5 region did not share a recognisable phenotype, but psychiatric disease was noted in 2 of 4 patients.

CONCLUSIONS

Overall, our findings broaden the phenotypic spectrum associated with 15q13.3 deletions and suggest that, in some individuals, deletion of 15q13.3 is not sufficient to cause disease. The existence of microdeletion syndromes, associated with an unpredictable and variable phenotypic outcome, will pose the clinician with diagnostic difficulties and challenge the commonly used paradigm in the diagnostic setting that aberrations inherited from a phenotypically normal parent are usually without clinical consequences.