SEARCH · Search NASA

Results for “differential expression”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

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

At least 91 records · Page 5

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Transcriptome and secretome profiling of sensory neurons reveals sex differences in pathways relevant to insulin sensing and insulin secretion

Sensory neurons in the dorsal root ganglia (DRG) convey somatosensory and metabolic cues to the central nervous system and release substances from stimulated terminal endings in peripheral organs. Sex-biased variations driven by the sex chromosome complement (XX and XY) have been implicated in the sensory–islet crosstalk. However, the molecular underpinnings of these male–female differences are not known. Here, we aim to characterize the molecular repertoire and the secretome profile of the lower thoracic spinal sensory neurons and to identify molecules with sex-biased insulin sensing- and/or insulin secretion-modulating activity that are encoded independently of circulating gonadal sex hormones. We used transcriptomics and proteomics to uncover differentially expressed genes and secreted molecules in lower thoracic T5-12 DRG sensory neurons derived from sexually immature 3-week-old male and female C57BL/6J mice. Comparative transcriptome and proteome analyses revealed differential gene expression and protein secretion in DRG neurons in males and females. Furthermore, the transcriptome analysis identified, among others, higher insulin signaling/sensing capabilities in female DRG neurons; secretome screening uncovered several sex-specific candidate molecules with potential regulatory functions in pancreatic β cells. Together, these data suggest a putative role of sensory interoception of insulin in the DRG–islet crosstalk with implications in sensory feedback loops in the regulation of β-cell activity in a sex-biased manner. Finally, we provide a valuable resource of molecular and secretory targets that can be leveraged for understanding insulin interoception and insulin secretion and inform the development of novel studies/approaches to fathom the role of the sensory–islet axis in the regulation of energy balance in males and females.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2020

Intra‐species genetic variability drives carbon metabolism and symbiotic host interactions in the ectomycorrhizal fungus Pisolithus microcarpus

Summary Ectomycorrhizal (ECM) fungi are integral to boreal and temperate forest ecosystem functioning and nutrient cycling. ECM fungi, however, originate from diverse saprotrophic lineages and the impacts of genetic variation across species, and especially within a given ECM species, on function and interactions with the environment is not well understood. Here, we explore the extent of intra‐species variation between four isolates of the ECM fungus Pisolithus microcarpus , in terms of gene regulation, carbon metabolism and growth, and interactions with a host, Eucalyptus grandis . We demonstrate that, while a core response to the host is maintained by all of the isolates tested, they have distinct patterns of gene expression and carbon metabolism, resulting in the differential expression of isolate‐specific response pathways in the host plant. Together, these results highlight the importance of using a wider range of individuals within a species to understand the broader ecological roles of ECM fungi and their host interactions.

Plett, Krista L.↗

DOE OSTI2021

The Transcriptional Response of Soil Bacteria to Long-Term Warming and Short-Term Seasonal Fluctuations in a Terrestrial Forest

Terrestrial ecosystems are an important carbon store, and this carbon is vulnerable to microbial degradation with climate warming. After 30 years of experimental warming, carbon stocks in a temperate mixed deciduous forest were observed to be reduced by 30% in the heated plots relative to the controls. In addition, soil respiration was seasonal, as was the warming treatment effect. We therefore hypothesized that long-term warming will have higher expressions of genes related to carbohydrate and lipid metabolism due to increased utilization of recalcitrant carbon pools compared to controls. Because of the seasonal effect of soil respiration and the warming treatment, we further hypothesized that these patterns will be seasonal. We used RNA sequencing to show how the microbial community responds to long-term warming (~30 years) in Harvard Forest, MA. Total RNA was extracted from mineral and organic soil types from two treatment plots (+5°C heated and ambient control), at two time points (June and October) and sequenced using Illumina NextSeq technology. Treatment had a larger effect size on KEGG annotated transcripts than on CAZymes, while soil types more strongly affected CAZymes than KEGG annotated transcripts, though effect sizes overall were small. Although, warming showed a small effect on overall CAZymes expression, several carbohydrate-associated enzymes showed increased expression in heated soils (~68% of all differentially expressed transcripts). Further, exploratory analysis using an unconstrained method showed increased abundances of enzymes related to polysaccharide and lipid metabolism and decomposition in heated soils. Compared to long-term warming, we detected a relatively small effect of seasonal variation on community gene expression. Together, these results indicate that the higher carbohydrate degrading potential of bacteria in heated plots can possibly accelerate a self-reinforcing carbon cycle-temperature feedback in a warming climate.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Global analysis of SBP gene family in Brachypodium distachyon reveals its association with spike development

