Transcriptome sequencing reveals novel biomarkers and immune cell infiltration in esophageal tumorigenesis

BACKGROUND Esophageal squamous cell carcinoma(ESCC)is one of the most common malignancies worldwide,and its development comprises a multistep process from intraepithelial neoplasia(IN)to carcinoma(CA).However,the critical regulators and underlying molecular mechanisms remain largely unknown.AIM To explore the genes and infiltrating immune cells in the microenvironment that are associated with the multistage progression of ESCC to facilitate diagnosis and early intervention.METHODS A mouse model mimicking the multistage development of ESCC was established by providing warter containing 4-nitroquinoline 1-oxide(4NQO)to C57BL/6 mice.Moreover,we established a control group without 4NQO treatment of mice.Then,transcriptome sequencing was performed for esophageal tissues from patients with different pathological statuses,including low-grade IN(LGIN),high-grade IN(HGIN),and CA,and controlled normal tissue(NOR)samples.Differentially expressed genes(DEGs)were identified in the LGIN,HGIN,and CA groups,and the biological functions of the DEGs were analyzed via Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analyses.The CIBERSORT algorithm was used to detect the pattern of immune cell infilt-ration.Immunohistochemistry(IHC)was also conducted to validate our results.Finally,the Luminex multiplex cytokine analysis was utilized to measure the serum cytokine levels in the mice.RESULTS Compared with those in the NOR group,a total of 681541,and 840 DEGs were obtained in the LGIN,HGIN,and CA groups,respectively.Using the intersection of the three sets of DEGs,we identified 86 genes as key genes involved in the development of ESCC.Enrichment analysis revealed that these genes were enriched mainly in the keratinization,epidermal cell differentiation,and interleukin(IL)-17 signaling pathways.CIBERSORT analysis revealed that,compared with those in the NOR group,M0 and M1 macrophages in the 4NQO group showed stronger infiltration,which was validated by IHC.Serum cytokine analysis revealed that,compared with those in the NOR group,IL-1βand IL-6 were upregulated,while IL-10 was downregulated in the LGIN,HGIN,and CA groups.Moreover,the expression of the representative key genes,such as S100a8 and Krt6b,was verified in external human samples,and the results of immunohistochemical staining were consistent with the findings in mice.CONCLUSION We identified a set of key genes represented by S100a8 and Krt6b and investigated their potential biological functions.In addition,we found that macrophage infiltration and abnormal alterations in the levels of inflam-mation-associated cytokines,such as IL-1β,IL-6,and IL-10,in the peripheral blood may be closely associated with the development of ESCC.

Evaluation of genetic response of mesenchymal stem cells to nanosecond pulsed electric fields by whole transcriptome sequencing

BACKGROUND Mesenchymal stem cells(MSCs)modulated by various exogenous signals have been applied extensively in regenerative medicine research.Notably,nanosecond pulsed electric fields(nsPEFs),characterized by short duration and high strength,significantly influence cell phenotypes and regulate MSCs differentiation via multiple pathways.Consequently,we used transcriptomics to study changes in messenger RNA(mRNA),long noncoding RNA(lncRNA),microRNA(miRNA),and circular RNA expression during nsPEFs application.AIM To explore gene expression profiles and potential transcriptional regulatory mechanisms in MSCs pretreated with nsPEFs.METHODS The impact of nsPEFs on the MSCs transcriptome was investigated through whole transcriptome sequencing.MSCs were pretreated with 5-pulse nsPEFs(100 ns at 10 kV/cm,1 Hz),followed by total RNA isolation.Each transcript was normalized by fragments per kilobase per million.Fold change and difference significance were applied to screen the differentially expressed genes(DEGs).Gene Ontology and Kyoto Encyclopedia of Genes and Genomes analyses were performed to elucidate gene functions,complemented by quantitative polymerase chain reaction verification.RESULTS In total,263 DEGs were discovered,with 92 upregulated and 171 downregulated.DEGs were predominantly enriched in epithelial cell proliferation,osteoblast differentiation,mesenchymal cell differentiation,nuclear division,and wound healing.Regarding cellular components,DEGs are primarily involved in condensed chromosome,chromosomal region,actin cytoskeleton,and kinetochore.From aspect of molecular functions,DEGs are mainly involved in glycosaminoglycan binding,integrin binding,nuclear steroid receptor activity,cytoskeletal motor activity,and steroid binding.Quantitative real-time polymerase chain reaction confirmed targeted transcript regulation.CONCLUSION Our systematic investigation of the wide-ranging transcriptional pattern modulated by nsPEFs revealed the differential expression of 263 mRNAs,2 miRNAs,and 65 lncRNAs.Our study demonstrates that nsPEFs may affect stem cells through several signaling pathways,which are involved in vesicular transport,calcium ion transport,cytoskeleton,and cell differentiation.

