Integration of human stem cell-derived in vitro systems and mouse preclinical models identifies complex pathophysiologic mechanisms in retinal dystrophy.

Rare DRAM2 coding variants cause retinal dystrophy with early macular involvement via unknown mechanisms. We found that DRAM2 is ubiquitously expressed in the human eye and expression changes were observed in eyes with more common maculopathy such as Age-related Macular Degeneration (AMD). To gain insights into pathogenicity of DRAM2-related retinopathy, we used a combination of in vitro and in vivo models. We found that DRAM2 loss in human pluripotent stem cell (hPSC)-derived retinal organoids caused the presence of additional mesenchymal cells. Interestingly, Dram2 loss in mice also caused increased proliferation of cells from the choroid in vitro and exacerbated choroidal neovascular lesions in vivo. Furthermore, we observed that DRAM2 loss in human retinal pigment epithelial (RPE) cells resulted in increased susceptibility to stress-induced cell death in vitro and that Dram2 loss in mice caused age-related photoreceptor degeneration. This highlights the complexity of DRAM2 function, as its loss in choroidal cells provided a proliferative advantage, whereas its loss in post-mitotic cells, such as photoreceptor and RPE cells, increased degeneration susceptibility. Different models such as human pluripotent stem cell-derived systems and mice can be leveraged to study and model human retinal dystrophies; however, cell type and species-specific expression must be taken into account when selecting relevant systems.

A systems biology approach uncovers novel disease mechanisms in age-related macular degeneration.

Age-related macular degeneration (AMD) is a leading cause of blindness, affecting 200 million people worldwide. To identify genes that could be targeted for treatment, we created a molecular atlas at different stages of AMD. Our resource is comprised of RNA sequencing (RNA-seq) and DNA methylation microarrays from bulk macular retinal pigment epithelium (RPE)/choroid of clinically phenotyped normal and AMD donor eyes (n = 85), single-nucleus RNA-seq (164,399 cells), and single-nucleus assay for transposase-accessible chromatin (ATAC)-seq (125,822 cells) from the retina, RPE, and choroid of 6 AMD and 7 control donors. We identified 23 genome-wide significant loci differentially methylated in AMD, over 1,000 differentially expressed genes across different disease stages, and an AMD Müller state distinct from normal or gliosis. Chromatin accessibility peaks in genome-wide association study (GWAS) loci revealed putative causal genes for AMD, including HTRA1 and C6orf223. Our systems biology approach uncovered molecular mechanisms underlying AMD, including regulators of WNT signaling, FRZB and TLE2, as mechanistic players in disease.

Dose escalation randomised study of efmarodocokin alfa in healthy volunteers and patients with ulcerative colitis.

BACKGROUND: The interleukin-22 cytokine (IL-22) has demonstrated efficacy in preclinical colitis models with non-immunosuppressive mechanism of action. Efmarodocokin alfa (UTTR1147A) is a fusion protein agonist that links IL-22 to the crystallisable fragment (Fc) of human IgG4 for improved pharmacokinetic characteristics, but with a mutation to minimise Fc effector functions. METHODS: This randomised, phase 1b study evaluated the safety, tolerability, pharmacokinetics and pharmacodynamics of repeat intravenous dosing of efmarodocokin alfa in healthy volunteers (HVs; n=32) and patients with ulcerative colitis (n=24) at 30-90 µg/kg doses given once every 2 weeks or monthly (every 4 weeks) for 12 weeks (6:2 active:placebo per cohort). RESULTS: The most common adverse events (AEs) were on-target, reversible, dermatological effects (dry skin, erythema and pruritus). Dose-limiting non-serious dermatological AEs (severe dry skin, erythema, exfoliation and discomfort) were seen at 90 µg/kg once every 2 weeks (HVs, n=2; patients, n=1). Pharmacokinetics were generally dose-proportional across the dose levels, but patients demonstrated lower drug exposures relative to HVs at the same dose. IL-22 serum biomarkers and IL-22-responsive genes in colon biopsies were induced with active treatment, and microbiota composition changed consistent with a reversal in baseline dysbiosis. As a phase 1b study, efficacy endpoints were exploratory only. Clinical response was observed in 7/18 active-treated and 1/6 placebo-treated patients; clinical remission was observed in 5/18 active-treated and 0/6 placebo-treated patients. CONCLUSION: Efmarodocokin alfa had an adequate safety and pharmacokinetic profile in HVs and patients. Biomarker data confirmed IL-22R pathway activation in the colonic epithelium. Results support further investigation of this non-immunosuppressive potential inflammatory bowel disease therapeutic. TRIAL REGISTRATION NUMBER: NCT02749630.

