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.

Massively parallel reporter assays and variant scoring identified functional variants and target genes for melanoma loci and highlighted cell-type specificity.

The most recent genome-wide association study (GWAS) of cutaneous melanoma identified 54 risk-associated loci, but functional variants and their target genes for most have not been established. Here, we performed massively parallel reporter assays (MPRAs) by using malignant melanoma and normal melanocyte cells and further integrated multi-layer annotation to systematically prioritize functional variants and susceptibility genes from these GWAS loci. Of 1,992 risk-associated variants tested in MPRAs, we identified 285 from 42 loci (78% of the known loci) displaying significant allelic transcriptional activities in either cell type (FDR < 1%). We further characterized MPRA-significant variants by motif prediction, epigenomic annotation, and statistical/functional fine-mapping to create integrative variant scores, which prioritized one to six plausible candidate variants per locus for the 42 loci and nominated a single variant for 43% of these loci. Overlaying the MPRA-significant variants with genome-wide significant expression or methylation quantitative trait loci (eQTLs or meQTLs, respectively) from melanocytes or melanomas identified candidate susceptibility genes for 60% of variants (172 of 285 variants). CRISPRi of top-scoring variants validated their cis-regulatory effect on the eQTL target genes, MAFF (22q13.1) and GPRC5A (12p13.1). Finally, we identified 36 melanoma-specific and 45 melanocyte-specific MPRA-significant variants, a subset of which are linked to cell-type-specific target genes. Analyses of transcription factor availability in MPRA datasets and variant-transcription-factor interaction in eQTL datasets highlighted the roles of transcription factors in cell-type-specific variant functionality. In conclusion, MPRAs along with variant scoring effectively prioritized plausible candidates for most melanoma GWAS loci and highlighted cellular contexts where the susceptibility variants are functional.

Cell-type-specific meQTLs extend melanoma GWAS annotation beyond eQTLs and inform melanocyte gene-regulatory mechanisms.

Although expression quantitative trait loci (eQTLs) have been powerful in identifying susceptibility genes from genome-wide association study (GWAS) findings, most trait-associated loci are not explained by eQTLs alone. Alternative QTLs, including DNA methylation QTLs (meQTLs), are emerging, but cell-type-specific meQTLs using cells of disease origin have been lacking. Here, we established an meQTL dataset by using primary melanocytes from 106 individuals and identified 1,497,502 significant cis-meQTLs. Multi-QTL colocalization with meQTLs, eQTLs, and mRNA splice-junction QTLs from the same individuals together with imputed methylome-wide and transcriptome-wide association studies identified candidate susceptibility genes at 63% of melanoma GWAS loci. Among the three molecular QTLs, meQTLs were the single largest contributor. To compare melanocyte meQTLs with those from malignant melanomas, we performed meQTL analysis on skin cutaneous melanomas from The Cancer Genome Atlas (n = 444). A substantial proportion of meQTL probes (45.9%) in primary melanocytes is preserved in melanomas, while a smaller fraction of eQTL genes is preserved (12.7%). Integration of melanocyte multi-QTLs and melanoma meQTLs identified candidate susceptibility genes at 72% of melanoma GWAS loci. Beyond GWAS annotation, meQTL-eQTL colocalization in melanocytes suggested that 841 unique genes potentially share a causal variant with a nearby methylation probe in melanocytes. Finally, melanocyte trans-meQTLs identified a hotspot for rs12203592, a cis-eQTL of a transcription factor, IRF4, with 131 candidate target CpGs. Motif enrichment and IRF4 ChIP-seq analysis demonstrated that these target CpGs are enriched in IRF4 binding sites, suggesting an IRF4-mediated regulatory network. Our study highlights the utility of cell-type-specific meQTLs.

A UVB-responsive common variant at chromosome band 7p21.1 confers tanning response and melanoma risk via regulation of the aryl hydrocarbon receptor, AHR.

