Search

Human cell line models, including the neuronal precursor line LUHMES, are important for investigating developmental transcriptional dynamics within imprinted regions, particularly the 15q11-q13 Angelman (AS) and Prader-Willi (PWS) syndrome locus. AS results from loss of maternal UBE3A in neurons, where the paternal allele is silenced by a convergent antisense transcript UBE3A-ATS, a lncRNA that terminates at PWAR1 in non-neurons. qRT-PCR analysis confirmed the exclusive and progressive increase in UBE3A-ATS in differentiating LUHMES neurons, validating their use for studying UBE3A silencing. Genome-wide transcriptome analyses revealed changes to 11 834 genes during neuronal differentiation, including the upregulation of most genes within the 15q11-q13 locus. To identify dynamic changes in chromatin loops linked to transcriptional activity, we performed a HiChIP validated by 4C, which identified two neuron-specific CTCF loops between MAGEL2-SNRPN and PWAR1-UBE3A. To determine if allele-specific differentially methylated regions (DMR) may be associated with CTCF loop anchors, whole genome long-read nanopore sequencing was performed. We identified a paternally hypomethylated DMR near the SNRPN upstream loop anchor exclusive to neurons and a paternally hypermethylated DMR near the PWAR1 CTCF anchor exclusive to undifferentiated cells, consistent with increases in neuronal transcription. Additionally, DMRs near CTCF loop anchors were observed in both cell types, indicative of allele-specific differences in chromatin loops regulating imprinted transcription. These results provide an integrated view of the 15q11-q13 epigenetic landscape during LUHMES neuronal differentiation, underscoring the complex interplay of transcription, chromatin looping, and DNA methylation. They also provide insights for future therapeutic approaches for AS and PWS.

Automated annotations of protein functions are error-prone because of our lack of knowledge of protein functions. For example, it is often impossible to predict the correct substrate for an enzyme or a transporter. Furthermore, much of the knowledge that we do have about the functions of proteins is missing from the underlying databases. We discuss how to use interactive tools to quickly find different kinds of information relevant to a protein's function. Many of these tools are available via PaperBLAST (http://papers.genomics.lbl.gov). Combining these tools often allows us to infer a protein's function. Ideally, accurate annotations would allow us to predict a bacterium's capabilities from its genome sequence, but in practice, this remains challenging. We describe interactive tools that infer potential capabilities from a genome sequence or that search a genome to find proteins that might perform a specific function of interest. Database URL: http://papers.genomics.lbl.gov.

BackgroundPrenatal folate supplementation has been consistently associated with a reduced risk of childhood acute lymphoblastic leukemia (ALL). Previous germline genetic studies examining the one carbon (folate) metabolism pathway were limited in sample size, scope, and population diversity and led to inconclusive results.MethodsWe evaluated whether ∼2,900 single-nucleotide polymorphisms (SNP) within 46 candidate genes involved in the folate metabolism pathway influence the risk of childhood ALL, using genome-wide data from nine case-control studies in the Childhood Cancer and Leukemia International Consortium (n = 9,058 cases including 4,510 children of European ancestry, 3,018 Latinx, and 1,406 Asians, and 92,364 controls). Each study followed a standardized protocol for quality control and imputation of genome-wide data and summary statistics were meta-analyzed for all children combined and by major ancestry group using METAL software.ResultsNone of the selected SNPs reached statistical significance, overall and for major ancestry groups (using adjusted Bonferroni P-value of 5 × 10-6 and less-stringent P-value of 3.5 × 10-5 accounting for the number of "independent" SNPs). None of the 10 top (nonsignificant) SNPs and corresponding genes overlapped across ancestry groups.ConclusionsThis large meta-analysis of original data does not reveal associations between many common genetic variants in the folate metabolism pathway and childhood ALL in various ancestry groups.ImpactGenetic variants in the folate pathway alone do not appear to substantially influence childhood acute lymphoblastic leukemia risk. Other mechanisms such as gene-folate interaction, DNA methylation, or maternal genetic effects may explain the observed associations with self-reported prenatal folate intake.

CRISPR-based high-throughput genome-wide loss-of-function screens are a valuable approach to functional genetics and strain engineering. The yeast Komagataella phaffii is a host of particular interest in the biopharmaceutical industry and as a metabolic engineering host for proteins and metabolites. Here, we design and validate a highly active 6-fold coverage genome-wide sgRNA library for this biotechnologically important yeast containing 30,848 active sgRNAs targeting over 99% of its coding sequences. Conducting fitness screens in the absence of functional non-homologous end joining (NHEJ), the dominant DNA repair mechanism in K. phaffii, provides a quantitative means to assess the activity of each sgRNA in the library. This approach allows for the experimental validation of each guide's targeting activity, leading to more precise screening outcomes. We used this approach to conduct growth screens with glucose as the sole carbon source and identify essential genes. Comparative analysis of the called gene sets identified a core set of K. phaffii essential genes, many of which relate to metabolic engineering targets, including protein production, secretion, and glycosylation. The high activity, genome-wide CRISPR library developed here enables functional genomic screening in K. phaffii, applied here to gene essentiality classification, and promises to enable other genetic screens.

The American black bear, Ursus americanus, is a widespread and ecologically important species in North America. In California, the black bear plays an important role in a variety of ecosystems and serves as an important species for recreational hunting. While research suggests that the populations in California are currently healthy, continued monitoring is critical, with genomic analyses providing an important surveillance tool. Here we report a high-quality, near chromosome-level genome assembly from a U. americanus sample from California. The primary assembly has a total length of 2.5 Gb contained in 316 scaffolds, a contig N50 of 58.9 Mb, a scaffold N50 of 67.6 Mb, and a BUSCO completeness score of 96%. This U. americanus genome assembly will provide an important resource for the targeted management of black bear populations in California, with the goal of achieving an appropriate balance between the recreational value of black bears and the maintenance of viable populations. The high quality of this genome assembly will also make it a valuable resource for comparative genomic analyses among black bear populations and among bear species.

The California Pipevine, Aristolochia californica Torr., is the only endemic California species within the cosmopolitan birthwort family Aristolochiaceae. It occurs as an understory vine in riparian and chaparral areas and in forest edges and windrows. The geographic range of this plant species almost entirely overlaps with that of its major specialized herbivore, the California Pipevine Swallowtail Butterfly Battus philenor hirsuta. While this species pair is a useful, ecologically well-understood system to study co-evolution, until recently, genomic resources for both have been lacking. Here, we report a new, chromosome-level assembly of A. californica as part of the California Conservation Genomics Project (CCGP). Following the sequencing and assembly strategy of the CCGP, we used Pacific Biosciences HiFi long reads and Hi-C chromatin proximity sequencing technology to produce a de novo assembled genome. Our genome assembly, the first for any species in the genus, contains 531 scaffolds spanning 661 megabase (Mb) pairs, with a contig N50 of 6.53 Mb, a scaffold N50 of 42.2 Mb, and BUSCO complete score of 98%. In combination with the recently published B. philenor hirsuta reference genome assembly, the A. californica reference genome assembly will be a powerful tool for studying co-evolution in a rapidly changing California landscape.

