Evidence has shown that differential transcriptomic profiles among human populations from diverse ancestries, supporting the role of genetic architecture in regulating gene expression alongside environmental stimuli. Genetic variants that regulate gene expression, known as expression quantitative trait loci (eQTL), are primarily shaped by human migration history and evolutionary forces, likewise, regulation of gene expression in principle could have been influenced by these events. Therefore, a comprehensive understanding of how human evolution impacts eQTL offers important insights into how phenotypic diversity is shaped. Recent studies, however, suggest that eQTL is enriched in genes that are selectively constrained. Whether eQTL is minimally affected by selective pressures remains an open question and requires comprehensive investigations. In addition, such studies are primarily dominated by the major populations of European ancestry, leaving many marginalized populations underrepresented. These observations indicate there exists a fundamental knowledge gap in the role of genomics variation on phenotypic diversity, which potentially hinders precision medicine. This article aims to revisit the abundance of eQTL across diverse populations and provide an overview of their impact from the population and evolutionary genetics perspective, subsequently discuss their influence on phenomics, as well as challenges and opportunities in the applications to precision medicine.

During the formation of kidney stones, the interaction between crystals and tubular epithelial cells (TECs) leads to tubular injury and dysfunction, which in turn promote stone formation. However, the molecular mechanisms underlying these changes in TECs remain elusive. Drug screening revealed that JQ1 inhibited the adhesion of calcium oxalate (CaOx) crystals to TECs. Its therapeutic effect is further confirmed in a glyoxylic acid-induced CaOx crystal deposition mouse model. Utilizing epigenomic and transcriptomic profiling, dynamic enhancer landscape and gene expression program associated with nephrolithiasis are charted. Bioinformatic analysis pinpointing the RXRα as a central transcription factor (TF) modulating enhancer activity. Importantly, the animal studies revealed that RXRα deletion promoted the CaOx crystal deposition, while its activation by Bexarotene (Bex), an FDA-approved drug, mitigated this progression. Mechanistically, under normal circumstances, RXRα inhibited nephrolithiasis-promoting genes by recruiting the HDAC3/SMART complex to repress enhancer activity. Yet, with the progression of CaOx crystal deposition, RXRα expression decreased, leading to enhancer activation and subsequent upregulation of nephrolithiasis-promoting genes. In summary, the work illustrates an epigenetic mechanism underlying TECs fate transition during CaOx crystal deposition and highlights the therapeutic potential of JQ1 and Bex in managing kidney stone diseases.

Most eukaryotic genes express more than one mature mRNA, defined as transcript variants. This complex phenomenon arises from various mechanisms, such as using alternative transcription start sites and alternative post-transcriptional processing events. The resulting transcript variants can lead to synthesizing proteins that possess distinct functional domains or may even generate noncoding RNAs, each with unique roles in cellular processes. The generation of these transcript variants is not merely a random occurrence; it is cell-type specific and varies with developmental stages, aging processes, or pathogenesis of diseases. This highlights the biological significance of transcript variants in regulating gene expression and their potential impact on cellular functionality. Despite the biological importance, investigating transcript variants has been hampered by challenges associated with detecting their expression. This review article addresses the advancements in molecular techniques in detecting transcript variants. Traditional methods such as RT-PCR and RT-qPCR can easily detect known transcript variants using primers that target unique exons associated with the variants. Other techniques like RACE-PCR and hybridization-based methods, including Northern blotting, RNase protection assays, and microarrays, have also been utilized to detect transcript variants. Nevertheless, RNA sequencing (RNA-Seq) has emerged as a powerful technique for identifying transcript variants, especially those with previously unknown sequences. The effectiveness of RNA sequencing in transcript variant detection depends on the specific sequencing approach and the precision of data analysis. By understanding the strengths and weaknesses of each laboratory technique, researchers can develop more effective strategies for detecting mRNA transcript variants. This ability will be crucial for our comprehensive understanding of gene regulation and the implications of transcript diversity in various biological contexts.

This study aimed to use the real-time RT-PCR method to detect the gene expression cytokines IL-10, TNFα, IFN-γ, and TGF-β in the serum of Alzheimer's patients.

This study was conducted on 100 serum samples of Alzheimer's patients. DNA extraction was performed on the samples with the Cinnaclone kit and PCR techniques were used to detect the presence of Helicobacter pylori, Chlamydia trachomatis, and Varicella zoster virus. Real-time RT-PCR was performed to measure the expression of TNFα, IFNɤ, TGFβ, and IL-10 genes with a Smobio kit.

The relative changes in the expression of TNFα, IFNɤ, TGFβ, and IL-10 genes were observed in Alzheimer's patients compared to the control samples without microorganisms, and a significant increase was observed (P < 0.05).

This study showed that the cytokines TNFα, IFNɤ, TGFβ, and IL-10, have an increase in Alzheimer's patients(P < 0.05). Therefore, the presence of the microorganisms accompanied by the rise and inducing the expression of cytokines compared to the groups without the mentioned microorganisms causes a significant increase in the production of cytokines effective in the occurrence or exacerbation of Alzheimer's disease.

Viral nervous necrosis (VNN) is an important viral disease threatening global aquaculture sustainability and affecting over 50 farmed and ecologically important fish species. A major QTL for resistance to VNN has been previously detected in European sea bass, but the underlying causal gene(s) and mutation(s) remain unknown. To identify the mechanisms and genetic factors underpinning resistance to VNN in European sea bass, we employed integrative analyses of multiple functional genomics assays in European sea bass.

The estimated heritability of VNN resistance was high (h2 ~ 0.40), and a major QTL explaining up to 38% of the genetic variance of the trait was confirmed on chromosome 3, with individuals with the resistant QTL genotype showing a 90% survivability against a VNN outbreak. Whole-genome resequencing analyses narrowed the location of this QTL to a small region containing 4 copies of interferon alpha inducible protein 27-like 2A (IFI27L2A) genes, and one copy of the interferon alpha inducible protein 27-like 2 (IFI27L2) gene. RNA sequencing revealed a clear association between the QTL genotype and the expression of two of the IFI27L2A genes, and the IFI27L2 gene. Integration with chromatin accessibility and histone modification data pinpointed two SNPs in active regulatory regions of two of these genes (IFI27L2A and IFI27L2), and transcription factor binding site gains for the resistant alleles were predicted. These alleles, particularly the SNP variant CHR3:10,077,301, exhibited higher frequencies (0.55 to 0.77) in Eastern Mediterranean Sea bass populations, which show considerably higher levels of resistance to VNN, as compared to susceptible West Mediterranean and Atlantic populations (0.15-0.25).

The SNP variant CHR3:10,077,301, through modulation of IFI27L2 and IFI27L2A genes, is likely the causative mutation underlying resistance to VNN in European sea bass. This is one of the first causative mutations discovered for disease resistance traits in fish and paves the way for marker-assisted selection as well as biotechnological approaches to enhance resistance to VNN in European sea bass and other susceptible species.

