SEARCH · Search PubMed

Results for “transcriptome annotation”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 379 records · Page 21Linked to original sources

PubMed2005

Identification of hypertension-related genes through an integrated genomic-transcriptomic approach.

In search for the genetic basis of hypertension, we applied an integrated genomic-transcriptomic approach to identify genes involved in the pathogenesis of hypertension in the Sabra rat model of salt-susceptibility. In the genomic arm of the project, we previously detected in male rats two salt-susceptibility QTLs on chromosome 1, SS1a (D1Mgh2-D1Mit11; span 43.1 cM) and SS1b (D1Mit11-D1Mit4; span 18 cM). In the transcriptomic arm, we studied differential gene expression in kidneys of SBH/y and SBN/y rats that had been fed regular diet or salt-loaded. We used the Affymetrix Rat Genome RAE230 GeneChip and probed >30,000 transcripts. The research algorithm called for an initial genome-wide screen for differentially expressed transcripts between the study groups. This step was followed by cluster analysis based on 2x2 ANOVA to identify transcripts that were of relevance specifically to salt-sensitivity and hypertension and to salt-resistance. The two arms of the project were integrated by identifying those differentially expressed transcripts that showed an allele-specific hypertensive effect on salt-loading and that mapped within the defined boundaries of the salt-susceptibility QTLs on chromosome 1. The differentially expressed transcripts were confirmed by RT-PCR. Of the 2933 genes annotated to rat chromosome 1, 1102 genes were identified within the boundaries of the two blood pressure QTLs. The microarray identified 2470 transcripts that were differentially expressed between the study groups. Cluster analysis identified genome-wide 192 genes that were relevant to salt-susceptibility and/or hypertension, 19 of which mapped to chromosome 1. Eight of these genes mapped within the boundaries of QTLs SS1a and SS1b. RT-PCR confirmed 7 genes, leaving TcTex1, Myadm, Lisch7, Axl-like, Fah, PRC1-like, and Serpinh1. None of these genes has been implicated in hypertension before. These genes become henceforth targets for our continuing search for the genetic basis of hypertension.

Algorithms↗

PubMed2004

Identification and functional analysis of 'hypothetical' genes expressed in Haemophilus influenzae.

The progress in genome sequencing has led to a rapid accumulation in GenBank submissions of uncharacterized 'hypothetical' genes. These genes, which have not been experimentally characterized and whose functions cannot be deduced from simple sequence comparisons alone, now comprise a significant fraction of the public databases. Expression analyses of Haemophilus influenzae cells using a combination of transcriptomic and proteomic approaches resulted in confident identification of 54 'hypothetical' genes that were expressed in cells under normal growth conditions. In an attempt to understand the functions of these proteins, we used a variety of publicly available analysis tools. Close homologs in other species were detected for each of the 54 'hypothetical' genes. For 16 of them, exact functional assignments could be found in one or more public databases. Additionally, we were able to suggest general functional characterization for 27 more genes (comprising approximately 80% total). Findings from this analysis include the identification of a pyruvate-formate lyase-like operon, likely to be expressed not only in H.influenzae but also in several other bacteria. Further, we also observed three genes that are likely to participate in the transport and/or metabolism of sialic acid, an important component of the H.influenzae lipo-oligosaccharide. Accurate functional annotation of uncharacterized genes calls for an integrative approach, combining expression studies with extensive computational analysis and curation, followed by eventual experimental verification of the computational predictions.

Amino Acid Sequence↗

PubMed2026

Contrasting regulation of protein-coding genes and lncRNA homeologs in allotetraploid Coffea arabica.

