Tissue-specific Gene Expression Research Articles

Genomic structural variation contributes to evolved changes in gene expression in high-altitude Tibetan sheep

Tibetan sheep were introduced to the Qinghai Tibet plateau roughly 3,000 B.P., making this species a good model for investigating genetic mechanisms of high-altitude adaptation over a relatively short timescale. Here, we characterize genomic structural variants (SVs) that distinguish Tibetan sheep from closely related, low-altitude Hu sheep, and we examine associated changes in tissue-specific gene expression. We document differentiation between the two sheep breeds in frequencies of SVs associated with genes involved in cardiac function and circulation. In Tibetan sheep, we identified high-frequency SVs in a total of 462 genes, including EPAS1, PAPSS2, and PTPRD. Single-cell RNA-Seq data and luciferase reporter assays revealed that the SVs had cis-acting effects on the expression levels of these three genes in specific tissues and cell types. In Tibetan sheep, we identified a high-frequency chromosomal inversion that exhibited modified chromatin architectures relative to the noninverted allele that predominates in Hu sheep. The inversion harbors several genes with altered expression patterns related to heart protection, brown adipocyte proliferation, angiogenesis, and DNA repair. These findings indicate that SVs represent an important source of genetic variation in gene expression and may have contributed to high-altitude adaptation in Tibetan sheep.

Cortisol, Stress, and Disease-Bidirectional Associations; Role for Corticosteroid-Binding Globulin?

Selye described stress as a unified neurohormonal mechanism maintaining homeostasis. Acute stress system activation is adaptive through neurocognitive, catecholaminergic, and immunomodulation mechanisms, followed by a reset via cortisol. Stress system components, the sympathoadrenomedullary system, hypothalamic-pituitary-adrenal axis, and limbic structures are implicated in many chronic diseases by establishing an altered homeostatic state, allostasis. Consequent "primary stress system disorders" were popularly accepted, with phenotypes based on conditions such as Cushing syndrome, pheochromocytoma, and adrenal insufficiency. Cardiometabolic and major depressive disorders are candidates for hypercortisolemic etiology, contrasting the "hypocortisolemic symptom triad" of stress sensitivity, chronic fatigue, and pain. However, acceptance of chronic stress etiology requires cause-and-effect associations, and practical utility such as therapeutics altering stress system function. Inherent predispositions to stress system perturbations may be relevant. Glucocorticoid receptor (GR) variants have been associated with metabolic/neuropsychological states. The SERPINA6 gene encoding corticosteroid-binding globulin (CBG), was the sole genetic factor in a single-nucleotide variation-genome-wide association study linkage study of morning plasma cortisol, a risk factor for cardiovascular disease, with alterations in tissue-specific GR-related gene expression. Studies showed genetically predicted high cortisol concentrations are associated with hypertension and anxiety, and low CBG concentrations/binding affinity, with the hypocortisolemic triad. Acquired CBG deficiency in septic shock results in 3-fold higher mortality when hydrocortisone administration produces equivocal results, consistent with CBG's role in spatiotemporal cortisol delivery. We propose some stress system disorders result from constitutional stress system variants rather than stressors themselves. Altered CBG:cortisol buffering may influence interstitial cortisol ultradian surges leading to pathological tissue effects, an example of stress system variants contributing to stress-related disorders.

TALE gene family: identification, evolutionary and expression analysis under various exogenous hormones and waterlogging stress in Cucumis sativus L.

BackgroundThree Amino acid Loop Extension (TALE) belongs to the homeobox group of genes that are important constituents of plant systems. The TALE gene family is instrumental not only in growth and development but also plays an essential role in regulating plant response to environmental adversaries.ResultsIn the present study, we isolated 21 CsTALE genes from the cucumber (Cucumis sativus L.) genome database. Bioinformatics tools were put in place to understand the structural and functional components of the CsTALE gene family. The evolutionary analysis dissected them into seven subclades (KNOX-I, KNOX-II, and BELL-I to BELL-V). The cis-acting elements in the promoter region of CsTALE genes disclosed that they are key regulators of hormonal and stress-related processes. Additionally, the STRING database advocated the concerting role of CsTALE proteins with other key transcription factors potent in plant developmental biology. The CsmiR319 and CsmiR167a-3p targeting the CsTALE15 and CsTALE16, respectively, further assert the importance of the CsTALE gene family posttranscriptional-related processes. Tissue-specific gene expression unfolded the fundamental involvement of CsTALE genes as they were expressed throughout the developmental stages. Under waterlogging stress, the CsTALE17 expressed significantly higher values in WL, WL-NAA, and WL-ETH but not in WL-MeJA-treated samples.ConclusionsThe present study reveals the evolution and functions of the CsTALE gene family in cucumber. Our work will provide a platform that will help future researchers address the issue of waterlogging stress in the Yangtze River Delta.

