PubMed Health⌕ Search

SEARCH · PubMed Health

Results for “spatial genomics”

Explore indexed PubMed citations for clinical trials, systematic reviews and public health research. Read source abstracts and follow each citation to its original PubMed record.

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

At least 1,369 records · Page 76Linked to original sources

Perseus: Lineage-Aware Refinement of Kraken2 Taxonomic Classification for Long Read Metagenomes.

MOTIVATION: Long-read metagenomic sequencing improves assembly contiguity and enables genome-resolved analysis of complex microbial communities, but accurate taxonomic classification of long reads and assembled contigs remains challenging. Highly scalable k-mer-based classifiers such as Kraken2 frequently over-assign fine-rank taxonomic labels when applied to long-read data, producing high false positive classification rates driven by sparse or localized k-mer matches, particularly in microbiomes with extensive taxonomic novelty. RESULTS: We present Perseus, a lineage-aware confidence estimation framework for taxonomic classification that models the spatial distribution and hierarchical consistency of k-mer evidence along sequences. This formulation reframes taxonomic classification as a hierarchical confidence estimation problem rather than a single-rank prediction task. Perseus refines k-mer-level taxonomic signals from Kraken2 using a multi-headed convolutional neural network that estimates calibrated confidence scores for taxonomic correctness at each canonical rank. Using these estimates, Perseus confirms assignments, backs off to higher taxonomic ranks, or abstains when evidence is insufficient, prioritizing correctness and lineage consistency over overly specific assignments. Across simulations of taxonomic novelty and real-world metagenomic datasets, Perseus consistently and substantially reduces the false assignment rate while improving precision and lineage-consistent accuracy. These improvements are most pronounced for long reads and assembled contigs, where spatial context enables reliable discrimination between consistent taxonomic signal and spurious matches. AVAILABILITY AND IMPLEMENTATION: Perseus integrates with existing Kraken2 workflows and is available at https://github.com/matnguyen/perseus.

Journal Article↗

The influence of maternal lineages on social affiliations among humpback whales (Megaptera novaeangliae) on their feeding grounds in the southern gulf of Maine.

Humpback whales on their feeding grounds in the Gulf of Maine typically form fluid fission/fusion groups of two to three individuals characterized by noncompetitive and, at times, cooperative behavior. Here we test the hypothesis that, despite the apparent absence of close kinship bonds, the fluid associations between feeding whales are influenced by "maternal lineages" as represented by mtDNA haplotypes. Using skin samples collected with a biopsy dart, variation in the hypervariable segment of the mtDNA control region identified 17 unique haplotypes among 159 individually identified whales from the southern Gulf of Maine. The haplotypes of a further 143 individuals were inferred from known direct maternal (cow-calf) relationships. The frequencies of associations among these 302 individuals were calculated from 21,617 sighting records collected from 1980 to 1995, excluding associations between a cow and her dependent calf. For groups of two where the haplotypes of both individuals were known (n = 3,151), individuals with the same haplotype were together significantly more often (26%) than expected by random association (20%). To account for different group sizes and associations with individuals of unknown haplotype and sex, we used Monte Carlo simulations to test for nonrandom associations in the full data set, as well as known female-only (n = 1,512), male-only (n = 730), and mixed-sex (n = 2,745) groups. Within-haplotype associations were significantly more frequent than expected at random for all groups (P = .002) and female-only groups (P = .011) but not male-only groups, while mixed-sex groups approached significance (P = .062). A Mantel test of individual pairwise association indices and haplotype identity confirmed that within-haplotype associations were more frequent than expected for all sex combinations except male-male associations, with females forming within-haplotype associations 1.7 times more often than expected by random assortment. Partial matrix correlations and permutation analyses indicated that the skew toward within-haplotype associations could not be accounted for by short-term temporal co-occurrence or fine-scale spatial distributions of individuals with shared haplotypes. While the mechanism by which individuals with a common mtDNA haplotype assort remains unknown, our results strongly suggest an influence of maternal lineages on the social organization of humpback whales within a regional feeding ground.

Animals↗

Copy-dependent and position-independent expression of rat aldolase A gene.

