Paternal dietary macronutrient balance and energy intake drive metabolic and behavioral differences among offspring.

Paternal diet can influence the phenotype of the next generation, yet, the dietary components inducing specific responses in the offspring are not identified. Here, we use the Nutritional Geometry Framework to determine the effects of pre-conception paternal dietary macronutrient balance on offspring metabolic and behavioral traits in mice. Ten isocaloric diets varying in the relative proportion of protein, fats, and carbohydrates are fed to male mice prior to mating. Dams and offspring are fed standard chow and never exposed to treatment diets. Body fat in female offspring is positively associated with the paternal consumption of fat, while in male offspring, an anxiety-like phenotype is associated to paternal diets low in protein and high in carbohydrates. Our study uncovers that the nature and the magnitude of paternal effects are driven by interactions between macronutrient balance and energy intake and are not solely the result of over- or undernutrition.

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.

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.

Splicing across adipocyte differentiation is highly dynamic and impacted by metabolic phenotype.

Adipose tissue dysfunction underlies many of the metabolic complications associated with obesity. A better understanding of the gene regulation differences present in metabolically unhealthy adipose tissue can provide insights into the mechanisms underlying adipose tissue dysfunction. Here, we used RNA-seq data collected from a differentiation time course of lean, obese, and obese with type 2 diabetes (T2D) individuals to characterize the role of alterative splicing in adipocyte differentiation and function. We found that splicing was highly dynamic across adipocyte differentiation in all three cohorts, and that the dynamics of splicing were significantly impacted by metabolic phenotype. We also found that there was very little overlap between genes that were differentially spliced in adipocyte differentiation and those that were differentially expressed, positioning alternative splicing as a largely independent gene regulatory mechanism whose impact would be missed when looking at gene expression changes alone. To assess the impact of alternative splicing across adipocyte differentiation on genetic risk for metabolic diseases, we integrated the differential splicing results generated here with genome-wide association study results for body mass index and T2D, and found that variants associated with T2D were enriched in regions that were differentially spliced in early differentiation. These findings provide insight into the role of alternative splicing in adipocyte differentiation and can serve as a resource to guide future variant-to-function studies.

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.

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.

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.

Genetics of sexually dimorphic adipose distribution in humans.

Obesity-associated morbidity is exacerbated by abdominal obesity, which can be measured as the waist-to-hip ratio adjusted for the body mass index (WHRadjBMI). Here we identify genes associated with obesity and WHRadjBMI and characterize allele-sensitive enhancers that are predicted to regulate WHRadjBMI genes in women. We found that several waist-to-hip ratio-associated variants map within primate-specific Alu retrotransposons harboring a DNA motif associated with adipocyte differentiation. This suggests that a genetic component of adipose distribution in humans may involve co-option of retrotransposons as adipose enhancers. We evaluated the role of the strongest female WHRadjBMI-associated gene, SNX10, in adipose biology. We determined that it is required for human adipocyte differentiation and function and participates in diet-induced adipose expansion in female mice, but not males. Our data identify genes and regulatory mechanisms that underlie female-specific adipose distribution and mediate metabolic dysfunction in women.

A pleiotropic hypoxia-sensitive EPAS1 enhancer is disrupted by adaptive alleles in Tibetans.

In Tibetans, noncoding alleles in EPAS1-whose protein product hypoxia-inducible factor 2α (HIF-2α) drives the response to hypoxia-carry strong signatures of positive selection; however, their functional mechanism has not been systematically examined. Here, we report that high-altitude alleles disrupt the activity of four EPAS1 enhancers in one or more cell types. We further characterize one enhancer (ENH5) whose activity is both allele specific and hypoxia dependent. Deletion of ENH5 results in down-regulation of EPAS1 and HIF-2α targets in acute hypoxia and in a blunting of the transcriptional response to sustained hypoxia. Deletion of ENH5 in mice results in dysregulation of gene expression across multiple tissues. We propose that pleiotropic adaptive effects of the Tibetan alleles in EPAS1 underlie the strong selective signal at this gene.

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.

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.

Asthma-associated genetic variants induce IL33 differential expression through an enhancer-blocking regulatory region.

Genome-wide association studies (GWAS) have implicated the IL33 locus in asthma, but the underlying mechanisms remain unclear. Here, we identify a 5 kb region within the GWAS-defined segment that acts as an enhancer-blocking element in vivo and in vitro. Chromatin conformation capture showed that this 5 kb region loops to the IL33 promoter, potentially regulating its expression. We show that the asthma-associated single nucleotide polymorphism (SNP) rs1888909, located within the 5 kb region, is associated with IL33 gene expression in human airway epithelial cells and IL-33 protein expression in human plasma, potentially through differential binding of OCT-1 (POU2F1) to the asthma-risk allele. Our data demonstrate that asthma-associated variants at the IL33 locus mediate allele-specific regulatory activity and IL33 expression, providing a mechanism through which a regulatory SNP contributes to genetic risk of asthma.

A functional genomics pipeline identifies pleiotropy and cross-tissue effects within obesity-associated GWAS loci.