Optimal variable identification for accurate detection of causal expression Quantitative Trait Loci with applications in heart-related diseases

Gene expression plays a pivotal role in various diseases, contributing significantly to their mechanisms. Most GWAS risk loci are in non-coding regions, potentially affecting disease risk by altering gene expression in specific tissues. This expression is notably tissue-specific, with genetic variants substantially influencing it. However, accurately detecting the expression Quantitative Trait Loci (eQTL) is challenging due to limited heritability in gene expression, extensive linkage disequilibrium (LD), and multiple causal variants. The single variant association approach in eQTL analysis is limited by its susceptibility to capture the combined effects of multiple variants, and a bias towards common variants, underscoring the need for a more robust method to accurately identify causal eQTL variants. To address this, we developed an algorithm, CausalEQTL, which integrates L0 +L1 penalized regression with an ensemble approach to localize eQTL, thereby enhancing prediction performance precisely. Our results demonstrate that CausalEQTL outperforms traditional models, including LASSO, Elastic Net, Ridge, in terms of power and overall performance. Furthermore, analysis of heart tissue data from the GTEx project revealed that eQTL sites identified by our algorithm provide deeper insights into heart-related tissue eQTL detection. This advancement in eQTL mapping promises to improve our understanding of the genetic basis of tissue-specific gene expression and its implications in disease. The source code and identified causal eQTLs for CausalEQTL are available on GitHub: https://github.com/zhc-moushang/CausalEQTL.

Primary Human Cell-Derived Extracellular Matrix from Decellularized Fibroblast Microtissues with Tissue-Dependent Composition and Microstructure.

Human extracellular matrix (ECM) exhibits complex protein composition and architecture depending on tissue and disease state, which remains challenging to reverse engineer. One promising approach is based on cell-secreted ECM from primary human fibroblasts that can be decellularized into acellular biomaterials. However, fibroblasts cultured on rigid culture plastic or biomaterial scaffolds can experience aberrant mechanical cues that perturb the biochemical, mechanical, and the efficiency of ECM production. Here, we demonstrate a method for preparing decellularized ECM using primary human fibroblasts with tissue and disease-specific features with two case studies: (1) cardiac fibroblasts; (2) lung fibroblasts from healthy or diseased donors. Cells aggregate into engineered microtissues in low adhesion microwells that deposited ECM and can be decellularized. We systematically investigate microtissue morphology, matrix architecture, and mechanical properties, along with transcriptomic and proteomic analysis. Microtissues exhibited tissue-specific gene expression and proteomics profiling, with ECM complexity similar to native tissues. Healthy lung microtissues exhibited web-like fibrillar collagen compared to dense patches in healthy heart microtissues. Diseased lung exhibited more disrupted collagen architecture than healthy. Decellularized microtissues had tissue-specific mechanical stiffness that was physiologically relevant. Importantly, decellularized microtissues supported viability and proliferation of human cells. We show that engineered microtissues of primary human fibroblasts seeded in low-adhesion microwells can be decellularized to produce human, tissue and disease-specific ECM. This approach should be widely applicable for generating personalized matrix that recapitulate tissues and disease states, relevant for culturing patient cells ex vivo as well as implantation for therapeutic treatments. The online version contains supplementary material available at 10.1007/s12195-024-00809-y.

Winners vs. losers: Schistosoma mansoni intestinal and liver eggs exhibit striking differences in gene expression and immunogenicity.