In order to understand the molecular mechanisms of the temporal and spatial differences of gene expression in higher organisms, rat aldolase A gene carrying two distinct promoters was introduced into fertilized eggs and the resulting transgenic mice were analyzed. The transgene expression is tissue-specific and is developmentally regulated. In addition, the expression is regulated in a copy-dependent manner irrespective of where the transgene is integrated, suggesting that a mechanism excluding the effect of the integration site exists within the transgene itself. To explore the conformational change of this gene in the genome, the DNase I hypersensitive sites of the gene were examined. Three sites (DHS-1,2, and 3) were identified upstream and downstream of the gene and these sites were retained in the transgene as well as in the gene observed endogenously.

Animals↗

Olfactory coding in a compound nose. Coexpression of odorant-binding proteins in Drosophila.

Odorant-binding proteins (OBPs) are small, soluble proteins present in the aqueous medium surrounding olfactory receptor neurons. Their function in olfaction is unknown: they have been proposed to facilitate the transit of hydrophobic molecules to olfactory receptors, to deactivate the odorant stimulus, and/or to play a role in chemosensory coding. We have examined the genomic organization and expression patterns of two olfactory-specific genes (OS-E and OS-F) of Drosophila melanogaster, the products of which are members of a protein family in Drosophila sharing sequence similarity with moth OBPs. We found that the OS-E and OS-F transcription units are located < 1 kb apart. They are oriented in the same direction and display a similar intron-exon organization. Expression of both OS-E and OS-F proteins is spatially restricted to the ventrolateral region of the Drosophila antenna. Within this region, both OS-E and OS-F proteins are expressed within two different types of sensory hairs: in most, if not all, sensilla trichodea and in approximately 40% of the interspersed small sensilla basiconica. We consistently observe that OS-E and OS-F are coexpressed, indicating that an individual sensillum can contain more than one odorant-binding protein. This finding has potential implications for the roles of odorant-binding proteins in olfactory coding.

Animals↗

Structural organization of poliovirus RNA replication is mediated by viral proteins of the P2 genomic region.

Transcriptionally active replication complexes bound to smooth membrane vesicles were isolated from poliovirus-infected cells. In electron microscopic, negatively stained preparations, the replication complex appeared as an irregularly shaped, oblong structure attached to several virus-induced vesicles of a rosettelike arrangement. Electron microscopic immunocytochemistry of such preparations demonstrated that the poliovirus replication complex contains the proteins coded by the P2 genomic region (P2 proteins) in a membrane-associated form. In addition, the P2 proteins are also associated with viral RNA, and they can be cross-linked to viral RNA by UV irradiation. Guanidine hydrochloride prevented the P2 proteins from becoming membrane bound but did not change their association with viral RNA. The findings allow the conclusion that the protein 2C or 2C-containing precursor(s) is responsible for the attachment of the viral RNA to the vesicular membrane and for the spatial organization of the replication complex necessary for its proper functioning in viral transcription. A model for the structure of the viral replication complex and for the function of the 2C-containing P2 protein(s) and the vesicular membranes is proposed.

Antibodies, Monoclonal↗

Multiple mechanisms promote the retained expression of gene duplicates in the tetraploid frog Xenopus laevis.

Gene duplication provides a window of opportunity for biological variants to persist under the protection of a co-expressed copy with similar or redundant function. Duplication catalyzes innovation (neofunctionalization), subfunction degeneration (subfunctionalization), and genetic buffering (redundancy), and the genetic survival of each paralog is triggered by mechanisms that add, compromise, or do not alter protein function. We tested the applicability of three types of mechanisms for promoting the retained expression of duplicated genes in 290 expressed paralogs of the tetraploid clawed frog, Xenopus laevis. Tests were based on explicit expectations concerning the ka/ks ratio, and the number and location of nonsynonymous substitutions after duplication. Functional constraints on the majority of paralogs are not significantly different from a singleton ortholog. However, we recover strong support that some of them have an asymmetric rate of nonsynonymous substitution: 6% match predictions of the neofunctionalization hypothesis in that (1) each paralog accumulated nonsynonymous substitutions at a significantly different rate and (2) the one that evolves faster has a higher ka/ks ratio than the other paralog and than a singleton ortholog. Fewer paralogs (3%) exhibit a complementary pattern of substitution at the protein level that is predicted by enhancement or degradation of different functional domains, and the remaining 13% have a higher average ka/ks ratio in both paralogs that is consistent with altered functional constraints, diversifying selection, or activity-reducing mutations after duplication. We estimate that these paralogs have been retained since they originated by genome duplication between 21 and 41 million years ago. Multiple mechanisms operate to promote the retained expression of duplicates in the same genome, in genes in the same functional class, over the same period of time following duplication, and sometimes in the same pair of paralogs. None of these paralogs are superfluous; degradation or enhancement of different protein subfunctions and neofunctionalization are plausible hypotheses for the retained expression of some of them. Evolution of most X. laevis paralogs, however, is consistent with retained expression via mechanisms that do not radically alter functional constraints, such as selection to preserve post-duplication stoichiometry or temporal, quantitative, or spatial subfunctionalization.

