Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “Preprint”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 235 records · Page 13Linked to original sources

Extrusion fountains are hallmarks of chromosome organization emerging upon zygotic genome activation.

The initiation of gene expression during development, known as zygotic genome activation (ZGA), is accompanied by massive changes in chromosome organization. However, the earliest events of chromosome folding and their functional roles remain unclear. Using Hi-C on zebrafish embryos, we discovered that chromosome folding begins early in development with the formation of "fountains", a novel element of chromosome organization. Emerging preferentially at enhancers, fountains exhibit an initial accumulation of cohesin, which later redistributes to CTCF sites at TAD borders. Knockouts of pioneer transcription factors driving ZGA enhancers result in the specific loss of fountains, establishing a causal link between enhancer activation and fountain formation. Polymer simulations demonstrate that fountains may arise as sites of facilitated cohesin loading, requiring two-sided but desynchronized loop extrusion, potentially caused by cohesin collisions with obstacles or internal switching. Moreover, we detected similar fountain patterns at enhancers in mouse cells. Fountains disappear upon acute cohesin depletion, as well as during mitosis, and reappear with cohesin loading in early G1. Altogether, fountains represent the first known enhancer-specific elements of chromosome organization and constitute starting points for chromosome folding during development, likely through facilitated cohesin loading.

Journal Article↗

Mega-Enhancer Bodies Organize Neuronal Long Genes in the Cerebellum.

Dynamic regulation of gene expression plays a key role in establishing the diverse neuronal cell types in the brain. Recent findings in genome biology suggest that three-dimensional (3D) genome organization has important, but mechanistically poorly understood functions in gene transcription. Beyond local genomic interactions between promoters and enhancers, we find that cerebellar granule neurons undergoing differentiation in vivo exhibit striking increases in long-distance genomic interactions between transcriptionally active genomic loci, which are separated by tens of megabases within a chromosome or located on different chromosomes. Among these interactions, we identify a nuclear subcompartment enriched for near-megabase long enhancers and their associated neuronal long genes encoding synaptic or signaling proteins. Neuronal long genes are differentially recruited to this enhancer-dense subcompartment to help shape the transcriptional identities of granule neuron subtypes in the cerebellum. SPRITE analyses of higher-order genomic interactions, together with IGM-based 3D genome modeling and imaging approaches, reveal that the enhancer-dense subcompartment forms prominent nuclear structures, which we term mega-enhancer bodies. These novel nuclear bodies reside in the nuclear periphery, away from other transcriptionally active structures, including nuclear speckles located in the nuclear interior. Together, our findings define additional layers of higher-order 3D genome organization closely linked to neuronal maturation and identity in the brain.

Journal Article↗

Oligomerization and positive feedback on membrane recruitment encode dynamically stable PAR-3 asymmetries in the C. elegans zygote.

Studies of PAR polarity have emphasized a paradigm in which mutually antagonistic PAR proteins form complementary polar domains in response to transient cues. A growing body of work suggests that the oligomeric scaffold PAR-3 can form unipolar asymmetries without mutual antagonism, but how it does so is largely unknown. Here we combine single molecule analysis and modeling to show how the interplay of two positive feedback loops promotes dynamically stable unipolar PAR-3 asymmetries in early C. elegans embryos. First, the intrinsic dynamics of PAR-3 membrane binding and oligomerization encode negative feedback on PAR-3 dissociation. Second, membrane-bound PAR-3 promotes its own recruitment through a mechanism that requires the anterior polarity proteins PAR-6 and PKC-3. Using a kinetic model tightly constrained by our experimental measurements, we show that these two feedback loops are individually required and jointly sufficient to encode dynamically stable and locally inducible unipolar PAR-3 asymmetries in the absence of posterior inhibition. Given the central role of PAR-3, and the conservation of PAR-3 membrane-binding, oligomerization, and core interactions with PAR-6/PKC-3, these results have widespread implications for PAR-mediated polarity in metazoa.

Journal Article↗

Active learning of enhancer and silencer regulatory grammar in photoreceptors.