A chromosome-level Bourbon assembly revealed that protein-coding homeologs are predominantly co-regulated between subgenomes. In contrast, intergenic lncRNAs display a modest, but statistically consistent bias toward subgenome E across diverse developmental and stress contexts. Coffea arabica is an allotetraploid species derived from natural hybridization between C. canephora and C. eugenioides, which contributed the C and E subgenomes, respectively. This genomic origin poses major challenges for genome assembly, annotation, and the interpretation of gene regulation. In this study, a high-quality genome assembly of C. arabica was generated and annotated, with particular emphasis on identifying protein-coding genes and intergenic long non-coding RNAs (lincRNAs). Homeologous relationships between genes from the C and E subgenomes were established, providing a robust framework to investigate subgenomic conservation and regulatory divergence. Using an extensive collection of publicly available RNA-seq libraries spanning multiple developmental stages, tissues, and environmental conditions, the relative transcriptional contribution of each subgenome was evaluated. On a global scale, gene expression was largely balanced between subgenomes, with no consistent evidence of subgenome dominance. While protein-coding genes showed comparable regulatory behavior across subgenomes, lincRNAs exhibited a more asymmetric expression pattern, suggesting higher subgenome-specific expression that is interpreted here as a consistent directional tendency rather than as evidence of subgenome dominance. Together, these results provide new insights into the regulatory architecture of the C. arabica genome and establish a foundational genomic and transcriptomic resource for future functional studies and crop improvement efforts.

Coffea↗

PubMed2005

The expanding role of microarrays in the investigation of macrophage responses to pathogens.

In the last few years, microarray technology has emerged as the method of choice for large-scale gene expression studies. It provides an efficient and rapid method to investigate the entire transcriptome of a cell. No research field has benefited more from microarray technology than the study of the exquisite interplay between pathogens and hosts. Numerous microarray studies have now been published in this field, which have provided insights into the mechanisms of host defence and the tactics employed by pathogens to circumvent these protection strategies. These studies have led to a more comprehensive understanding of the host immune response and identified new avenues of research for potential control strategies against pathogens. In the past, research has concentrated on human and mouse microarrays to investigate host-pathogen interactions, regardless of the host species. This trend is changing with the ever-expanding sequence resources now available for many pathogen and host species, including livestock animals. The use of species-specific microarrays has furthered our understanding of host-pathogen interactions for particular organisms and aided in the annotation of unknown genes. Macrophages play a central role in the host's innate and adaptive immune responses to pathogens. These cells are in the first line of defence and interact with a wide range of pathogens; many of which have evolved strategies to circumvent the macrophage defence mechanisms and survive within these cells. In this report, we review the wealth of studies using microarray technology to investigate the response of macrophages to pathogens. These studies illustrate how microarray technology has expanded our understanding of the dialogue between macrophage and pathogen and provide examples of the benefits and pitfalls of using this technique. Furthermore, we discuss the resources available to use microarray analysis to study the immune response of a non-human, non-rodent species, the cow.

Animals↗

PubMed2026

The machine-learning classifier ALLCatchR2 identifies 20 T-ALL subtypes across cohorts and age groups.

T-cell acute lymphoblastic leukemia (T-ALL) comprises molecularly diverse subtypes, but robust cross-cohort validations and operational gene-expression definitions are lacking. To establish a gene-expression-anchored framework for T-ALL subtyping, we aggregated 2314 transcriptomes (15 cohorts, age: 0.8-90.8 years). An extended unsupervised approach defined 17 main clusters and 3 subclusters in samples with high blast fractions. Supervised analyses added an overarching immature T-ALL (early T cell precursor [ETP]-like) definition and resolved the LMO2 γδ-like subtype. All clusters contained samples from at least two cohorts. Characteristic genomic driver enrichments were consistent across cohorts, while gene-expression clusters did not correspond exclusively to single driver events but also reflected developmental origins. A machine-learning classifier based on ALLCatchR, our B-cell acute lymphoblastic leukemia (B-ALL) classifier, identified these 20 transcriptomic subtypes and the immature T-ALL (ETP-like) signature with 0.995-1.0 accuracy in a validation set (n = 203). Testing the classifier on a second hold-out data set (n = 265 samples) showed that 92.7% of predictions matched with corresponding driver alterations. Across all samples, 83.2% of cases received high-confidence predictions, 7.3% candidate predictions, and 9.5% remained unclassified, largely because of low blast fractions. We identified a novel gene-expression cluster markedly enriched (P < 0.001) for clonal hematopoiesis mutations (IDH2 R140Q, DNMT3A) and a stem-/progenitor cell-like gene expression. This novel clonal hematopoiesis-related T-ALL subtype was observed in six cohorts and accounted for 8.9% of adults and 39.5% of patients aged >50 years. We extended ALLCatchR into ALLCatchR2, a free R package that now enables B-/T-lineage separation, gene-expression subtyping, blast estimation, and developmental annotation to harmonize T-ALL classification across studies and clinical contexts.

