Molecular characterization of Richter syndrome identifies de novo diffuse large B-cell lymphomas with poor prognosis.

Richter syndrome (RS) is the transformation of chronic lymphocytic leukemia (CLL) into aggressive lymphoma, most commonly diffuse large B-cell lymphoma (DLBCL). We characterize 58 primary human RS samples by genome-wide DNA methylation and whole-transcriptome profiling. Our comprehensive approach determines RS DNA methylation profile and unravels a CLL epigenetic imprint, allowing CLL-RS clonal relationship assessment without the need of the initial CLL tumor DNA. DNA methylation- and transcriptomic-based classifiers were developed, and testing on landmark DLBCL datasets identifies a poor-prognosis, activated B-cell-like DLBCL subset in 111/1772 samples. The classification robustly identifies phenotypes very similar to RS with a specific genomic profile, accounting for 4.3-8.3% of de novo DLBCLs. In this work, RS multi-omics characterization determines oncogenic mechanisms, establishes a surrogate marker for CLL-RS clonal relationship, and provides a clinically relevant classifier for a subset of primary "RS-type DLBCL" with unfavorable prognosis.

Correlation between DNA Methylation and Cell Proliferation Identifies New Candidate Predictive Markers in Meningioma.

Meningiomas are the most common primary tumors of the central nervous system. Based on the 2021 WHO classification, they are classified into three grades reflecting recurrence risk and aggressiveness. However, the WHO's histopathological criteria defining these grades are somewhat subjective. Together with reliable immunohistochemical proliferation indices, other molecular markers such as those studied with genome-wide epigenetics promise to revamp the current prognostic classification. In this study, 48 meningiomas of various grades were randomly included and explored for DNA methylation with the Infinium MethylationEPIC microarray over 850k CpG sites. We conducted differential and correlative analyses on grade and several proliferation indices and markers, such as mitotic index and Ki-67 or MCM6 immunohistochemistry. We also set up Cox proportional hazard models for extensive associations between CpG methylation and survival. We identified loci highly correlated with cell growth and a targeted methylation signature of regulatory regions persistently associated with proliferation, grade, and survival. Candidate genes under the control of these regions include SMC4, ESRRG, PAX6, DOK7, VAV2, OTX1, and PCDHA-PCDHB-PCDHG, i.e., the protocadherin gene clusters. This study highlights the crucial role played by epigenetic mechanisms in shaping dysregulated cellular proliferation and provides potential biomarkers bearing prognostic and therapeutic value for the clinical management of meningioma.

Epimutations in both the TESK2 and MMACHC promoters in the Epi-cblC inherited disorder of intracellular metabolism of vitamin B12.

BACKGROUND: epi-cblC is a recently discovered inherited disorder of intracellular vitamin B12 metabolism associating hematological, neurological, and cardiometabolic outcomes. It is produced by an epimutation at the promoter common to CCDC163P and MMACHC, which results from an aberrant antisense transcription due to splicing mutations in the antisense PRDX1 gene neighboring MMACHC. We studied whether the aberrant transcription produced a second epimutation by encompassing the CpG island of the TESK2 gene neighboring CCDC163P. METHODS: We unraveled the methylome architecture of the CCDC163P-MMACHC CpG island (CpG:33) and the TESK2 CpG island (CpG:51) of 17 epi-cblC cases. We performed an integrative analysis of the DNA methylome profiling, transcriptome reconstruction of RNA-sequencing (RNA-seq), chromatin immunoprecipitation sequencing (ChIP-Seq) of histone H3, and transcription expression of MMACHC and TESK2. RESULTS: The PRDX1 splice mutations and activation of numerous cryptic splice sites produced antisense readthrough transcripts encompassing the bidirectional MMACHC/CCDC163P promoter and the TESK2 promoter, resulting in the silencing of both the MMACHC and TESK2 genes through the deposition of SETD2-dependent H3K36me3 marks and the generation of epimutations in the CpG islands of the two promoters. CONCLUSIONS: The antisense readthrough transcription of the mutated PRDX1 produces an epigenetic silencing of MMACHC and TESK2. We propose using the term 'epi-digenism' to define this epigenetic disorder that affects two genes. Epi-cblC is an entity that differs from cblC. Indeed, the PRDX1 and TESK2 altered expressions are observed in epi-cblC but not in cblC, suggesting further evaluating the potential consequences on cancer risk and spermatogenesis.

