Whole genome sequencing unravels cryptic circulation of divergent dengue virus lineages in the rainforest region of Nigeria

Dengue is often misclassified and underreported in Africa due to inaccurate differential diagnoses of nonspecific febrile illnesses such as malaria, sparsity of diagnostic testing and poor clinical and genomic surveillance. There are limited reports on the seroprevalence and genetic diversity of dengue virus (DENV) in humans and vectors in Nigeria. In this study, we investigated the epidemiology and genetic diversity of dengue in the rainforest region of Nigeria. We screened 515 febrile patients who tested negative for malaria and typhoid fever in three hospitals in Oyo and Ekiti States in southern Nigeria with a combination of anti-dengue IgG/IgM/NS1 rapid test kits and metagenomic sequencing. We found that approximately 28% of screened patients had previous DENV exposure, with the highest prevalence in persons over sixty. Approximately 8% of the patients showed evidence of recent or current infection, and 2.7% had acute infection. Following sequencing of sixty samples, we assembled twenty DENV-1 genomes (3 complete and 17 partial). We found that all assembled genomes belonged to DENV-1 genotype III. Our phylogenetic analyses showed evidence of prolonged cryptic circulation of divergent DENV lineages in Oyo state. We were unable to resolve the source of DENV in Nigeria owing to limited sequencing data from the region. However, our sequences clustered closely with sequences in Tanzania and sequences reported in Chinese with travel history to Tanzania in 2019. This may reflect the wider unsampled bidirectional transmission of DENV-1 in Africa, which strongly emphasizes the importance of genomic surveillance in monitoring ongoing DENV transmission in Africa.

Epigenome-wide association study of long-term psychosocial stress in older adults

Long-term psychosocial stress is strongly associated with negative physical and mental health outcomes, as well as adverse health behaviours; however, little is known about the role that stress plays on the epigenome. One proposed mechanism by which stress affects DNA methylation is through health behaviours. We conducted an epigenome-wide association study (EWAS) of cumulative psychosocial stress (n = 2,689) from the Health and Retirement Study (mean age = 70.4 years), assessing DNA methylation (Illumina Infinium HumanMethylationEPIC Beadchip) at 789,656 CpG sites. For identified CpG sites, we conducted a formal mediation analysis to examine whether smoking, alcohol use, physical activity, and body mass index (BMI) mediate the relationship between stress and DNA methylation. Nine CpG sites were associated with psychosocial stress (all p < 9E-07; FDR q < 0.10). Additionally, health behaviours and/or BMI mediated 9.4% to 21.8% of the relationship between stress and methylation at eight of the nine CpGs. Several of the identified CpGs were in or near genes associated with cardiometabolic traits, psychosocial disorders, inflammation, and smoking. These findings support our hypothesis that psychosocial stress is associated with DNA methylation across the epigenome. Furthermore, specific health behaviours mediate only a modest percentage of this relationship, providing evidence that other mechanisms may link stress and DNA methylation.

Genetic variant classification by predicted protein structure: A case study on IRF6

Next-generation genome sequencing has revolutionized genetic testing, identifying numerous rare disease-associated gene variants. However, to impute pathogenicity, computational approaches remain inadequate and functional testing of gene variant is required to provide the highest level of evidence. The emergence of AlphaFold2 has transformed the field of protein structure determination, and here we outline a strategy that leverages predicted protein structure to enhance genetic variant classification. We used the gene IRF6 as a case study due to its clinical relevance, its critical role in cleft lip/palate malformation, and the availability of experimental data on the pathogenicity of IRF6 gene variants through phenotype rescue experiments in irf6-/- zebrafish. We compared results from over 30 pathogenicity prediction tools on 37 IRF6 missense variants. IRF6 lacks an experimentally derived structure, so we used predicted structures to explore associations between mutational clustering and pathogenicity. We found that among these variants, 19 of 37 were unanimously predicted as deleterious by computational tools. Comparing in silico predictions with experimental findings, 12 variants predicted as pathogenic were experimentally determined as benign. Even with the recently published AlphaMissense model, 15/18 (83%) of the predicted pathogenic variants were experimentally determined as benign. In comparison, mapping variants to the protein revealed deleterious mutation clusters around the protein binding domain, whereas N-terminal variants tend to be benign, suggesting the importance of structural information in determining pathogenicity of mutations in this gene. In conclusion, incorporating gene-specific structural features of known pathogenic/benign mutations may provide meaningful insights into pathogenicity predictions in a gene-specific manner and facilitate the interpretation of variant pathogenicity.

Lassa virus in novel hosts: insights into the epidemiology of lassa virus infections in southern Nigeria

Identification of the diverse animal hosts responsible for spill-over events from animals to humans is crucial for comprehending the transmission patterns of emerging infectious diseases, which pose significant public health risks. To better characterize potential animal hosts of Lassa virus (LASV), we assessed domestic and non-domestic animals from 2021-2022 in four locations in southern Nigeria with reported cases of Lassa fever (LF). Birds, lizards, and domestic mammals (dogs, pigs, cattle and goats) were screened using RT-qPCR, and whole genome sequencing was performed for lineage identification on selected LASV positive samples. Animals were also screened for exposure to LASV by enzyme-linked immunosorbent assay (ELISA). Among these animals, lizards had the highest positivity rate by PCR. Genomic sequencing of samples in most infected animals showed sub-lineage 2 g of LASV. Seropositivity was highest among cattle and lowest in pigs. Though the specific impact these additional hosts may have in the broader virus-host context are still unknown - specifically relating to pathogen diversity, evolution, and transmission - the detection of LASV in non-rodent hosts living in proximity to confirmed human LF cases suggests their involvement during transmission as potential reservoirs. Additional epidemiological data comparing viral genomes from humans and animals, as well as those circulating within the environment will be critical in understanding LASV transmission dynamics and will ultimately guide the development of countermeasures for this zoonotic health threat.

Methylation patterns associated with C-reactive protein in racially and ethnically diverse populations

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 < 1.2E-7). The downstream pathways affected by DNA methylation included the identification of IFI16 and IRF7 CpG-gene transcript pairs which contributed to the innate immune response gene enrichment pathway along with NLRC5, NOD2, and AIM2. Gene enrichment analysis also identified the nuclear factor-kappaB transcription pathway. Using two-sample Mendelian randomization (MR) we inferred methylation at three CpG sites as causal for CRP levels using both White and Black or African American MR instrument variables. Overall, we identified novel CpG sites and gene transcripts that could be valuable in understanding the specific cellular processes and pathogenic mechanisms involved in inflammation.

Asymmetric genetic attributions for one's own prosocial versus antisocial behavior

People tend to rate prosocial or positive behavior as more strongly influenced by the actor's genes than antisocial or negative behavior. The current study tested whether people would show a similar asymmetry when rating the role of genes in their own behavior, and if so, what variables might mediate this difference. Participants were prompted to think about an example of their own behavior from the past year that was either prosocial or antisocial. Those in the prosocial condition rated the role of genetics in causing the behavior as significantly greater than did those in the antisocial condition. A mediation analysis suggested that this asymmetry could be accounted for by a tendency to view prosocial behavior as more natural and more aligned with one's true self than antisocial behavior. These findings add to a growing body of evidence suggesting that people's reasoning about genetics may be influenced by evaluative judgments.

Reflections on 'common' genetic medical history questions: Time to examine the what, why, and how

OBJECTIVE

A central goal of patient-centered care is to establish a therapeutic relationship. While remaining in tune with patient emotions, genetics providers must ask questions to understand medical histories that will inform the differential diagnosis, evaluation plan, and potential treatments.

METHODS

