Search

Introduction: Mu opioid receptor gene (OPRM1) variant rs1799971 introduces a CpG site, which may influence DNA methylation (DNAm) and opioid/pain outcomes. Objectives: In this nested analysis, we investigated both OPRM1 A118G genotype and promoter/immediate downstream blood DNAm sequencing data for associations with opioid effects and chronic postsurgical pain (CPSP) in a surgical cohort. Methods: Prospectively recruited opioid naïve patients undergoing Nuss procedure or spinal fusion with rs1799971 genotypes (Illumina arrays), DNAm (next generation enzymatic methylation sequencing at Chr6:154,039,209-154,039,803) and outcomes—opioid analgesia (integrated opioid use 1 pain over postoperative days 0 and 1 normalized to surgery type), safety—respiratory depression (RD) in high opioid use groups, and CPSP (Numerical Rating Scale .3/10 2-12 months postsurgery)—were included. Linear and logistic regression were performed to test genetic and epigenetic associations, adjusted for sociodemographics, cell types, and analgesics. Results: In this cohort (N 5 112; 15.3 6 2.0 years, 50% female, 83% White, 55% had CPSP, 13% had RD), DNAm at Chr6: 154039216-154039217 was associated with CPSP (odds ratio [OR], 1.26; 95% confidence interval [CI], 1.00-1.57; P 5 0.03), Chr6: 154039661-154039662 with acute integrated pain (b 5 220.9, 95% CI, 240.70 to 21.10, P 5 0.04), Chr6:154039520-154039521 (OR, 1.49; 95% CI, 1.09-2.03; P 5 0.01), and Chr6:154039571-154039572 (OR, 1.47; 95% CI, 1.08-2.01; P 5 0.02) with RD. Significant CpG sites were located in Repressed Polycomb chromatin states. Genotype was not associated with DNAm or outcomes. Conclusion: Our analyses support OPRM1 DNAm as predictors of acute and chronic pain/opioid outcomes in children after painful surgery. Study limitations included absent GG genotype, low sequencing coverage, and lack of correction for multiple testing.

Profound functional switch of key regulatory factors may play a major role in homeostasis and disease. Dysregulation of circadian rhythm (CR) is strongly implicated in cancer with mechanisms poorly understood. We report here that the function of REV-ERBα, a major CR regulator of the orphan nuclear receptor subfamily, is dramatically altered in tumors in both its genome binding and functional mode. Loss of CR is linked to a functional inversion of REV-ERBα from a repressor in control of CR and metabolic gene programs in normal tissues to a strong activator in different cancers. Through changing its association from NCoR/HDAC3 corepressor complex to BRD4/p300 coactivators, REV-ERBα directly activates thousands of genes including tumorigenic programs such as MAPK and PI3K-Akt signaling. Functioning as a master transcriptional activator, REV-ERBα partners with pioneer factor FOXA1 and directly stimulates a large number of signaling genes, including multiple growth factors, receptor tyrosine kinases, RASs, AKTs, and MAPKs. Moreover, elevated REV-ERBα reprograms FOXA1 to bind new targets through a BRD4-mediated increase in local chromatin accessibility. Pharmacological targeting with SR8278 diminishes the function of both REV-ERBα and FOXA1 and synergizes with BRD4 inhibitor in effective suppression of tumorigenic programs and tumor growth. Thus, our study revealed a functional inversion by a CR regulator in driving gene reprogramming as an unexpected paradigm of tumorigenesis mechanism and demonstrated a high effectiveness of therapeutic targeting such switch.