Journal Article↗

PubMed2005

cGMP-independent nitric oxide signaling and regulation of the cell cycle.

BACKGROUND: Regulatory functions of nitric oxide (NO*) that bypass the second messenger cGMP are incompletely understood. Here, cGMP-independent effects of NO* on gene expression were globally examined in U937 cells, a human monoblastoid line that constitutively lacks soluble guanylate cyclase. Differentiated U937 cells (>80% in G0/G1) were exposed to S-nitrosoglutathione, a NO* donor, or glutathione alone (control) for 6 h without or with dibutyryl-cAMP (Bt2cAMP), and then harvested to extract total RNA for microarray analysis. Bt2cAMP was used to block signaling attributable to NO*-induced decreases in cAMP. RESULTS: NO* regulated 110 transcripts that annotated disproportionately to the cell cycle and cell proliferation (47/110, 43%) and more frequently than expected contained AU-rich, post-transcriptional regulatory elements (ARE). Bt2cAMP regulated 106 genes; cell cycle gene enrichment did not reach significance. Like NO*, Bt2cAMP was associated with ARE-containing transcripts. A comparison of NO* and Bt2cAMP effects showed that NO* regulation of cell cycle genes was independent of its ability to interfere with cAMP signaling. Cell cycle genes induced by NO* annotated to G1/S (7/8) and included E2F1 and p21/Waf1/Cip1; 6 of these 7 were E2F target genes involved in G1/S transition. Repressed genes were G2/M associated (24/27); 8 of 27 were known targets of p21. E2F1 mRNA and protein were increased by NO*, as was E2F1 binding to E2F promoter elements. NO* activated p38 MAPK, stabilizing p21 mRNA (an ARE-containing transcript) and increasing p21 protein; this increased protein binding to CDE/CHR promoter sites of p21 target genes, repressing key G2/M phase genes, and increasing the proportion of cells in G2/M. CONCLUSION: NO* coordinates a highly integrated program of cell cycle arrest that regulates a large number of genes, but does not require signaling through cGMP. In humans, antiproliferative effects of NO* may rely substantially on cGMP-independent mechanisms. Stress kinase signaling and alterations in mRNA stability appear to be major pathways by which NO* regulates the transcriptome.

Blotting, Western↗

PubMed2005

Endosperm-preferred expression of maize genes as revealed by transcriptome-wide analysis of expressed sequence tags.

The transcriptome-wide endosperm-preferred expression of maize genes was addressed by analyzing a large database of expressed sequence tags (ESTs). We generated 30,531 high quality sequence-reads from the 5'-ends of cDNA libraries from maize endosperm harvested at 10, 15, and 20 days after pollination. A further 196,900 maize sequence-reads retrieved from public databases were added to this endosperm collection to generate MAIZEST, a database with tools for data storage and analysis. MAIZEST contains 227,431 ESTs, one third of which represents developing endosperm and the remaining two-thirds represent transcripts from 49 cDNA libraries constructed from different organs and tissues. Assembling the MAIZEST ESTs generated 29,206 putative transcripts, of which a set of 4032 assembled sequences was composed exclusively of sequences derived from endosperm cDNA libraries. After sequence analysis using overlapping parameters, a sub-set of 2403 assembled sequences was functionally annotated and revealed a wide variety of putative new genes involved in endosperm development and metabolism.

Expressed Sequence Tags↗

PubMed2005

Combined transcriptome profiling reveals a novel family of arbuscular mycorrhizal-specific Medicago truncatula lectin genes.