195 audio-recorded conversations between providers and caregivers of pediatric patients with suspected genetic conditions were coded and analyzed. Coders identified sensitive history-taking questions asked by providers related to exposures and complications during pregnancy; ancestry and consanguinity; educational attainment of the caregiver; and family structure.

RESULTS

We highlight examples of providers: using stigmatizing language about conception or consanguinity; not clarifying the intent behind questions related to caregivers' educational attainment and work history; and making presumptions or assumptions about caregivers' race and ethnicity, family structure, and exposures during pregnancy.

CONCLUSION

Some questions and phrasing considered routine by genetics providers may interfere with patient-centered care by straining attempts to establish a therapeutic, trusting relationship. Additional research is needed to assess how question asking and phrasing impact rapport building and patient experience during genetics encounters.

PRACTICE IMPLICATIONS

Review of the purpose and need for medical history questions common to genetics practice could serve to improve patient-centered care.

A comparative study of structural variant calling in WGS from Alzheimer's disease families

Detecting structural variants (SVs) in whole-genome sequencing poses significant challenges. We present a protocol for variant calling, merging, genotyping, sensitivity analysis, and laboratory validation for generating a high-quality SV call set in whole-genome sequencing from the Alzheimer's Disease Sequencing Project comprising 578 individuals from 111 families. Employing two complementary pipelines, Scalpel and Parliament, for SV/indel calling, we assessed sensitivity through sample replicates (N = 9) with in silico variant spike-ins. We developed a novel metric, D-score, to evaluate caller specificity for deletions. The accuracy of deletions was evaluated by Sanger sequencing. We generated a high-quality call set of 152,301 deletions of diverse sizes. Sanger sequencing validated 114 of 146 detected deletions (78.1%). Scalpel excelled in accuracy for deletions ≤100 bp, whereas Parliament was optimal for deletions >900 bp. Overall, 83.0% and 72.5% of calls by Scalpel and Parliament were validated, respectively, including all 11 deletions called by both Parliament and Scalpel between 101 and 900 bp. Our flexible protocol successfully generated a high-quality deletion call set and a truth set of Sanger sequencing-validated deletions with precise breakpoints spanning 1-17,000 bp.

Altered myocardial lipid regulation in junctophilin-2-associated familial cardiomyopathies

Myocardial lipid metabolism is critical to normal heart function, whereas altered lipid regulation has been linked to cardiac diseases including cardiomyopathies. Genetic variants in the JPH2 gene can cause hypertrophic cardiomyopathy (HCM) and, in some cases, dilated cardiomyopathy (DCM). In this study, we tested the hypothesis that JPH2 variants identified in patients with HCM and DCM, respectively, cause distinct alterations in myocardial lipid profiles. Echocardiography revealed clinically significant cardiac dysfunction in both knock-in mouse models of cardiomyopathy. Unbiased myocardial lipidomic analysis demonstrated significantly reduced levels of total unsaturated fatty acids, ceramides, and various phospholipids in both mice with HCM and DCM, suggesting a common metabolic alteration in both models. On the contrary, significantly increased di- and triglycerides, and decreased co-enzyme were only found in mice with HCM. Moreover, mice with DCM uniquely exhibited elevated levels of cholesterol ester. Further in-depth analysis revealed significantly altered metabolites from all the lipid classes with either similar or opposing trends in JPH2 mutant mice with HCM or DCM. Together, these studies revealed, for the first time, unique alterations in the cardiac lipid composition-including distinct increases in neutral lipids and decreases in polar membrane lipids-in mice with HCM and DCM were caused by distinct JPH2 variants. These studies may aid the development of novel biomarkers or therapeutics for these inherited disorders.

Data from electronic healthcare records expand our understanding of X-linked genetic diseases

Disease specific cohort studies have reported details on X linked (XL) disorders affecting females. We investigated the spectrum and penetrance of XL disorders seen in electronic health records (EHR). We generated a cohort of individuals diagnosed with XL disorders at Vanderbilt University Medical Center over 20 years. Our cohort included 477 males and 203 females diagnosed with 108 different XL genetic disorders. We found large differences between the female/male (F/M) ratios for various XL disorders regardless of their OMIM annotated mode of inheritance. We identified four XL recessive disorders affecting women previously only described in men. Biomarkers for XL disease had unique gender-specific patterns differing between modes of inheritance. EHRs provide large cohorts of XL genetic disorders that give new insights compared to the literature. Differences in the F/M ratios and biomarkers of XL disorders observed likely result from disease specific and sex dependent penetrance. We conclude that observed gender ratios associated with specific XL disorders may be more useful than those predicted by Mendelian genetics provided by OMIM. Our findings of a gender specific penetrance and severity for XL disorders show unexpected differences from Mendelian predictions. Further work is required to validate our findings in larger combined EHR cohorts.

Novel molecular mechanism in Malan syndrome uncovered through genome sequencing reanalysis, exon-level Array, and RNA sequencing

The NFIX gene encodes a DNA-binding protein belonging to the nuclear factor one (NFI) family of transcription factors. Pathogenic variants of NFIX are associated with two autosomal dominant Mendelian disorders, Malan syndrome (MIM 614753) and Marshall-Smith syndrome (MIM 602535), which are clinically distinct due to different disease-causing mechanisms. NFIX variants associated with Malan syndrome are missense variants mostly located in exon 2 encoding the N-terminal DNA binding and dimerization domain or are protein-truncating variants that trigger nonsense-mediated mRNA decay (NMD) resulting in NFIX haploinsufficiency. NFIX variants associated with Marshall-Smith syndrome are protein-truncating and are clustered between exons 6 and 10, including a recurrent Alu-mediated deletion of exons 6 and 7, which can escape NMD. The more severe phenotype of Marshall-Smith syndrome is likely due to a dominant-negative effect of these protein-truncating variants that escape NMD. Here, we report a child with clinical features of Malan syndrome who has a de novo NFIX intragenic duplication. Using genome sequencing, exon-level microarray analysis, and RNA sequencing, we show that this duplication encompasses exons 6 and 7 and leads to NFIX haploinsufficiency. To our knowledge, this is the first reported case of Malan Syndrome caused by an intragenic NFIX duplication.

High-Throughput Single-Cell, Single-Mitochondrial DNA Assay Using Hydrogel Droplet Microfluidics

There is growing interest in understanding the biological implications of single cell heterogeneity and heteroplasmy of mitochondrial DNA (mtDNA), but current methodologies for single-cell mtDNA analysis limit the scale of analysis to small cell populations. Although droplet microfluidics have increased the throughput of single-cell genomic, RNA, and protein analysis, their application to sub-cellular organelle analysis has remained a largely unsolved challenge. Here, we introduce an agarose-based droplet microfluidic approach for single-cell, single-mtDNA analysis, which allows simultaneous processing of hundreds of individual mtDNA molecules within >10,000 individual cells. Our microfluidic chip encapsulates individual cells in agarose beads, designed to have a sufficiently dense hydrogel network to retain mtDNA after lysis and provide a robust scaffold for subsequent multi-step processing and analysis. To mitigate the impact of the high viscosity of agarose required for mtDNA retention on the throughput of microfluidics, we developed a parallelized device, successfully achieving ~95 % mtDNA retention from single cells within our microbeads at >700,000 drops/minute. To demonstrate utility, we analyzed specific regions of the single-mtDNA using a multiplexed rolling circle amplification (RCA) assay. We demonstrated compatibility with both microscopy, for digital counting of individual RCA products, and flow cytometry for higher throughput analysis.

Time-to-Event Genome-Wide Association Study for Incident Cardiovascular Disease in People With Type 2 Diabetes

OBJECTIVE

To identify genetic risk factors for incident cardiovascular disease (CVD) among people with type 2 diabetes (T2D).