Metabolic RNA labeling-based time-resolved single-cell RNA sequencing (scRNA-seq) has provided unprecedented tools to dissect the temporal dynamics and the complex gene regulatory networks of gene expression. However, this technology fails to reveal the spatial organization of cells in tissues, which also regulates the gene expression by intercellular communication. Herein, it is demonstrated that integrating time-resolved scRNA-seq with spatial transcriptomics is a new paradigm for spatiotemporal analysis. Metabolic RNA labeling-based time-resolved Well-TEMP-seq is first applied to profile the transcriptional dynamics of glioblastoma (GBM) cells and discover two potential pathways of EZH2-mediated mesenchymal transition in GBM. With spatial transcriptomics, it is further revealed that the crosstalk between CCL2+ malignant cells and IL10+ tumor-associated macrophages in the tumor microenvironment through an EZH2-FOSL2-CCL2 axis contributes to the mesenchymal transition in GBM. These discoveries show the power of integrative spatiotemporal scRNA-seq to elucidate the complex gene regulatory mechanism and advance the understanding of cellular processes in disease.

Variations in blood protein levels have been associated with a broad spectrum of complex diseases, including atherosclerotic cardiovascular disease (ACVD). These associations highlight the intricate interplay between local (e.g., cardiovascular) and systemic (non-cardiovascular) factors for the development of ACVD, emphasizing the need for a comprehensive, systems-level understanding of its etiology. To accomplish this, we developed a causal network inference framework by analyzing one of the largest serum proteomics studies to date, the Age, Gene/Environment Susceptibility-Reykjavik Study (AGES), a prospective population-based study of 7,523 serum proteins measured in 5,376 older adults. To reconstruct a causal network of serum proteins, we used cis -acting protein quantitative trait loci (pQTLs) as instrumental variables to infer causal relationships between protein pairs, while accounting for potential unobserved confounding factors. We identified 185 causal protein subnetworks (FDR = 1%, n ≥ 10 members), which collectively interacted with 5,611 target proteins, offering valuable biological insights and an overview of systemic homeostasis. Several subnetworks, many of which interact to establish a hierarchy of directional relationships, were significantly associated with future myocardial infarction and/or its long-term complications like heart failure, as well as with key cardiometabolic traits that contribute to the onset of ACVD.

Primary aldosteronism (PA), characterized by autonomous aldosterone overproduction, is a major cause of secondary hypertension with significant cardiovascular complications. Current treatments mainly focus on symptom management rather than addressing underlying mechanisms. This study aims to discover novel therapeutic targets for PA using integrated bioinformatics and experimental validation approaches. We employed a systematic approach combining: gene identification through transcriptome-wide association studies (TWAS); causal inference using summary data-based Mendelian randomization (SMR) and two-sample Mendelian randomization (MR) analyses; additional analyses included phenome-wide association analysis, enrichment analysis, protein-protein interaction (PPI) networks, drug repurposing, molecular docking and clinical validation through aldosterone-producing adenomas (APAs) tissue. Through systematic screening and prioritization, we identified 163 PA-associated genes, of which seven emerged as potential drug targets: CEP104, HIP1, TONSL, ZNF100, SHMT1, and two long non-coding RNAs (AC006369.2 and MRPL23-AS1). SHMT1 was identified as the most promising target, showing significantly elevated expression in APAs compared to adjacent non-tumorous tissues. Drug repurposing analysis identified four potential SHMT1-targeting compounds (Mimosine, Pemetrexed, Leucovorin, and Irinotecan), supported by molecular docking studies. The integration of multiple bioinformatics methods and experimental validation successfully identified novel drug targets for hyperaldosteronism. SHMT1, in particular, represents a promising candidate for future therapeutic development. These findings provide new opportunities for developing causative treatments for PA, though further clinical validation is warranted.

With the development of high-throughput profiling techniques, gene expressions have drawn significant attention due to their important biological implications, widespread data availability, and promising biological findings. The complex interactions and regulations among genes naturally lead to a network structure, which can provide a global view of molecular mechanisms and biological processes. This chapter provides a selective overview of constructing gene expression networks and utilizing them in downstream analysis. It also includes a demonstrating example.

Transcriptome-wide association studies (TWAS) have been widely used to identify thousands of likely causal genes for diseases and complex traits using predicted expression models. However, most existing TWAS methods rely on gene expression alone and overlook other regulatory mechanisms of gene expression, including DNA methylation and splicing, that contribute to the genetic basis of these complex traits and diseases. Here we introduce a multi-omics method that integrates gene expression, DNA methylation, and splicing data to improve the identification of associated genes with our traits of interest. Through simulations and by analyzing genome-wide association study (GWAS) summary statistics for 24 complex traits, we show that our integrated method, which leverages these complementary omics biomarkers, achieves higher statistical power, and improves the accuracy of likely causal gene identification in blood tissues over individual omics methods. Finally, we apply our integrated model to a lung cancer GWAS data set, demonstrating the integrated models improved identification of prioritized genes for lung cancer risk.

Gene regulatory networks (GRNs) are crucial for understanding gene regulation and cellular processes. Inferring GRNs helps uncover regulatory pathways, shedding light on the regulation and development of cellular processes. With the rise of high-throughput sequencing and advancements in computational technology, computational models have emerged as cost-effective alternatives to traditional experimental studies. Moreover, the surge in ChIP-seq data for TF-DNA binding has catalyzed the development of graph neural network (GNN)-based methods, greatly advancing GRN inference capabilities. However, most existing GNN-based methods suffer from the inability to capture long-distance structural semantic correlations due to transitive interactions. In this paper, we introduce a novel GNN-based model named Hierarchical Graph Transformer with Contrastive Learning for GRN (HGTCGRN) inference. HGTCGRN excels at capturing structural semantics using a hierarchical graph Transformer, which introduces a series of gene family nodes representing gene functions as virtual nodes to interact with nodes in the GRNS. These semantic-aware virtual-node embeddings are aggregated to produce node representations with varying emphasis. Additionally, we leverage gene ontology information to construct gene interaction networks for contrastive learning optimization of GRNs. Experimental results demonstrate that HGTCGRN achieves superior performance in GRN inference.

Multivariate network-based analytic methods such as weighted gene co-expression network analysis are frequently applied to human and animal gene-expression data to estimate the first principal component of a module, or module eigengene (ME). MEs are interpreted as multivariate summaries of correlated gene-expression patterns and network connectivity across genes within a module. As such, they have the potential to elucidate the mechanisms by which molecular genomic variation contributes to individual differences in complex traits. Although increasingly used to test for associations between modules and complex traits, the genetic and environmental etiology of MEs has not been empirically established. It is unclear if, and to what degree, individual differences in blood-derived MEs reflect random variation versus familial aggregation arising from heritable or shared environmental influences. We used biometrical genetic analyses to estimate the contribution of genetic and environmental influences on MEs derived from blood lymphocytes collected on a sample of N = 661 older male twins from the Vietnam Era Twin Study of Aging (VETSA) whose mean age at assessment was 67.7 years (SD = 2.6 years, range = 62-74 years). Of the 26 detected MEs, 14 (56%) had statistically significant additive genetic variation with an average heritability of 44% (SD = 0.08, range = 35%-64%). Despite the relatively small sample size, this demonstration of significant family aggregation including estimates of heritability in 14 of the 26 MEs suggests that blood-based MEs are reliable and merit further exploration in terms of their associations with complex traits and diseases.

