Comparative genomics: grasping human transcriptome evolution: what does it all mean?
Explore the source record for details and available documents.
SEARCH · Search PubMed
Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Explore the source record for details and available documents.
Recent molecular genetic studies have revealed that two major types of genomic instabilities, chromosomal instability (CIN) and microsatellite instability (MSI), exist in colorectal carcinomas. In order to clarify the molecular signature related to the CIN and MSI in colorectal carcinomas, we performed transcriptomic expression analysis on eight microsatellite instability-high (MSI-H) colorectal carcinomas and compared the results obtained with that of nine microsatellite stable (MSS) colorectal carcinomas using oligonucleotide microarrays containing 17 334 known genes and 1331 unknown genes or expression sequence tags (ESTs). Unsupervised two-way hierarchical clustering with 5724 genes successfully classified tumors from normal mucosa, and displayed a distinctive MSI-H carcinomas subgroup. Based on intensive filtering, 57 known genes and eight ESTs were found to be highly relevant to the differentiation of MSI-H and MSS colorectal carcinomas. These genes successfully distinguish the new test set of six MSI-H and five MSS colorectal carcinomas. Many up- and downregulated genes in MSI-H colorectal carcinomas were related to the previously reported phenotypic characteristics; increased mucin production and intense peritumoral immune response in MSI-H carcinomas. Some of these differences were confirmed by semiquantitative reverse transcription-PCR and immunohistochemical analysis. Our findings indicate that there are many different genetic and transcriptomic characteristics between MSI-H and MSS colorectal carcinomas, and some of these differently expressed genes can be used as diagnostic or prognostic markers.
BACKGROUND: Clear cell renal cell carcinoma (ccRCC) exhibits significant heterogeneity due to morphological changes and tumor microenvironment dynamics, influencing systemic therapy responses. While the role of high-mobility group AT-hook 2 (HMGA2) in tumor progression has been implicated in other cancers, its significance in ccRCC remains unclear. This study investigates the role of HMGA2 in these processes and its clinical impact. METHODS: Spatial transcriptomics (ST) was performed on primary ccRCC samples to investigate expression trajectories associated with HMGA2 expression and morphological evolution. In metastatic ccRCC cohorts treated with systemic therapy, immunohistochemistry and bulk RNA sequencing data were analyzed to evaluate molecular and clinical features in relation to HMGA2. Single-cell RNA sequencing (scRNA-seq) data were used to explore immune cell populations and their interactions. Based on these findings, multiplex immunohistochemistry (mIHC) assessed spatial distribution, cell-cell interactions, and pathological responses of key immune populations. RESULTS: HMGA2 expression was associated with aggressive morphological patterns, such as solid sheets and rhabdoid/sarcomatoid. ST revealed a progressive increase in HMGA2 expression along the morphological trajectory, marked by a shift from clear to eosinophilic cytoplasm, with eccentric nuclei and prominent nucleoli, and loss of vascular architecture. HMGA2-high tumors exhibited aggressive phenotypes driven by cell cycle, epithelial-mesenchymal transition, and inflammatory signaling pathways. Clinically, patients with high HMGA2 had worse progression-free survival but responded better to immune checkpoint inhibitor combination (Combo-ICI) therapy than to tyrosine kinase inhibitor monotherapy. To assess the immune landscape, scRNA-seq data revealed that HMGA2-high tumors were enriched with progenitor exhausted CD8+ T cells (Tpex), along with increased frequencies of conventional dendritic cell type 1 (cDC1) and inflammatory cDC type 2, which were found to interact with Tpex via ICAM-1. mIHC confirmed that Tpex were enriched among Combo-ICI responders in HMGA2-high tumors, with higher densities and closer proximity to ICAM-1+ cDC1. CONCLUSIONS: These findings suggest that dynamic HMGA2 expression contributes to morphological evolution and modulates immune responses through enhanced Tpex-cDCs engagement, serving as a potential marker for systemic therapy response in ccRCC. However, additional experimental studies are required to validate these mechanisms.
SRAP and TRAP are two kinds of newly developed molecular marker systems with the advantages of simplicity, high throughput, numerous co-dominant makers, highly reproducibility and ready to sequence, especially preferential targeting ORFs. The principle and protocol of SRAP and TRAP, which were developed in Brassica and Helianthus crops firstly, are introduced in the paper. Primer design is a key step for SRAP and TRAP, and some F-and R-primers were developed. Unlike SRAP technique, the TRAP requires cDNA or EST sequence information for fixed primer development. The annealing temperature is 35 degrees in the first 5 cycles and 50 degrees in the subsequent 30~35 cycles. The amplicons can be separated by polyacrylamide or agarose gel, and detected by autoradiography, silver or EB staining. SRAP and TRAP are applied in germplasm genetic diversity analysis, genetic map including transcriptome map construction, important trait gene tagging and gene cloning in many crops.
The human cerebral cortex is specialized into regions, but little is known about how human cellular lineages shape cortical regional variation and neuronal cell-type distribution during development. Here, we map single-cell lineages of human cortical regions and neuronal subtypes using >1,000 somatic single-nucleotide variants (sSNVs) identified from deep bulk whole-genome sequencing and analyzed over 25 regions and >72,000 single cells. In the fronto-parietal cortex, sSNVs are rarely restricted, marking neuron-generating clones that disperse into neighboring regions. In contrast, the primary visual cortex harbors 30%-70% more sSNVs than the neighboring secondary visual cortex. Clones at this border exhibit more restricted dispersion, suggesting late developmental lineage segregation. Single-nucleus sSNV and whole-transcriptome analysis reveal glutamatergic neuron clones with modest regional restrictions that share low-mosaic sSNVs with some GABAergic neurons, suggesting a recent dorsal cortical progenitor. Our analysis reveals human-specific cortical lineage patterns, regional differences in clonal patterns, and late divergence of some glutamatergic/GABAergic lineages.
INTRODUCTION: Alterations in the gut microbiota accompanied by intestinal inflammation are early features of Parkinson's disease (PD). Mutations in the leucine-rich repeat kinase 2 (LRRK2) gene represent a common genetic risk factor for PD and inflammatory bowel disease. Parabacteroides goldsteinii has been reported to alleviate intestinal and systemic inflammation. However, whether modulation of the gut microenvironment at early disease stage can attenuate PD progression remains unclear. OBJECTIVE: To investigate the impact of P. goldsteinii colonization prior to the onset of motor dysfunction on PD progression. METHODS: We established a germ-free PD mouse model carrying the LRRK2 G2019S mutation and administered P. goldsteinii orally at the pre-symptomatic stage to evaluate its effects on motor performance and PD-related neuropathology. Spatial and bulk RNA transcriptomic analyses of brain tissue, together with cytokine profiling, were conducted to assess central changes. To investigate gut immunomodulatory mechanisms, we performed intestinal bulk and single-cell RNA sequencing, spectral flow cytometry as well as cellular bioenergetic analyses. RESULTS: Germ-free conditions partially alleviated PD-like phenotypes in LRRK2 G2019S mice. Colonization with P. goldsteinii at 5-months of age, prior to motor symptom onset, further improved locomotor performance, reduced neuronal α-synuclein aggregations, and mitigated microglial activation and dopaminergic neurodegeneration. Neuroprotection was mediated through enhanced noncanonical neuronal IL-12 receptor-dependent neurotrophic support without activating the canonical STAT4 phosphorylation pathway, along with suppression of microglial activation and downregulation of LRRK2 kinase activity. At the intestinal level, P. goldsteinii suppressed TLR4-driven inflammation, expanded anti-inflammatory intraepithelial CD4+CD8αα+ T cells, promoted dendritic cell and macrophage differentiation, upregulated epithelial tight-junction genes, and improved mitochondrial bioenergetics in intestinal cells. CONCLUSION: P. goldsteinii colonization attenuates the progression of LRRK2-associated parkinsonism by restoring intestinal homeostasis and reducing neuroinflammation. These findings underscore the therapeutic potential of modulating the gut-immune-brain axis during the prodromal stage of PD.
Pre-harvest sprouting (PHS), where seeds germinate on panicles before harvest under humid conditions, is a serious global issue in cereal crop production, including rice. Fine-mapping of the previously reported chromosome 4 locus identified OsERF94 as a strong candidate gene for functional validation. In this study, we investigated the role of OsERF94 in PHS using CRISPR/Cas9 gene editing. The CRISPR/Cas9-mediated mutagenesis of OsERF94 induced frameshift mutations, resulting in a loss-of-function of OsERF94 in the 1-I-ET and 2-D-ET lines. The 1-I-ET and 2-D-ET lines exhibited significantly higher germination rates under PHS conditions compared to the wild type, indicating increased susceptibility to PHS. Whole-genome re-sequencing confirmed that few or no mutations could be detected at off-target candidate sites in both edited lines, ensuring the precision of the CRISPR/Cas9 gene editing. A transcriptome analysis revealed altered expression patterns of several GA-related genes, including OsLOL1, OsKO3, OsGA3ox2, and OsGA2ox5 in the OsERF94 mutant lines. The up-regulation of GA biosynthetic genes and the down-regulation of GA deactivation genes observed in both the OsERF94 mutant lines suggest possible alterations in GA metabolism during the early stages of PHS. Transient luciferase reporter assays using a single-luciferase system suggested that OsERF94 may be associated with changes in the promoter activities of several GA- and ethylene-related genes. These findings suggest that OsERF94 may contribute to the regulation of PHS, potentially through moderation of GA- and ethylene-related pathways. Overall, this study improves our understanding of the molecular role of OsERF94 in PHS and highlights its potential as a target for the genetic improvement of PHS resistance in rice-breeding programs.
Adrenocortical carcinoma (ACC) is a rare but aggressive malignancy with poor survival and limited treatment options. To comprehensively characterize its molecular landscape and identify clinically relevant subtypes, we performed an integrated genomic analysis - including whole-exome sequencing, RNA sequencing, and copy number variation profiling - on 61 Chinese patients with ACC. We identified recurrent mutations in TP53 (25%), CTNNB1 (15%), ZNRF3 (10%), and MEN1 (8%). Unsupervised clustering of transcriptomic data revealed four distinct molecular subtypes: cortisol-driven (CD, 14%), immune-suppressed (IS, 40%), cell cycle-altered (CCA, 22%), and immunomodulatory (IM, 24%). The CD subtype exhibited steroidogenic pathway activation; the IS subtype showed T cell receptor downregulation and the worst disease-free survival; the CCA subtype was marked by chromosomal instability and cell cycle gene overexpression; and the IM subtype displayed enriched immune signaling and favorable outcomes. Copy number analysis further uncovered focal amplifications (e.g. TERT, CDK4) and HLA-II deletions. This study establishes a novel molecular classification of ACC, providing a framework for subtype-specific therapeutic strategies, such as CDK4/6 inhibition for CCA and immunotherapy for IM tumors, while highlighting the clinical challenges of immune-cold IS tumors.
The envelope is an essential structure of the bacterial cell, and maintaining its integrity is a prerequisite for survival. To ensure proper function, transmembrane signal-transducing systems, such as two-component systems (TCS) and extracytoplasmic function (ECF) sigma factors, closely monitor its condition and respond to harmful perturbations. Both systems consist of a transmembrane sensor protein (histidine kinase or anti-sigma factor, respectively) and a corresponding cytoplasmic transcriptional regulator (response regulator or sigma factor, respectively) that mediates the cellular response through differential gene expression. The regulatory network of the cell envelope stress response is well studied in the gram-positive model organism Bacillus subtilis. It consists of at least two ECF sigma factors and four two-component systems. In this study, we describe the corresponding network in a close relative, Bacillus licheniformis. Based on sequence homology, domain architecture, and genomic context, we identified five TCS and eight ECF sigma factors as potential candidate regulatory systems mediating cell envelope stress response in this organism. We characterized the corresponding regulatory network by comparative transcriptomics and regulon mining as an initial screening tool. Subsequent in-depth transcriptional profiling was applied to define the inducer specificity of each identified cell envelope stress sensor. A total of three TCS and seven ECF sigma factors were shown to be induced by cell envelope stress in B. licheniformis. We noted a number of significant differences, indicative of a regulatory divergence between the two Bacillus species, in addition to the expected overlap in the respective responses.
Using small interfering RNA (siRNA) to induce sequence specific gene silencing is fast becoming a standard tool in functional genomics. As siRNAs in some cases tolerate mismatches with the mRNA target, knockdown of genes other than the intended target could make results difficult to interpret. In an investigation of 359 published siRNA sequences, we have found that about 75% of them have a risk of eliciting non-specific effects. A possible cause for this is the popular BLAST search engine, which is inappropriate for such short oligos as siRNAs. Furthermore, we used new special purpose hardware to do a transcriptome-wide screening of all possible siRNAs, and show that many unique siRNAs exist per target even if several mismatches are allowed. Hence, we argue that the risk of off-target effects is unnecessary and should be avoided in future siRNA design.
The Hrp type III protein secretion system is essential for pathogenicity of the Gram-negative plant pathogen Xanthomonas campestris pv. vesicatoria. Expression of the hrp gene cluster is controlled by HrpG, a two-component response regulator, and HrpX, an AraC-type transcriptional activator. Using the cDNA-AFLP technique, 30 hrpG-induced (hgi) and five hrpG-repressed (hgr) cDNA fragments were identified, defining a large hrpG-regulon in X. campestris pv. vesicatoria. Expression of most genes in the hrpG-regulon was dependent on hrpX. Seven cDNA fragments map to the known hrp gene cluster and flanking regions. All other genes appear to be scattered over the chromosome and endogenous plasmids. Sequence analysis identified genes encoding putative extracellular proteases, a putative transcriptional regulator and XopJ and XopB (Xanthomonas outer proteins), homologues of YopJ from Yersinia spp. and the avirulence protein AvrPphD of Pseudomonas syringae respectively. XopB is secreted by the Hrp type III secretion system. Analysis of deletion mutants in several hgi genes revealed a new virulence locus. This study demonstrates that cDNA-AFLP is a powerful tool to study prokaryotic transcriptomes and to identify genes contributing to Xanthomonas virulence and putative effector proteins.
Chiari 1 malformation (CM1) is the most common congenital malformation of the human hindbrain. Although prior studies have implicated chromatin-remodeling genes in CM1, the de novo genetic architecture and underlying neurodevelopmental mechanisms remain incompletely defined. To investigate the molecular genetics of a novel familial form of CM1 linked with syringomyelia and tethered cord and determine whether rare, damaging de novo variants (DNVs) contribute to sporadic CM1 risk with gene- and pathway-level resolution, we performed whole-exome sequencing in an ultra-rare multigenerational family with CM1 and associated spinal pathology, and in the largest assembled trio-based cohort to date, comprising 1,585 proband-parent trios with sporadic, idiopathic CM1 (2017-2025). The comparison cohort included 1,798 unaffected control siblings. Clinical phenotyping was by systematic medical record review. Structural domain mapping, in silico modeling, and integration with single-cell transcriptomic data from developing human cerebellum was conducted to assess biological plausibility. A heterozygous loss-of-function variant in CHD3 segregated with CM1 and syringomyelia in a multigenerational family. In the trio-based cohort, rare protein-altering DNVs were significantly enriched across multiple chromodomain helicase DNA-binding (CHD) genes, including CHD1, CHD3, CHD4, and CHD8, exceeding gene-specific mutation expectations (protein-damaging variants: P = 1.3 × 10-9; predicted loss-of-function variants: P = 8.6 × 10-5). CHD1 contained two pathogenic DNVs (p.A999D and p.E984K). CHD4 (p.D744N, p.T1813P, and p.I1102T) and CHD8 (p.R1402X, p.R1472X, and p.R2035X) each contained three new DNVs. Variants clustered within conserved ATPase, helicase, and chromodomain regions essential for chromatin remodeling, and these patients frequently had comorbid developmental delay and related neurodevelopmental features. Single-cell transcriptomic analyses demonstrated enrichment in Purkinje cells and inhibitory neurons of midgestational cerebellum, where CHD gene products form a coherent chromatin-regulatory network. Rare, large-effect DNVs that disrupt chromatin-remodeling programs contribute to sporadic CM1, implicating genetically encoded dysregulation of cerebellar development as a central disease mechanism. Exome sequencing may complement surgical evaluation of children with sporadic CM1, particularly when accompanied by neurodevelopmental concerns, informing prognosis and family counseling.
BACKGROUND: Seizures are prevalent in glioma patients, especially in those with low-grade gliomas. The interaction between gliomas and epilepsy involves complex biological mechanisms that are not fully understood. METHODS: We collected Genome-Wide Association Study data for epilepsy and gliomas, performed differential expression analysis, and conducted Gene Ontology (GO) enrichment analysis on the identified genes. Single-cell RNA sequencing data (scRNA-seq) from GSE221534 dataset in Gene Expression Omnibus (GEO) were used to analyze cell-cell interactions within glioma samples from patients with and without epilepsy. RESULTS: Mendelian Randomization (MR) analysis revealed significant associations between genetic variants related to epilepsy and glioma risk, suggesting a potential causal relationship, especially in astrocytomas. Differential expression analysis identified epilepsy-related genes that were significantly upregulated in astrocytoma tissues compared to normal brain tissues. GO enrichment analysis indicated that these genes are involved in critical biological processes such as neurogenesis and cellular signaling. The scRNA-seq analysis showed, compared to non-epileptic samples, glioma stem cells, microglia, and NK cells are increased in the core regions of astrocytomas in epileptic patients. Additionally, intercellular communication between tumor cells and other non-tumor cells is markedly enhanced in astrocytoma samples from epileptic patients. CONCLUSION: This study provides evidence of a genetic association between epilepsy and gliomas and elucidates the biological mechanisms through which epilepsy may influence glioma progression.
Microarray technology for mammalian cells has been utilized mainly for humans, mouse, and rat gene expression analysis. In this approach the feasibility of cross-species hybridization experiments using Chinese hamster ovary (CHO) cells was evaluated. Sequence alignments of available data for CHO were performed against mouse and rat transcripts to determine the homology between the investigated species. We implemented a probability model based on this homology in order to estimate the chance for successful hybridization using Agilent's 60-mer oligonucleotide platform. Heat-shock expression data from CHO, mouse 3T3, and rat A10 cells were generated to determine intraspecies variability, reproducibility, and specificity in order to assess the accuracy of this method. Detected signature genes, in particular from studies with the mouse arrays, showed a reliable similarity between these two rodents and were confirmed by quantitative RT-PCR. Our findings provide evidence that cross-species analysis can be a useful tool to study gene expression profiles of related organisms for which species-specific microarrays are not available.
The female Aedes aegypti mosquito's remarkable ability to hunt humans and transmit pathogens relies on her unique biology. Here, we present the Aedes aegypti Mosquito Cell Atlas, a comprehensive single-nucleus RNA sequencing dataset of more than 367,000 nuclei from 19 dissected tissues of adult female and male Aedes aegypti, providing cellular-level resolution of mosquito biology. We identify novel cell types and expand our understanding of sensory neuron organization of chemoreceptors to all sensory tissues. Our analysis uncovers male-specific cells and sexually dimorphic gene expression in the antenna and brain. In female mosquitoes, we find that glial cells in the brain, rather than neurons, undergo the most extensive transcriptional changes following blood feeding. Our findings provide insights into the cellular basis of mosquito behavior and sexual dimorphism. The Aedes aegypti Mosquito Cell Atlas resource enables systematic investigation of cell type-specific expression across all mosquito tissues.
Upland cotton (Gossypium hirsutum) produces the most widely used natural fibers, yet the regulatory mechanisms governing fiber cell elongation are not well understood. Through sequencing of a cotton fiber cDNA library and subsequent microarray analysis, we found that ethylene biosynthesis is one of the most significantly upregulated biochemical pathways during fiber elongation. The 1-Aminocyclopropane-1-Carboxylic Acid Oxidase1-3 (ACO1-3) genes responsible for ethylene production were expressed at significantly higher levels during this growth stage. The amount of ethylene released from cultured ovules correlated with ACO expression and the rate of fiber growth. Exogenously applied ethylene promoted robust fiber cell expansion, whereas its biosynthetic inhibitor l-(2-aminoethoxyvinyl)-glycine (AVG) specifically suppressed fiber growth. The brassinosteroid (BR) biosynthetic pathway was modestly upregulated during this growth stage, and treatment with BR or its biosynthetic inhibitor brassinazole (BRZ) also promoted or inhibited, respectively, fiber growth. However, the effect of ethylene treatment was much stronger than that of BR, and the inhibitory effect of BRZ on fiber cells could be overcome by ethylene, but the AVG effect was much less reversed by BR. These results indicate that ethylene plays a major role in promoting cotton fiber elongation. Furthermore, ethylene may promote cell elongation by increasing the expression of sucrose synthase, tubulin, and expansin genes.
Two collections of oligonucleotides have been designed for preparing pangenomic human and mouse microarrays. A total of 148,993 and 121,703 oligonucleotides were designed against human and mouse transcripts. Quality scores were created in order to select 25,342 human and 24,109 mouse oligonucleotides. They correspond to: (i) a BLAST-specificity score; (ii) the number of expressed sequence tags matching each probe; (iii) the distance to the 3' end of the target mRNA. Scores were also used to compare in silico the two microarrays with commercial microarrays. The sets described here, called RNG/MRC collections, appear at least as specific and sensitive as those from the commercial platforms. The RNG/MRC collections have now been used by an Anglo-French consortium to distribute more than 3500 microarrays to the academic community. Ad hoc identification of tissue-specific transcripts and a approximately 80% correlation with hybridizations performed on Affymetrix GeneChiptrade mark suggest that the RNG/MRC microarrays perform well. This work provides a comprehensive open resource for investigators working on human and mouse transcriptomes, as well as a generic method to generate new microarray collections in other organisms. All information related to these probes, as well as additional information about commercial microarrays have been stored in a freely-accessible database called MEDIANTE.
Little is known about the transcriptomic plasticity and adaptive mechanisms of circulating tumor cells (CTCs) during hematogeneous dissemination. Here we interrogate the transcriptome of 113 single CTCs from 4 different vascular sites, including hepatic vein (HV), peripheral artery (PA), peripheral vein (PV) and portal vein (PoV) using single-cell full-length RNA sequencing in hepatocellular carcinoma (HCC) patients. We reveal that the transcriptional dynamics of CTCs were associated with stress response, cell cycle and immune-evasion signaling during hematogeneous transportation. Besides, we identify chemokine CCL5 as an important mediator for CTC immune evasion. Mechanistically, overexpression of CCL5 in CTCs is transcriptionally regulated by p38-MAX signaling, which recruites regulatory T cells (Tregs) to facilitate immune escape and metastatic seeding of CTCs. Collectively, our results reveal a previously unappreciated spatial heterogeneity and an immune-escape mechanism of CTC, which may aid in designing new anti-metastasis therapeutic strategies in HCC.