TMPRSS2 and glycan receptors synergistically facilitate coronavirus entry.

The entry of coronaviruses is initiated by spike recognition of host cellular receptors, involving proteinaceous and/or glycan receptors. Recently, TMPRSS2 was identified as the proteinaceous receptor for HCoV-HKU1 alongside sialoglycan as a glycan receptor. However, the underlying mechanisms for viral entry remain unknown. Here, we investigated the HCoV-HKU1C spike in the inactive, glycan-activated, and functionally anchored states, revealing that sialoglycan binding induces a conformational change of the NTD and promotes the neighboring RBD of the spike to open for TMPRSS2 recognition, exhibiting a synergistic mechanism for the entry of HCoV-HKU1. The RBD of HCoV-HKU1 features an insertion subdomain that recognizes TMPRSS2 through three previously undiscovered interfaces. Furthermore, structural investigation of HCoV-HKU1A in combination with mutagenesis and binding assays confirms a conserved receptor recognition pattern adopted by HCoV-HKU1. These studies advance our understanding of the complex viral-host interactions during entry, laying the groundwork for developing new therapeutics against coronavirus-associated diseases.

Structural RNA components supervise the sequential DNA cleavage in R2 retrotransposon.

Retroelements are the widespread jumping elements considered as major drivers for genome evolution, which can also be repurposed as gene-editing tools. Here, we determine the cryo-EM structures of eukaryotic R2 retrotransposon with ribosomal DNA target and regulatory RNAs. Combined with biochemical and sequencing analysis, we reveal two essential DNA regions, Drr and Dcr, required for recognition and cleavage. The association of 3' regulatory RNA with R2 protein accelerates the first-strand cleavage, blocks the second-strand cleavage, and initiates the reverse transcription starting from the 3'-tail. Removing 3' regulatory RNA by reverse transcription allows the association of 5' regulatory RNA and initiates the second-strand cleavage. Taken together, our work explains the DNA recognition and RNA supervised sequential retrotransposition mechanisms by R2 machinery, providing insights into the retrotransposon and application reprogramming.

Broad neutralization of SARS-CoV-2 variants by an inhalable bispecific single-domain antibody.

The effectiveness of SARS-CoV-2 vaccines and therapeutic antibodies have been limited by the continuous emergence of viral variants and by the restricted diffusion of antibodies from circulation into the sites of respiratory virus infection. Here, we report the identification of two highly conserved regions on the Omicron variant receptor-binding domain recognized by broadly neutralizing antibodies. Furthermore, we generated a bispecific single-domain antibody that was able to simultaneously and synergistically bind these two regions on a single Omicron variant receptor-binding domain as revealed by cryo-EM structures. We demonstrated that this bispecific antibody can be effectively delivered to lung via inhalation administration and exhibits exquisite neutralization breadth and therapeutic efficacy in mouse models of SARS-CoV-2 infections. Importantly, this study also deciphered an uncommon and highly conserved cryptic epitope within the spike trimeric interface that may have implications for the design of broadly protective SARS-CoV-2 vaccines and therapeutics.

In vivo structural characterization of the SARS-CoV-2 RNA genome identifies host proteins vulnerable to repurposed drugs.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the cause of the ongoing coronavirus disease 2019 (COVID-19) pandemic. Understanding of the RNA virus and its interactions with host proteins could improve therapeutic interventions for COVID-19. By using icSHAPE, we determined the structural landscape of SARS-CoV-2 RNA in infected human cells and from refolded RNAs, as well as the regulatory untranslated regions of SARS-CoV-2 and six other coronaviruses. We validated several structural elements predicted in silico and discovered structural features that affect the translation and abundance of subgenomic viral RNAs in cells. The structural data informed a deep-learning tool to predict 42 host proteins that bind to SARS-CoV-2 RNA. Strikingly, antisense oligonucleotides targeting the structural elements and FDA-approved drugs inhibiting the SARS-CoV-2 RNA binding proteins dramatically reduced SARS-CoV-2 infection in cells derived from human liver and lung tumors. Our findings thus shed light on coronavirus and reveal multiple candidate therapeutics for COVID-19 treatment.

Cryo-EM Structure of Human Dicer and Its Complexes with a Pre-miRNA Substrate.

