Global endometrial DNA methylation analysis reveals insights into mQTL regulation and associated endometriosis disease risk and endometrial function

Epigenomics; Genetic predisposition to disease; Urogenital reproductive disordersEpigenómica; Predisposición genética a la enfermedad; Trastornos reproductivos urogenitalesEpigenòmica; Predisposició genètica a la malaltia; Trastorns reproductius urogenitalsEndometriosis is a leading cause of pain and infertility affecting millions of women globally. Herein, we characterize variation in DNA methylation (DNAm) and its association with menstrual cycle phase, endometriosis, and genetic variants through analysis of genotype data and methylation in endometrial samples from 984 deeply-phenotyped participants. We estimate that 15.4% of the variation in endometriosis is captured by DNAm and identify significant differences in DNAm profiles associated with stage III/IV endometriosis, endometriosis sub-phenotypes and menstrual cycle phase, including opening of the window for embryo implantation. Menstrual cycle phase was a major source of DNAm variation suggesting cellular and hormonally-driven changes across the cycle can regulate genes and pathways responsible for endometrial physiology and function. DNAm quantitative trait locus (mQTL) analysis identified 118,185 independent cis-mQTLs including 51 associated with risk of endometriosis, highlighting candidate genes contributing to disease risk. Our work provides functional evidence for epigenetic targets contributing to endometriosis risk and pathogenesis. Data generated serve as a valuable resource for understanding tissue-specific effects of methylation on endometrial biology in health and disease.This work has been supported by the National Institutes of Health (NIH) NICHD R01 HD089511. It was also supported, in part, by funding from Wellbeing of Women (through sponsorship from PwC) (R42533) and the Medical Research Council (MR/N024524/1 and MR/N022556/1) and NIH HD094842 (Harvard/MSU). K.K. was supported by NIH NCI R37 CA233774. A.F.M. was supported by an Australian Research Council Future Fellowship (FT200100837). G.W.M. was supported by NHMRC Fellowship (GNT1177194)

Global Analysis of Transcription Start Sites and Enhancers in Endometrial Stromal Cells and Differences Associated with Endometriosis.

Identifying tissue-specific molecular signatures of active regulatory elements is critical to understanding gene regulatory mechanisms. In this study, transcription start sites (TSS) and enhancers were identified using Cap analysis of gene expression (CAGE) across endometrial stromal cell (ESC) samples obtained from women with (n = 4) and without endometriosis (n = 4). ESC TSSs and enhancers were compared to those reported in other tissue and cell types in FANTOM5 and were integrated with RNA-seq and ATAC-seq data from the same samples for regulatory activity and network analyses. CAGE tag count differences between women with and without endometriosis were statistically tested and tags within close proximity to genetic variants associated with endometriosis risk were identified. Over 90% of tag clusters mapping to promoters were observed in cells and tissues in FANTOM5. However, some potential cell-type-specific promoters and enhancers were also observed. Regions of open chromatin identified using ATAC-seq provided further evidence of the active transcriptional regions identified by CAGE. Despite the small sample number, there was evidence of differences associated with endometriosis at 210 consensus clusters, including IGFBP5, CALD1 and OXTR. ESC TSSs were also located within loci associated with endometriosis risk from genome-wide association studies. This study provides novel evidence of transcriptional differences in endometrial stromal cells associated with endometriosis and provides a valuable cell-type specific resource of active TSSs and enhancers in endometrial stromal cells

Altered differentiation of endometrial mesenchymal stromal fibroblasts is associated with endometriosis susceptibility.

Cellular development is tightly regulated as mature cells with aberrant functions may initiate pathogenic processes. The endometrium is a highly regenerative tissue, shedding and regenerating each month. Endometrial stromal fibroblasts are regenerated each cycle from mesenchymal stem cells and play a pivotal role in endometriosis, a disease characterised by endometrial cells that grow outside the uterus. Why the cells of some women are more capable of developing into endometriosis lesions is not clear. Using isolated, purified and cultured endometrial cells of mesenchymal origin from 19 women with (n = 10) and without (n = 9) endometriosis we analysed the transcriptome of 33,758 individual cells and compared these to clinical characteristics and in vitro growth profiles. We show purified mesenchymal cell cultures include a mix of mesenchymal stem cells and two endometrial stromal fibroblast subtypes with distinct transcriptomic signatures indicative of varied progression through the differentiation processes. The fibroblast subgroup characterised by incomplete differentiation was predominantly (81%) derived from women with endometriosis and exhibited an altered in vitro growth profile. These results uncover an inherent difference in endometrial cells of women with endometriosis and highlight the relevance of cellular differentiation and its potential to contribute to disease susceptibility

Genotype-free demultiplexing of pooled single-cell RNA-seq.

A variety of methods have been developed to demultiplex pooled samples in a single cell RNA sequencing (scRNA-seq) experiment which either require hashtag barcodes or sample genotypes prior to pooling. We introduce scSplit which utilizes genetic differences inferred from scRNA-seq data alone to demultiplex pooled samples. scSplit also enables mapping clusters to original samples. Using simulated, merged, and pooled multi-individual datasets, we show that scSplit prediction is highly concordant with demuxlet predictions and is highly consistent with the known truth in cell-hashing dataset. scSplit is ideally suited to samples without external genotype information and is available at: https://github.com/jon-xu/scSplit

The genetic basis of endometriosis and comorbidity with other pain and inflammatory conditions