Analyses of Chicken Tenderness Traits Based on Transcriptome Sequencing

The calpain system is ubiquitous in cells, mainly comprising calpains and calpain inhibitors, and is a widespread calcium-dependent cysteine protease in organisms that is involved in many cellular processes such as muscle degradation in vivo and affects the tenderness of meat after animal slaughter. The study found 128 DEGs that probably regulated tenderness traits were selected from 16 significantly enriched GO terms by transcriptome sequencing analysis, and found that the developmental changes in the expression levels of the CAPN1 gene in the pectoral and leg muscles were significantly positively correlated ( P <0.05) with the cumulative growth values of live weight and comb weight. The developmental changes in the expression levels of the CAST gene in the pectoral and leg muscles were not significantly correlated with the cumulative growth values of live weight and comb weight. Our results helped demonstrate the potential molecular mechanisms of tenderness in chickens and provide valuable information for chicken breeding.

Use of transcriptome sequencing to explore the effect of CSRP3 on chicken myoblasts

The mechanisms that regulate the specificity and maintenance of chicken muscle fiber types remain largely unknown. In mammals, CSRP3 has been shown to play a vital role in the maintenance of typical muscle structure and function. This study investigated the role that CSRP3 plays in chicken skeletal muscle. First, the antibody against chicken CSRP3 protein was prepared, and the expression levels of the mRNA and protein of the CSRP3 gene in four chicken skeletal muscles with different myofiber compositions were compared. Then the effects of CSRP3 silencing on the expression profile of chicken myoblast transcriptomes were analyzed. The results showed that the expression levels of the mRNA and protein of the CSRP3 gene were both associated with the composition of fiber types in chicken skeletal muscles. A total of 650 genes with at least 1.5-fold differences(Q<0.05) were identified, of which 255 genes were upregulated and 395 genes were downregulated by CSRP3 silencing. Functional enrichment showed that several pathways, including adrenergic signaling in cardiomyocytes, adipocytokine signaling pathway and apelin signaling pathway, were significantly(P<0.05) enriched both in differentially expressed genes and all expressed genes. The co-expressed gene network suggested that CSRP3 silencing caused a compensatory upregulation(Q<0.05) of genes related to the assembly of myofibrils, muscle differentiation, and contraction. Meanwhile, two fast myosin heavy chain genes(MyH1B and MyH1E)were upregulated(Q<0.05) upon CSRP3 silencing. These results suggested that CSRP3 plays a crucial role in chicken myofiber composition, and affects the distribution of chicken myofiber types, probably by regulating the expression of MyH1B and MyH1E.

Transcriptome sequencing analysis of ursolic acid-mediated proliferation suppression on cutaneous T-cell lymphoma cells

Background:Ursolic acid is a triterpenoid compound found in natural plants that exhibits antiproliferative effects in various cancer cells.Our study is the first to demonstrate the strong inhibitory effects of ursolic acid on the proliferation of cutaneous T-cell lymphoma(CTCL)cells.We aimed to further investigate the underlying mechanism of the proliferation inhibition induced by ursolic acid in CTCL cells using transcriptome sequencing.Methods:Cell counting kit-8 assays were used to observe the effects of six traditional medicine monomers on the proliferation of CTCL cells.Transcriptome sequencing was used to identify differentially expressed genes after ursolic acid treatment.Bioinformatics analysis was performed to determine the potential mechanism.Real-time quantitative PCR and western blotting analyses were performed to confirm the sequencing results and verify the possible mechanisms of ursolic acid-mediated proliferation inhibition in CTCL cells.Results:Ursolic acid exhibited the strongest inhibitory effect on the proliferation of CTCL cells among the six traditional medicine monomers.Transcriptome sequencing analysis showed that 2,466 genes were significantly altered.Combined with Kyoto Encyclopedia of Genes and Genomes functional enrichment analysis and protein-protein interaction network analysis,the interaction of various pathways and signaling molecules,such as tumor necrosis factor-α,NLR family pyrin domain containing 1,c-Jun N-terminal kinase,and melanoma differentiation-associated gene 5,accounted for the anti-tumor effects of ursolic acid in CTCL cells.Conclusion:Ursolic acid significantly inhibited the proliferation of CTCL cells,and our study laid a theoretical foundation for the future treatment of CTCL using ursolic acid.