Human Dicer (hDicer) is a multi-domain protein belonging to the RNase III family. It plays pivotal roles in small RNA biogenesis during the RNA interference (RNAi) pathway by processing a diverse range of double-stranded RNA (dsRNA) precursors to generate ∼22 nt microRNA (miRNA) or small interfering RNA (siRNA) products for sequence-directed gene silencing. In this work, we solved the cryoelectron microscopy (cryo-EM) structure of hDicer in complex with its cofactor protein TRBP and revealed the precise spatial arrangement of hDicer's multiple domains. We further solved structures of the hDicer-TRBP complex bound with pre-let-7 RNA in two distinct conformations. In combination with biochemical analysis, these structures reveal a property of the hDicer-TRBP complex to promote the stability of pre-miRNA's stem duplex in a pre-dicing state. These results provide insights into the mechanism of RNA processing by hDicer and illustrate the regulatory role of hDicer's N-terminal helicase domain.

IL-13 secreted by ILC2s promotes the self-renewal of intestinal stem cells through circular RNA circPan3.

Intestinal stem cells (ISCs) are maintained by stemness signaling for precise modulation of self-renewal and differentiation under homeostasis. However, the way in which intestinal immune cells regulate the self-renewal of ISCs remains elusive. Here we found that mouse and human Lgr5+ ISCs showed high expression of the immune cell-associated circular RNA circPan3 (originating from the Pan3 gene transcript). Deletion of circPan3 in Lgr5+ ISCs impaired their self-renewal capacity and the regeneration of gut epithelium in a manner dependent on immune cells. circPan3 bound mRNA encoding the cytokine IL-13 receptor subunit IL-13Rα1 (Il13ra1) in ISCs to increase its stability, which led to the expression of IL-13Rα1 in ISCs. IL-13 produced by group 2 innate lymphoid cells in the crypt niche engaged IL-13Rα1 on crypt ISCs and activated signaling mediated by IL-13‒IL-13R, which in turn initiated expression of the transcription factor Foxp1. Foxp1 is associated with ß-catenin in rendering its nuclear translocation, which caused activation of the ß-catenin pathway and the maintenance of Lgr5+ ISCs.

CXCR5(+) follicular cytotoxic T cells control viral infection in B cell follicles.

During unresolved infections, some viruses escape immunological control and establish a persistant reservoir in certain cell types, such as human immunodeficiency virus (HIV), which persists in follicular helper T cells (TFH cells), and Epstein-Barr virus (EBV), which persists in B cells. Here we identified a specialized group of cytotoxic T cells (TC cells) that expressed the chemokine receptor CXCR5, selectively entered B cell follicles and eradicated infected TFH cells and B cells. The differentiation of these cells, which we have called 'follicular cytotoxic T cells' (TFC cells), required the transcription factors Bcl6, E2A and TCF-1 but was inhibited by the transcriptional regulators Blimp1, Id2 and Id3. Blimp1 and E2A directly regulated Cxcr5 expression and, together with Bcl6 and TCF-1, formed a transcriptional circuit that guided TFC cell development. The identification of TFC cells has far-reaching implications for the development of strategies to control infections that target B cells and TFH cells and to treat B cell-derived malignancies.

Coherent interfaces govern direct transformation from graphite to diamond.

Understanding the direct transformation from graphite to diamond has been a long-standing challenge with great scientific and practical importance. Previously proposed transformation mechanisms1-3, based on traditional experimental observations that lacked atomistic resolution, cannot account for the complex nanostructures occurring at graphite-diamond interfaces during the transformation4,5. Here we report the identification of coherent graphite-diamond interfaces, which consist of four basic structural motifs, in partially transformed graphite samples recovered from static compression, using high-angle annular dark-field scanning transmission electron microscopy. These observations provide insight into possible pathways of the transformation. Theoretical calculations confirm that transformation through these coherent interfaces is energetically favoured compared with those through other paths previously proposed1-3. The graphite-to-diamond transformation is governed by the formation of nanoscale coherent interfaces (diamond nucleation), which, under static compression, advance to consume the remaining graphite (diamond growth). These results may also shed light on transformation mechanisms of other carbon materials and boron nitride under different synthetic conditions.

Coordination of germ layer lineage choice by TET1 during primed pluripotency.

Gastrulation in the early postimplantation stage mammalian embryo begins when epiblast cells ingress to form the primitive streak or develop as the embryonic ectoderm. The DNA dioxygenase Tet1 is highly expressed in the epiblast and yet continues to regulate lineage specification during gastrulation when its expression is diminished. Here, we show how Tet1 plays a pivotal role upstream of germ layer lineage bifurcation. During the transition from naive pluripotency to lineage priming, a global reconfiguration redistributes Tet1 from Oct4-cobound promoters to distal regulatory elements at lineage differentiation genes, which are distinct from high-affinity sites engaged by Oct4. An altered chromatin landscape in Tet1-deficient primed epiblast-like cells is associated with enhanced Oct4 expression and binding to Nodal and Wnt target genes, resulting in collaborative signals that enhance mesendodermal and inhibit neuroectodermal gene expression during lineage segregation. A permissive role for Tet1 in neural fate induction involves Zic2-dependent engagement at neural target genes at lineage priming, is dependent on the signaling environment during gastrulation, and impacts neural tube closure after gastrulation. Our findings provide mechanistic information for epigenetic integration of pluripotency and signal-induced differentiation cues.