Gene expression regulation is a sophisticated, multi-stage process, and its robustness is critical to normal cell function and the survival of an organism. Previous studies indicate that differential gene expression at the RNA level is typically attenuated at the protein level through translational regulation. However, how post-transcriptional regulation (PTR) influences expression change during the RNA maturation process remains unclear. In this study, we investigated this by quantifying the magnitude of expression change in precursor RNA and mature RNA across a vast range of different biological conditions. We analyzed bulk tissue RNA sequencing data from 4689 samples, including healthy and diseased tissues from human, chimpanzee, rhesus macaque, and murine sources. We demonstrated that PTR tends to support homeostatic expression of mature RNA by amplifying normal tissue-specific expression of precursor RNA, while reducing expression change of precursor RNA in disease contexts. Our study provides insight into the general influence of PTR on gene expression homeostasis. Our analysis also suggests that intronic reads in RNA-seq studies may contain under-utilized information about disease associations. Additionally, our findings may assist in identifying new disease biomarkers and more effective ways of altering gene expression as a therapeutic strategy.

Systemic study of pathogenic pathways and interrelationships underlying genes associated with Alzheimer's disease (AD) facilitates the identification of new targets for effective treatments. Recently available large-scale multiomics datasets provide opportunities to use computational approaches for such studies. Here, we devised a novel disease gene identification (digID) computational framework that consists of a semi-supervised deep learning classifier to predict AD-associated genes and a protein-protein interaction (PPI) network-based analysis to prioritize the importance of these predicted genes in AD. digID predicted 1,529 AD-associated genes and revealed potentially new AD molecular mechanisms and therapeutic targets including GNAI1 and GNB1, two G-protein subunits that regulate cell signaling, and KNG1, an upstream modulator of CDC42 small G-protein signaling and mediator of inflammation and candidate coregulator of amyloid precursor protein (APP). Analysis of mRNA expression validated their dysregulation in AD brains but further revealed the significant spatial patterns in different brain regions as well as among different subregions of the frontal cortex and hippocampi. Super-resolution STochastic Optical Reconstruction Microscopy (STORM) further demonstrated their subcellular colocalization and molecular interactions with APP in a transgenic mouse model of both sexes with AD-like mutations. These studies support the predictions made by digID while highlighting the importance of concurrent biological validation of computationally identified gene clusters as potential new AD therapeutic targets.

The genetic analysis of potato is hampered by the complexity of tetrasomic inheritance. An ongoing effort aims to transform the clonally propagated tetraploid potato into a seed-propagated diploid crop, which would make genetic analyses much easier owing to disomic inheritance. Here, we construct and report the large-scale genetic and heterotic characteristics of a diploid F2 potato population derived from the cross of two highly homozygous inbred lines. We investigate 20,382 traits generated from multi-omics dataset and identify 25,770 quantitative trait loci (QTLs). Coupled with gene expression data, we construct a systems-genetics network for gene discovery in potatoes. Importantly, we explore the genetic basis of heterosis in this population, especially for yield and male fertility heterosis. We find that positive heterotic effects of yield-related QTLs and negative heterotic effects of metabolite QTLs (mQTLs) contribute to yield heterosis. Additionally, we identify a PME gene with a dominance heterotic effect that plays an important role in male fertility heterosis. This study provides genetic resources for the potato community and will facilitate the application of heterosis in diploid potato breeding.

A fundamental challenge in biomedicine is understanding the mechanisms predisposing individuals to disease. While previous research has suggested that switch-like gene expression is crucial in driving biological variation and disease susceptibility, a systematic analysis across multiple tissues is still lacking. By analyzing transcriptomes from 943 individuals across 27 tissues, we identified 1,013 switch-like genes. We found that only 31 (3.1%) of these genes exhibit switch-like behavior across all tissues. These universally switch-like genes appear to be genetically driven, with large exonic genomic structural variants explaining five (~18%) of them. The remaining switch-like genes exhibit tissue-specific expression patterns. Notably, tissue-specific switch-like genes tend to be switched on or off in unison within individuals, likely under the influence of tissue-specific master regulators, including hormonal signals. Among our most significant findings, we identified hundreds of concordantly switched-off genes in the stomach and vagina that are linked to gastric cancer (41-fold, p<10-4) and vaginal atrophy (44-fold, p<10-4), respectively. Experimental analysis of vaginal tissues revealed that low systemic levels of estrogen lead to a significant reduction in both the epithelial thickness and the expression of the switch-like gene ALOX12. We propose a model wherein the switching off of driver genes in basal and parabasal epithelium suppresses cell proliferation therein, leading to epithelial thinning and, therefore, vaginal atrophy. Our findings underscore the significant biomedical implications of switch-like gene expression and lay the groundwork for potential diagnostic and therapeutic applications.

There are no effective pharmacological treatments for sarcopenia. We aim to identify potential therapeutic targets for sarcopenia by integrating various publicly available datasets.

We integrated druggable genome data, cis-eQTL/cis-pQTL from human blood and skeletal muscle tissue, and GWAS summary data of sarcopenia-related traits to analyse the potential causal relationships between drug target genes and sarcopenia using the Mendelian Randomization (MR) method. Sensitivity analyses and Bayesian colocalization were employed to validate the causal relationships. We also assessed the side effects or additional indications of the identified drug targets using a phenome-wide MR (Phe-MR) approach and investigated actionable drugs for target genes using available databases.

MR analysis identified 17 druggable genes with potential causation to sarcopenia in human blood or skeletal muscle tissue. Six of them (HP, HLA-DRA, MAP 3K3, MFGE8, COL15A1, and AURKA) were further confirmed by Bayesian colocalization (PPH4 > 90%). The up-regulation of HP [higher ALM (beta: 0.012, 95% CI: 0.007-0.018, P = 1.2*10-5) and higher grip strength (OR: 0.96, 95% CI: 0.94-0.98, P = 4.2*10-5)], MAP 3K3 [higher ALM (beta: 0.24, 95% CI: 0.21-0.26, P = 1.8*10-94), higher grip strength (OR: 0.82, 95% CI: 0.75-0.90, P = 2.1*10-5), and faster walking pace (beta: 0.03, 95% CI: 0.02-0.05, P = 8.5*10-6)], and MFGE8 [higher ALM (muscle eQTL, beta: 0.09, 95% CI: 0.06-0.11, P = 6.1*10-13; blood pQTL, beta: 0.05, 95% CI: 0.03-0.07, P = 3.8*10-09)], as well as the down-regulation of HLA-DRA [lower ALM (beta: -0.09, 95% CI: -0.11 to -0.08, P = 5.4*10-36) and lower grip strength (OR: 1.13, 95% CI: 1.07-1.20, P = 1.8*10-5)] and COL15A1 [higher ALM (muscle eQTL, beta: -0.07, 95% CI: -0.10 to -0.04, P = 3.4*10-07; blood pQTL, beta: -0.05, 95% CI: -0.06 to -0.03, P = 1.6*10-07)], decreased the risk of sarcopenia. AURKA in blood (beta: -0.16, 95% CI: -0.22 to -0.09, P = 2.1*10-06) and skeletal muscle (beta: 0.03, 95% CI: 0.02 to 0.05, P = 5.3*10-05) tissues showed an inverse relationship with sarcopenia risk. The Phe-MR indicated that the six potential therapeutic targets for sarcopenia had no significant adverse effects. Drug repurposing analysis supported zinc supplementation and collagenase clostridium histolyticum might be potential therapeutics for sarcopenia by activating HP and inhibiting COL15A1, respectively.

