Identification of Wheat Glutamate Synthetase Gene Family and Expression Analysis under Nitrogen Stress.

Nitrogen (N), as the main component of biological macromolecules, maintains the basic process of plant growth and development. GOGAT, as a key enzyme in the N assimilation process, catalyzes α-ketoglutaric acid and glutamine to form glutamate. In this study, six GOGAT genes in wheat (Triticum aestivum L.) were identified and classified into two subfamilies, Fd-GOGAT (TaGOGAT2s) and NADH-GOGAT (TaGOGAT3s), according to the type of electron donor. Subcellular localization prediction showed that TaGOGAT3-D was localized in mitochondria and that the other five TaGOGATs were localized in chloroplasts. Via the analysis of promoter elements, many binding sites related to growth and development, hormone regulation and plant stress resistance regulations were found on the TaGOGAT promoters. The tissue-specificity expression analysis showed that TaGOGAT2s were mainly expressed in wheat leaves and flag leaves, while TaGOGAT3s were highly expressed in roots and leaves. The expression level of TaGOGATs and the enzyme activity of TaGOGAT3s in the leaves and roots of wheat seedlings were influenced by the treatment of N deficiency. This study conducted a systematic analysis of wheat GOGAT genes, providing a theoretical basis not only for the functional analysis of TaGOGATs, but also for the study of wheat nitrogen use efficiency (NUE).

GIGANTEA adjusts the response to shade at dusk by directly impinging on PHYTOCHROME INTERACTING FACTOR 7 function.

For plants adapted to bright light, a decrease in the amount of light received can be detrimental to their growth and survival. Consequently, in response to shade from surrounding vegetation, they initiate a suite of molecular and morphological changes known as the shade avoidance response through which stems and petioles elongate in search for light. Under sunlight-night cycles, the plant's responsiveness to shade varies across the day, being maximal at dusk time. While a role for the circadian clock in this regulation has long been proposed, mechanistic understanding of how it is achieved is incomplete. Here, we show that the clock component GIGANTEA (GI) directly interacts with the transcriptional regulator PHYTOCHROME INTERACTING FACTOR 7 (PIF7), a key player in the response to shade. GI represses PIF7 transcriptional activity and the expression of its target genes in response to shade, thereby fine-tuning the magnitude of the response to limiting light conditions. We find that under light/dark cycles, this function of GI is required to adequately modulate the gating of the response to shade at dusk. Importantly, we also show that this circuit primarily operates in epidermal cells, highlighting the relevance of tissue-specific clock-output connections for the regulation of plant development in resonance with the environment.

The cacao gene atlas: a transcriptome developmental atlas reveals highly tissue-specific and dynamically-regulated gene networks in Theobroma cacao L.

BACKGROUND: Theobroma cacao, the cocoa tree, is a tropical crop grown for its highly valuable cocoa solids and fat which are the basis of a 200-billion-dollar annual chocolate industry. However, the long generation time and difficulties associated with breeding a tropical tree crop have limited the progress of breeders to develop high-yielding disease-resistant varieties. Development of marker-assisted breeding methods for cacao requires discovery of genomic regions and specific alleles of genes encoding important traits of interest. To accelerate gene discovery, we developed a gene atlas composed of a large dataset of replicated transcriptomes with the long-term goal of progressing breeding towards developing high-yielding elite varieties of cacao. RESULTS: We describe the creation of the Cacao Transcriptome Atlas, its global characterization and define sets of genes co-regulated in highly organ- and temporally-specific manners. RNAs were extracted and transcriptomes sequenced from 123 different tissues and stages of development representing major organs and developmental stages of the cacao lifecycle. In addition, several experimental treatments and time courses were performed to measure gene expression in tissues responding to biotic and abiotic stressors. Samples were collected in replicates (3-5) to enable statistical analysis of gene expression levels for a total of 390 transcriptomes. To promote wide use of these data, all raw sequencing data, expression read mapping matrices, scripts, and other information used to create the resource are freely available online. We verified our atlas by analyzing the expression of genes with known functions and expression patterns in Arabidopsis (ACT7, LEA19, AGL16, TIP13, LHY, MYB2) and found their expression profiles to be generally similar between both species. We also successfully identified tissue-specific genes at two thresholds in many tissue types represented and a set of genes highly conserved across all tissues. CONCLUSION: The Cacao Gene Atlas consists of a gene expression browser with graphical user interface and open access to raw sequencing data files as well as the unnormalized and CPM normalized read count data mapped to several cacao genomes. The gene atlas is a publicly available resource to allow rapid mining of cacao gene expression profiles. We hope this resource will be used to help accelerate the discovery of important genes for key cacao traits such as disease resistance and contribute to the breeding of elite varieties to help farmers increase yields.