Genome-wide association studies (GWASs) have identified a melanoma-associated locus on chromosome band 7p21.1 with rs117132860 as the lead SNP and a secondary independent signal marked by rs73069846. rs117132860 is also associated with tanning ability and cutaneous squamous cell carcinoma (cSCC). Because ultraviolet radiation (UVR) is a key environmental exposure for all three traits, we investigated the mechanisms by which this locus contributes to melanoma risk, focusing on cellular response to UVR. Fine-mapping of melanoma GWASs identified four independent sets of candidate causal variants. A GWAS region-focused Capture-C study of primary melanocytes identified physical interactions between two causal sets and the promoter of the aryl hydrocarbon receptor (AHR). Subsequent chromatin state annotation, eQTL, and luciferase assays identified rs117132860 as a functional variant and reinforced AHR as a likely causal gene. Because AHR plays critical roles in cellular response to dioxin and UVR, we explored links between this SNP and AHR expression after both 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) and ultraviolet B (UVB) exposure. Allele-specific AHR binding to rs117132860-G was enhanced following both, consistent with predicted weakened AHR binding to the risk/poor-tanning rs117132860-A allele, and allele-preferential AHR expression driven from the protective rs117132860-G allele was observed following UVB exposure. Small deletions surrounding rs117132860 introduced via CRISPR abrogates AHR binding, reduces melanocyte cell growth, and prolongs growth arrest following UVB exposure. These data suggest AHR is a melanoma susceptibility gene at the 7p21.1 risk locus and rs117132860 is a functional variant within a UVB-responsive element, leading to allelic AHR expression and altering melanocyte growth phenotypes upon exposure.

Genetic testing in severe aplastic anemia is required for optimal hematopoietic cell transplant outcomes.

Patients with severe aplastic anemia (SAA) can have an unrecognized inherited bone marrow failure syndrome (IBMFS) because of phenotypic heterogeneity. We curated germline genetic variants in 104 IBMFS-associated genes from exome sequencing performed on 732 patients who underwent hematopoietic cell transplant (HCT) between 1989 and 2015 for acquired SAA. Patients with pathogenic or likely pathogenic (P/LP) variants fitting known disease zygosity patterns were deemed unrecognized IBMFS. Carriers were defined as patients with a single P/LP variant in an autosomal recessive gene or females with an X-linked recessive P/LP variant. Cox proportional hazard models were used for survival analysis with follow-up until 2017. We identified 113 P/LP single-nucleotide variants or small insertions/deletions and 10 copy number variants across 42 genes in 121 patients. Ninety-one patients had 105 in silico predicted deleterious variants of uncertain significance (dVUS). Forty-eight patients (6.6%) had an unrecognized IBMFS (33% adults), and 73 (10%) were carriers. No survival difference between dVUS and acquired SAA was noted. Compared with acquired SAA (no P/LP variants), patients with unrecognized IBMFS, but not carriers, had worse survival after HCT (IBMFS hazard ratio [HR], 2.13; 95% confidence interval[CI], 1.40-3.24; P = .0004; carriers HR, 0.96; 95% CI, 0.62-1.50; P = .86). Results were similar in analyses restricted to patients receiving reduced-intensity conditioning (n = 448; HR IBMFS = 2.39; P = .01). The excess mortality risk in unrecognized IBMFS attributed to death from organ failure (HR = 4.88; P < .0001). Genetic testing should be part of the diagnostic evaluation for all patients with SAA to tailor therapeutic regimens. Carriers of a pathogenic variant in an IBMFS gene can follow HCT regimens for acquired SAA.

Joint associations between established genetic susceptibility loci, pesticide exposures, and risk of prostate cancer.

