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 253 records · Page 14Linked to original sources

Transient Zn2+ deficiency induces replication stress and compromises daughter cell proliferation.

Cells must replicate their genome quickly and accurately, and they require metabolites and cofactors to do so. Ionic zinc (Zn2+) is an essential micronutrient that is required for hundreds of cellular processes, including DNA synthesis and adequate proliferation. Deficiency in this micronutrient impairs DNA synthesis and inhibits proliferation, but the mechanism is unknown. Using fluorescent reporters to track single cells via long-term live-cell imaging, we find that Zn2+ is required at the G1/S transition and during S-phase for timely completion of S-phase. A short pulse of Zn2+ deficiency impairs DNA synthesis and increases markers of replication stress. These markers of replication stress are reversed upon resupply of Zn2+. Finally, we find that if Zn2+ is removed during the mother cell's S-phase, daughter cells enter a transient quiescent state, maintained by sustained expression of p21, which disappears upon reentry into the cell cycle. In summary, short pulses of mild Zn2+ deficiency in S-phase specifically induce replication stress, which causes downstream proliferation impairments in daughter cells.

Biological sciences↗

Cyclin A/Cdk1 promotes chromosome alignment and timely mitotic progression.

To ensure genomic fidelity a series of spatially and temporally coordinated events are executed during prometaphase of mitosis, including bipolar spindle formation, chromosome attachment to spindle microtubules at kinetochores, the correction of erroneous kinetochore-microtubule (k-MT) attachments, and chromosome congression to the spindle equator. Cyclin A/Cdk1 kinase plays a key role in destabilizing k-MT attachments during prometaphase to promote correction of erroneous k-MT attachments. However, it is unknown if Cyclin A/Cdk1 kinase regulates other events during prometaphase. Here, we investigate additional roles of Cyclin A/Cdk1 in prometaphase by using an siRNA knockdown strategy to deplete endogenous Cyclin A from human cells. We find that depleting Cyclin A significantly extends mitotic duration, specifically prometaphase, because chromosome alignment is delayed. Unaligned chromosomes display erroneous monotelic, syntelic, or lateral k-MT attachments suggesting that bioriented k-MT attachment formation is delayed in the absence of Cyclin A. Mechanistically, chromosome alignment is likely impaired because the localization of the kinetochore proteins BUB1 kinase, KNL1, and MPS1 kinase are reduced in Cyclin A-depleted cells. Moreover, we find that Cyclin A promotes BUB1 kinetochore localization independently of its role in destabilizing k-MT attachments. Thus, Cyclin A/Cdk1 facilitates chromosome alignment during prometaphase to support timely mitotic progression.

Preprint↗

Genomic signatures of innovation and selection in the extremotolerant yeast Kluyveromyces marxianus.

Extremophiles can be the product of millions of years of evolutionary engineering and refinement. The underlying mechanisms can be quite distinct from the ones operating at earlier stages of trait innovation. In this work, we have developed the compost yeast Kluyveromyces marxianus, which diverged from its closest relative >20 million years ago, as a model for interspecies comparative biology and genomics. We applied a battery of growth assays to species of the Kluyveromyces genus and found that K. marxianus outperformed its relatives in a battery of heat and chemical stress conditions. We then generated and analyzed genomes from across the genus, to find derived genetic features associated with, and potentially causal for, K. marxianus traits. We found robust expansions in gene families in the K. marxianus genome, most notably among genes annotated as transmembrane transporters and in metabolism. In molecular-evolution tests, we identified adaptive protein variants at hundreds of genes, among which plasma membrane transporters were over-represented. Together, these signals enable a model for the molecular mechanisms and evolutionary pressures underlying K. marxianus traits, including gains in transporter function mediating stress resistance, and metabolic variants contributing to its capacity for rapid growth in challenging conditions. Such oligogenic architectures may be the rule rather than the exception in phenotypes that have evolved over long timescales.

Journal Article↗

Rhinovirus infection of airway epithelial cells uncovers the non-ciliated subset as a likely driver of genetic susceptibility to childhood-onset asthma.

