Abstract 3155: Gene expression profiles reveal distinct regulatory activities of transcription factors GATA1 and TAL1 upon AML relapse

The purpose of this study is to identify key regulatory pathways that potentially drive abnormal gene expression program in relapsed Acute Myeloid Leukemia (AML) patients, by integrative computational analyses on multi-omics molecular profiles. Relapsed AML remains a clinical challenge. Epigenetic heterogeneity may contribute to transcriptional dysregulation and disease progression in AML. However, what specific transcriptional programs and potential regulatory mechanisms contribute to disease relapse are not yet well understood. To characterize the global transcriptional landscapes in relapsed AML, we integrated genomics data from two cohorts of matched diagnosis and relapse patient specimens. We identified 5,416 differentially expressed genes (DEGs) between diagnosis and relapse in Cohort I. Unsupervised clustering yielded three distinct DEG groups: group A, B and C genes that were predominantly (88%) down-regulated, divergently regulated, or predominantly (65%) up-regulated, respectively, upon relapse. The expression pattern of all DEGs separated the patients into two clusters, most robustly by Group B genes. Interestingly, the majority of DEGs did not associate with changes in gene promoter methylation. Similar patterns were observed in Cohort II. We used Binding Analysis for Regulation of Transcription (BART) to identify transcriptional regulators (TRs) that potentially regulated the DEGs not associated with DNA methylation changes, and assessed the differential expression of identified TRs during disease progression. PU.1 was identified as a potential TR for Group A genes and was down-regulated upon relapse. GATA1 and TAL1 were identified as regulating Group B genes and were up-regulated in patient cluster1 and down-regulated in cluster2, consistent with the expression pattern of Group B genes. RBBP5 was a top predicted TR for Group C genes and was up-regulated upon relapse. We next validated the potential functionality of those predicted factors. In NSG mice transplanted with a human AML specimen, TAL1 and GATA1 were downregulated in AML cells collected four weeks after chemotherapy treatment, and were inferred as TRs for the down-regulated genes, similar to the patient data. PU.1 was inferred as regulating the up-regulated genes. Furthermore, we found that the level of differential expression of TAL1, GATA1, and PU.1 in each patient specimen associated with the correlation of DEG profiles between the patient specimen and TR perturbation in human-derived hematopoietic cell lines. Our results support the possibility that in some AML patients, TRs with roles in hematopoiesis and leukemia might contribute to disease relapse. Further mechanistic studies deciphering the molecular and phenotypic events facilitated by these TRs will yield significant insight into disease biology and possible therapeutic targeting approaches in relapsed AML.

