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.

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.

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.

Recent developments in application of single-cell RNA sequencing in the tumour immune microenvironment and cancer therapy

The advent of single-cell RNA sequencing(scRNA-seq)has provided insight into the tumour immune microenvironment(TIME).This review focuses on the application of scRNA-seq in investigation of the TIME.Over time,scRNA-seq methods have evolved,and components of the TIME have been deciphered with high resolution.In this review,we first introduced the principle of scRNA-seq and compared different sequencing approaches.Novel cell types in the TIME,a continuous transitional state,and mutual intercommunication among TIME components present potential targets for prognosis prediction and treatment in cancer.Thus,we concluded novel cell clusters of cancerassociated fibroblasts(CAFs),T cells,tumour-associated macrophages(TAMs)and dendritic cells(DCs)discovered after the application of scRNA-seq in TIME.We also proposed the development of TAMs and exhausted T cells,as well as the possible targets to interrupt the process.In addition,the therapeutic interventions based on cellular interactions in TIME were also summarized.For decades,quantification of the TIME components has been adopted in clinical practice to predict patient survival and response to therapy and is expected to play an important role in the precise treatment of cancer.Summarizing the current findings,we believe that advances in technology and wide application of single-cell analysis can lead to the discovery of novel perspectives on cancer therapy,which can subsequently be implemented in the clinic.Finally,we propose some future directions in the field of TIME studies that can be aided by scRNA-seq technology.

Single-cell RNA sequencing in cornea research:Insights into limbal stem cells and their niche regulation

The corneal epithelium is composed of stratified squamous epithelial cells on the outer surface of the eye,which acts as a protective barrier and is critical for clear and stable vision.Its continuous renewal or wound healing depends on the proliferation and differentiation of limbal stem cells(LSCs),a cell population that resides at the limbus in a highly regulated niche.Dysfunction of LSCs or their niche can cause limbal stem cell deficiency,a disease that is manifested by failed epithelial wound healing or even blindness.Nevertheless,compared to stem cells in other tissues,little is known about the LSCs and their niche.With the advent of single-cell RNA sequencing,our understanding of LSC characteristics and their microenvironment has grown considerably.In this review,we summarized the current findings from single-cell studies in the field of cornea research and focused on important advancements driven by this technology,including the heterogeneity of the LSC population,novel LSC markers and regulation of the LSC niche,which will provide a reference for clinical issues such as corneal epithelial wound healing,ocular surface reconstruction and interventions for related diseases.

Single-nuclei RNA sequencing uncovers heterogenous transcriptional signatures in Parkinson's disease associated with nuclear receptor-related factor 1 defect

Previous studies have found that deficiency in nuclear receptor-related factor 1(Nurr1),which participates in the development,differentiation,survival,and degeneration of dopaminergic neurons,is associated with Parkinson s disease,but the mechanism of action is perplexing.Here,we first asce rtained the repercussion of knocking down Nurr1 by pe rforming liquid chromatography coupled with tandem mass spectrometry.We found that 231 genes were highly expressed in dopaminergic neurons with Nurr1 deficiency,14 of which were linked to the Parkinson’s disease pathway based on Kyoto Encyclopedia of Genes and Genomes analysis.To better understand how Nurr1 deficiency autonomously invokes the decline of dopaminergic neurons and elicits Parkinson’s disease symptoms,we performed single-nuclei RNA sequencing in a Nurr1 LV-shRNA mouse model.The results revealed cellular heterogeneity in the substantia nigra and a number of activated genes,the preponderance of which encode components of the major histocompatibility Ⅱ complex.Cd74,H2-Ab1,H2-Aα,H2-Eb1,Lyz2,Mrc1,Slc6α3,Slc47α1,Ms4α4b,and Ptprc2 were the top 10 diffe rentially expressed genes.Immunofluorescence staining showed that,after Nurr1knockdown,the number of CD74-immunoreactive cells in mouse brain tissue was markedly increased.In addition,Cd74 expression was increased in a mouse model of Parkinson’s disease induced by treatment with 6-hydroxydopamine.Ta ken togethe r,our res ults suggest that Nurr1 deficiency results in an increase in Cd74 expression,thereby leading to the destruction of dopaminergic neuro ns.These findings provide a potential therapeutic target for the treatment of Parkinson’s disease.

