SARS-CoV-2 501Y.V2 variants lack higher infectivity but do have immune escape.

The 501Y.V2 variants of SARS-CoV-2 containing multiple mutations in spike are now dominant in South Africa and are rapidly spreading to other countries. Here, experiments with 18 pseudotyped viruses showed that the 501Y.V2 variants do not confer increased infectivity in multiple cell types except for murine ACE2-overexpressing cells, where a substantial increase in infectivity was observed. Notably, the susceptibility of the 501Y.V2 variants to 12 of 17 neutralizing monoclonal antibodies was substantially diminished, and the neutralization ability of the sera from convalescent patients and immunized mice was also reduced for these variants. The neutralization resistance was mainly caused by E484K and N501Y mutations in the receptor-binding domain of spike. The enhanced infectivity in murine ACE2-overexpressing cells suggests the possibility of spillover of the 501Y.V2 variants to mice. Moreover, the neutralization resistance we detected for the 501Y.V2 variants suggests the potential for compromised efficacy of monoclonal antibodies and vaccines.

The Impact of Mutations in SARS-CoV-2 Spike on Viral Infectivity and Antigenicity.

The spike protein of SARS-CoV-2 has been undergoing mutations and is highly glycosylated. It is critically important to investigate the biological significance of these mutations. Here, we investigated 80 variants and 26 glycosylation site modifications for the infectivity and reactivity to a panel of neutralizing antibodies and sera from convalescent patients. D614G, along with several variants containing both D614G and another amino acid change, were significantly more infectious. Most variants with amino acid change at receptor binding domain were less infectious, but variants including A475V, L452R, V483A, and F490L became resistant to some neutralizing antibodies. Moreover, the majority of glycosylation deletions were less infectious, whereas deletion of both N331 and N343 glycosylation drastically reduced infectivity, revealing the importance of glycosylation for viral infectivity. Interestingly, N234Q was markedly resistant to neutralizing antibodies, whereas N165Q became more sensitive. These findings could be of value in the development of vaccine and therapeutic antibodies.

Attenuated replication and pathogenicity of SARS-CoV-2 B.1.1.529 Omicron.

The Omicron (B.1.1.529) variant of SARS-CoV-2 emerged in November 2021 and is rapidly spreading among the human population1. Although recent reports reveal that the Omicron variant robustly escapes vaccine-associated and therapeutic neutralization antibodies2-10, the pathogenicity of the virus remains unknown. Here we show that the replication of Omicron is substantially attenuated in human Calu3 and Caco2 cells. Further mechanistic investigations reveal that Omicron is inefficient in its use of transmembrane serine protease 2 (TMPRSS2) compared with wild-type SARS-CoV-2 (HKU-001a) and previous variants, which may explain its reduced replication in Calu3 and Caco2 cells. The replication of Omicron is markedly attenuated in both the upper and lower respiratory tracts of infected K18-hACE2 mice compared with that of the wild-type strain and Delta (B.1.617.2) variant, resulting in its substantially ameliorated lung pathology. Compared with wild-type SARS-CoV-2 and the Alpha (B.1.1.7), Beta (1.351) and Delta variants, infection by Omicron causes the lowest reduction in body weight and the lowest mortality rate. Overall, our study demonstrates that the replication and pathogenicity of the Omicron variant of SARS-CoV-2 in mice is attenuated compared with the wild-type strain and other variants.

Engineering of Cas12a nuclease variants with enhanced genome-editing specificity.