Next-generation sequencing and genotype association studies reveal the association of HLA-DRB3*02:02 with delayed hypersensitivity to penicillins.

BACKGROUND: Nonimmediate (delayed)-allergic reactions to penicillins are common and some of them can be life-threatening. The genetic factors influencing these reactions are unknown/poorly known/poorly understood. We assessed the genetic predictors of a delayed penicillin allergy that cover the HLA loci. METHODS: Using next-generation sequencing (NGS), we genotyped the MHC region in 24 patients with delayed hypersensitivity compared with 20 patients with documented immediate hypersensitivity to penicillins recruited in Italy. Subsequently, we analyzed in silico Illumina Immunochip genotyping data that covered the HLA loci in 98 Spanish patients with delayed hypersensitivity and 315 with immediate hypersensitivity compared to 1,308 controls. RESULTS: The two alleles DRB3*02:02:01:02 and DRB3*02:02:01:01 were reported in twenty cases with delayed reactions (83%) and ten cases with immediate reactions (50%), but not in the Allele Frequency Net Database. Bearing at least one of the two alleles increased the risk of delayed reactions compared to immediate reactions, with an OR of 8.88 (95% CI, 3.37-23.32; p < .0001). The haplotype (ACAA) from rs9268835, rs6923504, rs6903608, and rs9268838 genetic variants of the HLA-DRB3 genomic region was significantly associated with an increased risk of delayed hypersensitivity to penicillins (OR, 1.7; 95% CI: 1.06-1.92; p = .001), but not immediate hypersensitivity. CONCLUSION: We showed that the HLA-DRB3 locus is strongly associated with an increased risk of delayed penicillin hypersensitivity, at least in Southwestern Europe. The determination of HLA-DRB3*02:02 alleles in the risk management of severe delayed hypersensitivity to penicillins should be evaluated further in larger population samples of different origins.

Low-frequency Coding Variants Associated With Body Mass Index Affect the Success of Bariatric Surgery.

CONTEXT: A recent study identified 14 low-frequency coding variants associated with body mass index (BMI) in 718 734 individuals predominantly of European ancestry. OBJECTIVE: We investigated the association of 2 genetic scores (GS) with i) the risk of severe/morbid obesity, ii) BMI variation before weight-loss intervention, iii) BMI change in response to an 18-month lifestyle/behavioral intervention program, and iv) BMI change up to 24 months after bariatric surgery. METHODS: The 14 low-frequency coding variants were genotyped or sequenced in 342 French adults with severe/morbid obesity and 574 French adult controls from the general population. We built risk and protective GS based on 6 BMI-increasing and 5 BMI-decreasing low-frequency coding variants that were polymorphic in our study. RESULTS: While the risk GS was not associated with severe/morbid obesity status, BMI-decreasing low-frequency coding variants were significantly less frequent in patients with severe/morbid obesity than in French adults from the general population. Neither the risk nor the protective GS was associated with BMI before intervention in patients with severe/morbid obesity, nor did they affect BMI change in response to a lifestyle/behavioral modification program. The protective GS was associated with a greater BMI decrease following bariatric surgery. The risk and protective GS were associated with a higher and lower risk of BMI regain after bariatric surgery. CONCLUSION: Our data indicate that in populations of European descent, low-frequency coding variants associated with BMI in the general population also affect the outcomes of bariatric surgery in patients with severe/morbid obesity.

Ionizing radiations induce shared epigenomic signatures unraveling adaptive mechanisms of cancerous cell lines with or without methionine dependency.

