Sex can be an important determinant of cancer phenotype, and exploring sex-biased tumor biology holds promise for identifying novel therapeutic targets and new approaches to cancer treatment. In an established isogenic murine model of glioblastoma (GBM), we discovered correlated transcriptome-wide sex differences in gene expression, H3K27ac marks, large Brd4-bound enhancer usage, and Brd4 localization to Myc and p53 genomic binding sites. These sex-biased gene expression patterns were also evident in human glioblastoma stem cells (GSCs). These observations led us to hypothesize that Brd4-bound enhancers might underlie sex differences in stem cell function and tumorigenicity in GBM. We found that male and female GBM cells exhibited sex-specific responses to pharmacological or genetic inhibition of Brd4. Brd4 knockdown or pharmacologic inhibition decreased male GBM cell clonogenicity and in vivo tumorigenesis while increasing both in female GBM cells. These results were validated in male and female patient-derived GBM cell lines. Furthermore, analysis of the Cancer Therapeutic Response Portal of human GBM samples segregated by sex revealed that male GBM cells are significantly more sensitive to BET (bromodomain and extraterminal) inhibitors than are female cells. Thus, Brd4 activity is revealed to drive sex differences in stem cell and tumorigenic phenotypes, which can be abrogated by sex-specific responses to BET inhibition. This has important implications for the clinical evaluation and use of BET inhibitors.

Simultaneous profiling of multiomic modalities within a single cell is a grand challenge for single-cell biology. While there have been impressive technical innovations demonstrating feasibility—for example, generating paired measurements of single-cell transcriptome (single-cell RNA sequencing [scRNA-seq]) and chromatin accessibility (single-cell assay for transposase-accessible chromatin using sequencing [scATAC-seq])—widespread application of joint profiling is challenging due to its experimental complexity, noise, and cost. Here, we introduce BABEL, a deep learning method that translates between the transcriptome and chromatin profiles of a single cell. Leveraging an interoperable neural network model, BABEL can predict single-cell expression directly from a cell’s scATAC-seq and vice versa after training on relevant data. This makes it possible to computationally synthesize paired multiomic measurements when only one modality is experimentally available. Across several paired single-cell ATAC and gene expression datasets in human and mouse, we validate that BABEL accurately translates between these modalities for individual cells. BABEL also generalizes well to cell types within new biological contexts not seen during training. Starting from scATAC-seq of patient-derived basal cell carcinoma (BCC), BABEL generated single-cell expression that enabled fine-grained classification of complex cell states, despite having never seen BCC data. These predictions are comparable to analyses of experimental BCC scRNA-seq data for diverse cell types related to BABEL’s training data. We further show that BABEL can incorporate additional single-cell data modalities, such as protein epitope profiling, thus enabling translation across chromatin, RNA, and protein. BABEL offers a powerful approach for data exploration and hypothesis generation.

Methylation of histone H3K9 is a hallmark of epigenetic silencing in eukaryotes. Nucleosome modifications often rely on positive feedback where enzymes are recruited by modified nucleosomes. A combination of local and global feedbacks has been proposed to account for some dynamic properties of heterochromatin, but the range at which the global feedbacks operate and the exact mode of heterochromatin propagation are not known. We investigated these questions in fission yeast. Guided by mathematical modeling, we incrementally increased the size of the mating-type region and profiled heterochromatin establishment over time. We observed exponential decays in the proportion of cells with active reporters, with rates that decreased with domain size. Establishment periods varied from a few generations in wild type to >200 generations in the longest region examined, and highly correlated silencing of two reporters located outside the nucleation center was observed. On a chromatin level, this indicates that individual regions are silenced in sudden bursts. Mathematical modeling accounts for these bursts if heterochromatic nucleosomes facilitate a deacetylation or methylation reaction at long range, in a distance-independent manner. A likely effector of three-dimensional interactions is the evolutionarily conserved Swi6^HP1 H3K9me reader, indicating the bursting behavior might be a general mode of heterochromatin propagation.

Interactions between genetic variants—epistasis—is pervasive in model systems and can profoundly impact evolutionary adaption, population disease dynamics, genetic mapping, and precision medicine efforts. In this work, we develop a model for structured polygenic epistasis, called coordinated epistasis (CE), and prove that several recent theories of genetic architecture fall under the formal umbrella of CE. Unlike standard epistasis models that assume epistasis and main effects are independent, CE captures systematic correlations between epistasis and main effects that result from pathway-level epistasis, on balance skewing the penetrance of genetic effects. To test for the existence of CE, we propose the even-odd (EO) test and prove it is calibrated in a range of realistic biological models. Applying the EO test in the UK Biobank, we find evidence of CE in 18 of 26 traits spanning disease, anthropometric, and blood categories. Finally, we extend the EO test to tissue-specific enrichment and identify several plausible tissue–trait pairs. Overall, CE is a dimension of genetic architecture that can capture structured, systemic forms of epistasis in complex human traits.