The eggs of the blood fluke Schistosoma mansoni are the main cause of the clinical manifestations of chronic schistosomiasis. After laying, the egg "winners" attach to the endothelium of the mesenteric vein and, after a period of development, induce the growth of a small granuloma, which facilitates their passage to the intestinal lumen. Egg "losers" carried by the bloodstream to non-specific tissues also undergo full development and induce large granuloma formation, but their life ends there. Although these trapped eggs represent a dead end in the parasite life cycle, the vast majority of studies attempting to describe the biology of the S. mansoni eggs have studied these liver-trapped "losers" instead of migrating intestinal "winners". This raises the fundamental question of how these eggs differ. With robust comparative transcriptomic analysis performed on S. mansoni eggs isolated 7 weeks post infection, we show that gene expression is critically dependent on tissue localization, both in the early and late stages of development. While mitochondrial genes and venom allergen-like proteins are significantly upregulated in mature intestinal eggs, well-described egg immunomodulators IPSE/alpha-1 and omega-1, together with micro-exon genes, are predominantly expressed in liver eggs. In addition, several proteases and protease inhibitors previously implicated in egg-host interactions display clear tissue-specific gene expression patterns. These major differences in gene expression could be then reflected in the observed different ability of liver and intestinal soluble egg antigens to elicit host immune responses and in the shorter viability of miracidia hatched from liver eggs. Our comparative analysis provides a new perspective on the biology of parasite's eggs in the context of their development and tissue localization. These findings could contribute to a broader and more accurate understanding of parasite eggs interactions with the host, which have historically been often restricted to liver eggs and sometimes inaccurately generalized.

The circadian clock in the choroid plexus drives rhythms in multiple cellular processes under the control of the suprachiasmatic nucleus

Choroid plexus (ChP), the brain structure primarily responsible for cerebrospinal fluid production, contains a robust circadian clock, whose role remains to be elucidated. The aim of our study was to [1] identify rhythmically controlled cellular processes in the mouse ChP and [2] assess the role and nature of signals derived from the master clock in the suprachiasmatic nuclei (SCN) that control ChP rhythms. To accomplish this goal, we used various mouse models (WT, mPer2Luc, ChP-specific Bmal1 knockout) and combined multiple experimental approaches, including surgical lesion of the SCN (SCNx), time-resolved transcriptomics, and single cell luminescence microscopy. In ChP of control (Ctrl) mice collected every 4 h over 2 circadian cycles in darkness, we found that the ChP clock regulates many processes, including the cerebrospinal fluid circadian secretome, precisely times endoplasmic reticulum stress response, and controls genes involved in neurodegenerative diseases (Alzheimer’s disease, Huntington’s disease, and frontotemporal dementia). In ChP of SCNx mice, the rhythmicity detected in vivo and ex vivo was severely dampened to a comparable extent as in mice with ChP-specific Bmal1 knockout, and the dampened cellular rhythms were restored by daily injections of dexamethasone in mice. Our data demonstrate that the ChP clock controls tissue-specific gene expression and is strongly dependent on the presence of a functional connection with the SCN. The results may contribute to the search for a novel link between ChP clock disruption and impaired brain health.

Genome-wide identification of the plant homeodomain-finger family in rye and ScPHD5 functions in cold tolerance and flowering time.

111 PHD genes were newly identified in rye genome and ScPHD5's role in regulating cold tolerance and flowering time was suggested. Plant homeodomain (PHD)-finger proteins regulate the physical properties of chromatin and control plant development and stress tolerance. Although rye (Secale cereale L.) is a major winter crop, PHD-finger proteins in rye have not been studied. Here, we identified 111 PHD genes in the rye genome that exhibited diverse gene and protein sequence structures. Phylogenetic tree analysis revealed that PHDs were genetically close in monocots and diverged from those in dicots. Duplication and synteny analyses demonstrated that ScPHDs have undergone several duplications during evolution and that high synteny is conserved among the Triticeae species. Tissue-specific and abiotic stress-responsive gene expression analyses indicated that ScPHDs were highly expressed in spikelets and developing seeds and were responsive to cold and drought stress. One of these genes, ScPHD5, was selected for further functional characterization. ScPHD5 was highly expressed in the spike tissues and was localized in the nuclei of rye protoplasts and tobacco leaves. ScPHD5-overexpressing Brachypodium was more tolerant to freezing stress than wild-type (WT), with increased CBF and COR gene expression. Additionally, these transgenic plants displayed an extremely early flowering phenotype that flowered more than two weeks earlier than the WT, and vernalization genes, rather than photoperiod genes, were increased in the WT. RNA-seq analysis revealed that diverse stress response genes, including HSPs, HSFs, LEAs, and MADS-box genes, were also upregulated in transgenic plants. Our study will help elucidate the roles of PHD genes in plant development and abiotic stress tolerance in rye.

