SARS‐CoV‐2 is a newly discovered beta coronavirus at the end of 2019, which is highly pathogenic and poses a serious threat to human health. In this paper, 1875 SARS‐CoV‐2 whole genome sequences and the sequence coding spike protein (S gene) sampled from the United States were used for bioinformatics analysis to study the molecular evolutionary characteristics of its genome and spike protein. The MCMC method was used to calculate the evolution rate of the whole genome sequence and the nucleotide mutation rate of the S gene. The results showed that the nucleotide mutation rate of the whole genome was 6.677 × 10⁻⁴ substitution per site per year, and the nucleotide mutation rate of the S gene was 8.066 × 10⁻⁴ substitution per site per year, which was at a medium level compared with other RNA viruses. Our findings confirmed the scientific hypothesis that the rate of evolution of the virus gradually decreases over time. We also found 13 statistically significant positive selection sites in the SARS‐CoV‐2 genome. In addition, the results showed that there were 101 nonsynonymous mutation sites in the amino acid sequence of S protein, including seven putative harmful mutation sites. This paper has preliminarily clarified the evolutionary characteristics of SARS‐CoV‐2 in the United States, providing a scientific basis for future surveillance and prevention of virus variants.

Colorectal cancer (CRC) is one of the most common malignancies and morbidity and mortality are increasing rapidly. Increasing evidence showed the close correlation between aberrant expression of certain RNAs and the occurrence and development of CRC. However, comprehensive analyses of differentially expressed profiles of linRNA in CRC based on large sample size have been lacking. In the present study, based on RNA-seq data obtained from the TCGA (The Cancer Genome Atlas) database, we identified 1176 lncRNAs, 245 miRNAs and 2083 mRNAs whichaberrantly expressed in the colorectal cancer tissues compared with the adjacent non-tumorous tissues. A Kaplan-Meier curve analysis was used to study the overall survival rate of the three RNA-related CRC patients. After constructing the ceRNA network, we performed the KEGG enrichment pathway analysis on ceRNA-related differentially expressed mRNAs and found that these mRNAs were remarkably enriched in the pathways associated with CRC. Combining the differentially expressed lncRNAs with clinical pathological variables of CRC patients, we also found that LINC00400 and LINC00355 not only contribute to the regulation of ceRNA network, but also show significantchanges in its expression in multiple CRC pathological stages, indicating that LINC00400 and LINC00355 can be considered as promising therapeutic targets for CRC.

Erythropoietin (EPO) is a glycoprotein hormone that plays a principal regulatory role in erythropoiesis and initiates cell homeostatic responses to environmental challenges. The Qinghai-Tibet Plateau is a natural laboratory for hypoxia adaptation. Gymnocypris dobula is a highly specialized plateau schizothoracine fish that is restricted to > 4500 m high-altitude freshwater rivers and ponds in the Qinghai-Tibet Plateau. The role of EPO in the adaptation of schizothoracine fish to hypoxia is unknown. The EPO and EPO receptor genes from G. dobula and four other schizothoracine fish from various altitudinal habitats were characterized. Schizothoracine EPOs are predicted to possess 2–3 N-glycosylation (NGS) sites, 4–5 casein kinase II phosphorylation (CK2) sites, 1–2 protein kinase C (PKC) phosphorylation sites, and four conserved cysteine residues within four helical domains, with variations in the numbers of NGS and CK2 sites in G. dobula. PAML analysis indicated a d N/d S value (ω) = 1.112 in the G. dobula lineage, and a few amino acids potentially under lineage-specific positive selection were detected within the G. dobula EPO. Similarly, EPO receptors of the two high-altitude schizothoracines (G. dobula and Ptychobarbus kaznakovi), were found to be statistically on the border of positive selection using the branch-site model (P-value = 0.096), and some amino acids located in the ligand-binding domain and the fibronectin type III domain were identified as potentially positive selection sites. Tissue EPO expression profiling based on transcriptome sequencing of three schizothoracines (G. dobula, Schizothorax nukiangensis Tsao, and Schizothorax prenanti) showed significant upregulation of EPO expression in the brain and less significantly in the gill of G. dobula. The elevated expression together with the rapid evolution of the EPO gene in G. dobula suggested a possible role for EPO in adaptation to hypoxia. To test this hypothesis, Gd-EPO and Sp-EPO were cloned into an expression vector and transfected into the cultured cell line 293 T. Significantly higher cell viability was observed in cells transfected with Gd-EPO than cells harboring Sp-EPO when challenged by hypoxia. The deduced EPO proteins of the schizothoracine fish contain characteristic structures and important domains similar to EPOs from other taxa. The presence of potentially positive selection sites in both EPO and EPOR in G. dobula suggest possible adaptive evolution in the ligand-receptor binding activity of the EPO signaling cascade in G. dobula. Functional study indicated that the EPO from high-altitude schizothoracine species demonstrated features of hypoxic adaptation by reducing toxic effects or improving cell survival when expressed in cultured cells, providing evidence of molecular adaptation to hypoxic conditions in the Qinghai-Tibet Plateau.