Endometriosis is a common condition associated with debilitating pelvic pain and infertility. A genome-wide association study meta-analysis, including 60,674 cases and 701,926 controls of European and East Asian descent, identified 42 genome-wide significant loci comprising 49 distinct association signals. Effect sizes were largest for stage 3/4 disease, driven by ovarian endometriosis. Identified signals explained up to 5.01% of disease variance and regulated expression or methylation of genes in endometrium and blood, many of which were associated with pain perception/maintenance (SRP14/BMF, GDAP1, MLLT10, BSN and NGF). We observed significant genetic correlations between endometriosis and 11 pain conditions, including migraine, back and multisite chronic pain (MCP), as well as inflammatory conditions, including asthma and osteoarthritis. Multitrait genetic analyses identified substantial sharing of variants associated with endometriosis and MCP/migraine. Targeted investigations of genetically regulated mechanisms shared between endometriosis and other pain conditions are needed to aid the development of new treatments and facilitate early symptomatic intervention

The genetic basis of endometriosis and comorbidity with other pain and inflammatory conditions

Endometriosis is a common condition associated with debilitating pelvic pain and infertility. A genome-wide association study meta-analysis, including 60,674 cases and 701,926 controls of European and East Asian descent, identified 42 genome-wide significant loci comprising 49 distinct association signals. Effect sizes were largest for stage 3/4 disease, driven by ovarian endometriosis. Identified signals explained up to 5.01% of disease variance and regulated expression or methylation of genes in endometrium and blood, many of which were associated with pain perception/maintenance (SRP14/BMF, GDAP1, MLLT10, BSN and NGF). We observed significant genetic correlations between endometriosis and 11 pain conditions, including migraine, back and multisite chronic pain (MCP), as well as inflammatory conditions, including asthma and osteoarthritis. Multitrait genetic analyses identified substantial sharing of variants associated with endometriosis and MCP/migraine. Targeted investigations of genetically regulated mechanisms shared between endometriosis and other pain conditions are needed to aid the development of new treatments and facilitate early symptomatic intervention

Genotype Data for Bullmastiff Dogs

Genome-wide SNP genotypes for 194 Bullmastiff dog

Genomic diversity and lymphoma predisposition in bullmastiffs

Lymphoma, cancer of the lymphocytes, is one of the most common malignancies seen in dogs. Studies have reported a high incidence of lymphoma in particular breeds suggesting a genetic component to the disease. The aim of this thesis was to investigate genetic factors influencing lymphoma predisposition in Bullmastiffs, one of the most at risk breeds, and assess the level of genetic diversity in the breed. Analysis of gene expression data from resting and stimulated lymphocytes in dogs characterized genes, and potential biomarkers, involved in the early response of T-cells to activation in dogs. Genealogical and molecular methods were used to assess the level of diversity in the Bullmastiff population. The small effective population size, popular sire use and presence of subpopulations highlight the need to maintain diversity and monitor use of breeding individuals. The application of NetView, a population analysis pipeline, in dogs was also trialed in German shepherd dogs, validating the pipelines ability to provide a more easily interpretable visualisation of the fine scale population structure. A genome-wide association analysis and fine-scale mapping detected an association between a 1.2Mb region on CFA13 containing five risk haplotypes and four candidate genes (MYC, miR-1204, miR-1205 and miR-1206) and early onset lymphoma in Bullmastiffs. A difference in haplotype frequency and homozygosity across two risk haplotypes was identified between high and low risk breeds, highlighting these two regions as candidates. Additional sources of variation in the Bullmastiff were identified from genotype and whole genome sequence data using copy number variant detection and variant discovery algorithms. Together this information has expanded understanding of the genomic factors influencing lymphoma risk in dogs and may be used by breeders to develop breeding strategies to reduce the incidence of lymphoma in the Bullmastiff population whilst maintaining diversity

T-cell activation and early gene response in dogs

T-cells play a crucial role in canine immunoregulation and defence against invading pathogens. Proliferation is fundamental to T-cell differentiation, homeostasis and immune response. Initiation of proliferation following receptor mediated stimuli requires a temporally programmed gene response that can be identified as immediate-early, mid- and late phases. The immediate-early response genes in T-cell activation engage the cell cycle machinery and promote subsequent gene activation events. Genes involved in this immedi-ate-early response in dogs are yet to be identified. The present study was undertaken to characterise the early T-cell gene response in dogs to improve understanding of the genetic mechanisms regulating immune function. Gene expression profiles were characterised using canine gene expression microarrays and quantitative reverse transcription PCR (qRT-PCR), and paired samples from eleven dogs. Significant functional annotation clusters were identified following stimulation with phytohemagluttinin (PHA) (5 μg/ml), including the Toll-like receptor signaling pathway and phosphorylation pathways. Using strict statistical criteria, 13 individual genes were found to be differentially expressed, nine of which have ontologies that relate to proliferation and cell cycle control. These included, prosta-glandin-endoperoxide synthase 2 (PTGS2/COX2), early growth response 1 (EGR1), growth arrest and DNA damage-inducible gene (GADD45B ), phorbol-12-myristate-13-acetateinduced protein 1 (PMAIP1), V-FOS FBJ murine osteosarcoma viral oncogene homolog (FOS), early growth response 2 (EGR2), hemogen (HEMGN ), polo-like kinase 2 (PLK2) and polo-like kinase 3 (PLK3). Differential gene expression was re-examined using qRT-PCR, which confirmed that EGR1, EGR2, PMAIP1, PTGS2, FOS and GADD45B were significantly upregulated in stimulated cells and ALAS2 downregulated. PTGS2 and EGR1 showed the highest levels of response in these dogs. Both of these genes are involved in cell cycle regulation. This study provides a comprehensive analysis of the early T-cell gene response to activation in dogs