PurposeThe purpose of this study was to evaluate RB1 expression and survival across ovarian carcinoma histotypes and how co-occurrence of BRCA1 or BRCA2 (BRCA) alterations and RB1 loss influences survival in tubo-ovarian high-grade serous carcinoma (HGSC).Experimental designRB1 protein expression was classified by immunohistochemistry in ovarian carcinomas of 7,436 patients from the Ovarian Tumor Tissue Analysis consortium. We examined RB1 expression and germline BRCA status in a subset of 1,134 HGSC, and related genotype to overall survival (OS), tumor-infiltrating CD8+ lymphocytes, and transcriptomic subtypes. Using CRISPR-Cas9, we deleted RB1 in HGSC cells with and without BRCA1 alterations to model co-loss with treatment response. We performed whole-genome and transcriptome data analyses on 126 patients with primary HGSC to characterize tumors with concurrent BRCA deficiency and RB1 loss.ResultsRB1 loss was associated with longer OS in HGSC but with poorer prognosis in endometrioid ovarian carcinoma. Patients with HGSC harboring both RB1 loss and pathogenic germline BRCA variants had superior OS compared with patients with either alteration alone, and their median OS was three times longer than those without pathogenic BRCA variants and retained RB1 expression (9.3 vs. 3.1 years). Enhanced sensitivity to cisplatin and paclitaxel was seen in BRCA1-altered cells with RB1 knockout. Combined RB1 loss and BRCA deficiency correlated with transcriptional markers of enhanced IFN response, cell-cycle deregulation, and reduced epithelial-mesenchymal transition. CD8+ lymphocytes were most prevalent in BRCA-deficient HGSC with co-loss of RB1.ConclusionsCo-occurrence of RB1 loss and BRCA deficiency was associated with exceptionally long survival in patients with HGSC, potentially due to better treatment response and immune stimulation.

SummaryHarmonizing variant indexing and allele assignments across datasets is crucial for data integrity in cross-dataset studies such as multi-cohort genome-wide association studies, meta-analyses, and the development, validation, and application of polygenic risk scores. Ensuring this indexing and allele consistency is a laborious, time-consuming, and error-prone process requiring a certain degree of computational proficiency. Here, we introduce GRIEVOUS, a command-line tool for cross-dataset variant homogenization. By means of an internal database and a custom indexing methodology, GRIEVOUS identifies, formats, and aligns all biallelic single nucleotide polymorphisms (SNPs) across all summary statistic and genotype files of interest. Upon completion of dataset harmonization, GRIEVOUS can also be used to extract the maximal set of biallelic SNPs common to all datasets.Availability and implementationGRIEVOUS and all supporting documentation and tutorials can be found at https://github.com/jvtalwar/GRIEVOUS. It is freely and publicly available under the MIT license and can be installed via pip.

Next-generation sequencing (NGS) has proven clinical utility on disease management and serves as an important tool for genomic surveillance. Currently, hurdles surrounding its implementation, namely the complex and demanding analytical workflows, have impeded its widespread use in many laboratories. To address this challenge, the UCLA Molecular Microbiology and Pathogen Genomics Laboratory evaluated the performance of the Tecan MagicPrep NGS system, a commercial automated solution for library preparation for clinical whole-genome sequencing assays, against the Illumina Nextera DNA Flex Library Prep. Using 35 unique organisms (28 bacteria and 7 fungi) for various clinical applications, including microbial identification and genomic characterization, we compared the quantity and quality of the prepared libraries and the resulting sequences, and concordance of the overall results. We also assessed the impact of its implementation on laboratory workflow. The MagicPrep NGS produced higher library concentrations with smaller sizes, and correspondingly, higher molarity. Quality metrics of the sequences, however, demonstrated no significant impact on the overall results, producing 100% concordance with the reference method. Importantly, workflow analysis showed 5 hours less hands-on time per run with more flexibility. This evaluation study indicates that performance of the MagicPrep NGS is comparable to the Nextera DNA Flex with the added benefit of improving workflow efficiency and reducing labor for performing routine clinical microbial whole-genome sequencing tests.

Delay discounting refers to the behavioral tendency to devalue rewards as a function of their delay in receipt. Heightened delay discounting has been associated with substance use disorders and multiple co-occurring psychopathologies. Human and animal genetic studies have established that delay discounting is heritable, but only a few associated genes have been identified. We aimed to identify novel genetic loci associated with delay discounting through a genome-wide association study (GWAS) using Heterogeneous Stock (HS) rats, a genetically diverse outbred population derived from eight inbred founder strains. We assessed delay discounting in 650 male and female HS rats using an adjusting amount procedure in which rats chose between smaller immediate sucrose rewards or a larger reward at various delays. Preference switch points were calculated and both exponential and hyperbolic functions were fitted to these indifference points. Area under the curve (AUC) and the discounting parameter k of both functions were used as delay discounting measures. GWAS for AUC, exponential k, and one indifference point identified significant loci on chromosomes 20 and 14. The gene Slc35f1, which encodes a member of the solute carrier family, was the sole gene within the chromosome 20 locus. That locus also contained an eQTL for Slc35f1, suggesting that heritable differences in the expression might be responsible for the association with behavior. Adgrl3, which encodes a latrophilin subfamily G-protein coupled receptor, was the sole gene within the chromosome 14 locus. These findings implicate novel genes in delay discounting and highlight the need for further exploration.

Schizophyllum commune is a mushroom-forming fungus notable for its distinctive fruiting bodies with split gills. It is used as a model organism to study mushroom development, lignocellulose degradation and mating type loci. It is a hypervariable species with considerable genetic and phenotypic diversity between the strains. In this study, we systematically phenotyped 16 dikaryotic strains for aspects of mushroom development and 18 monokaryotic strains for lignocellulose degradation. There was considerable heterogeneity among the strains regarding these phenotypes. The majority of the strains developed mushrooms with varying morphologies, although some strains only grew vegetatively under the tested conditions. Growth on various carbon sources showed strain-specific profiles. The genomes of seven monokaryotic strains were sequenced and analyzed together with six previously published genome sequences. Moreover, the related species Schizophyllum fasciatum was sequenced. Although there was considerable genetic variation between the genome assemblies, the genes related to mushroom formation and lignocellulose degradation were well conserved. These sequenced genomes, in combination with the high phenotypic diversity, will provide a solid basis for functional genomics analyses of the strains of S. commune.

Genome analysis of individuals affected by retinitis pigmentosa (RP) identified two rare nucleotide substitutions at the same genomic location on chromosome 11 (g.61392563 [GRCh38]), 69 base pairs upstream of the start codon of the ciliopathy gene TMEM216 (c.-69G>A, c.-69G>T [GenBank: NM_001173991.3]), in individuals of South Asian and African ancestry, respectively. Genotypes included 71 homozygotes and 3 mixed heterozygotes in trans with a predicted loss-of-function allele. Haplotype analysis showed single-nucleotide variants (SNVs) common across families, suggesting ancestral alleles within the two distinct ethnic populations. Clinical phenotype analysis of 62 available individuals from 49 families indicated a similar clinical presentation with night blindness in the first decade and progressive peripheral field loss thereafter. No evident systemic ciliopathy features were noted. Functional characterization of these variants by luciferase reporter gene assay showed reduced promotor activity. Nanopore sequencing confirmed the lower transcription of the TMEM216 c.-69G>T allele in blood-derived RNA from a heterozygous carrier, and reduced expression was further recapitulated by qPCR, using both leukocytes-derived RNA of c.-69G>T homozygotes and total RNA from genome-edited hTERT-RPE1 cells carrying homozygous TMEM216 c.-69G>A. In conclusion, these variants explain a significant proportion of unsolved cases, specifically in individuals of African ancestry, suggesting that reduced TMEM216 expression might lead to abnormal ciliogenesis and photoreceptor degeneration.

