Altered mechanisms of genital development identified through integration of DNA methylation and genomic measures in hypospadias.

Hypospadias is a common birth defect where the urethral opening forms on the ventral side of the penis. We performed integrative methylomic, genomic, and transcriptomic analyses to characterize sites of DNA methylation that influence genital development. In case-control and case-only epigenome-wide association studies (EWAS) of preputial tissue we identified 25 CpGs associated with hypospadias characteristics and used one-sample two stage least squares Mendelian randomization (2SLS MR) to show a causal relationship for 21 of the CpGs. The largest difference was 15.7% lower beta-value at cg14436889 among hypospadias cases than controls (EWAS P = 5.4e-7) and is likely causal (2SLS MR P = 9.8e-15). Integrative annotation using two-sample Mendelian randomization of these methylation regions highlight potentially causal roles of genes involved in germ layer differentiation (WDHD1, DNM1L, TULP3), beta-catenin signaling (PKP2, UBE2R2, TNKS), androgens (CYP4A11, CYP4A22, CYP4B1, CYP4X1, CYP4Z2P, EPHX1, CD33/SIGLEC3, SIGLEC5, SIGLEC7, KLK5, KLK7, KLK10, KLK13, KLK14), and reproductive traits (ACAA1, PLCD1, EFCAB4B, GMCL1, MKRN2, DNM1L, TEAD4, TSPAN9, KLK family). This study identified CpGs that remained differentially methylated after urogenital development and used the most relevant tissue sample available to study hypospadias. We identified multiple methylation sites and candidate genes that can be further evaluated for their roles in regulating urogenital development.

Inherited genetic susceptibility to acute lymphoblastic leukemia in Down syndrome.

Children with Down syndrome (DS) have a 20-fold increased risk of acute lymphoblastic leukemia (ALL) and distinct somatic features, including CRLF2 rearrangement in ∼50% of cases; however, the role of inherited genetic variation in DS-ALL susceptibility is unknown. We report the first genome-wide association study of DS-ALL, comprising a meta-analysis of 4 independent studies, with 542 DS-ALL cases and 1192 DS controls. We identified 4 susceptibility loci at genome-wide significance: rs58923657 near IKZF1 (odds ratio [OR], 2.02; Pmeta = 5.32 × 10-15), rs3731249 in CDKN2A (OR, 3.63; Pmeta = 3.91 × 10-10), rs7090445 in ARID5B (OR, 1.60; Pmeta = 8.44 × 10-9), and rs3781093 in GATA3 (OR, 1.73; Pmeta = 2.89 × 10-8). We performed DS-ALL vs non-DS ALL case-case analyses, comparing risk allele frequencies at these and other established susceptibility loci (BMI1, PIP4K2A, and CEBPE) and found significant association with DS status for CDKN2A (OR, 1.58; Pmeta = 4.1 × 10-4). This association was maintained in separate regression models, both adjusting for and stratifying on CRLF2 overexpression and other molecular subgroups, indicating an increased penetrance of CDKN2A risk alleles in children with DS. Finally, we investigated functional significance of the IKZF1 risk locus, and demonstrated mapping to a B-cell super-enhancer, and risk allele association with decreased enhancer activity and differential protein binding. IKZF1 knockdown resulted in significantly higher proliferation in DS than non-DS lymphoblastoid cell lines. Our findings demonstrate a higher penetrance of the CDKN2A risk locus in DS and serve as a basis for further biological insights into DS-ALL etiology.

The role of genetic variation in DGKK on moderate and severe hypospadias.

BACKGROUND: Recent genome-wide association studies of hypospadias have implicated the role of genetic variants in or near the diacylglycerol kinase kappa (DGKK) gene. However, these variants are largely identified among samples of mild and moderate hypospadias cases. Therefore, we evaluated previously identified DGKK variants among second- and third-degree hypospadias cases and controls recruited in Arkansas, a state characterized by a high birth prevalence of hypospadias. METHODS: Second- and third-degree hypospadias non-Hispanic white cases (n = 36 and n = 9, respectively) and controls (n = 45) were recruited at Arkansas Children's Hospital. Preputial tissue was collected on cases and controls between 2013 and 2017. Cases and controls were genotyped using the Illumina Infinium Global Screening Array. We used logistic regression models to assess the association of genotyped and imputed genetic variants mapped to the DGKK region with second- and third-degree hypospadias. RESULTS: All families self-reported as non-Hispanic white and genetic principal component analyses did not demonstrate evidence of population stratification. Five DGKK variants previously reported as associated with hypospadias were identified in the genotype data. None of the variants were associated with second- or third-degree hypospadias (range of odds ratios = 0.7-0.9, all p > .05). CONCLUSIONS: In our analyses, genetic variation in DGKK does not play a role in the development of moderate and severe hypospadias. Our findings provide support to the etiologic heterogeneity of hypospadias by all classifications of severity.

Genome-wide analysis of DNA methylation in buccal cells: a study of monozygotic twins and mQTLs.

BACKGROUND: DNA methylation arrays are widely used in epigenome-wide association studies and methylation quantitative trait locus (mQTL) studies. Here, we performed the first genome-wide analysis of monozygotic (MZ) twin correlations and mQTLs on data obtained with the Illumina MethylationEPIC BeadChip (EPIC array) and compared the performance of the EPIC array to the Illumina HumanMethylation450 BeadChip (HM450 array) for buccal-derived DNA. RESULTS: Good-quality EPIC data were obtained for 102 buccal-derived DNA samples from 49 MZ twin pairs (mean age = 7.5 years, range = 1-10). Differences between MZ twins in the cellular content of buccal swabs were a major driver for differences in their DNA methylation profiles, highlighting the importance to adjust for cellular composition in DNA methylation studies of buccal-derived DNA. After adjusting for cellular composition, the genome-wide mean correlation (r) between MZ twins was 0.21 for the EPIC array, and cis mQTL analysis in 84 twins identified 1,296,323 significant associations (FDR 5%), encompassing 33,749 methylation sites and 616,029 genetic variants. MZ twin correlations were slightly larger (p < 2.2 × 10-16) for novel EPIC probes (N = 383,066, mean r = 0.22) compared to probes that are also present on HM450 (N = 406,822, mean r = 0.20). In line with this observation, a larger percentage of novel EPIC probes was associated with genetic variants (novel EPIC probes with significant mQTL 4.7%, HM450 probes with mQTL 3.9%, p < 2.2 × 10-16). Methylation sites with a large MZ correlation and sites associated with mQTLs were most strongly enriched in epithelial cell DNase I hypersensitive sites (DHSs), enhancers, and histone mark H3K4me3. CONCLUSIONS: We conclude that the contribution of familial factors to individual differences in DNA methylation and the effect of mQTLs are larger for novel EPIC probes, especially those within regulatory elements connected to active regions specific to the investigated tissue.