A large number of SARS-related coronaviruses (SARSr-CoV) have been detected in horseshoe bats since 2005 in different areas of China. However, these bat SARSr-CoVs show sequence differences from SARS coronavirus (SARS-CoV) in different genes (S, ORF8, ORF3, etc) and are considered unlikely to represent the direct progenitor of SARS-CoV. Herein, we report the findings of our 5-year surveillance of SARSr-CoVs in a cave inhabited by multiple species of horseshoe bats in Yunnan Province, China. The full-length genomes of 11 newly discovered SARSr-CoV strains, together with our previous findings, reveals that the SARSr-CoVs circulating in this single location are highly diverse in the S gene, ORF3 and ORF8. Importantly, strains with high genetic similarity to SARS-CoV in the hypervariable N-terminal domain (NTD) and receptor-binding domain (RBD) of the S1 gene, the ORF3 and ORF8 region, respectively, were all discovered in this cave. In addition, we report the first discovery of bat SARSr-CoVs highly similar to human SARS-CoV in ORF3b and in the split ORF8a and 8b. Moreover, SARSr-CoV strains from this cave were more closely related to SARS-CoV in the non-structural protein genes ORF1a and 1b compared with those detected elsewhere. Recombination analysis shows evidence of frequent recombination events within the S gene and around the ORF8 between these SARSr-CoVs. We hypothesize that the direct progenitor of SARS-CoV may have originated after sequential recombination events between the precursors of these SARSr-CoVs. Cell entry studies demonstrated that three newly identified SARSr-CoVs with different S protein sequences are all able to use human ACE2 as the receptor, further exhibiting the close relationship between strains in this cave and SARS-CoV. This work provides new insights into the origin and evolution of SARS-CoV and highlights the necessity of preparedness for future emergence of SARS-like diseases.

Transcriptional dynamic in response to environmental and developmental cues are fundamental to biology, yet many mechanistic aspects are poorly understood. One such example is fungal plant pathogens, which use secreted proteins and small molecules, termed effectors, to suppress host immunity and promote colonization. Effectors are highly expressed in planta but remain transcriptionally repressed ex planta, but our mechanistic understanding of these transcriptional dynamics remains limited. We tested the hypothesis that repressive histone modification at H3-Lys27 underlies transcriptional silencing ex planta, and that exchange for an active chemical modification contributes to transcription of in planta induced genes. Using genetics, chromatin immunoprecipitation and sequencing and RNA-sequencing, we determined that H3K27me3 provides significant local transcriptional repression. We detail how regions that lose H3K27me3 gain H3K27ac, and these changes are associated with increased transcription. Importantly, we observed that many in planta induced genes were marked by H3K27me3 during axenic growth, and detail how altered H3K27 modification influences transcription. ChIP-qPCR during in planta growth suggests that H3K27 modifications are generally stable, but can undergo dynamics at specific genomic locations. Our results support the hypothesis that dynamic histone modifications at H3K27 contributes to fungal genome regulation and specifically contributes to regulation of genes important during host infection.

Glioblastoma (GBM) ranks among the most lethal of human cancers, containing glioma stem cells (GSCs) that display therapeutic resistance. Here, we report that the lncRNA INHEG is highly expressed in GSCs compared to differentiated glioma cells (DGCs) and promotes GSC self-renewal and tumorigenicity through control of rRNA 2'-O-methylation. INHEG induces the interaction between SUMO2 E3 ligase TAF15 and NOP58, a core component of snoRNP that guides rRNA methylation, to regulate NOP58 sumoylation and accelerate the C/D box snoRNP assembly. INHEG activation enhances rRNA 2'-O-methylation, thereby increasing the expression of oncogenic proteins including EGFR, IGF1R, CDK6 and PDGFRB in glioma cells. Taken together, this study identifies a lncRNA that connects snoRNP-guided rRNA 2'-O-methylation to upregulated protein translation in GSCs, supporting an axis for potential therapeutic targeting of gliomas.