Cis-regulatory elements (CREs) direct gene expression in health and disease, and models that can accurately predict their activities from DNA sequences are crucial for biomedicine. Deep learning represents one emerging strategy to model the regulatory grammar that relates CRE sequence to function. However, these models require training data on a scale that exceeds the number of CREs in the genome. We address this problem using active machine learning to iteratively train models on multiple rounds of synthetic DNA sequences assayed in live mammalian retinas. During each round of training the model actively selects sequence perturbations to assay, thereby efficiently generating informative training data. We iteratively trained a model that predicts the activities of sequences containing binding motifs for the photoreceptor transcription factor Cone-rod homeobox (CRX) using an order of magnitude less training data than current approaches. The model's internal confidence estimates of its predictions are reliable guides for designing sequences with high activity. The model correctly identified critical sequence differences between active and inactive sequences with nearly identical transcription factor binding sites, and revealed order and spacing preferences for combinations of motifs. Our results establish active learning as an effective method to train accurate deep learning models of cis-regulatory function after exhausting naturally occurring training examples in the genome.

Journal Article↗

Rare variation in noncoding regions with evolutionary signatures contributes to autism spectrum disorder risk.

Little is known about the role of noncoding regions in the etiology of autism spectrum disorder (ASD). We examined three classes of noncoding regions: Human Accelerated Regions (HARs), which show signatures of positive selection in humans; experimentally validated neural Vista Enhancers (VEs); and conserved regions predicted to act as neural enhancers (CNEs). Targeted and whole genome analysis of >16,600 samples and >4900 ASD probands revealed that likely recessive, rare, inherited variants in HARs, VEs, and CNEs substantially contribute to ASD risk in probands whose parents share ancestry, which enriches for recessive contributions, but modestly, if at all, in simplex family structures. We identified multiple patient variants in HARs near IL1RAPL1 and in a VE near SIM1 and showed that they change enhancer activity. Our results implicate both human-evolved and evolutionarily conserved noncoding regions in ASD risk and suggest potential mechanisms of how changes in regulatory regions can modulate social behavior.

Human Accelerated Regions↗

Natural variation of immune epitopes reveals intrabacterial antagonism.

Plants and animals detect biomolecules termed Microbe-Associated Molecular Patterns (MAMPs) and induce immunity. Agricultural production is severely impacted by pathogens which can be controlled by transferring immune receptors. However, most studies use a single MAMP epitope and the impact of diverse multi-copy MAMPs on immune induction is unknown. Here we characterized the epitope landscape from five proteinaceous MAMPs across 4,228 plant-associated bacterial genomes. Despite the diversity sampled, natural variation was constrained and experimentally testable. Immune perception in both Arabidopsis and tomato depended on both epitope sequence and copy number variation. For example, Elongation Factor Tu is predominantly single copy and 92% of its epitopes are immunogenic. Conversely, 99.9% of bacterial genomes contain multiple Cold Shock Proteins and 46% carry a non-immunogenic form. We uncovered a new mechanism for immune evasion, intrabacterial antagonism, where a non-immunogenic Cold Shock Protein blocks perception of immunogenic forms encoded in the same genome. These data will lay the foundation for immune receptor deployment and engineering based on natural variation.

comparative genomics↗

Asymmetric cortical projections to striatal direct and indirect pathways distinctly control actions.

The striatal direct and indirect pathways constitute the core for basal ganglia function in action control. Although both striatal D1- and D2-spiny projection neurons (SPNs) receive excitatory inputs from the cerebral cortex, whether or not they share inputs from the same cortical neurons, and how pathway-specific corticostriatal projections control behavior remain largely unknown. Here using a G-deleted rabies system in mice, we found that more than two-thirds of excitatory inputs to D2-SPNs also target D1-SPNs, while only one-third do so vice versa. Optogenetic stimulation of striatal D1- vs. D2-SPN-projecting cortical neurons differently regulate locomotion, reinforcement learning and sequence behavior, implying the functional dichotomy of pathway-specific corticostriatal subcircuits. These results reveal the partially segregated yet asymmetrically overlapping cortical projections on striatal D1- vs. D2-SPNs, and that the pathway-specific corticostriatal subcircuits distinctly control behavior. It has important implications in a wide range of neurological and psychiatric diseases affecting cortico-basal ganglia circuitry.