Uncovering Small RNAs in Penicillium digitatum by Transcriptome Sequencing

Small RNAs in Penicillium digitatum were identified and analyzed via transcriptome sequencing on the BGISEQ-500 platform. A total of 15 predicted miRNAs and 10718 novel siRNAs were found. Their length distribution, sequence, predicted construction, base bias, expression levels and potential targets were determined as well. Through pathway and KEGG enrichment analysis, the miRNA target genes were mostly involved in carbohydrate metabolism, transport and catabolism, translation and amino acid metabolism. The target genes involved in aflatoxin biosynthesis and proteasome had a higher rich factor value. The results will provide a theoretical foundation for understanding the developmental and pathogenic mechanisms of P. digitatum at the transcriptional level.

Transcriptome Sequencing and Analysis of Agaricus blazei Mycelia Under Cadmium Stress

To explore the mechanism of Agaricus blazei Murrill enrichment of the heavy metal cadmium,we employed Illumina high-throughput sequencing to analyze the transcriptomes of A.blazei mycelia treated with and without exogenous cadmium addition,and then the differentially expressed genes(DEGs)of the strains with high and low cadmium enrichment between the control and cadmium treatment were screened out.The results showed that the DEGs were mainly involved in steroid biosynthesis,antibiotic biosynthesis,protein processing in endoplasmic reticulum,glutathione metabolism and other pathways.Carbon metabolism and glutathione metabolism may play an important role in the response of A.blazei mycelium to cadmium stress.

Transcriptome sequencing of essential marine brown and red algal species in China and its significance in algal biology and phylogeny

Most phaeophytes （brown algae） and rhodophytes （red algae） dwell exclusively in marine habitats and play important roles in marine ecology and biodiversity. Many of these brown and red algae are also important resources for industries such as food, medicine and materials due to their unique metabolisms and me-tabolites. However, many fundamental questions surrounding their origins, early diversification, taxonomy, and special metabolisms remain unsolved because of poor molecular bases in brown and red algal study. As part of the 1 000 Plant Project, the marine macroalgal transcriptomes of 19 Phaeophyceae species and 21 Rhodophyta species from China＇s coast were sequenced, covering a total of 2 phyla, 3 classes, 11 orders, and 19 families. An average of 2 Gb per sample and a total 87.3 Gb of RNA-seq raw data were generated. Approxi-mately 15 000 to 25 000 unigenes for each brown algal sample and 5 000 to 10 000 unigenes for each red algal sample were annotated and analyzed. The annotation results showed obvious differences in gene expres-sion and genome characteristics between red algae and brown algae;these differences could even be seen between multicellular and unicellular red algae. The results elucidate some fundamental questions about the phylogenetic taxonomy within phaeophytes and rhodophytes, and also reveal many novel metabolic pathways. These pathways include algal CO2 fixation and particular carbohydrate metabolisms, and related gene/gene family characteristics and evolution in brown and red algae. These findings build on known algal genetic information and significantly improve our understanding of algal biology, biodiversity, evolution, and potential utilization of these marine algae.

Characterization of Chiton Ischnochiton hakodadensis Foot Based on Transcriptome Sequencing

Chiton(Ischnochiton hakodadensis) is one of marine mollusks well known for its eight separate shell plates. I. hakodadensis is important, which plays a vital role in the ecosystems it inhabits. So far, the genetic studies on the chiton are scarce due in part to insufficient genomic resources available for this species. In this study, we investigated the transcriptome of the chiton foot using Illumina sequencing technology. The reads were assembled and clustered into 256461 unigenes, of which 42247 were divided into diverse functional categories by Gene Ontology(GO) annotation terms, and 17256 mapped onto 365 pathways by KEGG pathway mapping. Meanwhile, a set of differentially expressed genes(DEGs) between distal and proximal muscles were identified as the foot adhesive locomotion associated, thus were useful for our future studies. Moreover, up to 679384 high-quality single nucleotide polymorphisms(SNPs) and 19814 simple sequence repeats(SSRs) were identified in this study, which are valuable for subsequent studies on genetic diversity and variation. The transcriptomic resource obtained in this study should aid to future genetic and genomic studies of chiton.

