Neutrophil extracellular trap formation and gene programs distinguish TST/IGRA sensitization outcomes among Mycobacterium tuberculosis exposed persons living with HIV.

Persons living with HIV (PLWH) have an increased risk for tuberculosis (TB). After prolonged and repeated exposure, some PLWH never develop TB and show no evidence of immune sensitization to Mycobacterium tuberculosis (Mtb) as defined by persistently negative tuberculin skin tests (TST) and interferon gamma release assays (IGRA). This group has been identified and defined as HIV+ persistently TB, tuberculin and IGRA negative (HITTIN). To investigate potential innate mechanisms unique to individuals with the HITTIN phenotype we compared their neutrophil Mtb infection response to that of PLWH, with no TB history, but who test persistently IGRA positive, and tuberculin positive (HIT). Neutrophil samples from 17 HITTIN (PMNHITTIN) and 11 HIT (PMNHIT) were isolated and infected with Mtb H37Rv for 1h and 6h. RNA was extracted and used for RNAseq analysis. Since there was no significant differential transcriptional response at 1h between infected PMNHITTIN and PMNHIT, we focused on the 6h timepoint. When compared to uninfected PMN, PMNHITTIN displayed 3106 significantly upregulated and 3548 significantly downregulated differentially expressed genes (DEGs) (absolute cutoff of a log2FC of 0.2, FDR < 0.05) whereas PMNHIT demonstrated 3816 significantly upregulated and 3794 significantly downregulated DEGs following 6h Mtb infection. Contrasting the log2FC 6h infection response to Mtb from PMNHITTIN against PMNHIT, 2285 genes showed significant differential response between the two groups. Overall PMNHITTIN had a lower fold change response to Mtb infection compared to PMNHIT. According to pathway enrichment, Apoptosis and NETosis were differentially regulated between HITTIN and HIT PMN responses after 6h Mtb infection. To corroborate the blunted NETosis transcriptional response measured among HITTIN, fluorescence microscopy revealed relatively lower neutrophil extracellular trap formation and cell loss in PMNHITTIN compared to PMNHIT, showing that PMNHITTIN have a distinct response to Mtb.

A systems biology approach identifies candidate drugs to reduce mortality in severely ill patients with COVID-19.

Despite the availability of highly efficacious vaccines, coronavirus disease 2019 (COVID-19) caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) lacks effective drug treatment, which results in a high rate of mortality. To address this therapeutic shortcoming, we applied a systems biology approach to the study of patients hospitalized with severe COVID. We show that, at the time of hospital admission, patients who were equivalent on the clinical ordinal scale displayed significant differential monocyte epigenetic and transcriptomic attributes between those who would survive and those who would succumb to COVID-19. We identified messenger RNA metabolism, RNA splicing, and interferon signaling pathways as key host responses overactivated by patients who would not survive. Those pathways are prime drug targets to reduce mortality of critically ill patients with COVID-19, leading us to identify tacrolimus, zotatifin, and nintedanib as three strong candidates for treatment of severely ill patients at the time of hospital admission.

Identification of new ETV6 modulators through a high-throughput functional screening.

ETV6 transcriptional activity is critical for proper blood cell development in the bone marrow. Despite the accumulating body of evidence linking ETV6 malfunction to hematological malignancies, its regulatory network remains unclear. To uncover genes that modulate ETV6 repressive transcriptional activity, we performed a specifically designed, unbiased genome-wide shRNA screen in pre-B acute lymphoblastic leukemia cells. Following an extensive validation process, we identified 13 shRNAs inducing overexpression of ETV6 transcriptional target genes. We showed that the silencing of AKIRIN1, COMMD9, DYRK4, JUNB, and SRP72 led to an abrogation of ETV6 repressive activity. We identified critical modulators of the ETV6 function which could participate in cellular transformation through the ETV6 transcriptional network.

Alveolar macrophages from persons living with HIV show impaired epigenetic response to Mycobacterium tuberculosis.

