Leveraging inter-individual transcriptional correlation structure to infer discrete signaling mechanisms across metabolic tissues.

Inter-organ communication is a vital process to maintain physiologic homeostasis, and its dysregulation contributes to many human diseases. Given that circulating bioactive factors are stable in serum, occur naturally, and are easily assayed from blood, they present obvious focal molecules for therapeutic intervention and biomarker development. Recently, studies have shown that secreted proteins mediating inter-tissue signaling could be identified by 'brute force' surveys of all genes within RNA-sequencing measures across tissues within a population. Expanding on this intuition, we reasoned that parallel strategies could be used to understand how individual genes mediate signaling across metabolic tissues through correlative analyses of gene variation between individuals. Thus, comparison of quantitative levels of gene expression relationships between organs in a population could aid in understanding cross-organ signaling. Here, we surveyed gene-gene correlation structure across 18 metabolic tissues in 310 human individuals and 7 tissues in 103 diverse strains of mice fed a normal chow or high-fat/high-sucrose (HFHS) diet. Variation of genes such as FGF21, ADIPOQ, GCG, and IL6 showed enrichments which recapitulate experimental observations. Further, similar analyses were applied to explore both within-tissue signaling mechanisms (liver PCSK9) and genes encoding enzymes producing metabolites (adipose PNPLA2), where inter-individual correlation structure aligned with known roles for these critical metabolic pathways. Examination of sex hormone receptor correlations in mice highlighted the difference of tissue-specific variation in relationships with metabolic traits. We refer to this resource as gene-derived correlations across tissues (GD-CAT) where all tools and data are built into a web portal enabling users to perform these analyses without a single line of code (gdcat.org). This resource enables querying of any gene in any tissue to find correlated patterns of genes, cell types, pathways, and network architectures across metabolic organs.

HyperCVAD versus pegaspargase-containing regimens for Hispanic adults with newly diagnosed B-cell acute lymphoblastic leukemia.

OBJECTIVE: There are significant disparities in outcomes among Hispanic patients with acute lymphoblastic leukemia (ALL). Recent studies have demonstrated favorable outcomes of pegaspargase-containing ALL regimens (PEG-CAR) in young adults however, outcomes in Hispanic ethnicity continue to be underreported. METHODS: We evaluated outcomes of newly diagnosed, adult B-cell ALL Hispanic and non-Hispanic patients consecutively treated with a PEG-CAR or HyperCVAD between January 2011 and November 2022. The primary endpoint was event-free survival (EFS) while secondary endpoints included cumulative incidence of relapse and overall survival (OS). RESULTS: Among 105 included patients, 48 (45.7%) were treated with a PEG-CAR and 57 (54.3%) with HyperCVAD. Median age was 38 years (range, 18-75 years), 61% were Hispanic, and 35.2% had poor-genetic risk. Hispanic patients demonstrated significantly worse 5-year EFS with a PEG-CAR compared to that seen with HyperCVAD (HR, 2.58; 95% CI, 1.32-5.04; p = .006) whereas non-Hispanic patients had better outcomes with PIR (52.4% vs. 42.0%). Hispanic ethnicity (p = .015) and male sex (p = .019) were independent predictors for poor OS. CONCLUSIONS: Hispanic patients with B-cell ALL had worse EFS with a PEG-CAR as compared with HyperCVAD. Future studies will aim to confirm these findings and establish a tailored treatment approach for this high-risk population.

Increased AID Results in Mutations at the CRLF2 Locus Implicated in Latin American ALL Health Disparities.

Activation-induced cytidine deaminase (AID) is a B cell-specific base editor required during class switch recombination and somatic hypermutation for B cell maturation and antibody diversification. However, it has also been implicated as a factor in the etiology of several B cell malignancies. Evaluating the AID-induced mutation load in patients at-risk for certain types of blood cancers is critical in assessing disease severity and treatment options. Here, we have developed a digital PCR (dPCR) assay that allows us to track the mutational landscape resulting from AID modification or DNA double-strand break (DSB) formation and repair at sites known to be prone to DSBs. Implementation of this new assay showed that increased AID levels in immature B cells increases genome instability at loci linked to translocation formation. This included the CRLF2 locus that is often involved in chromosomal translocations associated with a subtype of acute lymphoblastic leukemia (ALL) that disproportionately affects Latin Americans (LAs). To support this LA-specific identification of AID mutation signatures, we characterized DNA from immature B cells isolated from the bone marrow of ALL patients. Our ability to detect and quantify these mutation signatures will potentiate future risk identification, early detection of cancers, and reduction of associated cancer health disparities.

Disruption of the circadian clock drives Apc loss of heterozygosity to accelerate colorectal cancer.