Journal Article↗

Plasticity of Human Microglia and Brain Perivascular Macrophages in Aging and Alzheimer's Disease.

The complex roles of myeloid cells, including microglia and perivascular macrophages, are central to the neurobiology of Alzheimer's disease (AD), yet they remain incompletely understood. Here, we profiled 832,505 human myeloid cells from the prefrontal cortex of 1,607 unique donors covering the human lifespan and varying degrees of AD neuropathology. We delineated 13 transcriptionally distinct myeloid subtypes organized into 6 subclasses and identified AD-associated adaptive changes in myeloid cells over aging and disease progression. The GPNMB subtype, linked to phagocytosis, increased significantly with AD burden and correlated with polygenic AD risk scores. By organizing AD-risk genes into a regulatory hierarchy, we identified and validated MITF as an upstream transcriptional activator of GPNMB, critical for maintaining phagocytosis. Through cell-to-cell interaction networks, we prioritized APOE-SORL1 and APOE-TREM2 ligand-receptor pairs, associated with AD progression. In both human and mouse models, TREM2 deficiency disrupted GPNMB expansion and reduced phagocytic function, suggesting that GPNMB's role in neuroprotection was TREM2-dependent. Our findings clarify myeloid subtypes implicated in aging and AD, advancing the mechanistic understanding of their role in AD and aiding therapeutic discovery.

Journal Article↗

Evaluating the impact of modeling choices on the performance of integrated genetic and clinical models.

The value of genetic information for improving the performance of clinical risk prediction models has yielded variable conclusions. Many methodological decisions have the potential to contribute to differential results across studies. Here, we performed multiple modeling experiments integrating clinical and demographic data from electronic health records (EHR) and genetic data to understand which decision points may affect performance. Clinical data in the form of structured diagnostic codes, medications, procedural codes, and demographics were extracted from two large independent health systems and polygenic risk scores (PRS) were generated across all patients with genetic data in the corresponding biobanks. Crohn's disease was used as the model phenotype based on its substantial genetic component, established EHR-based definition, and sufficient prevalence for model training and testing. We investigated the impact of PRS integration method, as well as choices regarding training sample, model complexity, and performance metrics. Overall, our results show that including PRS resulted in higher performance by some metrics but the gain in performance was only robust when combined with demographic data alone. Improvements were inconsistent or negligible after including additional clinical information. The impact of genetic information on performance also varied by PRS integration method, with a small improvement in some cases from combining PRS with the output of a clinical model (late-fusion) compared to its inclusion an additional feature (early-fusion). The effects of other modeling decisions varied between institutions though performance increased with more compute-intensive models such as random forest. This work highlights the importance of considering methodological decision points in interpreting the impact on prediction performance when including PRS information in clinical models.

Preprint↗

One-carbon unit supplementation fuels tumor-infiltrating T cells and augments checkpoint blockade.

Nucleotides perform important metabolic functions, carrying energy and feeding nucleic acid synthesis. Here, we use isotope tracing-mass spectrometry to quantitate the contributions to purine nucleotides of salvage versus de novo synthesis. We further explore the impact of augmenting a key precursor for purine synthesis, one-carbon (1C) units. We show that tumors and tumor-infiltrating T cells (relative to splenic T cells) synthesize purines de novo. Purine synthesis requires two 1C units, which come from serine catabolism and circulating formate. Shortage of 1C units is a potential bottleneck for anti-tumor immunity. Elevating circulating formate drives its usage by tumor-infiltrating T cells. Orally administered methanol functions as a formate pro-drug, with deuteration enabling control of formate-production kinetics. In MC38 tumors, safe doses of methanol raise formate levels and augment anti-PD-1 checkpoint blockade, tripling durable regressions. Thus, 1C deficiency can gate antitumor immunity and this metabolic checkpoint can be overcome with pharmacological 1C supplementation.

Preprint↗