Worldwide trends to delay childbearing have increased parental ages at birth. Older parental age may harm offspring health, but mechanisms remain unclear. Alterations in offspring DNA methylation (DNAm) patterns could play a role as aging has been associated with methylation changes in gametes of older individuals. We meta-analyzed epigenome-wide associations of parental age with offspring blood DNAm of over 9500 newborns and 2000 children (5-10 years old) from the Pregnancy and Childhood Epigenetics consortium. In newborns, we identified 33 CpG sites in 13 loci with DNAm associated with maternal age (PFDR

Precision medicine holds great promise for improving cancer outcomes. Yet, there are large inequities in the demographics of patients from whom genomic data and models, including patient-derived xenografts (PDX), are developed and for whom treatments are optimized. In this study, we developed a genetic ancestry pipeline for the Cancer Genomics Cloud, which we used to assess the diversity of models currently available in the National Cancer Institute-supported PDX Development and Trial Centers Research Network (PDXNet). We showed that there is an under-representation of models derived from patients of non-European ancestry, consistent with other cancer model resources. We discussed these findings in the context of disparities in cancer incidence and outcomes among demographic groups in the US, as well as power analyses for biomarker discovery, to highlight the immediate need for developing models from minority populations to address cancer health equity in precision medicine. Our analyses identified key priority disparity-associated cancer types for which new models should be developed.SignificanceUnderstanding whether and how tumor genetic factors drive differences in outcomes among U.S. minority groups is critical to addressing cancer health disparities. Our findings suggest that many additional models will be necessary to understand the genome-driven sources of these disparities.

Recent studies have highlighted the essential role of RNA splicing, a key mechanism of alternative RNA processing, in establishing connections between genetic variations and disease. Genetic loci influencing RNA splicing variations show considerable influence on complex traits, possibly surpassing those affecting total gene expression. Dysregulated RNA splicing has emerged as a major potential contributor to neurological and psychiatric disorders, likely due to the exceptionally high prevalence of alternatively spliced genes in the human brain. Nevertheless, establishing direct associations between genetically altered splicing and complex traits has remained an enduring challenge. We introduce Spliced-Transcriptome-Wide Associations (SpliTWAS) to integrate alternative splicing information with genome-wide association studies to pinpoint genes linked to traits through exon splicing events. We applied SpliTWAS to two schizophrenia (SCZ) RNA-sequencing datasets, BrainGVEX and CommonMind, revealing 137 and 88 trait-associated exons (in 84 and 67 genes), respectively. Enriched biological functions in the associated gene sets converged on neuronal function and development, immune cell activation, and cellular transport, which are highly relevant to SCZ. SpliTWAS variants impacted RNA-binding protein binding sites, revealing potential disruption of RNA-protein interactions affecting splicing. We extended the probabilistic fine-mapping method FOCUS to the exon level, identifying 36 genes and 48 exons as putatively causal for SCZ. We highlight VPS45 and APOPT1, where splicing of specific exons was associated with disease risk, eluding detection by conventional gene expression analysis. Collectively, this study supports the substantial role of alternative splicing in shaping the genetic basis of SCZ, providing a valuable approach for future investigations in this area.

The contrast between the disruption of genome topology after cohesin loss and the lack of downstream gene expression changes instigates intense debates regarding the structure-function relationship between genome and gene regulation. Here, by analyzing transcriptome and chromatin accessibility at the single-cell level, we discover that, instead of dictating population-wide gene expression levels, cohesin supplies a general function to neutralize stochastic coexpression tendencies of cis-linked genes in single cells. Notably, cohesin loss induces widespread gene coactivation and chromatin co-opening tens of million bases apart in cis. Spatial genome and protein imaging reveals that cohesin prevents gene co-bursting along the chromosome and blocks spatial mixing of transcriptional hubs. Single-molecule imaging shows that cohesin confines the exploration of diverse enhancer and core promoter binding transcriptional regulators. Together, these results support that cohesin arranges nuclear topology to control gene coexpression in single cells.

IGHMBP2 is a nonessential, superfamily 1 DNA/RNA helicase that is mutated in patients with rare neuromuscular diseases SMARD1 and CMT2S. IGHMBP2 is implicated in translational and transcriptional regulation via biochemical association with ribosomal proteins, pre-rRNA processing factors, and tRNA-related species. To uncover the cellular consequences of perturbing IGHMBP2, we generated full and partial IGHMBP2 deletion K562 cell lines. Using polysome profiling and a nascent protein synthesis assay, we found that IGHMBP2 deletion modestly reduces global translation. We performed Ribo-seq and RNA-seq and identified diverse gene expression changes due to IGHMBP2 deletion, including ATF4 up-regulation. With recent studies showing the integrated stress response (ISR) can contribute to tRNA metabolism-linked neuropathies, we asked whether perturbing IGHMBP2 promotes ISR activation. We generated ATF4 reporter cell lines and found IGHMBP2 knockout cells demonstrate basal, chronic ISR activation. Our work expands upon the impact of IGHMBP2 in translation and elucidates molecular mechanisms that may link mutant IGHMBP2 to severe clinical phenotypes.

San Fernando Valley spineflower (Chorizanthe parryi var. fernandina [S. Watson] Jeps.) (Polygonaceae) is an herbaceous annual plant, endemic to California, and until rediscovered in 1999 had been thought to be extinct for almost seven decades. Historically documented at 10 locations, it currently persists at 2, separated by approximately 27 km. State listed as endangered, a description of its genetic diversity and structure is of conservation interest. After determining a lack of variation in ploidy, we examined genetic variation from samples within both populations: a common garden study for potentially adaptive genetic variation in selected growth and phenological traits and analysis of single nucleotide polymorphisms identified through restriction-site associated DNA sequencing. Both measures indicated that this highly restricted taxon nevertheless harbors substantial levels of genetic diversity and has significant between- and within-population genetic structure. Combining approaches from population genomics and common garden studies provided more insight into the patterns and basis of genetic diversity than is typical for studies of non-model species. Although local adaptation was not specifically studied (i.e., via reciprocal transplant studies), the differences determined from these two independent lines of evidence indicate that mixing gene pools between populations is not recommended at this time. Further, with significant differences revealed among subpopulations, we caution against mixing genotypes across subpopulations for the most part, and without much more evidence that this would not pose a risk of outbreeding depression. The importance of supporting pollinator health and diversity is highlighted. With genetic diversity—particularly with an annual species—being dynamic, fluctuating with the usual processes and with contributions from the soil seedbank, we recommend periodic resampling to monitor genetic diversity and structure. Climate change is anticipated to contribute to this variability.

The insights into interactions between host genetics and gut microbiome (GM) in colorectal tumor susceptibility (CTS) remains lacking. We used Collaborative Cross mouse population model to identify genetic and microbial determinants of Azoxymethane-induced CTS. We identified 4417 CTS-associated single nucleotide polymorphisms (SNPs) containing 334 genes that were transcriptionally altered in human colorectal cancers (CRCs) and consistently clustered independent human CRC cohorts into two subgroups with different prognosis. We discovered a set of genera in early-life associated with CTS and defined a 16-genus signature that accurately predicted CTS, the majority of which were correlated with human CRCs. We identified 547 SNPs associated with abundances of these genera. Mediation analysis revealed GM as mediators partially exerting the effect of SNP UNC3869242 within Duox2 on CTS. Intestine cell-specific depletion of Duox2 altered GM composition and contribution of Duox2 depletion to CTS was significantly influenced by GM. Our findings provide potential novel targets for personalized CRC prevention and treatment.

Systemic low-grade inflammation is a feature of chronic disease. C-reactive protein (CRP) is a common biomarker of inflammation and used as an indicator of disease risk; however, the role of inflammation in disease is not completely understood. Methylation is an epigenetic modification in the DNA which plays a pivotal role in gene expression. In this study we evaluated differential DNA methylation patterns associated with blood CRP level to elucidate biological pathways and genetic regulatory mechanisms to improve the understanding of chronic inflammation. The racially and ethnically diverse participants in this study were included as 50% White, 41% Black or African American, 7% Hispanic or Latino/a, and 2% Native Hawaiian, Asian American, American Indian, or Alaska Native (total n = 13,433) individuals. We replicated 113 CpG sites from 87 unique loci, of which five were novel (CADM3, NALCN, NLRC5, ZNF792, and cg03282312), across a discovery set of 1,150 CpG sites associated with CRP level (p

The complex pathology of mild traumatic brain injury (mTBI) is a main contributor to the difficulties in achieving a successful therapeutic regimen. Thyroxine (T4) administration has been shown to prevent the cognitive impairments induced by mTBI in mice but the mechanism is poorly understood. To understand the underlying mechanism, we carried out a single cell transcriptomic study to investigate the spatiotemporal effects of T4 on individual cell types in the hippocampus and frontal cortex at three post-injury stages in a mouse model of mTBI. We found that T4 treatment altered the proportions and transcriptomes of numerous cell types across tissues and timepoints, particularly oligodendrocytes, astrocytes, and microglia, which are crucial for injury repair. T4 also reversed the expression of mTBI-affected genes such as Ttr, mt-Rnr2, Ggn12, Malat1, Gnaq, and Myo3a, as well as numerous pathways such as cell/energy/iron metabolism, immune response, nervous system, and cytoskeleton-related pathways. Cell-type specific network modeling revealed that T4 mitigated select mTBI-perturbed dynamic shifts in subnetworks related to cell cycle, stress response, and RNA processing in oligodendrocytes. Cross cell-type ligand-receptor networks revealed the roles of App, Hmgb1, Fn1, and Tnf in mTBI, with the latter two ligands having been previously identified as TBI network hubs. mTBI and/or T4 signature genes were enriched for human genome-wide association study (GWAS) candidate genes for cognitive, psychiatric and neurodegenerative disorders related to mTBI. Our systems-level single cell analysis elucidated the temporal and spatial dynamic reprogramming of cell-type specific genes, pathways, and networks, as well as cell-cell communications as the mechanisms through which T4 mitigates cognitive dysfunction induced by mTBI.

Patterns of transcriptional activity are encoded in our genome through regulatory elements such as promoters or enhancers that, paradoxically, contain similar assortments of sequence-specific transcription factor (TF) binding sites1-3. Knowledge of how these sequence motifs encode multiple, often overlapping, gene expression programs is central to understanding gene regulation and how mutations in non-coding DNA manifest in disease4,5. Here, by studying gene regulation from the perspective of individual transcription start sites (TSSs), using natural genetic variation, perturbation of endogenous TF protein levels and massively parallel analysis of natural and synthetic regulatory elements, we show that the effect of TF binding on transcription initiation is position dependent. Analysing TF-binding-site occurrences relative to the TSS, we identified several motifs with highly preferential positioning. We show that these patterns are a combination of a TF's distinct functional profiles-many TFs, including canonical activators such as NRF1, NFY and Sp1, activate or repress transcription initiation depending on their precise position relative to the TSS. As such, TFs and their spacing collectively guide the site and frequency of transcription initiation. More broadly, these findings reveal how similar assortments of TF binding sites can generate distinct gene regulatory outcomes depending on their spatial configuration and how DNA sequence polymorphisms may contribute to transcription variation and disease and underscore a critical role for TSS data in decoding the regulatory information of our genome.

Pigmentation is an excellent trait to examine patterns of evolutionary change because it is often under natural selection. Benthic and limnetic threespine stickleback (Gasterosteus aculeatus) exhibit distinct pigmentation phenotypes, likely an adaptation to occupation of divergent niches. The genetic architecture of pigmentation in vertebrates appears to be complex. Prior QTL mapping of threespine stickleback pigmentation phenotypes has identified several candidate loci. However-relative to other morphological phenotypes (e.g., spines or lateral plates)-the genetic architecture of threespine stickleback pigmentation remains understudied. Here, we performed QTL mapping for two melanic pigmentation traits (melanophore density and lateral barring) using benthic-limnetic F2 crosses. The two traits mapped to different chromosomes, suggesting a distinct genetic basis. The resulting QTLs were additive, but explained a relatively small fraction of the total variance (~6%). QTLs maps differed by F1 family, suggesting variation in genetic architecture or ability to detect loci of small effect. Functional analysis identified enriched pathways for candidate loci. Several of the resulting candidate loci for pigmentation, including three loci in enriched pathways (bco1, sulf1, and tyms) have been previously indicated to affect pigmentation in other vertebrates. These findings add to a growing body of evidence suggesting pigmentation is often polygenic.

Glioblastoma (GBM) is the deadliest form of primary brain tumor with limited treatment options. Recent studies have profiled GBM tumor heterogeneity, revealing numerous axes of variation that explain the molecular and spatial features of the tumor. Here, we seek to bridge descriptive characterization of GBM cell type heterogeneity with the functional role of individual populations within the tumor. Our lens leverages a gene program-centric meta-atlas of published transcriptomic studies to identify commonalities between diverse tumors and cell types in order to decipher the mechanisms that drive them. This approach led to the discovery of a tumor-derived stem cell population with mixed vascular and neural stem cell features, termed a neurovascular progenitor (NVP). Following in situ validation and molecular characterization of NVP cells in GBM patient samples, we characterized their function in vivo. Genetic depletion of NVP cells resulted in altered tumor cell composition, fewer cycling cells, and extended survival, underscoring their critical functional role. Clonal analysis of primary patient tumors in a human organoid tumor transplantation system demonstrated that the NVP has dual potency, generating both neuronal and vascular tumor cells. Although NVP cells comprise a small fraction of the tumor, these clonal analyses demonstrated that they strongly contribute to the total number of cycling cells in the tumor and generate a defined subset of the whole tumor. This study represents a paradigm by which cell type-specific interrogation of tumor populations can be used to study functional heterogeneity and therapeutically targetable vulnerabilities of GBM.

One of the justifiable criticisms of human genetic studies is the underrepresentation of participants from diverse populations. Lack of inclusion must be addressed at-scale to identify causal disease factors and understand the genetic causes of health disparities. We present genome-wide associations for 2068 traits from 635,969 participants in the Department of Veterans Affairs Million Veteran Program, a longitudinal study of diverse United States Veterans. Systematic analysis revealed 13,672 genomic risk loci; 1608 were only significant after including non-European populations. Fine-mapping identified causal variants at 6318 signals across 613 traits. One-third (n = 2069) were identified in participants from non-European populations. This reveals a broadly similar genetic architecture across populations, highlights genetic insights gained from underrepresented groups, and presents an extensive atlas of genetic associations.

Drosophila prolongata is a member of the melanogaster species group and rhopaloa subgroup native to the subtropical highlands of southeast Asia. This species exhibits an array of recently evolved male-specific morphological, physiological, and behavioral traits that distinguish it from its closest relatives, making it an attractive model for studying the evolution of sexual dimorphism and testing theories of sexual selection. The lack of genomic resources has impeded the dissection of the molecular basis of sex-specific development and behavior in this species. To address this, we assembled the genome of D. prolongata using long-read sequencing and Hi-C scaffolding, resulting in a highly complete and contiguous (scaffold N50 2.2Mb) genome assembly of 220Mb. The repetitive content of the genome is 24.6%, the plurality of which are LTR retrotransposons (33.2%). Annotations based on RNA-seq data and homology to related species revealed a total of 19,330 genes, of which 16,170 are protein-coding. The assembly includes 98.5% of Diptera BUSCO genes, including 93.8% present as a single copy. Despite some likely regional duplications, the completeness of this genome suggests that it can be readily used for gene expression, GWAS, and other genomic analyses.

The evaluation of epigenetic features such as DNA methylation is becoming increasingly important in many biochemical processes like gene expression and transcription as well as in several diseases like schizophrenia or diabetes. Here, we report that self-assembled nanomechanical resonators entirely composed of DNA molecules can be used to explore gross changes in DNA methylation levels (0-25-50%), while careful control of tensile stress is needed to reduce the variability of resonance frequency for rigorous quantification. The effect of the tensile stress retained by the suspended DNA nanoresonators on the application of the technique is extensively explored using a combination of laser Doppler vibrometry and atomic force spectroscopy. DNA nanoresonators are real-time, label-free sensors and could avoid chemical functionalization and sample amplification. Therefore, they may represent in the future a key complementary routine tool for global DNA methylation analysis needed to evaluate the consequences of environmental stresses on the human genome.

The jaguar (Panthera onca) is the largest living cat species native to the Americas and one of few large American carnivorans to have survived into the Holocene. However, the extent to which jaguar diversity declined during the end-Pleistocene extinction event remains unclear. For example, Pleistocene jaguar fossils from North America are notably larger than the average extant jaguar, leading to hypotheses that jaguars from this continent represent a now-extinct subspecies (Panthera onca augusta) or species (Panthera augusta). Here, we used a hybridization capture approach to recover an ancient mitochondrial genome from a large, late Pleistocene jaguar from Kingston Saltpeter Cave, Georgia, United States, which we sequenced to 26-fold coverage. We then estimated the evolutionary relationship between the ancient jaguar mitogenome and those from other extinct and living large felids, including multiple jaguars sampled across the species' current range. The ancient mitogenome falls within the diversity of living jaguars. All sampled jaguar mitogenomes share a common mitochondrial ancestor ~400 thousand years ago, indicating that the lineage represented by the ancient specimen dispersed into North America from the south at least once during the late Pleistocene. While genomic data from additional and older specimens will continue to improve understanding of Pleistocene jaguar diversity in the Americas, our results suggest that this specimen falls within the variation of extant jaguars despite the relatively larger size and geographic location and does not represent a distinct taxon.

Up to 80% of Parkinson's disease patients develop dementia, but time to dementia varies widely from motor symptom onset. Dementia with Lewy bodies presents with clinical features similar to Parkinson's disease dementia, but cognitive impairment precedes or coincides with motor onset. It remains controversial whether dementia with Lewy bodies and Parkinson's disease dementia are distinct conditions or represent part of a disease spectrum. The biological mechanisms underlying disease heterogeneity, in particular the development of dementia, remain poorly understood, but will likely be the key to understanding disease pathways and, ultimately, therapy development. Previous genome-wide association studies in Parkinson's disease and dementia with Lewy bodies/Parkinson's disease dementia have identified risk loci differentiating patients from controls. We collated data for 7804 patients of European ancestry from Tracking Parkinson's, The Oxford Discovery Cohort, and Accelerating Medicine Partnership-Parkinson's Disease Initiative. We conducted a discrete phenotype genome-wide association study comparing Lewy body diseases with and without dementia to decode disease heterogeneity by investigating the genetic drivers of dementia in Lewy body diseases. We found that risk allele rs429358 tagging APOEe4 increases the odds of developing dementia, and that rs7668531 near the MMRN1 and SNCA-AS1 genes and an intronic variant rs17442721 tagging LRRK2 G2019S on chromosome 12 are protective against dementia. These results should be validated in autopsy-confirmed cases in future studies.

Transcription factors can promote gene expression through activation domains. Whole-genome screens have systematically mapped activation domains in transcription factors but not in non-transcription factor proteins (e.g., chromatin regulators and coactivators). To fill this knowledge gap, we employed the activation domain predictor PADDLE to analyze the proteomes of Arabidopsis thaliana and Saccharomyces cerevisiae. We screened 18,000 predicted activation domains from >800 non-transcription factor genes in both species, confirming that 89% of candidate proteins contain active fragments. Our work enables the annotation of hundreds of nuclear proteins as putative coactivators, many of which have never been ascribed any function in plants. Analysis of peptide sequence compositions reveals how the distribution of key amino acids dictates activity. Finally, we validated short, "universal" activation domains with comparable performance to state-of-the-art activation domains used for genome engineering. Our approach enables the genome-wide discovery and annotation of activation domains that can function across diverse eukaryotes.

The early 2-factor (E2F) family of transcription factors, including E2F1 through 8, plays a critical role in apoptosis, metabolism, proliferation, and angiogenesis within glioblastoma (GBM). However, the specific functions of E2F transcription factors (E2Fs) and their impact on the malignancy of Bevacizumab (BVZ)-responsive GBM subtypes remain unclear. This study used data from The Cancer Genome Atlas (TCGA), Chinese Glioma Genome Atlas (CGGA), European Molecular Biology Laboratory's European Bioinformatics Institute (EMBL-EBI), and Gene Expression Omnibus (GEO) to explore the impact of eight E2F family members on the clinical characteristics of BVZ-responsive GBM subtypes and possible mechanisms of recurrence after BVZ treatment. Using machine learning algorithms, including TreeBagger and deep neural networks, we systematically predicted and validated GBM patient survival terms based on the expression profiles of E2Fs across BVZ-responsive GBM subtypes. Our bioinformatics analyses suggested that a significant increase in E2F8 post-BVZ treatment may enhance the function of angiogenesis and stem cell proliferation, implicating this factor as a candidate mechanism of GBM recurrence after treatment. In addition, BVZ treatment in unresponsive GBM patients may potentially worsen disease progression. These insights underscore that E2F family members play important roles in GBM malignancy and BVZ treatment response, highlighting their potential as prognostic biomarkers, therapeutic targets, and recommending precision BVZ treatment to individual GBM patients.

Glioblastoma (GBM) is the most common malignant primary brain tumor and remains incurable. Previous work has shown that systemic administration of Decitabine (DAC) induces sufficient expression of cancer-testis antigens (CTA) in GBM for targeting by adoptive T-cell therapy in vivo. However, the mechanisms by which DAC enhances immunogenicity in GBM remain to be elucidated. Using New York esophageal squamous cell carcinoma 1 (NY-ESO-1) as a representative inducible CTA, we demonstrate in patient tissue, immortalized glioma cells, and primary patient-derived gliomaspheres that basal CTA expression is restricted by promoter hypermethylation in gliomas. DAC treatment of glioma cells specifically inhibits DNA methylation silencing to render NY-ESO-1 and other CTA into inducible tumor antigens at single-cell resolution. Functionally, NY-ESO-1 T-cell receptor-engineered effector cell targeting of DAC-induced antigen in primary glioma cells promotes specific and polyfunctional T-cell cytokine profiles. In addition to induction of CTA, DAC concomitantly reactivates tumor-intrinsic human endogenous retroviruses, interferon response signatures, and MHC-I. Overall, we demonstrate that DAC induces targetable tumor antigen and enhances T-cell functionality against GBM, ultimately contributing to the improvement of targeted immune therapies in glioma.SignificanceThis study dissects the tumor-intrinsic epigenetic and transcriptional mechanisms underlying enhanced T-cell functionality targeting decitabine-induced cancer-testis antigens in glioma. Our findings demonstrate concomitant induction of tumor antigens, reactivation of human endogenous retroviruses, and stimulation of interferon signaling as a mechanistic rationale to epigenetically prime human gliomas to immunotherapeutic targeting.

Consortia of multicellular magnetotactic bacteria (MMB) are currently the only known example of bacteria without a unicellular stage in their life cycle. Because of their recalcitrance to cultivation, most previous studies of MMB have been limited to microscopic observations. To study the biology of these unique organisms in more detail, we use multiple culture-independent approaches to analyze the genomics and physiology of MMB consortia at single-cell resolution. We separately sequenced the metagenomes of 22 individual MMB consortia, representing 8 new species, and quantified the genetic diversity within each MMB consortium. This revealed that, counter to conventional views, cells within MMB consortia are not clonal. Single consortia metagenomes were then used to reconstruct the species-specific metabolic potential and infer the physiological capabilities of MMB. To validate genomic predictions, we performed stable isotope probing (SIP) experiments and interrogated MMB consortia using fluorescence in situ hybridization (FISH) combined with nanoscale secondary ion mass spectrometry (NanoSIMS). By coupling FISH with bioorthogonal noncanonical amino acid tagging (BONCAT), we explored their in situ activity as well as variation of protein synthesis within cells. We demonstrate that MMB consortia are mixotrophic sulfate reducers and that they exhibit metabolic differentiation between individual cells, suggesting that MMB consortia are more complex than previously thought. These findings expand our understanding of MMB diversity, ecology, genomics, and physiology, as well as offer insights into the mechanisms underpinning the multicellular nature of their unique lifestyle.

Plant-microbe symbioses require intense interaction and genetic coordination to successfully establish in specific cell types of the host and symbiont. Traditional RNA-seq methodologies lack the cellular resolution to fully capture these complexities, but single-cell and spatial transcriptomics (ST) are now allowing scientists to probe symbiotic interactions at an unprecedented level of detail. Here, we discuss the advantages that novel spatial and single-cell transcriptomic technologies provide in studying plant-microbe endosymbioses and highlight key recent studies. Finally, we consider the remaining limitations of applying these approaches to symbiosis research, which are mainly related to the simultaneous capture of both plant and microbial transcripts within the same cells.

IntroductionEmerging evidence links changes in the gut microbiome to late-onset Alzheimer's disease (LOAD), necessitating examination of AD mouse models with consideration of the microbiome.MethodsWe used shotgun metagenomics and untargeted metabolomics to study the human amyloid beta knock-in (hAβ-KI) murine model for LOAD compared to both wild-type (WT) mice and a model for early-onset AD (3xTg-AD).ResultsEighteen-month female (but not male) hAβ-KI microbiomes were distinct from WT microbiomes, with AD genotype accounting for 18% of the variance by permutational multivariate analysis of variance (PERMANOVA). Metabolomic diversity differences were observed in females, however no individual metabolites were differentially abundant. hAβ-KI mice microbiomes were distinguishable from 3xTg-AD animals (81% accuracy by random forest modeling), with separation primarily driven by Romboutsia ilealis and Turicibacter species. Microbiomes were highly cage specific, with cage assignment accounting for more than 40% of the PERMANOVA variance between the groups.DiscussionThese findings highlight a sex-dependent variation in the microbiomes of hAβ-KI mice and underscore the importance of considering the microbiome when designing studies that use murine models for AD.HighlightsMicrobial diversity and the abundance of several species differed in human amyloid beta knock-in (hAβ-KI) females but not males. Correlations to Alzheimer's disease (AD) genotype were stronger for the microbiome than the metabolome. Microbiomes from hAβ-KI mice were distinct from 3xTg-AD mice. Cage effects accounted for most of the variance in the microbiome and metabolome.

The Mec1/ATR kinase is crucial for genome stability, yet the mechanism by which it prevents gross chromosomal rearrangements (GCRs) remains unknown. Here we find that in cells with deficient Mec1 signaling, GCRs accumulate due to the deregulation of multiple steps in homologous recombination (HR). Mec1 primarily suppresses GCRs through its role in activating the canonical checkpoint kinase Rad53, which ensures the proper control of DNA end resection. Upon loss of Rad53 signaling and resection control, Mec1 becomes hyperactivated and triggers a salvage pathway in which the Sgs1 helicase is recruited to sites of DNA lesions via the 911-Dpb11 scaffolds and phosphorylated by Mec1 to favor heteroduplex rejection and limit HR-driven GCR accumulation. Fusing an ssDNA recognition domain to Sgs1 bypasses the requirement of Mec1 signaling for GCR suppression and nearly eliminates D-loop formation, thus preventing non-allelic recombination events. We propose that Mec1 regulates multiple steps of HR to prevent GCRs while ensuring balanced HR usage when needed for promoting tolerance to replication stress.

IntroductionAltered immune signatures are emerging as a central theme in neurodegenerative disease, yet little is known about immune responses in early-onset Alzheimer's disease (EOAD).MethodsWe examined single-cell RNA-sequencing (scRNA-seq) data from peripheral blood mononuclear cells (PBMCs) and droplet digital polymerase chain reaction (ddPCR) data from CD4 T cells from participants with EOAD and clinically normal controls.ResultsWe analyzed PBMCs from 16 individuals by scRNA-seq and discovered increased interferon signaling-associated gene (ISAG) expression and striking expansion of antiviral-like ISAGhi T cells in EOAD. Isolating CD4 T cells from 19 individuals, including four cases analyzed by scRNA-seq, we confirmed increased expression of ISAGhi marker genes. Publicly available cerebrospinal fluid leukocyte scRNA-seq data from late-onset mild cognitive impairment and AD also revealed increased expression of interferon-response genes.DiscussionAntiviral-like ISAGhi T cells are expanded in EOAD. Additional research into these cells and the role of heightened peripheral IFN signaling in neurodegeneration is warranted.HighlightsInterferon-responsive T cells expanded in early-onset Alzheimer's disease (AD). Increased interferon-associated gene expression present in early- and late-onset AD. Peripheral immune changes in T and NK cells driven by females with early-onset AD.

Accurate diagnosis of lung cancer is important for treatment decision-making. Tumor biopsy and histologic examination are the standard for determining histologic lung cancer subtypes. Liquid biopsy, particularly cell-free DNA (cfDNA), has recently shown promising results in cancer detection and classification. In this study, we investigate the potential of cfDNA methylome for the noninvasive classification of lung cancer histologic subtypes. We focused on the two most prevalent lung cancer subtypes, lung adenocarcinoma and lung squamous cell carcinoma. Using a fragment-based marker discovery approach, we identified robust subtype-specific methylation markers from tumor samples. These markers were successfully validated in independent cohorts and associated with subtype-specific transcriptional activity. Leveraging these markers, we constructed a subtype classification model using cfDNA methylation profiles, achieving an AUC of 0.808 in cross-validation and an AUC of 0.747 in the independent validation. Tumor copy-number alterations inferred from cfDNA methylome analysis revealed potential for treatment selection. In summary, our study demonstrates the potential of cfDNA methylome analysis for noninvasive lung cancer subtyping, offering insights for cancer monitoring and early detection.SignificanceThis study explores the use of cfDNA methylomes for the classification of lung cancer subtypes, vital for effective treatment. By identifying specific methylation markers in tumor tissues, we developed a robust classification model achieving high accuracy for noninvasive subtype detection. This cfDNA methylome approach offers promising avenues for early detection and monitoring.

The high heritability of amyotrophic lateral sclerosis (ALS) contrasts with its low molecular diagnosis rate post-genetic testing, pointing to potential undiscovered genetic factors. To aid the exploration of these factors, we introduced EpiOut, an algorithm to identify chromatin accessibility outliers that are regions exhibiting divergent accessibility from the population baseline in a single or few samples. Annotation of accessible regions with histone chromatin immunoprecipitation sequencing and Hi-C indicates that outliers are concentrated in functional loci, especially among promoters interacting with active enhancers. Across different omics levels, outliers are robustly replicated, and chromatin accessibility outliers are reliable predictors of gene expression outliers and aberrant protein levels. When promoter accessibility does not align with gene expression, our results indicate that molecular aberrations are more likely to be linked to post-transcriptional regulation rather than transcriptional regulation. Our findings demonstrate that the outlier detection paradigm can uncover dysregulated regions in rare diseases. EpiOut is available at github.com/uci-cbcl/EpiOut.

IntroductionIdentifying reliable biomarkers that reflect cancer survivorship symptoms remains a challenge for researchers. DNA methylation (DNAm) measurements reflecting epigenetic changes caused by anti-cancer therapy may provide needed insights. Given lack of consensus describing utilization of DNAm data to predict survivorship issues, a review evaluating the current landscape is warranted.ObjectiveProvide an overview of current studies examining associations of DNAm with survivorship burdens in cancer survivors.MethodsA literature review was conducted including studies if they focused on cohorts of cancer survivors, utilized peripheral blood cell DNAm data, and evaluated the associations of DNAm and survivorship issues.ResultsA total of 22 studies were identified, with majority focused on breast (n = 7) or childhood cancer (n = 9) survivors, and half studies included less than 100 patients (n = 11). Survivorship issues evaluated included those related to neurocognition (n = 5), psychiatric health (n = 3), general wellness (n = 9), chronic conditions (n = 5), and treatment specific toxicities (n = 4). Studies evaluated epigenetic age metrics (n = 10) and DNAm levels at individual CpG sites or regions (n = 12) for their associations with survivorship issues in cancer survivors along with relevant confounding factors. Significant associations of measured DNAm in the peripheral blood samples of cancer survivors and survivorship issues were identified.Discussion/conclusionStudies utilizing epigenetic age metrics and differential methylation analysis demonstrated significant associations of DNAm measurements with survivorship burdens. Associations were observed encompassing diverse survivorship outcomes and timeframes relative to anti-cancer therapy initiation. These findings underscore the potential of these measurements as useful biomarkers in survivorship care and research.

Whole-genome bisulfite sequencing (BS-Seq) measures cytosine methylation changes at single-base resolution and can be used to profile cell-free DNA (cfDNA). In plasma, ultrashort single-stranded cfDNA (uscfDNA, ∼50 nt) has been identified together with 167 bp double-stranded mononucleosomal cell-free DNA (mncfDNA). However, the methylation profile of uscfDNA has not been described. Conventional BS-Seq workflows may not be helpful because bisulfite conversion degrades larger DNA into smaller fragments, leading to erroneous categorization as uscfDNA. We describe the '5mCAdpBS-Seq' workflow in which pre-methylated 5mC (5-methylcytosine) single-stranded adapters are ligated to heat-denatured cfDNA before bisulfite conversion. This method retains only DNA fragments that are unaltered by bisulfite treatment, resulting in less biased uscfDNA methylation analysis. Using 5mCAdpBS-Seq, uscfDNA had lower levels of DNA methylation (∼15%) compared to mncfDNA and was enriched in promoters and CpG islands. Hypomethylated uscfDNA fragments were enriched in upstream transcription start sites (TSSs), and the intensity of enrichment was correlated with expressed genes of hemopoietic cells. Using tissue-of-origin deconvolution, we inferred that uscfDNA is derived primarily from eosinophils, neutrophils, and monocytes. As proof-of-principle, we show that characteristics of the methylation profile of uscfDNA can distinguish non-small cell lung carcinoma from non-cancer samples. The 5mCAdpBS-Seq workflow is recommended for any cfDNA methylation-based investigations.

Juvenile dermatomyositis (JDM) is one of several childhood-onset autoimmune disorders characterized by a type I IFN response and autoantibodies. Treatment options are limited due to an incomplete understanding of how the disease emerges from dysregulated cell states across the immune system. We therefore investigated the blood of patients with JDM at different stages of disease activity using single-cell transcriptomics paired with surface protein expression. By immunophenotyping peripheral blood mononuclear cells, we observed skewing of the B cell compartment toward an immature naive state as a hallmark of JDM at diagnosis. Furthermore, we find that these changes in B cells are paralleled by T cell signatures suggestive of Th2-mediated inflammation that persist despite disease quiescence. We applied network analysis to reveal that hyperactivation of the type I IFN response in all immune populations is coordinated with previously masked cell states including dysfunctional protein processing in CD4+ T cells and regulation of cell death programming in NK cells, CD8+ T cells, and γδ T cells. Together, these findings unveil the coordinated immune dysregulation underpinning JDM and provide insight into strategies for restoring balance in immune function.

Acute lymphoblastic leukemia (ALL) is the most common cancer in children, yet few environmental risk factors have been identified. We previously found an association between early-life tobacco smoke exposure and frequency of somatic deletions of 8 leukemia driver genes among childhood ALL patients in the California Childhood Leukemia Study. To expand analysis genome-wide and examine potential mechanisms, we conducted tumor whole-genome sequencing in 35 ALL patients, including 18 with high prenatal tobacco exposure and 17 with low exposure as determined by established epigenetic biomarkers. High tobacco exposure patients had significantly more structural variants (P < .001) and deletions (P = .001) genome-wide than low exposure patients. Investigation of off-target RAG recombination revealed that 41% of deletions in the high tobacco exposure patients were putatively RAG-mediated (full RAG motif identified at one or both breakpoints) compared with only 21% in the low exposure group (P = .001). In a multilevel model, deletions in high tobacco exposure patients were 2.44-fold (95% CI:1.13-5.38) more likely to be putatively RAG-mediated than deletions in low exposure patients. No point mutational signatures were associated with prenatal tobacco exposure. Our findings suggest that early-life tobacco smoke exposure may promote leukemogenesis by driving development of somatic deletions in pre-leukemic lymphocytes via off-target RAG recombination.

Some gene polymorphisms can lead to monogenic diseases, whereas other polymorphisms may confer beneficial traits. A well-characterized example is congenital erythrocytosis-the non-pathogenic hyper-production of red blood cells-that is caused by a truncated erythropoietin receptor. Here we show that Cas9-mediated genome editing in CD34+ human haematopoietic stem and progenitor cells (HSPCs) can recreate the truncated form of the erythropoietin receptor, leading to substantial increases in erythropoietic output. We also show that combining the expression of the cDNA of a truncated erythropoietin receptor with a previously reported genome-editing strategy to fully replace the HBA1 gene with an HBB transgene in HSPCs (to restore normal haemoglobin production in cells with a β-thalassaemia phenotype) gives the edited HSPCs and the healthy red blood cell phenotype a proliferative advantage. Combining knowledge of human genetics with precise genome editing to insert natural human variants into therapeutic cells may facilitate safer and more effective genome-editing therapies for patients with genetic diseases.

Cowpea is a highly drought-adapted leguminous crop with great promise for improving agricultural sustainability and food security. Here, we report analyses derived from array-based genotyping of 2,021 accessions constituting a core subset of the worlds largest cowpea collection, held at the International Institute of Tropical Agriculture (IITA) in Ibadan, Nigeria. We used this dataset to examine genetic variation and population structure in worldwide cowpea. We confirm that the primary pattern of population structure is two geographically defined subpopulations originating in West and East Africa, respectively, and that population structure is associated with shifts in phenotypic distribution. Furthermore, we establish the cowpea core collection as a resource for genome-wide association studies by mapping the genetic basis of several phenotypes, with a focus on seed coat pigmentation patterning and color. We anticipate that the genotyped IITA Cowpea Core Collection will serve as a powerful tool for mapping complex traits, facilitating the acceleration of breeding programs to enhance the resilience of this crop in the face of rapid global climate change.

Clonal hematopoiesis of indeterminate potential (CHIP), whereby somatic mutations in hematopoietic stem cells confer a selective advantage and drive clonal expansion, not only correlates with age but also confers increased risk of morbidity and mortality. Here, we leverage genetically predicted traits to identify factors that determine CHIP clonal expansion rate. We used the passenger-approximated clonal expansion rate method to quantify the clonal expansion rate for 4,370 individuals in the National Heart, Lung, and Blood Institute (NHLBI) Trans-Omics for Precision Medicine (TOPMed) cohort and calculated polygenic risk scores for DNA methylation aging, inflammation-related measures and circulating protein levels. Clonal expansion rate was significantly associated with both genetically predicted and measured epigenetic clocks. No associations were identified with inflammation-related lab values or diseases and CHIP expansion rate overall. A proteome-wide search identified predicted circulating levels of myeloid zinc finger 1 and anti-Müllerian hormone as associated with an increased CHIP clonal expansion rate and tissue inhibitor of metalloproteinase 1 and glycine N-methyltransferase as associated with decreased CHIP clonal expansion rate. Together, our findings identify epigenetic and proteomic patterns associated with the rate of hematopoietic clonal expansion.

Activating signal co-integrator complex 1 (ASCC1) acts with ASCC-ALKBH3 complex in alkylation damage responses. ASCC1 uniquely combines two evolutionarily ancient domains: nucleotide-binding K-Homology (KH) (associated with regulating splicing, transcriptional, and translation) and two-histidine phosphodiesterase (PDE; associated with hydrolysis of cyclic nucleotide phosphate bonds). Germline mutations link loss of ASCC1 function to spinal muscular atrophy with congenital bone fractures 2 (SMABF2). Herein analysis of The Cancer Genome Atlas (TCGA) suggests ASCC1 RNA overexpression in certain tumors correlates with poor survival, Signatures 29 and 3 mutations, and genetic instability markers. We determined crystal structures of Alvinella pompejana (Ap) ASCC1 and Human (Hs) PDE domain revealing high-resolution details and features conserved over 500 million years of evolution. Extending our understanding of the KH domain Gly-X-X-Gly sequence motif, we define a novel structural Helix-Clasp-Helix (HCH) nucleotide binding motif and show ASCC1 sequence-specific binding to CGCG-containing RNA. The V-shaped PDE nucleotide binding channel has two His-Φ-Ser/Thr-Φ (HXT) motifs (Φ being hydrophobic) positioned to initiate cyclic phosphate bond hydrolysis. A conserved atypical active-site histidine torsion angle implies a novel PDE substrate. Flexible active site loop and arginine-rich domain linker appear regulatory. Small-angle X-ray scattering (SAXS) revealed aligned KH-PDE RNA binding sites with limited flexibility in solution. Quantitative evolutionary bioinformatic analyses of disease and cancer-associated mutations support implied functional roles for RNA binding, phosphodiesterase activity, and regulation. Collective results inform ASCC1's roles in transactivation and alkylation damage responses, its targeting by structure-based inhibitors, and how ASCC1 mutations may impact inherited disease and cancer.

The inability to scalably and precisely measure the activity of developmental cis-regulatory elements (CREs) in multicellular systems is a bottleneck in genomics. Here we develop a dual RNA cassette that decouples the detection and quantification tasks inherent to multiplex single-cell reporter assays. The resulting measurement of reporter expression is accurate over multiple orders of magnitude, with a precision approaching the limit set by Poisson counting noise. Together with RNA barcode stabilization via circularization, these scalable single-cell quantitative expression reporters provide high-contrast readouts, analogous to classic in situ assays but entirely from sequencing. Screening >200 regions of accessible chromatin in a multicellular in vitro model of early mammalian development, we identify 13 (8 previously uncharacterized) autonomous and cell-type-specific developmental CREs. We further demonstrate that chimeric CRE pairs generate cognate two-cell-type activity profiles and assess gain- and loss-of-function multicellular expression phenotypes from CRE variants with perturbed transcription factor binding sites. Single-cell quantitative expression reporters can be applied in developmental and multicellular systems to quantitatively characterize native, perturbed and synthetic CREs at scale, with high sensitivity and at single-cell resolution.

Following publication of the original article [1], the authors reported two errors in Fig. 1B and Fig. 2B. 1. The two EP + siGLO and DS + siGLO brightfield/β-gal images from Fig. 2A were accidentally duplicated in Fig. 1B. The images in Fig. 2A are correct. The images in Fig. 1B are replaced with representative images of untransfected EP and DS cells. The incorrect and the correct figure are given below. The incorrect Fig. 1: (Figure presented.) The correct Fig. 1: (Figure presented.) The replacement of these images does not change any of the other data in this figure, nor the conclusion of this figure. 2. The authors also noticed that the panels in Fig. 2B labelled with “DS” were incorrectly labelled. The bar charts and figure legend correctly state that these images are “DS + siGLO” images. These panels are now relabelled as “DS + siGLO”. This proposed relabelling does not change any of the other data in this figure, nor the conclusion of Fig. 2B. The incorrect and the correct figure are given below. The incorrect Fig. 2: (Figure presented.) The correct Fig. 2: (Figure presented.) In both instances, these errors do not change any data associated with these figures, nor do they change the conclusions of the corresponding figures and the paper.

Many autism spectrum disorder (ASD)-associated genes act as transcriptional regulators (TRs). Chromatin immunoprecipitation sequencing (ChIP-seq) was used to identify the regulatory targets of ARID1B, BCL11A, FOXP1, TBR1, and TCF7L2, ASD-associated TRs in the developing human and mouse cortex. These TRs shared substantial overlap in the binding sites, especially within open chromatin. The overlap within a promoter region, 1-2,000 bp upstream of the transcription start site, was highly predictive of brain-expressed genes. This signature was observed in 96 out of 102 ASD-associated genes. In vitro CRISPRi against ARID1B and TBR1 delineated downstream convergent biology in mouse cortical cultures. After 8 days, NeuN+ and CALB+ cells were decreased, GFAP+ cells were increased, and transcriptomic signatures correlated with the postmortem brain samples from individuals with ASD. We suggest that functional convergence across five ASD-associated TRs leads to shared neurodevelopmental outcomes of haploinsufficient disruption.