Single‑cell RNA sequencing opens a new era for cotton genomic research and gene functional analysis

Single-cell RNA sequencing(scRNA-seq)is one of the most advanced sequencing technologies for studying transcriptome landscape at the single-cell revolution.It provides numerous advantages over traditional RNA-seq.Since it was first used to profile single-cell transcriptome in plants in 2019,it has been extensively employed to perform different research in plants.Recently,scRNA-seq was also quickly adopted by the cotton research community to solve lots of scientific questions which have been never solved.In this comment,we highlighted the significant progress in employing scRNA-seq to cotton genetic and genomic study and its future potential applications.

Single-cell RNA sequencing analysis of the retina under acute high intraocular pressure

High intraocular pressure causes retinal ganglion cell injury in primary and secondary glaucoma diseases,yet the molecular landscape characteristics of retinal cells under high intraocular pressure remain unknown.Rat models of acute hypertension ocular pressure were established by injection of cross-linked hyaluronic acid hydrogel(Healaflow■).Single-cell RNA sequencing was then used to describe the cellular composition and molecular profile of the retina following high intraocular pressure.Our results identified a total of 12 cell types,namely retinal pigment epithelial cells,rod-photoreceptor cells,bipolar cells,Müller cells,microglia,cone-photoreceptor cells,retinal ganglion cells,endothelial cells,retinal progenitor cells,oligodendrocytes,pericytes,and fibroblasts.The single-cell RNA sequencing analysis of the retina under acute high intraocular pressure revealed obvious changes in the proportions of various retinal cells,with ganglion cells decreased by 23%.Hematoxylin and eosin staining and TUNEL staining confirmed the damage to retinal ganglion cells under high intraocular pressure.We extracted data from retinal ganglion cells and analyzed the retinal ganglion cell cluster with the most distinct expression.We found upregulation of the B3gat2 gene,which is associated with neuronal migration and adhesion,and downregulation of the Tsc22d gene,which participates in inhibition of inflammation.This study is the first to reveal molecular changes and intercellular interactions in the retina under high intraocular pressure.These data contribute to understanding of the molecular mechanism of retinal injury induced by high intraocular pressure and will benefit the development of novel therapies.

New challenges in hepatocellular carcinoma:A role for PIWIinteracting RNAs?

Hepatocellular carcinoma(HCC)is the most common and deadliest subtype of liver cancer worldwide and,therefore,poses an enormous threat to global health.Understanding the molecular mechanisms underlying the development and progression of HCC is central to improving our clinical approaches.PIWIinteracting RNAs(piRNAs)are a class of small non-coding RNAs that bind to PIWI family proteins to regulate gene expression at transcriptional and posttranscriptional levels.A growing body of work shows that the dysregulation of piRNAs plays a crucial role in the progression of various human cancers.In this editorial,we report on the current knowledge of HCC-associated piRNAs and their potential clinical utility.Based on the editorial by Papadopoulos and Trifylli,on the role and clinical evaluation of exosomal circular RNAs in HCC,we highlight this other emerging class of non-coding RNAs.

Integrated analysis of single-cell and bulk RNA-seq establishes a novel signature for prediction in gastric cancer