Causal Associations of Sleep Apnea with Alzheimer's Disease and Cardiovascular Disease: a Bidirectional Mendelian Randomization Analysis.

BACKGROUND: Sleep apnea (SA) has been linked to an increased risk of dementia in numerous observational studies; whether this is driven by neurodegenerative, vascular or other mechanisms is not clear. We sought to examine the bidirectional causal relationships between SA, Alzheimer's disease (AD), coronary artery disease (CAD), and ischemic stroke using Mendelian randomization (MR). METHODS: Using summary statistics from four recent, large genome-wide association studies of SA (n=523,366), AD (n=64,437), CAD (n=1,165,690), and stroke (n=1,308,460), we conducted bidirectional two-sample MR analyses. Our primary analytic method was fixed-effects inverse variance weighted MR; diagnostics tests and sensitivity analyses were conducted to verify the robustness of the results. RESULTS: We identified a significant causal effect of SA on the risk of CAD (odds ratio (OR IVW ) =1.35 per log-odds increase in SA liability, 95% confidence interval (CI) =1.25-1.47) and stroke (OR IVW =1.13, 95% CI =1.01-1.25). These associations were somewhat attenuated after excluding single-nucleotide polymorphisms associated with body mass index (BMI) (OR IVW =1.26, 95% CI =1.15-1.39 for CAD risk; OR IVW =1.08, 95% CI =0.96-1.22 for stroke risk). SA was not causally associated with a higher risk of AD (OR IVW =1.14, 95% CI =0.91-1.43). We did not find causal effects of AD, CAD, or stroke on risk of SA. CONCLUSIONS: These results suggest that SA increased the risk of CAD, and the identified causal association with stroke risk may be confounded by BMI. Moreover, no causal effect of SA on AD risk was found. Future studies are warranted to investigate cardiovascular pathways between sleep disorders, including SA, and dementia.

Preprint↗

MitoTracer facilitates the identification of informative mitochondrial mutations for precise lineage reconstruction.

Mitochondrial (MT) mutations serve as natural genetic markers for inferring clonal relationships using single cell sequencing data. However, the fundamental challenge of MT mutation-based lineage tracing is automated identification of informative MT mutations. Here, we introduced an open-source computational algorithm called "MitoTracer", which accurately identified clonally informative MT mutations and inferred evolutionary lineage from scRNA-seq or scATAC-seq samples. We benchmarked MitoTracer using the ground-truth experimental lineage sequencing data and demonstrated its superior performance over the existing methods measured by high sensitivity and specificity. MitoTracer is compatible with multiple single cell sequencing platforms. Its application to a cancer evolution dataset revealed the genes related to primary BRAF-inhibitor resistance from scRNA-seq data of BRAF-mutated cancer cells. Overall, our work provided a valuable tool for capturing real informative MT mutations and tracing the lineages among cells.

Journal Article↗

Human mutations in high-confidence Tourette disorder genes affect sensorimotor behavior, reward learning, and striatal dopamine in mice.