Navigating spatio-temporal microbiome dynamics: Environmental factors and trace elements shape the symbiont community of an invasive marine species.

The proliferation of marine invasive species is a mounting concern. While the role of microbial communities in invasive ascidian species is recognized, the role of seasonal shifts in microbiome composition remains largely unexplored. We sampled five individuals of the invasive ascidian Styela plicata quarterly from January 2020 to October 2021 in two harbours, examining gills, tunics, and surrounding water. By analysing Amplicon Sequence Variants (ASVs) and seawater trace elements, we found that compartment (seawater, tunic, or gills) was the primary differentiating factor, followed by harbour. Clear seasonal patterns were evident in seawater bacteria, less so in gills, and absent in tunics. We identified compartment-specific bacteria, as well as seasonal indicator ASVs and ASVs correlated with trace element concentrations. Among these bacteria, we found that Endozoicomonas, Hepatoplasma and Rhodobacteraceae species had reported functions which might be necessary for overcoming seasonality and trace element shifts. This study contributes to understanding microbiome dynamics in invasive holobiont systems, and the patterns found indicate a potential role in adaptation and invasiveness.

Understanding the multi-functionality and tissue-specificity of decellularized dental pulp matrix hydrogels for endodontic regeneration.

Dental pulp is the only soft tissue in the tooth which plays a crucial role in maintaining intrinsic multi-functional behaviors of the dentin-pulp complex. Nevertheless, the restoration of fully functional pulps after pulpitis or pulp necrosis, termed endodontic regeneration, remained a major challenge for decades. Therefore, a bioactive and in-situ injectable biomaterial is highly desired for tissue-engineered pulp regeneration. Herein, a decellularized matrix hydrogel derived from porcine dental pulps (pDDPM-G) was prepared and characterized through systematic comparison against the porcine decellularized nerve matrix hydrogel (pDNM-G). The pDDPM-G not only exhibited superior capabilities in facilitating multi-directional differentiation of dental pulp stem cells (DPSCs) during 3D culture, but also promoted regeneration of pulp-like tissues after DPSCs encapsulation and transplantation. Further comparative proteomic and transcriptome analyses revealed the differential compositions and potential mechanisms that endow the pDDPM-G with highly tissue-specific properties. Finally, it was realized that the abundant tenascin C (TNC) in pDDPM served as key factor responsible for the activation of Notch signaling cascades and promoted DPSCs odontoblastic differentiation. Overall, it is believed that pDDPM-G is a sort of multi-functional and tissue-specific hydrogel-based material that holds great promise in endodontic regeneration and clinical translation. STATEMENT OF SIGNIFICANCE: Functional hydrogel-based biomaterials are highly desirable for endodontic regeneration treatments. Decellularized extracellular matrix (dECM) preserves most extracellular matrix components of its native tissue, exhibiting unique advantages in promoting tissue regeneration and functional restoration. In this study, we prepared a porcine dental pulp-derived dECM hydrogel (pDDPM-G), which exhibited superior performance in promoting odontogenesis, angiogenesis, and neurogenesis of the regenerating pulp-like tissue, further showed its tissue-specificity compared to the peripheral nerve-derived dECM hydrogel. In-depth proteomic and transcriptomic analyses revealed that the activation of tenascin C-Notch axis played an important role in facilitating odontogenic regeneration. This biomaterial-based study validated the great potential of the dental pulp-specific pDDPM-G for clinical applications, and provides a springboard for research strategies in ECM-related regenerative medicine.