RESEARCH DESIGN AND METHODS

We conducted a multiancestry time-to-event genome-wide association study for incident CVD among people with T2D. We also tested 204 known coronary artery disease (CAD) variants for association with incident CVD.

RESULTS

Among 49,230 participants with T2D, 8,956 had incident CVD events (event rate 18.2%). We identified three novel genetic loci for incident CVD: rs147138607 (near CACNA1E/ZNF648, hazard ratio [HR] 1.23, P = 3.6 × 10-9), rs11444867 (near HS3ST1, HR 1.89, P = 9.9 × 10-9), and rs335407 (near TFB1M/NOX3, HR 1.25, P = 1.5 × 10-8). Among 204 known CAD loci, 5 were associated with incident CVD in T2D (multiple comparison-adjusted P < 0.00024, 0.05/204). A standardized polygenic score of these 204 variants was associated with incident CVD with HR 1.14 (P = 1.0 × 10-16).

CONCLUSIONS

The data point to novel and known genomic regions associated with incident CVD among individuals with T2D.

Iterative oxidation by TET1 is required for reprogramming of imprinting control regions and patterning of mouse sperm hypomethylated regions

Ten-eleven translocation (TET) enzymes iteratively oxidize 5-methylcytosine (5mC) to generate 5-hydroxymethylcytosine (5hmC), 5-formylcytosine, and 5-carboxylcytosine to facilitate active genome demethylation. Whether these bases are required to promote replication-coupled dilution or activate base excision repair during mammalian germline reprogramming remains unresolved due to the inability to decouple TET activities. Here, we generated two mouse lines expressing catalytically inactive TET1 (Tet1-HxD) and TET1 that stalls oxidation at 5hmC (Tet1-V). Tet1 knockout and catalytic mutant primordial germ cells (PGCs) fail to erase methylation at select imprinting control regions and promoters of meiosis-associated genes, validating the requirement for the iterative oxidation of 5mC for complete germline reprogramming. TET1V and TET1HxD rescue most hypermethylation of Tet1-/- sperm, suggesting the role of TET1 beyond its oxidative capability. We additionally identify a broader class of hypermethylated regions in Tet1 mutant mouse sperm that depend on TET oxidation for reprogramming. Our study demonstrates the link between TET1-mediated germline reprogramming and sperm methylome patterning.

Minimum information and guidelines for reporting a multiplexed assay of variant effect

Multiplexed assays of variant effect (MAVEs) have emerged as a powerful approach for interrogating thousands of genetic variants in a single experiment. The flexibility and widespread adoption of these techniques across diverse disciplines have led to a heterogeneous mix of data formats and descriptions, which complicates the downstream use of the resulting datasets. To address these issues and promote reproducibility and reuse of MAVE data, we define a set of minimum information standards for MAVE data and metadata and outline a controlled vocabulary aligned with established biomedical ontologies for describing these experimental designs.

Genetically guided precision medicine clinical decision support tools: a systematic review

OBJECTIVES

Patient care using genetics presents complex challenges. Clinical decision support (CDS) tools are a potential solution because they provide patient-specific risk assessments and/or recommendations at the point of care. This systematic review evaluated the literature on CDS systems which have been implemented to support genetically guided precision medicine (GPM).

MATERIALS AND METHODS

A comprehensive search was conducted in MEDLINE and Embase, encompassing January 1, 2011-March 14, 2023. The review included primary English peer-reviewed research articles studying humans, focused on the use of computers to guide clinical decision-making and delivering genetically guided, patient-specific assessments, and/or recommendations to healthcare providers and/or patients.

RESULTS

The search yielded 3832 unique articles. After screening, 41 articles were identified that met the inclusion criteria. Alerts and reminders were the most common form of CDS used. About 27 systems were integrated with the electronic health record; 2 of those used standards-based approaches for genomic data transfer. Three studies used a framework to analyze the implementation strategy.

DISCUSSION

Findings include limited use of standards-based approaches for genomic data transfer, system evaluations that do not employ formal frameworks, and inconsistencies in the methodologies used to assess genetic CDS systems and their impact on patient outcomes.

CONCLUSION

We recommend that future research on CDS system implementation for genetically GPM should focus on implementing more CDS systems, utilization of standards-based approaches, user-centered design, exploration of alternative forms of CDS interventions, and use of formal frameworks to systematically evaluate genetic CDS systems and their effects on patient care.

Measuring, visualizing, and diagnosing reference bias with biastools

Many bioinformatics methods seek to reduce reference bias, but no methods exist to comprehensively measure it. Biastools analyzes and categorizes instances of reference bias. It works in various scenarios: when the donor's variants are known and reads are simulated; when donor variants are known and reads are real; and when variants are unknown and reads are real. Using biastools, we observe that more inclusive graph genomes result in fewer biased sites. We find that end-to-end alignment reduces bias at indels relative to local aligners. Finally, we use biastools to characterize how T2T references improve large-scale bias.

Clinical case report of intractable paroxysmal sympathetic hyperactivity in TANGO2 deficiency disorder

TANGO2 deficiency disorder (TDD) is a neurodegenerative disease characterized by a broad and variable spectrum of clinical manifestations, even among individuals sharing the same pathogenic variants. Here, we report a severely affected individual with TDD presenting with intractable paroxysmal sympathetic hyperactivity (PSH). While progressive brain atrophy has been observed in TDD, PSH has not been reported. Despite comprehensive workup for an acute trigger, no definite cause was identified, and pharmacological interventions were ineffective to treat PSH. Ultimately care was redirected to comfort measures. This article expands the clinical phenotype of patients with TDD, highlights the possibility of PSH in these patients, and the need for continued research for better treatments of TDD.

Massively parallel screen uncovers many rare 3' UTR variants regulating mRNA abundance of cancer driver genes

Understanding the function of rare non-coding variants represents a significant challenge. Using MapUTR, a screening method, we studied the function of rare 3' UTR variants affecting mRNA abundance post-transcriptionally. Among 17,301 rare gnomAD variants, an average of 24.5% were functional, with 70% in cancer-related genes, many in critical cancer pathways. This observation motivated an interrogation of 11,929 somatic mutations, uncovering 3928 (33%) functional mutations in 155 cancer driver genes. Functional MapUTR variants were enriched in microRNA- or protein-binding sites and may underlie outlier gene expression in tumors. Further, we introduce untranslated tumor mutational burden (uTMB), a metric reflecting the amount of somatic functional MapUTR variants of a tumor and show its potential in predicting patient survival. Through prime editing, we characterized three variants in cancer-relevant genes (MFN2, FOSL2, and IRAK1), demonstrating their cancer-driving potential. Our study elucidates the function of tens of thousands of non-coding variants, nominates non-coding cancer driver mutations, and demonstrates their potential contributions to cancer.

Tracing genetic diversity captures the molecular basis of misfolding disease

Genetic variation in human populations can result in the misfolding and aggregation of proteins, giving rise to systemic and neurodegenerative diseases that require management by proteostasis. Here, we define the role of GRP94, the endoplasmic reticulum Hsp90 chaperone paralog, in managing alpha-1-antitrypsin deficiency on a residue-by-residue basis using Gaussian process regression-based machine learning to profile the spatial covariance relationships that dictate protein folding arising from sequence variants in the population. Covariance analysis suggests a role for the ATPase activity of GRP94 in controlling the N- to C-terminal cooperative folding of alpha-1-antitrypsin responsible for the correction of liver aggregation and lung-disease phenotypes of alpha-1-antitrypsin deficiency. Gaussian process-based spatial covariance profiling provides a standard model built on covariant principles to evaluate the role of proteostasis components in guiding information flow from genome to proteome in response to genetic variation, potentially allowing us to intervene in the onset and progression of complex multi-system human diseases.