The large majority of plants are capable of undergoing a tight symbiosis with arbuscular mycorrhizal (AM) fungi. During this symbiosis, highly specialized new structures called arbuscules are formed within the host cells, indicating that, during interaction with AM fungi, plants express AM-specific genetic programs. Despite increasing efforts, the number of genes known to be induced in the AM symbiosis is still low. In order to identify novel AM-induced genes which have not been listed before, 5,646 expressed sequence tags (ESTs) were generated from two Medicago truncatula cDNA libraries: a random cDNA library (MtAmp) and a suppression subtractive hybridization (SSH) library (MtGim), the latter being designed to enhance the cloning of mycorrhiza-upregulated genes. In silico expression analysis was applied to identify those tentative consensus sequences (TCs) of The Institute for Genomic Research M. truncatula gene index (MtGI) that are composed exclusively of ESTs deriving from the MtGim or MtAmp library, but not from any other cDNA library of the MtGI. This search revealed 115 MtAmp- or MTGim-specific TCs. For the majority of these TCs with sequence similarities to plant genes, the AM-specific expression was verified by quantitative reverse-transcription polymerase chain reaction. Annotation of the novel genes induced in mycorrhizal roots suggested their involvement in different transport as well as signaling processes and revealed a novel family of AM-specific lectin genes. The expression of reporter gene fusions in transgenic roots revealed an arbuscule-related expression of two members of the lectin gene family, indicating a role for AM-specific lectins during arbuscule formation or functioning.

Amino Acid Sequence↗

PubMed2026

Exploring the Effect of Whole-Genome Duplication on Salmonid LincRNA Repertoire.

Long intergenic non-coding RNAs (lincRNAs) are key epigenetic regulators of genome function, yet their evolutionary dynamics following whole-genome duplication (WGD) events remain poorly understood. Salmonids, which underwent a lineage-specific autotetraploidization (salmonid-specific WGD, ~88-100 million years ago), provide an excellent model to investigate the retention, divergence, and functional potential of recently duplicated non-coding elements. LincRNA repertoires were compared across five genome-annotated salmonids (Oncorhynchus tshawytscha, O. kisutch, O. mykiss, Salmo salar, and S. trutta) and their closest non-duplicated relative, northern pike (Esox lucius). LincRNAs represented ~5-7% of annotated genes in all salmonids except S. salar (18%). Sequence conservation was low relative to coding genes, with only 11-68 highly similar (e-value < 1 × 10-30; similarity > 70% and alignments > 100 nucleotides) putative orthologues shared between salmonids and northern pike, and 161-338 among salmonids alone. Synteny conservation was modest in lincRNAs, with lower conservation in putative orthologues (8-16%) compared to putative ohnologues (8-33%). Secondary structure conservation was associated with sequence similarity (ρ = -0.45; p = 2.2 × 10-16), and the association was stronger among WGD ohnologues than orthologues. In S. salar and O. mykiss, lincRNA putative ohnologues showed weaker expression correlations than coding genes, suggesting widespread regulatory divergence, possibly through neo- and subfunctionalisation. Conserved salmonid lincRNAs showed enriched predicted interactions with miRNAs involved in tumour suppression, brain, bone, and muscle development (e.g., miR-455, miR-365, miR124, miR-133a, miR-140, and miR-9), a finding supported by limited transcriptomic data. Although salmonid WGD expanded lincRNA repertoires, lincRNAs have undergone rapid sequence and transcriptional divergence, with limited conservation across species based on sequence similarity, chromosomal position, synteny, and secondary structure. A subset of conserved lincRNAs retains structural features and regulatory signatures consistent with roles as miRNA sponges in brain, skeletal, and muscle development and tumour suppression, potentially acting within conserved regulatory networks. These findings provide new insights into lincRNA evolution following genome duplication and highlight the need for experimental validation of their regulatory functions.

Animals↗

PubMed2006

Deciphering B-ZIP transcription factor interactions in vitro and in vivo.