BACKGROUND Single-cell sequencing technology provides the capability to analyze changes in specific cell types during the progression of disease.However,previous single-cell sequencing studies on gastric cancer(GC)have largely focused on immune cells and stromal cells,and further elucidation is required regarding the alterations that occur in gastric epithelial cells during the development of GC.AIM To create a GC prediction model based on single-cell and bulk RNA sequencing(bulk RNA-seq)data.METHODS In this study,we conducted a comprehensive analysis by integrating three singlecell RNA sequencing(scRNA-seq)datasets and ten bulk RNA-seq datasets.Our analysis mainly focused on determining cell proportions and identifying differentially expressed genes(DEGs).Specifically,we performed differential expression analysis among epithelial cells in GC tissues and normal gastric tissues(NAGs)and utilized both single-cell and bulk RNA-seq data to establish a prediction model for GC.We further validated the accuracy of the GC prediction model in bulk RNA-seq data.We also used Kaplan–Meier plots to verify the correlation between genes in the prediction model and the prognosis of GC.RESULTS By analyzing scRNA-seq data from a total of 70707 cells from GC tissue,NAG,and chronic gastric tissue,10 cell types were identified,and DEGs in GC and normal epithelial cells were screened.After determining the DEGs in GC and normal gastric samples identified by bulk RNA-seq data,a GC predictive classifier was constructed using the Least absolute shrinkage and selection operator(LASSO)and random forest methods.The LASSO classifier showed good performance in both validation and model verification using The Cancer Genome Atlas and Genotype-Tissue Expression(GTEx)datasets[area under the curve(AUC)_min=0.988,AUC_1se=0.994],and the random forest model also achieved good results with the validation set(AUC=0.92).Genes TIMP1,PLOD3,CKS2,TYMP,TNFRSF10B,CPNE1,GDF15,BCAP31,and CLDN7 were identified to have high importance values in multiple GC predictive models,and KM-PLOTTER analysis showed their relevance to GC prognosis,suggesting their potential for use in GC diagnosis and treatment.CONCLUSION A predictive classifier was established based on the analysis of RNA-seq data,and the genes in it are expected to serve as auxiliary markers in the clinical diagnosis of GC.

Applications of single-cell RNA sequencing in spermatogenesis and molecular evolution

Spermatogenic cell heterogeneity is determined by the complex process of spermatogenesis differentiation.However,effectively revealing the regulatory mechanisms underlying mammalian spermatogenic cell development and differentiation via traditional methods is difficult.Advances in technology have led to the emergence of many single-cell transcriptome sequencing protocols,which have partially addressed these challenges.In this review,we detail the principles of 10x Genomics technology and summarize the methods for downstream analysis of single-cell transcriptome sequencing data.Furthermore,we explore the role of single-cell transcriptome sequencing in revealing the heterogeneity of testicular ecological niche cells,delineating the establishment and disruption of testicular immune homeostasis during human spermatogenesis,investigating abnormal spermatogenesis in humans,and,ultimately,elucidating the molecular evolution of mammalian spermatogenesis.

Single-Nucleus RNA Sequencing Reveals Cardiac Macrophage Landscape in Hypoplastic Left Heart Syndrome

Background:Hypoplastic left heart syndrome(HLHS)is one of the most challenging congenital heart diseases in clinical treatment.In cardiac tissues,resident macrophages fulfill critical functions in maintaining a stable cardiac state and have strong regenerative capacity and organ specificity.However,the molecular mechanisms of macro-phages in HLHS remained unclear.Methods:Single-nucleus RNA sequencing(snRNA-seq)data of HLHS and healthy control(donors)samples obtained from the Gene Expression Omnibus(GEO)database were normalized and clustered using the Seurat package.The“FindMarkers”function was used to screen differentially expressed genes(DEGs)between the HLHS and donor groups and to analyze the functional enrichment of the set of genes of interest.Finally,cell-cell communication,pseudotime,and single-cell regulatory network inference and cluster-ing(SCENIC)analyses were used to study the mechanisms of macrophages in HLHS.Results:Based on the snRNA-seq data of HLHS and donors,we identified a total of 9 cell clusters,among which the proportion of macrophages was significantly less in the HLHS group than in the control group.Subdivision of macrophage subpopulations(Macrophages 1,2,and 3)showed that Macrophages 1 was mainly involved in nervous system development,angiogenesis,and apoptotic processes.In addition,analysis of communication between Macro-phages 1 and cardiomyocytes revealed that ligand-acceptor pairs such as GAS6/AXL,IL6,IGF1,THY1,and L1CAM were present only in the donor group.Finally,pesudotime and SCENIC analyses demonstrated that FOXO3 and ELF2 played a critical role for Macrophages 1 to maintain cardiac function in patients with HLHS.Conclusion:Our study improved the current understanding of the molecular mechanisms of macrophage devel-opment in HLHS,showing that manipulating the regulatory role of macrophages in the heart can be a novel treat-ment for HLHS.

Analysis of RNA-Dependent RNA Polymerase Sequence of Infectious Flacherie Virus Isolated in China and Its Expression in BmN Cells