SQUAMOSA-promoter binding like proteins (SBPs/SPLs) are plant specific transcription factors targeted by miR156 and involved in various biological pathways, playing multi-faceted developmental roles. This gene family is not well characterized in Brachypodium. We identified a total of 18 SBP genes in B. distachyon genome. Phylogenetic analysis revealed that SBP gene family in Brachypodium expanded through large scale duplication. A total of 10 BdSBP genes were identified as targets of miR156. Transcript cleavage analysis of selected BdSBPs by miR156 confirmed their antagonistic connection. Alternative splicing was observed playing an important role in BdSBPs and miR156 interaction. Characterization of T-DNA Bdsbp9 mutant showed reduced plant growth and spike length, reflecting its involvement in the spike development. Expression of a majority of BdSBPs elevated during spikelet initiation. Specifically, BdSBP1 and BdSBP3 differentially expressed in response to vernalization. Differential transcript abundance of BdSBP1, BdSBP3, BdSBP8, BdSBP9, BdSBP14, BdSBP18 and BdSBP23 genes was observed during the spike development under high temperature. Co-expression network, protein–protein interaction and biological pathway analysis indicate that BdSBP genes mainly regulate transcription, hormone, RNA and transport pathways. Our work reveals the multi-layered control of SBP genes and demonstrates their association with spike development and temperature sensitivity in Brachypodium.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2024

Multilevel analysis between Physcomitrium patens and Mortierellaceae endophytes explores potential long‐standing interaction among land plants and fungi

SUMMARY The model moss species Physcomitrium patens has long been used for studying divergence of land plants spanning from bryophytes to angiosperms. In addition to its phylogenetic relationships, the limited number of differential tissues, and comparable morphology to the earliest embryophytes provide a system to represent basic plant architecture. Based on plant–fungal interactions today, it is hypothesized these kingdoms have a long‐standing relationship, predating plant terrestrialization. Mortierellaceae have origins diverging from other land fungi paralleling bryophyte divergence, are related to arbuscular mycorrhizal fungi but are free‐living, observed to interact with plants, and can be found in moss microbiomes globally. Due to their parallel origins, we assess here how two Mortierellaceae species, Linnemannia elongata and Benniella erionia , interact with P. patens in coculture. We also assess how Mollicute ‐related or Burkholderia ‐related endobacterial symbionts (MRE or BRE) of these fungi impact plant response. Coculture interactions are investigated through high‐throughput phenomics, microscopy, RNA‐sequencing, differential expression profiling, gene ontology enrichment, and comparisons among 99 other P. patens transcriptomic studies. Here we present new high‐throughput approaches for measuring P. patens growth, identify novel expression of over 800 genes that are not expressed on traditional agar media, identify subtle interactions between P. patens and Mortierellaceae, and observe changes to plant–fungal interactions dependent on whether MRE or BRE are present. Our study provides insights into how plants and fungal partners may have interacted based on their communications observed today as well as identifying L. elongata and B. erionia as modern fungal endophytes with P. patens .

59 BASIC BIOLOGICAL SCIENCES↗

NASA NTRSDate not supplied

Optimizing a Small RNAseq Analysis Pipeline for NASA GeneLab Using Open-Source Tools and Libraries