Genome-Wide Epistatic Network Analyses of Semantic Fluency in Older Adults.

Semantic fluency impairment has been attributed to a wide range of neurocognitive and psychiatric conditions, especially in the older population. Moderate heritability estimates on semantic fluency were obtained from both twin and family-based studies suggesting genetic contributions to the observed variation across individuals. Currently, effort in identifying the genetic variants underlying the heritability estimates for this complex trait remains scarce. Using the semantic fluency scale and genome-wide SNP genotype data from the Long Life Family Study (LLFS), we performed a genome-wide association study (GWAS) and epistasis network analysis on semantic fluency in 2289 individuals aged over 60 years from the American LLFS cohorts and replicated the findings in 1129 individuals aged over 50 years from the Danish LLFS cohort. In the GWAS, two SNPs with genome-wide significance (rs3749683, p = 2.52 × 10-8; rs880179, p = 4.83 × 10-8) mapped to the CMYAS gene on chromosome 5 were detected. The epistasis network analysis identified five modules as significant (4.16 × 10-5 < p < 7.35 × 10-3), of which two were replicated (p < 3.10 × 10-3). These two modules revealed significant enrichment of tissue-specific gene expression in brain tissues and high enrichment of GWAS catalog traits, e.g., obesity-related traits, blood pressure, chronotype, sleep duration, and brain structure, that have been reported to associate with verbal performance in epidemiological studies. Our results suggest high tissue specificity of genetic regulation of gene expression in brain tissues with epistatic SNP networks functioning jointly in modifying individual verbal ability and cognitive performance.

Production of cleavage-resistant phytase transgenic pigs by handmade cloning

With the rapid development of intensive animal husbandry in the livestock industry, large quantities of manure waste containing phytate phosphorus are being generated. Phytase can effectively solve the problem of high phosphorus pollution in the feces of monogastric animals. Enviropig, which produces phytase in the salivary glands and secretes the enzyme in the saliva, were first generated in 1999. However, phytase is easily inactivated during digestion. To address this problem, cleavage-resistant phytase transgenic pigs were generated using handmade cloning in this study. Transgene construction was improved and three cell lines carrying Cafp were obtained. In total, 810 blastocysts were generated and 712 good-quality were transferred into six recipients. Fourteen piglets were born, of which six survived after weaning. Polymerase chain reaction and sequencing results showed that seven (three live and four dead) of the fourteen piglets carried Cafp. Phytase activity in the saliva of the six live cloned pigs was tested at four months of age, and only one pig had 0.155 FTU/mL enzyme activity. The other five pigs may not have been activated in the transgenic parotid gland. Among all the transgenic pigs, the highest phosphorus digestion rate was 59.2% of intake, representing a 25.4% decrease in fecal emission compared to the average of controls. Immunohistochemical results on the three Cafp-positive pigs that died after six months of age showed that the transgene was only expressed in parotid glands, confirming tissue-specific gene expression. In conclusion, cleavage-resistant phytase transgenic pigs were successfully produced through handmade cloning. The cloned pigs offer a unique biological approach to managing phosphorus nutrition and environmental pollution in animal husbandry.

UDP-glycosyltransferase gene family expansion and functional perspectives in five tenebrionid beetles

UDP-glycosyltransferases (UGTs) are widely distributed biotransformation enzymes crucial for various physiological processes. Through genomic analysis, 216 UGT genes were identified in five tenebrionoid beetles: Tebebrio molitor (46), Tribolium castaneum (43), Asbolus verrucosus (30), Hycleus cichorii (60), and Hycleus phaleratus (37), classified into 14 families, with 13 being coleopteran-specific. Lineage-specific expansion of the UGT321, UGT323, UGT324, and UGT347 families was observed, with evidence of gene duplications within these families. Notably, phytophagous blister beetles showed higher UGT321 gene counts compared to omnivorous darkling beetles, while UGT324 family expansion was noted in the latter. Investigation of temporal and tissue-specific UGT gene expression in T. molitor revealed over 60% of TmUGTs were highly expressed in the salivary gland, fat body, midgut, and Malpighian tubule, indicating roles in metabolism and detoxification. Particularly, UGT321 family members exhibited high expression across various tissues, suggesting diverse functions, including potential involvement in olfaction and detoxification. These findings emphasize the multifunctionality of UGTs and provide a basis for further research, utilizing T. molitor as a model organism.