Full gene sequence of RNA-dependent RNA polymerase （RdRp） from Bombyx mori infectious flacherie virus isolated in Zhejiang Province, China （Zhejiang01/CHN/2002） was cloned. The sequence was 1 920 nucleotides in length coding 639 amino acid residues. Sequences comparison of RdRp showed Zhejiang01/CHN/2002 was 99.7% nucleotide sequence and 99.1% amino acids sequence homology with Japanese strain. The RdRp sequence was aligned with 8 representative picorna（-like） viruses and 8 highly conserved regions were detected. The result indicated their relevance function. Phylogenetic tree of 14 picorna（-like） viruses which RdRp presumed protein sequences revealed that the viruses from Iflavirus genus formed an independent clade. The RdRp was successfully expressed in BmN cells using Bac-to-Bac expression system.

Identification, characterization and full-length sequence analysis of a novel endornavirus in common sunflower (Helianthus annuus L.)

To identify the possible quarantine viruses in seven common sunflower varieties imported from the United States of America and the Netherlands, we tested total RNAs extracted from the leaf tissues using next-generation sequencing of small RNAs. After analysis of small RNA sequencing data, no any quarantine virus was found, but a double-stranded RNA（dsRNA） molecule showing typical genomic features of endornavirus was detected in two varieties, X3939 and SH1108. Full-length sequence and phylogenetic analysis showed that it is a novel endornavirus, temporarily named as Helianthus annuus alphaendornavirus（HaEV）. Its full genome corresponds to a 14 662-bp dsRNA segment, including a 21-nt 5′ untranslated region（UTR）, 3＇ UTR ending with the unique sequence CCCCCCCC and lacking a poly（A） tail. An open reading frame（ORF） that encodes a deduced 4 867 amino acids（aa） polyprotein with three domains： RdRP, Hel and UGT（UDP-glycosyltransferase）. HaEV mainly distributed in the cytoplasm but less in the nucleus of leaf cells by fluorescence in situ hybridization（FISH） experiment. This virus has a high seed infection rate in the five varieties, X3907, X3939, A231, SH1108 and SR1320. To our knowledge, this is the first report about the virus of the family Endornaviridae in the common sunflower.

Identification and validation of a pyroptosis-related prognostic model for colorectal cancer based on bulk and single-cell RNA sequencing data

BACKGROUND Pyroptosis impacts the development of malignant tumors,yet its role in colorectal cancer(CRC)prognosis remains uncertain.AIM To assess the prognostic significance of pyroptosis-related genes and their association with CRC immune infiltration.METHODS Gene expression data were obtained from The Cancer Genome Atlas(TCGA)and single-cell RNA sequencing dataset GSE178341 from the Gene Expression Omnibus(GEO).Pyroptosis-related gene expression in cell clusters was analyzed,and enrichment analysis was conducted.A pyroptosis-related risk model was developed using the LASSO regression algorithm,with prediction accuracy assessed through K-M and receiver operating characteristic analyses.A nomo-gram predicting survival was created,and the correlation between the risk model and immune infiltration was analyzed using CIBERSORTx calculations.Finally,the differential expression of the 8 prognostic genes between CRC and normal samples was verified by analyzing TCGA-COADREAD data from the UCSC database.RESULTS An effective pyroptosis-related risk model was constructed using 8 genes-CHMP2B,SDHB,BST2,UBE2D2,GJA1,AIM2,PDCD6IP,and SEZ6L2(P<0.05).Seven of these genes exhibited differential expression between CRC and normal samples based on TCGA database analysis(P<0.05).Patients with higher risk scores demonstrated increased death risk and reduced overall survival(P<0.05).Significant differences in immune infiltration were observed between low-and high-risk groups,correlating with pyroptosis-related gene expression.CONCLUSION We developed a pyroptosis-related prognostic model for CRC,affirming its correlation with immune infiltration.This model may prove useful for CRC prognostic evaluation.

单细胞测序与Bulk RNA测序联合分析构建膀胱癌风险预后模型