Mitotic chromosomes are self-entangled and disentangle through a topoisomerase-II-dependent two-stage exit from mitosis

The topological state of chromosomes determines their mechanical properties, dynamics, and function. Recent work indicated that interphase chromosomes are largely free of entanglements. Here, we use Hi-C, polymer simulations, and multi-contact 3C and find that, by contrast, mitotic chromosomes are self-entangled. We explore how a mitotic self-entangled state is converted into an unentangled interphase state during mitotic exit. Most mitotic entanglements are removed during anaphase/telophase, with remaining ones removed during early G1, in a topoisomerase-II-dependent process. Polymer models suggest a two-stage disentanglement pathway: first, decondensation of mitotic chromosomes with remaining condensin loops produces entropic forces that bias topoisomerase II activity toward decatenation. At the second stage, the loops are released, and the formation of new entanglements is prevented by lower topoisomerase II activity, allowing the establishment of unentangled and territorial G1 chromosomes. When mitotic entanglements are not removed in experiments and models, a normal interphase state cannot be acquired.

The Brain Gene Registry: a data snapshot

Monogenic disorders account for a large proportion of population-attributable risk for neurodevelopmental disabilities. However, the data necessary to infer a causal relationship between a given genetic variant and a particular neurodevelopmental disorder is often lacking. Recognizing this scientific roadblock, 13 Intellectual and Developmental Disabilities Research Centers (IDDRCs) formed a consortium to create the Brain Gene Registry (BGR), a repository pairing clinical genetic data with phenotypic data from participants with variants in putative brain genes. Phenotypic profiles are assembled from the electronic health record (EHR) and a battery of remotely administered standardized assessments collectively referred to as the Rapid Neurobehavioral Assessment Protocol (RNAP), which include cognitive, neurologic, and neuropsychiatric assessments, as well as assessments for attention deficit hyperactivity disorder (ADHD) and autism spectrum disorder (ASD). Co-enrollment of BGR participants in the Clinical Genome Resource's (ClinGen's) GenomeConnect enables display of variant information in ClinVar. The BGR currently contains data on 479 participants who are 55% male, 6% Asian, 6% Black or African American, 76% white, and 12% Hispanic/Latine. Over 200 genes are represented in the BGR, with 12 or more participants harboring variants in each of these genes: CACNA1A, DNMT3A, SLC6A1, SETD5, and MYT1L. More than 30% of variants are de novo and 43% are classified as variants of uncertain significance (VUSs). Mean standard scores on cognitive or developmental screens are below average for the BGR cohort. EHR data reveal developmental delay as the earliest and most common diagnosis in this sample, followed by speech and language disorders, ASD, and ADHD. BGR data has already been used to accelerate gene-disease validity curation of 36 genes evaluated by ClinGen's BGR Intellectual Disability (ID)-Autism (ASD) Gene Curation Expert Panel. In summary, the BGR is a resource for use by stakeholders interested in advancing translational research for brain genes and continues to recruit participants with clinically reported variants to establish a rich and well-characterized national resource to promote research on neurodevelopmental disorders.

Linking genetic counseling communication skills to patient outcomes and experiences using a community-engagement and provider-engagement approach: research protocol for the GC-PRO mixed methods sequential explanatory study

INTRODUCTION

In over 50 years since the genetic counseling (GC) profession began, a systematic study of GC communication skills and patient-reported outcomes in actual sessions across multiple clinical specialties has never been conducted. To optimize GC quality and improve efficiency of care, the field must first be able to comprehensively measure GC skills and determine which skills are most critical to achieving positive patient experiences and outcomes. This study aims to characterise GC communication skills using a novel and pragmatic measure and link variations in communication skills to patient-reported outcomes, across clinical specialties and with patients from diverse backgrounds in the USA. Our community-engagement and provider-engagement approach is crucial to develop recommendations for quality, culturally informed GC care, which are greatly needed to improve GC practice.

METHODS AND ANALYSIS

A mixed methods, sequential explanatory design will be used to collect and analyze: audio-recorded GC sessions in cancer, cardiac, and prenatal/reproductive genetic indications; pre-visit and post-visit quantitative surveys capturing patient experiences and outcomes and post-visit qualitative interview data. A novel, practical checklist will measure GC communication skills. Coincidence analysis will identify patterns of GC skills that are consistent with high scores on patient-reported measures. Two-level, multilevel models will be used to evaluate how GC communication skills and other session/patient characteristics predict patient-reported outcomes. Four community advisory boards (CABs) and a genetic counselor advisory board will inform the study design and analysis.

ETHICS AND DISSEMINATION

This study has been approved by the single Institutional Review Board of the University of Minnesota. This research poses no greater than minimal risk to participants. Results from this study will be shared through national and international conferences and through community-based dissemination as guided by the study's CABs. A lay summary will also be disseminated to all participants.

Characterizing glucokinase variant mechanisms using a multiplexed abundance assay

BACKGROUND

Amino acid substitutions can perturb protein activity in multiple ways. Understanding their mechanistic basis may pinpoint how residues contribute to protein function. Here, we characterize the mechanisms underlying variant effects in human glucokinase (GCK) variants, building on our previous comprehensive study on GCK variant activity.

RESULTS

Using a yeast growth-based assay, we score the abundance of 95% of GCK missense and nonsense variants. When combining the abundance scores with our previously determined activity scores, we find that 43% of hypoactive variants also decrease cellular protein abundance. The low-abundance variants are enriched in the large domain, while residues in the small domain are tolerant to mutations with respect to abundance. Instead, many variants in the small domain perturb GCK conformational dynamics which are essential for appropriate activity.

CONCLUSIONS

In this study, we identify residues important for GCK metabolic stability and conformational dynamics. These residues could be targeted to modulate GCK activity, and thereby affect glucose homeostasis.

Effectiveness of Clopidogrel vs Alternative P2Y12 Inhibitors Based on the ABCD-GENE Score

BACKGROUND

An ABCD-GENE (age, body mass index, chronic kidney disease, diabetes, and CYP2C19 genetic variants) score ≥10 predicts reduced clopidogrel effectiveness, but its association with response to alternative therapy remains unclear.

OBJECTIVES

The aim of this study was to evaluate the association between ABCD-GENE score and the effectiveness of clopidogrel vs alternative P2Y12 inhibitor (prasugrel or ticagrelor) therapy after percutaneous coronary intervention (PCI).

METHODS

A total of 4,335 patients who underwent PCI, CYP2C19 genotyping, and P2Y12 inhibitor treatment were included. The primary outcome was major atherothrombotic events (MAE) within 1 year after PCI. Cox regression was performed to assess event risk in clopidogrel-treated (reference) vs alternatively treated patients, with stabilized inverse probability weights derived from exposure propensity scores after stratifying by ABCD-GENE score and further by CYP2C19 loss-of-function (LOF) genotype.

RESULTS

Among patients with scores <10 (n = 3,200), MAE was not different with alternative therapy vs clopidogrel (weighted HR: 0.89; 95% CI: 0.65-1.22; P = 0.475). The risk for MAE also did not significantly differ by treatment among patients with scores ≥10 (n = 1,135; weighted HR: 0.75; 95% CI: 0.51-1.11; P = 0.155). Among CYP2C19 LOF allele carriers, MAE risk appeared lower with alternative therapy in both the group with scores <10 (weighted HR: 0.50; 95% CI: 0.25-1.01; P = 0.052) and the group with scores ≥10 (weighted HR: 0.48; 95% CI: 0.29-0.80; P = 0.004), while there was no difference in the group with scores <10 and no LOF alleles (weighted HR: 1.03; 95% CI: 0.70-1.51; P = 0.885).