Application of single-cell RNA sequencing in head and neck squamous cell carcinoma

Single-cell RNA sequencing has been broadly applied to head and neck squamous cell carcinoma(HNSCC) for characterizing the heterogeneity and genomic mutations of HNSCC benefiting from the advantage of single-cell resolution. We summarized most of the current studies and aimed to explore their research methods and ideas, as well as how to transform them into clinical applications. Through single-cell RNA sequencing, we found the differences in tumor cells’ expression programs and differentiation tracks. The studies of immune microenvironment allowed us to distinguish immune cell subpopulations, the extensive expression of immune checkpoints, and the complex crosstalk network between immune cells and non-immune cells. For cancerassociated fibroblasts(CAFs), single-cell RNA sequencing had made an irreplaceable contribution to the exploration of their differentiation status, specific CAFs markers, and the interaction with tumor cells and immune cells. In addition, we demonstrated in detail how single-cell RNA sequencing explored the HNSCC epithelial-tomesenchymal transition(EMT) model and the mechanism of drug resistance, as well as its clinical value.

Single-cell RNA sequencing reveals the dynamics of hepatic non-parenchymal cells in autoprotection against acetaminophen-induced hepatotoxicity

Gaining a better understanding of autoprotection against drug-induced liver injury(DILI)may provide new strategies for its prevention and therapy.However,little is known about the underlying mechanisms of this phenomenon.We used single-cell RNA sequencing to characterize the dynamics and functions of hepatic non-parenchymal cells(NPCs)in autoprotection against DILI,using acetaminophen(APAP)as a model drug.Autoprotection was modeled through pretreatment with a mildly hepatotoxic dose of APAP in mice,followed by a higher dose in a secondary challenge.NPC subsets and dynamic changes were identified in the APAP(hepatotoxicity-sensitive)and APAP-resistant(hepatotoxicity-resistant)groups.A chemokine(C-C motif)ligand 2^(+)endothelial cell subset almost disappeared in the APAP-resistant group,and an R-spondin 3^(+)endothelial cell subset promoted hepatocyte proliferation and played an important role in APAP autoprotection.Moreover,the dendritic cell subset DC-3 may protect the liver from APAP hepatotoxicity by inducing low reactivity and suppressing the autoimmune response and occurrence of inflammation.DC-3 cells also promoted angiogenesis through crosstalk with endothelial cells via vascular endothelial growth factor-associated ligand-receptor pairs and facilitated liver tissue repair in the APAP-resistant group.In addition,the natural killer cell subsets NK-3 and NK-4 and the Sca-1^(-)CD62L^(+)natural killer T cell subset may promote autoprotection through interferon-γ-dependent pathways.Furthermore,macrophage and neutrophil subpopulations with anti-inflammatory phenotypes promoted tolerance to APAP hepatotoxicity.Overall,this study reveals the dynamics of NPCs in the resistance to APAP hepatotoxicity and provides novel insights into the mechanism of autoprotection against DILI at a high resolution.

Single cell RNA sequencing reveals mesenchymal heterogeneity and critical functions of Cd271 in tooth development