Our research indicated MAP 3K3, MFGE8, COL15A1, HP, and HLA-DRA may serve as promising targets for sarcopenia, while the effectiveness of zinc supplementation and collagenase clostridium histolyticum for sarcopenia requires further validation.

Fish gut is a versatile organ serving as the primary pathway for invasion by pathogens, particularly parasites, playing a crucial role in modulating the intestinal adaptive immune response. This study aimed to investigate the cellular-mediated reaction, mucosal acidity, and the expression of proliferating cell nuclear antigen (PCNA), vascular endothelial growth factor (VEGF), and CD68 in the intestines of catfish, Clarias gariepinus, naturally infected with helminths. Forty catfish were collected from the Nile River and examined for intestinal parasites. The intestinal tissues of the control and infected fish were fixed for histochemical and immunohistochemical studies. Two groups of helminths were found: cestodes Tetracampos ciliotheca and Polyonchobothrium clarias, and nematodes Paracamallanus cyathopharynx, with a prevalence rate of 63.63%, 18.0%, and 18.0%, respectively. Our results showed that the infected fish had a statistically significant rise in the activity of immune cells, including mast cells, eosinophil granular cells, and dendritic cells. This correlated with upregulation in the expressions of PCNA, VEGF, and CD68. Histochemical analyses demonstrated a marked increase in acidic mucus production, Sudan black B, and bromophenol mercury blue. This study enriches our understanding of the evolution of vertebrate immunity in combating intestinal parasitic infections and the host's adaptive responses.

Leaf scald, caused by Xanthomonas albilineans, is a severe disease affecting sugarcane worldwide. One of the most practical ways to control it is by developing resistant sugarcane cultivars. It is essential to identify genes associated with the response to leaf scald. A panel of 170 sugarcane genotypes was evaluated for resistance to leaf scald in field conditions for 2 years, followed by a 1-year greenhouse experiment. The phenotypic evaluation data showed a wide continuous distribution, with heritability values ranging from 0.58 to 0.84. Thirteen single nucleotide polymorphisms (SNPs) were identified, significantly associated with leaf scald resistance. Among these, eight were stable across multiple environments and association models. The candidate genes identified and validated based on RNA-seq and qRT-PCR included two genes that encode NB-ARC leucine-rich repeat (LRR)-containing domain disease-resistance protein. These findings provide a basis for developing marker-assisted selection strategies in sugarcane breeding programs.

Understanding the intricate relationships between diseases is critical for both prevention and recovery. However, there is a lack of suitable methodologies for exploring the precedence relationships within multiple censored time-to-event data, resulting in decreased analytical accuracy. This study introduces the Censored Event Precedence Analysis (CEPA), which is a nonparametric Bayesian approach suitable for understanding the precedence relationships in censored multivariate events. CEPA aims to analyze the precedence relationships between events to predict subsequent occurrences effectively. We applied CEPA to neonatal data from the National Health Insurance Service, identifying the precedence relationships among the seven most commonly diagnosed diseases categorized by the International Classification of Diseases. This analysis revealed a typical diagnostic sequence, starting with respiratory diseases, followed by skin, infectious, digestive, ear, eye, and injury-related diseases. Furthermore, simulation studies were conducted to demonstrate CEPA suitability for censored multivariate datasets compared to traditional models. The performance accuracy reached 76% for uniform distribution and 65% for exponential distribution, showing superior performance in all four tested environments. Therefore, the statistical approach based on CEPA enhances our understanding of disease interrelationships beyond competitive methodologies. By identifying disease precedence with CEPA, we can preempt subsequent disease occurrences and propose a healthcare system based on these relationships.

Powdery mildew disease (PMD) is caused by the obligate biotrophic fungus Microsphaera diffusa Cooke & Peck (M. diffusa) and results in significant yield losses in soybean (Glycine max (L.) Merr.) crops. By identifying disease-resistant genes and breeding soybean accessions with enhanced resistance, we can effectively mitigate the detrimental impact of PMD on soybeans. We analyzed PMD resistance in a diversity panel of 315 soybean accessions in two locations over 3 years, and candidate genes associated with PMD resistance were identified through genome-wide association studies (GWAS), haplotype analysis, qRT-PCR, and EMS mutant analysis. Based on the GWAS approach, we identified a region on chromosome 16 (Chr16) in which 21 genes form a gene cluster that is highly correlated with PMD resistance. In order to validate and refine these findings, we conducted haplotype analysis of 21 candidate genes and indicated there are single nucleotide polymorphisms (SNPs) and insertion-deletions (InDels) variations of Glyma.16G214000, Glyma.16G214200, Glyma.16G215100 and Glyma.16G215300 within the coding and promoter regions that exhibit a strong association with resistance against PMD. Subsequent structural analysis of candidate genes within this cluster revealed that in 315 accessions, the majority of accessions exhibited resistance to PMD when Glyma.16G214300, Glyma.16G214800 and Glyma.16G215000 were complete; however, they demonstrated susceptibility to PMD when these genes were incomplete. Quantitative real-time PCR assays (qRT-PCR) of possible candidate genes showed that 14 candidate genes (Glyma.16G213700, Glyma.16G213800, Glyma.16G213900, Glyma.16G214000, Glyma.16G214200, Glyma.16G214300, Glyma.16G214500, Glyma.16G214585, Glyma.16G214669, Glyma.16G214700, Glyma.16G214800, Glyma.16G215000, Glyma.16G215100 and Glyma.16G215300) were involved in PMD resistance. Finally, we evaluated the PMD resistance of mutant lines from the Williams 82 EMS mutations library, which revealed that mutants of Glyma.16G214000, Glyma.16G214200, Glyma.16G214300, Glyma.16G214800, Glyma.16G215000, Glyma.16G215100 and Glyma.16G215300, exhibited sensitivity to PMD. Combined with the analysis results of GWAS, haplotypes, qRT-PCR and mutants, the genes Glyma.16G214000, Glyma.16G214200, Glyma.16G214300, Glyma.16G214800, Glyma.16G215000, Glyma.16G215100 and Glyma.16G215300 were identified as highly correlated with PMD resistance. The candidate genes identified above are all NLR family genes, and these discoveries deepen our understanding of the molecular basis of PMD resistance in soybeans and will be useful for guiding breeding strategies.