Small RNA sequencing (small RNAseq) is a powerful tool for studying the regulation of gene expression in various organisms. Small RNAseq has been leveraged in space biology research to study how expression of small RNAs, e.g. micro RNAs (miRNAs), small interfering RNAs (siRNAs), and piwi-interacting RNAs (piRNAs), change upon exposure to the space environment. NASA GeneLab currently hosts small RNAseq raw data derived from space-relevant experiments on the Open Science Data Repository (OSDR). To maximize the accessibility of these data to the scientific community, in addition to hosting raw data, which is only interpretable by bioinformaticians, GeneLab plans to process all small RNAseq datasets and make those processed data available to the scientific community via the OSDR. In this study, we present the development of the GeneLab standardized pipeline for processing small RNAseq datasets. Using human, plant, and synthetic small RNAseq datasets, we interrogate various open-source software and publicly available databases to evaluate their accuracy and reproducibility in each step of the pipeline. For quality control and adapter detection and trimming, we evaluated TrimGalore!, FASTX, SeqKit, and DNApi methods to optimize alignment to reference genomes. We compared BWA, Bowtie, and Bowtie2 to determine the optimal alignment tool. For each alignment tool we also assessed various reference databases, including Ensembl reference genomes and different types of small RNA reference databases, including genome, hairpin, and miRNA references from the miRbase and MirGeneDB databases. To quantify the aligned data, we compared SAMtools, HTSeq, and RSEM for counting alignment events from each alignment tool used. Finally, we evaluated various tools, including DESeq2 and EdgeR, for data normalization and subsequent differential expression analysis. We will present the results from our comparative analyses for each pipeline step and propose a consensus pipeline for processing small RNAseq data derived from various organisms exposed to the space environment.

SmallRNAseq, NASA GeneLab, quality control, adapte↗

DOE OSTI2025

Genetic_and_epigenetic_signatures_of_Populus_Trichocarpa_in_response_to_abiotic_stress

Plants respond to abiotic stresses such as drought, heat, and salinity through both shared and stress-specific regulatory pathways. However, the role of epigenetic mechanisms, particularly DNA methylation, in modulating these responses is still underexplored. Here, we integrated transcriptome profiling with whole-genome bisulfite sequencing (WGBS) to investigate the dynamic relationship between gene expression and DNA methylation in Populus trichocarpa during brief early-stage exposure to abiotic stress. Each stress elicited distinct transcriptional and methylation signatures; however, a subset of conserved stress-responsive genes was commonly regulated across treatments. Approximately 5% of differentially expressed genes also displayed differential methylation patterns, suggesting a coordinated role for DNA methylation in regulating gene expression. Motif enrichment analysis of differentially methylated regions revealed binding sites of key transcription factor families, including ERF, bHLH, and ABF, highlighting potential role for methylation in modulating transcription factor targeting. Furthermore, we identified stress-inducible molecular markers with potential applications in early stress detection and functional dissection of gene regulatory pathways. Together, these findings provide new insights into the coordinated genetic and epigenetic responses to abiotic stress and provide a foundation for developing biosensors and breeding strategies to enhance stress resilience woody plants.

CBI↗

DOE OSTI2026

Genetic and epigenetic signatures of Populus trichocarpa in response to abiotic stresses

Plants respond to abiotic stresses such as drought, heat, and salinity through both shared and stress-specific regulatory pathways. However, the role of epigenetic mechanisms, particularly DNA methylation, in modulating these responses is still underexplored. Here, we integrated transcriptome profiling with whole-genome bisulfite sequencing (WGBS) to investigate the dynamic relationship between gene expression and DNA methylation in Populus trichocarpa during brief early-stage exposure to abiotic stress. Each stress elicited distinct transcriptional and methylation signatures; however, a subset of conserved stress-responsive genes was commonly regulated across treatments. We found that ∼5-6% of differentially expressed genes also displayed differential methylation patterns, suggesting a coordinated role for DNA methylation in regulating gene expression. Motif enrichment analysis of differentially methylated regions revealed binding sites of key transcription factor families, including ERF, bHLH, and ABF, highlighting potential role for methylation in modulating transcription factor targeting. Furthermore, we identified stress-inducible molecular markers with potential applications in early stress detection and functional dissection of gene regulatory pathways. Together, these findings provide new insights into the coordinated genetic and epigenetic responses to abiotic stress and provide a foundation for developing biosensors and breeding strategies to enhance stress resilience woody plants.

DNA methylation↗

NASA NTRS1995

Expression of S100 beta in sensory and secretory cells of the vertebrate inner ear