目的本研究致力于开发一种基于膀胱癌核心细胞差异性表达基因的预后风险评估模型。通过在多个数据集中验证其效力,旨在为膀胱癌患者提供一个新的临床应用工具,用于预后风险评估。方法本研究综合利用了单细胞和Bulk RNA测序数据。我们首先从GEO数据库下载并分析了相关的膀胱癌单细胞和芯片RNA数据集。通过生物信息学方法,我们鉴定了核心细胞的差异性表达基因,并对其进行了功能及通路富集分析。基于这些分析,我们使用单变量和多变量Cox回归方法筛选出与膀胱癌预后显著相关的关键基因,并据此构建了一个预后风险评估模型。该模型在TCGA-BLCA数据集中进一步进行了效力验证。结果经过全面的生物信息学分析,我们鉴定出了223个核心细胞的差异性表达基因。这些基因在细胞外基质的结构和功能方面发挥着重要作用。构建的预后风险评估模型包括5个独立的预后相关基因(MFAP5、PDE4D、ISG15、ADAMTS1和FGL2)。在GEO和TCGA-BLCA数据集中的验证结果表明,该模型具有良好的预测效力,为膀胱癌患者的预后评估提供新型生物学标志工具。结论本研究成功开发了一个基于5个关键基因标志物的膀胱癌预后风险评估模型并具有良好的预测效力。此模型的开发为膀胱癌的生物学研究和临床预后评估提供了一个新的工具,有助于更好地理解膀胱癌的生物学特性并指导患者的个性化治疗。

Expression changes of mi RNA-regulated genes associated with the formation of the leafy head in cabbage

The vegetative development of cabbage(Brassica oleracea var.capitata)passes through seedling,rosette,folding and heading stages.Leaves that form the rosette are large and mostly flat.In the following developmental stages,the plants produce leaves that curve inward to produce the leafy head.Many microRNAs and their target genes have been described participating in leaf development and leaf curvature.The aim of this study is to investigate the role of miRNA-regulated genes in the transition from the rosette to the heading stage.We compared the mi RNA and gene abundances between emerging rosette and heading leaves.To remove transcripts(miRNAs and genes)whose regulation was most likely associated with plant age rather than the change from rosette to heading stage,we utilized a non-heading collard green(B.oleracea var.acephala)morphotype as control.This resulted in 33 DEMs and 1998 DEGs with likely roles in the transition from rosette to heading stage in cabbage.Among these 1998 DEGs,we found enriched GO terms related to DNA-binding transcription factor activity,transcription regulator activity,iron ion binding,and photosynthesis.We predicted the target genes of these 33 DEMs and focused on the subset that was differentially expressed(1998DEGs)between rosette and heading stage leaves to construct mi RNA-target gene interaction networks.Our main finding is a role for miR396b-5p targeting two Arabidopsis thaliana orthologues of GROWTH REGULATING FACTORs 3(GRF3)and 4(GRF4)in pointed cabbage head formation.

Fine-needle aspiration technique under endoscopic ultrasound guidance:A technical approach for RNA profiling of pancreatic neoplasms

BACKGROUND Early diagnosis of pancreatic ductal adenocarcinoma(PDAC)has been a longstanding challenge.The prognosis of patients with PDAC depends on the stage at diagnosis.It is necessary to identify biomarkers for the detection and differentiation of pancreatic tumors and optimize PDAC sample preparation procedures for DNA and RNA analysis.Most molecular studies are done using paraffin-embedded blocks;however,the integrity of DNA and RNA is often compromised in this format.Moreover,RNA isolated from human pancreatic tissue samples is generally of low quality,in part,because of the high concentration of endogenous pancreatic RNAse activity present.AIM To assess the potential of endoscopic ultrasound-guided fine-needle aspiration(EUS-FNA)to obtain specimens from pancreatic neoplasms for subsequent RNA molecular profiling,including next-generation sequencing(NGS).METHODS Thirty-four EUS-FNA samples were included in this study:PDAC(n=15),chronic pancreatitis(n=5),pancreatic cysts(n=14),mucinous cysts(mucinous cystic neoplasia/intraductal papillary mucinous neoplasia)n=7,serous cystic neoplasms n=5,and pseudocysts n=2.Cyst material consisted of cyst fluid and cyst wall samples obtained by through-the-needle biopsy(TTNB).Samples were stored at -80℃ until analysis.RNA purity(A260/230,A260/280 ratios),concentration,and integrity(RIN)were assessed.Real-time polymerase chain reaction was conducted on all samples,and small RNA libraries were prepared from solid mass samples.RESULTS RNA was successfully extracted from 29/34(85%)EUS-FNA samples:100% pancreatic adenocarcinoma samples,100% chronic pancreatitis samples,70% pancreatic fluid cyst samples,and 50%TTNB samples.The relative expression of GAPDH and HPRT were obtained for all successfully extracted RNA samples(n=29)including lowquality RNA specimens.Low concentration and nonoptimal RIN values(no less than 3)of RNA extracted from EUS-FNA samples did not prevent NGS library preparation.The suitability of cyst fluid samples for RNA profiling varied.The quality of RNA extracted from mucinous cyst fluid had a median RIN of 7.7(5.0-8.2),which was compatible with that from solid neoplasms[6.2(0-7.8)],whereas the quality of the RNA extracted from all fluids of serous cystic neoplasms and TTNB samples had a RIN of 0.CONCLUSION The results demonstrate the high potential of EUS-FNA material for RNA profiling of various pancreatic lesions,including low-quality RNA specimens.

