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.

Background: Recent studies showed that circRNAs are involved in the biological process of some human cancers. However, little is known about their functions in colorectal cancer (CRC). Methods: Here we first revealed the expression profiles of circRNAs in the CRC tissues and the adjacent non-tumorous tissues using high-throughput sequencing. The sequence feature, chromosome location, alternative splicing and other characteristics of the circRNAs were also explored. The miRNA and mRNA expression profiles were then obtained by analyzing relevant CRC data retrived from the TCGA database. We obtained and analyzed the competing endogenous RNA (ceRNA) network of the top three pairs of the largest up-regulated and down-regulated circRNAs. Results: In this study, we obtained 50,410 circRNAs in the CRC tissue and the adjacent non-tumor tissues, of which 33.7% (16,975) were new, and revealed differential changes in circRNA expression during colorectal carcinogenesis. We have identified six potential key circRNAs (circPIEZO1-3, hsa_circ_0067163, hsa_circ_0140188, hsa_circ_0002632, hsa_circ_0001998 and hsa_circ_0023990) associated with CRC, which play important roles in carcinogenesis as ceRNA for regulation of miRNA-mRNA network. In the subsequent KEGG analysis, several CRC-related pathways were found. Conclusions: Our findings advance the understanding of the pathogenesis of CRC from the perspective of circRNAs and provide some circRNAs as candidate diagnostic biomarkers or potential therapeutic targets.

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.

Cancer immunotherapy has become an important means of cancer treatment; however, the complex composition and heterogeneity of the colorectal cancer (CRC) microenvironment pose a huge challenge to cancer immunotherapy. Using data downloaded from The Cancer Genome Atlas database, the differences in the microenvironment between cases with low and high immune scores were examined at the multiomics level using bioinformatics approaches. It was revealed that the samples with high immune scores had good cytolytic immune responses and relatively abundant stromal cells, as well as significant infiltration of 22 immune cell subsets and a high non‑synonymous mutation burden and neoantigen burden. All of these characteristics contribute to a good prognosis. To better understand the impact of immune‑related genes on prognosis, differentially expressed genes between the low and high immune score samples were identified and it was concluded that serpin family Emember 1 (SERPINE1) and ubiquitin C‑terminal hydrolase L1 (UCHL1) may be potential therapeutic targets. The relationship between the immune score and the infiltration of 22 immune cells and the difference in SERPINE1 expression were verified by analyzing the GSE17536 and GSE21510 data sets downloaded from the Gene Expression Omnibus database.The present study analyzed the unique properties of immune cells in the CRC microenvironment, which are of great significance for understanding CRC immune mechanism and may also provide novel ideas for the targeted design of cancer immunotherapy.