Systems genomics in age-related macular degeneration.

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

Dual targeting of lymphocyte homing and retention through α4ß7 and αEß7 inhibition in inflammatory bowel disease.

Anti-integrins are therapeutically effective for inflammatory bowel disease, yet the relative contribution of α4ß7 and αEß7 to gut lymphocyte trafficking is not fully elucidated. Here, we evaluate the effect of α4ß7 and αEß7 blockade using a combination of murine models of gut trafficking and longitudinal gene expression analysis in etrolizumab-treated patients with Crohn's disease (CD). Dual blockade of α4ß7 and αEß7 reduces CD8+ T cell accumulation in the gut to a greater extent than blockade of either integrin alone. Anti-αEß7 reduces epithelial:T cell interactions and promotes egress of activated T cells from the mucosa into lymphatics. Inflammatory gene expression is greater in human intestinal αEß7+ T cells. Etrolizumab-treated patients with CD display a treatment-specific reduction in inflammatory and cytotoxic intraepithelial lymphocytes (IEL) genes. Concurrent blockade of α4ß7 and αEß7 promotes reduction of cytotoxic IELs and inflammatory T cells in the gut mucosa through a stepwise inhibition of intestinal tissue entry and retention.

Integration of eQTL and a Single-Cell Atlas in the Human Eye Identifies Causal Genes for Age-Related Macular Degeneration.

Age-related macular degeneration (AMD) is a leading cause of vision loss. To better understand disease pathogenesis and identify causal genes in GWAS loci for AMD risk, we present a comprehensive database of human retina and retinal pigment epithelium (RPE). Our database comprises macular and non-macular RNA sequencing (RNA-seq) profiles from 129 donors, a genome-wide expression quantitative trait loci (eQTL) dataset that includes macula-specific retina and RPE/choroid, and single-nucleus RNA-seq (NucSeq) from human retina and RPE with subtype resolution from more than 100,000 cells. Using NucSeq, we find enriched expression of AMD candidate genes in RPE cells. We identify 15 putative causal genes for AMD on the basis of co-localization of genetic association signals for AMD risk and eye eQTL, including the genes TSPAN10 and TRPM1. These results demonstrate the value of our human eye database for elucidating genetic pathways and potential therapeutic targets for ocular diseases.

Diesel exhaust particles dysregulate multiple immunological pathways in murine macrophages: Lessons from microarray and scRNA-seq technologies.

Exposure to ambient particulate matter has been shown to promote a variety of disorders, including cardiovascular diseases predominantly of ischemic etiology. However, the mechanisms linking inhaled particulates with systemic vascular effects, resulting in worsened atherosclerosis, are not well defined. We assessed the potential role of macrophages in translating these effects by analyzing gene expression patterns in response to diesel exhaust particles (DEP) at the average cell level, using Affymetrix microarrays in peritoneal macrophages in culture (in vitro), and at the individual cell level, using single-cell RNA sequencing (scRNA-seq) in alveolar macrophages collected from exposed mice (in vivo). Peritoneal macrophages were harvested from C57BL/6J mice and treated with 25 µg/mL of a DEP methanol extract (DEPe). These cells exhibited significant (FDR < 0.05) differential expression of a large number of genes and enrichment in pathways, especially engaged in immune responses and antioxidant defense. DEPe led to marked upregulation of heme oxygenase 1 (Hmox1), the most significantly upregulated gene (FDR = 1.75E-06), and several other antioxidant genes. For the in vivo work, C57BL/6J mice were subjected to oropharyngeal aspiration of 200 µg of whole DEP. The gene expression profiles of the alveolar macrophages harvested from these mice were analyzed at the single-cell level using scRNA-seq, which showed significant dysregulation of a broad number of genes enriched in immune system pathways as well, but with a large heterogeneity in how individual alveolar macrophages responded to DEP exposures. Altogether, DEP pollutants dysregulated immunological pathways in macrophages that may mediate the development of pulmonary and systemic vascular effects.