The intra-S phase checkpoint targets Dna2 to prevent stalled replication forks from reversing.

When replication forks stall at damaged bases or upon nucleotide depletion, the intra-S phase checkpoint ensures they are stabilized and can restart. In intra-S checkpoint-deficient budding yeast, stalling forks collapse, and ∼10% form pathogenic chicken foot structures, contributing to incomplete replication and cell death (Lopes et al., 2001; Sogo et al., 2002; Tercero and Diffley, 2001). Using fission yeast, we report that the Cds1(Chk2) effector kinase targets Dna2 on S220 to regulate, both in vivo and in vitro, Dna2 association with stalled replication forks in chromatin. We demonstrate that Dna2-S220 phosphorylation and the nuclease activity of Dna2 are required to prevent fork reversal. Consistent with this, Dna2 can efficiently cleave obligate precursors of fork regression-regressed leading or lagging strands-on model replication forks. We propose that Dna2 cleavage of regressed nascent strands prevents fork reversal and thus stabilizes stalled forks to maintain genome stability during replication stress.

Better together against genetic heterogeneity: A sex-combined joint main and interaction analysis of 290 quantitative traits in the UK Biobank.

Genetic effects can be sex-specific, particularly for traits such as testosterone, a sex hormone. While sex-stratified analysis provides easily interpretable sex-specific effect size estimates, the presence of sex-differences in SNP effect implies a SNP×sex interaction. This suggests the usage of the often overlooked joint test, testing for an SNP's main and SNP×sex interaction effects simultaneously. Notably, even without individual-level data, the joint test statistic can be derived from sex-stratified summary statistics through an omnibus meta-analysis. Utilizing the available sex-stratified summary statistics of the UK Biobank, we performed such omnibus meta-analyses for 290 quantitative traits. Results revealed that this approach is robust to genetic effect heterogeneity and can outperform the traditional sex-stratified or sex-combined main effect-only tests. Therefore, we advocate using the omnibus meta-analysis that captures both the main and interaction effects. Subsequent sex-stratified analysis should be conducted for sex-specific effect size estimation and interpretation.

Canadian COVID-19 host genetics cohort replicates known severity associations.

The HostSeq initiative recruited 10,059 Canadians infected with SARS-CoV-2 between March 2020 and March 2023, obtained clinical information on their disease experience and whole genome sequenced (WGS) their DNA. We analyzed the WGS data for genetic contributors to severe COVID-19 (considering 3,499 hospitalized cases and 4,975 non-hospitalized after quality control). We investigated the evidence for replication of loci reported by the International Host Genetics Initiative (HGI); analyzed the X chromosome; conducted rare variant gene-based analysis and polygenic risk score testing. Population stratification was adjusted for using meta-analysis across ancestry groups. We replicated two loci identified by the HGI for COVID-19 severity: the LZTFL1/SLC6A20 locus on chromosome 3 and the FOXP4 locus on chromosome 6 (the latter with a variant significant at P < 5E-8). We found novel significant associations with MRAS and WDR89 in gene-based analyses, and constructed a polygenic risk score that explained 1.01% of the variance in severe COVID-19. This study provides independent evidence confirming the robustness of previously identified COVID-19 severity loci by the HGI and identifies novel genes for further investigation.

Comprehensive whole-genome analyses of the UK Biobank reveal significant sex differences in both genotype missingness and allele frequency on the X chromosome.

The UK Biobank is the most used dataset for genome-wide association studies (GWAS). GWAS of sex, essentially sex differences in minor allele frequencies (sdMAF), has identified autosomal SNPs with significant sdMAF, including in the UK Biobank, but the X chromosome was excluded. Our recent report identified multiple regions on the X chromosome with significant sdMAF, using short-read sequencing of other datasets. We performed a whole genome sdMAF analysis, with ~410 k white British individuals from the UK Biobank, using array genotyped, imputed or exome sequencing data. We observed marked sdMAF on the X chromosome, particularly at the boundaries between the pseudo-autosomal regions (PAR) and the non-PAR (NPR), as well as throughout the NPR, consistent with our earlier report. A small fraction of autosomal SNPs also showed significant sdMAF. Using the centrally imputed data, which relied mostly on low-coverage whole genome sequence, resulted in 2.1% of NPR SNPs with significant sdMAF. The whole exome sequencing also displays sdMAF on the X chromosome, including some NPR SNPs with heterozygous genotype calls in males. Genotyping, sequencing and imputation of X chromosomal SNPs requires further attention to ensure the integrity for downstream association analysis.