Amino Acid Substitution↗

Activin A induced expression of a fork head related gene in posterior chordamesoderm (notochord) of Xenopus laevis embryos.

A gene family encoding the fork head/HNF-3 domain has been identified in the South African clawed frog, Xenopus laevis. Screening of genomic DNA and gastrula stage derived cDNA libraries with a DNA fragment encoding the Drosophila fork head domain led to the isolation of a number of different clones encoding this motif. While one of the Xenopus fork head sequences, XFD-3, represents the Xenopus counterpart to rat HNF-3 beta, all other sequences encode novel types of fork head related proteins. Here we report on XFD-1, a DNA binding protein which can bind to the HNF-3 alpha target sequence. Analysis of temporal and spatial expression revealed that the gene is activated at blastula stage and that transcripts are localized in a rather thin stripe of cells invaginating the dorsal blastopore lip (organizer) during gastrulation. XFD-1 mRNA is localized within the notochord and, by the end of neurulation, is no longer detectable. In the animal cap assay the gene is activated by incubation with the vegetalizing factor (activin A) but not with bFGF.

Activins↗

A scale invariant clustering of genes on human chromosome 7.

BACKGROUND: Vertebrate genes often appear to cluster within the background of nontranscribed genomic DNA. Here an analysis of the physical distribution of gene structures on human chromosome 7 was performed to confirm the presence of clustering, and to elucidate possible underlying statistical and biological mechanisms. RESULTS: Clustering of genes was confirmed by virtue of a variance of the number of genes per unit physical length that exceeded the respective mean. Further evidence for clustering came from a power function relationship between the variance and mean that possessed an exponent of 1.51. This power function implied that the spatial distribution of genes on chromosome 7 was scale invariant, and that the underlying statistical distribution had a Poisson-gamma (PG) form. A PG distribution for the spatial scattering of genes was validated by stringent comparisons of both the predicted variance to mean power function and its cumulative distribution function to data derived from chromosome 7. CONCLUSION: The PG distribution was consistent with at least two different biological models: In the microrearrangement model, the number of genes per unit length of chromosome represented the contribution of a random number of smaller chromosomal segments that had originated by random breakage and reconstruction of more primitive chromosomes. Each of these smaller segments would have necessarily contained (on average) a gamma distributed number of genes. In the gene cluster model, genes would be scattered randomly to begin with. Over evolutionary timescales, tandem duplication, mutation, insertion, deletion and rearrangement could act at these gene sites through a stochastic birth death and immigration process to yield a PG distribution. On the basis of the gene position data alone it was not possible to identify the biological model which best explained the observed clustering. However, the underlying PG statistical model implicated neutral evolutionary mechanisms as the basis for this clustering.

Algorithms↗

Coexpression of two odorant-binding protein homologs in Drosophila: implications for olfactory coding.

Odorant-binding proteins (OBPs) are small soluble proteins present in the aqueous medium surrounding olfactory receptor neurons. Their function in olfaction is still unknown: they have been proposed to facilitate the transit of hydrophobic molecules to olfactory receptors, to deactivate the odorant stimulus, and/or to play a role in chemosensory coding. In this study we examine the genomic organization and expression patterns of two olfactory-specific genes (OS-E and OS-F) of Drosophila melanogaster, the products of which are members of a protein family in Drosophila sharing sequence similarity with moth OBPs. We show that the OS-E and OS-F transcription units are located <1 kb apart. They are oriented in the same direction and display a similar intron-exon organization. Expression of both OS-E and OS-F proteins is restricted spatially to the ventrolateral region of the Drosophila antenna. Within this region both OS-E and OS-F proteins are expressed within two different types of sensory hairs: in most, if not all, sensilla trichodea and in approximately 40% of the interspersed small sensilla basiconica. We consistently observe that OS-E and OS-F are coexpressed, indicating that an individual sensillum can contain more than one odorant-binding protein. The functional significance of the observed expression pattern and its implications for olfactory coding are discussed.