Asthma is a complex disease caused by genetic and environmental factors. Epidemiological studies have shown that in children, wheezing during rhinovirus infection (a cause of the common cold) is associated with asthma development during childhood. This has led scientists to hypothesize there could be a causal relationship between rhinovirus infection and asthma or that RV-induced wheezing identifies individuals at increased risk for asthma development. However, not all children who wheeze when they have a cold develop asthma. Genome-wide association studies (GWAS) have identified hundreds of genetic variants contributing to asthma susceptibility, with the vast majority of likely causal variants being non-coding. Integrative analyses with transcriptomic and epigenomic datasets have indicated that T cells drive asthma risk, which has been supported by mouse studies. However, the datasets ascertained in these integrative analyses lack airway epithelial cells. Furthermore, large-scale transcriptomic T cell studies have not identified the regulatory effects of most non-coding risk variants in asthma GWAS, indicating there could be additional cell types harboring these "missing regulatory effects". Given that airway epithelial cells are the first line of defense against rhinovirus, we hypothesized they could be mediators of genetic susceptibility to asthma. Here we integrate GWAS data with transcriptomic datasets of airway epithelial cells subject to stimuli that could induce activation states relevant to asthma. We demonstrate that epithelial cultures infected with rhinovirus significantly upregulate childhood-onset asthma-associated genes. We show that this upregulation occurs specifically in non-ciliated epithelial cells. This enrichment for genes in asthma risk loci, or 'asthma heritability enrichment' is also significant for epithelial genes upregulated with influenza infection, but not with SARS-CoV-2 infection or cytokine activation. Additionally, cells from patients with asthma showed a stronger heritability enrichment compared to cells from healthy individuals. Overall, our results suggest that rhinovirus infection is an environmental factor that interacts with genetic risk factors through non-ciliated airway epithelial cells to drive childhood-onset asthma.

Preprint↗

Distinguishing different psychiatric disorders using DDx-PRS.

Despite great progress on methods for case-control polygenic prediction (e.g. schizophrenia vs. control), there remains an unmet need for a method that genetically distinguishes clinically related disorders (e.g. schizophrenia (SCZ) vs. bipolar disorder (BIP) vs. depression (MDD) vs. control); such a method could have important clinical value, especially at disorder onset when differential diagnosis can be challenging. Here, we introduce a method, Differential Diagnosis-Polygenic Risk Score (DDx-PRS), that jointly estimates posterior probabilities of each possible diagnostic category (e.g. SCZ=50%, BIP=25%, MDD=15%, control=10%) by modeling variance/covariance structure across disorders, leveraging case-control polygenic risk scores (PRS) for each disorder (computed using existing methods) and prior clinical probabilities for each diagnostic category. DDx-PRS uses only summary-level training data and does not use tuning data, facilitating implementation in clinical settings. In simulations, DDx-PRS was well-calibrated (whereas a simpler approach that analyzes each disorder marginally was poorly calibrated), and effective in distinguishing each diagnostic category vs. the rest. We then applied DDx-PRS to Psychiatric Genomics Consortium SCZ/BIP/MDD/control data, including summary-level training data from 3 case-control GWAS ( N =41,917-173,140 cases; total N =1,048,683) and held-out test data from different cohorts with equal numbers of each diagnostic category (total N =11,460). DDx-PRS was well-calibrated and well-powered relative to these training sample sizes, attaining AUCs of 0.66 for SCZ vs. rest, 0.64 for BIP vs. rest, 0.59 for MDD vs. rest, and 0.68 for control vs. rest. DDx-PRS produced comparable results to methods that leverage tuning data, confirming that DDx-PRS is an effective method. True diagnosis probabilities in top deciles of predicted diagnosis probabilities were considerably larger than prior baseline probabilities, particularly in projections to larger training sample sizes, implying considerable potential for clinical utility under certain circumstances. In conclusion, DDx-PRS is an effective method for distinguishing clinically related disorders.

Journal Article↗

Unique growth and morphology properties of Clade 5 Clostridioides difficile strains revealed by single-cell time-lapse microscopy.