Molecular adaptations in response to exercise training are associated with tissue-specific transcriptomic and epigenomic signatures.

Regular exercise has many physical and brain health benefits, yet the molecular mechanisms mediating exercise effects across tissues remain poorly understood. Here we analyzed 400 high-quality DNA methylation, ATAC-seq, and RNA-seq datasets from eight tissues from control and endurance exercise-trained (EET) rats. Integration of baseline datasets mapped the gene location dependence of epigenetic control features and identified differing regulatory landscapes in each tissue. The transcriptional responses to 8 weeks of EET showed little overlap across tissues and predominantly comprised tissue-type enriched genes. We identified sex differences in the transcriptomic and epigenomic changes induced by EET. However, the sex-biased gene responses were linked to shared signaling pathways. We found that many G protein-coupled receptor-encoding genes are regulated by EET, suggesting a role for these receptors in mediating the molecular adaptations to training across tissues. Our findings provide new insights into the mechanisms underlying EET-induced health benefits across organs.

Tissue-specific differences in Ca2+ sensitivity of the mitochondrial permeability transition pore (PTP). Experiments in male rat liver and heart.

Permeability transition pore (PTP) opening dissipates ion and electron gradients across the internal mitochondrial membrane (IMM), including excess Ca2+ in the mitochondrial matrix. After opening, immediate PTP closure must follow to prevent outer membrane disruption, loss of cytochrome c, and eventual apoptosis. Flickering, defined as the rapid alternative opening/closing of PTP, has been reported in heart, which undergoes frequent, large variations in Ca2+. In contrast, in tissues that undergo depolarization events less often, such as the liver, PTP would not need to be as dynamic and thus these tissues would not be as resistant to stress. To evaluate this idea, it was decided to follow the reversibility of the permeability transition (PT) in isolated murine mitochondria from two different tissues: the very dynamic heart, and the liver, which suffers depolarizations less frequently. It was observed that in heart mitochondria PT remained reversible for longer periods and at higher Ca2+ loads than in liver mitochondria. In all cases, Ca2+ uptake was inhibited by ruthenium red and PT was delayed by Cyclosporine A. Characterization of this phenomenon included measuring the rate of oxygen consumption, organelle swelling and Ca2+ uptake and retention. Results strongly suggest that there are tissue-specific differences in PTP physiology, as it resists many more Ca2+ additions before opening in a highly active organ such as the heart than in an organ that seldom suffers Ca2+ loading, such as the liver.

Comparative transcriptomic analysis delineates adaptation strategies of Rana kukunoris toward cold stress on the Qinghai-Tibet Plateau.

BACKGROUND: Cold hardiness is fundamental for amphibians to survive during the extremely cold winter on the Qinghai-Tibet plateau. Exploring the gene regulation mechanism of freezing-tolerant Rana kukunoris could help us to understand how the frogs survive in winter. RESULTS: Transcriptome of liver and muscle of R. kukunoris collected in hibernation and spring were assisted by single molecule real-time (SMRT) sequencing technology. A total of 10,062 unigenes of R. kukunoris were obtained, and 9,924 coding sequences (CDS) were successfully annotated. Our examination of the mRNA response to whole body freezing and recover in the frogs revealed key genes concerning underlying antifreeze proteins and cryoprotectants (glucose and urea). Functional pathway analyses revealed differential regulated pathways of ribosome, energy supply, and protein metabolism which displayed a freeze-induced response and damage recover. Genes related to energy supply in the muscle of winter frogs were up-regulated compared with the muscle of spring frogs. The liver of hibernating frogs maintained modest levels of protein synthesis in the winter. In contrast, the liver underwent intensive high levels of protein synthesis and lipid catabolism to produce substantial quantity of fresh proteins and energy in spring. Differences between hibernation and spring were smaller than that between tissues, yet the physiological traits of hibernation were nevertheless passed down to active state in spring. CONCLUSIONS: Based on our comparative transcriptomic analyses, we revealed the likely adaptive mechanisms of R. kukunoris. Ultimately, our study expands genetic resources for the freezing-tolerant frogs.