Copper is one of a handful of biologically necessary heavy metals that is also a common environmental pollutant. Under normal conditions, copper ions are required for many key physiological processes. However, in excess, copper results in cell and tissue damage ranging in severity from temporary injury to permanent neurological damage. Because of its biological relevance, and because many conserved copper-responsive genes respond to nonessential heavy metal pollutants, copper resistance in Drosophila melanogaster is a useful model system with which to investigate the genetic control of the heavy metal stress response. Because heavy metal toxicity has the potential to differently impact specific tissues, we genetically characterized the control of the gene expression response to copper stress in a tissue-specific manner in this study. We assessed the copper stress response in head and gut tissue of 96 inbred strains from the Drosophila Synthetic Population Resource using a combination of differential expression analysis and expression quantitative trait locus mapping. Differential expression analysis revealed clear patterns of tissue-specific expression. Tissue and treatment specific responses to copper stress were also detected using expression quantitative trait locus mapping. Expression quantitative trait locus associated with MtnA, Mdr49, Mdr50, and Sod3 exhibited both genotype-by-tissue and genotype-by-treatment effects on gene expression under copper stress, illuminating tissue- and treatment-specific patterns of gene expression control. Together, our data build a nuanced description of the roles and interactions between allelic and expression variation in copper-responsive genes, provide valuable insight into the genomic architecture of susceptibility to metal toxicity, and highlight candidate genes for future functional characterization.

Human death is a complex, time-governed phenomenon that leads to the irreversible cessation of all bodily functions. Recent molecular and genetic studies have revealed remarkable experimental evidence of genetically programmed cellular death characterized by several physiological processes; however, the basic physiological function that occurs during the immediate postmortem period remains inadequately described. There is a paucity of knowledge connecting necrotic pathologies occurring in human organ tissues to complete functional loss of the human organism. Cells, tissues, organs, and organ systems show a range of differential resilience and endurance responses that occur during organismal death. Intriguingly, a persistent ambiguity in the study of postmortem physiological systems is the determination of the trajectory of a complex multicellular human body, far from life-sustaining homeostasis, following the gradual or sudden expiry of its regulatory systems. Recent groundbreaking investigations have resulted in a paradigm shift in understanding the cell biology and physiology of death. Two significant findings are that (i) most cells in the human body are microbial, and (ii) microbial cell abundance significantly increases after death. By addressing the physiological as well as the microbiological aspects of death, future investigations are poised to reveal innovative insights into the enigmatic biological activities associated with death and human decomposition. Understanding the elaborate crosstalk of abiotic and biotic factors in the context of death has implications for scientific discoveries important to informing translational knowledge regarding the transition from living to the non-living. There are important and practical needs for a transformative reestablishment of accepted models of biological death (i.e., artificial intelligence, AI) for more precise determinations of when the regulatory mechanisms for homeostasis of a living individual have ceased. In this review, we summarize mechanisms of physiological, genetic, and microbiological processes that define the biological changes and pathways associated with human organismal death and decomposition.

Transcriptomic data is a treasure trove in modern molecular biology, as it offers a comprehensive viewpoint into the intricate nuances of gene expression dynamics underlying biological systems. This genetic information must be utilized to infer biomolecular interaction networks that can provide insights into the complex regulatory mechanisms underpinning the dynamic cellular processes. Gene regulatory networks and protein-protein interaction networks are two major classes of such networks. This chapter thoroughly investigates the wide range of methodologies used for distilling insightful revelations from transcriptomic data that include association-based methods (based on correlation among expression vectors), probabilistic models (using Bayesian and Gaussian models), and interologous methods. We reviewed different approaches for evaluating the significance of interactions based on the network topology and biological functions of the interacting molecules and discuss various strategies for the identification of functional modules. The chapter concludes with highlighting network-based techniques of prioritizing key genes, outlining the centrality-based, diffusion- based, and subgraph-based methods. The chapter provides a meticulous framework for investigating transcriptomic data to uncover assembly of complex molecular networks for their adaptable analyses across a broad spectrum of biological domains.

Expression quantitative trait locus (eQTL) analysis is a useful tool to identify genetic loci that are associated with gene expression levels. Large collaborative efforts such as the Genotype-Tissue Expression (GTEx) project provide valuable resources for eQTL analysis in different tissues. Most existing methods, however, either focus on one tissue at a time, or analyze multiple tissues to identify eQTLs jointly present in multiple tissues. There is a lack of powerful methods to identify eQTLs in a target tissue while effectively borrowing strength from auxiliary tissues. In this paper, we propose a novel statistical framework to improve the eQTL detection efficacy in the tissue of interest with auxiliary information from other tissues. This framework can enhance the power of the hypothesis test for eQTL effects by incorporating shared and specific effects from multiple tissues into the test statistics. We also devise data-driven and distributed computing approaches for efficient implementation of eQTL detection when the number of tissues is large. Numerical studies in simulation demonstrate the efficacy of the proposed method, and the real data analysis of the GTEx example provides novel insights into eQTL findings in different tissues.

We have postulated that common changes in gene expression after treatment with different therapeutic classes of psychotropic drugs contribute to their common therapeutic mechanisms of action.

To test this hypothesis, we measured levels of cortical coding and non-coding RNA using GeneChip® Rat Exon 1.0 ST Array after treatment with vehicle (chow only), chow containing 1.8 g lithium carbonate/kg (n = 10) or chow containing 12 g sodium valproate/kg (n = 10) for 28 days. Differences in levels of RNA were identified using JMP Genomics 13 and the Panther Gene Ontology Classification System was used to identify potential consequences of RNA.

Compared to vehicle treatment, levels of cortical RNA for 543 and 583 coding and non-coding RNAs were different after treatment with valproate and lithium, respectively. Moreover, levels of 323 coding and non-coding RNAs were altered in a highly correlated way by treatment with valproate and lithium, changes that would impact on cholinergic, glutamatergic, serotonergic and dopaminergic neurotransmission as well as on voltage gated ion channels.

Our study suggests that treating with mood stabilisers cause many common changes in levels of RNA which will impact on CNS function, particularly affecting post-synaptic muscarinic receptor functioning and the release of multiple neurotransmitters.

Fusarium oxysporum is a prevalent fungal pathogen that diminishes soybean yield through seedling disease and root rot. Preventing Fusarium oxysporum root rot (FORR) damage entails on the identification of resistance genes and developing resistant cultivars. Therefore, conducting fine mapping and marker development for FORR resistance genes is of great significance for fostering the cultivation of resistant varieties. In this study, 350 soybean germplasm accessions, mainly from Northeast China, underwent genotyping using the SoySNP50K Illumina BeadChip, which includes 52,041 single nucleotide polymorphisms (SNPs). Their resistance to FORR was assessed in a greenhouse. Genome-wide association studies utilizing the general linear model, mixed linear model, compressed mixed linear model, and settlement of MLM under progressively exclusive relationship models were conducted to identify marker-trait associations while effectively controlling for population structure.