CONCLUSIONS

These data support the use of alternative therapy over clopidogrel in CYP2C19 LOF allele carriers after PCI, regardless of ABCD-GENE score, while clopidogrel is as effective as alternative therapy in non-LOF patients with scores <10.

DeLIVR: a deep learning approach to IV regression for testing nonlinear causal effects in transcriptome-wide association studies

Transcriptome-wide association studies (TWAS) have been increasingly applied to identify (putative) causal genes for complex traits and diseases. TWAS can be regarded as a two-sample two-stage least squares method for instrumental variable (IV) regression for causal inference. The standard TWAS (called TWAS-L) only considers a linear relationship between a gene's expression and a trait in stage 2, which may lose statistical power when not true. Recently, an extension of TWAS (called TWAS-LQ) considers both the linear and quadratic effects of a gene on a trait, which however is not flexible enough due to its parametric nature and may be low powered for nonquadratic nonlinear effects. On the other hand, a deep learning (DL) approach, called DeepIV, has been proposed to nonparametrically model a nonlinear effect in IV regression. However, it is both slow and unstable due to the ill-posed inverse problem of solving an integral equation with Monte Carlo approximations. Furthermore, in the original DeepIV approach, statistical inference, that is, hypothesis testing, was not studied. Here, we propose a novel DL approach, called DeLIVR, to overcome the major drawbacks of DeepIV, by estimating a related but different target function and including a hypothesis testing framework. We show through simulations that DeLIVR was both faster and more stable than DeepIV. We applied both parametric and DL approaches to the GTEx and UK Biobank data, showcasing that DeLIVR detected additional 8 and 7 genes nonlinearly associated with high-density lipoprotein (HDL) cholesterol and low-density lipoprotein (LDL) cholesterol, respectively, all of which would be missed by TWAS-L, TWAS-LQ, and DeepIV; these genes include BUD13 associated with HDL, SLC44A2 and GMIP with LDL, all supported by previous studies.

An ensemble penalized regression method for multi-ancestry polygenic risk prediction

Great efforts are being made to develop advanced polygenic risk scores (PRS) to improve the prediction of complex traits and diseases. However, most existing PRS are primarily trained on European ancestry populations, limiting their transferability to non-European populations. In this article, we propose a novel method for generating multi-ancestry Polygenic Risk scOres based on enSemble of PEnalized Regression models (PROSPER). PROSPER integrates genome-wide association studies (GWAS) summary statistics from diverse populations to develop ancestry-specific PRS with improved predictive power for minority populations. The method uses a combination ofL 1 (lasso) andL 2 (ridge) penalty functions, a parsimonious specification of the penalty parameters across populations, and an ensemble step to combine PRS generated across different penalty parameters. We evaluate the performance of PROSPER and other existing methods on large-scale simulated and real datasets, including those from 23andMe Inc., the Global Lipids Genetics Consortium, and All of Us. Results show that PROSPER can substantially improve multi-ancestry polygenic prediction compared to alternative methods across a wide variety of genetic architectures. In real data analyses, for example, PROSPER increased out-of-sample prediction R2 for continuous traits by an average of 70% compared to a state-of-the-art Bayesian method (PRS-CSx) in the African ancestry population. Further, PROSPER is computationally highly scalable for the analysis of large SNP contents and many diverse populations.

Stroke occurrence by hypertension treatment status in Ghana and Nigeria: A case-control study

BACKGROUND

Hypertension is preeminent among the vascular risk factors for stroke occurrence. The wide gaps in awareness, detection, treatment, and control rates of hypertension are fueling an epidemic of stroke in sub-Saharan Africa.

PURPOSE

To quantify the contribution of untreated, treated but uncontrolled, and controlled hypertension to stroke occurrence in Ghana and Nigeria.

METHODS

The Stroke Investigative Research and Educational Network (SIREN) is a case-control study across 16 study sites in Ghana and Nigeria. Cases were acute stroke (n = 3684) with age- and sex-matched stroke-free controls (n = 3684). We evaluated the associations of untreated hypertension, treated but uncontrolled hypertension, and controlled hypertension at BP of <140/90 mmHg with risk of stroke occurrence. We assessed the adjusted odds ratio and population-attributable risk of hypertension treatment control status associated with stroke occurrence.

RESULTS

The frequencies of no hypertension, untreated hypertension, treated but uncontrolled hypertension and controlled hypertension among stroke cases were 4.0%, 47.7%, 37.1%, and 9.2% vs 40.7%, 34.9%, 15.9%, and 7.7% respectively among stroke-free controls, p < 0.0001. The aOR and PAR (95% CI) for untreated hypertension were 6.58 (5.15-8.41) and 35.4% (33.4-37.4); treated but uncontrolled hypertension was 9.95 (7.60-13.02) and 35.9% (34.2-37.5); and controlled hypertension 5.37 (3.90-7.41) and 8.5% (7.6-9.5) respectively. Untreated hypertension contributed a PAR of 47.5% to the occurrence of intracerebral hemorrhage vs 29.5% for ischemic stroke. The aOR of untreated hypertension for stroke occurrence was 13.31 (7.64-23.19) for <50 years; 7.14 (4.51-11.31) for 50-64 years; and 3.48 (2.28-5.30) for 65 years or more.

CONCLUSION

The contribution of untreated hypertension and treated but uncontrolled hypertension to stroke occurrence among indigenous Africans is substantial. Implementing targeted interventions that address gaps in hypertension prevention and treatment, involving the local population, healthcare providers, and policymakers, can potentially substantially reduce the escalating burden of strokes in Africa.

Machine learning-driven identification of air toxic combinations associated with asthma symptoms among elementary school children in Spokane, Washington, USA

Air toxics are atmospheric pollutants with hazardous effects on health and the environment. Although methodological constraints have limited the number of air toxics assessed for associations with health and disease, advances in machine learning (ML) enable the assessment of a much larger set of environmental exposures. We used ML methods to conduct a retrospective study to identify combinations of 109 air toxics associated with asthma symptoms among 269 elementary school students in Spokane, Washington. Data on the frequency of asthma symptoms for these children were obtained from Spokane Public Schools. Their exposure to air toxics was estimated by using the Environmental Protection Agency's Air Toxics Screening Assessment and National Air Toxics Assessment. We defined three exposure periods: the most recent year (2019), the last three years (2017-2019), and the last five years (2014-2019). We analyzed the data using the ML-based Data-driven ExposurE Profile (DEEP) extraction method. DEEP identified 25 air toxic combinations associated with asthma symptoms in at least one exposure period. Three combinations (1,1,1-trichloroethane, 2-nitropropane, and 2,4,6-trichlorophenol) were significantly associated with asthma symptoms in all three exposure periods. Four air toxics (1,1,1-trichloroethane, 1,1,2,2-tetrachloroethane, BIS (2-ethylhexyl) phthalate (DEHP), and 2,4-dinitrophenol) were associated only in combination with other toxics, and would not have been identified by traditional statistical methods. The application of DEEP also identified a vulnerable subpopulation of children who were exposed to 13 of the 25 significant combinations in at least one exposure period. On average, these children experienced the largest number of asthma symptoms in our sample. By providing evidence on air toxic combinations associated with childhood asthma, our findings may contribute to the regulation of these toxics to improve children's respiratory health.

Estimating the overall fraction of phenotypic variance attributed to high-dimensional predictors measured with error