BACKGROUND Accumulating evidence suggests that the maxillary process,to which cranial crest cells migrate,is essential to tooth development.Emerging studies indicate that Cd271 plays an essential role in odontogenesis.However,the underlying mechanisms have yet to be elucidated.AIM To establish the functionally heterogeneous population in the maxillary process,elucidate the effects of Cd271 deficiency on gene expression differences.METHODS p75NTR knockout(Cd271-/-)mice(from American Jackson laboratory)were used to collect the maxillofacial process tissue of p75NTR knockout mice,and the wildtype maxillofacial process of the same pregnant mouse wild was used as control.After single cell suspension,the cDNA was prepared by loading the single cell suspension into the 10x Genomics Chromium system to be sequenced by NovaSeq6000 sequencing system.Finally,the sequencing data in Fastq format were obtained.The FastQC software is used to evaluate the quality of data and CellRanger analyzed the data.The gene expression matrix is read by R software,and Seurat is used to control and standardize the data,reduce the dimension and cluster.We search for marker genes for subgroup annotation by consulting literature and database;explore the effect of p75NTR knockout on mesenchymal stem cells(MSCs)gene expression and cell proportion by cell subgrouping,differential gene analysis,enrichment analysis and protein-protein interaction network analysis;understand the interaction between MSCs cells and the differentiation trajectory and gene change characteristics of p75NTR knockout MSCs by cell communication analysis and pseudo-time analysis.Last we verified the findings single cell sequencing in vitro.RESULTS We identified 21 cell clusters,and we re-clustered these into three subclusters.Importantly,we revealed the cell–cell communication networks between clusters.We clarified that Cd271 was significantly associated with the regulation of mineralization.CONCLUSION This study provides comprehensive mechanistic insights into the maxillary-process-derived MSCs and demonstrates that Cd271 is significantly associated with the odontogenesis in mesenchymal populations.

Differential mRNA expression profiling of oral squamous cell carcinoma by high-throughput RNA sequencing

Differentially expressed genes are thought to regulate the development and progression of oral squamous cell carcinomas （OSCC）. The purpose of this study was to screen differentially expressed mRNAs in OSCC and matched paraneoplastic normal tissues, and to explore the intrinsic mechanism of OSCC development and progres- sion. We obtained the differentially expressed mRNA expression profiles in 10 pairs of fresh-frozen OSCC tissue specimens and matched paraneoplastic normal tissue specimens by high-throughput RNA sequencing. By using Gene Ontology enrichment analysis and Kyoto Encyclopedia of Genes and Genomes （KEGG） pathway enrichment analyses, the functional significance of the differentially expressed genes were analyzed. We identified 1,120 sig- nificantly up-regulated mRNAs and 178 significantly down-regulated mRNAs in OSCC, compared to normal tissue. The differentially expressed mRNAs were involved in 20 biological processes and 68 signal pathways. Compared to adjacent normal tissue, the expression of MAGEAll was up-regulated; TCHH was down-regulated. These find- ings were verified by real-time PCR. These differentially expressed mRNAs may function as oncogenes or tumor suppressors in the development and progression of OSCC. This study provides novel insights into OSCC. However, further work is needed to determine if these differentially expressed mRNAs have potential roles as diagnostic bio- markers and candidate therapeutic targets for OSCC.

Genome-wide identification and comparative analysis of drought related genes in roots of two maize inbred lines with contrasting drought tolerance by RNA sequencing

Drought is one of the most important abiotic stresses affecting maize growth and development and therefore resulting in yield loss.Thus it is essential to understand molecular mechanisms of drought stress responses in maize for drought tolerance improvement.The root plays a critical role in plants sensing water deficit.In the present study,two maize inbred lines,H082183,a drought-tolerant line,and Lv28,a drought-sensitive line,were grown in the field and treated with different water conditions(moderate drought,severe drought,and well-watered conditions)during vegetative stage.The transcriptomes of their roots were investigated by RNA sequencing.There were 1428 and 512 drought-responsive genes(DRGs)in Lv28,688 and 3363 DRGs in H082183 under moderate drought and severe drought,respectively.A total of 31 Gene Ontology(GO)terms were significantly over-represented in the two lines,13 of which were enriched only in the DRGs of H082183.Based on results of Kyoto encyclopedia of genes and genomes(KEGG)enrichment analysis,"plant hormone signal transduction"and"starch and sucrose metabolism"were enriched in both of the two lines,while"phenylpropanoid biosynthesis"was only enriched in H082183.Further analysis revealed the different expression patterns of genes related to abscisic acid(ABA)signal pathway,trehalose biosynthesis,reactive oxygen scavenging,and transcription factors might contribute to drought tolerance in maize.Our results contribute to illustrating drought-responsive molecular mechanisms and providing gene resources for maize drought improvement.

Recent advances of single-cell RNA sequencing technology in mesenchymal stem cell research