The clustered regularly interspaced short palindromic repeat (CRISPR)-Cas12a system is a powerful tool in gene editing; however, crRNA-DNA mismatches might induce unwanted cleavage events, especially at the distal end of the PAM. To minimize this limitation, we engineered a hyper fidelity AsCas12a variant carrying the mutations S186A/R301A/T315A/Q1014A/K414A (termed HyperFi-As) by modifying amino acid residues interacting with the target DNA and crRNA strand. HyperFi-As retains on-target activities comparable to wild-type AsCas12a (AsCas12aWT) in human cells. We demonstrated that HyperFi-As has dramatically reduced off-target effects in human cells, and HyperFi-As possessed notably a lower tolerance to mismatch at the position of the PAM-distal region compared with the wild type. Further, a modified single-molecule DNA unzipping assay at proper constant force was applied to evaluate the stability and transient stages of the CRISPR/Cas ribonucleoprotein (RNP) complex. Multiple states were sensitively detected during the disassembly of the DNA-Cas12a-crRNA complexes. On off-target DNA substrates, the HyperFi-As-crRNA was harder to maintain the R-loop complex state compared to the AsCas12aWT, which could explain exactly why the HyperFi-As has low off-targeting effects in human cells. Our findings provide a novel version of AsCas12a variant with low off-target effects, especially capable of dealing with the high off-targeting in the distal region from the PAM. An insight into how the AsCas12a variant behaves at off-target sites was also revealed at the single-molecule level and the unzipping assay to evaluate multiple states of CRISPR/Cas RNP complexes might be greatly helpful for a deep understanding of how CRISPR/Cas behaves and how to engineer it in future.

Engineering of a compact, high-fidelity EbCas12a variant that can be packaged with its crRNA into an all-in-one AAV vector delivery system.

The CRISPR-associated endonuclease Cas12a has become a powerful genome-editing tool in biomedical research due to its ease of use and low off-targeting. However, the size of Cas12a severely limits clinical applications such as adeno-associated virus (AAV)-based gene therapy. Here, we characterized a novel compact Cas12a ortholog, termed EbCas12a, from the metagenome-assembled genome of a currently unclassified Erysipelotrichia. It has the PAM sequence of 5'-TTTV-3' (V = A, G, C) and the smallest size of approximately 3.47 kb among the Cas12a orthologs reported so far. In addition, enhanced EbCas12a (enEbCas12a) was also designed to have comparable editing efficiency with higher specificity to AsCas12a and LbCas12a in mammalian cells at multiple target sites. Based on the compact enEbCas12a, an all-in-one AAV delivery system with crRNA for Cas12a was developed for both in vitro and in vivo applications. Overall, the novel smallest high-fidelity enEbCas12a, this first case of the all-in-one AAV delivery for Cas12a could greatly boost future gene therapy and scientific research.

Ubiquitin-induced RNF168 condensation promotes DNA double-strand break repair.

Rapid accumulation of repair factors at DNA double-strand breaks (DSBs) is essential for DSB repair. Several factors involved in DSB repair have been found undergoing liquid-liquid phase separation (LLPS) at DSB sites to facilitate DNA repair. RNF168, a RING-type E3 ubiquitin ligase, catalyzes H2A.X ubiquitination for recruiting DNA repair factors. Yet, whether RNF168 undergoes LLPS at DSB sites remains unclear. Here, we identified K63-linked polyubiquitin-triggered RNF168 condensation which further promoted RNF168-mediated DSB repair. RNF168 formed liquid-like condensates upon irradiation in the nucleus while purified RNF168 protein also condensed in vitro. An intrinsically disordered region containing amino acids 460-550 was identified as the essential domain for RNF168 condensation. Interestingly, LLPS of RNF168 was significantly enhanced by K63-linked polyubiquitin chains, and LLPS largely enhanced the RNF168-mediated H2A.X ubiquitination, suggesting a positive feedback loop to facilitate RNF168 rapid accumulation and its catalytic activity. Functionally, LLPS deficiency of RNF168 resulted in delayed recruitment of 53BP1 and BRCA1 and subsequent impairment in DSB repair. Taken together, our finding demonstrates the pivotal effect of LLPS in RNF168-mediated DSB repair.

Genome-wide analysis of copy-number variation in humans with cleft lip and/or cleft palate identifies COBLL1, RIC1, and ARHGEF38 as clefting genes.