Even in well-characterized genomes, many transcripts are considered noncoding RNAs (ncRNAs) simply due to the absence of large open reading frames (ORFs). However, it is now becoming clear that many small ORFs (smORFs) produce peptides with important biological functions. In the process of characterizing the ribosome-bound transcriptome of an important cell type of the seminal fluid-producing accessory gland of Drosophila melanogaster, we detected an RNA, previously thought to be noncoding, called male-specific abdominal (msa). Notably, msa is nested in the HOX gene cluster of the Bithorax complex and is known to contain a micro-RNA within one of its introns. We find that this RNA encodes a “micropeptide” (9 or 20 amino acids, MSAmiP) that is expressed exclusively in the secondary cells of the male accessory gland, where it seems to accumulate in nuclei. Importantly, loss of function of this micropeptide causes defects in sperm competition. In addition to bringing insights into the biology of a rare cell type, this work underlines the importance of small peptides, a class of molecules that is now emerging as important actors in complex biological processes.

The synaptonemal complex is a tripartite proteinaceous ultrastructure that forms between homologous chromosomes during prophase I of meiosis in the majority of eukaryotes. It is characterized by the coordinated installation of transverse filament proteins between two lateral elements and is required for wild-type levels of crossing over and meiotic progression. We have generated null mutants of the duplicated Arabidopsis transverse filament genes zyp1a and zyp1b using a combination of T-DNA insertional mutants and targeted CRISPR/Cas mutagenesis. Cytological and genetic analysis of the zyp1 null mutants reveals loss of the obligate chiasma, an increase in recombination map length by 1.3- to 1.7-fold and a virtual absence of cross-over (CO) interference, determined by a significant increase in the number of double COs. At diplotene, the numbers of HEI10 foci, a marker for Class I interference-sensitive COs, are twofold greater in the zyp1 mutant compared to wild type. The increase in recombination in zyp1 does not appear to be due to the Class II interference-insensitive COs as chiasmata were reduced by ∼52% in msh5/zyp1 compared to msh5. These data suggest that ZYP1 limits the formation of closely spaced Class I COs in Arabidopsis. Our data indicate that installation of ZYP1 occurs at ASY1-labeled axial bridges and that loss of the protein disrupts progressive coalignment of the chromosome axes.

When addressing a genomic question, having a reliable and adequate reference genome is of utmost importance. This drives the necessity to refine and customize reference genomes (RGs). Our laboratory has recently developed a strategy, the Perfect Match Genomic Landscape (PMGL), to detect variation between genomes [K. Palacios-Flores et al.. Genetics 208, 1631–1641 (2018)]. The PMGL is precise and sensitive and, in contrast to most currently used algorithms, is nonstatistical in nature. Here we demonstrate the power of PMGL to refine and customize RGs. As a proof-of-concept, we refined different versions of the Saccharomyces cerevisiae RG. We applied the automatic PMGL pipeline to refine the genomes of microorganisms belonging to the three domains of life: the archaea Methanococcus maripaludis and Pyrococcus furiosus; the bacteria Escherichia coli, Staphylococcus aureus, and Bacillus subtilis; and the eukarya Schizosaccharomyces pombe, Aspergillus oryzae, and several strains of Saccharomyces paradoxus. We analyzed the reference genome of the virus SARS-CoV-2 and previously published viral genomes from patients’ samples with COVID-19. We performed a mutation-accumulation experiment in E. coli and show that the PMGL strategy can detect specific mutations generated at any desired step of the whole procedure. We propose that PMGL can be used as a final step for the refinement and customization of any haploid genome, independently of the strategies and algorithms used in its assembly.

