Chromosomal microarray analysis supplements exome sequencing to diagnose children with suspected inborn errors of immunity.

Purpose: Though copy number variants (CNVs) have been suggested to play a significant role in inborn errors of immunity (IEI), the precise nature of this role remains largely unexplored. We sought to determine the diagnostic contribution of CNVs using genome-wide chromosomal microarray analysis (CMA) in children with IEI. Methods: We performed exome sequencing (ES) and CMA for 332 unrelated pediatric probands referred for evaluation of IEI. The analysis included primary, secondary, and incidental findings. Results: Of the 332 probands, 134 (40.4%) received molecular diagnoses. Of these, 116/134 (86.6%) were diagnosed by ES alone. An additional 15/134 (11.2%) were diagnosed by CMA alone, including two likely de novo changes. Three (2.2%) participants had diagnostic molecular findings from both ES and CMA, including two compound heterozygotes and one participant with two distinct diagnoses. Half of the participants with CMA contribution to diagnosis had CNVs in at least one non-immune gene, highlighting the clinical complexity of these cases. Overall, CMA contributed to 18/134 diagnoses (13.4%), increasing the overall diagnostic yield by 15.5% beyond ES alone. Conclusion: Pairing ES and CMA can provide a comprehensive evaluation to clarify the complex factors that contribute to both immune and non-immune phenotypes. Such a combined approach to genetic testing helps untangle complex phenotypes, not only by clarifying the differential diagnosis, but in some cases by identifying multiple diagnoses contributing to the overall clinical presentation.

DNA methylation age in paired tumor and adjacent normal breast tissue in Chinese women with breast cancer.

BACKGROUND: Few studies have examined epigenetic age acceleration (AA), the difference between DNA methylation (DNAm) predicted age and chronological age, in relation to somatic genomic features in paired cancer and normal tissue, with less work done in non-European populations. In this study, we aimed to examine DNAm age and its associations with breast cancer risk factors, subtypes, somatic genomic profiles including mutation and copy number alterations and other aging markers in breast tissue of Chinese breast cancer (BC) patients from Hong Kong. METHODS: We performed genome-wide DNA methylation profiling of 196 tumor and 188 paired adjacent normal tissue collected from Chinese BC patients in Hong Kong (HKBC) using Illumina MethylationEPIC array. The DNAm age was calculated using Horvath's pan-tissue clock model. Somatic genomic features were based on data from RNA sequencing (RNASeq), whole-exome sequencing (WES), and whole-genome sequencing (WGS). Pearson's correlation (r), Kruskal-Wallis test, and regression models were used to estimate associations of DNAm AA with somatic features and breast cancer risk factors. RESULTS: DNAm age showed a stronger correlation with chronological age in normal (Pearson r = 0.78, P < 2.2e-16) than in tumor tissue (Pearson r = 0.31, P = 7.8e-06). Although overall DNAm age or AA did not vary significantly by tissue within the same individual, luminal A tumors exhibited increased DNAm AA (P = 0.004) while HER2-enriched/basal-like tumors exhibited markedly lower DNAm AA (P = < .0001) compared with paired normal tissue. Consistent with the subtype association, tumor DNAm AA was positively correlated with ESR1 (Pearson r = 0.39, P = 6.3e-06) and PGR (Pearson r = 0.36, P = 2.4e-05) gene expression. In line with this, we found that increasing DNAm AA was associated with higher body mass index (P = 0.039) and earlier age at menarche (P = 0.035), factors that are related to cumulative exposure to estrogen. In contrast, variables indicating extensive genomic instability, such as TP53 somatic mutations, high tumor mutation/copy number alteration burden, and homologous repair deficiency were associated with lower DNAm AA. CONCLUSIONS: Our findings provide additional insights into the complexity of breast tissue aging that is associated with the interaction of hormonal, genomic, and epigenetic mechanisms in an East Asian population.