Cleft lip with or without cleft palate (CL/P) is a common birth defect with a complex, heterogeneous etiology. It is well established that common and rare sequence variants contribute to the formation of CL/P, but the contribution of copy-number variants (CNVs) to cleft formation remains relatively understudied. To fill this knowledge gap, we conducted a large-scale comparative analysis of genome-wide CNV profiles of 869 individuals from the Philippines and 233 individuals of European ancestry with CL/P with three primary goals: first, to evaluate whether differences in CNV number, amount of genomic content, or amount of coding genomic content existed within clefting subtypes; second, to assess whether CNVs in our cohort overlapped with known Mendelian clefting loci; and third, to identify unestablished Mendelian clefting genes. Significant differences in CNVs across cleft types or in individuals with non-syndromic versus syndromic clefts were not observed; however, several CNVs in our cohort overlapped with known syndromic and non-syndromic Mendelian clefting loci. Moreover, employing a filtering strategy relying on population genetics data that rare variants are on the whole more deleterious than common variants, we identify several CNV-associated gene losses likely driving non-syndromic clefting phenotypes. By prioritizing genes deleted at a rare frequency across multiple individuals with clefts yet enriched in our cohort of individuals with clefts compared to control subjects, we identify COBLL1, RIC1, and ARHGEF38 as clefting genes. CRISPR-Cas9 mutagenesis of these genes in Xenopus laevis and Danio rerio yielded craniofacial dysmorphologies, including clefts analogous to those seen in human clefting disorders.

Mechanistic study of acupuncture on the pterygopalatine ganglion to improve allergic rhinitis: analysis of multi-target effects based on bioinformatics/network topology strategie.

One of the prevalent chronic inflammatory disorders of the nasal mucosa, allergic rhinitis (AR) has become more widespread in recent years. Acupuncture pterygopalatine ganglion (aPPG) is an emerging alternative therapy that is used to treat AR, but the molecular mechanisms underlying its anti-inflammatory effects are unclear. This work methodically demonstrated the multi-target mechanisms of aPPG in treating AR based on bioinformatics/topology using techniques including text mining, bioinformatics, and network topology, among others. A total of 16 active biomarkers and 108 protein targets related to aPPG treatment of AR were obtained. A total of 345 Gene Ontology terms related to aPPG of AR were identified, and 135 pathways were screened based on Kyoto Encyclopedia of Genes and Genomes analysis. Our study revealed for the first time the multi-targeted mechanism of action of aPPG in the treatment of AR. In animal experiments, aPPG ameliorated rhinitis symptoms in OVA-induced AR rats; decreased serum immunoglobulin E, OVA-sIgE, and substance P levels; elevated serum neuropeptide Y levels; and modulated serum Th1/Th2/Treg/Th17 cytokine expression by a mechanism that may be related to the inhibition of activation of the TLR4/NF-κB/NLRP3 signaling pathway. In vivo animal experiments once again validated the results of the bioinformatics analysis. This study revealed a possible multi-target mechanism of action between aPPG and AR, provided new insights into the potential pathogenesis of AR, and proved that aPPG was a promising complementary alternative therapy for the treatment of AR.

Enhanced polygenic risk score incorporating gene-environment interaction suggests the association of major depressive disorder with cardiac and lung function.

BACKGROUND: Depression has been linked to an increased risk of cardiovascular and respiratory diseases; however, its impact on cardiac and lung function remains unclear, especially when accounting for potential gene-environment interactions. METHODS: We developed a novel polygenic and gene-environment interaction risk score (PGIRS) integrating the major genetic effect and gene-environment interaction effect of depression-associated loci. The single nucleotide polymorphisms (SNPs) demonstrating major genetic effect or environmental interaction effect were obtained from genome-wide SNP association and SNP-environment interaction analyses of depression. We then calculated the depression PGIRS for non-depressed individuals, using smoking and alcohol consumption as environmental factors. Using linear regression analysis, we assessed the associations of PGIRS and conventional polygenic risk score (PRS) with lung function (N = 42 886) and cardiac function (N = 1791) in the subjects with or without exposing to smoking and alcohol drinking. RESULTS: We detected significant associations of depression PGIRS with cardiac and lung function, contrary to conventional depression PRS. Among smokers, forced vital capacity exhibited a negative association with PGIRS (ß = -0.037, FDR = 1.00 × 10-8), contrasting with no significant association with PRS (ß = -0.002, FDR = 0.943). In drinkers, we observed a positive association between cardiac index with PGIRS (ß = 0.088, FDR = 0.010), whereas no such association was found with PRS (ß = 0.040, FDR = 0.265). Notably, in individuals who both smoked and drank, forced expiratory volume in 1-second demonstrated a negative association with PGIRS (ß = -0.042, FDR = 6.30 × 10-9), but not with PRS (ß = -0.003, FDR = 0.857). CONCLUSIONS: Our findings underscore the profound impact of depression on cardiac and lung function, highlighting the enhanced efficacy of considering gene-environment interactions in PRS-based studies.

