Screening biomarkers for spinal cord injury using weighted gene co-expression network analysis and machine learning

Immune changes and inflammatory responses have been identified as central events in the pathological process of spinal co rd injury.They can greatly affect nerve regeneration and functional recovery.However,there is still limited understanding of the peripheral immune inflammato ry response in spinal cord inju ry.In this study.we obtained microRNA expression profiles from the peripheral blood of patients with spinal co rd injury using high-throughput sequencing.We also obtained the mRNA expression profile of spinal cord injury patients from the Gene Expression Omnibus(GEO)database(GSE151371).We identified 54 differentially expressed microRNAs and 1656 diffe rentially expressed genes using bioinformatics approaches.Functional enrichment analysis revealed that various common immune and inflammation-related signaling pathways,such as neutrophil extracellular trap formation pathway,T cell receptor signaling pathway,and nuclear factor-κB signal pathway,we re abnormally activated or inhibited in spinal cord inju ry patient samples.We applied an integrated strategy that combines weighted gene co-expression network analysis,LASSO logistic regression,and SVM-RFE algorithm and identified three biomarke rs associated with spinal cord injury:ANO10,BST1,and ZFP36L2.We verified the expression levels and diagnostic perfo rmance of these three genes in the original training dataset and clinical samples through the receiver operating characteristic curve.Quantitative polymerase chain reaction results showed that ANO20 and BST1 mRNA levels were increased and ZFP36L2 mRNA was decreased in the peripheral blood of spinal cord injury patients.We also constructed a small RNA-mRNA interaction network using Cytoscape.Additionally,we evaluated the proportion of 22 types of immune cells in the peripheral blood of spinal co rd injury patients using the CIBERSORT tool.The proportions of naive B cells,plasma cells,monocytes,and neutrophils were increased while the proportions of memory B cells,CD8^(+)T cells,resting natural killer cells,resting dendritic cells,and eosinophils were markedly decreased in spinal cord injury patients increased compared with healthy subjects,and ANO10,BST1 and ZFP26L2we re closely related to the proportion of certain immune cell types.The findings from this study provide new directions for the development of treatment strategies related to immune inflammation in spinal co rd inju ry and suggest that ANO10,BST2,and ZFP36L2 are potential biomarkers for spinal cord injury.The study was registe red in the Chinese Clinical Trial Registry(registration No.ChiCTR2200066985,December 12,2022).

Molecular Cloning and Sequence Analysis of Class Ⅱ Chitinase Gene in Leymus chinensis

[ Objective] The aim of this study was to clone Class Ⅱ chitinase gene in Leymus chinensis grown in saline land in Heilongjiang Province and analyze its sequence, which provided a foundation for further study on the biological function and application of chitinasa gene. [ Method] cDNA library of Leymus chinensis leaves were constructed, and its DNA sequence was determined or analyzed, while the homology of chitinasa gene and amino acid sequence was compared with that in GenBank. [ Result] One full length cDNA fragment with length of 996 bp was cloned from cDNA library of Leymus chinensis leaves. The length of ORF was 768 bp encoding 225 amino acids （GenBank accession number： EU344908）. The encoding products lacked CBD and C-terminal extension region from the view of structure, but had structural characters of Class Ⅱ chitinase gene, which indicated that amino acid sequence had high homology compared with Class Ⅱ chitinase gene of rye and wheat. The constructed recombinant vector pQE-LcChi2 could express a protein of 27 kD through induction, which was consistent with the deduced encoding product of pQE-LcChi2 gene. [ Conclusion] LcChl2 gene is an expression gene, which can express in E. coll.

Sequence Analysis of HA Genes from Three H9N2 Subtype Avian Influenza Viruses