An Allosteric Anti-tryptase Antibody for the Treatment of Mast Cell-Mediated Severe Asthma.

Severe asthma patients with low type 2 inflammation derive less clinical benefit from therapies targeting type 2 cytokines and represent an unmet need. We show that mast cell tryptase is elevated in severe asthma patients independent of type 2 biomarker status. Active ß-tryptase allele count correlates with blood tryptase levels, and asthma patients carrying more active alleles benefit less from anti-IgE treatment. We generated a noncompetitive inhibitory antibody against human ß-tryptase, which dissociates active tetramers into inactive monomers. A 2.15 Å crystal structure of a ß-tryptase/antibody complex coupled with biochemical studies reveal the molecular basis for allosteric destabilization of small and large interfaces required for tetramerization. This anti-tryptase antibody potently blocks tryptase enzymatic activity in a humanized mouse model, reducing IgE-mediated systemic anaphylaxis, and inhibits airway tryptase in Ascaris-sensitized cynomolgus monkeys with favorable pharmacokinetics. These data provide a foundation for developing anti-tryptase as a clinical therapy for severe asthma.

Epigenome-wide association in adipose tissue from the METSIM cohort.

Most epigenome-wide association studies to date have been conducted in blood. However, metabolic syndrome is mediated by a dysregulation of adiposity and therefore it is critical to study adipose tissue in order to understand the effects of this syndrome on epigenomes. To determine if natural variation in DNA methylation was associated with metabolic syndrome traits, we profiled global methylation levels in subcutaneous abdominal adipose tissue. We measured association between 32 clinical traits related to diabetes and obesity in 201 people from the Metabolic Syndrome in Men cohort. We performed epigenome-wide association studies between DNA methylation levels and traits, and identified associations for 13 clinical traits in 21 loci. We prioritized candidate genes in these loci using expression quantitative trait loci, and identified 18 high confidence candidate genes, including known and novel genes associated with diabetes and obesity traits. Using methylation deconvolution, we examined which cell types may be mediating the associations, and concluded that most of the loci we identified were specific to adipocytes. We determined whether the abundance of cell types varies with metabolic traits, and found that macrophages increased in abundance with the severity of metabolic syndrome traits. Finally, we developed a DNA methylation-based biomarker to assess type 2 diabetes risk in adipose tissue. In conclusion, our results demonstrate that profiling DNA methylation in adipose tissue is a powerful tool for understanding the molecular effects of metabolic syndrome on adipose tissue, and can be used in conjunction with traditional genetic analyses to further characterize this disorder.

Altered methylation of specific DNA loci in the liver of Bhmt-null mice results in repression of Iqgap2 and F2rl2 and is associated with development of preneoplastic foci.

Folate B12-dependent remethylation of homocysteine is important, but less is understood about the importance of the alternative betaine-dependent methylation pathway-catalyzed by betaine-homocysteine methyltransferase (BHMT)-for establishing and maintaining adequate DNA methylation across the genome. We studied C57Bl/6J Bhmt (betaine-homocysteine methyltransferase)-null mice at age 4, 12, 24, and 52 wk (N = 8) and observed elevation of S-adenosylhomocysteine concentrations and development of preneoplastic foci in the liver (increased placental glutathione S-transferase and cytokeratin 8-18 activity; starting at 12 wk). At 4 wk, we identified 63 differentially methylated CpGs (DMCs; false discovery rate < 5%) proximal to 81 genes (across 14 chromosomes), of which 18 were differentially expressed. Of these DMCs, 52% were located in one 15.5-Mb locus on chromosome 13, which encompassed the Bhmt gene and defined a potentially sensitive region with mostly decreased methylation. Analyzing Hybrid Mouse Diversity Panel data, which consisted of 100 inbred strains of mice, we identified 97 DMCs that were affected by Bhmt genetic variation in the same region, with 7 overlapping those found in Bhmt-null mice (P < 0.001). At all time points, we found a hypomethylated region mapping to Iqgap2 (IQ motif-containing GTPase activating protein 2) and F2rl2 (proteinase-activated receptor-3), 2 genes that were also silenced and underexpressed, respectively.-Lupu, D. S., Orozco, L. D., Wang, Y., Cullen, J. M., Pellegrini, M., Zeisel, S. H. Altered methylation of specific DNA loci in the liver of Bhmt-null mice results in repression of Iqgap2 and F2rl2 and is associated with development of preneoplastic foci.