Targeted long-read sequencing of the Ewing sarcoma 6p25.1 susceptibility locus identifies germline-somatic interactions with EWSR1-FLI1 binding.

Ewing sarcoma (EwS) is a rare bone and soft tissue malignancy driven by chromosomal translocations encoding chimeric transcription factors, such as EWSR1-FLI1, that bind GGAA motifs forming novel enhancers that alter nearby expression. We propose that germline microsatellite variation at the 6p25.1 EwS susceptibility locus could impact downstream gene expression and EwS biology. We performed targeted long-read sequencing of EwS blood DNA to characterize variation and genomic features important for EWSR1-FLI1 binding. We identified 50 microsatellite alleles at 6p25.1 and observed that EwS-affected individuals had longer alleles (>135 bp) with more GGAA repeats. The 6p25.1 GGAA microsatellite showed chromatin features of an EWSR1-FLI1 enhancer and regulated expression of RREB1, a transcription factor associated with RAS/MAPK signaling. RREB1 knockdown reduced proliferation and clonogenic potential and reduced expression of cell cycle and DNA replication genes. Our integrative analysis at 6p25.1 details increased binding of longer GGAA microsatellite alleles with acquired EWSR-FLI1 to promote Ewing sarcomagenesis by RREB1-mediated proliferation.

Immunogenetics associated with severe coccidioidomycosis.

Disseminated coccidioidomycosis (DCM) is caused by Coccidioides, pathogenic fungi endemic to the southwestern United States and Mexico. Illness occurs in approximately 30% of those infected, less than 1% of whom develop disseminated disease. To address why some individuals allow dissemination, we enrolled patients with DCM and performed whole-exome sequencing. In an exploratory set of 67 patients with DCM, 2 had haploinsufficient STAT3 mutations, and defects in ß-glucan sensing and response were seen in 34 of 67 cases. Damaging CLEC7A and PLCG2 variants were associated with impaired production of ß-glucan-stimulated TNF-α from PBMCs compared with healthy controls. Using ancestry-matched controls, damaging CLEC7A and PLCG2 variants were overrepresented in DCM, including CLEC7A Y238* and PLCG2 R268W. A validation cohort of 111 patients with DCM confirmed the PLCG2 R268W, CLEC7A I223S, and CLEC7A Y238* variants. Stimulation with a DECTIN-1 agonist induced DUOX1/DUOXA1-derived hydrogen peroxide [H2O2] in transfected cells. Heterozygous DUOX1 or DUOXA1 variants that impaired H2O2 production were overrepresented in discovery and validation cohorts. Patients with DCM have impaired ß-glucan sensing or response affecting TNF-α and H2O2 production. Impaired Coccidioides recognition and decreased cellular response are associated with disseminated coccidioidomycosis.

Clinical exome sequencing of 1000 families with complex immune phenotypes: Toward comprehensive genomic evaluations.

BACKGROUND: Prospective genetic evaluation of patients at this referral research hospital presents clinical research challenges. OBJECTIVES: This study sought not only a single-gene explanation for participants' immune-related presentations, but viewed each participant holistically, with the potential to have multiple genetic contributions to their immune phenotype and other heritable comorbidities relevant to their presentation and health. METHODS: This study developed a program integrating exome sequencing, chromosomal microarray, phenotyping, results return with genetic counseling, and reanalysis in 1505 individuals from 1000 families with suspected or known inborn errors of immunity. RESULTS: Probands were 50.8% female, 71.5% were ≥18 years, and had diverse immune presentations. Overall, 327 of 1000 probands (32.7%) received 361 molecular diagnoses. These included 17 probands with diagnostic copy number variants, 32 probands with secondary findings, and 31 probands with multiple molecular diagnoses. Reanalysis added 22 molecular diagnoses, predominantly due to new disease-gene associations (9 of 22, 40.9%). One-quarter of the molecular diagnoses (92 of 361) did not involve immune-associated genes. Molecular diagnosis was correlated with younger age, male sex, and a higher number of organ systems involved. This program also facilitated the discovery of new gene-disease associations such as SASH3-related immunodeficiency. A review of treatment options and ClinGen actionability curations suggest that at least 251 of 361 of these molecular diagnoses (69.5%) could translate into ≥1 management option. CONCLUSIONS: This program contributes to our understanding of the diagnostic and clinical utility whole exome analysis on a large scale.