Hand, foot, and mouth disease (HFMD) is an infectious disease that mainly affects infants and children, causing considerable morbidity and mortality worldwide. HFMD is commonly caused by enterovirus 71 (EV71) and coxsackieviruses A16 (CVA16), A6 (CVA6), and A10 (CVA10). Formalin-inactivated EV71 vaccines are currently available in China; however, these vaccines fail to confer cross-protection against infections by other HFMD-causing enteroviruses, highlighting the necessity of developing a multivalent HFMD vaccine. Our previous studies demonstrated that recombinant virus-like particles (VLP) of EV71, CVA16, and CVA6 are capable of inducing protective immunity against homologous virus challenges in mice. In this study, we generated CVA10-VLP using a baculovirus-insect cell expression system and then combined CVA10-VLP with EV71-VLP, CVA16-VLP, and CVA6-VLP to formulate a tetravalent VLP vaccine. Immunogenicity and protective efficacy of tetravalent VLP vaccine was compared with that of monovalent VLP vaccines. Mouse immunization studies revealed that the tetravalent vaccine elicited antigen-specific and long-lasting serum antibody responses comparable to those elicited by its corresponding monovalent vaccines. Moreover, tetravalent vaccine immune sera strongly neutralized EV71, CVA16, CVA10, and CVA6 strains with neutralization titers similar to those of their monovalent counterparts, indicating a good compatibility among the four antigens in the combination vaccine. Importantly, passively transferred tetravalent vaccine-immunized sera conferred efficient protection against single or mixed infections with EV71, CVA16, CVA10, and CVA6 viruses in mice, whereas the monovalent vaccines could only protect mice against homotypic virus infections but not heterotypic challenges. These results demonstrate that the tetravalent VLP vaccine represents a promising broad-spectrum HFMD vaccine candidate.

Genetic leukoencephalopathies (gLEs) are a group of heterogeneous disorders with white matter abnormalities affecting the central nervous system (CNS). The causative mutation in ~50% of gLEs is unknown. Using whole exome sequencing (WES), we identified homozygosity for a missense variant, VPS11: c.2536T>G (p.C846G), as the genetic cause of a leukoencephalopathy syndrome in five individuals from three unrelated Ashkenazi Jewish (AJ) families. All five patients exhibited highly concordant disease progression characterized by infantile onset leukoencephalopathy with brain white matter abnormalities, severe motor impairment, cortical blindness, intellectual disability, and seizures. The carrier frequency of the VPS11: c.2536T>G variant is 1:250 in the AJ population (n = 2,026). VPS11 protein is a core component of HOPS (homotypic fusion and protein sorting) and CORVET (class C core vacuole/endosome tethering) protein complexes involved in membrane trafficking and fusion of the lysosomes and endosomes. The cysteine 846 resides in an evolutionarily conserved cysteine-rich RING-H2 domain in carboxyl terminal regions of VPS11 proteins. Our data shows that the C846G mutation causes aberrant ubiquitination and accelerated turnover of VPS11 protein as well as compromised VPS11-VPS18 complex assembly, suggesting a loss of function in the mutant protein. Reduced VPS11 expression leads to an impaired autophagic activity in human cells. Importantly, zebrafish harboring a vps11 mutation with truncated RING-H2 domain demonstrated a significant reduction in CNS myelination following extensive neuronal death in the hindbrain and midbrain. Thus, our study reveals a defect in VPS11 as the underlying etiology for an autosomal recessive leukoencephalopathy disorder associated with a dysfunctional autophagy-lysosome trafficking pathway.

The Raf kinase inhibitory protein (RKIP) is down-regulated in multiple types of human cancers. Decreased RKIP transcription activity may be one of the major mechanisms responsible for the downregulation of RKIP expression in human diseases. To test this hypothesis, we need to gain basic knowledge of the transcriptional regulation of RKIP. To achieve this objective, we made a systematic effort to identify cis-acting elements and trans-acting factors that control RKIP promoter activity. We found that full RKIP promoter activity requires the region -56 to +261 relative to the transcription start site. Within the full promoter region, there are two motifs rich in G/C that responded to transcription factor Sp1, one cAMP-responsive element that responded to the transcription factor CREB, and one docking site for the histone acetylase p300. In human melanoma A375 cells and human cervical cancer HeLa cells, mutation or deletion of each of these cis-acting elements decreased promoter activity. In A375 cells, knockdown of the corresponding transcription factors Sp1, CREB, or p300 decreased RKIP promoter activity, whereas overexpression of CREB and p300 increased RKIP promoter activity. The results obtained with HeLa cells also supported the idea that Sp1 and CREB play positive roles in the regulation of RKIP transcription. These findings suggest that regulators of the expression or activity of Sp1, CREB, and p300 are involved in regulating RKIP transcription.