UNLABELLED: Tourette disorder (TD) is poorly understood, despite affecting 1/160 children. A lack of animal models possessing construct, face, and predictive validity hinders progress in the field. We used CRISPR/Cas9 genome editing to generate mice with mutations orthologous to human de novo variants in two high-confidence Tourette genes, CELSR3 and WWC1 . Mice with human mutations in Celsr3 and Wwc1 exhibit cognitive and/or sensorimotor behavioral phenotypes consistent with TD. Sensorimotor gating deficits, as measured by acoustic prepulse inhibition, occur in both male and female Celsr3 TD models. Wwc1 mice show reduced prepulse inhibition only in females. Repetitive motor behaviors, common to Celsr3 mice and more pronounced in females, include vertical rearing and grooming. Sensorimotor gating deficits and rearing are attenuated by aripiprazole, a partial agonist at dopamine type II receptors. Unsupervised machine learning reveals numerous changes to spontaneous motor behavior and less predictable patterns of movement. Continuous fixed-ratio reinforcement shows Celsr3 TD mice have enhanced motor responding and reward learning. Electrically evoked striatal dopamine release, tested in one model, is greater. Brain development is otherwise grossly normal without signs of striatal interneuron loss. Altogether, mice expressing human mutations in high-confidence TD genes exhibit face and predictive validity. Reduced prepulse inhibition and repetitive motor behaviors are core behavioral phenotypes and are responsive to aripiprazole. Enhanced reward learning and motor responding occurs alongside greater evoked dopamine release. Phenotypes can also vary by sex and show stronger affection in females, an unexpected finding considering males are more frequently affected in TD. SIGNIFICANCE STATEMENT: We generated mouse models that express mutations in high-confidence genes linked to Tourette disorder (TD). These models show sensorimotor and cognitive behavioral phenotypes resembling TD-like behaviors. Sensorimotor gating deficits and repetitive motor behaviors are attenuated by drugs that act on dopamine. Reward learning and striatal dopamine is enhanced. Brain development is grossly normal, including cortical layering and patterning of major axon tracts. Further, no signs of striatal interneuron loss are detected. Interestingly, behavioral phenotypes in affected females can be more pronounced than in males, despite male sex bias in the diagnosis of TD. These novel mouse models with construct, face, and predictive validity provide a new resource to study neural substrates that cause tics and related behavioral phenotypes in TD.

Preprint↗

Experimental Test of Evolutionary Safety of a CRISPR-Cas9 Gene-Drive Element.

CRISPR-Cas9 gene drives (CCGDs) are powerful tools for genetic control of wild populations, with applications from disease eradication to species conservation. However, Cas9 alone and in a complex with gRNA can cause double-stranded DNA breaks at off-target sites, which could increase the mutational load and lead to unintended loss-of-heterozygosity (LOH) events. These undesired effects raise potential concerns about the long-term evolutionary safety of CCGDs, but the magnitude of these effects is unknown. To measure how the presence of a CCGD or a Cas9 alone in the genome affects the rates of LOH events and de novo mutations, we carried out a mutation accumulation experiment in yeast Saccharomyces cerevisiae. We found no detectable effects on the genome-wide rates of mutations or LOH events. Our power calculations suggest that CCGD or Cas9 affect these rates by less than 30%, which is much less than natural variation for these traits in yeast. A more detailed examination shows that CCGD or Cas9 may alter the lengths and genomic distributions of LOH events, but the statistical support for these effects is weak. Thus, our results demonstrate that CCGDs impose at most a weak additional mutational burden in the yeast model. Although mutagenic effects of gene drives need to be further evaluated in other systems, our results add credence to the proposition that the evolutionary risks posed by well designed gene drives are likely acceptable.

Journal Article↗

Neotelomeres and Telomere-Spanning Chromosomal Arm Fusions in Cancer Genomes Revealed by Long-Read Sequencing.

Alterations in the structure and location of telomeres are key events in cancer genome evolution. However, previous genomic approaches, unable to span long telomeric repeat arrays, could not characterize the nature of these alterations. Here, we applied both long-read and short-read genome sequencing to assess telomere repeat-containing structures in cancers and cancer cell lines. Using long-read genome sequences that span telomeric repeat arrays, we defined four types of telomere repeat variations in cancer cells: neotelomeres where telomere addition heals chromosome breaks, chromosomal arm fusions spanning telomere repeats, fusions of neotelomeres, and peri-centromeric fusions with adjoined telomere and centromere repeats. Analysis of lung adenocarcinoma genome sequences identified somatic neotelomere and telomere-spanning fusion alterations. These results provide a framework for systematic study of telomeric repeat arrays in cancer genomes, that could serve as a model for understanding the somatic evolution of other repetitive genomic elements.

Telomere↗

Pervasive correlations between causal disease effects of proximal SNPs vary with functional annotations and implicate stabilizing selection.