Age-related DNA methylation in paired normal and tumour breast tissue in Chinese breast cancer patients.

Age-related DNA methylation is a potential mechanism contributing to breast cancer development. Studies of primarily Caucasian women have identified many CpG sites of age-related methylation in non-diseased breast tissue possibly driving cancer development over time. There is a paucity of studies involving Asian women whose ages at breast cancer onset are usually younger than Caucasians. We identified the 181 most consistent age-related methylation events in non-diseased breast tissue across published studies. Age-related methylation events were measured in adjacent normal and breast tumour tissue in an exclusively Asian population at the previously identified age-related methylation sites. Age-related methylation was found in 118 probes in adjacent normal breast tissue. Methylation of 99% of these sites was increased with age and predominantly located on CpG islands in promoter regions. To ascertain biological relevance to breast cancer, we focused on the 37 sites with overall higher methylation in tumour compared to adjacent normal samples. Some sites positively related to age, including AQP5 and CORO6, inversely correlated with gene expression. Several others have known involvement in suppression of carcinogenesis including GPC5 and SST, suggesting that perturbation of epigenetic regulation at these sites due to ageing may contribute to the progression of carcinogenesis. This study highlights an age-related methylation landscape in non-tumour tissue, consistent not just across studies, but also across different populations. We present candidate age-related methylation sites warranting further investigation as potential epigenetic drivers of breast cancer. They may serve as potential targets of site-specific demethylation intervention strategies for the prevention of age-related breast cancer.

Low-frequency variation near common germline susceptibility loci are associated with risk of Ewing sarcoma.

BACKGROUND: Ewing sarcoma (EwS) is a rare, aggressive solid tumor of childhood, adolescence and young adulthood associated with pathognomonic EWSR1-ETS fusion oncoproteins altering transcriptional regulation. Genome-wide association studies (GWAS) have identified 6 common germline susceptibility loci but have not investigated low-frequency inherited variants with minor allele frequencies below 5% due to limited genotyped cases of this rare tumor. METHODS: We investigated the contribution of rare and low-frequency variation to EwS susceptibility in the largest EwS genome-wide association study to date (733 EwS cases and 1,346 unaffected controls of European ancestry). RESULTS: We identified two low-frequency variants, rs112837127 and rs2296730, on chromosome 20 that were associated with EwS risk (OR = 0.186 and 2.038, respectively; P-value < 5×10-8) and located near previously reported common susceptibility loci. After adjusting for the most associated common variant at the locus, only rs112837127 remained a statistically significant independent signal (OR = 0.200, P-value = 5.84×10-8). CONCLUSIONS: These findings suggest rare variation residing on common haplotypes are important contributors to EwS risk. IMPACT: Motivate future targeted sequencing studies for a comprehensive evaluation of low-frequency and rare variation around common EwS susceptibility loci.

Subsequent Neoplasm Risk Associated With Rare Variants in DNA Damage Response and Clinical Radiation Sensitivity Syndrome Genes in the Childhood Cancer Survivor Study.