The insecticide dichlorodiphenyltrichloroethane (DDT) and its persistent metabolite, dichlorodiphenyldichloroethylene (DDE), have been associated with increased adiposity and obesity in multiple generations of rodents and humans. These lipophilic pollutants accumulate in adipose tissue and appear to decrease energy expenditure through the impairment of thermogenesis in brown adipose tissue (BAT). We hypothesized that impaired thermogenesis is due to persistent epigenetic modifications of BAT. To address this, we exposed C57BL/6 J mice to DDT or DDE from gestational day (GD) 11.5 to postnatal day (PND) 5, evaluated longitudinal body temperature, and performed reduced representation bisulfite sequencing and RNA sequencing of BAT from infant and adult offspring. Exposure to DDT or DDE reduced core body temperature in adult mice, and differential methylation at the pathway and gene level was persistent from infancy to adulthood. Furthermore, thermogenesis and biological pathways essential for thermogenic function, such as oxidative phosphorylation and mechanistic target of rapamycin kinase (mTOR) signaling, were enriched with differential methylation and RNA transcription in adult mice exposed to DDT or DDE. PAZ6 human brown preadipocytes were differentiated in the presence of DDT or DDE to understand the brown adipocyte-autonomous effect of these pollutants. In vitro exposure led to limited changes in RNA expression; however, mitochondrial membrane potential was decreased in vitro with 0.1 µM and 1 µM doses of DDT or DDE. These results demonstrate that concentrations of DDT and DDE relevant to human exposure have a significant effect on thermogenesis, the transcriptome, and DNA methylome of mouse BAT and the mitochondrial function of human brown adipocytes.

The human hippocampus and prefrontal cortex play critical roles in learning and cognition1,2, yet the dynamic molecular characteristics of their development remain enigmatic. Here we investigated the epigenomic and three-dimensional chromatin conformational reorganization during the development of the hippocampus and prefrontal cortex, using more than 53,000 joint single-nucleus profiles of chromatin conformation and DNA methylation generated by single-nucleus methyl-3C sequencing (snm3C-seq3)3. The remodelling of DNA methylation is temporally separated from chromatin conformation dynamics. Using single-cell profiling and multimodal single-molecule imaging approaches, we have found that short-range chromatin interactions are enriched in neurons, whereas long-range interactions are enriched in glial cells and non-brain tissues. We reconstructed the regulatory programs of cell-type development and differentiation, finding putatively causal common variants for schizophrenia strongly overlapping with chromatin loop-connected, cell-type-specific regulatory regions. Our data provide multimodal resources for studying gene regulatory dynamics in brain development and demonstrate that single-cell three-dimensional multi-omics is a powerful approach for dissecting neuropsychiatric risk loci.

Our ability to predict, control, or design biological function is fundamentally limited by poorly annotated gene function. This can be particularly challenging in non-model systems. Accordingly, there is motivation for new high-throughput methods for accurate functional annotation. Here, we used complementation of auxotrophs and DNA barcode sequencing (Coaux-Seq) to enable high-throughput characterization of protein function. Fragment libraries from eleven genetically diverse bacteria were tested in twenty different auxotrophic strains of Escherichia coli to identify genes that complement missing biochemical activity. We recovered 41% of expected hits, with effectiveness ranging per source genome, and observed success even with distant E. coli relatives like Bacillus subtilis and Bacteroides thetaiotaomicron. Coaux-Seq provided the first experimental validation for 53 proteins, of which 11 are less than 40% identical to an experimentally characterized protein. Among the unexpected function identified was a sulfate uptake transporter, an O-succinylhomoserine sulfhydrylase for methionine synthesis, and an aminotransferase. We also identified instances of cross-feeding wherein protein overexpression and nearby non-auxotrophic strains enabled growth. Altogether, Coaux-Seq's utility is demonstrated, with future applications in ecology, health, and engineering.

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.2 Mb) genome assembly of 220 Mb. The repetitive content of the genome is 24.6%, the plurality of which are long terminal repeats 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, genome-wide association studies (GWAS), and other genomic analyses.

Alterations in the nuclear envelope are linked to a variety of rare diseases termed laminopathies. A single amino acid substitution at position 12 (A12T) of the human nuclear envelope protein BAF (Barrier to Autointegration Factor) causes Néstor-Guillermo Progeria Syndrome (NGPS). This premature ageing condition leads to growth retardation and severe skeletal defects, but the underlying mechanisms are unknown. Here, we have generated a novel in vivo model for NGPS by modifying the baf-1 locus in C. elegans to mimic the human NGPS mutation. These baf-1(G12T) mutant worms displayed multiple phenotypes related to fertility, lifespan, and stress resistance. Importantly, nuclear morphology deteriorated faster during aging in baf-1(G12T) compared to wild-type animals, recapitulating an important hallmark of cells from progeria patients. Although localization of BAF-1(G12T) was similar to wild-type BAF-1, lamin accumulation at the nuclear envelope was reduced in mutant worms. Tissue-specific chromatin binding and transcriptome analyses showed reduced BAF-1 association in most genes deregulated by the baf-1(G12T) mutation, suggesting that altered BAF chromatin association induces NGPS phenotypes via altered gene expression.