[ Objective] The study aimed to understand the genetic characters of H9N2 subtype avian influenza viruses isolated in Belling area. [ Method] HA genes of three H9N2 subtype avian influenza viruses A/Chicken/Beijing/xu/00, A/Chicken/Beijing/bei/00 and A/Chicken/Beijing/ liu/00 were amplified by RT-PCR and then sequenced. [ Result] The results of phylogenetic analysis showed that A/Chicken/Beijing/xu/00, A/ Chicken/Beijing/bei/00 and A/Chicken/Beijing/liu/00 shared the nucleotide homologies of 84.8% （ Dk/HK/Y439/97 ） -98.0% （ Ck/GX17/00 ）, 85.1% （Dk/HK/Y439/97） - 99.1% （ Ck/GXl 7/00）, 90.7% （ Ck/BJ/3/01 ） - 99.1% （Ck/GX17/00） with the isolates from Hongkong and other are- as of Chinese Mainland respectively. At the same time, the analysis of amino acid indicated that the three isolates belonged to low pathogenic H9N2 isolates of avian origin. The 226^th amino acid of them were L （ Leu）, suggesting their high binding affinity to human cells. There were seven glyco- sylation sites in HA protein, five from HA1 and two from HA2. [ Cenclusien] By analysis at molecular level, it could be concluded that A/Chicken/ Beijing/xu/00, A/Chicken/Beijing/bei/00 and A/Chicken/Beijing/liu/00 were low pathogenic H9N2 isolates of avian origin.

Determination and Analysis of Mitochondrial ND2 Gene Sequence of Anas platyrhynchos

[Objective] The study was to analyze the phylogenesis of Anas platyrhynchos. [Method] Complete sequence of mitochondrial ND2 gene of 4 Anas platyrhynchos was determined by direct DNA sequencing based on PCR products. Combined with ND2 gene sequences of the Anas Linnaeus accessed in GenBank, phylogenetic tree was constructed by Neighbor-joining and maximum parsimony methods. [Result] The ND2 gene sequences of 4 Anas platyrhynchos were identical(1 041 bp in length; the nucleotide contents of A, G, T, and C were 28.91%, 13.35%, 20.75% and 36.98% respectively; A+T content approximated to that of C+G). Sequences of ND2 gene of mallard were same as spotbill duck, and had high homology with others. The phylogenetic trees indicated mallard and spotbilled duck were close in genetic relationship, both shared a haplotype; then Philippine duck, green-winged teal and northern pintail fell into branch ''A". [Conclusion] The domestic duck may be domesticated from mallard and spotbilled duck.

Cloning and Sequence Analysis of 16S rRNA and COI Gene in Mitochondrial DNA of Scortum barcoo

[Objective] The aim was to provide molecular biological basis for the researches on the genetic resources,genetic relationship among species and phyletic evolution of S.barcoo.[Method] PCR amplification and sequencing were used to study the 16S rRNA and COI gene fragments.[Result] As for 16S rRNA gene fragments,nucleotide sequences of 791 bp were obtained,and the A,T,G and C contents in this fragment were 31.6%,21.4%,20.4% and 26.7%respectively.As for the COI gene fragments,the size was 631 bp and the A,T,G And C contents were 27.7%,23.6%,29.8% and 18.9% respectively.Among these two gene fragments,the content of GC was lower than AT,and AT/GC of these two fragments was 1.13 and 1.05 respectively.[Conclusion] The genetic characteristics of gene fragments of 16S rRNA and COI of S.barcoo suggested that the variation in the same species was relatively low.The sequences of 16S rRNA gene in three samples the same,while the sequences of COI gene was also the same,indicating that these two gene of S.barcoo were conservative.

Cloning and Sequence Analysis of Lactate Dehydrogenase C(LDH-C)Gene from Black-lipped Pika in Western Sichuan Plateau

[Objective] The aim was to lay a foundation for research of contraceptive rodenticide with LDH-C4 as target protein. [ Method] EST sequence of LDH-C gene from black-lipped pika was cloned by PCR with degenerate primers; then the full length open reading frame （ORF） and 3＇UTR sequence were cloned by RACE technique. [ Result] The full length cDNA was 1 498 bp containing an ORF of 996 bp and a 3＇UTR of 486 bp. The ORF encoded a polypeptide of 332 amino acids. The alignment of LDH-C gene ORF nucleotide sequences from different species showed that the gene was conserved even between large taxons. The phylogenic tree showed that black-lipped pika LDH-C was closer to prima- tes and artiodactyla than to rodents. [Conclusion] cDNA sequence of LDH-C gene from black-lipped pika was cloned successfully.