PURPOSE: Radiotherapy for childhood cancer is associated with elevated subsequent neoplasm (SN) risk, but the contribution of rare variants in DNA damage response and radiation sensitivity genes to SN risk is unknown. PATIENTS AND METHODS: We conducted whole-exome sequencing in a cohort of childhood cancer survivors originally diagnosed during 1970 to 1986 (mean follow-up, 32.7 years), with reconstruction of doses to body regions from radiotherapy records. We identified patients who developed SN types previously reported to be related to radiotherapy (RT-SNs; eg, basal cell carcinoma [BCC], breast cancer, meningioma, thyroid cancer, sarcoma) and matched controls (sex, childhood cancer type/diagnosis, age, SN location, radiation dose, survival). Conditional logistic regression assessed SN risk associated with potentially protein-damaging rare variants (SnpEff, ClinVar) in 476 DNA damage response or radiation sensitivity genes with exact permutation-based P values using a Bonferroni-corrected significance threshold of P < 8.06 × 10-5. RESULTS: Among 5,105 childhood cancer survivors of European descent, 1,108 (21.7%) developed at least 1 RT-SN. Out-of-field RT-SN risk, excluding BCC, was associated with homologous recombination repair (HRR) gene variants (patient cases, 23.2%; controls, 10.8%; odds ratio [OR], 2.6; 95% CI, 1.7 to 3.9; P = 4.79 × 10-5), most notably but nonsignificantly for FANCM (patient cases, 4.0%; matched controls, 0.6%; P = 9.64 × 10-5). HRR variants were not associated with likely in/near-field RT-SNs, excluding BCC (patient cases, 12.7%; matched controls, 12.9%; P = .92). Irrespective of radiation dose, risk for RT-SNs was also associated with EXO1 variants (patient cases, 1.8%; controls, 0.4%; P = 3.31 × 10-5), another gene implicated in DNA double-strand break repair. CONCLUSION: In this large-scale discovery study, we identified novel associations between RT-SN risk after childhood cancer and potentially protein-damaging rare variants in genes involved in DNA double-strand break repair, particularly HRR. With replication, these results could affect screening recommendations for childhood cancer survivors and risk-benefit assessments of treatment approaches.

Frequency of Pathogenic Germline Variants in Cancer-Susceptibility Genes in Patients With Osteosarcoma.

Importance: Osteosarcoma, the most common malignant bone tumor in children and adolescents, occurs in a high number of cancer predisposition syndromes that are defined by highly penetrant germline mutations. The germline genetic susceptibility to osteosarcoma outside of familial cancer syndromes remains unclear. Objective: To investigate the germline genetic architecture of 1244 patients with osteosarcoma. Design, Setting, and Participants: Whole-exome sequencing (n = 1104) or targeted sequencing (n = 140) of the DNA of 1244 patients with osteosarcoma from 10 participating international centers or studies was conducted from April 21, 2014, to September 1, 2017. The results were compared with the DNA of 1062 individuals without cancer assembled internally from 4 participating studies who underwent comparable whole-exome sequencing and 27 173 individuals of non-Finnish European ancestry who were identified through the Exome Aggregation Consortium (ExAC) database. In the analysis, 238 high-interest cancer-susceptibility genes were assessed followed by testing of the mutational burden across 736 additional candidate genes. Principal component analyses were used to identify 732 European patients with osteosarcoma and 994 European individuals without cancer, with outliers removed for patient-control group comparisons. Patients were subsequently compared with individuals in the ExAC group. All data were analyzed from June 1, 2017, to July 1, 2019. Main Outcomes and Measures: The frequency of rare pathogenic or likely pathogenic genetic variants. Results: Among 1244 patients with osteosarcoma (mean [SD] age at diagnosis, 16 [8.9] years [range, 2-80 years]; 684 patients [55.0%] were male), an analysis restricted to individuals with European ancestry indicated a significantly higher pathogenic or likely pathogenic variant burden in 238 high-interest cancer-susceptibility genes among patients with osteosarcoma compared with the control group (732 vs 994, respectively; P = 1.3 × 10-18). A pathogenic or likely pathogenic cancer-susceptibility gene variant was identified in 281 of 1004 patients with osteosarcoma (28.0%), of which nearly three-quarters had a variant that mapped to an autosomal-dominant gene or a known osteosarcoma-associated cancer predisposition syndrome gene. The frequency of a pathogenic or likely pathogenic cancer-susceptibility gene variant was 128 of 1062 individuals (12.1%) in the control group and 2527 of 27 173 individuals (9.3%) in the ExAC group. A higher than expected frequency of pathogenic or likely pathogenic variants was observed in genes not previously linked to osteosarcoma (eg, CDKN2A, MEN1, VHL, POT1, APC, MSH2, and ATRX) and in the Li-Fraumeni syndrome-associated gene, TP53. Conclusions and Relevance: In this study, approximately one-fourth of patients with osteosarcoma unselected for family history had a highly penetrant germline mutation requiring additional follow-up analysis and possible genetic counseling with cascade testing.