CpG Islands, Gene Expression and Pseudogenization: A Case for a Potential Trilogy.

BACKGROUND: The promoters of mammalian genes contain clusters of CG dinucleotides known as CpG islands. Most mammalian housekeeping genes predominantly contain CpG islands (CGIs), facilitating gene transcription. Numerous studies have explored the physiological implications of the relationship between CGIs and gene expression. However, the evolutionary implications of this relationship remain largely unexplored. Pseudogenes, in contrast, are genomic remnants that have lost their function over evolutionary time. METHODS: In our current research, we employed comparative genomic techniques to demonstrate a correlation between the absence of gene expression due to a lack of CGIs in the gene promoters and pseudogenization. RESULTS: We showed that there is a significant enrichment of tissue-specific genes in the functional orthologs of pseudogenes. We also found a significant correlation between the lack of CGIs and enriched tissue specificity in these functional orthologs of pseudogenes. CONCLUSIONS: We inferred that perhaps tissue-specific genes are more prone to the process of pseudogenization. In this way, because of their impact on gene expression, CGIs may affect the fate of a gene. To our knowledge, this is the first study to propose a connection between CGIs, gene expression, and the pseudogenization process and discuss the evolutionary implications of this potential trilogy.

Enriched atlas of lncRNA and protein-coding genes for the GRCg7b chicken assembly and its functional annotation across 47 tissues.

Gene atlases for livestock are steadily improving thanks to new genome assemblies and new expression data improving the gene annotation. However, gene content varies across databases due to differences in RNA sequencing data and bioinformatics pipelines, especially for long non-coding RNAs (lncRNAs) which have higher tissue and developmental specificity and are harder to consistently identify compared to protein coding genes (PCGs). As done previously in 2020 for chicken assemblies galgal5 and GRCg6a, we provide a new gene atlas, lncRNA-enriched, for the latest GRCg7b chicken assembly, integrating "NCBI RefSeq", "EMBL-EBI Ensembl/GENCODE" reference annotations and other resources such as FAANG and NONCODE. As a result, the number of PCGs increases from 18,022 (RefSeq) and 17,007 (Ensembl) to 24,102, and that of lncRNAs from 5789 (RefSeq) and 11,944 (Ensembl) to 44,428. Using 1400 public RNA-seq transcriptome representing 47 tissues, we provided expression evidence for 35,257 (79%) lncRNAs and 22,468 (93%) PCGs, supporting the relevance of this atlas. Further characterization including tissue-specificity, sex-differential expression and gene configurations are provided. We also identified conserved miRNA-hosting genes with human counterparts, suggesting common function. The annotated atlas is available at gega.sigenae.org.

ApCPS2 contributes to medicinal diterpenoid biosynthesis and defense against insect herbivore in Andrographis paniculata.