More than 200 genetic variants have been independently associated with prostate cancer risk. Studies among farmers have also observed increased prostate cancer risk associated with exposure to specific organophosphate (fonofos, terbufos, malathion, dimethoate) and organochlorine (aldrin, chlordane) insecticides. We examined the joint associations between these pesticides, established prostate cancer loci, and prostate cancer risk among 1,162 cases (588 aggressive) and 2,206 frequency-matched controls nested in the Agricultural Health Study cohort. History of lifetime pesticide use was combined with a polygenic risk score (PRS) generated using 256 established prostate cancer risk variants. Logistic regression models estimated the joint associations of the pesticides, the PRS, and the 256 individual genetic variants with risk of total and aggressive prostate cancer. Likelihood ratio tests assessed multiplicative interaction. We observed interaction between ever use of fonofos and the PRS in relation to total and aggressive prostate cancer risk. Compared to the reference group (never use, PRS < median), men with ever use of fonofos and PRS > median had elevated risks of total (OR 1.35 [1.06-1.73], p-interaction = 0.03) and aggressive (OR 1.49 [1.09-2.04], p-interaction = 0.19) prostate cancer. There was also suggestion of interaction between pesticides and individual genetic variants occurring in regions associated with DNA damage response (CDH3, EMSY genes) and with variants related to altered androgen receptor-driven transcriptional programs critical for prostate cancer. Our study provides evidence that men with greater genetic susceptibility to prostate cancer may be at higher risk if they are also exposed to pesticides and suggests potential mechanisms by which pesticides may increase prostate cancer risk.

Next-generation sequencing errors due to genetic variation in WRAP53 encoding TCAB1 on chromosome 17.

Next-generation sequencing (NGS) is a valuable tool, but has limitations in sequencing through repetitive runs of single nucleotides (homopolymers). Pathogenic germline variants in WRAP53 encoding telomere Cajal body protein 1 (TCAB1) are a known cause of dyskeratosis congenita. We identified a significant NGS error in WRAP53, c.1562dup, p.Ala522Glyfs*8 (rs755116516 G>-/GG/GGG) that did not validate by Sanger sequencing. This error occurs because rs755116516 G>-/GG/GGG (Chr17:7,606,714) is polymorphic, and variants at this site challenge the ability of NGS to accurately call the correct number of nucleotides in a homopolymer run. This was further complicated by the fact that chr17:7,606,721 (rs769202794) is multiallelic G>A, C, T, and that chr17:7,606,722 is also multiallelic (rs7640C>A/G/T and rs373064567C>delC). In addition to the expert interpretation of potentially clinically actionable variants, it recommended that all variants in regions of the genome with homopolymers be validated by Sanger sequencing before clinical action.

Rare germline variants in PALB2 and BRCA2 in familial and sporadic chordoma.

Chordoma is a rare bone tumor with genetic risk factors largely unknown. We conducted a whole-exome sequencing (WES) analysis of germline DNA from 19 familial chordoma cases in five pedigrees and 137 sporadic chordoma patients and identified 17 rare germline variants in PALB2 and BRCA2, whose products play essential roles in homologous recombination (HR) and tumor suppression. One PALB2 variant showed disease cosegregation in a family with four affected people or obligate gene carrier. Chordoma cases had a significantly increased burden of rare variants in both genes when compared to population-based controls. Four of the six PALB2 variants identified from chordoma patients modestly affected HR function and three of the 11 BRCA2 variants caused loss of function in experimental assays. These results, together with previous reports of abnormal morphology and Brachyury expression of the notochord in Palb2 knockout mouse embryos and genomic signatures associated with HR defect and HR gene mutations in advanced chordomas, suggest that germline mutations in PALB2 and BRCA2 may increase chordoma susceptibility. Our data shed light on the etiology of chordoma and support the previous finding that PARP-1 inhibitors may be a potential therapy for some chordoma patients.

The oral microbiome and breast cancer and nonmalignant breast disease, and its relationship with the fecal microbiome in the Ghana Breast Health Study.