Animals↗

Colorectal Liver Metastasis Pathomics Model: Integrating Single-Cell and Spatial Transcriptome Analysis With Pathomics for Predicting Liver Metastasis in Colorectal Cancer.

The liver is the primary target organ for hematologic metastasis of colorectal cancer (CRC), and CRC liver metastasis (CRLM) often precludes radical resection, making it the leading cause of death in patients with CRC. To improve the identification and prediction of liver metastasis risk, we identified a cell type of liver metastasis--triggering malignant cells (LMTMCs) through integrating single-cell RNA sequencing and spatial transcriptome analysis. Multiomics cell communication analysis indicated that the interaction between fibroblasts and LMTMCs through the COL1A1-CD44/SDC4 and LAMA4-CD44 signaling axes could promote CRLM. By applying the one-class logistic regression algorithm, we developed a CRLM scoring system in the bulk RNA-sequencing data according to the abundance of LMTMCs in each individual. Using the grouping labels derived from the CRLM scoring system in the bulk data and the corresponding whole-slide images without any manual annotations at the region or pixel level, processed via slide-level weakly supervised learning, a deep-learning model based on the ResNet18 architecture, called Colorectal Liver Metastasis Pathomics Model, was developed to predict the risk of liver metastasis in patients with CRC. The Colorectal Liver Metastasis Pathomics Model achieved an area under the curve of 0.84 at the internal test set of The Cancer Genome Atlas-CRC histology images. In the external independent validation sets, namely the Affiliated Hospital of Southwest Medical University and the Affiliated Traditional Chinese Medicine Hospital of Southwest Medical University cohorts, the areas under the curve were 0.89 and 0.72, respectively, indicating effective classification performances. This study provided new insights and tools for the early identification of CRLM and demonstrated the potential of combining multiomics with deep learning-based pathomics in cancer research.

Humans↗

Tissue-nonspecific alkaline phosphatase participates in the establishment and growth of feather germs in embryonic chick skin cultures.

Alkaline phosphatase activity is present in the mesoderm of embryonic chick skin and becomes spatially restricted to the dermal condensation of the developing feather germs. Inhibitors to tissue-nonspecific (liver/bone/kidney), but not intestinal alkaline phosphatase inhibit the establishment and growth of feather germs in cultured skins. A window of maximum sensitivity to the inhibitor was observed to be the first day of culture when early development and establishment of pattern takes place. The cDNA for the avian tissue-nonspecific alkaline phosphatase was cloned and sequenced, and Southern analysis revealed a single copy of this gene in the avian genome. Northern analysis revealed that a 2.8 kb transcript for this form of alkaline phosphatase is present in developing skin.

Alkaline Phosphatase↗

Histone lysine methylation patterns in human cell types are arranged in distinct three-dimensional nuclear zones.

The impact of histone lysine methylation as an essential epigenetic mechanism for gene regulation has been demonstrated by numerous studies where it was functionally and structurally linked to euchromatin and heterochromatin. Most of these data have been obtained by biochemical and two-dimensional (2D)-microscopic techniques providing little information about the global nuclear arrangement of histone modifications. We investigated the 3D architecture and spatial interrelationships of different histone lysine methylation sites (tri-H3K4, mono-H4K20, mono-H3K9, tri-H3K27, tri-H4K20 and tri-H3K9) in various human cell types. Immunofluorescence and confocal microscopy were used together with a quantitative evaluation of 3D images, to reveal spatial relations of specific methylation sites with either centromeres, nascent RNA or with each other. A close association with centromeres was found only for histone methylation sites previously linked to constitutively repressed chromatin. Differences observed in these sites in relation to the cell cycle emphasize the potential relevance of the dynamic properties of heterochromatin for nuclear functions. Nascent RNA was found associated, though to a different degree, with all histone methylation sites, supporting the increasing evidence that transcription occurs across a wide range of the human genome. Finally we demonstrated by simultaneous visualization of different histone lysine methylation sites that methylation patterns are organized in distinct nuclear zones with little apparent intermingling.