An alarming rise in young onset colorectal cancer (CRC) has been reported; however, the underlying molecular mechanism remains undefined. Suspected risk factors of young onset CRC include environmental aspects, such as lifestyle and dietary factors, which are known to affect the circadian clock. We find that both genetic disruption and environmental disruption of the circadian clock accelerate Apc-driven CRC pathogenesis in vivo. Using an intestinal organoid model, we demonstrate that clock disruption promotes transformation by driving Apc loss of heterozygosity, which hyperactivates Wnt signaling. This up-regulates c-Myc, a known Wnt target, which drives heightened glycolytic metabolism. Using patient-derived organoids, we show that circadian rhythms are lost in human tumors. Last, we identify that variance between core clock and Wnt pathway genes significantly predicts the survival of patients with CRC. Overall, our findings demonstrate a previously unidentified mechanistic link between clock disruption and CRC, which has important implications for young onset cancer prevention.

Diversity upon diversity: linking DNA double-strand break repair to blood cancer health disparities.

Chromosomal translocations arising from aberrant repair of multiple DNA double-strand breaks (DSBs) are a defining characteristic of many cancers. DSBs are an essential part of physiological processes in antibody-producing B cells. The B cell environment is poised to generate genome instability leading to translocations relevant to the pathology of blood cancers. These are a diverse set of cancers, but limited data from under-represented groups have pointed to health disparities associated with each. We focus on the DSBs that occur in developing B cells and propose the most likely mechanism behind the formation of translocations. We also highlight specific cancers in which these rearrangements occur and address the growing concern of health disparities associated with them.

Constitutively active Artemis nuclease recognizes structures containing single-stranded DNA configurations.

The Artemis nuclease recognizes and endonucleolytically cleaves at single-stranded to double-stranded DNA (ss/dsDNA) boundaries. It is also a key enzyme in the non-homologous end joining (NHEJ) DNA double-strand break repair pathway. Previously, a truncated form, Artemis-413, was developed that is constitutively active both in vitro and in vivo. Here, we use this constitutively active form of Artemis to detect DNA structures with ss/dsDNA boundaries that arise under topological stress. Topoisomerases prevent abnormal levels of torsional stress through modulation of positive and negative supercoiling. We show that overexpression of Artemis-413 in yeast cells carrying genetic mutations that ablate topoisomerase activity have an increased frequency of DNA double-strand breaks (DSBs). Based on the biochemical activity of Artemis, this suggests an increase in ss/dsDNA-containing structures upon increased torsional stress, with DSBs arising due to Artemis cutting at these ss/dsDNA structures. Camptothecin targets topoisomerase IB (Top1), and cells treated with camptothecin show increased DSBs. We find that expression of Artemis-413 in camptothecin-treated cells leads to a reduction in DSBs, the opposite of what we find with topoisomerase genetic mutations. This contrast between outcomes not only confirms that topoisomerase mutation and topoisomerase poisoning have distinct effects on cells, but also demonstrates the usefulness of Artemis-413 to study changes in DNA structure.

Concept of DNA Lesion Longevity and Chromosomal Translocations.

A subset of chromosomal translocations related to B cell malignancy in human patients arises due to DNA breaks occurring within defined 20-600 base pair (bp) zones. Several factors influence the breakage rate at these sites including transcription, DNA sequence, and topological tension. These factors favor non-B DNA structures that permit formation of transient single-stranded DNA (ssDNA), making the DNA more vulnerable to agents such as the enzyme activation-induced cytidine deaminase (AID) and reactive oxygen species (ROS). Certain DNA lesions created during the ssDNA state persist after the DNA resumes its normal duplex structure. We propose that factors favoring both formation of transient ssDNA and persistent DNA lesions are key in determining the DNA breakage mechanism.

Nonhomologous DNA end-joining for repair of DNA double-strand breaks.

Nonhomologous DNA end-joining (NHEJ) is the predominant double-strand break (DSB) repair pathway throughout the cell cycle and accounts for nearly all DSB repair outside of the S and G2 phases. NHEJ relies on Ku to thread onto DNA termini and thereby improve the affinity of the NHEJ enzymatic components consisting of polymerases (Pol µ and Pol λ), a nuclease (the Artemis·DNA-PKcs complex), and a ligase (XLF·XRCC4·Lig4 complex). Each of the enzymatic components is distinctive for its versatility in acting on diverse incompatible DNA end configurations coupled with a flexibility in loading order, resulting in many possible junctional outcomes from one DSB. DNA ends can either be directly ligated or, if the ends are incompatible, processed until a ligatable configuration is achieved that is often stabilized by up to 4 bp of terminal microhomology. Processing of DNA ends results in nucleotide loss or addition, explaining why DSBs repaired by NHEJ are rarely restored to their original DNA sequence. Thus, NHEJ is a single pathway with multiple enzymes at its disposal to repair DSBs, resulting in a diversity of repair outcomes.