Citation Format: Zhenjia Wang, Yaseswini Neelamraju, Cem Meydan, Nicholas Dunham, Jorge Gandara, Tak Lee, Subhash Prajapati, Franck Rapaport, Caroline Sheridan, Paul Zumbo, Michael Becker, Lars Bullinger, Martin Carroll, Richard D’Andrea, Richard Dillon, Ross Levine, Christopher E. Mason, Ari Melnick, Donna Neuberg, Stefan Bekiranov, Chongzhi Zang, Francine E. Garrett-Bakelman. Gene expression profiles reveal distinct regulatory activities of transcription factors GATA1 and TAL1 upon AML relapse [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 3155.

Validation of the CARPREG II risk stratification model and the WHOm scale in pregnant women with heart disease

 

The physiological changes of pregnancy imply a state of haemodynamic stress, which increases the risk of maternal-fetal complications in women with cardiac pathology. Risk stratification models allow clinical decisions to be guided and optimal diagnostic, follow-up and management strategies to be established according to each estimated category in order to reduce adverse outcomes.

Objective

To validate the CARPREG II and WHOm risk prediction models in a population of pregnant women with heart disease.

Materials and methods

Validation and comparison study of a risk prediction model in a prospective cohort of pregnant women with heart disease assessed between 2016 and 2019 by a cardiobstetric team in an average income country. The CARPREG II score was established and patients were classified according to the WHOm risk scale. We assessed cardiovascular and perinatal outcomes and determined the calibration and level of discrimination of these tools.

Results

In a cohort of 328 pregnant women (27±7 years), 33% (n=110) had congenital heart disease, followed by arrhythmias in 30% (n=98), valvular pathologies in 14% (n=46) and cardiomyopathies in 9% (n=29). In 56% of the pregnancies, the route of delivery was caesarean section, 84% of these by obstetric indication. A cardiac event occurred in 15%, with left heart failure (5.3%) and arrhythmias (2.3%) being the most frequent. The frequency of maternal death of cardiac origin during the study was 1.6%. Neonatal outcomes occurred in 37% of gestations (preterm delivery (16%) and low weight for gestational age (8.4%)) and obstetric events in 12.5%: pregnancy-induced hypertension (9%) and postpartum haemorrhage (2.3%).

NYHA functional class III-IV or cyanosis (OR 12 95% CI 3.1 - 46.4) and left ventricular dysfunction (LVEF <55%) (OR 3 95% CI 1 - 10.9) were the most statistically significant risk predictors. Discrimination of both models was adequate (AUC-ROC of 0.74 95% CI 0.64 - 0.84) for the CARPREG II risk index and 0.77 for the WHOm scale (95% CI 0.69 - 0.86) (Figure 2). Calibration is also good in the study population (Hosmer- Lemeshow goodness-of-fit 0.6 and 0.1, respectively). By including in the CARPREG II model the variables ejection fraction and pulmonary artery systolic pressure in their numerical and not dichotomised form, a discrete improvement in the predictive ability of the scale is evident (AUC-ROC 0.81 95% CI 0.71–0.91).

Conclusions

The CARPREG II and WHOm risk stratification models have good ability to discriminate the risk of adverse cardiac outcomes in pregnant women with heart disease and fit our population. To improve the predictive power of CARPREG II, the variables pulmonary hypertension and left ventricular dysfunction could be used numerically and not dichotomised as in the original model.

Funding Acknowledgement

Type of funding sources: Private hospital(s). Main funding source(s): San Vicente Foundation University Hospital, Cardio-obstetric Center

Microbiota-Induced TNF-like Ligand 1A Drives Group 3 Innate Lymphoid Cell-Mediated Barrier Protection and Intestinal T Cell Activation during Colitis

Inflammatory bowel disease (IBD) results from a dysregulated interaction between the microbiota and a genetically susceptible host. Genetic studies have linked TNFSF15 polymorphisms and its protein TNF-like ligand 1A (TL1A) with IBD, but the functional role of TL1A is not known. Here, we found that adherent IBD-associated microbiota induced TL1A release from CX3CR1+ mononuclear phagocytes (MNPs). Using cell-specific genetic deletion models, we identified an essential role for CX3CR1+MNP-derived TL1A in driving group 3 innate lymphoid cell (ILC3) production of interleukin-22 and mucosal healing during acute colitis. In contrast to this protective role in acute colitis, TL1A-dependent expression of co-stimulatory molecule OX40L in MHCII+ ILC3s during colitis led to co-stimulation of antigen-specific T cells that was required for chronic T cell colitis. These results identify a role for ILC3s in activating intestinal T cells and reveal a central role for TL1A in promoting ILC3 barrier immunity during colitis.

Liver transplant recipients have a higher prevalence of anal squamous intraepithelial lesions

Background:

Anal squamous intraepithelial lesions (ASIL) are precancerous lesions of anal squamous cell carcinoma, with a higher prevalence in immunosuppressed patients. There are some studies in kidney transplant recipients, but there is no information regarding prevalence in liver transplantation. Our aim was to evaluate the prevalence of ASIL in this setting.

Methods:

Prospective case–control study involving liver transplant recipients without any other known risk factor for ASIL (n=59), which were compared with a healthy control group (n=57). All were submitted to anal cytology and high-resolution anoscopy was performed in those with abnormal results.

Results:

Ten (17%) of liver transplant recipients had abnormal cytological results, seven patients had atypical squamous cells of undetermined significance (ASC-US), one patient had atypical squamous cells that cannot exclude high-grade (ASC-H) and two patients had high-grade squamous intraepithelial lesions (HSIL). In the control group, one patient (2%) had an ASC-US result (P=0.005). Anal squamous intraepithelial lesions were confirmed in 7 out of 10 of liver transplant patients and 0 out of 1 in the controls (P=0.013) by high-resolution anoscopy with biopsies. Current smoking was the only risk factor for abnormal cytology (odds ratio=5.87, 95% confidence intervals=1.22–28.12, P=0.027).

Conclusions:

Liver transplant patients have a higher risk of ASIL. Screening should be considered, especially in smokers.

Assessment of REPLI-g Multiple Displacement Whole Genome Amplification (WGA) Techniques for Metagenomic Applications

Amplification of minute quantities of DNA is a fundamental challenge in low-biomass metagenomic and microbiome studies because of potential biases in coverage, guanine-cytosine (GC) content, and altered species abundances. Whole genome amplification (WGA), although widely used, is notorious for introducing artifact sequences, either by amplifying laboratory contaminants or by nonrandom amplification of a sample's DNA. In this study, we investigate the effect of REPLI-g multiple displacement amplification (MDA; Qiagen, Valencia, CA, USA) on sequencing data quality and species abundance detection in 8 paired metagenomic samples and 1 titrated, mixed control sample. We extracted and sequenced genomic DNA (gDNA) from 8 environmental samples and compared the quality of the sequencing data for the MDA and their corresponding non-MDA samples. The degree of REPLI-g MDA bias was evaluated by sequence metrics, species composition, and cross-validating observed species abundance and species diversity estimates using the One Codex and MetaPhlAn taxonomic classification tools. Here, we provide evidence of the overall efficacy of REPLI-g MDA on retaining sequencing data quality and species abundance measurements while providing increased yields of high-fidelity DNA. We find that species abundance estimates are largely consistent across samples, even with REPLI-g amplification, as demonstrated by the Spearman's rank order coefficient (R² > 0.8). However, REPLI-g MDA often produced fewer classified reads at the species, genera, and family level, resulting in decreased species diversity. We also observed some areas with the PCR "jackpot effect," with varying input DNA values for the Metagenomics Research Group (MGRG) controls at specific genomic loci. We visualize this effect in whole genome coverage plots and with sequence composition analyses and note these caveats of the MDA method. Despite overall concordance of species abundance between the amplified and unamplified samples, these results demonstrate that amplification of DNA using the REPLI-g method has some limitations. These concerns could be addressed by future improvements in the enzymes or methods for REPLI-g to be considered a >99% robust method for increasing the amount of high-fidelity DNA from low-biomass samples or at the very least, accounted for during computational analysis of MDA samples.

Genome assembly and geospatial phylogenomics of the bed bug Cimex lectularius

The common bed bug (Cimex lectularius) has been a persistent pest of humans for thousands of years, yet the genetic basis of the bed bug's basic biology and adaptation to dense human environments is largely unknown. Here we report the assembly, annotation and phylogenetic mapping of the 697.9-Mb Cimex lectularius genome, with an N50 of 971 kb, using both long and short read technologies. A RNA-seq time course across all five developmental stages and male and female adults generated 36,985 coding and noncoding gene models. The most pronounced change in gene expression during the life cycle occurs after feeding on human blood and included genes from the Wolbachia endosymbiont, which shows a simultaneous and coordinated host/commensal response to haematophagous activity. These data provide a rich genetic resource for mapping activity and density of C. lectularius across human hosts and cities, which can help track, manage and control bed bug infestations.

Genomic DNA transposition induced by human PGBD5

Transposons are mobile genetic elements that are found in nearly all organisms, including humans. Mobilization of DNA transposons by transposase enzymes can cause genomic rearrangements, but our knowledge of human genes derived from transposases is limited. In this study, we find that the protein encoded by human PGBD5, the most evolutionarily conserved transposable element-derived gene in vertebrates, can induce stereotypical cut-and-paste DNA transposition in human cells. Genomic integration activity of PGBD5 requires distinct aspartic acid residues in its transposase domain, and specific DNA sequences containing inverted terminal repeats with similarity to piggyBac transposons. DNA transposition catalyzed by PGBD5 in human cells occurs genome-wide, with precise transposon excision and preference for insertion at TTAA sites. The apparent conservation of DNA transposition activity by PGBD5 suggests that genomic remodeling contributes to its biological function.

Mining the Archives: A Cross-Platform Analysis of Gene Expression Profiles in Archival Formalin-Fixed Paraffin-Embedded Tissues

Formalin-fixed paraffin-embedded (FFPE) tissue samples represent a potentially invaluable resource for transcriptomic research. However, use of FFPE samples in genomic studies has been limited by technical challenges resulting from nucleic acid degradation. Here we evaluated gene expression profiles derived from fresh-frozen (FRO) and FFPE mouse liver tissues preserved in formalin for different amounts of time using 2 DNA microarray protocols and 2 whole-transcriptome sequencing (RNA-seq) library preparation methodologies. The ribo-depletion protocol outperformed the other methods by having the highest correlations of differentially expressed genes (DEGs), and best overlap of pathways, between FRO and FFPE groups. The effect of sample time in formalin (18 h or 3 weeks) on gene expression profiles indicated that test article treatment, not preservation method, was the main driver of gene expression profiles. Meta- and pathway analyses indicated that biological responses were generally consistent for 18 h and 3 week FFPE samples compared with FRO samples. However, clear erosion of signal intensity with time in formalin was evident, and DEG numbers differed by platform and preservation method. Lastly, we investigated the effect of time in paraffin on genomic profiles. Ribo-depletion RNA-seq analysis of 8-, 19-, and 26-year-old control blocks resulted in comparable quality metrics, including expected distributions of mapped reads to exonic, untranslated region, intronic, and ribosomal fractions of the transcriptome. Overall, our results indicate that FFPE samples are appropriate for use in genomic studies in which frozen samples are not available, and that ribo-depletion RNA-seq is the preferred method for this type of analysis in archival and long-aged FFPE samples.

Abstract 1103: Human tumorigenesis induced by endogenous DNA transposase

Recent cancer genome surveys have revealed extremely low rates of coding gene mutations in distinct tumor subtypes, suggesting that alternative mechanisms must contribute to their pathogenesis. Transposons are mobile genetic elements that are found in all living organisms, including humans where they occupy nearly half of the genome. Their mobilization can cause structural rearrangements in normal and cancer cells. However, it remains unknown whether transposition is a cause of cellular transformation or merely a bystander effect of dysregulated gene expression. Here, we report that PGBD5, a recently characterized human gene related to the piggyBac transposase from the cabbage looper moth, is aberrantly expressed in rhabdoid tumors, medulloblastoma, acute leukemias, and some sarcomas and carcinomas. Ectopic expression of PGBD5 in non-transformed primary human cells is sufficient to induce anchorage independence in vitro and penetrant tumor formation in immunodeficient mice in vivo. PGBD5 expression is sufficient to induce genomic mobilization of engineered DNA transposons in human cells, and purified recombinant PGBD5 exhibits transposase domain-dependent endonuclease activity in vitro. Flanking-sequence exponential anchored PCR and massively parallel sequencing of DNA transposon integrations revealed distinct activity on piggyBac-like inverted terminal repeats, and preference for specific euchromatic human genomic loci. This enables mapping of structural rearrangements of endogenous human transposable elements in primary human tumor genomes, some of which target genes involved in cellular transformation. We find that PGBD5 transposase-induced cell transformation is associated with morphologic de-differentiation, induction of distinct Polycomb gene expression programs and structural chromatin remodeling, consistent with its epigenetic control. These findings reveal an unanticipated mechanism of human tumorigenesis, genomic plasticity and structural alterations of non-coding regulatory genomic loci in human cancer.

Citation Format: Anton Henssen, Amy Eisenberg, Eileen Jiang, Elizabeth Henaff, Richard Koche, Melissa Burns, Julianne R. Carson, Gouri Nanjangud, Eric Still, Jorge Gandara, Paolo Cifani, Avantika Dhabaria, Xiaodong Huang, Elisa de Stanchina, Elizabeth Mullen, Hanno Steen, Elizabeth Perlman, Jeffrey Dome, Cristina Antonescu, Cedric Feschotte, Christopher E. Mason, Alex Kentsis. Human tumorigenesis induced by endogenous DNA transposase. [abstract]. In: Proceedings of the 106th Annual Meeting of the American Association for Cancer Research; 2015 Apr 18-22; Philadelphia, PA. Philadelphia (PA): AACR; Cancer Res 2015;75(15 Suppl):Abstract nr 1103. doi:10.1158/1538-7445.AM2015-1103

Geospatial Resolution of Human and Bacterial Diversity with City-Scale Metagenomics

The panoply of microorganisms and other species present in our environment influence human health and disease, especially in cities, but have not been profiled with metagenomics at a city-wide scale. We sequenced DNA from surfaces across the entire New York City (NYC) subway system, the Gowanus Canal, and public parks. Nearly half of the DNA (48%) does not match any known organism; identified organisms spanned 1,688 bacterial, viral, archaeal, and eukaryotic taxa, which were enriched for harmless genera associated with skin (e.g., Acinetobacter). Predicted ancestry of human DNA left on subway surfaces can recapitulate U.S. Census demographic data, and bacterial signatures can reveal a station’s history, such as marine-associated bacteria in a hurricane-flooded station. Some evidence of pathogens was found (Bacillus anthracis), but a lack of reported cases in NYC suggests that the pathogens represent a normal, urban microbiome. This baseline metagenomic map of NYC could help long-term disease surveillance, bioterrorism threat mitigation, and health management in the built environment of cities.

Multi-platform assessment of transcriptome profiling using RNA-seq in the ABRF next-generation sequencing study

High-throughput RNA sequencing (RNA-seq) greatly expands the potential for genomics discoveries, but the wide variety of platforms, protocols and performance capabilitites has created the need for comprehensive reference data. Here we describe the Association of Biomolecular Resource Facilities next-generation sequencing (ABRF-NGS) study on RNA-seq. We carried out replicate experiments across 15 laboratory sites using reference RNA standards to test four protocols (poly-A-selected, ribo-depleted, size-selected and degraded) on five sequencing platforms (Illumina HiSeq, Life Technologies PGM and Proton, Pacific Biosciences RS and Roche 454). The results show high intraplatform (Spearman rank R > 0.86) and inter-platform (R > 0.83) concordance for expression measures across the deep-count platforms, but highly variable efficiency and cost for splice junction and variant detection between all platforms. For intact RNA, gene expression profiles from rRNA-depletion and poly-A enrichment are similar. In addition, rRNA depletion enables effective analysis of degraded RNA samples. This study provides a broad foundation for cross-platform standardization, evaluation and improvement of RNA-seq.

Experience with indeterminate mammograms

Most radiologists recommend follow-up mammography in 3 to 6 months for asymptomatic patients with mammographic lesions that do not appear suspicious enough to warrant immediate biopsy. In this retrospective cohort study, the medical records of 776 patients were audited 24 to 35 months after mammography to estimate the frequency of indeterminate lesions, the probability of malignancy, and the rate of patients' compliance with recommendations. Almost 2% of patients retested after an indeterminate mammogram were found to have breast cancer (95% confidence limits 0.2% and 5.9%). Only 26% of the patients retested during the audit period had the test within the suggested upper limit of 6 months. We conclude that timely follow-up after an indeterminate mammogram is advisable and that for the population we studied, compliance with this recommendation is low.