Yuzhe Sun

Neuroligin-3 (Nlgn3) is an autism-associated cell-adhesion molecule that interacts with neurexins and is robustly expressed in both neurons and astrocytes. Neuronal Nlgn3 is an essential regulator of synaptic transmission but the function of astrocytic Nlgn3 is largely unknown. Given the high penetrance of Nlgn3 mutations in autism and the emerging role of astrocytes in neuropsychiatric disorders, we here asked whether astrocytic Nlgn3 might shape neural circuit properties in the cerebellum similar to neuronal Nlgn3. Imaging of tagged Nlgn3 protein produced by CRISPR/Cas9-mediated genome editing showed that Nlgn3 is enriched in the cell body but not the fine processes of cerebellar astrocytes (Bergmann glia). Astrocyte-specific knockout of Nlgn3 did not detectably alter the number of synapses, synaptic transmission, or astrocyte morphology in mouse cerebellum. However, spatial transcriptomic analyses revealed a significant shift in gene expression among multiple cerebellar cell types after the deletion of astrocytic Nlgn3. Hence, in contrast to neuronal Nlgn3, astrocytic Nlgn3 in the cerebellum is not involved in shaping synapses but may modulate gene expression in specific brain areas.

Coronavirus 2019 (COVID-19) is a complex disease that affects billions of people worldwide. Currently, effective etiological treatment of COVID-19 is still lacking; COVID-19 also causes damages to various organs that affects therapeutics and mortality of the patients. Surveillance of the treatment responses and organ injury assessment of COVID-19 patients are of high clinical value. In this study, we investigated the characteristic fragmentation patterns and explored the potential in tissue injury assessment of plasma cell-free DNA in COVID-19 patients. Through recruitment of 37 COVID-19 patients, 32 controls and analysis of 208 blood samples upon diagnosis and during treatment, we report gross abnormalities in cfDNA of COVID-19 patients, including elevated GC content, altered molecule size and end motif patterns. More importantly, such cfDNA fragmentation characteristics reflect patient-specific physiological changes during treatment. Further analysis on cfDNA tissue-of-origin tracing reveals frequent tissue injuries in COVID-19 patients, which is supported by clinical diagnoses. Hence, our work demonstrates and extends the translational merit of cfDNA fragmentation pattern as valuable analyte for effective treatment monitoring, as well as tissue injury assessment in COVID-19.

Hepatocellular carcinoma (HCC) is the fourth main reason of cancer-related death. Codonopsis pilosula is a commonly used traditional Chinese medicine (TCM) for patients with HCC. However, its potential mechanism for treatment of HCC remains unclear. Here, we used transcriptomics and network pharmacology to explore the potential molecular mechanisms of Codonopsis pilosula. In our study, twelve differentially expressed genes (DEGs) (5 upregulated and 7 downregulated) of Codonopsis pilosula treating HepG2 cells (a kind of HCC cell) were identified. Among the 12 DEGs, HMOX1 may play an essential role. Codonopsis pilosula mainly affects the mineral absorption pathway in HCC. We acquired 2957, 1877, and 255 targets from TCMID, SymMap, and TCMSP, respectively. Codonopsis pilosula could upregulate HMOX1 via luteolin, capsaicin, and sulforaphane. Our study provided new understanding of the potential pharmacological mechanisms of Codonopsis pilosula in treating HCC and pointed out a direction for further experimental research.

The Coronavirus Disease 2019 (COVID-19) epidemic was first detected in late-December 2019. So far, it has caused 203,815,431 confirmed cases and 4,310,623 deaths in the world. We collected sequences from 150,659 COVID-19 patients. Based on the previous phylogenomic analysis, we found three major branches of the virus RNA genomic mutation located in Asia, America, and Europe which is consistent with other studies. We selected sites with high mutation frequencies from Asia, America, and Europe. There are only 13 common mutation sites in these three regions. It infers that the viral mutations are highly dependent on their location and different locations have specific mutations. Most mutations can lead to amino acid substitutions, which occurred in 3/5’UTR, S/N/M protein, and ORF1ab/3a/8/10. Thus, the mutations may affect the pathogenesis of the virus. In addition, we applied an ARIMA model to predict the short-term frequency change of these top mutation sites during the spread of the disease. We tested a variety of settings of the ARIMA model to optimize the prediction effect of three patterns. This model can provide good help for predicting short-term mutation frequency changes.