Mergeomics: multidimensional data integration to identify pathogenic perturbations to biological systems.

BACKGROUND: Complex diseases are characterized by multiple subtle perturbations to biological processes. New omics platforms can detect these perturbations, but translating the diverse molecular and statistical information into testable mechanistic hypotheses is challenging. Therefore, we set out to create a public tool that integrates these data across multiple datasets, platforms, study designs and species in order to detect the most promising targets for further mechanistic studies. RESULTS: We developed Mergeomics, a computational pipeline consisting of independent modules that 1) leverage multi-omics association data to identify biological processes that are perturbed in disease, and 2) overlay the disease-associated processes onto molecular interaction networks to pinpoint hubs as potential key regulators. Unlike existing tools that are mostly dedicated to specific data type or settings, the Mergeomics pipeline accepts and integrates datasets across platforms, data types and species. We optimized and evaluated the performance of Mergeomics using simulation and multiple independent datasets, and benchmarked the results against alternative methods. We also demonstrate the versatility of Mergeomics in two case studies that include genome-wide, epigenome-wide and transcriptome-wide datasets from human and mouse studies of total cholesterol and fasting glucose. In both cases, the Mergeomics pipeline provided statistical and contextual evidence to prioritize further investigations in the wet lab. The software implementation of Mergeomics is freely available as a Bioconductor R package. CONCLUSION: Mergeomics is a flexible and robust computational pipeline for multidimensional data integration. It outperforms existing tools, and is easily applicable to datasets from different studies, species and omics data types for the study of complex traits.

DNA Methylation Indicates Susceptibility to Isoproterenol-Induced Cardiac Pathology and Is Associated With Chromatin States.

RATIONALE: Only a small portion of the known heritability of cardiovascular diseases, such as heart failure, can be explained based on single-gene mutations. Chromatin structure and regulation provide a substrate through which genetic differences in noncoding regions may affect cellular function and response to disease, but the mechanisms are unknown. OBJECTIVE: We conducted genome-wide measurements of DNA methylation in different strains of mice that are susceptible and resistant to isoproterenol-induced dysfunction to test the hypothesis that this epigenetic mark may play a causal role in the development of heart failure. METHODS AND RESULTS: BALB/cJ and BUB/BnJ mice, determined to be susceptible and resistant to isoproterenol-induced heart failure, respectively, were administered the drug for 3 weeks via osmotic minipump. Reduced representational bisulfite sequencing was then used to compare the differences between the cardiac DNA methylomes in the basal state between strains and then after isoproterenol treatment. Single-base resolution DNA methylation measurements were obtained and revealed a bimodal distribution of methylation in the heart, enriched in lone intergenic CpGs and depleted from CpG islands around genes. Isoproterenol induced global decreases in methylation in both strains; however, the basal methylation pattern between strains shows striking differences that may be predictive of disease progression before environmental stress. The global correlation between promoter methylation and gene expression (as measured by microarray) was modest and revealed itself only with focused analyses of transcription start site and gene body regions (in contrast to when gene methylation was examined in toto). Modules of comethylated genes displayed correlation with other protein-based epigenetic marks, supporting the hypothesis that chromatin modifications act in a combinatorial manner to specify transcriptional phenotypes in the heart. CONCLUSIONS: This study provides the first single-base resolution map of the mammalian cardiac DNA methylome and the first case-control analysis of the changes in DNA methylation with heart failure. The findings demonstrate marked genetic differences in DNA methylation that are associated with disease progression.

Epigenome-wide association of liver methylation patterns and complex metabolic traits in mice.

Heritable epigenetic factors can contribute to complex disease etiology. Here we examine the contribution of DNA methylation to complex traits that are precursors to heart disease, diabetes, and osteoporosis. We profiled DNA methylation in the liver using bisulfite sequencing in 90 mouse inbred strains, genome-wide expression levels, proteomics, metabolomics, and 68 clinical traits and performed epigenome-wide association studies (EWAS). We found associations with numerous clinical traits including bone density, insulin resistance, expression, and protein and metabolite levels. A large proportion of associations were unique to EWAS and were not identified using GWAS. Methylation levels were regulated by genetics largely in cis, but we also found evidence of trans regulation, and we demonstrate that genetic variation in the methionine synthase reductase gene Mtrr affects methylation of hundreds of CpGs throughout the genome. Our results indicate that natural variation in methylation levels contributes to the etiology of complex clinical traits.