Over the last 15 years, numerous studies have addressed the structural rules that regulate dimerization stability and dimerization specificity of the leucine zipper, a dimeric parallel coiled-coil domain that can either homodimerize or heterodimerize. Initially, these studies were performed with a limited set of B-ZIP proteins, sequence-specific DNA binding proteins that dimerize using the leucine zipper domain to bind DNA. A global analysis of B-ZIP leucine zipper dimerization properties can be rationalized using a limited number of structural rules [J.R. Newman, A.E. Keating, Comprehensive identification of human bZIP interactions with coiled-coil arrays, Science 300 (2003) 2097-2101]. Today, however, access to the genomic sequences of many different organisms has made possible the annotation of all B-ZIP proteins from several species and has generated a bank of data that can be used to refine, and potentially expand, these rules. Already, a comparative analysis of the B-ZIP proteins from Arabidopsis thaliana and Homo sapiens has revealed that the same amino acids are used in different patterns to generate diverse B-ZIP dimerization patterns [C.D. Deppmann, A. Acharya, V. Rishi, B. Wobbes, S. Smeekens, E.J. Taparowsky, C. Vinson, Dimerization specificity of all 67 B-ZIP motifs in Arabidopsis thaliana: a comparison to Homo sapiens B-ZIP motifs, Nucleic Acids Res. 32 (2004) 3435-3445]. The challenge ahead is to investigate the biological significance of different B-ZIP protein-protein interactions. Gaining insight at this level will rely on ongoing investigations to (a) define the role of target DNA on modulating B-ZIP dimerization partners, (b) characterize the B-ZIP transcriptome in various cells and tissues through mRNA microarray analysis, (c) identify the genomic localization of B-ZIP binding at a genomic level using the chromatin immunoprecipitation assay, and (d) develop more sophisticated imaging technologies to visualize dimer dynamics in single cells and whole organisms. Studies of B-ZIP family leucine zipper dimerization and the regulatory mechanisms that control their biological activities could serve as a paradigm for deciphering the biophysical and biological parameters governing other well-characterized protein-protein interaction motifs. This review will focus on the dimerization specificity of coiled-coil proteins, particularly the human B-ZIP transcription family that consists of 53 proteins that use the leucine zipper coiled-coil as a dimerization motif.

Amino Acid Motifs↗

PubMed2006

Differential use of signal peptides and membrane domains is a common occurrence in the protein output of transcriptional units.

Membrane organization describes the orientation of a protein with respect to the membrane and can be determined by the presence, or absence, and organization within the protein sequence of two features: endoplasmic reticulum signal peptides and alpha-helical transmembrane domains. These features allow protein sequences to be classified into one of five membrane organization categories: soluble intracellular proteins, soluble secreted proteins, type I membrane proteins, type II membrane proteins, and multi-spanning membrane proteins. Generation of protein isoforms with variable membrane organizations can change a protein's subcellular localization or association with the membrane. Application of MemO, a membrane organization annotation pipeline, to the FANTOM3 Isoform Protein Sequence mouse protein set revealed that within the 8,032 transcriptional units (TUs) with multiple protein isoforms, 573 had variation in their use of signal peptides, 1,527 had variation in their use of transmembrane domains, and 615 generated protein isoforms from distinct membrane organization classes. The mechanisms underlying these transcript variations were analyzed. While TUs were identified encoding all pairwise combinations of membrane organization categories, the most common was conversion of membrane proteins to soluble proteins. Observed within our high-confidence set were 156 TUs predicted to generate both extracellular soluble and membrane proteins, and 217 TUs generating both intracellular soluble and membrane proteins. The differential use of endoplasmic reticulum signal peptides and transmembrane domains is a common occurrence within the variable protein output of TUs. The generation of protein isoforms that are targeted to multiple subcellular locations represents a major functional consequence of transcript variation within the mouse transcriptome.

Animals↗

PubMed2004

The transcriptome of adult female Anopheles darlingi salivary glands.

Anopheles (Nyssorhynchus) darlingi is an important malaria vector in South and Central America; however, little is known about molecular aspects of its biology. Genomic and proteomic analyses were performed on the salivary gland products of Anopheles darlingi. A total of 593 randomly selected, salivary gland-derived cDNAs were sequenced and assembled based on their similarities into 288 clusters. The putative translated proteins were classified into three categories: (S) secretory products, (H) housekeeping products and (U) products with unknown cell location and function. Ninety-three clusters encode putative secreted proteins and several of them, such as an anophelin, a thrombin inhibitor, apyrases and several new members of the D7 protein family, were identified as molecules involved in haematophagy. Sugar-feeding related enzymes (alpha-glucosidases and alpha-amylase) also were found among the secreted salivary products. Ninety-nine clusters encode housekeeping proteins associated with energy metabolism, protein synthesis, signal transduction and other cellular functions. Ninety-seven clusters encode proteins with no similarity with known proteins. Comparison of the sequence divergence of the S and H categories of proteins of An. darlingi and An. gambiae revealed that the salivary proteins are less conserved than the housekeeping proteins, and therefore are changing at a faster evolutionary rate. Tabular and supplementary material containing the cDNA sequences and annotations are available at http://www.ncbi.nlm.nih.gov/projects/Mosquito/A_darlingi_sialome/