Identifying, understanding, and correcting technical artifacts on the sex chromosomes in next-generation sequencing data.

BACKGROUND: Mammalian X and Y chromosomes share a common evolutionary origin and retain regions of high sequence similarity. Similar sequence content can confound the mapping of short next-generation sequencing reads to a reference genome. It is therefore possible that the presence of both sex chromosomes in a reference genome can cause technical artifacts in genomic data and affect downstream analyses and applications. Understanding this problem is critical for medical genomics and population genomic inference. RESULTS: Here, we characterize how sequence homology can affect analyses on the sex chromosomes and present XYalign, a new tool that (1) facilitates the inference of sex chromosome complement from next-generation sequencing data; (2) corrects erroneous read mapping on the sex chromosomes; and (3) tabulates and visualizes important metrics for quality control such as mapping quality, sequencing depth, and allele balance. We find that sequence homology affects read mapping on the sex chromosomes and this has downstream effects on variant calling. However, we show that XYalign can correct mismapping, resulting in more accurate variant calling. We also show how metrics output by XYalign can be used to identify XX and XY individuals across diverse sequencing experiments, including low- and high-coverage whole-genome sequencing, and exome sequencing. Finally, we discuss how the flexibility of the XYalign framework can be leveraged for other uses including the identification of aneuploidy on the autosomes. XYalign is available open source under the GNU General Public License (version 3). CONCLUSIONS: Sex chromsome sequence homology causes the mismapping of short reads, which in turn affects downstream analyses. XYalign provides a reproducible framework to correct mismapping and improve variant calling on the sex chromsomes.

Blood DNA methylation and breast cancer risk: a meta-analysis of four prospective cohort studies.

BACKGROUND: Environmental and genetic factors play an important role in the etiology of breast cancer. Several small blood-based DNA methylation studies have reported risk associations with methylation at individual CpGs and average methylation levels; however, these findings require validation in larger prospective cohort studies. To investigate the role of blood DNA methylation on breast cancer risk, we conducted a meta-analysis of four prospective cohort studies, including a total of 1663 incident cases and 1885 controls, the largest study of blood DNA methylation and breast cancer risk to date. METHODS: We assessed associations with methylation at 365,145 CpGs present in the HumanMethylation450 (HM450K) Beadchip, after excluding CpGs that did not pass quality controls in all studies. Each of the four cohorts estimated odds ratios (ORs) and 95% confidence intervals (CI) for the association between each individual CpG and breast cancer risk. In addition, each study assessed the association between average methylation measures and breast cancer risk, adjusted and unadjusted for cell-type composition. Study-specific ORs were combined using fixed-effect meta-analysis with inverse variance weights. Stratified analyses were conducted by age at diagnosis (< 50, ≥ 50), estrogen receptor (ER) status (+/-), and time since blood collection (< 5, 5-10, > 10 years). The false discovery rate (q value) was used to account for multiple testing. RESULTS: The average age at blood draw ranged from 52.2 to 62.2 years across the four cohorts. Median follow-up time ranged from 6.6 to 8.4 years. The methylation measured at individual CpGs was not associated with breast cancer risk (q value > 0.59). In addition, higher average methylation level was not associated with risk of breast cancer (OR = 0.94, 95% CI = 0.85, 1.05; P = 0.26; P for study heterogeneity = 0.86). We found no evidence of modification of this association by age at diagnosis (P = 0.17), ER status (P = 0.88), time since blood collection (P = 0.98), or CpG location (P = 0.98). CONCLUSIONS: Our data indicate that DNA methylation measured in the blood prior to breast cancer diagnosis in predominantly postmenopausal women is unlikely to be associated with substantial breast cancer risk on the HM450K array. Larger studies or with greater methylation coverage are needed to determine if associations exist between blood DNA methylation and breast cancer risk.