Molecular Cloning, Sequence Analysis and Prokaryotic Expression of Ovine Activin Receptor Type IIB(ActRIIB) Gene

Objective] This study aimed to clone ovine activin receptor type llB （Ac-tRIIB） gene, construct the prokaryotic expression vector and express the target gene in vitro, thus providing basis for further function verification. [Method] The template cDNA which was reversely transcribed from total RNA of sheep liver tissue, was subjected to polymerase chain reaction （PCR） using specific primers of ActRIIB. The ful-length cDNA of ovine ActRIIB was obtained by pMD18-T cloning and sequencing for bioinformatics analysis. Ovine ActRIIB encoding sequence was subcloned into prokaryotic expression vector pET41a with restriction sites BamHl/Notl, and then transformed into BL21 （DE3）. The induced products by lPTG were analyzed with SDS-PAGE and Western Blot. [Result] The amplified ful-length cDNA of ovine Ac-tRllB gene was 1 564 bp in length （Genbank accession number： JX422071.1） with an open reading frame of 1 539 bp, encoding 512 amino acides. Ovine ActRllB shared the highest homology （99.6%） with bovine ActRllB. ActRllB had highly ho-mologous C-terminal domains and belonged to the TGFβ family. After prokaryotic expression, an approximately 92 kD His-tagged ActRllB recombinant protein was obtained, which was consistent with the excepted result. [Conclusion] Cloning and successful expression of ovine ActRIIB laid solid foundation for further investigation of its biological function.

Cloning and Sequence Analysis of HN and F Protein Genes from a Strain of Goose Paramyxovirus

[ Objective] The study was to clone HN and F genes from GX1 strain of goose paramyxovirus and analyze their sequences. [ Method] According to the full nucleotide sequence of GPV-SF02 strain of goose paramyxovirus, two pairs of pdmers were designed to amplify the HN and F genes from GX1 strain of goose paramyxovirus isolated from diseased goose in Guangxi Zhuang Autonomous Region; the amplified products were ligated into pMD18-T vector and sequenced. [ Result ] HN and F genes of this strain tested were 1 716 and 1 662 bp in full nucleotide length, respectively; both showed the homologues of about 97.3% with GPV- SF02 strain, of 80.3% -97.5% with strains LaSota, F48E9 and JS, of just 84.8% with Miyadera strain. [ Conclusion] The results show that isolated strain BX1 matches to virulent APMV-1 strain, belonging to genotype Ⅶ of APMV-1 strain.

Construction of Midgut Tissue-Specific cDNA Library of Bombyx mandarina M. and Isolation and Sequence Analysis of Serine Protease Gene Fragment

[Objective] The aim of the study is to construct cDNA library of midgut tissue of wild silkworm and isolate the serine protease gene. [Method] The midgut tissue-specific cDNA library of wild silkworm was constructed via cDNA Library Construction Kit (TaKaRa), then the serine protease gene was cloned via sequencing of the yielded cDNA library. [Result] The titer of cDNA library reached 6.2×105 pfu/ml, average insert size was about 1.2 kb. The serine protease gene cDNA fragment was obtained from colony sequencing (Accession No: EU672968). The nucleotide sequence of the cloned 854 bp fragment encodes 284 amino acid residues. Homology analyses showed some homology between putative amino acid sequence of the cloned fragment and amino acid sequences of serine proteases from other ten insects. [Conclusion] The results may avail to reveal the resistance of silkworm and wild silkworm to exotic intrusion.

Cloning and Sequence Analysis of Sheeppox Virus RPO30 Gene