The Notch2 is a critical membrane receptor for B-cell functions, and also displays various biological roles in lymphoma pathogenesis. In this article, we reported that 3 of 69 (4.3%) diffuse large B-cell lymphomas (DLBCLs) exhibited a truncate NOTCH2 mutation at the nucleotide 7605 (G/A) in the cDNA sequence, which led to partial deletion of the C-terminal of PEST (proline-, glutamic acid-, serine- and threonine-rich) domain. The truncate Notch2 activated both the Notch2 and the NF-κB signals and promoted the proliferation of B-cell lymphoma cell lines, including DLBCL and Burkitt's lymphoma cell lines. Moreover, the ectopic proliferation was completely inhibited by ammonium pyrrolidinedithiocarbamate (PDTC), an NF-κB inhibitor. Simultaneously, PDTC also reduced the expression level of Notch2. Based on these results, we conclude that the Notch2 receptor with PEST domain truncation enhances cell proliferation which may be associated with the activation of the Notch2 and the NF-κB signaling. Our results are expected to provide a possible target for new DLBCL therapies by suppressing the Notch2 and the NF-κB signaling.

The heterotrimeric G-protein mediates growth and development by perceiving and transmitting signals in multiple organisms. Alternative splicing (AS), a vital process for regulating gene expression at the post-transcriptional level, plays a significant role in plant adaptation and evolution. Here, we identified five splicing variants of Gγ subunit gene TaGS3 (TaGS3.1 to TaGS3.5), which showed expression divergence during wheat polyploidization, and differential function in grain weight and size determination. TaGS3.1 overexpression significantly reduced grain weight by 5.89% and grain length by 5.04%, while TaGS3.2-3.4 overexpression did not significantly alter grain size compared to wild type. Overexpressing TaGS3.5 significantly increased the grain weight by 5.70% and grain length by 4.30%. Biochemical assays revealed that TaGS3 isoforms (TaGS3.1-3.4) with an intact OSR domain interact with WGB1 to form active Gβγ heterodimers that further interact with WGA1 to form inactive Gαβγ heterotrimers. Truncated isoforms TaGS3.2-3.4 , which lack the C-terminal Cys-rich region but have enhanced binding affinity to WGB1, antagonistically compete with TaGS3.1 to bind WGB1, while TaGS3.5 with an incomplete OSR domain does not interact with WGB1. Taking these observations together, we proposed that TaGS3 differentially regulates grain size via AS, providing a strategy by which the grain size is fine-tuned and regulated at the post-transcriptional level.

Purple corn is a rich source of anthocyanins. In the experiment, two anthocyanins-enriched purple corn lines Ha0414 and Ha6130 were identified. The anthocyanins were respectively accumulated in the pericarp of Ha0414 and the aleurone layer of Ha6130 with different composition and content. Transcriptome analysis of the two tissues in both lines identified 16 and 14 differentially expressed genes belonging to anthocyanin metabolism pathway in pericarp and the aleurone layer, individually. Of these genes, two genes encoding 2-oxoglutarate (2OG) and Fe (II)-dependent oxygenase superfamily proteins, and one gene annotated as UDP-glycosyltransferase superfamily protein exhibited increased transcript abundance in both the colored pericarp and aleurone layer. Otherwise, one gene annotated as flavonoid 3', 5'-hydroxylase, and another gene encoding flavonoid 3'-monooxygenase displayed increased transcript abundance in the aleurone layer of Ha6130. Moreover, 36 transcription factors were identified with increased transcript abundance in the pericarp of Ha0414, such as bHLH transcription factors, WRKY transcription factors, and HB transcription factors. And 79 transcription factors were isolated with an increased expression level in the aleurone layer of Ha6130, including MYB transcription factors, MYB-related transcription factors, and bHLH transcription factors. These genes expression may result in the tissue-specific accumulation of anthocyanins in pericarp and aleurone layer.