Chronic active Epstein-Barr virus disease originates from infected hematopoietic stem cells.

ABSTRACT: Chronic active Epstein-Barr virus (EBV) disease (CAEBV) is a lethal syndrome because of persistent EBV infection. When diagnosed as CAEBV, EBV infection was observed in multiple hematopoietic lineages, but the etiology of CAEBV is still elusive. Bone marrow and peripheral cells derived from 5 patients with CAEBV, 1 patient with EBV-associated hemophagocytic lymphohistiocytosis, and 2 healthy controls were analyzed. Multiple assays were applied to identify and characterize EBV-infected cells, including quantitative polymerase chain reaction, PrimeFlow, and single-cell RNA-sequencing (scRNA-seq). Based on scRNA-seq data, alterations in gene expression of particular cell types were analyzed between patients with CAEBV and controls, and between infected and uninfected cells. One patient with CAEBV was treated with allogeneic hematopoietic stem cell transplantation (HSCT), and the samples derived from this patient were analyzed again 6 months after HSCT. EBV infected the full spectrum of the hematopoietic system including both lymphoid and myeloid lineages, as well as the hematopoietic stem cells (HSCs) of the patients with CAEBV. EBV-infected HSCs exhibited a higher differentiation rate toward downstream lineages, and the EBV infection had an impact on both the innate and adaptive immunity, resulting in inflammatory symptoms. EBV-infected cells were thoroughly removed from the hematopoietic system after HSCT. Taken together, multiple lines of evidence presented in this study suggest that CAEBV disease originates from the infected HSCs, which might potentially lead to innovative therapy strategies for CAEBV.

CropGS-Hub: a comprehensive database of genotype and phenotype resources for genomic prediction in major crops.

The explosive amount of multi-omics data has brought a paradigm shift both in academic research and further application in life science. However, managing and reusing the growing resources of genomic and phenotype data points presents considerable challenges for the research community. There is an urgent need for an integrated database that combines genome-wide association studies (GWAS) with genomic selection (GS). Here, we present CropGS-Hub, a comprehensive database comprising genotype, phenotype, and GWAS signals, as well as a one-stop platform with built-in algorithms for genomic prediction and crossing design. This database encompasses a comprehensive collection of over 224 billion genotype data and 434 thousand phenotype data generated from >30 000 individuals in 14 representative populations belonging to 7 major crop species. Moreover, the platform implemented three complete functional genomic selection related modules including phenotype prediction, user model training and crossing design, as well as a fast SNP genotyper plugin-in called SNPGT specifically built for CropGS-Hub, aiming to assist crop scientists and breeders without necessitating coding skills. CropGS-Hub can be accessed at https://iagr.genomics.cn/CropGS/.

Chromosomal-scale genomes of two Rosa species provide insights into genome evolution and ascorbate accumulation.