Algorithms↗

Compositional heterogeneity within, and uniformity between, DNA sequences of yeast chromosomes.

The heterogeneity within, and similarities between, yeast chromosomes are studied. For the former, we show by the size distribution of domains, coding density, size distribution of open reading frames, spatial power spectra, and deviation from binomial distribution for C + G% in large moving windows that there is a strong deviation of the yeast sequences from random sequences. For the latter, not only do we graphically illustrate the similarity for the above mentioned statistics, but we also carry out a rigorous analysis of variance (ANOVA) test. The hypothesis that all yeast chromosomes are similar cannot be rejected by this test. We examine the two possible explanations of this interchromosomal uniformity: a common origin, such as genome-wide duplication (polyploidization), and a concerted evolutionary process.

Analysis of Variance↗

Spatial control of transgene expression in rice (Oryza sativa L.) using the GAL4 enhancer trapping system.

We used enhancer trapping with the GAL4 transcriptional activator from yeast to obtain spatial control of transgene expression in all organs of the model monocotyledonous species rice (Oryza sativa L. cv. Nipponbare). Our T-DNA enhancer trapping cassette consisted of two principle components: (1) the minimal promoter-equipped gal4 gene placed adjacent to the right border, and (2) the green fluorescent protein gene (gfp) fused to the upstream activation sequence element (UAS) to which GAL4 binds and activates expression, so that gfp expression corresponds to gal4 expression. Agrobacterium-mediated integration of the cassette into the rice genome often brings the gal4 gene under transcriptional control of local genomic enhancers and promoters, resulting in gal4/gfp expression patterns ranging in specificity from single-cell types to constitutive expression. We produced more than 13 000 enhancer trap lines with this cassette and screened T(0) adult plants (1982 lines), T(1) seed (2684 lines) and T(1) seedlings (2667 lines) for gfp expression. Approximately 30% of the lines produced GFP, and we identified lines with gfp expression in specific cell types of all major organs of the rice plant. Subsequently, using the GUS reporter gene (uidA), we demonstrated that UAS:geneX constructs can be transactivated in specific cell types where gal4 and gfp are expressed, thus providing an excellent system for the manipulation of gene expression and physiological function in specific cell types of rice.

DNA-Binding Proteins↗

No statistical support for correlation between the positions of protein interaction sites and alternatively spliced regions.

BACKGROUND: Alternative splicing is an efficient mechanism for increasing the variety of functions fulfilled by proteins in a living cell. It has been previously demonstrated that alternatively spliced regions often comprise functionally important and conserved sequence motifs. The objective of this work was to test the hypothesis that alternative splicing is correlated with contact regions of protein-protein interactions. RESULTS: Protein sequence spans involved in contacts with an interaction partner were delineated from atomic structures of transient interaction complexes and juxtaposed with the location of alternatively spliced regions detected by comparative genome analysis and spliced alignment. The total of 42 alternatively spliced isoforms were identified in 21 amino acid chains involved in biomolecular interactions. Using this limited dataset and a variety of sophisticated counting procedures we were not able to establish a statistically significant correlation between the positions of protein interaction sites and alternatively spliced regions. CONCLUSIONS: This finding contradicts a naïve hypothesis that alternatively spliced regions would correlate with points of contact. One possible explanation for that could be that all alternative splicing events change the spatial structure of the interacting domain to a sufficient degree to preclude interaction. This is indirectly supported by the observed lack of difference in the behaviour of relatively short regions affected by alternative splicing and cases when large portions of proteins are removed. More structural data on complexes of interacting proteins, including structures of alternative isoforms, are needed to test this conjecture.

Alternative Splicing↗

Monte Carlo simulation of the production of short DNA fragments by low-linear energy transfer radiation using higher-order DNA models.