The genomes of living lungfishes can inform on the molecular-developmental basis of the Devonian sarcopterygian fish-tetrapod transition. We de novo sequenced the genomes of the African (Protopterus annectens) and South American lungfishes (Lepidosiren paradoxa). The Lepidosiren genome (about 91 Gb, roughly 30 times the human genome) is the largest animal genome sequenced so far and more than twice the size of the Australian (Neoceratodus forsteri)1 and African2 lungfishes owing to enlarged intergenic regions and introns with high repeat content (about 90%). All lungfish genomes continue to expand as some transposable elements (TEs) are still active today. In particular, Lepidosiren's genome grew extremely fast during the past 100 million years (Myr), adding the equivalent of one human genome every 10 Myr. This massive genome expansion seems to be related to a reduction of PIWI-interacting RNAs and C2H2 zinc-finger and Krüppel-associated box (KRAB)-domain protein genes that suppress TE expansions. Although TE abundance facilitates chromosomal rearrangements, lungfish chromosomes still conservatively reflect the ur-tetrapod karyotype. Neoceratodus' limb-like fins still resemble those of their extinct relatives and remained phenotypically static for about 100 Myr. We show that the secondary loss of limb-like appendages in the Lepidosiren-Protopterus ancestor was probably due to loss of sonic hedgehog limb-specific enhancers.

Transcription factors (TFs) bind combinatorially to cis-regulatory elements, orchestrating transcriptional programs. Although studies of chromatin state and chromosomal interactions have demonstrated dynamic neurodevelopmental cis-regulatory landscapes, parallel understanding of TF interactions lags. To elucidate combinatorial TF binding driving mouse basal ganglia development, we integrated chromatin immunoprecipitation sequencing (ChIP-seq) for twelve TFs, H3K4me3-associated enhancer-promoter interactions, chromatin and gene expression data, and functional enhancer assays. We identified sets of putative regulatory elements with shared TF binding (TF-pRE modules) that orchestrate distinct processes of GABAergic neurogenesis and suppress other cell fates. The majority of pREs were bound by one or two TFs; however, a small proportion were extensively bound. These sequences had exceptional evolutionary conservation and motif density, complex chromosomal interactions, and activity as in vivo enhancers. Our results provide insights into the combinatorial TF-pRE interactions that activate and repress expression programs during telencephalon neurogenesis and demonstrate the value of TF binding toward modeling developmental transcriptional wiring.

The drivers of immune evasion are not entirely clear, limiting the success of cancer immunotherapies. Here we applied single-cell spatial and perturbational transcriptomics to delineate immune evasion in high-grade serous tubo-ovarian cancer. To this end, we first mapped the spatial organization of high-grade serous tubo-ovarian cancer by profiling more than 2.5 million cells in situ in 130 tumors from 94 patients. This revealed a malignant cell state that reflects tumor genetics and is predictive of T cell and natural killer cell infiltration levels and response to immune checkpoint blockade. We then performed Perturb-seq screens and identified genetic perturbations-including knockout of PTPN1 and ACTR8-that trigger this malignant cell state. Finally, we show that these perturbations, as well as a PTPN1/PTPN2 inhibitor, sensitize ovarian cancer cells to T cell and natural killer cell cytotoxicity, as predicted. This study thus identifies ways to study and target immune evasion by linking genetic variation, cell-state regulators and spatial biology.