Genome-wide identification and analyses of ZmAPY genes reveal their roles involved in maize development and abiotic stress responses

Apyrase is a class of enzyme that catalyzes the hydrolysis of nucleoside triphosphates/diphosphates (NTP/NDP), which widely involved in regulation of plant growth and stress responses. However, apyrase family genes in maize have not been identified, and their characteristics and functions are largely unknown. In this study, we identified 16 apyrases (named as ZmAPY1-ZmAPY16) in maize genome, and analyzed their phylogenetic relationships, gene structures, chromosomal distribution, upstream regulatory transcription factors and expression patterns. Analysis of the transcriptome database unveiled tissue-specific and abiotic stress-responsive expression of ZmAPY genes in maize. qPCR analysis further confirmed their responsiveness to drought, heat, and cold stresses. Association analyses indicated that variations of ZmAPY5 and ZmAPY16 may regulate maize agronomic traits and drought responses. Our findings shed light on the molecular characteristics and evolutionary history of maize apyrase genes, highlighting their roles in various biological processes and stress responses. This study forms a basis for further exploration of apyrase functions in maize.

The impact of exercise on gene regulation in association with complex trait genetics.

Endurance exercise training is known to reduce risk for a range of complex diseases. However, the molecular basis of this effect has been challenging to study and largely restricted to analyses of either few or easily biopsied tissues. Extensive transcriptome data collected across 15 tissues during exercise training in rats as part of the Molecular Transducers of Physical Activity Consortium has provided a unique opportunity to clarify how exercise can affect tissue-specific gene expression and further suggest how exercise adaptation may impact complex disease-associated genes. To build this map, we integrate this multi-tissue atlas of gene expression changes with gene-disease targets, genetic regulation of expression, and trait relationship data in humans. Consensus from multiple approaches prioritizes specific tissues and genes where endurance exercise impacts disease-relevant gene expression. Specifically, we identify a total of 5523 trait-tissue-gene triplets to serve as a valuable starting point for future investigations [Exercise; Transcription; Human Phenotypic Variation].

Seasonal tissue-specific gene expression in wild crown-of-thorns starfish reveals reproductive and stress-related transcriptional systems.

Animals are influenced by the season, yet we know little about the changes that occur in most species throughout the year. This is particularly true in tropical marine animals that experience relatively small annual temperature and daylight changes. Like many coral reef inhabitants, the crown-of-thorns starfish (COTS), well known as a notorious consumer of corals and destroyer of coral reefs, reproduces exclusively in the summer. By comparing gene expression in 7 somatic tissues procured from wild COTS sampled on the Great Barrier Reef, we identified more than 2,000 protein-coding genes that change significantly between summer and winter. COTS genes that appear to mediate conspecific communication, including both signalling factors released into the surrounding sea water and cell surface receptors, are up-regulated in external secretory and sensory tissues in the summer, often in a sex-specific manner. Sexually dimorphic gene expression appears to be underpinned by sex- and season-specific transcription factors (TFs) and gene regulatory programs. There are over 100 TFs that are seasonally expressed, 87% of which are significantly up-regulated in the summer. Six nuclear receptors are up-regulated in all tissues in the summer, suggesting that systemic seasonal changes are hormonally controlled, as in vertebrates. Unexpectedly, there is a suite of stress-related chaperone proteins and TFs, including HIFa, ATF3, C/EBP, CREB, and NF-κB, that are uniquely and widely co-expressed in gravid females. The up-regulation of these stress proteins in the summer suggests the demands of oogenesis in this highly fecund starfish affects protein stability and turnover in somatic cells. Together, these circannual changes in gene expression provide novel insights into seasonal changes in this coral reef pest and have the potential to identify vulnerabilities for targeted biocontrol.

Germline cis variant determines epigenetic regulation of the anti-cancer drug metabolism gene dihydropyrimidine dehydrogenase (DPYD)