Comparative analysis on transcriptome sequencings of six Sargassum species in China

Species of Sargassum are distributed worldwide, and are of great ecological and economic importance in marine ecosystems and bioresources. In this study, transcriptome sequencings of six Sargassum species were performed for the first time using an Illumina platform. For each sample, a total of 2.1-2.5 Gb of nucle-otides are collected and assembled into 69 871-116 790 scaffolds, with an average length of 410-550 bp and N50 length of 756-1 462 bp. A total of 20 512-28 684 unigenes of each sample were annotated and compared well with known gene sequences from nr database. Clusters of Orthologous Groups （COG）, gene ontology （GO） and Kyoto Encyclopedia of Genes and Genomes （KEGG） analyses were also performed for further un-derstanding of gene functions and regulation pathways. Gene expression levels were calculated based on RPKM values and compared among these species, especially for those genes related to carbohydrate metab-olism. Cluster analyses indicated that the differences of global gene expression between S. fusiforme, which was nominated as Hizikia fusiformis before, and other five species were not significant. Further phylogenet-ic analysis of 108 orthologous genes confirmed that S. fusiforme had closer relationship with S. hemiphyllum rather than S. horneri. These transcriptome data provided valuable information for better understanding of genome and gene characteristics of Sargassum algae and benefiting comparative and phylogenetic studies of Phaeophyceae species in future studies.

Transcriptome sequencing analysis of lncRNA expression in peripheral blood mononuclear cells from patients with COPD

Objective:To screen abnormally expressed lncRNAs in peripheral blood mononuclear cells from patients with COPD.Methods:The peripheral blood of 3 COPD patients and 3 normal controls were collected from our hospital,mononuclear cells were isolated,RNA was extracted and then transcriptome sequencing was performed.The expression difference between the two groups of samples was calculated based on p<0.05 and|log2FC|>1.Plot the difference lncRNA heat map and volcano map.The Lncpro database may predict mRNAs regulated by differential lncRNA,and perform the GO function and KEGG signaling clustering.Results:There were 67 lncRNAs between the COPD group and the control group that met the difference of p<0.05 and|log2FC|>1,of which 33 were up-regulated and 34 were down-regulated.Between the two groups.The target genes are mainly enriched in GO functions:regulatory functions of multicellular biological processes,regulatory functions of development processes,structured morphogenesis functions,system development functions,and development process functions.Target genes are mainly enriched in KEGG signaling pathways:multi-species apoptotic pathway,TGF-βsignaling pathway,complement and coagulation cascade pathway,colorectal cancer pathway and apoptosis pathway.Conclusion:Our results provide general information and possible regulatory functions and pathways of lncRNA expression changes in peripheral blood mononuclear cells of COPD,which may help clarify the underlying mechanism of COPD.

Screening of morphology-related genes based on predatorinduced transcriptome sequencing and the functional analysis of Dagcut gene in Daphnia galeata

High fish predation pressure can trigger"induced defense"in Daphnia species,resulting in phenotypic plasticity in morphology,behavior,or life-history traits.The molecular mechanisms of defense morphogenesis(e.g.,the tail spine and helmet)in Daphnia remain unclear.In the pres-ent study,the tail spine,helmet,and body of Daphnia galeata under fish and non-fish kairomones conditions were collected for transcriptome analysis.A total of 24 candidate genes related to the morphological defense of D.galeata were identified,including 2 trypsin,one cuticle protein,1 C1qDC protein,and 2 ferritin genes.The function of the Dagcut gene(D.galeata cuticle protein gene)in relation to tail spine morphology was assessed using RNA interference(RNAi).Compared with the EGFP(Enhanced green fluorescent protein)treatment,after RNAi,the expression levels of the Dagcut gene(D.galeata cuticle protein gene)showed a significant decrease.Correspondingly,the tail spines of the offspring pro-duced by D.galeata after RNAi of the Dagcut gene appeared curved during the experiment.In whole-mount in situ hybridization,a clear signal site was detected on the tail spine of D.galeata before RNAi which disappeared after RNAi.Our results suggest that the Dagcut gene may play an important role in tail spine formation of D.galeata,and will provide a theoretical basis for studying the molecular mechanisms of the morpho-logical plasticityin cladocera inthefuture.