BACKGROUND: Although radiation therapy represents a core cancer treatment modality, its efficacy is hampered by radioresistance. The effect of ionizing radiations (IRs) is well known regarding their ability to induce genetic alterations; however, their impact on the epigenome landscape in cancer, notably at the CpG dinucleotide resolution, remains to be further deciphered. In addition, no evidence is available regarding the effect of IRs on the DNA methylome profile according to the methionine dependency phenotype, which represents a hallmark of metabolic adaptation in cancer. METHODS: We used a case-control study design with a fractionated irradiation regimen on four cancerous cell lines representative of HCC (HepG2), melanoma (MeWo and MeWo-LC1, which exhibit opposed methionine dependency phenotypes), and glioblastoma (U251). We performed high-resolution genome-wide DNA methylome profiling using the MethylationEPIC BeadChip on baseline conditions, irradiated cell lines (cumulative dose of 10 Gy), and non-irradiated counterparts. We performed epigenome-wide association studies to assess the effect of IRs and methionine-dependency-oriented analysis by carrying out epigenome-wide conditional logistic regression. We looked for epigenome signatures at the locus and single-probe (CpG dinucleotide) levels and through enrichment analyses of gene ontologies (GO). The EpiMet project was registered under the ID#AAP-BMS_003_211. RESULTS: EWASs revealed shared GO annotation pathways associated with increased methylation signatures for several biological processes in response to IRs, including blood circulation, plasma membrane-bounded cell projection organization, cell projection organization, multicellular organismal process, developmental process, and animal organ morphogenesis. Epigenome-wide conditional logistic regression analysis on the methionine dependency phenotype highlighted several epigenome signatures related to cell cycle and division and responses to IR and ultraviolet light. CONCLUSIONS: IRs generated a variation in the methylation level of a high number of CpG probes with shared biological pathways, including those associated with cell cycle and division, responses to IRs, sustained angiogenesis, tissue invasion, and metastasis. These results provide insight on shared adaptive mechanisms of the epigenome in cancerous cell lines in response to IR. Future experiments should focus on the tryptic association between IRs, the initiation of a radioresistance phenotype, and their interaction with methionine dependency as a hallmark of metabolic adaptation in cancer.

Population and evolutionary genetics of the PAH locus to uncover overdominance and adaptive mechanisms in phenylketonuria: Results from a multiethnic study.

BACKGROUND: Phenylketonuria (PKU) is the most common inborn error of amino acid metabolism in Europe. The reasons underlying the high prevalence of heterozygous carriers are not clearly understood. We aimed to look for pathogenic PAH variant enrichment according to geographical areas and patients' ethnicity using a multiethnic nationwide cohort of patients with PKU in France. We subsequently appraised the population differentiation, balancing selection and the molecular evolutionary history of the PAH locus. METHODS: The French nationwide PKU study included patients who have been referred at the national level to the University Hospital of Nancy, and for whom a molecular diagnosis of phenylketonuria was made by Sanger sequencing. We performed enrichment analyses by comparing alternative allele frequencies using Fisher's exact test with Bonferroni adjustment. We estimated the amount of genetic differentiation among populations using Wright's fixation index (Fst). To estimate the molecular evolutionary history of the PAH gene, we performed phylogenetic and evolutionary analyses using whole-genome and exome-sequencing data from healthy individuals and non-PKU patients, respectively. Finally, we used exome-wide association study to decipher potential genetic loci associated with population divergence on PAH. FINDINGS: The study included 696 patients and revealed 132 pathogenic PAH variants. Three geographical areas showed significant enrichment for a pathogenic PAH variant: North of France (p.Arg243Leu), North-West of France (p.Leu348Val), and Mediterranean coast (p.Ala403Val). One PAH variant (p.Glu280Gln) was significantly enriched among North-Africans (OR = 23·23; 95% CI: 9·75-55·38). PAH variants exhibiting a strong genetic differentiation were significantly enriched in the 'Biopterin_H' domain (OR = 6·45; 95% CI: 1·99-20·84), suggesting a balancing selection pressure on the biopterin function of PAH. Phylogenetic and timetree analyses were consistent with population differentiation events on European-, African-, and Asian-ancestry populations. The five PAH variants most strongly associated with a high selection pressure were phylogenetically close and were located within the biopterin domain coding region of PAH or in its vicinity. Among the non-PAH loci potentially associated with population divergence, two reached exome-wide significance: SSPO (SCO-spondin) and DBH (dopamine beta-hydroxylase), involved in neuroprotection and metabolic adaptation, respectively. INTERPRETATION: Our data provide evidence on the combination of evolutionary and adaptive events in populations with distinct ancestries, which may explain the overdominance of some genetic variants on PAH. FUNDING: French National Institute of Health and Medical Research (INSERM) UMR_S 1256.