Genome-wide association studies (GWAS) have identified many disease-associated variants, yet mechanisms underlying these associations remain unclear. To understand obesity-associated variants, we generate gene regulatory annotations in adipocytes and hypothalamic neurons across cellular differentiation stages. We then test variants in 97 obesity-associated loci using a massively parallel reporter assay and identify putatively causal variants that display cell type specific or cross-tissue enhancer-modulating properties. Integrating these variants with gene regulatory information suggests genes that underlie obesity GWAS associations. We also investigate a complex genomic interval on 16p11.2 where two independent loci exhibit megabase-range, cross-locus chromatin interactions. We demonstrate that variants within these two loci regulate a shared gene set. Together, our data support a model where GWAS loci contain variants that alter enhancer activity across tissues, potentially with temporally restricted effects, to impact the expression of multiple genes. This complex model has broad implications for ongoing efforts to understand GWAS.

Establishment of human induced trophoblast stem-like cells from term villous cytotrophoblasts.

Human trophoblast stem cells (hTSC) can be isolated from first trimester placenta but not from term placenta. Here we demonstrate that villous cytotrophoblasts (vCTB) from term placenta can be reprogrammed into induced trophoblastic stem-like cells (iTSC) by introducing sets of transcription factors. The iTSCs express TSC markers such as GATA3, TEAD4 and ELF5, and are multipotent, validated by their differentiation into both extravillous trophoblasts (EVT) and syncytiotrophoblasts (STB) in vitro and in vivo. The iTSC can be passaged indefinitely in vitro without slowing of growth. The transcriptome profile of these cells closely resembles the profile of hTSC isolated from first trimester placentae but different from the term placental vCTB from which they originated. The ability to reprogram cells from term placenta into iTSC will allow study of early gestation events which impact placental function later in gestation, including preeclampsia and spontaneous preterm birth.

Extensive pleiotropism and allelic heterogeneity mediate metabolic effects of IRX3 and IRX5.

Whereas coding variants often have pleiotropic effects across multiple tissues, noncoding variants are thought to mediate their phenotypic effects by specific tissue and temporal regulation of gene expression. Here, we investigated the genetic and functional architecture of a genomic region within the FTO gene that is strongly associated with obesity risk. We show that multiple variants on a common haplotype modify the regulatory properties of several enhancers targeting IRX3 and IRX5 from megabase distances. We demonstrate that these enhancers affect gene expression in multiple tissues, including adipose and brain, and impart regulatory effects during a restricted temporal window. Our data indicate that the genetic architecture of disease-associated loci may involve extensive pleiotropy, allelic heterogeneity, shared allelic effects across tissues, and temporally restricted effects.

Genetic variants in LGALS3 are related to lower galectin-3 serum levels and clinical outcomes in systemic sclerosis patients: A case-control study.

INTRODUCTION: Systemic sclerosis (SSc) is a rare complex disease characterized by vascular damage, autoimmunity, and extensive skin and internal organs fibrosis. Galectin-3 (Gal-3) is encoded by gene LGALS3 (Lectin, Galactoside-Binding, Soluble, 3; 14q22.3) and it has been reported to play a central role in self-tolerance, inflammation, and fibrosis. OBJECTIVE: To investigate associations among LGALS3 single nucleotide polymorphisms (SNPs) and serum levels Gal-3 and SSc susceptibility and their clinical features. METHODS: A case-control study with 88 patients and 151 matched controls was performed. LGALS3 variants were analyzed by the TaqMan real-time polymerase chain reaction (PCR) system whereas Gal-3 serum levels were measured by sandwich enzyme linked immunosorbent assay (ELISA). Associations among genotypes, clinical features, and Gal-3 levels were performed by univariable and multivariable analysis through statistical packages. RESULTS: The LGALS3 rs4652 A/C genotype was more frequent in SSc patients than controls according to overdominant model [OR 1.89 (CI 95% 1.01 - 3.52); p = .046]. Also, LGALS3 rs4652 C/C polymorphic genotype was associated with lower patient Gal-3 levels (p = .03) and control group (p = 0.005), as noted by generalized linear model (GLM). The LGALS3 rs1009977 G/T controls showed higher Gal-3 levels than wild-type and polymorphic genotypes (p = .03); however, in SSc patients, no difference was found. None of the LGALS3 SNPs or Gal-3 levels was associated with clinical manifestations in SSc patients. Considering only the SSc group, GLM analysis pointed LGALS3 rs4652 and rs2075601, pulmonary arterial hypertension (PAH), myopathy, and health assessment questionnaire (HAQ) and scleroderma health assessment questionnaire (SHAQ) as important predictors for Gal-3 levels. CONCLUSION: The LGALS3 rs4652 A/C was more frequent in SSc patients and related to lower Gal-3 levels. These findings were corroborated through a GLM to estimate Gal-3 values. Also, by model equations, Gal-3 levels may be predicted by HAQ, SHAQ, PAH, myopathy, and LGALS3 rs4652 and rs2075601 factors. In these ways, we suggest that galectins may be promising biomarkers to identify susceptibility to SSc as well as to identify HAQ, SHAQ, PAH, and myopathy outcomes.

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.