Clostridioides difficile is a gastrointestinal pathogen of both humans and agricultural animals and thus a major One Health threat. The C. difficile species consists of five main clades, with Clade 5 currently undergoing speciation from Clades 1-4. Clade 5 strains are highly prevalent in agricultural animals and can cause zoonotic infections, suggesting that these strains have evolved phenotypes that distinguish them from Clade 1-4 strains. Here, we compare the growth properties of Clade 5 strains to those of Clade 1-4 strains using anaerobic time-lapse microscopy coupled with automated image analysis. Our analyses indicate that Clade 5 strains grow faster and are more likely to form long chains of cells than Clade 1-4 strains. Using comparative genomic and CRISPRi analyses, we show that the chaining phenotype of Clade 5 strains is driven by the orientation of the invertible cmr switch sequence, with chaining strains exhibiting a bias to the cmr-ON state. Interestingly, Clade 5 strains with a bias towards the cmr-ON state shifted to a largely cmr-OFF state during murine infection, suggesting that the cmr-OFF state is under positive selection during infection. Collectively, our data reveal that Clade 5 strains have distinct growth properties, which may allow them to inhabit diverse ecological niches.

Journal Article↗

Complete chromosome 21 centromere sequencing of families with Down syndrome reveals centromere size asymmetry.

Down syndrome, the most common form of human intellectual disability, is caused by nondisjunction and chromosome 21 trisomy (T21). Small centromeres have been hypothesized to contribute to its aetiology and studies on mammals suggest that larger centromeres are more efficiently transmitted, yet complete sequencing of chromosome 21 (chr21) centromeres has been particularly challenging. Using long-read sequencing, we sequenced and assembled the centromeres from eight families that include a child with free T21 (1 trio, 6 child-mother duos, and 1 singleton) all resulting from maternal meiosis I errors. Two of these families carry the smallest chr21 centromeres (143 and 181 kbp) observed in female individuals to date, exhibiting a ~10.7- and ~19.4-fold centromeric α-satellite higher-order repeat array size difference between the maternally inherited homologs, respectively. In both cases, the longer centromere harbors a poorly defined centromere dip region, marked by DNA hypomethylation, in the proband but not in the mother. A comparison of all proband chr21 centromeres (n=24) to those of controls (n=261) shows that small centromeres are not enriched in families with T21 (p-value=0.73); contrarily, chr21 extreme centromere size asymmetry (>10-fold) is unique of T21 (p-value=0.003), suggesting that this feature may represent a genetic risk factor for a subset of families with free T21. Additionally, phylogenetic reconstruction reveals that human chr21 has been particularly prone to such variation with some of the biggest size differences occurring over the last ~17 thousand years of human evolution.

Down syndrome↗

Functional overlap between the mammalian Sar1a and Sar1b paralogs in vivo.

Proteins carrying a signal peptide and/or a transmembrane domain enter the intracellular secretory pathway at the endoplasmic reticulum (ER) and are transported to the Golgi apparatus via COPII vesicles or tubules. SAR1 initiates COPII coat assembly by recruiting other coat proteins to the ER membrane. Mammalian genomes encode two SAR1 paralogs, SAR1A and SAR1B. While these paralogs exhibit ~90% amino acid sequence identity, it is unknown whether they perform distinct or overlapping functions in vivo. We now report that genetic inactivation of Sar1a in mice results in lethality during mid-embryogenesis. We also confirm previous reports that complete deficiency of murine Sar1b results in perinatal lethality. In contrast, we demonstrate that deletion of Sar1b restricted to hepatocytes is compatible with survival, though resulting in hypocholesterolemia that can be rescued by adenovirus-mediated overexpression of either SAR1A or SAR1B. To further examine the in vivo function of these 2 paralogs, we genetically engineered mice with the Sar1a coding sequence replacing that of Sar1b at the endogenous Sar1b locus. Mice homozygous for this allele survive to adulthood and are phenotypically normal, demonstrating complete or near-complete overlap in function between the two SAR1 protein paralogs in mice. These data also suggest upregulation of SAR1A gene expression as a potential approach for the treatment of SAR1B deficiency (chylomicron retention disease) in humans.

Preprint↗

Alternative pre-mRNA Splicing and Gene Expression Patterns in Midbrain Lineage Cells Carrying Familial Parkinson's Disease Mutations.