From single-cell RNA-sequencing (scRNA-seq) and spatial transcriptomics (ST), one can extract high-dimensional gene expression patterns that can be described by intercellular communication networks or decoupled gene modules. These two descriptions of information flow are often assumed to occur independently. However, intercellular communication drives directed flows of information that are mediated by intracellular gene modules, in turn triggering outflows of other signals. Methodologies to describe such intercellular flows are lacking. We present FlowSig, a method that infers communication-driven intercellular flows from scRNA-seq or ST data using graphical causal modeling and conditional independence. We benchmark FlowSig using newly generated experimental cortical organoid data and synthetic data generated from mathematical modeling. We demonstrate FlowSig's utility by applying it to various studies, showing that FlowSig can capture stimulation-induced changes to paracrine signaling in pancreatic islets, demonstrate shifts in intercellular flows due to increasing COVID-19 severity and reconstruct morphogen-driven activator-inhibitor patterns in mouse embryogenesis.

Advances in sequencing technologies have enabled the comparison of high-quality genomes of diverse primate species, revealing vast amounts of divergence due to structural variation. Given their large size, structural variants (SVs) can simultaneously alter the function and regulation of multiple genes. Studies estimate that collectively more than 3.5% of the genome is divergent in humans versus other great apes, impacting thousands of genes. Functional genomics and gene-editing tools in various model systems recently emerged as an exciting frontier - investigating the wide-ranging impacts of SVs on molecular, cellular, and systems-level phenotypes. This review examines existing research and identifies future directions to broaden our understanding of the functional roles of SVs on phenotypic innovations and diversity impacting uniquely human features, ranging from cognition to metabolic adaptations.

Diet is a robust entrainment cue that regulates diurnal rhythms of the gut microbiome. We and others have shown that disruption of the circadian clock drives the progression of colorectal cancer (CRC). While certain bacterial species have been suggested to play driver roles in CRC, it is unknown whether the intestinal clock impinges on the microbiome to accelerate CRC pathogenesis. To address this, genetic disruption of the circadian clock, in an Apc-driven mouse model of CRC, was used to define the impact on the gut microbiome. When clock disruption is combined with CRC, metagenomic sequencing identified dysregulation of many bacterial genera including Bacteroides, Helicobacter, and Megasphaera. We identify functional changes to microbial pathways including dysregulated nucleic acid, amino acid, and carbohydrate metabolism, as well as disruption of intestinal barrier function. Our findings suggest that clock disruption impinges on microbiota composition and intestinal permeability that may contribute to CRC pathogenesis.

Despite efforts to understand breast cancer biology, metastatic disease remains a clinical challenge. Identifying suppressors of breast cancer progression and mechanisms of transition to more invasive phenotypes could provide game changing therapeutic opportunities. Transcriptional deregulation is central to all malignancies, highlighted by the extensive reprogramming of regulatory elements that underlie oncogenic programs. Among these, super-enhancers (SEs) stand out due to their enrichment in genes controlling cancer hallmarks. To reveal novel breast cancer dependencies, we integrated the analysis of the SE landscape with master regulator activity inference for a series of breast cancer cell lines. As a result, we identified T - h elper-inducing Poxviruses and Zinc-finger ( PO Z)/ K rüppel-like factor (ThPOK, ZBTB7B ), a CD4 + cell lineage commitment factor, as a breast cancer master regulator that is recurrently associated with a SE. ThPOK expression is highest in luminal breast cancer but is significantly reduced in the basal subtype. Manipulation of ThPOK levels in cell lines shows that its repressive function restricts breast cancer cells to an epithelial phenotype by suppressing the expression of genes involved in the epithelial-mesenchymal transition (EMT), WNT/β-catenin target genes, and the pro-metastatic TGFβ pathway. Our study reveals ThPOK as a master transcription factor that restricts the acquisition of metastatic features in breast cancer cells.

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.