We evaluated anti-S100 beta expression in the chick (Gallus domesticus) inner ear and determined that: 1) the monomer anti-S100 beta is expressed differentially in the vestibular and auditory perikarya; 2) expression of S100 beta in the afferent nerve terminals is time-related to synapse and myelin formation; 3) the expression of the dimer anti-S100 alpha alpha beta beta and monomer anti-S100 beta overlaps in most inner ear cell types. Most S100 alpha alpha beta beta positive cells express S100 beta, but S100 beta positive cells do not always express S100 alpha alpha beta beta. 4) the expression of S100 beta is diffused over the perikaryal cytoplasm and nuclei of the acoustic ganglia but is concentrated over the nuclei of the vestibular perikarya. 6) S100 beta is expressed in secretory cells, and it is co-localized with GABA in sensory cells. 7) Color thresholding objective quantitation indicates that the amount of S100 beta was higher (mean 22, SD +/- 4) at E19 than at E9 (mean 34, SD +/- 3) in afferent axons. 8) Moreover, S100 beta was unchanged between E11-E19 in the perikaryal cytoplasm, but did change over the nuclei. At E9, 74%, and at E21, 5% of vestibular perikarya were positive. The data suggest that S100 beta may be physically associated with neuronal and ionic controlling cells of the vertebrate inner ear, where it could provide a dual ionic and neurotrophic modulatory function.

NASA Discipline Number 40-10↗

DOE OSTI2026

Longitudinal Plasma Proteomic Profiling Reveals Divergent Immune Responses in Durably Cured and Relapsed Pulmonary Tuberculosis

Background: Predicting the risk of tuberculosis (TB) relapse is vital to improving treatment outcomes. Although clinical risk factors of relapse are well characterized, the biological mechanisms driving relapse, particularly host immune responses, remain poorly understood. Elucidating these mechanisms is necessary to better predict relapse risk. Methods: We conducted a longitudinal, global proteomic study on 60 participants with active pulmonary TB, half who were durably cured and half who relapsed. Plasma was collected at seven time-points: at treatment initiation (baseline), during therapy, and 52 weeks post-baseline. Samples were analyzed by high-resolution LC-MS/MS. Results: 2,418 proteins were identified across all samples, with 1,756 being differentially expressed relative to baseline (unadjusted p < 0.05). 956 proteins were differentially abundant between cured and relapsed participants. Relapsed participants showed heightened humoral immunity throughout treatment, as well as upregulated complement activation and HDL particles. Cured participants exhibited elevated recovery-related pathways by week 4, including downregulated epithelial invasion and upregulated oxygen transport. Conclusions: Heightened humoral and innate immune responses were associated with relapse, whereas recovery signatures were associated with durable cure. These findings advance our understanding of host responses to treatment and provide a basis for developing blood-based biomarkers to identify patients at increased risk of relapse.

LC-MS/MS↗

DOE OSTI2023

The global dynamic of DNA methylation in response to heat stress revealed epigenetic mechanism of heat acclimation in Saccharina japonica

Saccharina japonica is an ecologically and economically important kelp in cold‐temperate regions. When it is cultivated on a large scale in the temperate and even subtropical zones, heat stress is a frequent abiotic stress. This study is the first attempt to reveal the regulatory mechanism of the response to heat stress from the perspective of DNA methylation in S. japonica . We firstly obtained the characteristics of variation in the methylome under heat stress, and observed that heat stress caused a slight increase in the overall methylation level and methylation rate, especially in the non‐coding regions of the genome. Secondly, we noted that methylation was probably one of factors affecting the expression of genes, and that methylation within the gene body was positively correlated with the gene expression (rho = 0.0784). Moreover, it was found that among the differentially expressed genes regulated by methylation, many genes were related to heat stress response, such as HSP gene family, genes of antioxidant enzymes, genes related to proteasome‐ubiquitination pathway, and plant cell signaling pathways. This study demonstrated that DNA methylation is involved in regulating the response to heat stress, laying a foundation for studying the acclimation and adaptation of S. japonica to heat stress from an epigenetic perspective.

Liu, Fuli↗

NASA NTRSDate not supplied

Translational Line of Sight: a multi-omics longitudinal study of the murine retinal response to spaceflight hazard analogs