Association of whole-person eigen-polygenic risk scores with Alzheimer's disease.

Late-Onset Alzheimer's Disease (LOAD) is a heterogeneous neurodegenerative disorder with complex etiology and high heritability. Its multifactorial risk profile and large portions of unexplained heritability suggest the involvement of yet unidentified genetic risk factors. Here we describe the "whole person" genetic risk landscape of polygenic risk scores for 2218 traits in 2044 elderly individuals and test if novel eigen-PRSs derived from clustered subnetworks of single-trait PRSs can improve the prediction of LOAD diagnosis, rates of cognitive decline, and canonical LOAD neuropathology. Network analyses revealed distinct clusters of PRSs with clinical and biological interpretability. Novel eigen-PRSs (ePRS) from these clusters significantly improved LOAD-related phenotypes prediction over current state-of-the-art LOAD PRS models. Notably, an ePRS representing clusters of traits related to cholesterol levels was able to improve variance explained in a model of the brain-wide beta-amyloid burden by 1.7% (likelihood ratio test P = 9.02 × 10-7). All associations of ePRS with LOAD phenotypes were eliminated by the removal of APOE-proximal loci. However, our association analysis identified modules characterized by PRSs of high cholesterol and LOAD. We believe this is due to the influence of the APOE region from both PRSs. We found significantly higher mean SNP effects for LOAD in the intersecting APOE region SNPs. Combining genetic risk factors for vascular traits and dementia could improve current single-trait PRS models of LOAD, enhancing the use of PRS in risk stratification. Our results are catalogued for the scientific community, to aid in generating new hypotheses based on our maps of clustered PRSs and associations with LOAD-related phenotypes.

eXclusionarY: 10 years later, where are the sex chromosomes in GWASs?

10 years ago, a detailed analysis showed that only 33% of genome-wide association study (GWAS) results included the X chromosome. Multiple recommendations were made to combat such exclusion. Here, we re-surveyed the research landscape to determine whether these earlier recommendations had been translated. Unfortunately, among the genome-wide summary statistics reported in 2021 in the NHGRI-EBI GWAS Catalog, only 25% provided results for the X chromosome and 3% for the Y chromosome, suggesting that the exclusion phenomenon not only persists but has also expanded into an exclusionary problem. Normalizing by physical length of the chromosome, the average number of studies published through November 2022 with genome-wide-significant findings on the X chromosome is ∼1 study/Mb. By contrast, it ranges from ∼6 to ∼16 studies/Mb for chromosomes 4 and 19, respectively. Compared with the autosomal growth rate of ∼0.086 studies/Mb/year over the last decade, studies of the X chromosome grew at less than one-seventh that rate, only ∼0.012 studies/Mb/year. Among the studies that reported significant associations on the X chromosome, we noted extreme heterogeneities in data analysis and reporting of results, suggesting the need for clear guidelines. Unsurprisingly, among the 430 scores sampled from the PolyGenic Score Catalog, 0% contained weights for sex chromosomal SNPs. To overcome the dearth of sex chromosome analyses, we provide five sets of recommendations and future directions. Finally, until the sex chromosomes are included in a whole-genome study, instead of GWASs, we propose such studies would more properly be referred to as "AWASs," meaning "autosome-wide scans."

Restoring neuronal iron homeostasis revitalizes neurogenesis after spinal cord injury.

Spinal cord injury (SCI) can lead to iron overloading and subsequent neuronal ferroptosis, which hinders the recovery of locomotor function. However, it is still unclear whether the maintenance of neuronal iron homeostasis enables to revitalize intrinsic neurogenesis. Herein, we report the regulation of cellular iron homeostasis after SCI via the chelation of excess iron ions and modulation of the iron transportation pathway using polyphenol-based hydrogels for the revitalization of intrinsic neurogenesis. The reversed iron overloading can promote neural stem/progenitor cell differentiation into neurons and elicit the regenerative potential of newborn neurons, which is accompanied by improved axon reinnervation and remyelination. Notably, polyphenol-based hydrogels significantly increase the neurological motor scores from ~8 to 18 (out of 21) and restore the transmission of sensory and motor electrophysiological signals after SCI. Maintenance of iron homeostasis at the site of SCI using polyphenol-based hydrogels provides a promising paradigm to revitalize neurogenesis for the treatment of iron accumulation-related nervous system diseases.