[Objective] The study aimed to clone RPO30 gene from Sheeppox virus （SPPV） and predict the structure and function of the sequence. [Method] RPO30 gene of SPPV was cloned with PCR, linked into pMD18-T simple vector and then transformed into E. coli DH5a. In blue-white screen, the white colonies were selected to prepare plasmids. The positive plasmids were selected by double digestion and PCR, and then sequenced. Finally, the structure and function of the sequence obtained were predicted by bioinformatics methods. [Results] The RPO30 gene was successfully obtained; its ORF was 585 bp, encoding 193 amino acids and containing a recognition site for Hind III. Moreover, the SPPV RPO30 gene shared different homologies with the RPO30 gene sequences of other pox virus strains from GenBank database. Further analysis by biological software showed that in RPO30 protein, amino acids 4-12, 18-26, 50- 61, 68- 92 and 176-190 had a high possibility to form the active center, and acting to these regions was likely to inactivate the enzyme encoded by the sequence, thus to inhibit viral replication efficiently. [Conclusion] This study will lay foundation for further study on the structure and function of RPO30.

Sequence Comparison and Analysis of HA Gene of Four H9N2 Avian Influenza Virus Isolates

[ Objective] To determine the HA gene sequences of four H9N2 Avian influenza virus （AIV） strains and carry out comparative analysis so as to understand the difference and variation pattern of each strain from the angle of molecular biology and to know the distribution and epidemic law of H9N2 AIV. [Method] One pair of primers was designed referring to HA gene sequences of H9N2 AIV. The HA genes of A/Chicken/Hebei/WD/98 （H9N2; WD98 for short）, A/Chicken/Hebei/ZD/04 （H9N2; ZD04 for short））, A/Chicken/Beijing/MY/06 （H9N2; MY06 for short） ）, and A/Chicken/Beijing/PG/08 （H9N2; PG08 for short）） were amplified, cloned and sequenced. Then the HA gene sequences of these strains were compared with that of 10 H9N2 AIV stains in GenBank. [Result] The ORF of HA genes of the four strains was 1 683 bp in size, encoding 516 amino acids. The HA gene sequences of the four strains, WD98, MY06, PG08, and ZD04, were 82.6% -95.1%, 83.0% -99.0%, 82.7% -95.5%, and 81.3% -95.7% homologous to that of the 10 H9N2 AIV stains, respectively. And the homology of amino acid was respectively 86.6% -96.3%, 86.6% -97.9%, 87.0% -97.1%, and 86.9% -97.3%. [ Conclusion] The HA gene has greatly high homology among different strains.

Cloning and Sequence Analysis on IGF-1 Gene of Hubei White Swine

[Objective] The study aimed at cloning and analyzing the insulin-like growth factor-1 （IGF-1） gene from liver of Hubei white swine. [Method] The total RNA was extracted by using Trizol from the liver of Hubei white swine and used as template to amplify IGF-1 gene cDNA by RT-PCR. The cDNA product was cloned into pCRII vector, screened with blue-white colonies, digested with double enzymes and sequenced. [Result] The sequencing result indicated that the IGF-1 gene consisted of 607 nucleotides, containing 5＇-untranslated region at nucleotides 1-145, a complete ORF at nucleotides 146-538 encoding 130 amino acids, and 3＇-untranslated region at nucleotides 539-607. It shared 100% homology with the porcine IGF-1 gene reported by Muller et al. [Conclusion] The successful cloning and sequencing of the Hubei white swine IGF-1 gene confirmed that IGF-I gene was highly conserved, which provided technical basis for the use of transgenic technology for breeding of Hubei white swine.

Cloning,Sequence Analysis and Amino Acid Structure Prediction of TLR9 Gene from Wild Ovis ammon in Xinjiang