Parkinson's disease (PD) arises from genetic and environmental factors. Human genetics has identified mutations in ~20 inherited familial genes linked to monogenic forms of PD. To investigate the effects of individual familial PD mutations, human pluripotent embryonic stem cells (hPSCs) carrying 12 distinct familial PD mutations were differentiated into midbrain lineage cells, including dopaminergic (mDA) neurons. Global gene expression and pre-mRNA splicing patterns were analyzed in midbrain cultures carrying pathogenic PD mutations in the PRKN, SNCA, LRRK2, PINK1, DNAJC6, FBXO7, SYNJ1, DJ1, VPS13C, ATP13A2 and GBA1 genes. This analysis revealed that these familial PD mutations lead to pre-mRNA splicing changes linked to RNA splicing factors and to pathways controlling cell projections, cytoskeleton, GTPase regulation and others. Importantly, we have also shown that subsets of these splicing changes overlap with changes found in PD patient postmortem brains. Mutation-specific pre-mRNA isoforms may function as both diagnostic biomarkers for familial PD-associated genotypes and promising therapeutic targets.

Journal Article↗

A Massively Parallel CRISPR-Based Screening Platform for Modifiers of Neuronal Activity.

Understanding the complex interplay between gene expression and neuronal activity is crucial for unraveling the molecular mechanisms underlying cognitive function and neurological disorders. Here, we developed pooled screens for neuronal activity, using CRISPR interference (CRISPRi) and the fluorescent calcium integrator CaMPARI2. Using this screening method, we evaluated 1343 genes for their effect on excitability in human iPSC-derived neurons, revealing potential links to neurodegenerative and neurodevelopmental disorders. These genes include known regulators of neuronal excitability, such as TARPs and ion channels, as well as genes associated with autism spectrum disorder and Alzheimer's disease not previously described to affect neuronal excitability. This CRISPRi-based screening platform offers a versatile tool to uncover molecular mechanisms controlling neuronal activity in health and disease.

Journal Article↗

Nonlinear transcriptional responses to gradual modulation of transcription factor dosage.

Genomic loci associated with common traits and diseases are typically non-coding and likely impact gene expression, sometimes coinciding with rare loss-of-function variants in the target gene. However, our understanding of how gradual changes in gene dosage affect molecular, cellular, and organismal traits is currently limited. To address this gap, we induced gradual changes in gene expression of four genes using CRISPR activation and inactivation. Downstream transcriptional consequences of dosage modulation of three master trans-regulators associated with blood cell traits (GFI1B, NFE2, and MYB) were examined using targeted single-cell multimodal sequencing. We showed that guide tiling around the TSS is the most effective way to modulate cis gene expression across a wide range of fold-changes, with further effects from chromatin accessibility and histone marks that differ between the inhibition and activation systems. Our single-cell data allowed us to precisely detect subtle to large gene expression changes in dozens of trans genes, revealing that many responses to dosage changes of these three TFs are nonlinear, including non-monotonic behaviours, even when constraining the fold-changes of the master regulators to a copy number gain or loss. We found that the dosage properties are linked to gene constraint and that some of these nonlinear responses are enriched for disease and GWAS genes. Overall, our study provides a straightforward and scalable method to precisely modulate gene expression and gain insights into its downstream consequences at high resolution.

Journal Article↗

Targeted, Genome-scale Overexpression in Proteobacteria.

Targeted, genome-scale gene perturbation screens using Clustered Regularly Interspaced Short Palindromic Repeats interference (CRISPRi) and activation (CRISPRa) have revolutionized eukaryotic genetics, advancing medical, industrial, and basic research. Although CRISPRi knockdowns have been broadly applied in bacteria, options for genome-scale gene overexpression face key limitations. Here, we develop a facile approach for genome-scale overexpression in bacteria we call, "CRISPRtOE" (CRISPR transposition and OverExpression). We first create a platform for comprehensive gene targeting using CRISPR-associated transposons (CAST) and show that transposition occurs at a higher frequency in non-transcribed DNA. We then demonstrate that CRISPRtOE can upregulate gene expression in Proteobacteria with medical and industrial relevance by integrating synthetic promoters of varying strength upstream of target genes. Finally, we employ CRISPRtOE screening at the genome-scale in the model bacterium Escherichia coli and the non-model biofuel producer Zymomonas mobilis, recovering known and novel antibiotic and engineering targets. We envision that CRISPRtOE will be a valuable overexpression tool for antibiotic mode of action, industrial strain optimization, and gene function discovery in bacteria.