Folate can promote the methionine-dependent reprogramming of glioblastoma cells towards pluripotency.

Methionine dependency of tumor growth, although not well-understood, is detectable by 11C-methionine positron emission tomography and may contribute to the aggressivity of glioblastomas (GBM) and meningiomas. Cytosolic folate cycle is required for methionine synthesis. Its dysregulation may influence cell reprogramming towards pluripotency. We evaluated methionine-dependent growth of monolayer (ML) cells and stem cell-like tumor spheres (TS) derived from 4 GBM (U251, U87, LN299, T98G) and 1 meningioma (IOMM-LEE) cell lines. Our data showed that for all cell lines studied, exogenous methionine is required for TS formation but not for ML cells proliferation. Furthermore, for GBM cell lines, regardless of the addition of folate cycle substrates (folic acid and formate), the level of 3 folate isoforms, 5-methytetrahydrofolate, 5,10-methenyltetrahydrofolate, and 10-formyltetrahydrofolate, were all downregulated in TS relative to ML cells. Unlike GBM cell lines, in IOMM-LEE cells, 5-methyltetrahydrofolate was actually more elevated in TS than ML, and only 5,10-methenyltetrahydrofolate and 10-formyltetrahydrofolate were downregulated. The functional significance of this variation in folate cycle repression was revealed by the finding that Folic Acid and 5-methyltetrahydrofolate promote the growth of U251 TS but not IOMM-LEE TS. Transcriptome-wide sequencing of U251 cells revealed that DHFR, SHMT1, and MTHFD1 were downregulated in TS vs ML, in concordance with the low activity cytosolic folate cycle observed in U251 TS. In conclusion, we found that a repressed cytosolic folate cycle underlies the methionine dependency of GBM and meningioma cell lines and that 5-methyltetrahydrofolate is a key metabolic switch for glioblastoma TS formation. The finding that folic acid facilitates TS formation, although requiring further validation in diseased human tissues, incites to investigate whether excessive folate intake could promote cancer stem cells formation in GBM patients.

BRIP1 coding variants are associated with a high risk of hepatocellular carcinoma occurrence in patients with HCV- or HBV-related liver disease.