Mesenchymal stem cells(MSCs)are multipotent stromal cells with great potential for clinical applications.However,little is known about their cell heterogeneity at a single-cell resolution,which severely impedes the development of MSC therapy.In this review,we focus on advances in the identification of novel surface markers and functional subpopulations of MSCs made by single-cell RNA sequencing and discuss their participation in the pathophysiology of stem cells and related diseases.The challenges and future directions of single-cell RNA sequencing in MSCs are also addressed in this review.

Panoramic comparison between NK cells in healthy and cancerous liver through single-cell RNA sequencing

Objective:NK cells play crucial roles in the immune defense mechanisms against viral infections and transformed cells.However,the developmental progression,transcriptomic landscape,and functional subtypes of liver NK cells are not well defined.Hepatocellular carcinoma(HCC)accounts for approximately 80%of primary liver cancer worldwide,yet the biological characteristics of NK cells in the HCC environment are unclear.Therefore,we aimed to determine these cells’roles in tumorigenesis and prognosis.Methods:We compared the single-cell RNA sequencing profiles of NK cells purified from blood(n=1),healthy liver tissues(n=3),HCC tumor tissues(n=4),and peritumor liver tissues(n=1)to identify NK cell subsets.Furthermore,we performed bioinformatics analysis by using The Cancer Genome Atlas(TCGA)data to identify prognostic biomarkers simultaneously overexpressed in the blood and tumor tissues of patients with HCC.Results:Transcriptomic analysis revealed 5 NK cell subsets(L1-NK-CD56bright,L2-NK-CD56dim,L3-NK-HLA,L4-LrNK-FCGR3A,and L5-LrNK-XCL1)in the healthy liver tissues.However,the transitional L3 subset and the CXCR6+CD16+L4 subset with strong anti-tumor activity were absent in the HCC and peritumor liver tissues.Furthermore,4 common prognosis-associated genes(RHOB,TALDO1,HLA-DPA1,and TKT)were significantly overexpressed in the paired tumor tissue and blood.Conclusions:Our study revealed 5 specific subsets of NK cells in healthy human liver tissues.However,only 3 of the 5 NK cell subsets were present in HCC and peritumor tissues.The cytotoxic NK cell subsets were absent in HCC tissues.Furthermore,we identified 4 potential non-invasive prognostic biomarkers in patients with HCC.

Prudent application of single-cell RNA sequencing in understanding cellular features and functional phenotypes in cancer studies

This decade has seen remarkable advances in the field of high-throughput single cell techniques.Single-cell RNA sequencing(sc RNA-seq)has proven to be a powerful strategy to study the heterogeneity in clinical samples,providing an unbiased approach to uncover the characteristics in different cell subsets.To ensure the reproducibility and robustness of biological discoveries,researchers need to be aware of hidden caveats in tissue dissociation,cell capturing and transcripts measurement which may affect cell composition assessment and cellular function annotation.With measured interpretation of data and innovations in experimental and technical approaches,sc RNA-seq can greatly unravel the heterogeneity in complex system and improve our understandings in tissue homeostasis and cancer biology.

Landscape of cell heterogeneity and evolutionary trajectory in ulcerative colitis-associated colon cancer revealed by single-cell RNA sequencing

Objective:The goal of this study was to get preliminary insight on the intra-tumor heterogeneity in colitisassociated cancer(CAC)and to reveal a potential evolutionary trajectory from ulcerative colitis(UC)to CAC at the single-cell level.Methods:Fresh samples of tumor tissues and adjacent UC tissues from a CAC patient with pT3N1M0 stage cancer were examined by single-cell RNA sequencing(scRNA-seq).Data from The Cancer Genome Atlas(TCGA)and The Human Protein Atlas were used to confirm the different expression levels in normal and tumor tissues and to determine their relationships with patient prognosis.Results:Ultimately,4,777 single-cell transcriptomes(1,220 genes per cell)were examined,of which 2,250(47%)and 2,527(53%)originated from tumor and adjacent UC tissues,respectively.We defined the composition of cancer-associated stromal cells and identified six cell clusters,including myeloid,T and B cells,fibroblasts,endothelial and epithelial cells.Notable pathways and transcription factors involved in these cell clusters were analyzed and described.Moreover,the precise cellular composition and developmental trajectory from UC to UCassociated colon cancer were graphed,and it was predicted that CD74,CLCA1,and DPEP1 played a potential role in disease progression.Conclusions:scRNA-seq technology revealed intra-tumor cell heterogeneity in UC-associated colon cancer,and might provide a promising direction to identify novel potential therapeutic targets in the evolution from UC to CAC.