Study on compound Duzhong Jiangu granule in the treatment of knee osteoarthritis based on transcriptome analysis strategy

Duzhong Jiangu Granule in the treatment of primary osteoarthritis(POA)model of postmenopausal kidney deficiency type in Hartley female guinea pigs after ovariectomy,and the correlation between gene expression of bone marrow tissue,cartilage tissue,and knee osteoarthritis.Methods:383-months-old Hartley female guinea pigs after one week of adaptive feeding were weighed about 400 g±20 g,numbered,and sorted by ear tag.Six of them were selected as normal groups by looking up random number tables.The remainder were removed from the bilateral ovaries to construct the postmenopausal kidney deficiency model,and castrated guinea pigs were used to construct the postmenopausal kidney deficiency POA model.After the modeling cycle,a guinea pig from the blank group and a guinea pig from the model group were sacrificed and the right knee was observed.The model was established and the experiments continued.There were five guinea pigs in the blank group,and the remaining model guinea pigs were randomly divided into model control group,high-dose group of compound Duzhong Jiangu granules,middle-dose group of compound Duzhong Jiangu granules,low-dose group of compound Duzhong Jiangu granules and design group,with 5 guinea pigs in each group.Blanks and model groups were given a cellulose sodium solution by gavage.The guinea pigs were sacrificed after 30 days of intragastric administration.The left knee cartilage and bone marrow of the blank group,model group,middle dose group,and high dose group of compound Duzhong Jiangu granule were collected and applied to transcriptome sequencing,and the sequencing data were analyzed,including differential gene expression analysis,functional enrichment analysis of database established by Gene Ontology federation(GO)and Kyoto Encyclopedia of Gene and Genome(KEGG)pathway enrichment analysis.The complete specimens of the right knee joint were collected,and the morphological changes of the cartilage of the right knee joint in each group were observed by saffron rapid green staining,and the subchondral bone was quantitatively analyzed by Micro CT so that the expression of TRAF6,MIP-1βand IL-1βprotein in NF-kappa B signaling pathway was detected by Western Blot technique(WB).Results:The results of Safranin Fast Green staining showed that Compound Duzhong Jiangu Granules could effectively reduce the degree of morphological damage of articular cartilage in guinea pigs with the POA model.According to the analysis results of the subchondral bone structure under Micro CT,Compound Duzhong Jiangu Granules can improve the bone condition of the POA model,thus delaying the process of degenerative changes of the knee joint.From the results of transcriptome analysis,Compound Duzhong Jiangu Granules can inhibit the expression of related genes in POA model guinea pigs.According to the results of Wester Bolt verification,Compound Duzhong Jiangu Granules can effectively improve knee osteoarthritis.Conclusions:The effect of Compound Duzhong Jiangu Granules on OA is obvious,and its mechanism may be related to the expression of genes GZMK,Jchain,igkc,IGHV3-74,IGHV3-11,IGHV4-1,CCL5,and IGKV1–39.

Comprehensive profiling of EBV gene expression in nasopharyn- geal carcinoma through paired-end transcriptome sequencing

The latent expression pattern of Epstein-Barr Virus （EBV）igenes in nasopharyngeal carcinoma （NPC） has been extensively investigated, and the expression of several lytic genes in NPC has been reported. However, comprehensive information through EBV transcriptome analysis in NPC is limited. We performed paired-end RNA-seq to systematically and comprehensively characterize the expression of EBV genes in NPC tissue and C666-1 NPC cell line, which consistently carries EBV. In addition to the transcripts restricted to type II latency infection, the type Ⅲ latency EBNA3s genes and a substantial number of lytic genes, such as BZLF1, BRLF1, and BMRF1, were detected through RNA-seq and were further verified in C666-1 cells and NPC tissue through real- time PCR. We also performed clustering analysis to classify NPC patient groups in terms of EBV gene expression, which presented two subtypes of NPC samples. Results revealed interesting patterns of EBV gene expression in NPC patients. This clustering was correlated with many signaling pathways, such as those related to heterotrimeric G-protein signaling, inflammation mediated by chemokine and cytokine signaling, ribosomes, protein metabolism, influenza infection, and ECM-receptor interaction. Our combined findings suggested that the expression of EBV genes in NPC is restricted not only to type II latency genes but also to type Ⅲ latency and lyric genes. This study provided further insights into the potential role of EBV in the development of NPC.