Rosa roxburghii and Rosa sterilis, two species belonging to the Rosaceae family, are widespread in the southwest of China. These species have gained recognition for their remarkable abundance of ascorbate in their fresh fruits, making them an ideal vitamin C resource. In this study, we generated two high-quality chromosome-scale genome assemblies for R. roxburghii and R. sterilis, with genome sizes of 504 and 981.2 Mb, respectively. Notably, we present a haplotype-resolved, chromosome-scale assembly for diploid R. sterilis. Our results indicated that R. sterilis originated from the hybridization of R. roxburghii and R. longicuspis. Genome analysis revealed the absence of recent whole-genome duplications in both species and identified a series of duplicated genes that possibly contributing to the accumulation of flavonoids. We identified two genes in the ascorbate synthesis pathway, GGP and GalLDH, that show signs of positive selection, along with high expression levels of GDP-d-mannose 3', 5'-epimerase (GME) and GDP-l-galactose phosphorylase (GGP) during fruit development. Furthermore, through co-expression network analysis, we identified key hub genes (MYB5 and bZIP) that likely regulate genes in the ascorbate synthesis pathway, promoting ascorbate biosynthesis. Additionally, we observed the expansion of terpene synthase genes in these two species and tissue expression patterns, suggesting their involvement in terpenoid biosynthesis. Our research provides valuable insights into genome evolution and the molecular basis of the high concentration of ascorbate in these two Rosa species.

Signatures of Adaptation and Purifying Selection in Highland Populations of Dasiphora fruticosa.

High mountains harbor a considerable proportion of biodiversity, but we know little about how diverse plants adapt to the harsh environment. Here we finished a high-quality genome assembly for Dasiphora fruticosa, an ecologically important plant distributed in the Qinghai-Tibetan Plateau and lowland of the Northern Hemisphere, and resequenced 592 natural individuals to address how this horticulture plant adapts to highland. Demographic analysis revealed D. fruticosa underwent a bottleneck after Naynayxungla Glaciation. Selective sweep analysis of two pairs of lowland and highland populations identified 63 shared genes related to cell wall organization or biogenesis, cellular component organization, and dwarfism, suggesting parallel adaptation to highland habitats. Most importantly, we found that stronger purging of estimated genetic load due to inbreeding in highland populations apparently contributed to their adaptation to the highest mountain. Our results revealed how plants could tolerate the extreme plateau, which could provide potential insights for species conservation and crop breeding.

WMSA 2: a multiple DNA/RNA sequence alignment tool implemented with accurate progressive mode and a fast win-win mode combining the center star and progressive strategies.

Multiple sequence alignment is widely used for sequence analysis, such as identifying important sites and phylogenetic analysis. Traditional methods, such as progressive alignment, are time-consuming. To address this issue, we introduce StarTree, a novel method to fast construct a guide tree by combining sequence clustering and hierarchical clustering. Furthermore, we develop a new heuristic similar region detection algorithm using the FM-index and apply the k-banded dynamic program to the profile alignment. We also introduce a win-win alignment algorithm that applies the central star strategy within the clusters to fast the alignment process, then uses the progressive strategy to align the central-aligned profiles, guaranteeing the final alignment's accuracy. We present WMSA 2 based on these improvements and compare the speed and accuracy with other popular methods. The results show that the guide tree made by the StarTree clustering method can lead to better accuracy than that of PartTree while consuming less time and memory than that of UPGMA and mBed methods on datasets with thousands of sequences. During the alignment of simulated data sets, WMSA 2 can consume less time and memory while ranking at the top of Q and TC scores. The WMSA 2 is still better at the time, and memory efficiency on the real datasets and ranks at the top on the average sum of pairs score. For the alignment of 1 million SARS-CoV-2 genomes, the win-win mode of WMSA 2 significantly decreased the consumption time than the former version. The source code and data are available at https://github.com/malabz/WMSA2.

Bioinformatics toolbox for exploring target mutation-induced drug resistance.

Drug resistance is increasingly among the main issues affecting human health and threatening agriculture and food security. In particular, developing approaches to overcome target mutation-induced drug resistance has long been an essential part of biological research. During the past decade, many bioinformatics tools have been developed to explore this type of drug resistance, and they have become popular for elucidating drug resistance mechanisms in a low cost, fast and effective way. However, these resources are scattered and underutilized, and their strengths and limitations have not been systematically analyzed and compared. Here, we systematically surveyed 59 freely available bioinformatics tools for exploring target mutation-induced drug resistance. We analyzed and summarized these resources based on their functionality, data volume, data source, operating principle, performance, etc. And we concisely discussed the strengths, limitations and application examples of these tools. Specifically, we tested some predictive tools and offered some thoughts from the clinician's perspective. Hopefully, this work will provide a useful toolbox for researchers working in the biomedical, pesticide, bioinformatics and pharmaceutical engineering fields, and a good platform for non-specialists to quickly understand drug resistance prediction.