Persons living with HIV (PLWH) are at increased risk of tuberculosis (TB). HIV-associated TB is often the result of recent infection with Mycobacterium tuberculosis (M. tuberculosis) followed by rapid progression to disease. Alveolar macrophages (AMs) are the first cells of the innate immune system that engage M. tuberculosis, but how HIV and antiretroviral therapy (ART) affect the anti-mycobacterial response of AMs is not known. To investigate the impact of HIV and ART on the transcriptomic and epigenetic response of AMs to M. tuberculosis, we obtained AMs by bronchoalveolar lavage from 20 PLWH receiving ART, 16 control subjects who were HIV-free (HC), and 14 subjects who received ART as preexposure prophylaxis (PrEP) to prevent HIV infection. Following in vitro challenge with M. tuberculosis, AMs from each group displayed overlapping but distinct profiles of significantly up- and downregulated genes in response to M. tuberculosis. Comparatively, AMs isolated from both PLWH and PrEP subjects presented a substantially weaker transcriptional response. In addition, AMs from HC subjects challenged with M. tuberculosis responded with pronounced chromatin accessibility changes while AMs obtained from PLWH and PrEP subjects displayed no significant changes in their chromatin state. Collectively, these results revealed a stronger adverse effect of ART than HIV on the epigenetic landscape and transcriptional responsiveness of AMs.

Functional Analysis of Promoter Variants in Genes Involved in Sex Steroid Action, DNA Repair and Cell Cycle Control.

Genetic variants affecting the regulation of gene expression are among the main causes of human diversity. The potential importance of regulatory polymorphisms is underscored by results from Genome Wide Association Studies, which have already implicated such polymorphisms in the susceptibility to complex diseases such as breast cancer. In this study, we re-sequenced the promoter regions of 24 genes involved in pathways related to breast cancer including sex steroid action, DNA repair, and cell cycle control in 60 unrelated Caucasian individuals. We constructed haplotypes and assessed the functional impact of promoter variants using gene reporter assays and electrophoretic mobility shift assays. We identified putative functional variants within the promoter regions of estrogen receptor 1 (ESR1), ESR2, forkhead box A1 (FOXA1), ubiquitin interaction motif containing 1 (UIMC1) and cell division cycle 7 (CDC7). The functional polymorphism on CDC7, rs13447455, influences CDC7 transcriptional activity in an allele-specific manner and alters DNA⁻protein complex formation in breast cancer cell lines. Moreover, results from the Breast Cancer Association Consortium show a marginal association between rs13447455 and breast cancer risk (p=9.3x10-5), thus warranting further investigation. Furthermore, our study has helped provide methodological solutions to some technical difficulties that were encountered with gene reporter assays, particularly regarding inter-clone variability and statistical consistency.

Recurrent somatic BRAF insertion (p.V504_R506dup): a tumor marker and a potential therapeutic target in pilocytic astrocytoma.

Pilocytic astrocytoma (PA) is emerging as a tumor entity with dysregulated RAS/RAF/MEK/ERK signaling. In this study, we report the identification of a novel recurrent BRAF insertion (p.V504_R506dup) in five PA cases harboring exclusively this somatic tandem duplication. This recurrent alteration leads to an addition of three amino acids in the kinase domain of BRAF and has functional impact on activating MAPK phosphorylation. Importantly, we show that this mutation confers resistance to RAF inhibitors without changing effectiveness while downstream MEK inhibitors remain effective. Our results further emphasize the importance of BRAF alterations in PA and the need to characterize them in a given tumor as this can affect therapeutic strategies and their potential use as tumor marker in molecular diagnostics.

Mutational dynamics of early and late relapsed childhood ALL: rapid clonal expansion and long-term dormancy.

