Identification of cell-type-specific marker genes from co-expression patterns in tissue samples.
Ontology highlight
ABSTRACT: Motivation
Marker genes, defined as genes that are expressed primarily in a single cell type, can be identified from the single cell transcriptome; however, such data are not always available for the many uses of marker genes, such as deconvolution of bulk tissue. Marker genes for a cell type, however, are highly correlated in bulk data, because their expression levels depend primarily on the proportion of that cell type in the samples. Therefore, when many tissue samples are analyzed, it is possible to identify these marker genes from the correlation pattern.Results
To capitalize on this pattern, we develop a new algorithm to detect marker genes by combining published information about likely marker genes with bulk transcriptome data in the form of a semi-supervised algorithm. The algorithm then exploits the correlation structure of the bulk data to refine the published marker genes by adding or removing genes from the list.Availability and implementation
We implement this method as an R package markerpen, hosted on CRAN (https://CRAN.R-project.org/package=markerpen).Supplementary information
Supplementary data are available at Bioinformatics online.
SUBMITTER: Qiu Y
PROVIDER: S-EPMC8504631 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:BackgroundCellular senescence is a complex stress response that impacts cellular function and organismal health. Multiple developmental and environmental factors, such as intrinsic cellular cues, radiation, oxidative stress, oncogenes, and protein accumulation, activate genes and pathways that can lead to senescence. Enormous efforts have been made to identify and characterize senescence genes (SnGs) in stress and disease systems. However, the prevalence of senescent cells in healthy human tissues and the global SnG expression signature in different cell types are poorly understood.MethodsThis study performed an integrative gene network analysis of bulk and single-cell RNA-seq data in non-diseased human tissues to investigate SnG co-expression signatures and their cell-type specificity.ResultsThrough a comprehensive transcriptomic network analysis of 50 human tissues in the Genotype-Tissue Expression Project (GTEx) cohort, we identified SnG-enriched gene modules, characterized SnG co-expression patterns, and constructed aggregated SnG networks across primary tissues of the human body. Our network approaches identified 51 SnGs highly conserved across the human tissues, including CDKN1A (p21)-centered regulators that control cell cycle progression and the senescence-associated secretory phenotype (SASP). The SnG-enriched modules showed remarkable cell-type specificity, especially in fibroblasts, endothelial cells, and immune cells. Further analyses of single-cell RNA-seq and spatial transcriptomic data independently validated the cell-type specific SnG signatures predicted by the network analysis.ConclusionsThis study systematically revealed the co-regulated organizations and cell type specificity of SnGs in major human tissues, which can serve as a blueprint for future studies to map senescent cells and their cellular interactions in human tissues.
| S-EPMC8744330 | biostudies-literature
Project description:Cell-type specific genes were recognized by interrogating microarrays carrying Dictyostelium gene fragments with probes prepared from fractions enriched in prestalk and prespore cells. Cell-type specific accumulation of mRNA from 17 newly identified genes was confirmed by Northern analyses. DNA microarrays carrying 690 targets were used to determine expression profiles during development. The profiles were fit to a biologically based kinetic equation to extract the times of transcription onset and cessation. Although the majority of the genes that were cell-type enriched at the slug stage were first expressed as the prespore and prestalk cells sorted out in aggregates, some were found to be expressed earlier before the cells had even aggregated. These early genes may have been initially expressed in all cells and then preferentially turned over in one or the other cell type. Alternatively, cell type divergence may start soon after the initiation of development.
| S-EPMC59697 | biostudies-literature
Project description:BACKGROUND:Cellular development requires the precise control of gene expression states. Transcription factors are involved in this regulatory process through their combinatorial binding with DNA. Information about transcription factor binding sites can help determine which combinations of factors work together to regulate a gene, but it is unclear how far the binding data from one cell type can inform about regulation in other cell types. RESULTS:By integrating data on co-localized transcription factor binding sites in the K562 cell line with expression data across 38 distinct hematopoietic cell types, we developed regression models to describe the relationship between the expression of target genes and the transcription factors that co-localize nearby. With K562 binding sites identifying the predictors, the proportion of expression explained by the models is statistically significant only for monocytic cells (p-value< 0.001), which are closely related to K562. That is, cell type specific binding patterns are crucial for choosing the correct transcription factors for the model. Comparison of predictors obtained from binding sites in the GM12878 cell line with those from K562 shows that the amount of difference between binding patterns is directly related to the quality of the prediction. By identifying individual genes whose expression is predicted accurately by the binding sites, we are able to link transcription factors FOS, TAF1 and YY1 to a sparsely studied gene LRIG2. We also find that the activity of a transcription factor may be different depending on the cell type and the identity of other co-localized factors. CONCLUSION:Our approach shows that gene expression can be explained by a modest number of co-localized transcription factors, however, information on cell-type specific binding is crucial for understanding combinatorial gene regulation.
| S-EPMC3441573 | biostudies-literature
Project description:Transcriptome diversity provides the key to cellular identity. One important contribution to expression diversity is the use of alternative promoters, which creates mRNA isoforms by expanding the choice of transcription initiation sites of a gene. The proximity of the basal promoter to the transcription initiation site enables prediction of a promoter's location based on the gene annotations. We show that annotation of alternative promoters regulating expression of transcripts with distinct first exons enables a novel methodology to quantify expression levels and tissue specificity of mRNA isoforms.The use of distinct alternative first exons in 3,296 genes was examined using exon-microarray data from 11 human tissues. Comparing two transcripts from each gene we found that the activity of alternative promoters (i.e., P1 and P2) was not correlated through tissue specificity or level of expression. Furthermore neither P1 nor P2 conferred any bias for tissue-specific or ubiquitous expression. Genes associated with specific diseases produced transcripts whose limited expression patterns were consistent with the tissue affected in disease. Notably, genes that were historically designated as tissue-specific or housekeeping had alternative isoforms that showed differential expression. Furthermore, only a small number of alternative promoters showed expression exclusive to a single tissue indicating that "tissue preference" provides a better description of promoter activity than tissue specificity. When compared to gene expression data in public databases, as few as 22% of the genes had detailed information for more than one isoform, whereas the remainder collapsed the expression patterns from individual transcripts into one profile.We describe a computational pipeline that uses microarray data to assess the level of expression and breadth of tissue profiles for transcripts with distinct first exons regulated by alternative promoters. We conclude that alternative promoters provide individualized regulation that is confirmed through expression levels, tissue preference and chromatin modifications. Although the selective use of alternative promoters often goes uncharacterized in gene expression analyses, transcripts produced in this manner make unique contributions to the cell that requires further exploration.
| S-EPMC2923625 | biostudies-literature
Project description:BackgroundCaenorhabditis elegans has traditionally been used as a model for studying nematode biology, but its small size limits the ability for researchers to perform some experiments such as high-throughput tissue-specific gene expression studies. However, the dissection of individual tissues is possible in the parasitic nematode Ascaris suum due to its relatively large size. Here, we take advantage of the recent genome sequencing of Ascaris suum and the ability to physically dissect its separate tissues to produce a wide-scale tissue-specific nematode RNA-seq datasets, including data on three non-reproductive tissues (head, pharynx, and intestine) in both male and female worms, as well as four reproductive tissues (testis, seminal vesicle, ovary, and uterus). We obtained fundamental information about the biology of diverse cell types and potential interactions among tissues within this multicellular organism.Methodology/principal findingsOverexpression and functional enrichment analyses identified many putative biological functions enriched in each tissue studied, including functions which have not been previously studied in detail in nematodes. Putative tissue-specific transcriptional factors and corresponding binding motifs that regulate expression in each tissue were identified, including the intestine-enriched ELT-2 motif/transcription factor previously described in nematode intestines. Constitutively expressed and novel genes were also characterized, with the largest number of novel genes found to be overexpressed in the testis. Finally, a putative acetylcholine-mediated transcriptional network connecting biological activity in the head to the male reproductive system is described using co-expression networks, along with a similar ecdysone-mediated system in the female.Conclusions/significanceThe expression profiles, co-expression networks and co-expression regulation of the 10 tissues studied and the tissue-specific analysis presented here are a valuable resource for studying tissue-specific biological functions in nematodes.
| S-EPMC3916258 | biostudies-literature
Project description:The patterns of evolution and expression of tissue-specific genes are poorly understood beyond sex-specific genes. Accordingly, we identified 3,191 tissue-specific genes and 38,745 common genes using 22 RNA-seq datasets from cultivated peanut. The expression levels of tissue-specific genes were significantly lower than those of common genes. Further, the expression levels of sex-specific genes were significantly higher than those of somatic tissue-specific genes. Among sex-specific genes, the expression levels of gynoecium-specific genes were significantly higher than those of androecium-specific genes. Function-specific genes were lacking among tissue-specific genes, and tissue-specific gene annotations overlapped among different tissues. Duplicate gene pairs were classified as homogeneous pairs expressed within the same tissue or heterogeneous pairs expressed in different tissues. Heterogeneous gene pairs evolved more rapidly than homogeneous gene pairs. In addition, somatic tissue-specific genes evolved faster than sex-specific genes. Molecular signatures of selection indicated that somatic tissue-specific genes have mainly experienced relaxed selection, while sex-specific genes have been under stronger selective constraint. Somatic tissue-specific genes had higher codon usage bias than sex-specific genes. These contrasting patterns between somatic tissue-specific and sex-specific genes provide new insights into the basic biology and evolution of peanut, an important crop.
| S-EPMC5566475 | biostudies-literature
Project description:AimTo determine the gene expression profile data for the whole liver during development of dimethylni-trosamine (DMN)-induced hepatic fibrosis.MethodsMarker genes were identified for different types of hepatic cells, including hepatic stellate cells (HSCs), Kupffer cells (including other inflammatory cells), and hepatocytes, using independent temporal DNA microarray data obtained from isolated hepatic cells.ResultsThe cell-type analysis of gene expression gave several key results and led to formation of three hypotheses: (1) changes in the expression of HSC-specific marker genes during fibrosis were similar to gene expression data in in vitro cultured HSCs, suggesting a major role of the self-activating characteristics of HSCs in formation of fibrosis; (2) expression of mast cell-specific marker genes reached a peak during liver fibrosis, suggesting a possible role of mast cells in formation of fibrosis; and (3) abnormal expression of hepatocyte-specific marker genes was found across several metabolic pathways during fibrosis, including sulfur-containing amino acid metabolism, fatty acid metabolism, and drug metabolism, suggesting a mechanistic relationship between these abnormalities and symptoms of liver fibrosis.ConclusionAnalysis of marker genes for specific hepatic cell types can identify the key aspects of fibrogenesis. Sequential activation of inflammatory cells and the self-supporting properties of HSCs play an important role in development of fibrosis.
| S-EPMC4100637 | biostudies-literature
Project description:BackgroundThe genetic background of cancer remains complex and challenging to integrate. Many somatic mutations within genes are known to cause and drive cancer, while genome-wide association studies (GWAS) of cancer have revealed many germline risk factors associated with cancer. However, the overlap between known somatic driver genes and positional candidate genes from GWAS loci is surprisingly small. We hypothesised that genes from multiple independent cancer GWAS loci should show tissue-specific co-regulation patterns that converge on cancer-specific driver genes.ResultsWe studied recent well-powered GWAS of breast, prostate, colorectal and skin cancer by estimating co-expression between genes and subsequently prioritising genes that show significant co-expression with genes mapping within susceptibility loci from cancer GWAS. We observed that the prioritised genes were strongly enriched for cancer drivers defined by COSMIC, IntOGen and Dietlein et al. The enrichment of known cancer driver genes was most significant when using co-expression networks derived from non-cancer samples of the relevant tissue of origin.ConclusionWe show how genes within risk loci identified by cancer GWAS can be linked to known cancer driver genes through tissue-specific co-expression networks. This provides an important explanation for why seemingly unrelated sets of genes that harbour either germline risk factors or somatic mutations can eventually cause the same type of disease.
| S-EPMC11247850 | biostudies-literature
Project description:Identification of marker genes associated with a specific tissue/cell type is a fundamental challenge in genetic and cell research. Marker genes are of great importance for determining cell identity, and for understanding tissue specific gene function and the molecular mechanisms underlying complex diseases.We have developed a new bioinformatics tool called MGFM (Marker Gene Finder in Microarray data) to predict marker genes from microarray gene expression data. Marker genes are identified through the grouping of samples of the same type with similar marker gene expression levels. We verified our approach using two microarray data sets from the NCBI's Gene Expression Omnibus public repository encompassing samples for similar sets of five human tissues (brain, heart, kidney, liver, and lung). Comparison with another tool for tissue-specific gene identification and validation with literature-derived established tissue markers established functionality, accuracy and simplicity of our tool. Furthermore, top ranked marker genes were experimentally validated by reverse transcriptase-polymerase chain reaction (RT-PCR). The sets of predicted marker genes associated with the five selected tissues comprised well-known genes of particular importance in these tissues. The tool is freely available from the Bioconductor web site, and it is also provided as an online application integrated into the CellFinder platform ( http://cellfinder.org/analysis/marker ).MGFM is a useful tool to predict tissue/cell type marker genes using microarray gene expression data. The implementation of the tool as an R-package as well as an application within CellFinder facilitates its use.
| S-EPMC4552366 | biostudies-literature
Project description:Promoters play a crucial role in controlling the spatial and temporal expression of genes at transcriptional levels in the process of higher plant growth and development. The spatial, efficient, and correct regulation of exogenous genes expression, as desired, is the key point in plant genetic engineering research. Constitutive promoters widely used in plant genetic transformation are limited because, sometimes, they may cause potential negative effects. This issue can be solved, to a certain extent, by using tissue-specific promoters. Compared with constitutive promoters, a few tissue-specific promoters have been isolated and applied. In this study, based on the transcriptome data, a total of 288 tissue-specific genes were collected, expressed in seven tissues, including the leaves, stems, flowers, pods, seeds, roots, and nodules of soybean (Glycine max). KEGG pathway enrichment analysis was carried out, and 52 metabolites were annotated. A total of 12 tissue-specific genes were selected via the transcription expression level and validated through real-time quantitative PCR, of which 10 genes showed tissue-specific expression. The 3-kb 5' upstream regions of ten genes were obtained as putative promoters. Further analysis showed that all the 10 promoters contained many tissue-specific cis-elements. These results demonstrate that high-throughput transcriptional data can be used as effective tools, providing a guide for high-throughput novel tissue-specific promoter discovery.
| S-EPMC10298445 | biostudies-literature