The results demonstrated that these models effectively managed population structure. Eight SNP loci significantly associated with FORR resistance in soybean were detected, primarily located on Chromosome 6. Notably, there was a strong linkage disequilibrium between the large-effect SNPs ss715595462 and ss715595463, contributing substantially to phenotypic variation. Within the genetic interval encompassing these loci, 28 genes were present, with one gene Glyma.06G088400 encoding a protein kinase family protein containing a leucine-rich repeat domain identified as a potential candidate gene in the reference genome of Williams82. Additionally, quantitative real-time reverse transcription polymerase chain reaction analysis evaluated the gene expression levels between highly resistant and susceptible accessions, focusing on primary root tissues collected at different time points after F. oxysporum inoculation. Among the examined genes, only this gene emerged as the strongest candidate associated with FORR resistance.

The identification of this candidate gene Glyma.06G088400 improves our understanding of soybean resistance to FORR and the markers strongly linked to resistance can be beneficial for molecular marker-assisted selection in breeding resistant soybean accessions against F. oxysporum.

Molecular omics data is increasingly ubiquitous throughout medicine. In organ transplantation, recent large-scale research efforts are generating the 'transplant-ome' - the entire set of molecular omics data, including the genome, transcriptome, proteome, and metabolome. Importantly, early studies in anesthesiology have demonstrated how perioperative interventions alter molecular profiles in various patient populations. The next step for anesthesiologists and intensivists will be to tailor perioperative care to the transplant-ome of individual liver transplant patients.

In liver transplant patients, elements of the transplant-ome predict complications and point to novel interventions. Importantly, molecular profiles of both the donor organ and recipient contribute to this risk, and interventions like normothermic machine perfusion influence these profiles. As we can now measure various omics molecules simultaneously, we can begin to understand how these molecules interact to form molecular networks and emerging technologies offer noninvasive and continuous ways to measure these networks throughout the perioperative period. Molecules that regulate these networks are likely mediators of complications and actionable clinical targets throughout the perioperative period.

The transplant-ome can be used to tailor perioperative care to the individual liver transplant patient. Monitoring molecular networks continuously and noninvasively would provide new opportunities to optimize perioperative management.

The integration of data from multiple genomic assays from humans and non-human model organisms is an effective approach to identify genes involved in skeletal fragility and fracture risk due to osteoporosis and other conditions. This review summarizes genome-wide genetic variation and gene expression data resources relevant to the discovery of genes contributing to skeletal fragility and fracture risk.

Genome-wide association studies (GWAS) of osteoporosis-related traits are summarized, in addition to gene expression in bone tissues in humans and non-human organisms, with a focus on rodent models related to skeletal fragility and fracture risk. Gene discovery approaches using these genomic data resources are described. We also describe the Musculoskeletal Knowledge Portal (MSKKP) that integrates much of the available genomic data relevant to fracture risk. The available genomic resources provide a wealth of knowledge and can be analyzed to identify genes related to fracture risk. Genomic resources that would fill particular scientific gaps are discussed.

Powdery mildew (PMD), caused by the pathogen Microsphaera diffusa, leads to substantial yield decreases in susceptible soybean under favorable environmental conditions. Effective prevention of soybean PMD damage can be achieved by identifying resistance genes and developing resistant cultivars. In this study, we genotyped 331 soybean germplasm accessions, primarily from Northeast China, using the SoySNP50K BeadChip, and evaluated their resistance to PMD in a greenhouse setting. To identify marker-trait associations while effectively controlling for population structure, we conducted genome-wide association studies utilizing factored spectrally transformed linear mixed models, mixed linear models, efficient mixed-model association eXpedited, and compressed mixed linear models. The results revealed seven single nucleotide polymorphism (SNP) loci strongly associated with PMD resistance in soybean. Among these, one SNP was localized on chromosome (Chr) 14, and six SNPs with low linkage disequilibrium were localized near or in the region of previously mapped genes on Chr 16. In the reference genome of Williams82, we discovered 96 genes within the candidate region, including 17 resistance (R)-like genes, which were identified as potential candidate genes for PMD resistance. In addition, we performed quantitative real-time reverse transcription polymerase chain reaction analysis to evaluate the gene expression levels in highly resistant and susceptible genotypes, focusing on leaf tissues collected at different times after M. diffusa inoculation. Among the examined genes, three R-like genes, including Glyma.16G210800, Glyma.16G212300, and Glyma.16G213900, were identified as strong candidates associated with PMD resistance. This discovery can significantly enhance our understanding of soybean resistance to PMD. Furthermore, the significant SNPs strongly associated with resistance can serve as valuable markers for genetic improvement in breeding M. diffusa-resistant soybean cultivars.

Common human diseases result from the interplay of genes and their biologically associated pathways. Genetic pathway analyses provide more biological insight as compared to conventional gene-based analysis. In this article, we propose a framework combining genetic data into pathway structure and using an ensemble of convolutional neural networks (CNNs) along with a Canonical Correlation Regularizer layer for comprehensive prediction of disease risk. The novelty of our approach lies in our two-step framework: (i) utilizing the CNN's effectiveness to extract the complex gene associations within individual genetic pathways and (ii) fusing features from ensemble of CNNs through Canonical Correlation Regularization layer to incorporate the interactions between pathways which share common genes. During prediction, we also address the important issues of interpretability of neural network models, and identifying the pathways and genes playing an important role in prediction.

Implementation of our methodology into three real cancer genetic datasets for different prediction tasks validates our model's generalizability and robustness. Comparing with conventional models, our methodology provides consistently better performance with AUC improvement of 11% on predicting early/late-stage kidney cancer, 10% on predicting kidney versus liver cancer type and 7% on predicting survival status in ovarian cancer as compared to the next best conventional machine learning model. The robust performance of our deep learning algorithm indicates that disease prediction using neural networks in multiple functionally related genes across different pathways improves genetic data-based prediction and understanding molecular mechanisms of diseases.

https://github.com/divya031090/ReGeNNe.

Transcription Factors (TFs) influence gene expression by facilitating or disrupting the formation of transcription initiation machinery at particular genomic loci. Because genomic localization of TFs is in part driven by TF recognition of DNA sequence, variation in TF binding sites can disrupt TF-DNA associations and affect gene regulation. To identify variants that impact TF binding in human brain tissues, we quantified allele bias for 93 TFs analyzed with ChIP-seq experiments of multiple structural brain regions from two donors. Using graph genomes constructed from phased genomic sequence data, we compared ChIP-seq signal between alleles at heterozygous variants within each tissue sample from each donor. Comparison of results from different brain regions within donors and the same regions between donors provided measures of allele bias reproducibility. We identified thousands of DNA variants that show reproducible bias in ChIP-seq for at least one TF. We found that alleles that are rarer in the general population were more likely than common alleles to exhibit large biases, and more frequently led to reduced TF binding. Combining ChIP-seq with RNA-seq, we identified TF-allele interaction biases with RNA bias in a phased allele linked to 6,709 eQTL variants identified in GTEx data, 3,309 of which were found in neural contexts. Our results provide insights into the effects of both common and rare variation on gene regulation in the brain. These findings can facilitate mechanistic understanding of cis-regulatory variation associated with biological traits, including disease.