FMAlign2: a novel fast multiple nucleotide sequence alignment method for ultralong datasets.

MOTIVATION: In bioinformatics, multiple sequence alignment (MSA) is a crucial task. However, conventional methods often struggle with aligning ultralong sequences. To address this issue, researchers have designed MSA methods rooted in a vertical division strategy, which segments sequence data for parallel alignment. A prime example of this approach is FMAlign, which utilizes the FM-index to extract common seeds and segment the sequences accordingly. RESULTS: FMAlign2 leverages the suffix array to identify maximal exact matches, redefining the approach of FMAlign from searching for global chains to partial chains. By using a vertical division strategy, large-scale problem is deconstructed into manageable tasks, enabling parallel execution of subMSA. Furthermore, sequence-profile alignment and refinement are incorporated to concatenate subsets, yielding the final result seamlessly. Compared to FMAlign, FMAlign2 markedly augments the segmentation of sequences and significantly reduces the time while maintaining accuracy, especially on ultralong datasets. Importantly, FMAlign2 enhances existing MSA methods by conferring the capability to handle sequences reaching billions in length within an acceptable time frame. AVAILABILITY AND IMPLEMENTATION: Source code and datasets are available at https://github.com/malabz/FMAlign2 and https://zenodo.org/records/10435770.

An integrated multi-omics analysis of identifies distinct molecular characteristics in pulmonary infections of Pseudomonas aeruginosa.

Pseudomonas aeruginosa (P. aeruginosa) can cause severe acute infections, including pneumonia and sepsis, and cause chronic infections, commonly in patients with structural respiratory diseases. However, the molecular and pathophysiological mechanisms of P. aeruginosa respiratory infection are largely unknown. Here, we performed assays for transposase-accessible chromatin using sequencing (ATAC-seq), transcriptomics, and quantitative mass spectrometry-based proteomics and ubiquitin-proteomics in P. aeruginosa-infected lung tissues for multi-omics analysis, while ATAC-seq and transcriptomics were also examined in P. aeruginosa-infected mouse macrophages. To identify the pivotal factors that are involved in host immune defense, we integrated chromatin accessibility and gene expression to investigate molecular changes in P. aeruginosa-infected lung tissues combined with proteomics and ubiquitin-proteomics. Our multi-omics investigation discovered a significant concordance for innate immunological and inflammatory responses following P. aeruginosa infection between hosts and alveolar macrophages. Furthermore, we discovered that multi-omics changes in pioneer factors Stat1 and Stat3 play a crucial role in the immunological regulation of P. aeruginosa infection and that their downstream molecules (e.g., Fas) may be implicated in both immunosuppressive and inflammation-promoting processes. Taken together, these findings indicate that transcription factors and their downstream signaling molecules play a critical role in the mobilization and rebalancing of the host immune response against P. aeruginosa infection and may serve as potential targets for bacterial infections and inflammatory diseases, providing insights and resources for omics analyses.

Cell wall invertase 3 plays critical roles in providing sugars during pollination and fertilization in cucumber.

In plants, pollen-pistil interactions during pollination and fertilization mediate pollen hydration and germination, pollen tube growth, and seed set and development. Cell wall invertases (CWINs) help provide the carbohydrates for pollen development; however, their roles in pollination and fertilization have not been well established. In cucumber (Cucumis sativus), CsCWIN3 showed the highest expression in flowers, and we further examined CsCWIN3 for functions during pollination to seed set. Both CsCWIN3 transcript and CsCWIN3 protein exhibited similar expression patterns in the sepals, petals, stamen filaments, anther tapetum, and pollen of male flowers, as well as in the stigma, style, transmitting tract, and ovule funiculus of female flowers. Notably, repression of CsCWIN3 in cucumber did not affect the formation of parthenocarpic fruit but resulted in an arrested growth of stigma integuments in female flowers and a partially delayed dehiscence of anthers with decreased pollen viability in male flowers. Consequently, the pollen tube grew poorly in the gynoecia after pollination. In addition, CsCWIN3-RNA interference plants also showed affected seed development. Considering that sugar transporters could function in cucumber fecundity, we highlight the role of CsCWIN3 and a potential close collaboration between CWIN and sugar transporters in these processes. Overall, we used molecular and physiological analyses to determine the CsCWIN3-mediated metabolism during pollen formation, pollen tube growth, and plant fecundity. CsCWIN3 has essential roles from pollination and fertilization to seed set but not parthenocarpic fruit development in cucumber.