Peripheral blood RNA biomarkers can predict lesion severity in degenerative cervical myelopathy

Degenerative cervical myelopathy is a common cause of spinal cord injury,with longer symptom duration and higher myelopathy severity indicating a worse prognosis.While numerous studies have investigated serological biomarkers for acute spinal cord injury,few studies have explored such biomarkers for diagnosing degenerative cervical myelopathy.This study involved 30 patients with degenerative cervical myelopathy(51.3±7.3 years old,12 women and 18 men),seven healthy controls(25.7±1.7 years old,one woman and six men),and nine patients with cervical spondylotic radiculopathy(51.9±8.6 years old,three women and six men).Analysis of blood samples from the three groups showed clear differences in transcriptomic characteristics.Enrichment analysis identified 128 differentially expressed genes that were enriched in patients with neurological disabilities.Using least absolute shrinkage and selection operator analysis,we constructed a five-gene model(TBCD,TPM2,PNKD,EIF4G2,and AP5Z1)to diagnose degenerative cervical myelopathy with an accuracy of 93.5%.One-gene models(TCAP and SDHA)identified mild and severe degenerative cervical myelopathy with accuracies of 83.3%and 76.7%,respectively.Signatures of two immune cell types(memory B cells and memory-activated CD4^(+)T cells)predicted levels of lesions in degenerative cervical myelopathy with 80%accuracy.Our results suggest that peripheral blood RNA biomarkers could be used to predict lesion severity in degenerative cervical myelopathy.

The Influence of Aerial Exposure on Sea Anemones Aulactinia veratra Mucin Genes Expression Using the RNA Sequencing

Mucin genes are the main component of mucus. The sea anemone species, Aulactinia veratra (Phylum Cnidaria) contains different types of mucin genes. In the intertidal zone, A. veratra is found to be exposed to air during the low tide and produces large quantities of mucus as an external covering. The relation between low tide and mucus secretion is still unclear, and what is the role of mucin during arial exposure is not yet investigated. This study hypothesised that the mucin genes in A. veratra would have significantly high expression in response to aerial exposure. Therefore, the aim of current study was to examine and analyses the response of A. veratra mucins in response to an experiment involving three hours of aerial exposure. To achieve this, aim the RNA-sequencing and bioinformatics analyses were used to examine the expression profile of A. veratra mucin genes in response to aerial exposure. The generated results have shown that, Mucin4-like and mucin5B-like were up-regulated in response to the three hours of aerial exposure in A. veratra. This finding shows a significant role of mucin5B-like and mucin4-like genes in response to air stress at low tide. The data generated from this study could be used in conjunction with future mucin gene studies of sea anemones and other cnidarians to compare A. veratra mucin gene expression results across time, and to extend our understanding of mucin stress response in this phylum.

Prediction of candidate small non-coding RNAs inAgrobacterium by computational analysis