Anthocyanins are natural pigments that have important functions in plant growth and development. Radish taproots are rich in anthocyanins which confer different taproot colors and are potentially beneficial to human health. The crop differentially accumulates anthocyanin during various stages of growth, yet molecular mechanisms underlying this differential anthocyanin accumulation remains unknown. In the present study, transcriptome analysis was used to concisely identify putative genes involved in anthocyanin biosynthesis in radish. Spatial-temporal transcript expressions were then profiled in four color variant radish cultivars. From the total transcript sequences obtained through illumina sequencing, 102 assembled unigenes, and 20 candidate genes were identified to be involved in anthocyanin biosynthesis. Fifteen genomic sequences were isolated and sequenced from radish taproot. The length of these sequences was between 900 and 1,579 bp, and the unigene coverage to all of the corresponding cloned sequences was more than 93%. Gene structure analysis revealed that RsF3'H is intronless and anthocyanin biosynthesis genes (ABGs) bear asymmetrical exons, except RsSAM. Anthocyanin accumulation showed a gradual increase in the leaf of the red radish and the taproot of colored cultivars during development, with a rapid increase at 30 days after sowing (DAS), and the highest content at maturity. Spatial-temporal transcriptional analysis of 14 genes revealed detectable expressions of 12 ABGs in various tissues at different growth levels. The investigation of anthocyanin accumulation and gene expression in four color variant radish cultivars, at different stages of development, indicated that total anthocyanin correlated with transcript levels of ABGs, particularly RsUFGT, RsF3H, RsANS, RsCHS3 and RsF3'H1. Our results suggest that these candidate genes play key roles in phenotypic and spatial-temporal anthocyanin accumulation in radish through coordinated regulation and the major control point in anthocyanin biosynthesis in radish is RsUFGT. The present findings lend invaluable insights into anthocyanin biosynthesis and may facilitate genetic manipulation for enhanced anthocyanin content in radish.

Gastric cancer (GC) is a highly heterogeneous disease. To identify potential clinically actionable therapeutic targets that may inform individualized treatment strategies, we performed whole-exome sequencing on 78 GCs of differing histologies and anatomic locations, as well as whole-genome sequencing on two GC cases, each with three primary tumors and two matching lymph node metastases. The data showed two distinct GC subtypes with either high-clonality (HiC) or low-clonality (LoC). The HiC subtype of intratumoral heterogeneity was associated with older age, TP53 (tumor protein P53) mutation, enriched C > G transition, and significantly shorter survival, whereas the LoC subtype was associated with younger age, ARID1A (AT rich interactive domain 1A) mutation, and significantly longer survival. Phylogenetic tree analysis of whole-genome sequencing data from multiple samples of two patients supported the clonal evolution of GC metastasis and revealed the accumulation of genetic defects that necessitate combination therapeutics. The most recurrently mutated genes, which were validated in a separate cohort of 216 cases by targeted sequencing, were members of the homologous recombination DNA repair, Wnt, and PI3K-ERBB pathways. Notably, the drugable NRG1 (neuregulin-1) and ERBB4 (V-Erb-B2 avian erythroblastic leukemia viral oncogene homolog 4) ligand-receptor pair were mutated in 10% of GC cases. Mutations of the BRCA2 (breast cancer 2, early onset) gene, found in 8% of our cohort and validated in The Cancer Genome Atlas GC cohort, were associated with significantly longer survivals. These data define distinct clinicogenetic forms of GC in the Chinese population that are characterized by specific mutation sets that can be investigated for efficacy of single and combination therapies.