Prostate cancer is frequently treated with radiotherapy. Unfortunately, aggressive radioresistant relapses can arise, and the molecular underpinnings of radioresistance are unknown. Modern clinical radiotherapy is evolving to deliver higher doses of radiation in fewer fractions (hypofractionation). We therefore analyzed genomic, transcriptomic, and proteomic data to characterize prostate cancer radioresistance in cells treated with both conventionally fractionated and hypofractionated radiotherapy. Independent of fractionation schedule, resistance to radiotherapy involved massive genomic instability and abrogation of DNA mismatch repair. Specific prostate cancer driver genes were modulated at the RNA and protein levels, with distinct protein subcellular responses to radiotherapy. Conventional fractionation led to a far more aggressive biomolecular response than hypofractionation. Testing preclinical candidates identified in cell lines, we revealed POLQ (DNA Polymerase Theta) as a radiosensitizer. POLQ-modulated radioresistance in model systems and was predictive of it in large patient cohorts. The molecular response to radiation is highly multimodal and sheds light on prostate cancer lethality.SignificanceRadiation is standard of care in prostate cancer. Yet, we have little understanding of its failure. We demonstrate a new paradigm that radioresistance is fractionation specific and identified POLQ as a radioresistance modulator.

Understanding cell fate decision-making during complex biological processes is an open challenge that is now aided by high-resolution single-cell sequencing technologies. Specifically, it remains challenging to identify and characterize transition states corresponding to "tipping points" whereby cells commit to new cell states. Here, we present a computational method that takes advantage of single-cell transcriptomics data to infer the stability and gene regulatory networks (GRNs) along cell lineages. Our method uses the unspliced and spliced counts from single-cell RNA sequencing data and cell ordering along lineage trajectories to train an RNA splicing multivariate model, from which cell-state stability along the lineage is inferred based on spectral analysis of the model's Jacobian matrix. Moreover, the model infers the RNA cross-species interactions resulting in GRNs and their variation along the cell lineage. When applied to epithelial-mesenchymal transition in ovarian and lung cancer-derived cell lines, our model predicts a saddle-node transition between the epithelial and mesenchymal states passing through an unstable, intermediate cell state. Furthermore, we show that the underlying GRN controlling epithelial-mesenchymal transition rearranges during the transition, resulting in denser and less modular networks in the intermediate state. Overall, our method represents a flexible tool to study cell lineages with a combination of theory-driven modeling and single-cell transcriptomics data.

RNA structural switches are key regulators of gene expression in bacteria, but their characterization in Metazoa remains limited. Here, we present SwitchSeeker, a comprehensive computational and experimental approach for systematic identification of functional RNA structural switches. We applied SwitchSeeker to the human transcriptome and identified 245 putative RNA switches. To validate our approach, we characterized a previously unknown RNA switch in the 3' untranslated region of the RORC (RAR-related orphan receptor C) transcript. In vivo dimethyl sulfate (DMS) mutational profiling with sequencing (DMS-MaPseq), coupled with cryogenic electron microscopy, confirmed its existence as two alternative structural conformations. Furthermore, we used genome-scale CRISPR screens to identify trans factors that regulate gene expression through this RNA structural switch. We found that nonsense-mediated messenger RNA decay acts on this element in a conformation-specific manner. SwitchSeeker provides an unbiased, experimentally driven method for discovering RNA structural switches that shape the eukaryotic gene expression landscape.

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 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