Childhood acute lymphoblastic leukemia (cALL) is the most frequent pediatric cancer. Refractory/relapsed cALL presents a survival rate of ∼45% and is still one of the leading causes of death by disease among children. Mechanisms, such as clonal competition and evolutionary adaptation, govern treatment resistance. However, the underlying clonal dynamics leading to multiple relapses and differentiating early (<36 months postdiagnosis) from late relapse events remain elusive. Here, we use an integrative genome-based analysis combined with serial sampling of relapsed tumors (from primary tumor to ≤4 relapse events) from 19 pre-B-cell cALL patients (8 early and 11 late relapses) to assess the fitness of somatic mutations and infer their ancestral relationships. By quantifying both general clonal dynamics and newly acquired subclonal diversity, we show that 2 distinct evolutionary patterns govern early and late relapse: on one hand, a highly dynamic pattern, sustained by a putative defect of DNA repair processes, illustrating the quick emergence of fitter clones, and on the other hand, a quasi-inert evolution pattern, suggesting the escape from dormancy of leukemia stem cells likely spared from initial cytoreductive therapy. These results offer new insights into cALL relapse mechanisms and highlight the pressing need for adapted treatment strategies to circumvent resistance mechanisms.

SNooPer: a machine learning-based method for somatic variant identification from low-pass next-generation sequencing.

BACKGROUND: Next-generation sequencing (NGS) allows unbiased, in-depth interrogation of cancer genomes. Many somatic variant callers have been developed yet accurate ascertainment of somatic variants remains a considerable challenge as evidenced by the varying mutation call rates and low concordance among callers. Statistical model-based algorithms that are currently available perform well under ideal scenarios, such as high sequencing depth, homogeneous tumor samples, high somatic variant allele frequency (VAF), but show limited performance with sub-optimal data such as low-pass whole-exome/genome sequencing data. While the goal of any cancer sequencing project is to identify a relevant, and limited, set of somatic variants for further sequence/functional validation, the inherently complex nature of cancer genomes combined with technical issues directly related to sequencing and alignment can affect either the specificity and/or sensitivity of most callers. RESULTS: For these reasons, we developed SNooPer, a versatile machine learning approach that uses Random Forest classification models to accurately call somatic variants in low-depth sequencing data. SNooPer uses a subset of variant positions from the sequencing output for which the class, true variation or sequencing error, is known to train the data-specific model. Here, using a real dataset of 40 childhood acute lymphoblastic leukemia patients, we show how the SNooPer algorithm is not affected by low coverage or low VAFs, and can be used to reduce overall sequencing costs while maintaining high specificity and sensitivity to somatic variant calling. When compared to three benchmarked somatic callers, SNooPer demonstrated the best overall performance. CONCLUSIONS: While the goal of any cancer sequencing project is to identify a relevant, and limited, set of somatic variants for further sequence/functional validation, the inherently complex nature of cancer genomes combined with technical issues directly related to sequencing and alignment can affect either the specificity and/or sensitivity of most callers. The flexibility of SNooPer's random forest protects against technical bias and systematic errors, and is appealing in that it does not rely on user-defined parameters. The code and user guide can be downloaded at https://sourceforge.net/projects/snooper/ .

Genomic characterization of pediatric T-cell acute lymphoblastic leukemia reveals novel recurrent driver mutations.

T-cell acute lymphoblastic leukemia (T-ALL) is an aggressive hematologic malignancy with variable prognosis. It represents 15% of diagnosed pediatric ALL cases and has a threefold higher incidence among males. Many recurrent alterations have been identified and help define molecular subgroups of T-ALL, however the full range of events involved in driving transformation remain to be defined. Using an integrative approach combining genomic and transcriptomic data, we molecularly characterized 30 pediatric T-ALLs and identified common recurrent T-ALL targets such as FBXW7, JAK1, JAK3, PHF6, KDM6A and NOTCH1 as well as novel candidate T-ALL driver mutations including the p.R35L missense mutation in splicesome factor U2AF1 found in 3 patients and loss of function mutations in the X-linked tumor suppressor genes MED12 (frameshit mutation p.V167fs, splice site mutation g.chrX:70339329T>C, missense mutation p.R1989H) and USP9X (nonsense mutation p.Q117*). In vitro functional studies further supported the putative role of these novel T-ALL genes in driving transformation. U2AF1 p.R35L was shown to induce aberrant splicing of downstream target genes, and shRNA knockdown of MED12 and USP9X was shown to confer resistance to apoptosis following T-ALL relevant chemotherapy drug treatment in Jurkat leukemia cells. Interestingly, nearly 60% of novel candidate driver events were identified among immature T-ALL cases, highlighting the underlying genomic complexity of pediatric T-ALL, and the need for larger integrative studies to decipher the mechanisms that contribute to its various subtypes and provide opportunities to refine patient stratification and treatment.