In prospective genomic studies (e.g., DNA methylation, metagenomics, and transcriptomics), it is crucial to estimate the overall fraction of phenotypic variance (OFPV) attributed to the high-dimensional genomic variables, a concept similar to heritability analyses in genome-wide association studies (GWAS). Unlike genetic variants in GWAS, these genomic variables are typically measured with error due to technical limitation and temporal instability. While the existing methods developed for GWAS can be used, ignoring measurement error may severely underestimate OFPV and mislead the design of future studies. Assuming that measurement error variances are distributed similarly between causal and noncausal variables, we show that the asymptotic attenuation factor equals to the average intraclass correlation coefficients of all genomic variables, which can be estimated based on a pilot study with repeated measurements. We illustrate the method by estimating the contribution of microbiome taxa to body mass index and multiple allergy traits in the American Gut Project. Finally, we show that measurement error does not cause meaningful bias when estimating the correlation of effect sizes for two traits.

Fine-tuning Large Language Models for Rare Disease Concept Normalization

Objective

We aim to develop a novel method for rare disease concept normalization by fine-tuning Llama 2, an open-source large language model (LLM), using a domain-specific corpus sourced from the Human Phenotype Ontology (HPO).

Methods

We developed an in-house template-based script to generate two corpora for fine-tuning. The first (NAME) contains standardized HPO names, sourced from the HPO vocabularies, along with their corresponding identifiers. The second (NAME+SYN) includes HPO names and half of the concept's synonyms as well as identifiers. Subsequently, we fine-tuned Llama2 (Llama2-7B) for each sentence set and conducted an evaluation using a range of sentence prompts and various phenotype terms.

Results

When the phenotype terms for normalization were included in the fine-tuning corpora, both models demonstrated nearly perfect performance, averaging over 99% accuracy. In comparison, ChatGPT-3.5 has only ~20% accuracy in identifying HPO IDs for phenotype terms. When single-character typos were introduced in the phenotype terms, the accuracy of NAME and NAME+SYN is 10.2% and 36.1%, respectively, but increases to 61.8% (NAME+SYN) with additional typo-specific fine-tuning. For terms sourced from HPO vocabularies as unseen synonyms, the NAME model achieved 11.2% accuracy, while the NAME+SYN model achieved 92.7% accuracy.

Conclusion

Our fine-tuned models demonstrate ability to normalize phenotype terms unseen in the fine-tuning corpus, including misspellings, synonyms, terms from other ontologies, and laymen's terms. Our approach provides a solution for the use of LLM to identify named medical entities from the clinical narratives, while successfully normalizing them to standard concepts in a controlled vocabulary.

A multiplexed, single-cell sequencing screen identifies compounds that increase neurogenic reprogramming of murine Muller glia

Retinal degeneration in mammals causes permanent loss of vision, due to an inability to regenerate naturally. Some non-mammalian vertebrates show robust regeneration, via Muller glia (MG). We have recently made significant progress in stimulating adult mouse MG to regenerate functional neurons by transgenic expression of the proneural transcription factor Ascl1. While these results showed that MG can serve as an endogenous source of neuronal replacement, the efficacy of this process is limited. With the goal of improving this in mammals, we designed a small molecule screen using sci-Plex, a method to multiplex up to thousands of single nucleus RNA-seq conditions into a single experiment. We used this technology to screen a library of 92 compounds, identified, and validated two that promote neurogenesis in vivo. Our results demonstrate that high-throughput single-cell molecular profiling can substantially improve the discovery process for molecules and pathways that can stimulate neural regeneration and further demonstrate the potential for this approach to restore vision in patients with retinal disease.

Hierarchical joint analysis of marginal summary statistics-Part I: Multipopulation fine mapping and credible set construction

Recent advancement in genome-wide association studies (GWAS) comes from not only increasingly larger sample sizes but also the shift in focus towards underrepresented populations. Multipopulation GWAS increase power to detect novel risk variants and improve fine-mapping resolution by leveraging evidence and differences in linkage disequilibrium (LD) from diverse populations. Here, we expand upon our previous approach for single-population fine-mapping through Joint Analysis of Marginal SNP Effects (JAM) to a multipopulation analysis (mJAM). Under the assumption that true causal variants are common across studies, we implement a hierarchical model framework that conditions on multiple SNPs while explicitly incorporating the different LD structures across populations. The mJAM framework can be used to first select index variants using the mJAM likelihood with different feature selection approaches. In addition, we present a novel approach leveraging the ideas of mediation to construct credible sets for these index variants. Construction of such credible sets can be performed given any existing index variants. We illustrate the implementation of the mJAM likelihood through two implementations: mJAM-SuSiE (a Bayesian approach) and mJAM-Forward selection. Through simulation studies based on realistic effect sizes and levels of LD, we demonstrated that mJAM performs well for constructing concise credible sets that include the underlying causal variants. In real data examples taken from the most recent multipopulation prostate cancer GWAS, we showed several practical advantages of mJAM over other existing multipopulation methods.

Evaluating and improving the representation of bacterial contents in long-read metagenome assemblies

BACKGROUND

In the metagenomic assembly of a microbial community, abundant species are often thought to assemble well given their deeper sequencing coverage. This conjuncture is rarely tested or evaluated in practice. We often do not know how many abundant species are missing and do not have an approach to recover them.

RESULTS

Here, we propose k-mer based and 16S RNA based methods to measure the completeness of metagenome assembly. We show that even with PacBio high-fidelity (HiFi) reads, abundant species are often not assembled, as high strain diversity may lead to fragmented contigs. We develop a novel reference-free algorithm to recover abundant metagenome-assembled genomes (MAGs) by identifying circular assembly subgraphs. Complemented with a reference-free genome binning heuristics based on dimension reduction, the proposed method rescues many abundant species that would be missing with existing methods and produces competitive results compared to those state-of-the-art binners in terms of total number of near-complete genome bins.

CONCLUSIONS

Our work emphasizes the importance of metagenome completeness, which has often been overlooked. Our algorithm generates more circular MAGs and moves a step closer to the complete representation of microbial communities.

Heterozygous loss-of-function SMC3 variants are associated with variable growth and developmental features

Heterozygous missense variants and in-frame indels in SMC3 are a cause of Cornelia de Lange syndrome (CdLS), marked by intellectual disability, growth deficiency, and dysmorphism, via an apparent dominant-negative mechanism. However, the spectrum of manifestations associated with SMC3 loss-of-function variants has not been reported, leading to hypotheses of alternative phenotypes or even developmental lethality. We used matchmaking servers, patient registries, and other resources to identify individuals with heterozygous, predicted loss-of-function (pLoF) variants in SMC3, and analyzed population databases to characterize mutational intolerance in this gene. Here, we show that SMC3 behaves as an archetypal haploinsufficient gene: it is highly constrained against pLoF variants, strongly depleted for missense variants, and pLoF variants are associated with a range of developmental phenotypes. Among 14 individuals with SMC3 pLoF variants, phenotypes were variable but coalesced on low growth parameters, developmental delay/intellectual disability, and dysmorphism, reminiscent of atypical CdLS. Comparisons to individuals with SMC3 missense/in-frame indel variants demonstrated an overall milder presentation in pLoF carriers. Furthermore, several individuals harboring pLoF variants in SMC3 were nonpenetrant for growth, developmental, and/or dysmorphic features, and some had alternative symptomatologies with rational biological links to SMC3. Analyses of tumor and model system transcriptomic data and epigenetic data in a subset of cases suggest that SMC3 pLoF variants reduce SMC3 expression but do not strongly support clustering with functional genomic signatures of typical CdLS. Our finding of substantial population-scale LoF intolerance in concert with variable growth and developmental features in subjects with SMC3 pLoF variants expands the scope of cohesinopathies, informs on their allelic architecture, and suggests the existence of additional clearly LoF-constrained genes whose disease links will be confirmed only by multilayered genomic data paired with careful phenotyping.