Non-invasively evaluating gene expression products in human pre-implantation embryos remains a significant challenge. Here, we develop a non-invasive method for comprehensive characterization of the extracellular RNAs (exRNAs) in a single droplet of spent media that was used to culture human in vitro fertilization embryos. We generate the temporal extracellular transcriptome atlas (TETA) of human pre-implantation development. TETA consists of 245 exRNA sequencing datasets for five developmental stages. These data reveal approximately 4,000 exRNAs at each stage. The exRNAs of the developmentally arrested embryos are enriched with the genes involved in negative regulation of the cell cycle, revealing an exRNA signature of developmental arrest. Furthermore, a machine-learning model can approximate the morphology-based rating of embryo quality based on the exRNA levels. These data reveal the widespread presence of coding gene-derived exRNAs at every stage of human pre-implantation development, and these exRNAs provide rich information on the physiology of the embryo.

IL-2 plays a key role in the survival and proliferation of immune cells, especially T lymphocytes. Its expression is precisely regulated at transcriptional and posttranscriptional level. IL-2 is known to be regulated by RNA binding proteins, such as tristetraprolin (TTP), via an AU-rich element (ARE) in the 3'-untranslated region (3'UTR) to influence the stability of mRNA. MCPIP1, identified as a novel RNase, can degrade IL-6, IL-12 and TNF-α mRNA by an ARE-independent pathway in the activation of macrophages. Here, we reported that MCPIP1 was induced in the activation of T lymphocytes and negatively regulated IL-2 gene expression in both mouse and human primary T lymphocytes through destabilizing its mRNA. A set of Luciferase reporter assay demonstrated that a non-ARE conserved element in IL-2 3'UTR, which formed a stem-loop structure, responded to MCPIP1 activity.RNA immunoprecipitation and Biotin pulldown experiments further suggested that MCPIP1 could modestly bind to IL-2 mRNA. Taken together, these data demonstrate that MCPIP1 down-regulates IL-2 via an ARE-independent pathway.

Host response diffuse large B-cell lymphoma (HR DLBCL) shares features of histologically defined T-cell/histiocyte-rich B-cell lymphoma, including fewer genetic abnormalities, frequent splenic and bone marrow involvement, and younger age at presentation. HR DLBCL is inherently less responsive to the standard treatment for DLBCL. Moreover, the mechanism of infiltration of HR DLBCL with preexisting abundant T-cells and dendritic cells is unknown, and their associated underlying immune responses incompletely defined. Here, hub genes and pathogenesis associated with HR DLBCL were explored to reveal molecular mechanisms and treatment targets.

ALG-2, a 22-kDa penta-EF-hand protein, is involved in cell death, signal transduction, membrane trafficking, etc., by interacting with various proteins in mammalian cells in a Ca2+-dependent manner. Most known ALG-2-interacting proteins contain proline-rich regions in which either PPYPXnYP (type 1 motif) or PXPGF (type 2 motif) is commonly found. Previous X-ray crystal structural analysis of the complex between ALG-2 and an ALIX peptide revealed that the peptide binds to the two hydrophobic pockets. In the present study, we resolved the crystal structure of the complex between ALG-2 and a peptide of Sec31A (outer shell component of coat complex II, COPII; containing the type 2 motif) and found that the peptide binds to the third hydrophobic pocket (Pocket 3). While amino acid substitution of Phe85, a Pocket 3 residue, with Ala abrogated the interaction with Sec31A, it did not affect the interaction with ALIX. On the other hand, amino acid substitution of Tyr180, a Pocket 1 residue, with Ala caused loss of binding to ALIX, but maintained binding to Sec31A. We conclude that ALG-2 recognizes two types of motifs at different hydrophobic surfaces. Furthermore, based on the results of serial mutational analysis of the ALG-2-binding sites in Sec31A, the type 2 motif was newly defined.