The molecular mechanisms of hepatocellular carcinoma (HCC) carcinogenesis are still not fully understood. DNA repair defects may influence HCC risk. The aim of the study was to look for potential genetic variants of DNA repair genes associated with HCC risk among patients with alcohol- or viral-induced liver disease. We performed four case-control studies on 2,006 European- (Derivation#1 and #2 studies) and African-ancestry (Validation#1 and #2 studies) patients originating from several cohorts in order to assess the association between genetic variants on DNA repair genes and HCC risk using a custom array encompassing 94 genes. In the Derivation#1 study, the BRIP1 locus reached array-wide significance (Chi-squared SV-Perm, P=5.00×10-4) among the 253 haplotype blocks tested for their association with HCC risk, in patients with viral cirrhosis but not among those with alcoholic cirrhosis. The BRIP1 haplotype block included three exonic variants (rs4986763, rs4986764, rs4986765). The BRIP1 'AAA' haplotype was significantly associated with an increased HCC risk [odds ratio (OR), 2.01 (1.19-3.39); false discovery rate (FDR)-P=1.31×10-2]. In the Derivation#2 study, results were confirmed for the BRIP1 'GGG' haplotype [OR, 0.53 (0.36-0.79); FDR-P=3.90×10-3]. In both Validation#1 and #2 studies, BRIP1 'AAA' haplotype was significantly associated with an increased risk of HCC [OR, 1.71 (1.09-2.68); FDR-P=7.30×10-2; and OR, 6.45 (4.17-9.99); FDR-P=2.33×10-19, respectively]. Association between the BRIP1 locus and HCC risk suggests that impaired DNA mismatch repair might play a role in liver carcinogenesis, among patients with HCV- or HBV-related liver disease.

Cystathionine ß-synthase genetic variant rs2124459 is associated with a reduced risk of cleft palate in French and Belgian populations.

BACKGROUND: Orofacial cleft (OFC) is the most prevalent craniofacial birth defect. Genes involved in one-carbon, folate and vitamin B12 metabolisms have been associated with OFC but no study performed a concomitant assessment on genes involved in these three pathways. OBJECTIVE: We looked for potential genetic variants associated with OFC using an exhaustive gene panel of one-carbon metabolism. METHODS: We performed a case-control discovery study on children with OFC (236 cases, 145 controls) and their related mothers (186 cases, 127 controls). We performed a replication study on the top significant genetic variant in an independent group from Belgium (248 cases, 225 controls). RESULTS: In the discovery study on 'mothers', the CBS locus reached array-wide significance (p=9.13×10-6; Bonferroni p=4.77×10-3; OR 0.47 (0.33 to 0.66)) among the 519 haplotypes tested for their association with OFC risk. Within the CBS haplotype block (rs2124459, rs6586282, rs4920037, rs234705, rs234709), the rs2124459 was the most significantly associated with a reduced risk of OFC (p=1.77×10-4; Bonferroni p=2.00×10-2; OR 0.53 (0.38 to 0.74), minor allele). The rs2124459 was associated with a reduced risk of cleft palate (CP) (p=6.78×10-5; Bonferroni p=7.80×10-3; OR 0.40 (0.25 to 0.63)). In the 'children' group, the rs2124459 was associated with a reduced risk of CP (p=0.02; OR 0.61 (0.40 to 0.93), minor allele). The association between rs2124459 and reduced risk of CP was replicated in an independent children population from Belgium (p=0.02; OR 0.64 (0.44 to 0.93), minor allele). CONCLUSIONS: The CBS rs2124459 was associated with a reduced risk of CP in both French and Belgian populations. These results highlight the prominent involvement of the vitamin B6-dependent transsulfuration pathway of homocysteine in OFC risk and the interest for evaluating vitamin B6 status in further population studies.

Exome-Wide Association Study Identifies New Low-Frequency and Rare UGT1A1 Coding Variants and UGT1A6 Coding Variants Influencing Serum Bilirubin in Elderly Subjects: A Strobe Compliant Article.

Genome-wide association studies (GWASs) have identified loci contributing to total serum bilirubin level. However, no exome-wide approaches have been performed to address this question. Using exome-wide approach, we assessed the influence of protein-coding variants on unconjugated, conjugated, and total serum bilirubin levels in a well-characterized cohort of 773 ambulatory elderly subjects from Italy. Coding variants were replicated in 227 elderly subjects from the same area. We identified 4 missense rare (minor allele frequency, MAF < 0.5%) and low-frequency (MAF, 0.5%-5%) coding variants located in the first exon of the UGT1A1 gene, which encodes for the substrate-binding domain (rs4148323 [MAF = 0.06%; p.Gly71Arg], rs144398951 [MAF = 0.06%; p.Ile215Val], rs35003977 [MAF = 0.78%; p.Val225Gly], and rs57307513 [MAF = 0.06%; p.Ser250Pro]). These variants were in strong linkage disequilibrium with 3 intronic UGT1A1 variants (rs887829, rs4148325, rs6742078), which were significantly associated with total bilirubin level (P = 2.34 × 10(-34), P = 7.02 × 10(-34), and P = 8.27 × 10(-34)), as well as unconjugated, and conjugated bilirubin levels. We also identified UGT1A6 variants in association with total (rs6759892, p.Ser7Ala, P = 1.98 × 10(-26); rs2070959, p.Thr181Ala, P = 2.87 × 10(-27); and rs1105879, p.Arg184Ser, P = 3.27 × 10(-29)), unconjugated, and conjugated bilirubin levels. All UGT1A1 intronic variants (rs887829, rs6742078, and rs4148325) and UGT1A6 coding variants (rs6759892, rs2070959, and rs1105879) were significantly associated with gallstone-related cholecystectomy risk. The UGT1A6 variant rs2070959 (p.Thr181Ala) was associated with the highest risk of gallstone-related cholecystectomy (OR, 4.58; 95% CI, 1.58-13.28; P = 3.21 × 10(-3)). Using an exome-wide approach we identified coding variants on UGT1A1 and UGT1A6 genes in association with serum bilirubin level and hyperbilirubinemia risk in elderly subjects. UGT1A1 intronic single-nucleotide polymorphisms (SNPs) (rs6742078, rs887829, rs4148324) serve as proxy markers for the low-frequency and rare UGT1A1 variants, thereby providing mechanistic explanation to the relationship between UGT1A1 intronic SNPs and the UGT1A1 enzyme activity. UGT1A1 and UGT1A6 variants might be potentially associated with gallstone-related cholecystectomy risk.

HLA-DRA variants predict penicillin allergy in genome-wide fine-mapping genotyping.

BACKGROUND: Immediate reactions to ß-lactams are the most common causes of anaphylactic reactions and can be life-threatening. The few known genetic factors influencing these reactions suggest a link with atopy and inflammation. OBJECTIVE: We performed a fine-mapping genome-wide association study of the genetic predictors of ß-lactam allergy to better understand the underlying mechanisms. METHODS: We studied 387 patients with immediate allergic reactions to ß-lactams and 1124 paired control subjects from Spain. We replicated the results in 299 patients and 362 paired control subjects from Italy. RESULTS: We found significant associations with the single nucleotide polymorphisms rs4958427 of ZNF300 (c.64-471G>A, P = 9.9 × 10(-9)), rs17612 of C5 (c.4311A>C [p.Glu1437Asp], P = 7.5 × 10(-7)), rs7754768 and rs9268832 of the HLA-DRA | HLA-DRB5 interregion (P = 1.6 × 10(-6) and 4.9 × 10(-6)), and rs7192 of HLA-DRA (c.724T>G [p.Leu242Val], P = 7.4 × 10(-6)) in an allelic model, with similar results in an additive model. Single nucleotide polymorphisms of HLA-DRA and ZNF300 predicted skin test positivity to amoxicillin and other penicillins but not to cephalosporins. A haplotype block in HLA-DRA and the HLA-DRA | HLA-DRB5 interregion encompassed a motif involved in balanced expression of the α- and ß-chains of MHC class II, whereas rs7192 was predicted to influence α-chain conformation. HLA-DRA rs7192 and rs8084 were significantly associated with allergy to penicillins and amoxicillin (P = 6.0 × 10(-4) and P = 4.0 × 10(-4), respectively) but not to cephalosporins in the replication study. CONCLUSIONS: Gene variants of HLA-DRA and the HLA-DRA | HLA-DRB5 interregion were significant predictors of allergy to penicillins but not to cephalosporins. These data suggest complex gene-environment interactions in which genetic susceptibility of HLA type 2 antigen presentation plays a central role.