Nitrogen fixation in soybean consumes a tremendous amount of energy, leading to substantial differences in energy metabolism and mitochondrial activities between nodules and uninoculated roots. While C-to-U RNA editing and intron splicing of mitochondrial transcripts are common in plant species, their roles in relation to nodule functions are still elusive. In this study, we performed RNA-seq to compare transcript profiles and RNA editing of mitochondrial genes in soybean nodules and roots. A total of 631 RNA editing sites were identified on mitochondrial transcripts, with 12% or 74 sites differentially edited among the transcripts isolated from nodules, stripped roots, and uninoculated roots. Eight out of these 74 differentially edited sites are located on the matR transcript, of which the degrees of RNA editing were the highest in the nodule sample. The degree of mitochondrial intron splicing was also examined. The splicing efficiencies of several introns in nodules and stripped roots were higher than in uninoculated roots. These include nad1 introns 2/3/4, nad4 intron 3, nad5 introns 2/3, cox2 intron 1, and ccmFc intron 1. A greater splicing efficiency of nad4 intron 1, a higher NAD4 protein abundance, and a reduction in supercomplex I + III₂ were also observed in nodules, although the causal relationship between these observations requires further investigation.

Exosomes are cell-derived lipid bilayer particles which are abundant in biological fluids. Exosome is an emerging source of biomarkers to diagnose various human diseases. Sequencing based exosomal studies could provide a comprehensive view of exosomal RNA and protein. To extracted these inclusions, exosomes should be isolated from the plasma first. Several exosome isolation methods were introduced since the discover of exosome. To promote the clinical application of exosomal inclusions, different isolation methods should be compared. We isolated exosomes from human plasma by using user-friendly and commercially available kits, SBI ExoQuick and QIAGEN exoRNeasy. Subsequently, small RNA sequencing was performed with two groups of isolated exosome samples and one group of plasma samples. No fundamental differences of exRNA yield between SC and EQ were found. In RNA profile analysis, the small RNA aligned reads, miRNA pattern, sample clustering varied as a result of methodological differences. Small RNA isolated by ExoQuick presented better data quality and RNA profile than exoRNeasy. This study compared sRNA sequencing data generated from two exosome isolation kits, it provides a reference for future small RNA studies and biomarker prediction in human plasma exosome.

Previously, several long non-coding RNAs (lncRNAs) were characterized as regulators in phosphate (Pi) starvation responses. However, systematic studies of novel lncRNAs involved in the Pi starvation signaling pathways have not been reported. Here, we used a genome-wide sequencing and bioinformatics approach to identify both poly(A) + and poly(A)– lncRNAs that responded to Pi starvation in Arabidopsis thaliana. We sequenced shoot and root transcriptomes of the Arabidopsis seedlings grown under Pi-sufficient and Pi-deficient conditions, and predicted 1212 novel lncRNAs, of which 78 were poly(A)– lncRNAs. By employing strand-specific RNA libraries, we discovered many novel antisense lncRNAs for the first time. We further defined 309 lncRNAs that were differentially expressed between P+ and P– conditions in either shoots or roots. Through Gene Ontology enrichment of the associated protein-coding genes (co-expressed or close on the genome), we found that many lncRNAs were adjacent or co-expressed with the genes involved in several Pi starvation related processes, including cell wall organization and photosynthesis. In total, we identified 104 potential lncRNA targets of PHR1, a key regulator for transcriptional response to Pi starvation. Moreover, we identified 16 candidate lncRNAs as potential targets of miR399, another key regulator of plant Pi homeostasis. Altogether, our data provide a rich resource of candidate lncRNAs involved in the Pi starvation regulatory network.