Journal Article↗

Proteome-scale tissue mapping using mass spectrometry based on label-free and multiplexed workflows.

Multiplexed bimolecular profiling of tissue microenvironment, or spatial omics, can provide deep insight into cellular compositions and interactions in healthy and diseased tissues. Proteome-scale tissue mapping, which aims to unbiasedly visualize all the proteins in a whole tissue section or region of interest, has attracted significant interest because it holds great potential to directly reveal diagnostic biomarkers and therapeutic targets. While many approaches are available, however, proteome mapping still exhibits significant technical challenges in both protein coverage and analytical throughput. Since many of these existing challenges are associated with mass spectrometry-based protein identification and quantification, we performed a detailed benchmarking study of three protein quantification methods for spatial proteome mapping, including label-free, TMT-MS2, and TMT-MS3. Our study indicates label-free method provided the deepest coverages of ~3500 proteins at a spatial resolution of 50 µm and the highest quantification dynamic range, while TMT-MS2 method holds great benefit in mapping throughput at >125 pixels per day. The evaluation also indicates both label-free and TMT-MS2 provide robust protein quantifications in identifying differentially abundant proteins and spatially co-variable clusters. In the study of pancreatic islet microenvironment, we demonstrated deep proteome mapping not only enables the identification of protein markers specific to different cell types, but more importantly, it also reveals unknown or hidden protein patterns by spatial co-expression analysis.

Journal Article↗

The expected polygenic risk score (ePRS) framework: an equitable metric for quantifying polygenetic risk via modeling of ancestral makeup.

Polygenic risk scores (PRSs) depend on genetic ancestry due to differences in allele frequencies between ancestral populations. This leads to implementation challenges in diverse populations. We propose a framework to calibrate PRS based on ancestral makeup. We define a metric called "expected PRS" (ePRS), the expected value of a PRS based on one's global or local admixture patterns. We further define the "residual PRS" (rPRS), measuring the deviation of the PRS from the ePRS. Simulation studies confirm that it suffices to adjust for ePRS to obtain nearly unbiased estimates of the PRS-outcome association without further adjusting for PCs. Using the TOPMed dataset, the estimated effect size of the rPRS adjusting for the ePRS is similar to the estimated effect of the PRS adjusting for genetic PCs. Similarly, we applied the ePRS framework to six cardiovascular-related traits in the All of Us dataset, and the results are consistent with those from the TOPMed analysis. The ePRS framework can protect from population stratification in association analysis and provide an equitable strategy to quantify genetic risk across diverse populations.

Journal Article↗

Structure of a polymorphic repeat at the CACNA1C schizophrenia locus.

Genetic variation within intron 3 of the CACNA1C calcium channel gene is associated with schizophrenia and other neuropsychiatric disorders, but analysis of the causal variants and their effect is complicated by a nearby variable-number tandem repeat (VNTR). Here, we explored the structure and population variability of the CACNA1C intron 3 VNTR using 155 long-read genome assemblies from 78 diverse individuals. Based on sequence differences among repeat units, we clustered individual sequences into 7 VNTR structural alleles called Types. Three Types were related through large duplications, but the other Types diverged much earlier such that only 12 repeat units at the 5' end of the VNTR were shared across most Types. The most diverged Types were rare and present only in individuals with African ancestry, but a multiallelic structural polymorphism was present across populations at different frequencies, consistent with expansion of the VNTR preceding the emergence of early hominins. We demonstrated that this polymorphism was in complete linkage disequilibrium with fine-mapped schizophrenia variants from genome-wide association studies (GWAS), and that this risk haplotype was associated with decreased CACNA1C gene expression in the brain. Our work suggests that sequence variation within a human-specific VNTR affects gene expression, and provides a detailed characterization of new alleles at a flagship neuropsychiatric locus.

Variable-number tandem repeat↗

Comprehensive analyses of a large human gut Bacteroidales culture collection reveal species and strain level diversity and evolution.