CLIC5: a novel ETV6 target gene in childhood acute lymphoblastic leukemia.

The most common rearrangement in childhood precursor B-cell acute lymphoblastic leukemia is the t(12;21)(p13;q22) translocation resulting in the ETV6-AML1 fusion gene. A frequent concomitant event is the loss of the residual ETV6 allele suggesting a critical role for the ETV6 transcriptional repressor in the etiology of this cancer. However, the precise mechanism through which loss of functional ETV6 contributes to disease pathogenesis is still unclear. To investigate the impact of ETV6 loss on the transcriptional network and to identify new transcriptional targets of ETV6, we used whole transcriptome analysis of both pre-B leukemic cell lines and patients combined with chromatin immunoprecipitation. Using this integrative approach, we identified 4 novel direct ETV6 target genes: CLIC5, BIRC7, ANGPTL2 and WBP1L To further evaluate the role of chloride intracellular channel protein CLIC5 in leukemogenesis, we generated cell lines overexpressing CLIC5 and demonstrated an increased resistance to hydrogen peroxide-induced apoptosis. We further described the implications of CLIC5's ion channel activity in lysosomal-mediated cell death, possibly by modulating the function of the transferrin receptor with which it colocalizes intracellularly. For the first time, we showed that loss of ETV6 leads to significant overexpression of CLIC5, which in turn leads to decreased lysosome-mediated apoptosis. Our data suggest that heightened CLIC5 activity could promote a permissive environment for oxidative stress-induced DNA damage accumulation, and thereby contribute to leukemogenesis.

A novel somatic mutation in ACD induces telomere lengthening and apoptosis resistance in leukemia cells.

BACKGROUND: The identification of oncogenic driver mutations has largely relied on the assumption that genes that exhibit more mutations than expected by chance are more likely to play an active role in tumorigenesis. Major cancer sequencing initiatives have therefore focused on recurrent mutations that are more likely to be drivers. However, in specific genetic contexts, low frequency mutations may also be capable of participating in oncogenic processes. Reliable strategies for identifying these rare or even patient-specific (private) mutations are needed in order to elucidate more personalized approaches to cancer diagnosis and treatment. METHODS: Here we performed whole-exome sequencing on three cases of childhood pre-B acute lymphoblastic leukemia (cALL), representing three cytogenetically-defined subgroups (high hyperdiploid, t(12;21) translocation, and cytogenetically normal). We applied a data reduction strategy to identify both common and rare/private somatic events with high functional potential. Top-ranked candidate mutations were subsequently validated at high sequencing depth on an independent platform and in vitro expression assays were performed to evaluate the impact of identified mutations on cell growth and survival. RESULTS: We identified 6 putatively damaging non-synonymous somatic mutations among the three cALL patients. Three of these mutations were well-characterized common cALL mutations involved in constitutive activation of the mitogen-activated protein kinase pathway (FLT3 p.D835Y, NRAS p.G13D, BRAF p.G466A). The remaining three patient-specific mutations (ACD p.G223V, DOT1L p.V114F, HCFC1 p.Y103H) were novel mutations previously undescribed in public cancer databases. Cytotoxicity assays demonstrated a protective effect of the ACD p.G223V mutation against apoptosis in leukemia cells. ACD plays a key role in protecting telomeres and recruiting telomerase. Using a telomere restriction fragment assay, we also showed that this novel mutation in ACD leads to increased telomere length in leukemia cells. CONCLUSION: This study identified ACD as a novel gene involved in cALL and points to a functional role for ACD in enhancing leukemia cell survival. These results highlight the importance of rare/private somatic mutations in understanding cALL etiology, even within well-characterized molecular subgroups.