A genome-wide association study identifies candidate genes for sleep disturbances in depressed individuals.

OBJECTIVE: This study aimed to identify candidate loci and genes related to sleep disturbances in depressed individuals and clarify the co-occurrence of sleep disturbances and depression from the genetic perspective. METHODS: The study subjects (including 58,256 self-reported depressed individuals and 6,576 participants with PHQ-9 score ≥ 10, respectively) were collected from the UK Biobank, which were determined based on the Patient Health Questionnaire (PHQ-9) and self-reported depression status, respectively. Sleep related traits included chronotype, insomnia, snoring and daytime dozing. Genome-wide association studies (GWASs) of sleep related traits in depressed individuals were conducted by PLINK 2.0 adjusting age, sex, Townsend deprivation index and 10 principal components as covariates. The CAUSALdb database was used to explore the mental traits associated with the candidate genes identified by the GWAS. RESULTS: GWAS detected 15 loci significantly associated with chronotype in the subjects with self-reported depression, such as rs12736689 at RNASEL (P = 1.00 × 10- 09), rs509476 at RGS16 (P = 1.58 × 10- 09) and rs1006751 at RFX4 (P = 1.54 × 10- 08). 9 candidate loci were identified in the subjects with PHQ-9 ≥ 10, of which 2 loci were associated with insomnia such as rs115379847 at EVC2 (P = 3.50 × 10- 08), and 7 loci were associated with daytime dozing, such as rs140876133 at SMYD3 (P = 3.88 × 10- 08) and rs139156969 at ROBO2 (P = 3.58 × 10- 08). Multiple identified genes, such as RNASEL, RGS16, RFX4 and ROBO2 were reported to be associated with chronotype, depression or cognition in previous studies. CONCLUSION: Our study identified several candidate genes related to sleep disturbances in depressed individuals, which provided new clues for understanding the biological mechanism underlying the co-occurrence of depression and sleep disorders.

Xbp1 promotes odontoblastic differentiation through modulating mitochondrial homeostasis.

Odontoblast differentiation depends on the orderly recruitment of transcriptional factors (TFs) in the transcriptional regulatory network. The depletion of crucial TFs disturbs dynamic alteration of the chromatin landscape and gene expression profile, leading to developmental defects. Our previous studies have revealed that the basic leucine zipper (bZIP) TF family is crucial in odontoblastic differentiation, but the function of bZIP TF family member XBP1 is still unknown. Here, we showed the stage-specific expression patterns of the spliced form Xbp1s during tooth development. Elevated Xbp1 expression and nuclear translocation of XBP1S in mesenchymal stem cells (MSCs) were induced by differentiation medium in vitro. Diminution of Xbp1 expression impaired the odontogenic differentiation potential of MSCs. The further integration of ATAC-seq and RNA-seq identified Hspa9 as a direct downstream target, an essential mitochondrial chaperonin gene that modulated mitochondrial homeostasis. The amelioration of mitochondrial dysfunction rescued the impaired odontogenic differentiation potential of MSCs caused by the diminution of Xbp1. Furthermore, the overexpression of Hspa9 rescued Xbp1-deficient defects in odontoblastic differentiation. Our study illustrates the crucial role of Xbp1 in odontoblastic differentiation via modulating mitochondrial homeostasis and brings evidence to the therapy of mitochondrial diseases caused by genetic defects.