Genome-wide association study identifies multiple new loci associated with Ewing sarcoma susceptibility.

Ewing sarcoma (EWS) is a pediatric cancer characterized by the EWSR1-FLI1 fusion. We performed a genome-wide association study of 733 EWS cases and 1346 unaffected individuals of European ancestry. Our study replicates previously reported susceptibility loci at 1p36.22, 10q21.3 and 15q15.1, and identifies new loci at 6p25.1, 20p11.22 and 20p11.23. Effect estimates exhibit odds ratios in excess of 1.7, which is high for cancer GWAS, and striking in light of the rarity of EWS cases in familial cancer syndromes. Expression quantitative trait locus (eQTL) analyses identify candidate genes at 6p25.1 (RREB1) and 20p11.23 (KIZ). The 20p11.22 locus is near NKX2-2, a highly overexpressed gene in EWS. Interestingly, most loci reside near GGAA repeat sequences and may disrupt binding of the EWSR1-FLI1 fusion protein. The high locus to case discovery ratio from 733 EWS cases suggests a genetic architecture in which moderate risk SNPs constitute a significant fraction of risk.

Successful use of whole genome amplified DNA from multiple source types for high-density Illumina SNP microarrays.

BACKGROUND: The recommended genomic DNA input requirements for whole genome single nucleotide polymorphism microarrays can limit the scope of molecular epidemiological studies. We performed a large-scale evaluation of whole genome amplified DNA as input into high-density, whole-genome Illumina® Infinium® SNP microarray. RESULTS: Overall, 6622 DNA samples from 5970 individuals were obtained from three distinct biospecimen sources and genotyped using gDNA and/or wgaDNA inputs. When genotypes from the same individual were compared with standard, native gDNA input amount, we observed 99.94% mean concordance with wgaDNA input. CONCLUSIONS: Our results demonstrate that carefully conducted studies with wgaDNA inputs can yield high-quality genotyping results. These findings should enable investigators to consider expansion of ongoing studies using high-density SNP microarrays, currently challenged by small amounts of available DNA.

Genome-wide association study identifies the GLDC/IL33 locus associated with survival of osteosarcoma patients.

Survival rates for osteosarcoma, the most common primary bone cancer, have changed little over the past three decades and are particularly low for patients with metastatic disease. We conducted a multi-institutional genome-wide association study (GWAS) to identify germline genetic variants associated with overall survival in 632 patients with osteosarcoma, including 523 patients of European ancestry and 109 from Brazil. We conducted a time-to-event analysis and estimated hazard ratios (HR) and 95% confidence intervals (CI) using Cox proportional hazards models, with and without adjustment for metastatic disease. The results were combined across the European and Brazilian case sets using a random-effects meta-analysis. The strongest association after meta-analysis was for rs3765555 at 9p24.1, which was inversely associated with overall survival (HR = 1.76; 95% CI 1.41-2.18, p = 4.84 × 10-7 ). After imputation across this region, the combined analysis identified two SNPs that reached genome-wide significance. The strongest single association was with rs55933544 (HR = 1.9; 95% CI 1.5-2.4; p = 1.3 × 10-8 ), which localizes to the GLDC gene, adjacent to the IL33 gene and was consistent across both the European and Brazilian case sets. Using publicly available data, the risk allele was associated with lower expression of IL33 and low expression of IL33 was associated with poor survival in an independent set of patients with osteosarcoma. In conclusion, we have identified the GLDC/IL33 locus on chromosome 9p24.1 as associated with overall survival in patients with osteosarcoma. Further studies are needed to confirm this association and shed light on the biological underpinnings of this susceptibility locus.