Transcriptome Sequencing and de novo Analysis for Oviductus Ranae of Rana chensinensis Using Illumina RNA-Seq Technology

Oviductus Ranae is the dried oviduct of female Rana tem-poraria chensinensis （David）, distributed mainly in North- eastern China. Oviductus Ranae is one of the best-known and highly valued oriental foods and medicines. Traditional Chinese medicine holds that Oviductus Ranae can nourish yin, moisten lung and replenish the kidney essence. Meanwhile, activities of Oviductus Ranae such as anti-aging, anti-lipemic, anti-oxidation and anti-fatigue have also been demonstrated by modern phar-macological studies. Previous studies have shown that Oviductus Ranae is mainly composed of proteins, which are up to 50% or more.

Integrative Analyses of Lung Squamous Cell Carcinoma in Ten Chinese Patients with Transcriptome Sequencing

Few effective therapies have been developed for the treatment of lung squamous cell carcinoma （SQCC）, in part due to a lack of un- derstanding regarding the mechanisms underlying the initiation and development of this disease. Whole transcriptome sequencing not only provides insight into the expression of all transcribed genes, but offers an efficient approach for identifying genetic variations, including gene fusions, mutations and alternative splicing. In this study, we performed whole transcriptome sequencing of 10 patients with stage IIIA lung SQCC, and discovered a large number of single nucleotide variants （SNVs： mean of 12.2 SNVs/Mb）, with C〉T/G〉A and A〉G/T〉C transitions being the most frequently observed. Additionally, a total of 132 gene fusions were identified based upon TopHat alignments, 70.5% （93/132） of which occurred as a result of intra-chromosomal rearrangements. Based on the number of supporting reads for each fusion, we further validated 20 of the 26 top gene fusions by RT-PCR and Sanger sequencing. Taken together, these data provide an in-depth view of transcriptional alterations in lung SQCC patients, and may be useful for identification of new therapeutic targets.

Integrated transcriptome sequencing and RNA interference reveals molecular changes in Diaphorina citri after exposure to validamycin

Validamycin has been widely used as a specific competitive inhibitor of tre-halase.In our previous research,validamycin significantly inhibited trehalase activity and chitin synthesis in Diaphorina citri,resulting in abnormal phenotypes.However,the mechanism of validamycin's action on D.citri remains unclear.Here,using a comparative transcriptome analysis,464 differentially expressed genes(DEGs)in D.citri were identified after validamycin treatment.A Gene Ontology enrichment analysis revealed that these DEGs were mainly involved in“small molecule process”,“structural molecule ac-tivity”and“transition metal ion binding”.DEGs involved in chitin metabolism,cuticle synthesis and insecticide detoxification were validated by reverse transcription quantitative polymerase chain reaction.The RNA interference of D.citri chitinase-like protein EN03 and D.citri cuticle protein 7 genes significantly affected D.citri molting.Moreover,the recombinant chitinase-like protein EN03 exhibited a chitin-binding property,and an antimicrobial activity against Bacillus subtilis.This study provides a first insight into the molecular changes in D.citri after exposure to validamycin and identifies two effective RNA interference targets for D.citri control.

Transcriptome-based analysis of key genes and pathways affecting the linoleic acid content in chickens

Linoleic acid is an essential polyunsaturated fatty acid that cannot be synthesized by humans or animals themselves and can only be obtained externally.The amount of linoleic acid present has an impact on the quality and flavour of meat and indirectly affects consumer preference.However,the molecular mechanisms influencing the deposition of linoleic acid in organisms are not clear.As the molecular mechanisms of linoleic acid deposition are not well understood,to investigate the main effector genes affecting the linoleic acid content,this study aimed to screen for hub genes in slow-type yellow-feathered chickens by transcriptome sequencing(RNA-Seq)and weighted gene coexpression network analysis(WGCNA).We screened for candidate genes associated with the linoleic acid content in slow-type yellow-feathered broilers.A total of 399 Tiannong partridge chickens were slaughtered at 126 days of age,fatty acid levels were measured in pectoral muscle,and pectoral muscle tissue was collected for transcriptome sequencing.Transcriptome sequencing results were combined with phenotypes for WGCNA to screen for candidate genes.KEGG enrichment analysis was also performed on the genes that were significantly enriched in the modules with the highest correlation.A total of 13310 genes were identified after quality control of transcriptomic data from 399 pectoral muscle tissues.WGCNA was performed,and a total of 26 modules were obtained,eight of which were highly correlated with the linoleic acid content.Four key genes,namely,MDH2,ATP5B,RPL7A and PDGFRA,were screened according to the criteria|GS|>0.2 and|MM|>0.8.The functional enrichment results showed that the genes within the target modules were mainly enriched in metabolic pathways.In this study,a large-sample-size transcriptome analysis revealed that metabolic pathways play an important role in the regulation of the linoleic acid content in Tiannong partridge chickens,and MDH2,ATP5B,RPL7A and PDGFRA were screened as important candidate genes affecting the linoleic acid content.The results of this study provide a theoretical basis for selecting molecular markers and comprehensively understanding the molecular mechanism affecting the linoleic acid content in muscle,providing an important reference for the breeding of slow-type yellowfeathered broiler chickens.