Accessible and easy-to-use standardized bioinformatics workflows are necessary to advance microbiome research from observational studies to large-scale, data-driven approaches. Standardized multi-omics data enables comparative studies, data reuse, and applications of machine learning to model biological processes. To advance broad accessibility of standardized multi-omics bioinformatics workflows, the National Microbiome Data Collaborative (NMDC) has developed the Empowering the Development of Genomics Expertise (NMDC EDGE) resource, a user-friendly, open-source web application (https://nmdc-edge.org). Here, we describe the design and main functionality of the NMDC EDGE resource for processing metagenome, metatranscriptome, natural organic matter, and metaproteome data. The architecture relies on three main layers (web application, orchestration, and execution) to ensure flexibility and expansion to future workflows. The orchestration and execution layers leverage best practices in software containers and accommodate high-performance computing and cloud computing services. Further, we have adopted a robust user research process to collect feedback for continuous improvement of the resource. NMDC EDGE provides an accessible interface for researchers to process multi-omics microbiome data using production-quality workflows to facilitate improved data standardization and interoperability.

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.

The Human Genome Project was an enormous accomplishment, providing a foundation for countless explorations into the genetics and genomics of the human species. Yet for many years, the human genome reference sequence remained incomplete and lacked representation of human genetic diversity. Recently, two major advances have emerged to address these shortcomings: complete gap-free human genome sequences, such as the one developed by the Telomere-to-Telomere Consortium, and high-quality pangenomes, such as the one developed by the Human Pangenome Reference Consortium. Facilitated by advances in long-read DNA sequencing and genome assembly algorithms, complete human genome sequences resolve regions that have been historically difficult to sequence, including centromeres, telomeres, and segmental duplications. In parallel, pangenomes capture the extensive genetic diversity across populations worldwide. Together, these advances usher in a new era of genomics research, enhancing the accuracy of genomic analysis, paving the path for precision medicine, and contributing to deeper insights into human biology.

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.

The under-representation of non-European cohorts in neurodegenerative disease genome-wide association studies (GWAS) hampers precision medicine efforts. Despite the inherent genetic and phenotypic diversity in these diseases, GWAS research consistently exhibits a disproportionate emphasis on participants of European ancestry. This study reviews GWAS up to 2022, focusing on non-European or multi-ancestry neurodegeneration studies. We conducted a systematic review of GWAS results and publications up to 2022, focusing on non-European or multi-ancestry neurodegeneration studies. Rigorous article inclusion and quality assessment methods were employed. Of 123 neurodegenerative disease (NDD) GWAS reviewed, 82% predominantly featured European ancestry participants. A single European study identified over 90 risk loci, compared to a total of 50 novel loci in identified in all non-European or multi-ancestry studies. Notably, only six of the loci have been replicated. The significant under-representation of non-European ancestries in NDD GWAS hinders comprehensive genetic understanding. Prioritizing genomic diversity in future research is crucial for advancing NDD therapies and understanding. HIGHLIGHTS: Eighty-two percent of neurodegenerative genome-wide association studies (GWAS) focus on Europeans. Only 6 of 50 novel neurodegenerative disease (NDD) genetic loci have been replicated. Lack of diversity significantly hampers understanding of NDDs. Increasing diversity in NDD genetic research is urgently required. New initiatives are aiming to enhance diversity in NDD research.

A primary goal of psychiatry is to better understand the pathways that link genetic risk to psychiatric symptoms. Here, we tested association of diagnosis and endophenotypes with overall and neurotransmitter pathway-specific polygenic risk in patients with early-stage psychosis. Subjects included 205 demographically diverse cases with a psychotic disorder who underwent comprehensive psychiatric and neurological phenotyping and 115 matched controls. Following genotyping, we calculated polygenic scores (PGSs) for schizophrenia (SZ) and bipolar disorder (BP) using Psychiatric Genomics Consortium GWAS summary statistics. To test if overall genetic risk can be partitioned into affected neurotransmitter pathways, we calculated pathway PGSs (pPGSs) for SZ risk affecting each of four major neurotransmitter systems: glutamate, GABA, dopamine, and serotonin. Psychosis subjects had elevated SZ PGS versus controls; cases with SZ or BP diagnoses had stronger SZ or BP risk, respectively. There was no significant association within psychosis cases between individual symptom measures and overall PGS. However, neurotransmitter-specific pPGSs were moderately associated with specific endophenotypes; notably, glutamate was associated with SZ diagnosis and with deficits in cognitive control during task-based fMRI, while dopamine was associated with global functioning. Finally, unbiased endophenotype-driven clustering identified three diagnostically mixed case groups that separated on primary deficits of positive symptoms, negative symptoms, global functioning, and cognitive control. All clusters showed strong genome-wide risk. Cluster 2, characterized by deficits in cognitive control and negative symptoms, additionally showed specific risk concentrated in glutamatergic and GABAergic pathways. Due to the intensive characterization of our subjects, the present study was limited to a relatively small cohort. As such, results should be followed up with additional research at the population and mechanism level. Our study suggests pathway-based PGS analysis may be a powerful path forward to study genetic mechanisms driving psychiatric endophenotypes.

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.

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.

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.