[Objective] The study aimed at cloning and identifying the toll receptor gene 9（TLR9） of wild Ovis ammon in Xinjiang,and predicting its structure and function.[Method] The TLR9 complete sequence of wild Ovis ammon was cloned from its peripheral blood by PCR technology.Then the PCR products were purified by agarose gel electrophoresis and then sequenced.Finally,the structure and function of TLR9 sequence were predicted by molecular biological software.[Result] The complete sequence of TLR9 gene was 3 192 bp in length,encoding 1 064 amino acids with a signal peptide composed of 30 amino acids;the leucine percentage reached as high as 18.5%.The TLR9 amino acid possibly contained three hydrophobic regions,at amino acids 455-475,740-760 and 780-800.The 3-D structure of TLR9 was constructed by the extracellular LRR domain and intracellular TIL domain（Toll/IL IR）.[Conclusion] The characteristics of TLR9 provided theoretical basis for further study on the TLR9 gene of wild Ovis ammon in Xinjiang.

Cloning and Sequence Analysis of 5′ Flanking Region and Exon of Inhibin α Precursor Gene in Goats

[Objective] The aim of this study was to reveal the relationship between inihibin （INH） α precursor gene and seasonal reproduction of goats, and investigate the evolutionary conservation of INHα precursor gene. [ Method] Cloning and sequence analysis of 5＇ flanking region and exon of inihibinα （INHE） precursor gene in twenty ewes between non-seasonal estrous breed （Haimen goats） and seasonal estrous breed （Anhui white goats） was analyzed in this study. [ Result] Compared with Anhui white goats, INHα precursor gene in Haimen goats had three SNP but no amino acid change, while its nucleotide homology was 99.7% and amino acid homology was 100%. The nucleotide homology of INHα precursor gene in goat, cattle, pig, person, chicken, horse, rat and dog ranged from 12.7% to 96.5%. [ Conclusion] INHα precursor gene tends to be highly conserved in species, and any change of nucleotide and amino acid maybe directly influence the function of the whole gene coding and regulation.

Cloning and Sequence Analysis of CmNAC Gene from Cucurbita moschata

[Objective] The aim was to clone the CmNAC gene from Cucurbita moschata and analyze the sequence characteristics. [Method] A pair of degenerate primers was designed based on the conserved sequences of NAC gene from Brassica napus, Lycopersicon esculentum and Capsicum annuum. NAC transcription factor gene was amplified by RT-PCR from Cucurbita moschata leaves and cloned into pMD-19T vector; then the recombinant clones were sequenced. Finally, the sequences of nucleic acid and amino acid were analyzed using BLAST and DNAMAN software. [Result] The NAC transcription factor gene cloned from C. moschata included 442 bp encoding 147 amino acids, named CmNAC. The NAC gene fragment contained a conserved region like other plant NAC genes and belonged to the NAC family ATAF1/2 subfamily. [Conclusion] The stress resistance related gene NAC cloned from C. moschata is a foundation for further study on the biological function of the gene and plant genetic engineering.

Isolation and Analysis of α-Gliadin Gene Coding Sequences from Triticum durum

Three coding sequences of gliadins genes, designed as Gli2_Dul, Gli2_Du2 and Gli2_Du3, were isolated from the genomic DNA of Triticum durum accessions CItr5083. Gli2_Dul and Gli2_Du2 contain 945 and 864 bp, encoding the mature proteins with 314 and 287 amino acid residues, respectively. Gli2_Du3 is recognized as a pseudogene due to the stop codon occurring in the coding region. The pseudogenes, commonly occurring in gliadins family, are attributed to the single base change C→T. The amino acid sequences deduced from these gene sequences were characterized with the typical structure of α-gliadin proteins, including the toxic sequences （PSQQQP）. The peptide fraction PF（Y）PP（Q）is thought to be an extra unit of repetitive domain, slightly diverging from the previous report. Six cysteine residues were observed within two unique domains. Phylogenetic analysis showed Gli2_Du2 and Gli2_Du3 were closely related to the genes on chromosome 6A, whereas Gli2_Dul seems to be more homologous with the genes on chromosome 6B.

Isolation and Sequence Analysis of HMW Glutenin Subunit 1Dy10.1 Ecoding Gene from Xinjiang Wheat (Triticum petropavlovskyi Udacz. et Migusch)