Genome-Wide Association Study to Identify Susceptibility Loci That Modify Radiation-Related Risk for Breast Cancer After Childhood Cancer.

Background: Childhood cancer survivors treated with chest-directed radiotherapy have substantially elevated risk for developing breast cancer. Although genetic susceptibility to breast cancer in the general population is well studied, large-scale evaluation of breast cancer susceptibility after chest-directed radiotherapy for childhood cancer is lacking. Methods: We conducted a genome-wide association study of breast cancer in female survivors of childhood cancer, pooling two cohorts with detailed treatment data and systematic, long-term follow-up: the Childhood Cancer Survivor Study and St. Jude Lifetime Cohort. The study population comprised 207 survivors who developed breast cancer and 2774 who had not developed any subsequent neoplasm as of last follow-up. Genotyping and subsequent imputation yielded 16 958 466 high-quality variants for analysis. We tested associations in the overall population and in subgroups stratified by receipt of lower than 10 and 10 or higher gray breast radiation exposure. We report P values and pooled per-allele risk estimates from Cox proportional hazards regression models. All statistical tests were two-sided. Results: Among survivors who received 10 or higher gray breast radiation exposure, a locus on 1q41 was associated with subsequent breast cancer risk (rs4342822, nearest gene PROX1 , risk allele frequency in control subjects [RAF controls ] = 0.46, hazard ratio = 1.92, 95% confidence interval = 1.49 to 2.44, P = 7.09 × 10 -9 ). Two rare variants also showed potentially promising associations (breast radiation ≥10 gray: rs74949440, 11q23, TAGLN , RAF controls = 0.02, P = 5.84 × 10 -8 ; <10 gray: rs17020562, 1q32.3, RPS6KC1 , RAF controls = 0.0005, P = 6.68 × 10 -8 ). Associations were restricted to these dose subgroups, with consistent findings in the two survivor cohorts. Conclusions: Our study provides strong evidence that germline genetics outside high-risk syndromes could modify the effect of radiation exposure on breast cancer risk after childhood cancer.

Evaluating the Causal Link Between Malaria Infection and Endemic Burkitt Lymphoma in Northern Uganda: A Mendelian Randomization Study.

BACKGROUND: Plasmodium falciparum (Pf) malaria infection is suspected to cause endemic Burkitt Lymphoma (eBL), but the evidence remains unsettled. An inverse relationship between sickle cell trait (SCT) and eBL, which supports that between malaria and eBL, has been reported before, but in small studies with low power. We investigated this hypothesis in children in a population-based study in northern Uganda using Mendelian Randomization. METHODS: Malaria-related polymorphisms (SCT, IL10, IL1A, CD36, SEMA3C, and IFNAR1) were genotyped in 202 eBL cases and 624 controls enrolled during 2010-2015. We modeled associations between genotypes and eBL or malaria using logistic regression. FINDINGS: SCT was associated with decreased risk of eBL (adjusted odds ratio [OR] 0·37, 95% CI 0·21-0·66; p=0·0003). Decreased risk of eBL was associated with IL10 rs1800896-CT (OR 0·73, 95% CI 0·50-1·07) and -CC genotypes (OR 0·53, 95% CI 0·29-0·95, ptrend=0·019); IL1A rs2856838-AG (OR 0·56, 95% CI 0·39-0·81) and -AA genotype (OR 0·50, 95% CI 0·28-1·01, ptrend=0·0016); and SEMA3C rs4461841-CT or -CC genotypes (OR 0·57, 95% CI 0·35-0·93, p=0·0193). SCT and IL10 rs1800896, IL1A rs2856838, but not SEMA3C rs4461841, polymorphisms were associated with decreased risk of malaria in the controls. INTERPRETATION: Our results support a causal effect of malaria infection on eBL.