Glioblastoma (GBM) is the most common malignant primary brain tumor in adults. Currently, the prognosis of the patients with GBM is very poor and new molecular targets and treatment strategies are urgently needed to combat it. MicroRNA-148a (miR-148a) has been shown to be dysregulated in certain tumor types. However, the role of miR-148a in the pathogenesis of GBM is not fully understood. Here we comprehensively analyzed the roles of miR-148a, downstream DLGAP1, and their molecular pathways in GBM. We showed that miR-148a promote the proliferation and growth of GBM cells both in vitro and in vivo, and also induced the migration, invasion, and EMT (epithelial–mesenchymal transition) program of GBM cells by directly targeting DLGAP1. Furthermore, we identified 31 new miR-148a targets and found that miR-148a function was mainly involved in the cell adhesion signaling pathway and was associated with nervous system diseases. Our findings provide a new mechanism for miR-148a-mediated GBM cell invasion and reveal previously unreported targets of miR-148a as well as novel miR-148a-mediated regulatory networks in GBM. These results increase the understanding of the role of miR-148a in GBM and may lead to novel therapeutic strategies for GBM.

Recent studies have highlighted the value of microRNA-21 (miR-21) as a prognostic biomarker in gliomas. However, the role of miR-21 in predicting prognosis remains controversial. We performed a comprehensive study based upon a meta-analysis and The Cancer Genome Atlas (TCGA) glioma dataset validation to clarify the prognostic significance of miR-21 in glioma patients. In this study, we searched Embase, PubMed, Web of science, CNKI, SinoMed, and Wanfang databases for records up to May 2018. Relevant data were extracted to assess the correlation between miR-21 expression and survival in glioma patients. Pooled hazard ratios (HRs) with 95% confidence intervals (CIs) were used to describe association strength. We further used multivariate Cox regression analysis to assess miR-21 expression in the TCGA glioma dataset to validate the relationship between miR-21 expression and survival. Nine studies were included in the meta-analysis. Among them, eight studies provided data on overall survival (OS) with a pooled HR of 1.91 (95% CI: 1.34, 2.73), indicating that higher expression of miR-21 was significantly associated with worse OS in glioma patients; for the other study, which provided data on progression-free survival (PFS), no statistically significant HR was reported for PFS in the glioma patients (HR = 1.23, 95% CI: 0.41, 3.72). A multivariate Cox regression analysis of the miR-21 expression in the TCGA glioma dataset revealed that overexpression of miR-21 was a potential independent prognostic biomarker of poorer OS (HR = 1.27, 95% CI: 1.01, 1.59) and poorer PFS (HR = 1.46, 95% CI: 1.17, 1.82). Our findings suggest that higher expression of miR-21 is correlated with poorer glioma prognosis.

Background: Similar to humans, the zebrafish brain plays a central role in regulating sexual reproduction, maturation and sexual behavior. However, systematic studies of the dimorphic patterns of gene expression in the brain of male and female zebrafish are lacking. Results: In this study, the mRNA and lncRNA expression profiles were obtained from the brain tissue samples of the three male and three female zebrafish by high-throughput transcriptome sequencing. We identified a total of 108 mRNAs and 50 lncRNAs with sex-based differential expression. We randomly selected four differentially expressed genes for RT-qPCR verification and the results certified that the expression pattern showed a similar trend between RNA-seq and RT-qPCR results. Protein-protein interaction network analysis, Gene Ontology (GO) analysis, and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis were performed to obtain the biological significance of differentially expressed mRNA in the brain dimorphism of zebrafish. Finally, a Pearson correlation analysis was performed to construct the co-expression network of the mRNAs and lncRNAs. Conclusions: We found that 12 new lncRNAs not only have significant gender specificity in the brain of zebrafish, and this finding may provide a clue to further study of the functional difference between male and female zebrafish brain.