AID and Reactive Oxygen Species Can Induce DNA Breaks within Human Chromosomal Translocation Fragile Zones.

DNA double-strand breaks (DSBs) occurring within fragile zones of less than 200 base pairs account for the formation of the most common human chromosomal translocations in lymphoid malignancies, yet the mechanism of how breaks occur remains unknown. Here, we have transferred human fragile zones into S. cerevisiae in the context of a genetic assay to understand the mechanism leading to DSBs at these sites. Our findings indicate that a combination of factors is required to sensitize these regions. Foremost, DNA strand separation by transcription or increased torsional stress can expose these DNA regions to damage from either the expression of human AID or increased oxidative stress. This damage causes DNA lesions that, if not repaired quickly, are prone to nuclease cleavage, resulting in DSBs. Our results provide mechanistic insight into why human neoplastic translocation fragile DNA sequences are more prone to enzymes or agents that cause longer-lived DNA lesions.

Non-homologous DNA end joining and alternative pathways to double-strand break repair.

DNA double-strand breaks (DSBs) are the most dangerous type of DNA damage because they can result in the loss of large chromosomal regions. In all mammalian cells, DSBs that occur throughout the cell cycle are repaired predominantly by the non-homologous DNA end joining (NHEJ) pathway. Defects in NHEJ result in sensitivity to ionizing radiation and the ablation of lymphocytes. The NHEJ pathway utilizes proteins that recognize, resect, polymerize and ligate the DNA ends in a flexible manner. This flexibility permits NHEJ to function on a wide range of DNA-end configurations, with the resulting repaired DNA junctions often containing mutations. In this Review, we discuss the most recent findings regarding the relative involvement of the different NHEJ proteins in the repair of various DNA-end configurations. We also discuss the shunting of DNA-end repair to the auxiliary pathways of alternative end joining (a-EJ) or single-strand annealing (SSA) and the relevance of these different pathways to human disease.

RNA Polymerase Collision versus DNA Structural Distortion: Twists and Turns Can Cause Break Failure.

The twisting of DNA due to the movement of RNA polymerases is the basis of numerous classic experiments in molecular biology. Recent mouse genetic models indicate that chromosomal breakage is common at sites of transcriptional turbulence. Two key studies on this point mapped breakpoints to sites of either convergent or divergent transcription but arrived at different conclusions as to which is more detrimental and why. The issue hinges on whether DNA strand separation is the basis for the chromosomal instability or collision of RNA polymerases.

Dissecting the Roles of Divergent and Convergent Transcription in Chromosome Instability.

The interplay of transcription, topological tension, and chromosome breakage is a subject of intense interest, but, with so many facets to the problem, it is difficult to test. Here, we vary the orientation of promoters relative to one another in a yeast system that permits sensitive detection of chromosome breaks. Interestingly, convergent transcription that would direct RNA polymerases into one another does not increase chromosome breakage. In contrast, divergent transcription that would create underwound and potentially single-stranded DNA does cause a marked increase in chromosome breakage. Furthermore, we examine the role that topoisomerases are playing in preventing genome instability at these promoters and find that Top2 is required to prevent instability at converging promoters.

The repetitive portion of the Xenopus IgH Mu switch region mediates orientation-dependent class switch recombination.

Vertebrates developed immunoglobulin heavy chain (IgH) class switch recombination (CSR) to express different IgH constant regions. Most double-strand breaks for Ig CSR occur within the repetitive portion of the switch regions located upstream of each set of constant domain exons for the Igγ, Igα or Igϵ heavy chain. Unlike mammalian switch regions, Xenopus switch regions do not have a high G-density on the non-template DNA strand. In previous studies, when Xenopus Sµ DNA was moved to the genome of mice, it is able to support substantial CSR when it is used to replace the murine Sγ1 region. Here, we tested both the 2kb repetitive portion and the 4.6 kb full-length portions of the Xenopus Sµ in both their natural (forward) orientation relative to the constant domain exons, as well as the opposite (reverse) orientation. Consistent with previous work, we find that the 4.6 kb full-length Sµ mediates similar levels of CSR in both the forward and reverse orientations. Whereas, the forward orientation of the 2kb portion can restore the majority of the CSR level of the 4.6 kb full-length Sµ, the reverse orientation poorly supports R-looping and no CSR. The forward orientation of the 2kb repetitive portion has more GG dinucleotides on the non-template strand than the reverse orientation. The correlation of R-loop formation with CSR efficiency, as demonstrated in the 2kb repetitive fragment of the Xenopus switch region, confirms a role played by R-looping in CSR that appears to be conserved through evolution.

Complexities due to single-stranded RNA during antibody detection of genomic rna:dna hybrids.