Animals↗

PubMed2006

Differential brain transcriptome of beta4 nAChR subunit-deficient mice: is it the effect of the null mutation or the background strain?

Studies using mice with beta4 nicotinic acetylcholine receptor (nAChR) subunit deficiency (beta4-/- mice) helped reveal the roles of this subunit in bradycardiac response to vagal stimulation, nicotine-induced seizure activity and anxiety. To identify genes that might be related to beta4-containing nAChRs activity, we compared the mRNA expression profiles of brains from beta4-/- and wild-type mice using Affymetrix U74Av2 microarray. Seventy-seven genes significantly differentiated between these two experimental groups. Of them, the two most downregulated were spastic paraplegia 21 (human) homolog (Spg21) and 6-pyruvoyl-tetrahydropterin synthase (Pts) genes. Since the targeted mutagenesis of the beta4 nAChR subunit was done by using two mouse strains, 129SvEv and C57BL/6J, it is possible that the genes closely linked to the mutated beta4 gene represent the 129SvEv allele and not the control C57BL/6J-driven allele. We examined this possibility by using public database and quantitative RT-PCR. The expression levels of Spg21 and Pts genes that, like the beta4 gene, are localized on mouse chromosome 9, as well as the expression levels of other genes located on this chromosome, were dependent on the mouse background strain. The 67 differentially expressed genes that are not located on chromosome 9 were further analyzed for overrepresented functional annotations and transcription regulatory elements compared with the entire microarray. Genes encoding for proteins involved in tyrosine phosphatase activity, calcium ion binding, cell growth and/or maintenance, and chromosome organization were overrepresented. Our data enhance the understanding of the molecular interactions involved in the beta4 nAChR subunit function. They also emphasize the need for careful interpretation of expression microarray studies done on genetically manipulated animals.

Animals↗

PubMed2006

Integrative genomics identifies distinct molecular classes of neuroblastoma and shows that multiple genes are targeted by regional alterations in DNA copy number.

Neuroblastoma is remarkable for its clinical heterogeneity and is characterized by genomic alterations that are strongly correlated with tumor behavior. The specific genes that influence neuroblastoma biology and are targeted by genomic alterations remain largely unknown. We quantified mRNA expression in a highly annotated series of 101 prospectively collected diagnostic neuroblastoma primary tumors using an oligonucleotide-based microarray. Genomic copy number status at the prognostically relevant loci 1p36, 2p24 (MYCN), 11q23, and 17q23 was determined by PCR and was aberrant in 26, 20, 40, and 38 cases, respectively. In addition, 72 diagnostic neuroblastoma primary tumors assayed in a different laboratory were used as an independent validation set. Unsupervised hierarchical clustering showed that gene expression was highly correlated with genomic alterations and clinical markers of tumor behavior. The vast majority of samples with MYCN amplification and 1p36 loss of heterozygosity (LOH) clustered together on a terminal node of the sample dendrogram, whereas the majority of samples with 11q deletion clustered separately and both of these were largely distinct from the copy number neutral group of tumors. Genes involved in neurodevelopment were broadly overrepresented in the more benign tumors, whereas genes involved in RNA processing and cellular proliferation were highly represented in the most malignant cases. By combining transcriptomic and genomic data, we showed that LOH at 1p and 11q was associated with significantly decreased expression of 122 (61%) and 88 (27%) of the genes mapping to 1p35-36 and all of 11q, respectively, suggesting that multiple genes may be targeted by LOH events. A total of 71 of the 1p35-36 genes were also differentially expressed in the independent validation data set, providing a prioritized list of candidate neuroblastoma suppressor genes. Taken together, these data are consistent with the hypotheses that the neuroblastoma transcriptome is a sensitive marker of underlying tumor biology and that chromosomal deletion events in this cancer likely target multiple genes through alteration in mRNA dosage. Lead positional candidates for neuroblastoma suppressor genes can be inferred from these data, but the potential multiplicity of transcripts involved has significant implications for ongoing gene discovery strategies.