Phenotypic variation in human is the results of genetic variation and environmental influences. Understanding the contribution of genetic and environmental components to phenotypic variation is of great interest. The variance explained by genome-wide single nucleotide polymorphisms (SNPs) typically represents a small proportion of the phenotypic variance for complex traits, which may be because the genome is only a part of the whole biological process to shape the phenotypes. In this study, we propose to partition the phenotypic variance of three anthropometric traits, using gene expression levels and environmental variables from GTEx data. We use the gene expression of four tissues that are deemed relevant for the anthropometric traits (two adipose tissues, skeletal muscle tissue and blood tissue). Additionally, we estimate the transcriptome-environment correlation that partly underlies the phenotypes of the anthropometric traits. We found that genetic factors play a significant role in determining body mass index (BMI), with the proportion of phenotypic variance explained by gene expression levels of visceral adipose tissue being 0.68 (SE = 0.06). However, we also observed that environmental factors such as age, sex, ancestry, smoking status, and drinking alcohol status have a small but significant impact (0.005, SE = 0.001). Interestingly, we found a significant negative correlation between the transcriptomic and environmental effects on BMI (transcriptome-environment correlation = -0.54, SE = 0.14), suggesting an antagonistic relationship. This implies that individuals with lower genetic profiles may be more susceptible to the effects of environmental factors on BMI, while those with higher genetic profiles may be less susceptible. We also show that the estimated transcriptomic variance varies across tissues, e.g., the gene expression levels of whole blood tissue and environmental variables explain a lower proportion of BMI phenotypic variance (0.16, SE = 0.05 and 0.04, SE = 0.004 respectively). We observed a significant positive correlation between transcriptomic and environmental effects (1.21, SE = 0.23) for this tissue. In conclusion, phenotypic variance partitioning can be done using gene expression and environmental data even with a small sample size (n = 838 from GTEx data), which can provide insights into how the transcriptomic and environmental effects contribute to the phenotypes of the anthropometric traits.

Copper is one of a handful of biologically necessary heavy metals that is also a common environmental pollutant. Under normal conditions, copper ions are required for many key physiological processes. However, in excess, copper quickly results in cell and tissue damage that can range in severity from temporary injury to permanent neurological damage. Because of its biological relevance, and because many conserved copper-responsive genes also respond to other non-essential heavy metal pollutants, copper resistance in Drosophila melanogaster is a useful model system with which to investigate the genetic control of the response to heavy metal stress. Because heavy metal toxicity has the potential to differently impact specific tissues, we genetically characterized the control of the gene expression response to copper stress in a tissue-specific manner in this study. We assessed the copper stress response in head and gut tissue of 96 inbred strains from the Drosophila Synthetic Population Resource (DSPR) using a combination of differential expression analysis and expression quantitative trait locus (eQTL) mapping. Differential expression analysis revealed clear patterns of tissue-specific expression, primarily driven by a more pronounced gene expression response in gut tissue. eQTL mapping of gene expression under control and copper conditions as well as for the change in gene expression following copper exposure (copper response eQTL) revealed hundreds of genes with tissue-specific local cis-eQTL and many distant trans-eQTL. eQTL associated with MtnA, Mdr49, Mdr50, and Sod3 exhibited genotype by environment effects on gene expression under copper stress, illuminating several tissue- and treatment-specific patterns of gene expression control. Together, our data build a nuanced description of the roles and interactions between allelic and expression variation in copper-responsive genes, provide valuable insight into the genomic architecture of susceptibility to metal toxicity, and highlight many candidate genes for future functional characterization.

The continuous worldwide seagrasses decline calls for immediate actions in order to preserve this precious marine ecosystem. The main stressors that have been linked with decline in seagrasses are 1) the increasing ocean temperature due to climate change and 2) the continuous inputs of nutrients (eutrophication) associated with coastal human activities. To avoid the loss of seagrass populations, an "early warning" system is needed. We used Weighed Gene Co-expression Network Analysis (WGCNA), a systems biology approach, to identify potential candidate genes that can provide an early warning signal of stress in the Mediterranean iconic seagrass Posidonia oceanica, anticipating plant mortality. Plants were collected from both eutrophic (EU) and oligotrophic (OL) environments and were exposed to thermal and nutrient stress in a dedicated mesocosm. By correlating the whole-genome gene expression after 2-weeks exposure with the shoot survival percentage after 5-weeks exposure to stressors, we were able to identify several transcripts that indicated an early activation of several biological processes (BP) including: protein metabolic process, RNA metabolic process, organonitrogen compound biosynthetic process, catabolic process and response to stimulus, which were shared among OL and EU plants and among leaf and shoot apical meristem (SAM), in response to excessive heat and nutrients. Our results suggest a more dynamic and specific response of the SAM compared to the leaf, especially the SAM from plants coming from a stressful environment appeared more dynamic than the SAM from a pristine environment. A vast list of potential molecular markers is also provided that can be used as targets to assess field samples.

Despite decades of genetic studies on late-onset Alzheimer's disease, the underlying molecular mechanisms remain unclear. To better comprehend its complex etiology, we use an integrative approach to build robust predictive (causal) network models using two large human multi-omics datasets. We delineate bulk-tissue gene expression into single cell-type gene expression and integrate clinical and pathologic traits, single nucleotide variation, and deconvoluted gene expression for the construction of cell type-specific predictive network models. Here, we focus on neuron-specific network models and prioritize 19 predicted key drivers modulating Alzheimer's pathology, which we then validate by knockdown in human induced pluripotent stem cell-derived neurons. We find that neuronal knockdown of 10 of the 19 targets significantly modulates levels of amyloid-beta and/or phosphorylated tau peptides, most notably JMJD6. We also confirm our network structure by RNA sequencing in the neurons following knockdown of each of the 10 targets, which additionally predicts that they are upstream regulators of REST and VGF. Our work thus identifies robust neuronal key drivers of the Alzheimer's-associated network state which may represent therapeutic targets with relevance to both amyloid and tau pathology in Alzheimer's disease.

Technological advancement, data democratisation, and decreasing costs have led to a revolution in molecular biology in which the entire set of DNA, RNA, proteins, and various other molecules - the 'multi-omic' profile - can be measured in humans. Sequencing 1 million bases of human DNA now costs US$0.01, and emerging technologies soon promise to reduce the cost of sequencing the whole genome to US$100. These trends have made it feasible to sample the multi-omic profile of millions of people, much of which is publicly available for medical research. Can anaesthesiologists use these data to improve patient care? This narrative review brings together a rapidly growing literature in multi-omic profiling across numerous fields that points to the future of precision anaesthesiology. Here, we discuss how DNA, RNA, proteins, and other molecules interact in molecular networks that can be used for preoperative risk stratification, intraoperative optimisation, and postoperative monitoring. This literature provides evidence for four fundamental insights: (1) Clinically similar patients have different molecular profiles and, as a consequence, different outcomes. (2) Vast, publicly available, and rapidly growing molecular datasets have been generated in chronic disease patients and can be repurposed to estimate perioperative risk. (3) Multi-omic networks are altered in the perioperative period and influence postoperative outcomes. (4) Multi-omic networks can serve as empirical, molecular measurements of a successful postoperative course. With this burgeoning universe of molecular data, the anaesthesiologist-of-the-future will tailor their clinical management to an individual's multi-omic profile to optimise postoperative outcomes and long-term health.