Kalmegh (Andrographis paniculata) spatiotemporally produces medicinally-important ent-labdane-related diterpenoids (ent-LRDs); andrographolide (AD), 14-deoxy-11,12-didehydroandrographolide (DDAD), neoandrographolide (NAD). ApCPS1 and ApCPS2, the ent-copalyl pyrophosphate (ent-CPP)-producing class II diterpene synthases (diTPSs) were identified, but their contributions to ent-CPP precursor supply for ent-LRD biosynthesis were not well understood. Here, we characterized ApCPS4, an additional ent-CPP-forming diTPS. Further, we elucidated in planta function of the ent-CPP-producing diTPSs (ApCPS1,2,4) by integrating transcript-metabolite co-profiles, biochemical analysis and gene functional characterization. ApCPS1,2,4 localized to the plastids, where diterpenoid biosynthesis occurs in plants, but ApCPS1,2,4 transcript expression patterns and ent-LRD contents revealed a strong correlation of ApCPS2 expression and ent-LRD accumulation in kalmegh. ApCPS1,2,4 upstream sequences differentially activated ß-glucuronidase (GUS) in Arabidopsis and transiently-transformed kalmegh. Similar to higher expression of ApCPS1 in kalmegh stem, ApCPS1 upstream sequence activated GUS in stem/hypocotyl of Arabidopsis and kalmegh. However, ApCPS2,4 upstream sequences weakly activated GUS expression in Arabidopsis, which was not well correlated with ApCPS2,4 transcript expression in kalmegh tissues. Whereas, ApCPS2,4 upstream sequences could activate GUS expression at a considerable level in kalmegh leaf and roots/calyx, respectively, suggesting the involvement of transcriptional regulator(s) of ApCPS2,4 that might participate in kalmegh-specific diterpenoid pathway. Interestingly, ApCPS2-silenced kalmegh showed a drastic reduction in AD, DDAD and NAD contents and compromised defense against insect herbivore Spodoptera litura. However, ent-LRD contents and herbivore defense in ApCPS1 or ApCPS4-silenced plants remained largely unaltered. Overall, these results suggested an important role of ApCPS2 in producing ent-CPP for medicinal ent-LRD biosynthesis and defense against insect herbivore.

Genetic effects of sequence-conserved enhancer-like elements on human complex traits.

BACKGROUND: The vast majority of findings from human genome-wide association studies (GWAS) map to non-coding sequences, complicating their mechanistic interpretations and clinical translations. Non-coding sequences that are evolutionarily conserved and biochemically active could offer clues to the mechanisms underpinning GWAS discoveries. However, genetic effects of such sequences have not been systematically examined across a wide range of human tissues and traits, hampering progress to fully understand regulatory causes of human complex traits. RESULTS: Here we develop a simple yet effective strategy to identify functional elements exhibiting high levels of human-mouse sequence conservation and enhancer-like biochemical activity, which scales well to 313 epigenomic datasets across 106 human tissues and cell types. Combined with 468 GWAS of European (EUR) and East Asian (EAS) ancestries, these elements show tissue-specific enrichments of heritability and causal variants for many traits, which are significantly stronger than enrichments based on enhancers without sequence conservation. These elements also help prioritize candidate genes that are functionally relevant to body mass index (BMI) and schizophrenia but were not reported in previous GWAS with large sample sizes. CONCLUSIONS: Our findings provide a comprehensive assessment of how sequence-conserved enhancer-like elements affect complex traits in diverse tissues and demonstrate a generalizable strategy of integrating evolutionary and biochemical data to elucidate human disease genetics.

Epigenetic analyses in forensic medicine: future and challenges.

The possibility of using epigenetics in forensic investigation has gradually risen over the last few years. Epigenetic changes with their dynamic nature can either be inherited or accumulated throughout a lifetime and be reversible, prompting investigation of their use across various fields. In forensic sciences, multiple applications have been proposed, such as the discrimination of monozygotic twins, identifying the source of a biological trace left at a crime scene, age prediction, determination of body fluids and tissues, human behavior association, wound healing progression, and determination of the post-mortem interval (PMI). Despite all these applications, not all the studies considered the impact of PMI and post-sampling effects on the epigenetic modifications and the tissue-specificity of the epigenetic marks.This review aims to highlight the substantial forensic significance that epigenetics could support in various forensic investigations. First, basic concepts in epigenetics, describing the main epigenetic modifications and their functions, in particular, DNA methylation, histone modifications, and non-coding RNA, with a particular focus on forensic applications, were covered. For each epigenetic marker, post-mortem stability and tissue-specificity, factors that should be carefully considered in the study of epigenetic biomarkers in the forensic context, have been discussed. The advantages and limitations of using post-mortem tissues have been also addressed, proposing directions for these innovative strategies to analyze forensic specimens.

Comparative RNA-Seq Analysis Revealed Tissue-Specific Splicing Variations during the Generation of the PDX Model.