The transcriptome analysis of cleft lip/palate-related PTCH1 variants in GMSM-K cells show carcinogenic potential

Cancer progression involves the sonic hedgehog(SHH)pathway,in which the receptor PTCH1 actives the downstream pathways.Dysfunction of PTCH1 can lead to nevoid basal cell carcinoma Syndrome(NBCCs)including neoplastic disease and congenital disorder.To evaluate the relationship between PTCH1 and cancer,we applied the CRISPR/Cas9 system to knock out PTCH1 in oral nontumorous epithelial cells(GMSM-K).Then we screened six PTCH1 variants associated with cleft lip/palate(CL/P),one of the congenital disorders in NBCCs,and generated PTCH1 variant and wild-type recombinant PTCH1^(−/−)GMSM-K cell lines.Transcriptome sequencing was conducted in these cell lines.The results revealed that differentially expressed genes(DEGs)in PTCH1^(−/−)GMSM-K were enriched in extracellular compartments,contributing epithelial diseases by pathway enrichment analysis.RT-PCR confirmed that KRT34,KRT81,KRT86,PDGFB,and WNT10B genes,associated with extracellular compartments were highly expressed in PTCH1^(−/−).The Kyoto Encyclopedia of Genes and Genomes analysis also suggested that DEGs are closely related to focal adhesion,transcriptional misregulation,and proteoglycans in breast and gastric cancers.Comparative analysis of samples revealed that the CL/P-associated PTCH1 variants A443G and V908G are potentially carcinogenic.These findings provide new insights into the carcinogenic potential of PTCH1 dysfunction.

Transcriptome Analysis of a Wild Eggplant Germplasm M239 in Response to Verticillium dahliae Infection

In this study,wild eggplant germplasm No.M239,which is highly susceptible to Verticillium wilt,was used as the experimental material.The physiological and biochemical indices(SOD,PAL,MDA and soluble protein)of M239 roots were measured at different times(0,12,24,36,48,60 and 72 h)post inoculation with Verticillium dahliae,and the key time points for the M239 response to Verticillium wilt infection were screened.Then,RNA-Seq technology was used to screen the differentially expressed genes(DEGs)in M239 roots at 0,12 and 48 h post-inoculation(hpi).The transcriptional results of M239 were also compared with those resistance genes from some reported wild relative Solanum species(S.sisymbriifolium and S.aculeatissimum).Then some DEGs were chosen for validation by qRT–PCR.The results showed that 12 and 48 hpi were the turning points in the changes in all physiological and biochemical indices.A total of 6,783 DEGs were identified by RNA-Seq,including 6,141 DEGs(3,046 upregulated and 3,095 downregulated)at the M_12 h vs.M_0 h,1,903 DEGs(792 upregulated and 1,111 downregulated)at M_48 h vs.M_12 h,and 1,261 DEGs that appeared simultaneously in both stages.KEGG enrichment analysis showed that there were 5 metabolic pathways enriched from DEGs,which were mostly related to primary metabolism,such as glycolysis,amino acid and ribosome biogenesis.Compared with the NCBI non-redundant protein(NR)database,one Ve2 homologous gene and 8 PR protein-related genes were screened.Transcription factor analysis showed that there were a large number of DEGs,such as MYB,AP2-EREBP,bHLH,NAC and Orphans in the two stages.Compared with the reported Verticillium wilt-resistant wild eggplant species,it was found that there were fewer genes and enriched metabolic pathways in the M239 response to Verticillium wilt infection,and it also lacked the response of some known key resistance genes.These results proved that the above resistance genes and metabolic pathways played a key role in the wild eggplant response to V.dahliae infection.