Enhancers are critical for regulating tissue-specific gene expression, and genetic variants within enhancer regions have been suggested to contribute to various cancer-related processes, including therapeutic resistance. However, the precise mechanisms remain elusive. Using a well-defined drug-gene pair, we identified an enhancer region for dihydropyrimidine dehydrogenase (DPD, DPYD gene) expression that is relevant to the metabolism of the anti-cancer drug 5-fluorouracil (5-FU). Using reporter systems, CRISPR genome-edited cell models, and human liver specimens, we demonstrated in vitro and vivo that genotype status for the common germline variant (rs4294451; 27% global minor allele frequency) located within this novel enhancer controls DPYD transcription and alters resistance to 5-FU. The variant genotype increases recruitment of the transcription factor CEBPB to the enhancer and alters the level of direct interactions between the enhancer and DPYD promoter. Our data provide insight into the regulatory mechanisms controlling sensitivity and resistance to 5-FU.

Germline cis variant determines epigenetic regulation of the anti-cancer drug metabolism gene dihydropyrimidine dehydrogenase (DPYD).

Enhancers are critical for regulating tissue-specific gene expression, and genetic variants within enhancer regions have been suggested to contribute to various cancer-related processes, including therapeutic resistance. However, the precise mechanisms remain elusive. Using a well-defined drug-gene pair, we identified an enhancer region for dihydropyrimidine dehydrogenase (DPD, DPYD gene) expression that is relevant to the metabolism of the anti-cancer drug 5-fluorouracil (5-FU). Using reporter systems, CRISPR genome-edited cell models, and human liver specimens, we demonstrated in vitro and vivo that genotype status for the common germline variant (rs4294451; 27% global minor allele frequency) located within this novel enhancer controls DPYD transcription and alters resistance to 5-FU. The variant genotype increases recruitment of the transcription factor CEBPB to the enhancer and alters the level of direct interactions between the enhancer and DPYD promoter. Our data provide insight into the regulatory mechanisms controlling sensitivity and resistance to 5-FU.

Alternative splicing and its regulation in the malaria vector Anopheles gambiae

IntroductionAlternative splicing (AS) is a highly conserved mechanism that allows for the expansion of the coding capacity of the genome, through modifications of the way that multiple isoforms are expressed or used to generate different phenotypes. Despite its importance in physiology and disease, genome-wide studies of AS are lacking in most insects, including mosquitoes. Even for model organisms, chromatin associated processes involved in the regulation AS are poorly known.MethodsIn this study, we investigated AS in the mosquito Anopheles gambiae in the context of tissue-specific gene expression and mosquito responses to a Plasmodium falciparum infection, as well as the relationship between patterns of differential isoform expression and usage with chromatin accessibility changes. For this, we combined RNA-seq and ATAC-seq data from A. gambiae midguts and salivary glands, infected and non-infected.ResultsWe report differences between tissues in the expression of 392 isoforms and in the use of 247 isoforms. Secondly, we find a clear and significant association between chromatin accessibility states and tissue-specific patterns of AS. The analysis of differential accessible regions located at splicing sites led to the identification of several motifs resembling the binding sites of Drosophila transcription factors. Finally, the genome-wide analysis of tissue-dependent enhancer activity revealed that approximately 20% of A. gambiae transcriptional enhancers annotate to a differentially expressed or used isoform, and that their activation status is linked to AS differences between tissues.ConclusionThis research elucidates the role of AS in mosquito vector gene expression and identifies regulatory regions potentially involved in AS regulation, which could be important in the development of novel strategies for vector control.

Tissue-specific gene expression of genome-wide significant loci associated with major depressive disorder subtypes