Chromosome Aberrations↗

PubMed2003

A four-tiered transcriptional regulatory circuit controls flagellar biogenesis in Pseudomonas aeruginosa.

The single polar flagellum of Pseudomonas aeruginosa is an important virulence and colonization factor of this opportunistic pathogen. In this study, the annotation of the genes belonging to the fla regulon was updated and their organization was analysed in strains PAK and PAO1, representative type-a and type-b strains of P. aeruginosa respectively. The flagellar genes are clustered in three non-contiguous regions of the chromosome. A polymorphic locus flanked by flgJ and fleQ in Region I contains a glycosylation island in PAK. The expression and ordered assembly of the complex multicomponent flagellum is intricately regulated. Dedicated flagellar genes fleQ, fleS, fleR, fliA, flgM and fleN encode proteins that participate in the regulation of the flagellar transcriptional circuit. In addition, expression of the flagellum is coordinately regulated with other P. aeruginosa virulence factors by the alternative sigma factor sigma54, encoded by rpoN. In order to gain insight into the hierarchical regulation of flagellar genes, deletion mutations were constructed in fleQ, fleR, fliA and rpoN. The transcriptional impact of these mutations was examined by transcriptional profiling using a P. aeruginosa whole genome microarray. Analysis of the transcriptomes generated for each of these mutants indicates a four-tiered (Classes I-IV) hierarchy of transcriptional regulation. Class I genes are constitutively expressed and include the transcriptional regulator fleQ and the alternative sigma factor fliA (sigma28). Class II genes including fleSR, encoding a two-component regulatory system require FleQ and RpoN (sigma54) for their transcriptional activation. Class III genes are positively regulated by the activated response regulator FleR in concert with RpoN. The transcription of Class IV genes is dependent on the availability of free FliA following the export of the FliA specific antisigma factor FlgM through the basal body rod-hook structure (assembled from Class II and III gene products). Two previously uncharacterized genes, which are coordinately regulated with known flagellar genes have been identified by genome-wide analysis and their role in flagellar biogenesis was analysed.

Bacterial Proteins↗

PubMed2026

The current and future perspective of ChickenGTEx project and its applications in precision breeding.

The Chicken Genotype-Tissue Expression (ChickenGTEx) project was established to systematically characterize the regulatory landscape of the chicken genome and to accelerate the translation of functional genomics into precision breeding. By integrating whole-genome sequencing with multi-tissue transcriptomic profiling, ChickenGTEx provides a comprehensive atlas of gene expression regulation across diverse tissues and physiological systems. Current findings demonstrate that complex production traits are governed by coordinated regulatory networks rather than isolated loci, with substantial contributions from tissue-specific gene expression, structural variation, and genotype-by-sex interactions. Sex-dependent regulatory effects further refine the genetic architecture of metabolic, immune, and reproductive traits, highlighting the importance of incorporating sex as a biological variable in genomic analyses. Application of integrative omics frameworks within elite layer populations has revealed multilayer regulatory mechanisms underlying extended laying performance, feed efficiency, metabolic health, and eggshell quality. By partitioning phenotypic variance into genetic, regulatory, and host-microbiome components, these approaches move beyond association-based mapping toward causal inference and biological interpretation. Importantly, validated regulatory loci identified through ChickenGTEx and related analyses provide actionable markers for genomic selection and rational targets for precision genome modification. Looking forward, continued expansion of regulatory atlases, incorporation of single-cell and longitudinal data in diverse environmental conditions, and integration of functional annotation into breeding pipelines will further enhance prediction accuracy and sustainable genetic improvement. The ChickenGTEx project thus represents a foundational platform bridging functional genomics and practical poultry breeding.