The exceptionally rich fossil record available for the equid family has provided textbook examples of macroevolutionary changes. Horses, asses, and zebras represent three extant subgenera of Equus lineage, while the Sussemionus subgenus is another remarkable Equus lineage ranging from North America to Ethiopia in the Pleistocene. We sequenced 26 archaeological specimens from Northern China in the Holocene that could be assigned morphologically and genetically to Equus ovodovi, a species representative of Sussemionus. We present the first high-quality complete genome of the Sussemionus lineage, which was sequenced to 13.4× depth of coverage. Radiocarbon dating demonstrates that this lineage survived until ~3500 years ago, despite continued demographic collapse during the Last Glacial Maximum and the great human expansion in East Asia. We also confirmed the Equus phylogenetic tree and found that Sussemionus diverged from the ancestor of non-caballine equids ~2.3-2.7 million years ago and possibly remained affected by secondary gene flow post-divergence. We found that the small genetic diversity, rather than enhanced inbreeding, limited the species' chances of survival. Our work adds to the growing literature illustrating how ancient DNA can inform on extinction dynamics and the long-term resilience of species surviving in cryptic population pockets.

The axon initial segment (AIS) is a specialized neuronal compartment required for action potential generation and neuronal polarity. However, understanding the mechanisms regulating AIS structure and function has been hindered by an incomplete knowledge of its molecular composition. Here, using immuno-proximity biotinylation we further define the AIS proteome and its dynamic changes during neuronal maturation. Among the many AIS proteins identified, we show that SCRIB is highly enriched in the AIS both in vitro and in vivo, and exhibits a periodic architecture like the axonal spectrin-based cytoskeleton. We find that ankyrinG interacts with and recruits SCRIB to the AIS. However, loss of SCRIB has no effect on ankyrinG. This powerful and flexible approach further defines the AIS proteome and provides a rich resource to elucidate the mechanisms regulating AIS structure and function.

Histone acetyltransferases (HATs) play an important role in eukaryotic transcription. Eight HATs identified in rice (OsHATs) can be organized into four families, namely the CBP (OsHAC701, OsHAC703, and OsHAC704), TAFII250 (OsHAF701), GNAT (OsHAG702, OsHAG703, and OsHAG704), and MYST (OsHAM701) families. The biological functions of HATs in rice remain unknown, so a comprehensive protein sequence analysis of the HAT families was conducted to investigate their potential functions. In addition, the subcellular localization and expression patterns of the eight OsHATs were analyzed.

The camellia weevil (CW [Curculio chinensis]) is a notorious host-specific predator of the seeds of Camellia species in China, causing seed losses of up to 60%. The weevil is capable of overcoming host tree chemical defenses, while the mechanisms of how these beetles contend with the toxic compounds are still unknown. Here, we examined the interaction between the gut microbes of CW and camellia seed chemistry and found that beetle-associated bacterial symbionts mediate tea saponin degradation. We demonstrate that the gut microbial community profile of CW was significantly plant associated, and the gut bacterial community associated with CW feeding on Camellia oleifera seeds is enriched with genes involved in tea saponin degradation compared with those feeding on Camellia sinensis and Camellia reticulata seeds. Twenty-seven bacteria from the genera Enterobacter, Serratia, Acinetobacter, and Micrococcus subsisted on tea saponin as a sole source of carbon and nitrogen, and Acinetobacter species are identified as being involved in the degradation of tea saponin. Our results provide the first metagenome of gut bacterial communities associated with a specialist insect pest of Camellia trees, and the results are consistent with a potential microbial contribution to the detoxification of tree-defensive chemicals.IMPORTANCE The gut microbiome may play an important role in insect-plant interactions mediated by plant secondary metabolites, but the microbial communities and functions of toxic plant feeders are still poorly characterized. In the present study, we provide the first metagenome of gut bacterial communities associated with a specialist weevil feeding on saponin-rich and saponin-low camellia seeds, and the results reveal the correlation between bacterial diversity and plant allelochemicals. We also used cultured microbes to establish their saponin-degradative capacity outside the insect. Our results provide new experimental context to better understand how gut microbial communities are influenced by plant secondary metabolites and how the resistance mechanisms involving microbes have evolved to deal with the chemical defenses of plants.

Psilocybin is a classic psychedelics, which has been shown to have antidepressant effects by many studies in recent years. In this study, we aim to evaluate the efficacy, acceptability and tolerability of psilocybin in the treatment of primary (major depressive disorder) or secondary (experiencing distress related to life-threatening diagnoses and terminal illness) depression.