A realistic DNA target model has been developed and implemented in the biophysical simulation code PARTRAC. It describes five levels of the B-DNA structure (nucleotides, DNA helices, nucleosomes, chromatin fiber structure and chromatin fiber loops) on an atomic level for the whole genome inside a mammalian cell nucleus. The model is capable of describing regular solenoidal, crossed-linker or zigzag structures as well as repeating stochastic arrangements of nucleosomes in the chromatin fiber. Electron tracks resulting from monoenergetic electrons with energies up to 100 keV and from 220 kVp X rays, starting at random positions in the cell, were superimposed on four DNA target models with different chromatin fiber structures. The yields of SSBs, DSBs and short single- and double-stranded DNA fragments were determined from spatial coincidences with strand atoms. Two parameters of the model-the energy necessary to create an SSB and the distance between two breaks that would be scored as a DSB-were adapted to equate simulated and measured strand break yields after X irradiation of human fibroblast cells. The integral fractions of short single- and double-stranded fragments were rather similar for all condensed chromatin fiber structures; they agreed with experimental data for DNA fragments below 2 kbp. The simulated fragment size distributions in the range from 0.1 to 1.5 kbp reflected the fiber structure irrespective of strandedness or electron energy. The distributions using a stochastic arrangement of nucleosomes in the chromatin fiber were found to be in better accordance with experimental data than those obtained with regular fiber structures.

DNA↗

The Karyote physico-chemical genomic, proteomic, metabolic cell modeling system.

Modeling approaches to the dynamics of a living cell are presented that are strongly based on its underlying physical and chemical processes and its hierarchical spatio-temporal organization. Through the inclusion of a broad spectrum of processes and a rigorous analysis of the multiple scale nature of cellular dynamics, we are attempting to advance cell modeling and its applications. The presentation focuses on our cell modeling system, which integrates data archiving and quantitative physico-chemical modeling and information theory to provide a seamless approach to the modeling/data analysis endeavor. Thereby the rapidly growing mess of genomic, proteomic, metabolic, and cell physiological data can be automatically used to develop and calibrate a predictive cell model. The discussion focuses on the Karyote cell modeling system and an introduction to the CellX and VirusX models. The Karyote software system integrates three elements: (1) a model-building and data archiving module that allows one to define a cell type to be modeled through its reaction network, structure, and transport processes as well as to choose the surrounding medium and other parameters of the phenomenon to be modeled; (2) a genomic, proteomic, metabolic cell simulator that solves the equations of metabolic reaction, transcription/translation polymerization and the exchange of molecules between parts of the cell and with the surrounding medium; and (3) an information theory module (ITM) that automates model calibration and development, and integrates a variety of data types with the cell dynamic computations. In Karyote, reactions may be fast (equilibrated) or slow (finite rate), and the special effects of enzymes and other minority species yielding steady-state cycles of arbitrary complexities are accounted for. These features of the dynamics are handled via rigorous multiple scale analysis. A user interface allows for an automated generation and solution of the equations of multiple timescale, compartmented dynamics. Karyote is based on a fixed intracellular structure. However, cell response to changes in the host medium, damage, development or transformation to abnormality can involve dramatic changes in intracellular structure. As this changes the nature of the cellular dynamics, a new model, CellX, is being developed based on the spatial distribution of concentration and other variables. This allows CellX to capture the self-organizing character of cellular behavior. The self-assembly of organelles, viruses, and other subcellular bodies is being addressed in a second new model, VirusX, that integrates molecular mechanics and continuum theory. VirusX is designed to study the influence of a host medium on viral self-assembly, structural stability, infection of a single cell, and transmission of disease.

Animals↗

Axes, boundaries and coordinates: the ABCs of fly leg development.

Recent studies of gene expression in the developing fruitfly leg support a model--Meinhardt's Boundary Model--which seems to contradict the prevailing paradigm for pattern formation in the imaginal discs of Drosophila--the Polar Coordinate Model. Reasoning from geometric first principles, this article examines the strengths and weaknesses of these hypotheses, plus some baffling phenomena that neither model can comfortably explain. The deeper question at issue is: how does the fly's genome encode the three-dimensional anatomy of the adult? Does it demarcate territories and boundaries (as in a geopolitical map) and then use those boundaries and their points of intersection as a scaffolding on which to erect the anatomy (the Boundary Model)? Or does it assign cellular fates within a relatively seamless coordinate system (the Polar Coordinate Model)? The existence of hybrid Cartesian-polar models shows that the alternatives may not be so clear-cut: a single organ might utilize different systems that are spatially superimposed or temporally sequential.

Animals↗