A novel HMW glutenin subunit gene 1Dy10. 1 was isolated and characterized from Xinjiang wheat （Triticum petropavlovskyi. Udacz. et Migusch） accession Daomai 2. The complete open reading frame （ORF） of 1Dy10. 1 was 1 965 bp, encoding 655 amino acids. The numbers and distribution of cysteines in 1Dy 10.1 were similar to those of 1Dy10 and other y-type subunits. In the N-terminal of 1Dy 10.1, an amino acid was changed from L （leucine） to P （proline） at position 55. The repetitive domain of 1Dy10.1 differed from those of known HMW subunits by substitutions, insertions or/and deletions involving single or more amino acid residues. In the repetitive domain of subunit 1Dyl 0.1, the deletion of tripeptide GQQ in the consensus unit PGQGQQ resulted in the appearance of the motif PGQ that have not been observed in other known y-type HMW subunits. In comparison with the subunit 1Dy 12, a deletion of dipeptide GQ, which occurred in subunit IDy10, was also observed in subunit 1Dy10.1. The cloned IDylO. 1 gene had been successfully expressed in Escherichia coli, and the expressed protein had the identical mobility with the endogenous subunit IDy10.1 from seed.

Vitis amuerensis CBF3 Gene Isolation,Sequence Analysis and Expression

The full length cDNA sequence of CBF3 （CRT/DRE-binding factor） was cloned from Vitis amurensis by reverse transcription polymerase chain reaction （RT-PCR） using the primers designed based on CBF genes available in GenBank. Sequence analysis showed that the gene had 854 bp long and its coding sequence contained 720 bp, encoding a protein with 239 amino acids and an AP2 structural domain. The molecular mass of CBF3 was predicted to be 25.9 kDa and its theoretical isoelectric point was 7.02 and aliphatic index was 59.29. The average hydropathicity of the protein was -0.551. The tertiary structures of CBF3 were also analyzed. The prokaryotic expression vector pGEX-4T-CBF3 containing CBF3 gene was constructed and CBF3 fusion protein （52 kDa） was produced in Escherichia coli after induction with 1 mmol L-1 IPTG. Further studies are needed to evaluate its potential application for improving plant resistance to cold and other stress condition such as drought and salinity.

Cloning and Sequence Analysis of a Novel Cold-Adapted Lipase Gene from Strain lip35 (Pseudomonas sp.)

A combination method of the usual-PCR and reverse-PCR for the cloning of a novel lipase gene directly from the total genomic DNA of strain lip35 （Pseudomonas sp.） is described, whereby a lipase gene （lip） was cloned directly from genomic DNA. The sequence data have been deposited in the GenBank and EMBL data bank with the accession number EU414288. The nucleotide sequence showed a major open reading frame encoding a 59-kDa protein of 566 amino acid residues, which contained a lipase consensus sequence GXSXG. The lipase lip had 74 and 70% homologies with the lipases of an uncultured bacterium and P. fluorescens PfO-1, respectively, but it did not show any overall homology with lipases from other origins. The functional lipase was obtained when the lip gene was expressed in Pichia pastoris GS115.

Gene cloning and sequence analysis of the cold-adapted chaperones DnaK and DnaJ from deep-sea psychrotrophic bacterium Pseudoalteromonas sp. SM9913

Pseudoalteromonas sp. SM9913 is a phychrotmphic bacterium isolated from the deep-sea sediment. The genes encoding chaperones DnaJ and DnaK of P. sp. SM9913 were cloned by normal PCR and TAIL - PCR （GenBank accession Nos DQ640312, DQ504163 ）. The chaperones DnaJ and DnaK from the strain SM9913 contain such conserved domains as those of many other bacteria, and show some cold-adapted characteristics in their structures when compared with those from psychro-, meso-and themophilic bacteria. It is indicated that chaperones DnaJ and DnaK of P. sp. SM9913 may be adapted to low temperature in deep-sea and function well in assisting folding, assembling and translocation of proteins at low temperature. This research lays a foundation for the further study on the cold-adapted mechanism of chaperones DnaJ and DnaK of cold-adapted microorganisms.