Modulation of deep neural circuits with sonogenetics.

Noninvasive control of neuronal activity in the deep brain can be illuminating for probing brain function and treating dysfunctions. Here, we present a sonogenetic approach for controlling distinct mouse behavior with circuit specificity and subsecond temporal resolution. Targeted neurons in subcortical regions were made to express a mutant large conductance mechanosensitive ion channel (MscL-G22S), enabling ultrasound to trigger activity in MscL-expressing neurons in the dorsal striatum and increase locomotion in freely moving mice. Ultrasound stimulation of MscL-expressing neurons in the ventral tegmental area could activate the mesolimbic pathway to trigger dopamine release in the nucleus accumbens and modulate appetitive conditioning. Moreover, sonogenetic stimulation of the subthalamic nuclei of Parkinson's disease model mice improved their motor coordination and mobile time. Neuronal responses to ultrasound pulse trains were rapid, reversible, and repeatable. We also confirmed that the MscL-G22S mutant is more effective to sensitize neurons to ultrasound compared to the wild-type MscL. Altogether, we lay out a sonogenetic approach which can selectively manipulate targeted cells to activate defined neural pathways, affect specific behaviors, and relieve symptoms of neurodegenerative disease.

The mechanosensitive ion channel Piezo1 contributes to ultrasound neuromodulation.

Transcranial low-intensity ultrasound is a promising neuromodulation modality, with the advantages of noninvasiveness, deep penetration, and high spatiotemporal accuracy. However, the underlying biological mechanism of ultrasonic neuromodulation remains unclear, hindering the development of efficacious treatments. Here, the well-known Piezo1 was studied through a conditional knockout mouse model as a major mediator for ultrasound neuromodulation ex vivo and in vivo. We showed that Piezo1 knockout (P1KO) in the right motor cortex of mice significantly reduced ultrasound-induced neuronal calcium responses, limb movement, and muscle electromyogram (EMG) responses. We also detected higher Piezo1 expression in the central amygdala (CEA), which was found to be more sensitive to ultrasound stimulation than the cortex was. Knocking out the Piezo1 in CEA neurons showed a significant reduction of response under ultrasound stimulation, while knocking out astrocytic Piezo1 showed no-obvious changes in neuronal responses. Additionally, we excluded an auditory confound by monitoring auditory cortical activation and using smooth waveform ultrasound with randomized parameters to stimulate P1KO ipsilateral and contralateral regions of the same brain and recording evoked movement in the corresponding limb. Thus, we demonstrate that Piezo1 is functionally expressed in different brain regions and that it is an important mediator of ultrasound neuromodulation in the brain, laying the ground for further mechanistic studies of ultrasound.

A flexible summary statistics-based colocalization method with application to the mucin cystic fibrosis lung disease modifier locus.

Mucus obstruction is a central feature in the cystic fibrosis (CF) airways. A genome-wide association study (GWAS) of lung disease by the CF Gene Modifier Consortium (CFGMC) identified a significant locus containing two mucin genes, MUC20 and MUC4. Expression quantitative trait locus (eQTL) analysis using human nasal epithelia (HNE) from 94 CF-affected Canadians in the CFGMC demonstrated MUC4 eQTLs that mirrored the lung association pattern in the region, suggesting that MUC4 expression may mediate CF lung disease. Complications arose, however, with colocalization testing using existing methods: the locus is complex and the associated SNPs span a 0.2 Mb region with high linkage disequilibrium (LD) and evidence of allelic heterogeneity. We previously developed the Simple Sum (SS), a powerful colocalization test in regions with allelic heterogeneity, but SS assumed eQTLs to be present to achieve type I error control. Here we propose a two-stage SS (SS2) colocalization test that avoids a priori eQTL assumptions, accounts for multiple hypothesis testing and the composite null hypothesis, and enables meta-analysis. We compare SS2 to published approaches through simulation and demonstrate type I error control for all settings with the greatest power in the presence of high LD and allelic heterogeneity. Applying SS2 to the MUC20/MUC4 CF lung disease locus with eQTLs from CF HNE revealed significant colocalization with MUC4 (p = 1.31 × 10-5) rather than with MUC20. The SS2 is a powerful method to inform the responsible gene(s) at a locus and guide future functional studies. SS2 has been implemented in the application LocusFocus.