Evaluating and improving health equity and fairness of polygenic scores

Polygenic scores (PGSs) are quantitative metrics for predicting phenotypic values, such as human height or disease status. Some PGS methods require only summary statistics of a relevant genome-wide association study (GWAS) for their score. One such method is Lassosum, which inherits the model selection advantages of Lasso to select a meaningful subset of the GWAS single-nucleotide polymorphisms as predictors from their association statistics. However, even efficient scores like Lassosum, when derived from European-based GWASs, are poor predictors of phenotype for subjects of non-European ancestry; that is, they have limited portability to other ancestries. To increase the portability of Lassosum, when GWAS information and estimates of linkage disequilibrium are available for both ancestries, we propose Joint-Lassosum (JLS). In the simulation settings we explore, JLS provides more accurate PGSs compared to other methods, especially when measured in terms of fairness. In analyses of UK Biobank data, JLS was computationally more efficient but slightly less accurate than a Bayesian comparator, SDPRX. Like all PGS methods, JLS requires selection of predictors, which are determined by data-driven tuning parameters. We describe a new approach to selecting tuning parameters and note its relevance for model selection for any PGS. We also draw connections to the literature on algorithmic fairness and discuss how JLS can help mitigate fairness-related harms that might result from the use of PGSs in clinical settings. While no PGS method is likely to be universally portable, due to the diversity of human populations and unequal information content of GWASs for different ancestries, JLS is an effective approach for enhancing portability and reducing predictive bias.

An open source knowledge graph ecosystem for the life sciences

Translational research requires data at multiple scales of biological organization. Advancements in sequencing and multi-omics technologies have increased the availability of these data, but researchers face significant integration challenges. Knowledge graphs (KGs) are used to model complex phenomena, and methods exist to construct them automatically. However, tackling complex biomedical integration problems requires flexibility in the way knowledge is modeled. Moreover, existing KG construction methods provide robust tooling at the cost of fixed or limited choices among knowledge representation models. PheKnowLator (Phenotype Knowledge Translator) is a semantic ecosystem for automating the FAIR (Findable, Accessible, Interoperable, and Reusable) construction of ontologically grounded KGs with fully customizable knowledge representation. The ecosystem includes KG construction resources (e.g., data preparation APIs), analysis tools (e.g., SPARQL endpoint resources and abstraction algorithms), and benchmarks (e.g., prebuilt KGs). We evaluated the ecosystem by systematically comparing it to existing open-source KG construction methods and by analyzing its computational performance when used to construct 12 different large-scale KGs. With flexible knowledge representation, PheKnowLator enables fully customizable KGs without compromising performance or usability.

Return of polygenic risk scores in research: Stakeholders' views on the eMERGE-IV study

Research on polygenic risk scores (PRSs) for common, genetically complex chronic diseases aims to improve health-related predictions, tailor risk-reducing interventions, and improve health outcomes. Yet, the study and use of PRSs in clinical settings raise equity, clinical, and regulatory challenges that can be greater for individuals from historically marginalized racial, ethnic, and other minoritized communities. As part of the National Human Genome Research Institute-funded Electronic Medical Records and Genomics IV Network, we conducted online focus groups with patients/community members, clinicians, and members of institutional review boards to explore their views on key issues, including PRS research, return of PRS results, clinical translation, and barriers and facilitators to health behavioral changes in response to PRS results. Across stakeholder groups, our findings indicate support for PRS development and a strong interest in having PRS results returned to research participants. However, we also found multi-level barriers and significant differences in stakeholders' views about what is needed and possible for successful implementation. These include researcher-participant interaction formats, health and genomic literacy, and a range of structural barriers, such as financial instability, insurance coverage, and the absence of health-supporting infrastructure and affordable healthy food options in poorer neighborhoods. Our findings highlight the need to revisit and implement measures in PRS studies (e.g., incentives and resources for follow-up care), as well as system-level policies to promote equity in genomic research and health outcomes.

The impact of package selection and versioning on single-cell RNA-seq analysis

Standard single-cell RNA-sequencing analysis (scRNA-seq) workflows consist of converting raw read data into cell-gene count matrices through sequence alignment, followed by analyses including filtering, highly variable gene selection, dimensionality reduction, clustering, and differential expression analysis. Seurat and Scanpy are the most widely-used packages implementing such workflows, and are generally thought to implement individual steps similarly. We investigate in detail the algorithms and methods underlying Seurat and Scanpy and find that there are, in fact, considerable differences in the outputs of Seurat and Scanpy. The extent of differences between the programs is approximately equivalent to the variability that would be introduced in benchmarking scRNA-seq datasets by sequencing less than 5% of the reads or analyzing less than 20% of the cell population. Additionally, distinct versions of Seurat and Scanpy can produce very different results, especially during parts of differential expression analysis. Our analysis highlights the need for users of scRNA-seq to carefully assess the tools on which they rely, and the importance of developers of scientific software to prioritize transparency, consistency, and reproducibility for their tools.

Probing enzyme-dependent pseudouridylation using direct RNA sequencing to assess neuronal epitranscriptome plasticity

Chemical modifications in mRNAs such as pseudouridine (psi) can regulate gene expression, although our understanding of the functional impact of individual psi modifications, especially in neuronal cells, is limited. We apply nanopore direct RNA sequencing to investigate psi dynamics under cellular perturbations in SH-SY5Y cells. We assign sites to psi synthases using siRNA-based knockdown. A steady-state enzyme-substrate model reveals a strong correlation between psi synthase and mRNA substrate levels and psi modification frequencies. Next, we performed either differentiation or lead-exposure to SH-SY5Y cells and found that, upon lead exposure, not differentiation, the modification frequency is less dependent on enzyme levels suggesting translational control. Finally, we compared the plasticity of psi sites across cellular states and found that plastic sites can be condition-dependent or condition-independent; several of these sites fall within transcripts encoding proteins involved in neuronal processes. Our psi analysis and validation enable investigations into the dynamics and plasticity of RNA modifications.

An Ensemble Penalized Regression Method for Multi-ancestry Polygenic Risk Prediction

Great efforts are being made to develop advanced polygenic risk scores (PRS) to improve the prediction of complex traits and diseases. However, most existing PRS are primarily trained on European ancestry populations, limiting their transferability to non-European populations. In this article, we propose a novel method for generating multi-ancestry Polygenic Risk scOres based on enSemble of PEnalized Regression models (PROSPER). PROSPER integrates genome-wide association studies (GWAS) summary statistics from diverse populations to develop ancestry-specific PRS with improved predictive power for minority populations. The method uses a combination of ℒ 1 (lasso) and ℒ 2 (ridge) penalty functions, a parsimonious specification of the penalty parameters across populations, and an ensemble step to combine PRS generated across different penalty parameters. We evaluate the performance of PROSPER and other existing methods on large-scale simulated and real datasets, including those from 23andMe Inc., the Global Lipids Genetics Consortium, and All of Us. Results show that PROSPER can substantially improve multi-ancestry polygenic prediction compared to alternative methods across a wide variety of genetic architectures. In real data analyses, for example, PROSPER increased out-of-sample prediction R2 for continuous traits by an average of 70% compared to a state-of-the-art Bayesian method (PRS-CSx) in the African ancestry population. Further, PROSPER is computationally highly scalable for the analysis of large SNP contents and many diverse populations.

Widespread variation in molecular interactions and regulatory properties among transcription factor isoforms