Major depressive disorder (MDD) is a clinically and genetically heterogeneous disorder. To reduce heterogeneity, large-scale genome-wide association studies have recently identified genome-wide significant loci associated with seven MDD subtypes. However, it was unclear in which tissues the genes near those loci are specifically expressed. We investigated whether genes related to specific MDD subtypes would be preferably expressed in a specific tissue. At 14 novel subtype-specific loci related to seven MDD subtypes—(1) non-atypical-like features MDD, (2) early-onset MDD, (3) recurrent MDD, (4) MDD with suicidal thoughts, (5) MDD without suicidal thoughts, (6) MDD with moderate impairment, and (7) postpartum depression, we investigated whether 22 genome-wide significant genetic variant-mapped genes were tissue-specifically expressed in brain, female reproductive, male specific, cardiovascular, gastrointestinal, or urinary tissues in the Genotype-Tissue Expression (GTEx) subjects (n ≤ 948). To confirm the tissue-specific expression in the GTEx, we used independent Human Protein Atlas (HPA) RNA-seq subjects (n ≤ 95). Of 22 genes, nine and five genes were tissue-specifically expressed in brain and female reproductive tissues, respectively (p < 2.27 × 10−3). RTN1, ERBB4, and AMIGO1 related to early-onset MDD, recurrent MDD, or MDD with suicidal thoughts were highly expressed in brain tissues (d = 1.19–2.71), while OAS1, LRRC9, DHRS7, PSMA5, SYPL2, and GULP1 related to non-atypical-like features MDD, early-onset MDD, MDD with suicidal thoughts, or postpartum depression were expressed at low levels in brain tissues (d = −0.17−−1.48). DFNA5, CTBP2, PCNX4, SDCCAG8, and GULP1, which are related to early-onset MDD, MDD with moderate impairment, or postpartum depression, were highly expressed in female reproductive tissues (d = 0.80–2.08). Brain and female reproductive tissue-specific expression was confirmed in the HPA RNA-seq subjects. Our findings suggest that brain and female reproductive tissue-specific expression might contribute to the pathogenesis of MDD subtypes.

Emergence of enhancers at late DNA replicating regions

Enhancers are fast-evolving genomic sequences that control spatiotemporal gene expression patterns. By examining enhancer turnover across mammalian species and in multiple tissue types, we uncover a relationship between the emergence of enhancers and genome organization as a function of germline DNA replication time. While enhancers are most abundant in euchromatic regions, enhancers emerge almost twice as often in late compared to early germline replicating regions, independent of transposable elements. Using a deep learning sequence model, we demonstrate that new enhancers are enriched for mutations that alter transcription factor (TF) binding. Recently evolved enhancers appear to be mostly neutrally evolving and enriched in eQTLs. They also show more tissue specificity than conserved enhancers, and the TFs that bind to these elements, as inferred by binding sequences, also show increased tissue-specific gene expression. We find a similar relationship with DNA replication time in cancer, suggesting that these observations may be time-invariant principles of genome evolution. Our work underscores that genome organization has a profound impact in shaping mammalian gene regulation.

Analysis of long-range chromatin contacts, compartments and looping between mouse embryonic stem cells, lens epithelium and lens fibers

BackgroundNuclear organization of interphase chromosomes involves individual chromosome territories, “open” and “closed” chromatin compartments, topologically associated domains (TADs) and chromatin loops. The DNA- and RNA-binding transcription factor CTCF together with the cohesin complex serve as major organizers of chromatin architecture. Cellular differentiation is driven by temporally and spatially coordinated gene expression that requires chromatin changes of individual loci of various complexities. Lens differentiation represents an advantageous system to probe transcriptional mechanisms underlying tissue-specific gene expression including high transcriptional outputs of individual crystallin genes until the mature lens fiber cells degrade their nuclei.ResultsChromatin organization between mouse embryonic stem (ES) cells, newborn (P0.5) lens epithelium and fiber cells were analyzed using Hi-C. Localization of CTCF in both lens chromatins was determined by ChIP-seq and compared with ES cells. Quantitative analyses show major differences between number and size of TADs and chromatin loop size between these three cell types. In depth analyses show similarities between lens samples exemplified by overlaps between compartments A and B. Lens epithelium-specific CTCF peaks are found in mostly methylated genomic regions while lens fiber-specific and shared peaks occur mostly within unmethylated DNA regions. Major differences in TADs and loops are illustrated at the ~ 500 kb Pax6 locus, encoding the critical lens regulatory transcription factor and within a larger ~ 15 Mb WAGR locus, containing Pax6 and other loci linked to human congenital diseases. Lens and ES cell Hi-C data (TADs and loops) together with ATAC-seq, CTCF, H3K27ac, H3K27me3 and ENCODE cis-regulatory sites are shown in detail for the Pax6, Sox1 and Hif1a loci, multiple crystallin genes and other important loci required for lens morphogenesis. The majority of crystallin loci are marked by unexpectedly high CTCF-binding across their transcribed regions.ConclusionsOur study has generated the first data on 3-dimensional (3D) nuclear organization in lens epithelium and lens fibers and directly compared these data with ES cells. These findings generate novel insights into lens-specific transcriptional gene control, open new research avenues to study transcriptional condensates in lens fiber cells, and enable studies of non-coding genetic variants linked to cataract and other lens and ocular abnormalities.