Small non-coding RNAs with important regulatory roles are not confined to eukaryotes. Recent work has uncovered a growing number of bacterial small RNAs （sRNAs）, some of which have been shown to regulate critical cellular processes. Computational approaches, in combination with molecular experiments, have played an important role in the identification of these sRNAs. At present, there is no information on the presence of small non-coding RNAs and their genes in the Agrobacterium tumefaciens genome. To identify potential sRNAs in this important bacterium, deep sequencing of the short RNA populations isolated from Agrobacterium tumefaciens C58 was carried out. From a data set of more than 10,000 short sequences, 16 candidate sRNAs have been tentatively identified based on computational analysis. All of these candidates can form stem-loop structures by RNA folding predictions and the majority of the secondary structures are rich in GC base-pairs：：Some are followed by a short stretch of U residues, indicative of a rho-independent transcription terminator, whereas some of the short RNAs are found in the stem region of the hairpin, indicative of eukaryotic-like sRNAs. Experimental strategies will need to be used to verify these candidates. The study of an expanded list of candidate sRNAs in Agrobacterium will allow a more complete understanding of the range of roles played by regulatory RNAs in prokaryotes.

Analysis of long noncoding RNA-associated competing endogenous RNA network in glucagon-like peptide-1 receptor agonist-mediated protection inβcells

BACKGROUND Long noncoding RNAs(lncRNAs)and mRNAs are widely involved in various physiological and pathological processes.The use of glucagon-like peptide-1 receptor agonists(GLP-1RAs)is a novel therapeutic strategy that could promote insulin secretion and decrease the rate ofβ-cell apoptosis in type 2 diabetes mellitus(T2DM)patients.However,the specific lncRNAs and mRNAs and their functions in these processes have not been fully identified and elucidated.AIM To identify the lncRNAs and mRNAs that are involved in the protective effect of GLP-1RA inβcells,and their roles.METHODS Rat gene microarray was used to screen differentially expressed(DE)lncRNAs and mRNAs inβcells treated with geniposide,a GLP-1RA.Gene Ontology(GO)and Kyoto Encyclopedia of Genes and Genomes(KEGG)pathway enrichment analyses were performed to assess the underlying functions of DE mRNAs.Hub mRNAs were filtered using the STRING database and the Cytoscape plugin,CytoHubba.In order to reveal the regulatory relationship between lncRNAs and hub mRNAs,their co-expression network was constructed based on the Pearson coefficient of DE lncRNAs and mRNAs,and competing endogenous RNA(ceRNA)mechanism was explored through miRanda and TargetScan databases.RESULTS We identified 308 DE lncRNAs and 128 DE mRNAs with a fold change filter of≥1.5 and P value<0.05.GO and KEGG pathway enrichment analyses indicated that the most enriched terms were G-protein coupled receptor signaling pathway,inflammatory response,calcium signaling pathway,positive regulation of cell proliferation,and ERK1 and ERK2 cascade.Pomc,Htr2a,and Agtr1a were screened as hub mRNAs using the STRING database and the Cytoscape plugin,CytoHubba.This result was further verified using SwissTargetPrediction tool.Through the co-expression network and competing endogenous(ceRNA)mechanism,we identified seven lncRNAs(NONRATT027738,NONRATT027888,NONRATT030038,etc.)co-expressed with the three hub mRNAs(Pomc,Htr2a,and Agtr1a)based on the Pearson coefficient of the expression levels.These lncRNAs regulated hub mRNA functions by competing with six miRNAs(rno-miR-5132-3p,rno-miR-344g,rno-miR-3075,etc.)via the ceRNA mechanism.Further analysis indicated that lncRNA NONRATT027738 interacts with all the three hub mRNAs,suggesting that it is at a core position within the ceRNA network.CONCLUSION We have identified key lncRNAs and mRNAs,and highlighted here how they interact through the ceRNA mechanism to mediate the protective effect of GLP-1RA inβcells.

Circular RNA expression and the competitive endogenous RNA network in pathological,age-related macular degeneration events:A cross-platform normalization study