The space environment includes unique hazards like radiation and microgravity which adversely affect physiology and behavior of humans and rodent models. To better characterize the retinal response to spaceflight, we assessed a multi-omics NASA GeneLab dataset where 6-month-old female mice were gamma irradiated and/or hindlimb unloaded for 21 days followed by whole transcriptome shotgun sequencing (RNA-Seq) and reduced representation bisulfite sequencing (RRBS) of retina samples collected at 7 days, 1 month or 4 months post-exposure. We compared time-matched epigenomic and transcriptomic retinal profiles revealing a total of 4,178 differentially methylated loci or regions, and 457 differentially expressed genes. Highest correlation in methylation differences was seen across different conditions at the same time point (e.g., between radiation exposure and hindlimb unloaded at 7 days). Biological processes related to nucleotide metabolism were enriched in all groups with activation at 1 month and suppression at 7 days and 4 months. Genes and processes related to Notch and Wnt signaling showed alterations 4 months post-exposure. Interestingly, Notch3 and Lrg1 showed differential patterns in the NASA Twins Study in-flight samples and in response to stressors in the murine retina in the current study. A total of 23 genes were both differentially methylated and expressed, including genes involved in retinal disease or cataract development (Crybb3, Fgfr1, Pitpnm3, Sipa1l3, Sox9) and inflammatory response (B4galt6, Ppm1a, Sphk1). To our knowledge, the current multi-omics analysis is the first multi-omics study to interrogate the epigenomic and transcriptomic impacts of radiation and hindlimb unloading on the retina in isolation and in combination. The results provide an insight into the retinal response to individual spaceflight hazard analogs and their interplay at different post-exposure stages and contributes towards a mechanistic understanding of spaceflight-induced vision impairment using ground-based models.

Prachi Kothiyal↗

NASA NTRSDate not supplied

A Multi-omics Longitudinal Study of the Murine Retinal Response to Chronic Low-dose Irradiation and/or Simulated Microgravity

Prachi Kothiyal↗

DOE OSTI2020

Evolutionary Variation in MADS Box Dimerization Affects Floral Development and Protein Abundance in Maize

Interactions between MADS-box transcription factors are critical in the regulation of floral development, and shifting MADS-box protein-protein interactions are predicted to have influenced floral evolution. However, precisely how evolutionary variation in protein-protein interactions affects MADS-box protein function remains unknown. To assess the impact of changing MADS-box protein-protein interactions on transcription factor function, we turned to the grasses, where interactions between B-class MADS-box proteins vary. We tested the functional consequences of this evolutionary variability using maize as an experimental system. We found that differential B-class dimerization was associated with subtle, quantitative differences in stamen shape. In contrast, differential dimerization resulted in large-scale changes to downstream gene expression. Differential dimerization also affected B-class complex composition and abundance, independent of transcript levels. This indicates that differential B-class dimerization affects protein degradation, revealing an important consequence for evolutionary variability in MADS-box interactions. Our results highlight complexity in the evolution of developmental gene networks - changing protein-protein interactions could affect not only the composition of transcription factor complexes, but also their degradation and persistence in developing flowers. Here, our results also show how coding change in a pleiotropic master regulator could have small, quantitative effects on development.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2021

PD-L1 induces macrophage polarization toward the M2 phenotype via Erk/Akt/mTOR

Highlights: • PD-L1 treatment upregulates CD206 but not CD86 expression. • PD-L1-induced DEGs are related to the immune response, mitochondrial function, and metabolism. • Nivolumab, LY294002, U0126, and rapamycin inhibit the PD-L1-stimulated CD206 expression. • PD-L1 promotes M2 polarization via Erk/Akt/mTOR. PD-L1 (programmed death-ligand 1) is the ligand of PD-1 (programmed cell death protein 1) and regulates inhibitory immune responses. It is well known that PD-L1 suppresses T cell function via binding to PD-1. However, little is known about the role of the PD-1/PD-L1 axis in macrophage polarization. According to previous studies, the function of the PD-1/PD-L1 axis in macrophage polarization is controversial, and the underlying mechanism has not been fully elucidated. Thus, we treated THP-1-derived macrophages with human PD-L1 Fc to determine the role of the PD-1/PD-L1 axis in macrophage polarization. To further explore the mechanism, we performed RNA sequencing and used specific inhibitors to identify the implicated signalling pathways. In this study, we found that PD-L1 induces the upregulation of CD206 expression, which is inhibited by nivolumab, LY294002, U0126, and rapamycin. Evaluation of differentially expressed genes (DEGs) and bioinformatics analysis indicated that PD-L1 also induces the upregulation of the expression of genes that maintain mitochondrial function and mediate metabolic switching. In addition, we did not detect PD-L1-induced CD86 alterations, indicating that PD-L1 treatment has no significant influence on M1 polarization. Taken together, these results suggest that PD-L1 binds to PD-1 and promotes M2 polarization accompanied by mitochondrial function enhancement and metabolic reprogramming via Erk/Akt/mTOR. This study elucidates the role of PD-L1 in macrophage polarization and verifies the underlying mechanisms for the first time. Considering that aberrantly upregulated PD-L1 expression contributes to a wide variety of diseases, targeting PD-L1-mediated macrophage polarization is a prospective therapeutic strategy for both neoplastic and nonneoplastic diseases.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2020