Characteristics of alveolar macrophages in bronchioalveolar lavage fluids from active tuberculosis patients identified by single-cell RNA sequencing

Tuberculosis(TB),is an infectious disease caused by Mycobacterium tuberculosis(M.tuberculosis),and presents with high morbidity and mortality.Alveolar macrophages play an important role in TB pathogenesis although there is heterogeneity and functional plasticity.This study aimed to show the characteristics of alveolar macrophages from bronchioalveolar lavage fluid(BALF)in active TB patients.Single-cell RNA sequencing(scRNA-seq)was performed on BALF cells from three patients with active TB and additional scRNA-seq data from three healthy adults were established as controls.Transcriptional profiles were analyzed and compared by differential gene expression and functional enrichment analysis.We applied pseudo-temporal trajectory analysis to investigate correlations and heterogeneity within alveolar macrophage subclusters.Alveolar macrophages from active TB patients at the single-cell resolution are described.We found that TB patients have higher cellular percentages in five macrophage subclusters.Alveolar macrophage subclusters with increased percentages were involved in inflammatory signaling pathways as well as the basic macrophage functions.The TB-increased alveolar macrophage subclusters might be derived from M1-like polarization state,before switching to an M2-like polarization state with the development of M.tuberculosis infection.Cell-cell communications of alveolar macrophages also increased and enhanced in active TB patients.Overall,our study demonstrated the characteristics of alveolar macrophages from BALF in active TB patients by using scRNA-seq.

Single-cell RNA sequencing of mitoticarrested prospermatogonia with DAZL::GFP chickens and revealing unique epigenetic reprogramming of chickens

Background:Germ cell mitotic arrest is conserved in many vertebrates,including birds,although the time of entry or exit into quiescence phase differs.Mitotic arrest is essential for the normal differentiation of male germ cells into spermatogonia and accompanies epigenetic reprogramming and meiosis inhibition from embryonic development to post-hatch.However,mitotic arrest was not well studied in chickens because of the difficulty in obtaining pure germ cells from relevant developmental stage.Results:We performed single-cell RNA sequencing to investigate transcriptional dynamics of male germ cells during mitotic arrest in DAZL::GFP chickens.Using differentially expressed gene analysis and K-means clustering to analyze cells at different developmental stages(E12,E16,and hatch),we found that metabolic and signaling pathways were regulated,and that the epigenome was reprogrammed during mitotic arrest.In particular,we found that histone H3K9 and H3K14 acetylation(by HDAC2)and DNA demethylation(by DNMT3B and HELLS)led to a transcriptionally permissive chromatin state.Furthermore,we found that global DNA demethylation occurred gradually after the onset of mitotic arrest,indicating that the epigenetic-reprogramming schedule of the chicken genome differs from that of the mammalian genome.DNA hypomethylation persisted after hatching,and methylation was slowly re-established 3 weeks later.Conclusions:We found a unique epigenetic-reprogramming schedule of mitotic-arrested chicken prospermatogonia and prolonged hypomethylation after hatching.This will provide a foundation for understanding the process of germ-cell epigenetic regulation in several species for which this process is not clearly described.Our findings on the biological processes related to sex-specific differentiation of prospermatogonia could help studying germline development in vitro more elaborately.

Weighted single-step GWAS and RNA sequencing reveals key candidate genes associated with physiological indicators of heat stress in Holstein cattle