Recent findings regarding nicotinamide adenine dinucleotide (NAD⁺)-capped RNAs (NAD-RNAs) indicate that prokaryotes and eukaryotes employ noncanonical RNA capping to regulate gene expression. Two methods for transcriptome-wide analysis of NAD-RNAs, NAD captureSeq and NAD tagSeq, are based on copper-catalyzed azide-alkyne cycloaddition (CuAAC) click chemistry to label NAD-RNAs. However, copper ions can fragment/degrade RNA, interfering with the analyses. Here we report development of NAD tagSeq II, which uses copper-free, strain-promoted azide-alkyne cycloaddition (SPAAC) for labeling NAD-RNAs, followed by identification of tagged RNA by single-molecule direct RNA sequencing. We used this method to compare NAD-RNA and total transcript profiles of Escherichia coli cells in the exponential and stationary phases. We identified hundreds of NAD-RNA species in E. coli and revealed genome-wide alterations of NAD-RNA profiles in the different growth phases. Although no or few NAD-RNAs were detected from some of the most highly expressed genes, the transcripts of some genes were found to be primarily NAD-RNAs. Our study suggests that NAD-RNAs play roles in linking nutrient cues with gene regulation in E. coli.

Most glioblastomas (GBMs) achieve cellular immortality by acquiring a mutation in the telomerase reverse transcriptase (TERT) promoter. TERT promoter mutations create a binding site for a GA binding protein (GABP) transcription factor complex, whose assembly at the promoter is associated with TERT reactivation and telomere maintenance. Here, we demonstrate increased binding of a specific GABPB1L-isoform–containing complex to the mutant TERT promoter. Furthermore, we find that TERT promoter mutant GBM cells, unlike wild-type cells, exhibit a critical near-term dependence on GABPB1L for proliferation, notably also posttumor establishment in vivo. Up-regulation of the protein paralogue GABPB2, which is normally expressed at very low levels, can rescue this dependence. More importantly, when combined with frontline temozolomide (TMZ) chemotherapy, inducible GABPB1L knockdown and the associated TERT reduction led to an impaired DNA damage response that resulted in profoundly reduced growth of intracranial GBM tumors. Together, these findings provide insights into the mechanism of cancer-specific TERT regulation, uncover rapid effects of GABPB1L-mediated TERT suppression in GBM maintenance, and establish GABPB1L inhibition in combination with chemotherapy as a therapeutic strategy for TERT promoter mutant GBM.

It was recently shown that the major genetic risk factor associated with becoming severely ill with COVID-19 when infected by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is inherited from Neandertals. New, larger genetic association studies now allow additional genetic risk factors to be discovered. Using data from the Genetics of Mortality in Critical Care (GenOMICC) consortium, we show that a haplotype at a region on chromosome 12 associated with requiring intensive care when infected with the virus is inherited from Neandertals. This region encodes proteins that activate enzymes that are important during infections with RNA viruses. In contrast to the previously described Neandertal haplotype that increases the risk for severe COVID-19, this Neandertal haplotype is protective against severe disease. It also differs from the risk haplotype in that it has a more moderate effect and occurs at substantial frequencies in all regions of the world outside Africa. Among ancient human genomes in western Eurasia, the frequency of the protective Neandertal haplotype may have increased between 20,000 and 10,000 y ago and again during the past 1,000 y.

Duchenne muscular dystrophy (DMD) is an X-linked recessive disorder characterized by progressive muscle degeneration and weakness due to mutations in the dystrophin gene. The symptoms of DMD share similarities with those of accelerated aging. Recently, hydrogen sulfide (H₂S) supplementation has been suggested to modulate the effects of age-related decline in muscle function, and metabolic H₂S deficiencies have been implicated in affecting muscle mass in conditions such as phenylketonuria. We therefore evaluated the use of sodium GYY4137 (NaGYY), a H₂S-releasing molecule, as a possible approach for DMD treatment. Using the dys-1(eg33) Caenorhabditis elegans DMD model, we found that NaGYY treatment (100 µM) improved movement, strength, gait, and muscle mitochondrial structure, similar to the gold-standard therapeutic treatment, prednisone (370 µM). The health improvements of either treatment required the action of the kinase JNK-1, the transcription factor SKN-1, and the NAD-dependent deacetylase SIR-2.1. The transcription factor DAF-16 was required for the health benefits of NaGYY treatment, but not prednisone treatment. AP39 (100 pM), a mitochondria-targeted H₂S compound, also improved movement and strength in the dys-1(eg33) model, further implying that these improvements are mitochondria-based. Additionally, we found a decline in total sulfide and H₂S-producing enzymes in dystrophin/utrophin knockout mice. Overall, our results suggest that H₂S deficit may contribute to DMD pathology, and rectifying/overcoming the deficit with H₂S delivery compounds has potential as a therapeutic approach to DMD treatment.