Higher-than-expected population prevalence of potentially pathogenic germline TP53 variants in individuals unselected for cancer history.

Li-Fraumeni syndrome (LFS) is an autosomal-dominant cancer predisposition disorder associated with pathogenic germline variants in TP53, with a high penetrance over an individual's lifetime. The actual population prevalence of pathogenic germline TP53 mutations is still unclear, most likely due to biased selection of cancer affected families. The aim of this study was to estimate the population prevalence of potentially pathogenic TP53 exonic variants in three sequencing databases, totaling 63,983 unrelated individuals. Potential pathogenicity was defined using an original algorithm combining bioinformatic prediction tools, suggested clinical significance, and functional data. We identified 34 different potentially pathogenic TP53 variants in 131 out of 63,983 individuals (0.2%). Twenty-eight (82%) of these variants fell within the DNA-binding domain of TP53, with an enrichment for specific variants that were not previously identified as LFS mutation hotspots, such as the p.R290H and p.N235S variants. Our findings reveal that the population prevalence of potentially pathogenic TP53 variants may be up to 10 times higher than previously estimated from family-based studies. These results point to the need for further studies aimed at evaluating cancer penetrance modifiers as well as the risk associated between cancer and rare TP53 variants.

Whole exome sequencing in 75 high-risk families with validation and replication in independent case-control studies identifies TANGO2, OR5H14, and CHAD as new prostate cancer susceptibility genes.

Prostate cancer (PCa) susceptibility is defined by a continuum from rare, high-penetrance to common, low-penetrance alleles. Research to date has concentrated on identification of variants at the ends of that continuum. Taking an alternate approach, we focused on the important but elusive class of low-frequency, moderately penetrant variants by performing disease model-based variant filtering of whole exome sequence data from 75 hereditary PCa families. Analysis of 341 candidate risk variants identified nine variants significantly associated with increased PCa risk in a population-based, case-control study of 2,495 men. In an independent nested case-control study of 7,121 men, there was risk association evidence for TANGO2 p.Ser17Ter and the established HOXB13 p.Gly84Glu variant. Meta-analysis combining the case-control studies identified two additional variants suggestively associated with risk, OR5H14 p.Met59Val and CHAD p.Ala342Asp. The TANGO2 and HOXB13 variants co-occurred in cases more often than expected by chance and never in controls. Finally, TANGO2 p.Ser17Ter was associated with aggressive disease in both case-control studies separately. Our analyses identified three new PCa susceptibility alleles in the TANGO2, OR5H14 and CHAD genes that not only segregate in multiple high-risk families but are also of importance in altering disease risk for men from the general population. This is the first successful study to utilize sequencing in high-risk families for the express purpose of identifying low-frequency, moderately penetrant PCa risk mutations.

Female chromosome X mosaicism is age-related and preferentially affects the inactivated X chromosome.

To investigate large structural clonal mosaicism of chromosome X, we analysed the SNP microarray intensity data of 38,303 women from cancer genome-wide association studies (20,878 cases and 17,425 controls) and detected 124 mosaic X events >2 Mb in 97 (0.25%) women. Here we show rates for X-chromosome mosaicism are four times higher than mean autosomal rates; X mosaic events more often include the entire chromosome and participants with X events more likely harbour autosomal mosaic events. X mosaicism frequency increases with age (0.11% in 50-year olds; 0.45% in 75-year olds), as reported for Y and autosomes. Methylation array analyses of 33 women with X mosaicism indicate events preferentially involve the inactive X chromosome. Our results provide further evidence that the sex chromosomes undergo mosaic events more frequently than autosomes, which could have implications for understanding the underlying mechanisms of mosaic events and their possible contribution to risk for chronic diseases.