Cholangiocytes transit from quiescence to hyperproliferation during cystogenesis in polycystic liver disease (PLD), the severity of which displays prominent sex differences. Epigenetic regulation plays important roles in cell state transition. We aimed to investigate the sex-specific epigenetic basis of hepatic cystogenesis and to develop therapeutic strategies targeting epigenetic modifications for PLD treatment.

Normal and cystic primary cholangiocytes were isolated from wild-type and PLD mice of both sexes. Chromatin states were characterized by analyzing chromatin accessibility (ATAC sequencing) and multiple histone modifications (chromatin immunoprecipitation sequencing). Differential gene expression was determined by transcriptomic analysis (RNA sequencing). Pharmacologic inhibition of epigenetic modifying enzymes was undertaken in PLD model mice.

Through genome-wide profiling of chromatin dynamics, we revealed a profound increase of global chromatin accessibility during cystogenesis in both male and female PLD cholangiocytes. We identified a switch from H3K9me3 to H3K9ac on cis-regulatory DNA elements of cyst-associated genes and showed that inhibition of H3K9ac acetyltransferase or H3K9me3 demethylase slowed cyst growth in male, but not female, PLD mice. In contrast, we found that H3K27ac was specifically increased in female PLD mice and that genes associated with H3K27ac-gained regions were enriched for cyst-related pathways. In an integrated epigenomic and transcriptomic analysis, we identified an estrogen receptor alpha-centered transcription factor network associated with the H3K27ac-regulated cystogenic gene expression program in female PLD mice.

Our findings highlight the multi-layered sex-specific epigenetic dynamics underlying cholangiocyte state transition and reveal a potential epigenetic therapeutic strategy for male PLD patients.

In the present study, we elucidate a sex-specific epigenetic mechanism underlying the cholangiocyte state transition during hepatic cystogenesis and identify epigenetic drugs that effectively slow cyst growth in male PLD mice. These findings underscore the importance of sex difference in the pathogenesis of PLD and may guide researchers and physicians to develop sex-specific personalized approaches for PLD treatment.

Expanding a statistical approach called Mendelian randomization to include multiple variables may help researchers to identify new molecular causes of specific traits.

Single-cell RNA sequencing (scRNA-seq) reveals the transcriptional heterogeneity of cells, but the static snapshots fail to reveal the time-resolved dynamics of transcription. Herein, we develop Well-TEMP-seq, a high-throughput, cost-effective, accurate, and efficient method for massively parallel profiling the temporal dynamics of single-cell gene expression. Well-TEMP-seq combines metabolic RNA labeling with scRNA-seq method Well-paired-seq to distinguish newly transcribed RNAs marked by T-to-C substitutions from pre-existing RNAs in each of thousands of single cells. The Well-paired-seq chip ensures a high single cell/barcoded bead pairing rate (~80%) and the improved alkylation chemistry on beads greatly alleviates chemical conversion-induced cell loss (~67.5% recovery). We further apply Well-TEMP-seq to profile the transcriptional dynamics of colorectal cancer cells exposed to 5-AZA-CdR, a DNA-demethylating drug. Well-TEMP-seq unbiasedly captures the RNA dynamics and outperforms the splicing-based RNA velocity method. We anticipate that Well-TEMP-seq will be broadly applicable to unveil the dynamics of single-cell gene expression in diverse biological processes.

Global genetic networks provide additional information for the analysis of human diseases, beyond the traditional analysis that focuses on single genes or local networks. The Gaussian graphical model (GGM) is widely applied to learn genetic networks because it defines an undirected graph decoding the conditional dependence between genes. Many algorithms based on the GGM have been proposed for learning genetic network structures. Because the number of gene variables is typically far more than the number of samples collected, and a real genetic network is typically sparse, the graphical lasso implementation of GGM becomes a popular tool for inferring the conditional interdependence among genes. However, graphical lasso, although showing good performance in low dimensional data sets, is computationally expensive and inefficient or even unable to work directly on genome-wide gene expression data sets. In this study, the method of Monte Carlo Gaussian graphical model (MCGGM) was proposed to learn global genetic networks of genes. This method uses a Monte Carlo approach to sample subnetworks from genome-wide gene expression data and graphical lasso to learn the structures of the subnetworks. The learned subnetworks are then integrated to approximate a global genetic network. The proposed method was evaluated with a relatively small real data set of RNA-seq expression levels. The results indicate the proposed method shows a strong ability of decoding the interactions with high conditional dependences among genes. The method was then applied to genome-wide data sets of RNA-seq expression levels. The gene interactions with high interdependence from the estimated global networks show that most of the predicted gene-gene interactions have been reported in the literatures playing important roles in different human cancers. Also, the results validate the ability and reliability of the proposed method to identify high conditional dependences among genes in large-scale data sets.

Our knowledge of the cell-type-specific mechanisms of insulin resistance remains limited. To dissect the cell-type-specific molecular signatures of insulin resistance, we performed a multiscale gene network analysis of adipose and muscle tissues in African and European ancestry populations. In adipose tissues, a comparative analysis revealed ethnically conserved cell-type signatures and two adipocyte subtype-enriched modules with opposite insulin sensitivity responses. The modules enriched for adipose stem and progenitor cells as well as immune cells showed negative correlations with insulin sensitivity. In muscle tissues, the modules enriched for stem cells and fibro-adipogenic progenitors responded to insulin sensitivity oppositely. The adipocyte and muscle fiber-enriched modules shared cellular-respiration-related genes but had tissue-specific rearrangements of gene regulations in response to insulin sensitivity. Integration of the gene co-expression and causal networks further pinpointed key drivers of insulin resistance. Together, this study revealed the cell-type-specific transcriptomic networks and signaling maps underlying insulin resistance in major glucose-responsive tissues. A record of this paper's transparent peer review process is included in the supplemental information.

Expression quantitative trait locus detection has become increasingly important for understanding how noncoding variants contribute to disease susceptibility and complex traits. The major challenges in expression quantitative trait locus fine-mapping and causal variant discovery relate to the impact of linkage disequilibrium on signals due to one or multiple functional variants that lie within a credible set. We perform expression quantitative trait locus fine-mapping using the all-but-one approach, conditioning each signal on all others detected in an interval, on the Consortium for the Architecture of Gene Expression cohorts of microarray-based peripheral blood gene expression in 2,138 European-ancestry human adults. We contrast these results with traditional forward stepwise conditional analysis and a Bayesian localization method. All-but-one conditioning significantly modifies effect-size estimates for 51% of 2,351 expression quantitative trait locus peaks, but only modestly affects credible set size and location. On the other hand, both conditioning approaches result in unexpectedly low overlap with Bayesian credible sets, with just 57% peak concordance and between 50% and 70% SNP sharing, leading us to caution against the assumption that any one localization method is superior to another. We also cross reference our results with ATAC-seq data, cell-type-specific expression quantitative trait locus, and activity-by-contact-enhancers, leading to the proposal of a 5-tier approach to further reduce credible set sizes and prioritize likely causal variants for all known inflammatory bowel disease risk loci active in immune cells.