Tissue-specific gene expression generates fundamental differences in the function of each tissue and affects the characteristics of the tumors that are created as a result. However, it is unclear how much the tissue specificity is conserved during grafting of the primary tumor into an immune-compromised mouse model. Here, we performed a comparative RNA-seq analysis of four different primary-patient derived xenograft (PDX) tumors. The analysis revealed a conserved RNA biotype distribution of primary-PDX pairs, as revealed by previous works. Interestingly, we detected significant changes in the splicing pattern of PDX, which was mainly comprised of skipped exons. This was confirmed by splicing variant-specific RT-PCR analysis. On the other hand, the correlation analysis for the tissue-specific genes indicated overall strong positive correlations between the primary and PDX tumor pairs, with the exception of gastric cancer cases, which showed an inverse correlation. These data propose a tissue-specific change in splicing events during PDX formation as a variable factor that affects primary-PDX integrity.

Comparison of Alternative Splicing Landscapes Revealed by Long-Read Sequencing in Hepatocyte-Derived HepG2 and Huh7 Cultured Cells and Human Liver Tissue.

The long-read RNA sequencing developed by Oxford Nanopore Technologies provides a direct quantification of transcript isoforms, thereby making it possible to present alternative splicing (AS) profiles as arrays of single splice variants with different abundances. Additionally, AS profiles can be presented as arrays of genes characterized by the degree of alternative splicing (the DAS-the number of detected splice variants per gene). Here, we successfully utilized the DAS to reveal biological pathways influenced by the alterations in AS in human liver tissue and the hepatocyte-derived malignant cell lines HepG2 and Huh7, thus employing the mathematical algorithm of gene set enrichment analysis. Furthermore, analysis of the AS profiles as abundances of single splice variants by using the graded tissue specificity index τ provided the selection of the groups of genes expressing particular splice variants specifically in liver tissue, HepG2 cells, and Huh7 cells. The majority of these splice variants were translated into proteins products and appeal to be in focus regarding further insights into the mechanisms underlying cell malignization. The used metrics are intrinsically suitable for transcriptome-wide AS profiling using long-read sequencing.

Multitissue Circadian Proteome Atlas of WT and Per1-/-/Per2-/- Mice.

The molecular basis of circadian rhythm, driven by core clock genes such as Per1/2, has been investigated on the transcriptome level, but not comprehensively on the proteome level. Here we quantified over 11,000 proteins expressed in eight types of tissues over 46 h with an interval of 2 h, using WT and Per1/Per2 double knockout mouse models. The multitissue circadian proteome landscape of WT mice shows tissue-specific patterns and reflects circadian anticipatory phenomena, which are less obvious on the transcript level. In most peripheral tissues of double knockout mice, reduced protein cyclers are identified when compared with those in WT mice. In addition, PER1/2 contributes to controlling the anticipation of the circadian rhythm, modulating tissue-specific cyclers as well as key pathways including nucleotide excision repair. Severe intertissue temporal dissonance of circadian proteome has been observed in the absence of Per1 and Per2. The γ-aminobutyric acid might modulate some of these temporally correlated cyclers in WT mice. Our study deepens our understanding of rhythmic proteins across multiple tissues and provides valuable insights into chronochemotherapy. The data are accessible at https://prot-rhythm.prottalks.com/.

Discovery of Novel miRNAs in Colorectal Cancer: Potential Biological Roles and Clinical Utility.

Deregulated miRNAs are associated with colorectal cancer (CRC), with alterations depending on the tumor location. Novel tissue-specific miRNAs have been identified in different tumors and are associated with cancer. We used miRMaster to identify novel miRNAs in CRC from the TCGA and GEO data (discovery and validation groups). We used TCGA data from five tissues to analyze miRNA tissue specificity. miRDB was used to predict miRNA targets, and the UCSC Xena Browser was used to evaluate target expression. After successive analyses, we identified 15 novel miRNAs with the same expression patterns in CRC in both the discovery and validation groups. Four molecules (nov-miR-13844-5p, nov-miR-7154-5p, nov-miR-5035-3p, and nov-miR-590-5p) were differentially expressed in proximal and distal CRC. The nov-miR-3345-5p and nov-miR-13172-3p, which are upregulated in tumors, are only expressed in colorectal tissues. These molecules have been linked to a worse prognosis in right-sided colon and rectal carcinomas. An analysis revealed an association between eight novel miRNAs and 81 targets, mostly cancer-related genes, with varying expression based on tumor location. These findings provide new miRNAs with potential biological relevance, molecular biomarkers, and therapeutic targets for CRC treatment.