Reduced activity of insulin/insulin-like growth factor signaling (IIS) increases healthy lifespan among diverse animal species. Downstream of IIS, multiple evolutionarily conserved transcription factors (TFs) are required; however, distinct TFs are likely responsible for these effects in different tissues. Here we have asked which TFs can extend healthy lifespan within distinct cell types of the adult nervous system in Drosophila. Starting from published single-cell transcriptomic data, we report that forkhead (FKH) is endogenously expressed in neurons, whereas forkhead-box-O (FOXO) is expressed in glial cells. Accordingly, we find that neuronal FKH and glial FOXO exert independent prolongevity effects. We have further explored the role of neuronal FKH in a model of Alzheimer’s disease-associated neuronal dysfunction, where we find that increased neuronal FKH preserves behavioral function and reduces ubiquitinated protein aggregation. Finally, using transcriptomic profiling, we identify Atg17, a member of the Atg1 autophagy initiation family, as one FKH-dependent target whose neuronal overexpression is sufficient to extend healthy lifespan. Taken together, our results underscore the importance of cell type-specific mapping of TF activity to preserve healthy function with age.

Single-nucleotide variants (SNVs), pertinent to aging and disease, occur sporadically in the human genome, hence necessitating single-cell measurements. However, detection of single-cell SNVs suffers from false positives (FPs) due to intracellular single-stranded DNA damage and the process of whole-genome amplification (WGA). Here, we report a single-cell WGA method termed multiplexed end-tagging amplification of complementary strands (META-CS), which eliminates nearly all FPs by virtue of DNA complementarity, and achieved the highest accuracy thus far. We validated META-CS by sequencing kindred cells and human sperm, and applied it to other human tissues. Investigation of mature single human neurons revealed increasing SNVs with age and potentially unrepaired strand-specific oxidative guanine damage. We determined SNV frequencies along the genome in differentiated single human blood cells, and identified cell type-dependent mutational patterns for major types of lymphocytes.

Most eukaryotes possess two RecA-like recombinases (ubiquitous Rad51 and meiosis-specific Dmc1) to promote interhomolog recombination during meiosis. However, some eukaryotes have lost Dmc1. Given that mammalian and yeast Saccharomyces cerevisiae (Sc) Dmc1 have been shown to stabilize recombination intermediates containing mismatches better than Rad51, we used the Pezizomycotina filamentous fungus Trichoderma reesei to address if and how Rad51-only eukaryotes conduct interhomolog recombination in zygotes with high sequence heterogeneity. We applied multidisciplinary approaches (next- and third-generation sequencing technology, genetics, cytology, bioinformatics, biochemistry, and single-molecule biophysics) to show that T. reesei Rad51 (TrRad51) is indispensable for interhomolog recombination during meiosis and, like ScDmc1, TrRad51 possesses better mismatch tolerance than ScRad51 during homologous recombination. Our results also indicate that the ancestral TrRad51 evolved to acquire ScDmc1-like properties by creating multiple structural variations, including via amino acid residues in the L1 and L2 DNA-binding loops.

Nearly 50 hereditary diseases result from the inheritance of abnormally long repetitive DNA microsatellites. While it was originally believed that the size of inherited repeats is the key factor in disease development, it has become clear that somatic instability of these repeats throughout an individual’s lifetime strongly contributes to disease onset and progression. Importantly, somatic instability is commonly observed in terminally differentiated, postmitotic cells, such as neurons. To unravel the mechanisms of repeat instability in nondividing cells, we created an experimental system to analyze the mutability of Friedreich’s ataxia (GAA)_n repeats during chronological aging of quiescent Saccharomyces cerevisiae. Unexpectedly, we found that the predominant repeat-mediated mutation in nondividing cells is large-scale deletions encompassing parts, or the entirety, of the repeat and adjacent regions. These deletions are caused by breakage at the repeat mediated by mismatch repair (MMR) complexes MutSβ and MutLα and DNA endonuclease Rad1, followed by end-resection by Exo1 and repair of the resulting double-strand breaks (DSBs) via nonhomologous end joining. We also observed repeat-mediated gene conversions as a result of DSB repair via ectopic homologous recombination during chronological aging. Repeat expansions accrue during chronological aging as well—particularly in the absence of MMR-induced DSBs. These expansions depend on the processivity of DNA polymerase δ while being counteracted by Exo1 and MutSβ, implicating nick repair. Altogether, these findings show that the mechanisms and types of (GAA)_n repeat instability differ dramatically between dividing and nondividing cells, suggesting that distinct repeat-mediated mutations in terminally differentiated somatic cells might influence Friedreich’s ataxia pathogenesis.