Intergenerational genomic DNA methylation patterns in mouse hybrid strains.

BACKGROUND: DNA methylation is a contributing factor to both rare and common human diseases, and plays a major role in development and gene silencing. While the variation of DNA methylation among individuals has been partially characterized, the degree to which methylation patterns are preserved across generations is still poorly understood. To determine the extent of methylation differences between two generations of mice we examined DNA methylation patterns in the livers of eight parental and F1 mice from C57BL/6J and DBA/2J mouse strains using bisulfite sequencing. RESULTS: We find a large proportion of reproducible methylation differences between C57BL/6J and DBA/2J chromosomes in CpGs, which are highly heritable between parent and F1 mice. We also find sex differences in methylation levels in 396 genes, and 11% of these are differentially expressed between females and males. Using a recently developed approach to identify allelically methylated regions independently of genotypic differences, we identify 112 novel putative imprinted genes and microRNAs, and validate imprinting at the RNA level in 10 of these genes. CONCLUSIONS: The majority of DNA methylation differences among individuals are associated with genetic differences, and a much smaller proportion of these epigenetic differences are due to sex, imprinting or stochastic intergenerational effects. Epigenetic differences can be a determining factor in heritable traits and should be considered in association studies for molecular and clinical traits, as we observed that methylation differences in the mouse model are highly heritable and can have functional consequences on molecular traits such as gene expression.

Intrauterine calorie restriction affects placental DNA methylation and gene expression.

Maternal nutrient restriction causes the development of adult onset chronic diseases in the intrauterine growth restricted (IUGR) fetus. Investigations in mice have shown that either protein or calorie restriction during pregnancy leads to glucose intolerance, increased fat mass, and hypercholesterolemia in adult male offspring. Some of these phenotypes are shown to persist in successive generations. The molecular mechanisms underlying IUGR remain unclear. The placenta is a critical organ for mediating changes in the environment and the development of embryos. To shed light on molecular mechanisms that might affect placental responses to differing environments we examined placentas from mice that had been exposed to different diets. We measured gene expression and whole genome DNA methylation in both male and female placentas of mice exposed to either caloric restriction or ad libitum diets. We observed several differentially expressed pathways associated with IUGR phenotypes and, most importantly, a significant decrease in the overall methylation between these groups as well as sex-specific effects that are more pronounced in males. In addition, a set of significantly differentially methylated genes that are enriched for known imprinted genes were identified, suggesting that imprinted loci may be particularly susceptible to diet effects. Lastly, we identified several differentially methylated microRNAs that target genes associated with immunological, metabolic, gastrointestinal, cardiovascular, and neurological chronic diseases, as well as genes responsible for transplacental nutrient transfer and fetal development.

Maximal information component analysis: a novel non-linear network analysis method.

BACKGROUND: Network construction and analysis algorithms provide scientists with the ability to sift through high-throughput biological outputs, such as transcription microarrays, for small groups of genes (modules) that are relevant for further research. Most of these algorithms ignore the important role of non-linear interactions in the data, and the ability for genes to operate in multiple functional groups at once, despite clear evidence for both of these phenomena in observed biological systems. RESULTS: We have created a novel co-expression network analysis algorithm that incorporates both of these principles by combining the information-theoretic association measure of the maximal information coefficient (MIC) with an Interaction Component Model. We evaluate the performance of this approach on two datasets collected from a large panel of mice, one from macrophages and the other from liver by comparing the two measures based on a measure of module entropy, Gene Ontology (GO) enrichment, and scale-free topology (SFT) fit. Our algorithm outperforms a widely used co-expression analysis method, weighted gene co-expression network analysis (WGCNA), in the macrophage data, while returning comparable results in the liver dataset when using these criteria. We demonstrate that the macrophage data has more non-linear interactions than the liver dataset, which may explain the increased performance of our method, termed Maximal Information Component Analysis (MICA) in that case. CONCLUSIONS: In making our network algorithm more accurately reflect known biological principles, we are able to generate modules with improved relevance, particularly in networks with confounding factors such as gene by environment interactions.