Tissue-specific Gene Expression Research Articles

Ontogeny and tissue specific expression profiles of recombination activating genes (RAGs) during development in Nile tilapia, Oreochromisniloticus

Recombination activating genes (RAGs) mediates the process of rearrangement and somatic recombination (V(D)J) to generate different antibody repertoire. Studies on the expression pattern of adaptive immune genes during ontogenic development are crucial for the formulation of fish immunization strategy. In the present study, Nile tilapia was taken to explore the relative expression profile of RAG genes during their developmental stages. The developmental stages of Nile tilapia, i.e., unfertilized egg, 0, 2, 4, 6, 8, 10, 12, 14, 16, 18, 20, 22, 24, 26, 28 and 30 days post-hatch (dph) and kidney, blood, gill, liver and spleen tissues from adult fish were collected and the cDNA synthesis was carried out. Gene specific primers for RAG-1 and RAG-2 of Nile tilapia were designed and their annealing temperature (Tm) was optimized by gradient PCR. Consequently, PCR was performed to confirm the specific amplification of RAG-1 and RAG-2 genes. Quantitative real-time PCR (qRT-PCR) gene expression of RAG-1 and RAG-2 were noticed in all the developmental stages; however, a significant increase was observed after 12 dph and peaked at 24 dph, followed by a gradual decrease until 30 dph. Tissue-specific gene expression profiling revealed that the highest expression of RAG-1 and RAG-2 was observed in the kidney, followed by spleen, gill, liver and blood. The findings of the study explored the suitable timing of lymphoid maturation that could be technically used for the adoption of strategies to improve disease resistance of fish larvae for mitigating larval mortality.

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.

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.

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.

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.

Tissue specific tumor-gene link prediction through sampling based GNN using a heterogeneous network.

A tissue sample is a valuable resource for understanding a patient's symptoms and health status in relation to tumor growth. Recent research seeks to establish a connection between tissue-specific tumor samples and genetic markers (genes). This breakthrough has paved the way for personalized cancer therapies. With this motivation, the proposed model constructs a heterogeneous network based on tumor sample-gene relation data and gene-gene interaction data. This network also incorporates tissue-specific gene expression and primary site-based gene counts as features, enabling tissue-specific predictions. Graph neural networks (GNNs) have proven effective in modeling complex interactions and predicting links within this network. The proposed model has successfully predicted tumor-gene associations by leveraging sampling-based GNNs and link layer embedding. The model's performance metrics, such as AUC-ROC scores, reached approximately 94%, demonstrating the potential of this heterogeneous network in predicting tissue-specific tumor sample-gene links. This paper's findings highlight the importance of tissue-specific associations in cancer research.

Evolution of tissue-specific expression of ancestral genes across vertebrates and insects.

Regulation of gene expression is arguably the main mechanism underlying the phenotypic diversity of tissues within and between species. Here we assembled an extensive transcriptomic dataset covering 8 tissues across 20 bilaterian species and performed analyses using a symmetric phylogeny that allowed the combined and parallel investigation of gene expression evolution between vertebrates and insects. We specifically focused on widely conserved ancestral genes, identifying strong cores of pan-bilaterian tissue-specific genes and even larger groups that diverged to define vertebrate and insect tissues. Systematic inferences of tissue-specificity gains and losses show that nearly half of all ancestral genes have been recruited into tissue-specific transcriptomes. This occurred during both ancient and, especially, recent bilaterian evolution, with several gains being associated with the emergence of unique phenotypes (for example, novel cell types). Such pervasive evolution of tissue specificity was linked to gene duplication coupled with expression specialization of one of the copies, revealing an unappreciated prolonged effect of whole-genome duplications on recent vertebrate evolution.

The transcriptome assembly of the European freshwater mussel Unio elongatulus C. Pfeiffer, 1825

Freshwater mussels of the order Unionida are a global conservation concern. Species of this group are strictly freshwater, sessile, slow-growing animals and, extremely sensitive to environmental changes. Human-mediated changes in freshwater habitats are imposing enormous pressure on the survival of freshwater mussels. Although a few flagship species are protected in Europe, other highly imperilled species receive much less attention. Moreover, knowledge about biology, ecology, and evolution and proper conservation assessments of many European species are still sparse. This knowledge gap is further aggravated by the lack of genomic resources available, which are key tools for conservation. Here we present the transcriptome assembly of Unio elongatulus C. Pfeiffer, 1825, one of the least studied European freshwater mussels. Using the individual sequencing outputs from eight physiologically representative mussel tissues, we provide an annotated panel of tissue-specific Relative Gene Expression profiles. These resources are pivotal to studying the species’ biological and ecological features, as well as helping to understand its vulnerability to current and future threats.

Cis-regulatory interfaces reveal the molecular mechanisms underlying the notochord gene regulatory network of Ciona

Tissue-specific gene expression is fundamental in development and evolution, and is mediated by transcription factors and by the cis-regulatory regions (enhancers) that they control. Transcription factors and their respective tissue-specific enhancers are essential components of gene regulatory networks responsible for the development of tissues and organs. Although numerous transcription factors have been characterized from different organisms, the knowledge of the enhancers responsible for their tissue-specific expression remains fragmentary. Here we use Ciona to study the enhancers associated with ten transcription factors expressed in the notochord, an evolutionary hallmark of the chordate phylum. Our results illustrate how two evolutionarily conserved transcription factors, Brachyury and Foxa2, coordinate the deployment of other notochord transcription factors. The results of these detailed cis-regulatory analyses delineate a high-resolution view of the essential notochord gene regulatory network of Ciona, and provide a reference for studies of transcription factors, enhancers, and their roles in development, disease, and evolution.

Olfactory genes affect major depression in highly educated, emotionally stable, lean women: a bridge between animal models and precision medicine

Most current approaches to establish subgroups of depressed patients for precision medicine aim to rely on biomarkers that require highly specialized assessment. Our present aim was to stratify participants of the UK Biobank cohort based on three readily measurable common independent risk factors, and to investigate depression genomics in each group to discover common and separate biological etiology. Two-step cluster analysis was run separately in males (n = 149,879) and females (n = 174,572), with neuroticism (a tendency to experience negative emotions), body fat percentage, and years spent in education as input variables. Genome-wide association analyses were implemented within each of the resulting clusters, for the lifetime occurrence of either a depressive episode or recurrent depressive disorder as the outcome. Variant-based, gene-based, gene set-based, and tissue-specific gene expression test were applied. Phenotypically distinct clusters with high genetic intercorrelations in depression genomics were found. A two-cluster solution was the best model in each sex with some differences including the less important role of neuroticism in males. In females, in case of a protective pattern of low neuroticism, low body fat percentage, and high level of education, depression was associated with pathways related to olfactory function. While also in females but in a risk pattern of high neuroticism, high body fat percentage, and less years spent in education, depression showed association with complement system genes. Our results, on one hand, indicate that alteration of olfactory pathways, that can be paralleled to the well-known rodent depression models of olfactory bulbectomy, might be a novel target towards precision psychiatry in females with less other risk factors for depression. On the other hand, our results in multi-risk females may provide a special case of immunometabolic depression.