Animals↗

PubMed2026

Isoform-Level Analysis Reveals Reproducible Early Changes in Transcript Usage During Human Vaccine Responses.

Vaccine-induced transcriptional responses have been extensively characterized at the gene level, but whether vaccination also alters transcript isoform usage remains largely unexplored. Here, we reanalyzed longitudinal whole-blood RNA-seq data from a discovery cohort of mRNA COVID-19 vaccine recipients using the IsoformSwitchAnalyzeR framework and validated the findings in an independent cohort. Key findings were validated by full-length RNA long-read sequencing and extended to four additional vaccine cohorts covering distinct platforms and pathogens. mRNA vaccination induced a rapid and transient wave of differential transcript usage, peaking at 24 h post-vaccination with 131 isoforms significantly altered across 107 genes, before largely resolving by Day 14. Isoform switching events were reproducible across independent cohorts and confirmed by full-length RNA long-read sequencing. Structural annotation of switching transcripts, including RMI2, WARS1, and NT5C3A, revealed changes affecting predicted protein domains and signal peptides. Notably, highly concordant isoform switching patterns were observed across MVA-based SARS-CoV-2, influenza, and Ebola vaccine cohorts and showed dose-dependent modulation. Overall, differential transcript isoform usage is a rapid and transient feature of the early human immune response to vaccination that was observed across multiple vaccine platforms. These findings reveal an underappreciated layer of transcriptional regulation that complements conventional gene-level analyses and warrants integration into future vaccine immunogenicity studies.

Humans↗

PubMed2026

Exploring genomic regions regulating the liver transcriptome and energy homeostasis in pigs.

In pigs, energy homeostasis has an impact on meat quality and health. In a Duroc pig population, 30 quantitative trait locus (QTL) regions associated with fatty acid (FA) composition in adipose tissue, plasma, liver and muscle were previously identified. Mapping of expression quantitative trait locus (eQTL) regions will provide a molecular hypothesis for genotype-phenotype interactions and may allow the identification of shared causal variants, key to increasing our understanding of the genetic regulation of FA composition and energy homeostasis. However, gene expression is impacted by environmental factors, while individual-level allelic imbalance (AI) can be more reliable and can be surveyed via allelic-specific expression (ASE) analysis. Furthermore, treatment of ASE as a quantitative trait allows the identification of allele-specific expression quantitative trait loci (aseQTLs), which are variants whose heterozygosity is linked to the AI of a nearby single-nucleotide polymorphism (SNP), pointing to regulatory elements. In this study, liver was selected as a key metabolic hub with an important role in the regulation of energy homeostasis, and 310 liver RNA sequencing samples were analysed using a combination of (1) eQTL mapping, (2) ASE analysis, and (3) aseQTL mapping methods. A total of 2 188 eQTL regions were identified, mostly cis-eQTL regions (73.17%). ASE analysis reported 1 964 ASE SNPs, associated with 633 genes. Finally, aseQTL mapping reported 64 172 aseQTL, associated with the AI of 31 genes. Colocalisation analysis combined with ASE analysis showed that the expression of FADS1 and FADS2 genes is associated with the polyunsaturated FA composition in several tissues, where microRNA regulation may be present. Finally, in the DGAT2 gene, annotated in a QTL region associated with multiple FAs in adipose tissue, ASE revealed allelic imbalance in the 3' untranslated region (UTR) of this gene. Allelic differential expression can be caused by a 13-bp insertion affecting messenger RNA stability, previously described, exemplifying how allelic imbalance is caused by a post-transcriptional regulatory mechanism undetectable by eQTL mapping. Furthermore, aseQTLs were associated with this gene, linked to a previously identified copy number variant not yet associated with DGAT2 expression. These results demonstrate how ASE analysis and aseQTL mapping can complement eQTL mapping, as they resolved a complex region affected by allelic heterogeneity, a main confounding effect of QTL mapping. In conclusion, the combination of eQTL mapping, ASE analysis and aseQTL mapping allowed the characterisation of the regulation of liver gene expression, improving our understanding of the genetic determinism of energy homeostasis.

Allele-specific expression↗