Whole-exome sequencing of a rare case of familial childhood acute lymphoblastic leukemia reveals putative predisposing mutations in Fanconi anemia genes.

BACKGROUND: Acute lymphoblastic leukemia (ALL) is the most common pediatric cancer. While the multi-step model of pediatric leukemogenesis suggests interplay between constitutional and somatic genomes, the role of inherited genetic variability remains largely undescribed. Nonsyndromic familial ALL, although extremely rare, provides the ideal setting to study inherited contributions to ALL. Toward this goal, we sequenced the exomes of a childhood ALL family consisting of mother, father and two non-twinned siblings diagnosed with concordant pre-B hyperdiploid ALL and previously shown to have inherited a rare form of PRDM9, a histone H3 methyltransferase involved in crossing-over at recombination hotspots and Holliday junctions. We postulated that inheritance of additional rare disadvantaging variants in predisposing cancer genes could affect genomic stability and lead to increased risk of hyperdiploid ALL within this family. METHODS: Whole exomes were captured using Agilent's SureSelect kit and sequenced on the Life Technologies SOLiD System. We applied a data reduction strategy to identify candidate variants shared by both affected siblings. Under a recessive disease model, we focused on rare non-synonymous or frame-shift variants in leukemia predisposing pathways. RESULTS: Though the family was nonsyndromic, we identified a combination of rare variants in Fanconi anemia (FA) genes FANCP/SLX4 (compound heterozygote - rs137976282/rs79842542) and FANCA (rs61753269) and a rare homozygous variant in the Holliday junction resolvase GEN1 (rs16981869). These variants, predicted to affect protein function, were previously identified in familial breast cancer cases. Based on our in-house database of 369 childhood ALL exomes, the sibs were the only patients to carry this particularly rare combination and only a single hyperdiploid patient was heterozygote at both FANCP/SLX4 positions, while no FANCA variant allele carriers were identified. FANCA is the most commonly mutated gene in FA and is essential for resolving DNA interstrand cross-links during replication. FANCP/SLX4 and GEN1 are involved in the cleavage of Holliday junctions and their mutated forms, in combination with the rare allele of PRDM9, could alter Holliday junction resolution leading to nondisjunction of chromosomes and segregation defects. CONCLUSION: Taken together, these results suggest that concomitant inheritance of rare variants in FANCA, FANCP/SLX4 and GEN1 on the specific genetic background of this familial case, could lead to increased genomic instability, hematopoietic dysfunction, and higher risk of childhood leukemia.

Bacteriophage mv4 site-specific recombination: the central role of the P2 mv4Int-binding site.

The contributions of the five (mv4)Int- and two (mv4)Xis arm-binding sites to the spatial intasome organization of bacteriophage mv4 were found not to be equivalent. The 8-bp overlap region was mapped to the left extremity of the core region and is directly flanked by the P2 Int arm-binding site. These results and the absence of characteristic Int core-binding sites suggest that the P2 site is the determinant for integrase positioning and recognition of the core region.

Functional analysis of promoter variants in KU70 and their role in cancer susceptibility.

KU70 is involved in the DNA double-strand break repair pathway, which plays a critical role in maintaining genomic stability and preventing cancer. Genetic variation within the KU70 gene has been shown to be associated with increased risk of several types of cancers including breast cancer. Here, we used gene reporter and gel shift assays combined with site-directed mutagenesis to characterize genetic variation within the KU70 proximal promoter region and investigate the putative functional role of regulatory variation and altered KU70 expression in modulating an individual's susceptibility to disease. We show that the variant rs2267437C>G significantly influences KU70 transcriptional activity in an allele- specific manner and alters DNA-protein complex formation in breast cancer cell lines. Our data provide a possible molecular explanation for the associations observed between the KU70 regulatory variant rs2267437 and breast cancer risk.