The oral microbiome, like the fecal microbiome, may be related to breast cancer risk. Therefore, we investigated whether the oral microbiome was associated with breast cancer and nonmalignant breast disease, and its relationship with the fecal microbiome in a case-control study in Ghana. A total of 881 women were included (369 breast cancers, 93 nonmalignant cases and 419 population-based controls). The V4 region of the 16S rRNA gene was sequenced from oral and fecal samples. Alpha-diversity (observed amplicon sequence variants [ASVs], Shannon index and Faith's Phylogenetic Diversity) and beta-diversity (Bray-Curtis, Jaccard and weighted and unweighted UniFrac) metrics were computed. MiRKAT and logistic regression models were used to investigate the case-control associations. Oral sample alpha-diversity was inversely associated with breast cancer and nonmalignant breast disease with odds ratios (95% CIs) per every 10 observed ASVs of 0.86 (0.83-0.89) and 0.79 (0.73-0.85), respectively, compared to controls. Beta-diversity was also associated with breast cancer and nonmalignant breast disease compared to controls (P ≤ .001). The relative abundances of Porphyromonas and Fusobacterium were lower for breast cancer cases compared to controls. Alpha-diversity and presence/relative abundance of specific genera from the oral and fecal microbiome were strongly correlated among breast cancer cases, but weakly correlated among controls. Particularly, the relative abundance of oral Porphyromonas was strongly, inversely correlated with fecal Bacteroides among breast cancer cases (r = -.37, P ≤ .001). Many oral microbial metrics were strongly associated with breast cancer and nonmalignant breast disease, and strongly correlated with fecal microbiome among breast cancer cases, but not controls.

Shwachman Diamond syndrome: narrow genotypic spectrum and variable clinical features.

BACKGROUND AND OBJECTIVES: Shwachman Diamond syndrome (SDS) is an inherited bone marrow failure syndrome (IBMFS) associated with pancreatic insufficiency, neutropenia, and skeletal dysplasia. Biallelic pathogenic variants (PV) in SBDS account for >90% of SDS. We hypothesized that the SDS phenotype varies based on genotype and conducted a genotype-phenotype correlation study to better understand these complexities. METHODS: We reviewed records of all patients with SDS or SDS-like syndromes in the National Cancer Institute's (NCI) IBMFS study. Additional published SDS cohorts were reviewed and compared with the NCI cohort. RESULTS: PVs in SBDS were present in 32/47 (68.1%) participants. Biallelic inheritance of SBDS c.258 + 2T > C and c.183_184TA > CT was the most common genotype in our study (25/32, 78.1%) and published cohorts. Most patients had the SDS hallmark features of neutropenia (45/45, 100%), pancreatic insufficiency (41/43, 95.3%), and/or bony abnormalities (29/36, 80.6%). Developmental delay was common (20/34, 58.8%). Increased risk of hematologic malignancies at young ages and the rarity of solid malignancies was observed in both the NCI cohort and published studies. CONCLUSIONS: SDS is a complex childhood illness with a narrow genotypic spectrum. Patients may first present to primary care, gastroenterology, orthopedic, and/or hematology clinics. Coordinated multidisciplinary care is important for diagnosis and patient management. CLINICAL TRIAL REGISTRATION: ClinicalTrials.gov Identifier: NCT00027274. IMPACT: The clinical and genetic spectrum of Shwachman Diamond Syndrome was comprehensively evaluated, and the findings illustrate the importance of a multidisciplinary approach for these complex patients. Our work reveals: 1. a narrow genotypic spectrum in SDS; 2. a low risk of solid tumors in patients with SDS; 3. patients with SDS have clinical manifestations in multiple organ systems.

Aflatoxin levels and prevalence of TP53 aflatoxin-mutations in hepatocellular carcinomas in Mexico.

OBJECTIVE: To determine the exposure to aflatoxin B1 (AFB1) in southern Mexico and the presence of the aflatoxin signature mutation in hepatocellular carcinoma (HCC) tissue from patients from a cancer referral center. MATERIALS AND METHODS: We estimated the prevalence and distribution of AFB1 in a representative sample of 100 women and men from Chiapas using the National Health and Nutrition Survey 2018-19. We also examined the presence of the aflatoxin signature mutation in codon 249 (R249S), and other relevant mutations of the TP53 gene in HCC tissue blocks from 24 women and 26 men treated in a national cancer referral center. RESULTS: The prevalence of AFB1 in serum samples was 85.5% (95%CI 72.1-93.1) and the median AFB1 was 0.117 pg/µL (IQR, 0.050-0.350). We detected TP53 R249S in three of the 50 HCCs (6.0%) and observed four other G>T transversions potentially induced by AFB1. CONCLUSION: Our analysis provides evidence that AFB1 may have a relevant role on HCC etiology in Mexico.

Population Frequency of Fanconi Pathway Gene Variants and Their Association with Survival After Hematopoietic Cell Transplantation for Severe Aplastic Anemia.

Severe aplastic anemia (SAA) is most frequently immune-mediated; however, rare inherited bone marrow failure syndromes, such as Fanconi anemia (FA), may be causal and can present as aplastic anemia (AA). FA is primarily an autosomal recessive disorder caused by the presence of 2 pathogenic variants in a single FA/BRCA DNA repair pathway gene. Patients with SAA often undergo genetic testing during clinical evaluation that may identify single deleterious alleles in FA pathway genes. We quantified the rate of germline single deleterious alleles in 22 FA genes using both a general population database (3234 variants, 125,748 exomes) and in a cohort of patients with SAA undergoing hematopoietic cell transplantation (HCT) (21 variants in 730 patients). The variants were classified as deleterious using in silico tools (REVEL, MetaSVM, VEP) and database resources (ClinVar, LOVD-FA). We found similar rates of single deleterious alleles in FA genes in both groups (2.6% and 2.9%). The presence of a single deleterious variant in a gene for FA in SAA HCT recipients did not affect the overall survival after HCT (hazard ratio, 0.85; 95% CI, 0.37 to 1.95; P  = 0.71), or post-HCT cancer risk (P = 0.52). Our results demonstrate that the identification of a germline monoallelic deleterious variant in an FA gene in patients with idiopathic SAA does not influence the outcome of HCT. Our findings suggest that there is no need for special treatment considerations for patients with SAA and a single deleterious FA allele identified on genetic testing.

Pre-transplant short telomeres are associated with high mortality risk after unrelated donor haematopoietic cell transplant for severe aplastic anaemia.

Telomeres are essential for chromosomal stability and markers of biological age. We evaluated the effect of pre-transplant short (<10th percentile-for-age) or very short (<5th or <1st percentile-for-age) leucocyte telomere length on survival after unrelated donor haematopoietic cell transplantation (HCT) for acquired severe aplastic anaemia (SAA). Patient pre-transplant blood samples and clinical data were available at the Center for International Blood and Marrow Transplant Research. We used quantitative real time polymerase chain reaction to measure relative telomere length (RTL) in 490 SAA patients who received HCT between 1990 and 2013 (median age = 20 years). One hundred and twelve patients (22·86%) had pre-HCT RTL <10th percentile-for-age, with the majority below the 5th percentile (N = 80, 71·43%). RTL <10th percentile-for-age was associated with a higher risk of post-HCT mortality (hazard ratio [HR] = 1·78, 95% confidence interval [CI]=1·18-2·69, P = 0·006) compared with RTL ≥50th percentile; no survival differences were noted in longer RTL categories (P > 0·10). Time-dependent effects for post-HCT mortality were only observed in relation to very short RTL; HR comparing RTL <5th versus ≥5th percentile = 1·38, P = 0·15 for the first 12 months after HCT, and HR = 3·91, P < 0·0001, thereafter, P-heterogeneity = 0·008; the corresponding HRs for RTL <1st versus ≥1st percentile = 1·29, P = 0·41, and HR = 5·18, P < 0·0001, P-heterogeneity = 0·005. The study suggests a potential role for telomere length in risk stratification of SAA patients in regard to their HCT survival.

Immune gene expression profiling reveals heterogeneity in luminal breast tumors.

BACKGROUND: Heterogeneity of immune gene expression patterns of luminal breast cancer (BC), which is clinically heterogeneous and overall considered as low immunogenic, has not been well studied especially in non-European populations. Here, we aimed at characterizing the immune gene expression profile of luminal BC in an Asian population and associating it with patient characteristics and tumor genomic features. METHODS: We performed immune gene expression profiling of tumor and adjacent normal tissue in 92 luminal BC patients from Hong Kong using RNA-sequencing data and used unsupervised consensus clustering to stratify tumors. We then used luminal patients from The Cancer Genome Atlas (TCGA, N = 564) and a Korean breast cancer study (KBC, N = 112) as replication datasets. RESULTS: Based on the expression of 130 immune-related genes, luminal tumors were stratified into three distinct immune subtypes. Tumors in one subtype showed higher level of tumor-infiltrating lymphocytes (TILs), characterized by T cell gene activation, higher expression of immune checkpoint genes, higher nonsynonymous mutation burden, and higher APOBEC-signature mutations, compared with other luminal tumors. The high-TIL subtype was also associated with lower ESR1/ESR2 expression ratio and increasing body mass index. The comparison of the immune profile in tumor and matched normal tissue suggested a tumor-derived activation of specific immune responses, which was only seen in high-TIL patients. Tumors in a second subtype were characterized by increased expression of interferon-stimulated genes and enrichment for TP53 somatic mutations. The presence of three immune subtypes within luminal BC was replicated in TCGA and KBC, although the pattern was more similar in Asian populations. The germline APOBEC3B deletion polymorphism, which is prevalent in East Asian populations and was previously linked to immune activation, was not associated with immune subtypes in our study. This result does not support the hypothesis that the germline APOBEC3B deletion polymorphism is the driving force for immune activation in breast tumors in Asian populations. CONCLUSION: Our findings suggest that immune gene expression and associated genomic features could be useful to further stratify luminal BC beyond the current luminal A/B classification and a subset of luminal BC patients may benefit from checkpoint immunotherapy, at least in Asian populations.

Rare germline variants in known melanoma susceptibility genes in familial melanoma.

Known high-risk cutaneous malignant melanoma (CMM) genes account for melanoma risk in <40% of melanoma-prone families, suggesting the existence of additional high-risk genes or perhaps a polygenic mechanism involving multiple genetic modifiers. The goal of this study was to systematically characterize rare germline variants in 42 established melanoma genes among 144 CMM patients in 76 American CMM families without known mutations using data from whole-exome sequencing. We identified 68 rare (<0.1% in public and in-house control datasets) nonsynonymous variants in 25 genes. We technically validated all loss-of-function, inframe insertion/deletion, and missense variants predicted as deleterious, and followed them up in 1, 559 population-based CMM cases and 1, 633 controls. Several of these variants showed disease co-segregation within families. Of particular interest, a stopgain variant in TYR was present in five of six CMM cases/obligate gene carriers in one family and a single population-based CMM case. A start gain variant in the 5'UTR region of PLA2G6 and a missense variant in ATM were each seen in all three affected people in a single family, respectively. Results from rare variant burden tests showed that familial and population-based CMM patients tended to have higher frequencies of rare germline variants in albinism genes such as TYR, TYRP1, and OCA2 (P < 0.05). Our results suggest that rare nonsynonymous variants in low- or intermediate-risk CMM genes may influence familial CMM predisposition, warranting further investigation of both common and rare variants in genes affecting functionally important pathways (such as melanogenesis) in melanoma risk assessment.

Combined somatic mutation and copy number analysis in the survival of familial CLL.

Recurrent large-scale somatic copy number alterations (SCNAs), and somatic point mutations can be analysed to stratify patients with chronic lymphocytic leukaemia (CLL) into distinct prognostic groups. To investigate the relationship between SCNAs and somatic mutations, we performed whole-exome sequencing and single nucleotide polymorphism microarray analyses on 98 CLL patients from 40 families with a high burden of CLL. Overall, 69 somatic mutations in 29 CLL driver genes were detected among 45 subjects (46%), with the most frequently mutated genes being TP53 (8·2%), NOTCH1 (8·2%) and ATM (5·1%). Additionally, 142 SCNAs from 54 subjects (57%) were detected, including losses of chromosome 13q14 (28·9%), 11q (5·6%), 17p (2·1%), and gain of chromosome 12 (4·2%). We found that patients having both an adverse point mutation in a CLL driver gene and an unfavourable SCNA tended to have poorer survival (Hazard ratio [HR] = 3·17, 95% confidence interval [CI] = 0·97-10·35; P = 0·056) than patients having either a point mutation (HR = 1·34, 95%CI = 0·66-2·71; P = 0·42) or SCNAs (HR = 2·65, 95%CI = 0·77-9·13; P = 0·12). TP53 mutation carriers were associated with the poorest overall survival (HR = 4·39, 95%CI = 1·28-15·04; P = 0·018). Our study suggests that combining SCNA and mutational data could contribute to predicting outcome in familial CLL.

Germline mutations in Protection of Telomeres 1 in two families with Hodgkin lymphoma.

In a previous whole exome sequencing of patients from 41 families with Hodgkin lymphoma, we identified two families with distinct heterozygous rare coding variants in POT1 (D224N and Y36H), both in a highly conserved region of the gene. POT1 D224N mutant did not bind to a single-stranded telomere oligonucleotide in vitro suggesting the mutation perturbs POT1's ability to bind to the telomeric G-rich overhang. Human HT1080 cells expressing POT1 D224N and lymphoblastoid cells carrying Y36H both showed increased telomere length and fragility in comparison to wild type cells. This strongly suggests that mutant POT1 causes chromosome instability and may play a role in lymphomagenesis in these families.

Association of Aflatoxin and Gallbladder Cancer.

BACKGROUND & AIMS: Aflatoxin, which causes hepatocellular carcinoma, may also cause gallbladder cancer. We investigated whether patients with gallbladder cancer have higher exposure to aflatoxin than patients with gallstones. METHODS: We measured aflatoxin B1 (AFB1)-lysine adducts in plasma samples from the Shanghai Biliary Tract Cancer case-control study, conducted from 1997 through 2001. We calculated age- and sex-adjusted odds ratios (ORs) and 95% confidence intervals (CIs) and the population-attributable fraction for 209 patients with gallbladder cancer and gallstones vs 250 patients with gallstones without cancer (controls). In 54 patients with gallbladder cancer, tumor tissue was examined for the R249S mutation in TP53, associated with aflatoxin exposure, through targeted sequencing. RESULTS: The AFB1-lysine adduct was detected in 67 (32%) of 209 patients with gallbladder cancer and 37 (15%) of the 250 controls (χ2 P < .0001), almost threefold more patients with gallbladder cancer than controls (OR, 2.71; 95% CI, 1.70-4.33). Among participants with detectable levels of AFB1-lysine, the median level of AFB1-lysine was 5.4 pg/mg in those with gallbladder cancer, compared with 1.2 pg/mg in controls. For patients in the fourth quartile of AFB1-lysine level vs the first quartile, the OR for gallbladder cancer was 7.61 (95% CI, 2.01-28.84). None of the 54 gallbladder tumors sequenced were found to have the R249S mutation in TP53. The population-attributable fraction for cancer related to aflatoxin was 20% (95% CI, 15%-25%). CONCLUSIONS: In a case-control study of patients with gallbladder cancer and gallstones vs patients with gallstones without cancer, we associated exposure to aflatoxin (based on plasma level of AFB1-lysine) with gallbladder cancer. Gallbladder cancer does not appear associate with the R249S mutation in TP53. If aflatoxin is a cause of gallbladder cancer, it may have accounted for up to 20% of the gallbladder cancers in Shanghai, China, during the study period, and could account for an even higher proportion in high-risk areas. If our findings are verified, reducing aflatoxin exposure might reduce the incidence of gallbladder cancer.

Novel and known ribosomal causes of Diamond-Blackfan anaemia identified through comprehensive genomic characterisation.

BACKGROUND: Diamond-Blackfan anaemia (DBA) is an inherited bone marrow failure syndrome (IBMFS) characterised by erythroid hypoplasia. It is associated with congenital anomalies and a high risk of developing specific cancers. DBA is caused predominantly by autosomal dominant pathogenic variants in at least 15 genes affecting ribosomal biogenesis and function. Two X-linked recessive genes have been identified. OBJECTIVES: We aim to identify the genetic aetiology of DBA. METHODS: Of 87 families with DBA enrolled in an institutional review board-approved cohort study (ClinicalTrials.gov Identifier:NCT00027274), 61 had genetic testing information available. Thirty-five families did not have a known genetic cause and thus underwent comprehensive genomic evaluation with whole exome sequencing, deletion and CNV analyses to identify their disease-associated pathogenic variant. Controls for functional studies were healthy mutation-negative individuals enrolled in the same study. RESULTS: Our analyses uncovered heterozygous pathogenic variants in two previously undescribed genes in two families. One family had a non-synonymous variant (p.K77N) in RPL35; the second family had a non-synonymous variant (p. L51S) in RPL18. Both of these variants result in pre-rRNA processing defects. We identified heterozygous pathogenic variants in previously known DBA genes in 16 of 35 families. Seventeen families who underwent genetic analyses are yet to have a genetic cause of disease identified. CONCLUSIONS: Overall, heterozygous pathogenic variants in ribosomal genes were identified in 44 of the 61 families (72%). De novo pathogenic variants were observed in 57% of patients with DBA. Ongoing studies of DBA genomics will be important to understand this complex disorder.

Relative Telomere Length before Hematopoietic Cell Transplantation and Outcome after Unrelated Donor Hematopoietic Cell Transplantation for Acute Leukemia.

Telomeres are tandem nucleotide repeats and a protein complex located at the end of the chromosomes maintaining genomic stability. Their potential as a predictive biomarker for outcomes after allogeneic hematopoietic cell transplant (HCT) in hematologic malignancies is still unclear. From the Center for International Blood and Marrow Transplant Research we randomly selected 536 acute leukemia patients from those who underwent myeloablative 8/8 HLA-matched unrelated donor HCT between 2005 and 2012 and who had an available pre-HCT blood sample in the repository. Relative telomere length (RTL) was measured by real-time quantitative PCR. We used Kaplan-Meier and competing risk estimators to calculate survival probability and cumulative incidence, respectively, across patient RTL tertiles. Cox proportional hazard regression was used for adjusted analyses. The study included 396 acute myeloid leukemia (AML) and 140 acute lymphoblastic leukemia (ALL) patients. Median age at HCT was 41 years (range, .5 to 66), and median follow-up for survivors was 5.1 years (range, .4 to 8.3). Significant inverse correlations between age and RTL were observed in patients with AML (r = -.44, P < .0001) and ALL (r = -.48, P < .0001). Patients with ALL had longer RTL than those with AML (.48 versus .43, respectively); the difference was not statistically significant after adjusting for patient age (P = .96). Pre-HCT RTL in acute leukemia patients was not statistically significantly associated with overall survival (HR for longest RTL compared with shortest, .91; 95% CI, .65 to 1.28), disease-free survival (HR, .90; 95% CI, .64 to 1.25), transplant-related mortality (HR, .97; 95% CI, .60 to 1.59), incidence of relapse (HR, .89; 95% CI, .56 to 1.40), neutrophil engraftment (HR, 1.06; 95% CI, .85 to 1.32), or grades II to IV acute graft-versus-host disease (HR, 1.11; 95% CI, .81 to 1.53), grades III-IV acute graft-versus-host disease (HR, .92; 95% CI, .54 to 1.59), and chronic graft-versus-host disease (HR, 1.10; 95% CI, .81 to 1.50). In this study, recipient pre-HCT RTL had no prognostic role in post-transplant outcomes in acute leukemia patients.