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.

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.

Genomic landscape of allelic imbalance in premalignant atypical adenomatous hyperplasias of the lung.

BACKGROUND: Genomic investigation of atypical adenomatous hyperplasia (AAH), the only known precursor lesion to lung adenocarcinomas (LUAD), presents challenges due to the low mutant cell fractions. This necessitates sensitive methods for detection of chromosomal aberrations to better study the role of critical alterations in early lung cancer pathogenesis and the progression from AAH to LUAD. METHODS: We applied a sensitive haplotype-based statistical technique to detect chromosomal alterations leading to allelic imbalance (AI) from genotype array profiling of 48 matched normal lung parenchyma, AAH and tumor tissues from 16 stage-I LUAD patients. To gain insights into shared developmental trajectories among tissues, we performed phylogenetic analyses and integrated our results with point mutation data, highlighting significantly-mutated driver genes in LUAD pathogenesis. FINDINGS: AI was detected in nine AAHs (56%). Six cases exhibited recurrent loss of 17p. AI and the enrichment of 17p events were predominantly identified in patients with smoking history. Among the nine AAH tissues with detected AI, seven exhibited evidence for shared chromosomal aberrations with matched LUAD specimens, including losses harboring tumor suppressors on 17p, 8p, 9p, 9q, 19p, and gains encompassing oncogenes on 8q, 12p and 1q. INTERPRETATION: Chromosomal aberrations, particularly 17p loss, appear to play critical roles early in AAH pathogenesis. Genomic instability in AAH, as well as truncal chromosomal aberrations shared with LUAD, provide evidence for mutation accumulation and are suggestive of a cancerized field contributing to the clonal selection and expansion of these premalignant lesions. FUND: Supported in part by Cancer Prevention and Research Institute of Texas (CPRIT) grant RP150079 (PS and HK), NIH grant R01HG005859 (PS) and The University of Texas MD Anderson Cancer Center Core Support Grant.

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.