BACKGROUND: Long genomic R-loops in eukaryotes were first described at the immunoglobulin heavy chain locus switch regions using bisulfite sequencing and functional studies. A mouse monoclonal antibody called S9.6 has been used for immunoprecipitation (IP) to identify R-loops, based on the assumption that it is specific for RNA:DNA over other nucleic acid duplexes. However, recent work has demonstrated that a variable domain of S9.6 binds AU-rich RNA:RNA duplexes with a KD that is only 5.6-fold weaker than for RNA:DNA duplexes. Most IP protocols do not pre-clear the genomic nucleic acid with RNase A to remove free RNA. Fold back of ssRNA can readily generate RNA:RNA duplexes that may bind the S9.6 antibody, and adventitious binding of RNA may also create short RNA:DNA regions. Here we investigate whether RNase A is needed to obtain reliable IP with S9.6. FINDINGS: As our test locus, we chose the most well-documented site for kilobase-long mammalian genomic R-loops, the immunoglobulin heavy chain locus (IgH) class switch regions. The R-loops at this locus can be induced by using cytokines to stimulate transcription from germline transcript promoters. We tested IP using S9.6 with and without various RNase treatments. The RNase treatments included RNase H to destroy the RNA in an RNA:DNA duplex and RNase A to destroy single-stranded (ss) RNA to prevent it from binding S9.6 directly (as duplex RNA) and to prevent the ssRNA from annealing to the genome, resulting in adventitious RNA:DNA hybrids. We find that optimal detection of RNA:DNA duplexes requires removal of ssRNA using RNase A. Without RNase A treatment, known regions of R-loop formation containing RNA:DNA duplexes can not be reliably detected. With RNase A treatment, a signal can be detected over background, but only within a limited 2 or 3-fold range, even with a stable kilobase-long genomic R-loop. CONCLUSION: Any use of the S9.6 antibody must be preceded by RNase A treatment to remove free ssRNA that may compete for the S9.6 binding by forming RNA:RNA regions or short, transient RNA:DNA duplexes. Caution should be used when interpreting S9.6 data, and confirmation by independent structural and functional methods is essential.

The role of G-density in switch region repeats for immunoglobulin class switch recombination.

The boundaries of R-loops are well-documented at immunoglobulin heavy chain loci in mammalian B cells. Within primary B cells or B cell lines, the upstream boundaries of R-loops typically begin early in the repetitive portion of the switch regions. Most R-loops terminate within the switch repetitive zone, but the remainder can extend a few hundred base pairs further, where G-density on the non-template DNA strand gradually drops to the genome average. Whether the G-density determines how far the R-loops extend is an important question. We previously studied the role of G-clusters in initiating R-loop formation, but we did not examine the role of G-density in permitting the elongation of the R-loop, after it had initiated. Here, we vary the G-density of different portions of the switch region in a murine B cell line. We find that both class switch recombination (CSR) and R-loop formation decrease significantly when the overall G-density is reduced from 46% to 29%. Short 50 bp insertions with low G-density within switch regions do not appear to affect either CSR or R-loop elongation, whereas a longer (150 bp) insertion impairs both. These results demonstrate that G-density is an important determinant of the length over which mammalian genomic R-loops extend.

The strength of an Ig switch region is determined by its ability to drive R loop formation and its number of WGCW sites.

R loops exist at the murine IgH switch regions and possibly other locations, but their functional importance is unclear. In biochemical systems, R loop initiation requires DNA sequence regions containing clusters of G nucleotides, but cellular studies have not been done. Here, we vary the G-clustering, total switch region length, and the number of target sites (WGCW sites for the activation-induced deaminase) at synthetic switch regions in a murine B cell line to determine the effect on class switch recombination (CSR). G-clusters increase CSR regardless of their immediate proximity to the WGCW sites. This increase is accompanied by an increase in R loop formation. CSR efficiency correlates better with the absolute number of WGCW sites in the switch region rather than the total switch region length or density of WGCW sites. Thus, the overall strength of the switch region depends on G-clusters, which initiate R loop formation, and on the number of WGCW sites.

Non-homologous end joining often uses microhomology: implications for alternative end joining.

Artemis and PALF (also called APLF) appear to be among the primary nucleases involved in non-homologous end joining (NHEJ) and responsible for most nucleolytic end processing in NHEJ. About 60% of NHEJ events show an alignment of the DNA ends that use 1 or 2bp of microhomology (MH) between the two DNA termini. Thus, MH is a common feature of NHEJ. For most naturally occurring human chromosomal deletions (e.g., after oxidative damage or radiation) and translocations, such as those seen in human neoplasms and as well as inherited chromosomal structural variations, MH usage occurs at a frequency that is typical of NHEJ, and does not suggest major involvement of alternative pathways that require more extensive MH. Though we mainly focus on human NHEJ at double-strand breaks, comparison on these points to other eukaryotes, primarily S. cerevisiae, is informative.