Background: The study of molecular processes regulating heat stress response in dairy cattle is paramount for developing mitigation strategies to improve heat tolerance and animal welfare. Therefore, we aimed to identify quantitative trait loci(QTL) regions associated with three physiological indicators of heat stress response in Holstein cattle, including rectal temperature(RT), respiration rate score(RS), and drooling score(DS). We estimated genetic parameters for all three traits. Subsequently, a weighted single-step genome-wide association study(WssGWAS) was performed based on 3200 genotypes, 151,486 phenotypic records, and 38,101 animals in the pedigree file. The candidate genes located within the identified QTL regions were further investigated through RNA sequencing(RNA-seq) analyses of blood samples for four cows collected in April(non-heat stress group) and four cows collected in July(heat stress group).Results: The heritability estimates for RT, RS, and DS were 0.06, 0.04, and 0.03, respectively. Fourteen, 19, and 20 genomic regions explained 2.94%, 3.74%, and 4.01% of the total additive genetic variance of RT, RS, and DS, respectively. Most of these genomic regions are located in the Bos taurus autosome(BTA) BTA3, BTA6, BTA8, BTA12, BTA14, BTA21, and BTA24. No genomic regions overlapped between the three indicators of heat stress, indicating the polygenic nature of heat tolerance and the complementary mechanisms involved in heat stress response. For the RNA-seq analyses, 2627 genes were significantly upregulated and 369 downregulated in the heat stress group in comparison to the control group. When integrating the WssGWAS, RNA-seq results, and existing literature, the key candidate genes associated with physiological indicators of heat stress in Holstein cattle are: PMAIP1, SBK1, TMEM33, GATB, CHORDC1, RTN4IP1, and BTBD7.Conclusions: Physiological indicators of heat stress are heritable and can be improved through direct selection. Fifty-three QTL regions associated with heat stress indicators confirm the polygenic nature and complex genetic determinism of heat tolerance in dairy cattle. The identified candidate genes will contribute for optimizing genomic evaluation models by assigning higher weights to genetic markers located in these regions as well as to the design of SNP panels containing polymorphisms located within these candidate genes.

Single-cell RNA sequencing reveals the landscape of maize root tips and assists in identification of cell type-specific nitrate-response genes

The root system is fundamental for maize growth and yield.Characterizing its heterogeneity and cell type-specific response to nitrate at the single-cell level will shed light on root development and nutrient uptake.We profiled the transcriptomes of>7000 cells derived from root tips of maize seedlings grown on media with or without nitrate,and identified 11 major cell types or tissues and 85 cell type-specific nitrate-response genes,including several known nitrate metabolic genes.A pseudotime analysis showed a continuous pseudotime series with the beginning at meristematic zone cells and showed that the root hair cell was derived by differentiation of a subset of epidermal cells.Interspecies comparison of root cells between maize and rice revealed the conservation and divergence of the root cell types and identified 57,216,and 80 conserved orthologous genes in root hair,endodermis,and phloem cells respectively.This study provides a global view of maize root tip developmental processes and responses to nitrate at the single-cell level.The genes described in the present work could serve as targets for further genetic analyses and accurate regulation of gene expression and phenotypic variation in specific cell types or tissues.

Small RNA sequencing revealed aberrant piRNA expression profiles in deciduas of recurrent spontaneous abortion patients

Piwi-interacting RNAs(piRNAs)is a novel class of non-coding RNAs.However,changes in piRNA expression profiles in recurrent spontaneous abortion(RSA)have not yet been investigated.The aim of this study was to identify differentially expressed piRNAs in deciduas of RSA patients.Decidua tissues were collected by curettage from recruited RSA patients and normal early pregnant(NEP)women with their informed consent.Small RNA sequencing was used to evaluate the differences in piRNA expression profiles between RSA and NEP.The present results demonstrated that the counts of total piRNA reads in RSA samples were increased compared with those in NEP samples(0.21%vs.0.11%).Differential expression analysis identified 29 upregulated piRNAs and 18 downregulated piRNAs in RSA samples.RT-qPCR further confirmed that the expression levels of uniq-109625,uniq-89328,uniq-50651 and uniq-4569 were decreased in 8 RSA tissues,compared with 13 NEP tissues.Otherwise,pi-22628 and uniq-173406 were increased in 8 RSA tissues.Based on GO term and KEGG pathway analysis,we speculate that these piRNAs regulate RSA by targeting extracellular matrix component pathway,cell adhesion pathway and focal adhesion pathway.PiRNAs may be involved in RSA pathogenesis by target genes function on adhesion and extracellular matrix component.