The genetic architecture of human diseases and complex traits has been extensively studied, but little is known about the relationship of causal disease effect sizes between proximal SNPs, which have largely been assumed to be independent. We introduce a new method, LD SNP-pair effect correlation regression (LDSPEC), to estimate the correlation of causal disease effect sizes of derived alleles between proximal SNPs, depending on their allele frequencies, LD, and functional annotations; LDSPEC produced robust estimates in simulations across various genetic architectures. We applied LDSPEC to 70 diseases and complex traits from the UK Biobank (average N=306K), meta-analyzing results across diseases/traits. We detected significantly nonzero effect correlations for proximal SNP pairs (e.g., -0.37±0.09 for low-frequency positive-LD 0-100bp SNP pairs) that decayed with distance (e.g., -0.07±0.01 for low-frequency positive-LD 1-10kb), varied with allele frequency (e.g., -0.15±0.04 for common positive-LD 0-100bp), and varied with LD between SNPs (e.g., +0.12±0.05 for common negative-LD 0-100bp) (because we consider derived alleles, positive-LD and negative-LD SNP pairs may yield very different results). We further determined that SNP pairs with shared functions had stronger effect correlations that spanned longer genomic distances, e.g., -0.37±0.08 for low-frequency positive-LD same-gene promoter SNP pairs (average genomic distance of 47kb (due to alternative splicing)) and -0.32±0.04 for low-frequency positive-LD H3K27ac 0-1kb SNP pairs. Consequently, SNP-heritability estimates were substantially smaller than estimates of the sum of causal effect size variances across all SNPs (ratio of 0.87±0.02 across diseases/traits), particularly for certain functional annotations (e.g., 0.78±0.01 for common Super enhancer SNPs)-even though these quantities are widely assumed to be equal. We recapitulated our findings via forward simulations with an evolutionary model involving stabilizing selection, implicating the action of linkage masking, whereby haplotypes containing linked SNPs with opposite effects on disease have reduced effects on fitness and escape negative selection.

Journal Article↗

Cell type-specific loops linked to RNA polymerase II elongation in human neural differentiation.

DNA is folded into higher-order structures that shape and are shaped by genome function. The role for long-range loops in the establishment of new gene expression patterns during cell fate transitions remains poorly understood. Here, we investigate the link between cell-specific loops and RNA polymerase II (RNAPolII) during neural lineage commitment. We find thousands of loops decommissioned or gained de novo upon differentiation of human induced pluripotent stem cells (hiPSCs) to neural progenitors (NPCs) and post-mitotic neurons. During hiPSC-to-NPC and NPC-to-neuron transitions, genes changing from RNAPolII initiation to elongation are >4-fold more likely to anchor cell-specific loops than repressed genes. Elongated genes exhibit significant mRNA upregulation when connected in cell-specific promoter-enhancer loops but not invariant promoter-enhancer loops, promoter-promoter loops, or unlooped. Genes transitioning from repression to RNAPolII initiation exhibit slight mRNA increase independent of loop status. Our data link cell-specific loops and robust RNAPolII-mediated elongation during neural cell fate transitions.

Preprint↗

Inferring Metabolic States from Single Cell Transcriptomic Data via Geometric Deep Learning.

The ability to measure gene expression at single-cell resolution has elevated our understanding of how biological features emerge from complex and interdependent networks at molecular, cellular, and tissue scales. As technologies have evolved that complement scRNAseq measurements with things like single-cell proteomic, epigenomic, and genomic information, it becomes increasingly apparent how much biology exists as a product of multimodal regulation. Biological processes such as transcription, translation, and post-translational or epigenetic modification impose both energetic and specific molecular demands on a cell and are therefore implicitly constrained by the metabolic state of the cell. While metabolomics is crucial for defining a holistic model of any biological process, the chemical heterogeneity of the metabolome makes it particularly difficult to measure, and technologies capable of doing this at single-cell resolution are far behind other multiomics modalities. To address these challenges, we present GEFMAP (Gene Expression-based Flux Mapping and Metabolic Pathway Prediction), a method based on geometric deep learning for predicting flux through reactions in a global metabolic network using transcriptomics data, which we ultimately apply to scRNAseq. GEFMAP leverages the natural graph structure of metabolic networks to learn both a biological objective for each cell and estimate a mass-balanced relative flux rate for each reaction in each cell using novel deep learning models.

Preprint↗