Time-varying and tissue-dependent effects of adiposity on leptin levels: A Mendelian randomization study.

Background: Findings from Mendelian randomization (MR) studies are conventionally interpreted as lifelong effects, which typically do not provide insight into the molecular mechanisms underlying the effect of an exposure on an outcome. In this study, we apply two recently developed MR approaches (known as 'lifecourse' and 'tissue-partitioned' MR) to investigate lifestage-specific effects and tissues of action in the relationship between adiposity and circulating leptin levels. Methods: Genetic instruments for childhood and adult adiposity were incorporated into a multivariable MR (MVMR) framework to estimate lifestage-specific effects on leptin levels measured during early life (mean age: 10 y) in the Avon Longitudinal Study of Parents and Children and in adulthood (mean age: 55 y) using summary-level data from the deCODE Health study. This was followed by partitioning body mass index (BMI) instruments into those whose effects are putatively mediated by gene expression in either subcutaneous adipose or brain tissues, followed by using MVMR to simultaneously estimate their separate effects on childhood and adult leptin levels. Results: There was strong evidence that childhood adiposity has a direct effect on leptin levels at age 10 y in the lifecourse (ß = 1.10 SD change in leptin levels, 95% CI = 0.90-1.30, p=6 × 10-28), whereas evidence of an indirect effect was found on adulthood leptin along the causal pathway involving adulthood body size (ß = 0.74, 95% CI = 0.62-0.86, p=1 × 10-33). Tissue-partitioned MR analyses provided evidence to suggest that BMI exerts its effect on leptin levels during both childhood and adulthood via brain tissue-mediated pathways (ß = 0.79, 95% CI = 0.22-1.36, p=6 × 10-3 and ß = 0.51, 95% CI = 0.32-0.69, p=1 × 10-7, respectively). Conclusions: Our findings demonstrate the use of lifecourse MR to disentangle direct and indirect effects of early-life exposures on time-varying complex outcomes. Furthermore, by integrating tissue-specific data, we highlight the etiological importance of appetite regulation in the effect of adiposity on leptin levels. Funding: This work was supported by the Integrative Epidemiology Unit, which receives funding from the UK Medical Research Council and the University of Bristol (MC_UU_00011/1).

Sexual conflict, heterochrony and tissue specificity as evolutionary problems of adaptive plasticity in development.

Differential gene expression represents a fundamental cause and manifestation of phenotypic plasticity. Adaptive phenotypic plasticity in gene expression as a trait evolves when alleles that mediate gene regulation serve to increase organismal fitness by improving the alignment of variation in gene expression with variation in circumstances. Among the diverse circumstances that a gene encounters are distinct cell types, developmental stages and sexes, as well as an organism's extrinsic ecological environments. Consequently, adaptive phenotypic plasticity provides a common framework to consider diverse evolutionary problems by considering the shared implications of alleles that produce context-dependent gene expression. From this perspective, adaptive plasticity represents an evolutionary resolution to conflicts of interest that arise from any negatively pleiotropic effects of expression of a gene across ontogeny, among tissues, between the sexes, or across extrinsic environments. This view highlights shared properties within the general relation of fitness, trait expression and context that may nonetheless differ substantively in the grain of selection within and among generations to influence the likelihood of adaptive plasticity as an evolutionary response. Research programmes that historically have focused on these separate issues may use the insights from one another by recognizing their shared dependence on context-dependent gene regulatory evolution.