Most human Transcription factors (TFs) genes encode multiple protein isoforms differing in DNA binding domains, effector domains, or other protein regions. The global extent to which this results in functional differences between isoforms remains unknown. Here, we systematically compared 693 isoforms of 246 TF genes, assessing DNA binding, protein binding, transcriptional activation, subcellular localization, and condensate formation. Relative to reference isoforms, two-thirds of alternative TF isoforms exhibit differences in one or more molecular activities, which often could not be predicted from sequence. We observed two primary categories of alternative TF isoforms: "rewirers" and "negative regulators", both of which were associated with differentiation and cancer. Our results support a model wherein the relative expression levels of, and interactions involving, TF isoforms add an understudied layer of complexity to gene regulatory networks, demonstrating the importance of isoform-aware characterization of TF functions and providing a rich resource for further studies.

Longitudinal profiling of the microbiome at four body sites reveals core stability and individualized dynamics during health and disease

To understand the dynamic interplay between the human microbiome and host during health and disease, we analyzed the microbial composition, temporal dynamics, and associations with host multi-omics, immune, and clinical markers of microbiomes from four body sites in 86 participants over 6 years. We found that microbiome stability and individuality are body-site specific and heavily influenced by the host. The stool and oral microbiome are more stable than the skin and nasal microbiomes, possibly due to their interaction with the host and environment. We identify individual-specific and commonly shared bacterial taxa, with individualized taxa showing greater stability. Interestingly, microbiome dynamics correlate across body sites, suggesting systemic dynamics influenced by host-microbial-environment interactions. Notably, insulin-resistant individuals show altered microbial stability and associations among microbiome, molecular markers, and clinical features, suggesting their disrupted interaction in metabolic disease. Our study offers comprehensive views of multi-site microbial dynamics and their relationship with host health and disease.

Pfp-fm: an accelerated FM-index

FM-indexes are crucial data structures in DNA alignment, but searching with them usually takes at least one random access per character in the query pattern. Ferragina and Fischer [1] observed in 2007 that word-based indexes often use fewer random accesses than character-based indexes, and thus support faster searches. Since DNA lacks natural word-boundaries, however, it is necessary to parse it somehow before applying word-based FM-indexing. In 2022, Deng et al. [2] proposed parsing genomic data by induced suffix sorting, and showed that the resulting word-based FM-indexes support faster counting queries than standard FM-indexes when patterns are a few thousand characters or longer. In this paper we show that using prefix-free parsing-which takes parameters that let us tune the average length of the phrases-instead of induced suffix sorting, gives a significant speedup for patterns of only a few hundred characters. We implement our method and demonstrate it is between 3 and 18 times faster than competing methods on queries to GRCh38, and is consistently faster on queries made to 25,000, 50,000 and 100,000 SARS-CoV-2 genomes. Hence, it seems our method accelerates the performance of count over all state-of-the-art methods with a moderate increase in the memory. The source code for PFP - FM is available at https://github.com/AaronHong1024/afm .

Bayesian estimation of gene constraint from an evolutionary model with gene features

Measures of selective constraint on genes have been used for many applications including clinical interpretation of rare coding variants, disease gene discovery, and studies of genome evolution. However, widely-used metrics are severely underpowered at detecting constraint for the shortest ∼25% of genes, potentially causing important pathogenic mutations to be overlooked. We developed a framework combining a population genetics model with machine learning on gene features to enable accurate inference of an interpretable constraint metric, shet. Our estimates outperform existing metrics for prioritizing genes important for cell essentiality, human disease, and other phenotypes, especially for short genes. Our new estimates of selective constraint should have wide utility for characterizing genes relevant to human disease. Finally, our inference framework, GeneBayes, provides a flexible platform that can improve estimation of many gene-level properties, such as rare variant burden or gene expression differences.

DNA binding analysis of rare variants in homeodomains reveals homeodomain specificity-determining residues

Homeodomains (HDs) are the second largest class of DNA binding domains (DBDs) among eukaryotic sequence-specific transcription factors (TFs) and are the TF structural class with the largest number of disease-associated mutations in the Human Gene Mutation Database (HGMD). Despite numerous structural studies and large-scale analyses of HD DNA binding specificity, HD-DNA recognition is still not fully understood. Here, we analyze 92 human HD mutants, including disease-associated variants and variants of uncertain significance (VUS), for their effects on DNA binding activity. Many of the variants alter DNA binding affinity and/or specificity. Detailed biochemical analysis and structural modeling identifies 14 previously unknown specificity-determining positions, 5 of which do not contact DNA. The same missense substitution at analogous positions within different HDs often exhibits different effects on DNA binding activity. Variant effect prediction tools perform moderately well in distinguishing variants with altered DNA binding affinity, but poorly in identifying those with altered binding specificity. Our results highlight the need for biochemical assays of TF coding variants and prioritize dozens of variants for further investigations into their pathogenicity and the development of clinical diagnostics and precision therapies.

MUSSEL: Enhanced Bayesian polygenic risk prediction leveraging information across multiple ancestry groups

Polygenic risk scores (PRSs) are now showing promising predictive performance on a wide variety of complex traits and diseases, but there exists a substantial performance gap across populations. We propose MUSSEL, a method for ancestry-specific polygenic prediction that borrows information in summary statistics from genome-wide association studies (GWASs) across multiple ancestry groups via Bayesian hierarchical modeling and ensemble learning. In our simulation studies and data analyses across four distinct studies, totaling 5.7 million participants with a substantial ancestral diversity, MUSSEL shows promising performance compared to alternatives. For example, MUSSEL has an average gain in prediction R2 across 11 continuous traits of 40.2% and 49.3% compared to PRS-CSx and CT-SLEB, respectively, in the African ancestry population. The best-performing method, however, varies by GWAS sample size, target ancestry, trait architecture, and linkage disequilibrium reference samples; thus, ultimately a combination of methods may be needed to generate the most robust PRSs across diverse populations.

A molecular proximity sensor based on an engineered, dual-component guide RNA

One of the goals of synthetic biology is to enable the design of arbitrary molecular circuits with programmable inputs and outputs. Such circuits bridge the properties of electronic and natural circuits, processing information in a predictable manner within living cells. Genome editing is a potentially powerful component of synthetic molecular circuits, whether for modulating the expression of a target gene or for stably recording information to genomic DNA. However, programming molecular events such as protein-protein interactions or induced proximity as triggers for genome editing remains challenging. Here we demonstrate a strategy termed "P3 editing", which links protein-protein proximity to the formation of a functional CRISPRCas9 dual-component guide RNA. By engineering the crRNA:tracrRNA interaction, we demonstrate that various known protein-protein interactions, as well as the chemically-induced dimerization of protein domains, can be used to activate prime editing or base editing in human cells. Additionally, we explore how P3 editing can incorporate outputs from ADAR-based RNA sensors, potentially allowing specific RNAs to induce specific genome edits within a larger circuit. Our strategy enhances the controllability of CRISPR-based genome editing, facilitating its use in synthetic molecular circuits deployed in living cells.

A revamped rat reference genome improves the discovery of genetic diversity in laboratory rats

The seventh iteration of the reference genome assembly for Rattus norvegicus-mRatBN7.2-corrects numerous misplaced segments and reduces base-level errors by approximately 9-fold and increases contiguity by 290-fold compared with its predecessor. Gene annotations are now more complete, improving the mapping precision of genomic, transcriptomic, and proteomics datasets. We jointly analyzed 163 short-read whole-genome sequencing datasets representing 120 laboratory rat strains and substrains using mRatBN7.2. We defined ∼20.0 million sequence variations, of which 18,700 are predicted to potentially impact the function of 6,677 genes. We also generated a new rat genetic map from 1,893 heterogeneous stock rats and annotated transcription start sites and alternative polyadenylation sites. The mRatBN7.2 assembly, along with the extensive analysis of genomic variations among rat strains, enhances our understanding of the rat genome, providing researchers with an expanded resource for studies involving rats.