Single‐parent expression drives dynamic gene expression complementation in maize hybrids

SUMMARY Single‐parent expression (SPE) is defined as gene expression in only one of the two parents. SPE can arise from differential expression between parental alleles, termed non‐presence/absence (non‐PAV) SPE, or from the physical absence of a gene in one parent, termed PAV SPE. We used transcriptome data of diverse Zea mays (maize) inbreds and hybrids, including 401 samples from five different tissues, to test for differences between these types of SPE genes. Although commonly observed, SPE is highly genotype and tissue specific. A positive correlation was observed between the genetic distance of the two inbred parents and the number of SPE genes identified. Regulatory analysis showed that PAV SPE and non‐PAV SPE genes are mainly regulated by cis effects, with a small fraction under trans regulation. Polymorphic transposable element insertions in promoter sequences contributed to the high level of cis regulation for PAV SPE and non‐PAV SPE genes. PAV SPE genes were more frequently expressed in hybrids than non‐PAV SPE genes. The expression of parentally silent alleles in hybrids of non‐PAV SPE genes was relatively rare but occurred in most hybrids. Non‐PAV SPE genes with expression of the silent allele in hybrids are more likely to exhibit above high parent expression level than hybrids that do not express the silent allele, leading to non‐additive expression. This study provides a comprehensive understanding of the nature of non‐PAV SPE and PAV SPE genes and their roles in gene expression complementation in maize hybrids.

Li, Zhi↗

DOE OSTI2022

Transcriptome Mining Provides Insights into Cell Wall Metabolism and Fiber Lignification in Agave tequilana Weber

Resilience of growing in arid and semiarid regions and a high capacity of accumulating sugar-rich biomass with low lignin percentages have placed Agave species as an emerging bioenergy crop. Although transcriptome sequencing of fiber-producing agave species has been explored, molecular bases that control wall cell biogenesis and metabolism in agave species are still poorly understood. Here, through RNAseq data mining, we reconstructed the cellulose biosynthesis pathway and the phenylpropanoid route producing lignin monomers in A. tequilana, and evaluated their expression patterns in silico and experimentally. Most of the orthologs retrieved showed differential expression levels when they were analyzed in different tissues with contrasting cellulose and lignin accumulation. Phylogenetic and structural motif analyses of putative CESA and CAD proteins allowed to identify those potentially involved with secondary cell wall formation. RT-qPCR assays revealed enhanced expression levels of AtqCAD5 and AtqCESA7 in parenchyma cells associated with extraxylary fibers, suggesting a mechanism of formation of sclerenchyma fibers in Agave similar to that reported for xylem cells in model eudicots. Overall, our results provide a framework for understanding molecular bases underlying cell wall biogenesis in Agave species studying mechanisms involving in leaf fiber development in monocots.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2021

Age-related telomere attrition causes aberrant gene expression in sub-telomeric regions

Telomere attrition has been proposed as a biomarker and causal factor in aging. In addition to causing cellular senescence and apoptosis, telomere shortening has been found to affect gene expression in subtelomeric regions. Here, we analyzed the distribution of age-related differentially expressed genes from the GTEx RNA sequencing database of 54 tissue types from 979 human subjects and found significantly more upregulated than downregulated genes in subtelomeric regions as compared to the genome-wide average. Our data demonstrate spatial relationships between telomeres and gene expression in aging.

59 BASIC BIOLOGICAL SCIENCES↗