Species of the Bacteroidales order are among the most abundant and stable bacterial members of the human gut microbiome with diverse impacts on human health. While Bacteroidales strains and species are genomically and functionally diverse, order-wide comparative analyses are lacking. We cultured and sequenced the genomes of 408 Bacteroidales isolates from healthy human donors representing nine genera and 35 species and performed comparative genomic, gene-specific, mobile gene, and metabolomic analyses. Families, genera, and species could be grouped based on many distinctive features. However, we also show extensive DNA transfer between diverse families, allowing for shared traits and strain evolution. Inter- and intra-specific diversity is also apparent in the metabolomic profiling studies. This highly characterized and diverse Bacteroidales culture collection with strain-resolved genomic and metabolomic analyses can serve as a resource to facilitate informed selection of strains for microbiome reconstitution.

Preprint↗

CROPseq-multi: a universal solution for multiplexed perturbation in high-content pooled CRISPR screens.

Forward genetic screens seek to dissect complex biological systems by systematically perturbing genetic elements and observing the resulting phenotypes. While standard screening methodologies introduce individual perturbations, multiplexing perturbations improves the performance of single-target screens and enables combinatorial screens for the study of genetic interactions. Current tools for multiplexing perturbations are limited by technical challenges and do not offer compatibility across diverse screening methodologies, including enrichment, single-cell sequencing, and optical pooled screens. Here, we report the development of CROPseq-multi (CSM), a CROPseq1-inspired lentiviral system to multiplex Streptococcus pyogenes (Sp) Cas9-based perturbations with versatile readout compatibility and high performance for both perturbation and barcode identification. CSM has equivalent per-guide activity to CROPseq and low lentiviral recombination frequencies. Dual-guide CSM libraries are constructed in a single, facile molecular cloning step that facilitates the use of unique molecular identifiers. CSM is compatible with enrichment screening methodologies, single-cell RNA-sequencing readouts, and optical pooled screens. For optical pooled screens, an optimized and multiplexed in situ detection protocol improves barcode counts 10-fold (for mRNA detection), enables detection of recombination events, and reduces the number of sequencing cycles required for decoding by 3-fold relative to CROPseq. CROPseq-multi-v2 (CSMv2) adds compatibility for detection methods based on T7 RNA polymerase in vitro transcription2-5. CSM provides a single system for CRISPR screens that is compatible with individual and combinatorial perturbations, diverse SpCas9-based perturbation technologies, and multiple high-content, single-cell phenotypic readouts.

CRISPR Cas9↗

Natural history of eukaryotic DNA viruses with double jelly-roll major capsid proteins.

The phylum Preplasmiviricota (kingdom Bamfordvirae, realm Varidnaviria) is a broad assemblage of diverse viruses with comparatively short double-stranded DNA genomes (<50 kbp) that produce icosahedral capsids built from double jelly-roll major capsid proteins. Preplasmiviricots infect hosts from all cellular domains, testifying to their ancient origin and, in particular, are associated with six of the seven supergroups of eukaryotes. Preplasmiviricots comprise four major groups of viruses, namely, polintons, polinton-like viruses (PLVs), virophages, and adenovirids. We employed protein structure modeling and analysis to show that protein-primed DNA polymerases (pPolBs) of polintons, virophages, and cytoplasmic linear plasmids encompass an N-terminal domain homologous to the terminal proteins (TPs) of prokaryotic PRD1-like tectivirids and eukaryotic adenovirids that are involved in protein-primed replication initiation, followed by a viral ovarian tumor-like cysteine deubiquitinylase (vOTU) domain. The vOTU domain is likely responsible for the cleavage of the TP from the large pPolB polypeptide and is inactivated in adenovirids, in which TP is a separate protein. Many PLVs and transpovirons encode a distinct derivative of polinton-like pPolB that retains the TP, vOTU and pPolB polymerization palm domains but lacks the exonuclease domain and instead contains a supefamily 1 helicase domain. Analysis of the presence/absence and inactivation of the vOTU domains, and replacement of pPolB with other DNA polymerases in eukaryotic preplasmiviricots enabled us to outline a complete scenario for their origin and evolution.

adenovirus↗