Age-related macular degeneration(AMD)causes irreversible blindness in people aged over 50 worldwide.The dysfunction of the retinal pigment epithelium is the primary cause of atrophic AMD.In the current study,we used the ComBat and Training Distribution Matching method to integrate data obtained from the Gene Expression Omnibus database.We analyzed the integrated sequencing data by the Gene Set Enrichment Analysis.Peroxisome and tumor necrosis factor-α(TNF-α)signaling and nuclear factor kappa B(NF-κB)were among the top 10 pathways,and thus we selected them to construct AMD cell models to identify differentially expressed circular RNAs(circRNAs).We then constructed a competing endogenous RNA network,which is related to differentially expressed circRNAs.This network included seven circRNAs,15 microRNAs,and 82 mRNAs.The Kyoto Encyclopedia of Genes and Genomes analysis of mRNAs in this network showed that the hypoxia-inducible factor-1(HIF-1)signaling pathway was a common downstream event.The results of the current study may provide insights into the pathological processes of atrophic AMD.

Expression of long non-coding RNAs in complete transection spinal cord injury: a transcriptomic analysis

Long non-coding RNAs(lncRNAs)are abundantly expressed in the central nervous system and exert a critical role in gene regulation via multiple biological processes.To uncover the functional significance and molecular mechanisms of lncRNAs in spinal cord injury(SCI),the expression signatures of lncRNAs were profiled using RNA sequencing(RNA-seq)technology in a Sprague-Dawley rat model of the 10th thoracic vertebra complete transection SCI.Results showed that 116 of 14,802 detected lncRNAs were differentially expressed,among which 16—including eight up-regulated(H19,Vof16,Hmox2-ps1,LOC100910973,Ybx1-ps3,Nnat,Gcgr,LOC680254)and eight down-regulated(Rmrp,Terc,Ngrn,Ppp2r2b,Cox6a2,Rpl37a-ps1,LOC360231,Rpph1)—demonstrated fold changes>2 in response to transection SCI.A subset of these RNA-seq results was validated by quantitative real-time PCR.The levels of 821 mRNAs were also significantly altered post-SCI;592 mRNAs were up-regulated and 229 mRNAs were down-regulated by more than 2-fold.Gene Ontology(GO)and Kyoto Encyclopedia of Genes and Genomes(KEGG)analyses showed that differentially expressed mRNAs were related to GO biological processes and molecular functions such as injury and inflammation response,wound repair,and apoptosis,and were significantly enriched in 15 KEGG pathways,including cell phagocytosis,tumor necrosis factor alpha pathway,and leukocyte migration.Our results reveal the expression profiles of lncRNAs and mRNAs in the rat spinal cord of a complete transection model,and these differentially expressed lncRNAs and mRNAs represent potential novel targets for SCI treatment.We suggest that lncRNAs may play an important role in the early immuno-inflammatory response after spinal cord injury.This study was approved by the Administration Committee of Experimental Animals,Guangdong Province,China.

A new Multilocus Sequence Analysis Scheme for Mycobacterium tuberculosis

Objective Tuberculosis remains one of the most serious infectious diseases in the world. In this study, a scheme of Mycobacterium tuberculosis （M. tuberculosis） multilocus sequence analysis （MLSA） was established for the phylogenetic and epidemiology analysis. Methods To establish the scheme of M. tuberculosis MLSA, the genome of H37Rv, CCDC5079 and CCDC5180 were compared, and some variable genes were chosen to be the MLSA typing scheme. 44 M. tuberculosis clinical isolates were typed by MLSA, IS6110-RFLP, and soligotyping, to evaluate the MLSA methods. Results After comparison of the genome, seven high discrimination gene loci （recX, rpsL, rmlC, rpmG1, mprA, gcvH, ideR） were chosen to be the MLSA typing scheme finally. 11 variable SNP sites of those seven genes were found among the 44 M. tuberculosis isolate strains and 11 sequence types （STs） were identified. Based on the Hunter-Gaston Index （HGI）, MLSA typing was not as good for discrimination at the strain level as IS6110-RFLP, but the HGI was much better than that of spoligotyping. In addition, the MEGA analysis result of MLSA data was similar to spoligotyping/PGG lineage, showing a strong phylogenetic signal in the modern strains of M. tuberculosis. The MLSA data analysis by eBURST revealed that 4 sequence types （ST） came into a main cluster, showing the major clonal complexes in those 44 strains. Conclusion MLSA genotyping not only can be used for molecular typing, but also is an ideal method for the phylogenetic analysis for M. tuberculosis.