The role of regulatory variation in complex traits and disease

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

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

Crosstalk between epitranscriptomic and epigenomic modifications and its implication in human diseases

Crosstalk between N6-methyladenosine (m6A) and epigenomes is crucial for gene regulation, but its regulatory directionality and disease significance remain unclear. Here, we utilize quantitative trait loci (QTLs) as genetic instruments to delineate directional maps of crosstalk between m6A and two epigenomic traits, DNA methylation (DNAme) and H3K27ac. We identify 47 m6A-to-H3K27ac and 4,733 m6A-to-DNAme and, in the reverse direction, 106 H3K27ac-to-m6A and 61,775 DNAme-to-m6A regulatory loci, with differential genomic location preference observed for different regulatory directions. Integrating these maps with complex diseases, we prioritize 20 genome-wide association study (GWAS) loci for neuroticism, depression, and narcolepsy in brain; 1,767 variants for asthma and expiratory flow traits in lung; and 249 for coronary artery disease, blood pressure, and pulse rate in muscle. This study establishes disease regulatory paths, such as rs3768410-DNAme-m6A-asthma and rs56104944-m6A-DNAme-hypertension, uncovering locus-specific crosstalk between m6A and epigenomic layers and offering insights into regulatory circuits underlying human diseases.

Multi-omics analysis unravels chemical roadmap and genetic basis for peach fruit aroma improvement

Selection of fruits with enhanced health benefits and superior flavor is an important aspect of peach breeding. Understanding the genetic interplay between appearance and flavor chemicals remains a major challenge. We identify the most important volatiles contributing to consumer preferences for peach, thus establishing priorities for improving flavor quality. We quantify volatiles of a peach population consisting of 184 accessions and demonstrate major reductions in the important flavor volatiles linalool and Z-3-hexenyl acetate in red-fleshed accessions. We identify 474 functional gene regulatory networks (GRNs), among which GRN05 plays a crucial role in controlling both red flesh and volatile content through the NAM/ATAF1/2/CUC (NAC) transcription factor PpBL. Overexpressing PpBL results in reduced expression of PpNAC1, a positive regulator for Z-3-hexenyl acetate and linalool synthesis. Additionally, we identify haplotypes for three tandem PpAATs that are significantly correlated with reduced gene expression and ester content. We develop genetic resources for improvement of fruit quality.

Unraveling the molecular mechanisms of Ace2-mediated post-COVID-19 cognitive dysfunction through systems genetics approach

The dysregulation of Angiotensin-converting enzyme 2 (ACE2) in central nervous system is believed associates with COVID-19 induced cognitive dysfunction. However, the detailed mechanism remains largely unknown. In this study, we performed a comprehensive system genetics analysis on hippocampal ACE2 based on BXD mice panel. Expression quantitative trait loci (eQTLs) mapping showed that Ace2 was strongly trans-regulated, and the elevation of Ace2 expression level was significantly correlated with impaired cognitive functions. Further Gene co-expression analysis showed that Ace2 may be correlated with the membrane proteins in Calcium signaling pathway. Further, qRT-PCR confirmed that SARS-CoV-2 spike S1 protein upregulated ACE2 expression together with eight membrane proteins in Calcium Signaling pathway. Moreover, such elevation can be attenuated by recombinant ACE2. Collectively, our findings revealed a potential mechanism of Ace2 in cognitive dysfunction, which could be beneficial for COVID-19-induced cognitive dysfunction prevention and potential treatment.

Expression genome-wide association study identifies key regulatory variants enriched with metabolic and immune functions in four porcine tissues

BACKGROUND

Integration of high throughput DNA genotyping and RNA-sequencing data enables the discovery of genomic regions that regulate gene expression, known as expression quantitative trait loci (eQTL). In pigs, efforts to date have been mainly focused on purebred lines for traits with commercial relevance as such growth and meat quality. However, little is known on genetic variants and mechanisms associated with the robustness of an animal, thus its overall health status. Here, the liver, lung, spleen, and muscle transcriptomes of 100 three-way crossbred female finishers were studied, with the aim of identifying novel eQTL regulatory regions and transcription factors (TFs) associated with regulation of porcine metabolism and health-related traits.

RESULTS

An expression genome-wide association study with 535,896 genotypes and the expression of 12,680 genes in liver, 13,310 genes in lung, 12,650 genes in spleen, and 12,595 genes in muscle resulted in 4,293, 10,630, 4,533, and 6,871 eQTL regions for each of these tissues, respectively. Although only a small fraction of the eQTLs were annotated as cis-eQTLs, these presented a higher number of polymorphisms per region and significantly stronger associations with their target gene compared to trans-eQTLs. Between 20 and 115 eQTL hotspots were identified across the four tissues. Interestingly, these were all enriched for immune-related biological processes. In spleen, two TFs were identified: ERF and ZNF45, with key roles in regulation of gene expression.

CONCLUSIONS

This study provides a comprehensive analysis with more than 26,000 eQTL regions identified that are now publicly available. The genomic regions and their variants were mostly associated with tissue-specific regulatory roles. However, some shared regions provide new insights into the complex regulation of genes and their interactions that are involved with important traits related to metabolism and immunity.

Smooth muscle expression of RNA editing enzyme ADAR1 controls vascular integrity and progression of atherosclerosis

Mapping the genomic architecture of complex disease has been predicated on the understanding that genetic variants influence disease risk through modifying gene expression. However, recent discoveries have revealed that a significant burden of disease heritability in common autoinflammatory disorders and coronary artery disease is mediated through genetic variation modifying post-transcriptional modification of RNA through adenosine-to-inosine (A-to-I) RNA editing. This common RNA modification is catalyzed by ADAR enzymes, where ADAR1 edits specific immunogenic double stranded RNA (dsRNA) to prevent activation of the double strand RNA (dsRNA) sensor MDA5 ( IFIH1 ) and stimulation of an interferon stimulated gene (ISG) response. Multiple lines of human genetic data indicate impaired RNA editing and increased dsRNA sensing to be an important mechanism of coronary artery disease (CAD) risk. Here, we provide a crucial link between observations in human genetics and mechanistic cell biology leading to progression of CAD. Through analysis of human atherosclerotic plaque, we implicate the vascular smooth muscle cell (SMC) to have a unique requirement for RNA editing, and that ISG induction occurs in SMC phenotypic modulation, implicating MDA5 activation. Through culture of human coronary artery SMCs, generation of a conditional SMC specific Adar1 deletion mouse model on a pro-atherosclerosis background, and with incorporation of single cell RNA sequencing cellular profiling, we further show that Adar1 controls SMC phenotypic state, is required to maintain vascular integrity, and controls progression of atherosclerosis and vascular calcification. Through this work, we describe a fundamental mechanism of CAD, where cell type and context specific RNA editing and sensing of dsRNA mediates disease progression, bridging our understanding of human genetics and disease causality.

Single cell transcriptomes and multiscale networks from persons with and without Alzheimer's disease

The emergence of single nucleus RNA sequencing (snRNA-seq) offers to revolutionize the study of Alzheimer's disease (AD). Integration with complementary multiomics data such as genetics, proteomics and clinical data provides powerful opportunities to link cell subpopulations and molecular networks with a broader disease-relevant context. We report snRNA-seq profiles from superior frontal gyrus samples from 101 well characterized subjects from the Banner Brain and Body Donation Program in combination with whole genome sequences. We report findings that link common AD risk variants with CR1 expression in oligodendrocytes as well as alterations in hematological parameters. We observed an AD-associated CD83(+) microglial subtype with unique molecular networks and which is associated with immunoglobulin IgG4 production in the transverse colon. Our major observations were replicated in two additional, independent snRNA-seq data sets. These findings illustrate the power of multi-tissue molecular profiling to contextualize snRNA-seq brain transcriptomics and reveal disease biology.

Trans-eQTL hotspots shape complex traits by modulating cellular states

Regulatory genetic variation shapes gene expression, providing an important mechanism connecting DNA variation and complex traits. The causal relationships between gene expression and complex traits remain poorly understood. Here, we integrated transcriptomes and 46 genetically complex growth traits in a large cross between two strains of the yeast Saccharomyces cerevisiae. We discovered thousands of genetic correlations between gene expression and growth, suggesting potential functional connections. Local regulatory variation was a minor source of these genetic correlations. Instead, genetic correlations tended to arise from multiple independent trans-acting regulatory loci. Trans-acting hotspots that affect the expression of numerous genes accounted for particularly large fractions of genetic growth variation and of genetic correlations between gene expression and growth. Genes with genetic correlations were enriched for similar biological processes across traits, but with heterogeneous direction of effect. Our results reveal how trans-acting regulatory hotspots shape complex traits by altering cellular states.

Analysis of Evolutionary Conservation, Expression Level, and Genetic Association at a Genome-wide Scale Reveals Heterogeneity Across Polygenic Phenotypes

Understanding the expression level and evolutionary rate of associated genes with human polygenic diseases provides crucial insights into their disease-contributing roles. In this work, we leveraged genome-wide association studies (GWASs) to investigate the relationship between the genetic association and both the evolutionary rate (dN/dS) and expression level of human genes associated with the two polygenic diseases of schizophrenia and coronary artery disease. Our findings highlight a distinct variation in these relationships between the two diseases. Genes associated with both diseases exhibit a significantly greater variance in evolutionary rate compared to those implicated in monogenic diseases. Expanding our analyses to 4,756 complex traits in the GWAS atlas database, we unraveled distinct trait categories with a unique interplay among the evolutionary rate, expression level, and genetic association of human genes. In most polygenic traits, highly expressed genes were more associated with the polygenic phenotypes compared to lowly expressed genes. About 69% of polygenic traits displayed a negative correlation between genetic association and evolutionary rate, while approximately 30% of these traits showed a positive correlation between genetic association and evolutionary rate. Our results demonstrate the presence of a spectrum among complex traits, shaped by natural selection. Notably, at opposite ends of this spectrum, we find metabolic traits being more likely influenced by purifying selection, and immunological traits that are more likely shaped by positive selection. We further established the polygenic evolution portal (evopolygen.de) as a resource for investigating relationships and generating hypotheses in the field of human polygenic trait evolution.

Genomic determinants, architecture, and constraints in drought-related traits in Corymbia calophylla

BACKGROUND

Drought adaptation is critical to many tree species persisting under climate change, however our knowledge of the genetic basis for trees to adapt to drought is limited. This knowledge gap impedes our fundamental understanding of drought response and application to forest production and conservation. To improve our understanding of the genomic determinants, architecture, and trait constraints, we assembled a reference genome and detected ~ 6.5 M variants in 432 phenotyped individuals for the foundational tree Corymbia calophylla.

RESULTS

We found 273 genomic variants determining traits with moderate heritability (h2SNP = 0.26-0.64). Significant variants were predominantly in gene regulatory elements distributed among several haplotype blocks across all chromosomes. Furthermore, traits were constrained by frequent epistatic and pleiotropic interactions.

CONCLUSIONS

Our results on the genetic basis for drought traits in Corymbia calophylla have several implications for the ability to adapt to climate change: (1) drought related traits are controlled by complex genomic architectures with large haplotypes, epistatic, and pleiotropic interactions; (2) the most significant variants determining drought related traits occurred in regulatory regions; and (3) models incorporating epistatic interactions increase trait predictions. Our findings indicate that despite moderate heritability drought traits are likely constrained by complex genomic architecture potentially limiting trees response to climate change.

Identifying Novel Proteins for Chronic Pain: Integration of Human Brain Proteomes and Genome-wide Association Data

Numerous genome-wide association studies have identified risk genes for chronic pain, yet the mechanisms by which genetic variants modify susceptibility have remained elusive. We sought to identify key genes modulating chronic pain risk by regulating brain protein expression. We integrated brain proteomic data with the largest genome-wide dataset for multisite chronic pain (N = 387,649) in a proteome-wide association study (PWAS) using discovery and confirmatory proteomic datasets (N = 376 and 152) from the dorsolateral prefrontal cortex. Leveraging summary data-based Mendelian randomization and Bayesian colocalization analysis, we pinpointed potential causal genes, while a transcriptome-wide association study integrating 452 human brain transcriptomes investigated whether cis-effects on protein abundance extended to the transcriptome. Single-cell RNA-sequencing data and single-nucleus transcriptomic data revealed cell-type-specific expression patterns for identified causal genes in the dorsolateral prefrontal cortex and dorsal root ganglia (DRG), complemented by RNA microarray analysis of expression profiles in other pain-related brain regions. Of the 22 genes cis-regulating protein abundance identified by the discovery PWAS, 18 (82%) were deemed causal by summary data-based Mendelian randomization or Bayesian colocalization analysis analyses, with 7 of these 18 genes (39%) replicating in the confirmatory PWAS, including guanosine diphosphate-mannose pyrophosphorylase B, which also associated at the transcriptome level. Several causal genes exhibited selective expression in excitatory and inhibitory neurons, oligodendrocytes, and astrocytes, while most identified genes were expressed across additional pain-related brain regions. This integrative proteogenomic approach identified 18 high-confidence causal genes for chronic pain, regulated by cis-effects on brain protein levels, suggesting promising avenues for treatment research and indicating a contributory role for the DRG. PERSPECTIVE: The current post genome-wide association study analyses identified 18 high-confidence causal genes regulating chronic pain risk via cis-modulation of brain protein abundance, suggesting promising avenues for future chronic pain therapies. Additionally, the significant expression of these genes in the DRG indicated a potential contributory role, warranting further investigation.

Plasma proteome analysis implicates novel proteins as potential therapeutic targets for chronic kidney disease: A proteome-wide association study

Chronic kidney disease (CKD) is prevalent globally with limited therapeutic drugs available. To systemically identify novel proteins involved in the pathogenesis of CKD and possible therapeutic targets, we integrated human plasma proteomes with the genome-wide association studies (GWASs) of CKD, estimated glomerular filtration rate (eGFR) and blood urea nitrogen (BUN) to perform proteome-wide association study (PWAS), Mendelian Randomization and Bayesian colocalization analyses. The single-cell RNA sequencing data of healthy human and mouse kidneys were analyzed to explore the cell-type specificity of identified genes. Functional enrichment analysis was conducted to investigate the involved signaling pathways. The PWAS identified 22 plasma proteins significantly associated with CKD. Of them, the significant associations of three proteins (INHBC, LMAN2, and SNUPN) were replicated in the GWASs of eGFR, and BUN. Mendelian Randomization analyses showed that INHBC and SNUPN were causally associated with CKD, eGFR, and BUN. The Bayesian colocalization analysis identified shared causal variants for INHBC in CKD, eGFR, and BUN (all PP4 > 0.75). The single-cell RNA sequencing revealed that the INHBC gene was sparsely scattered within the kidney cells. This proteomic study revealed that INHBC, LMAN2, and SNUPN may be involved in the pathogenesis of CKD, which represent novel therapeutic targets and warrant further exploration in future research.

Machine learning integrative approaches to advance computational immunology

The study of immunology, traditionally reliant on proteomics to evaluate individual immune cells, has been revolutionized by single-cell RNA sequencing. Computational immunologists play a crucial role in analysing these datasets, moving beyond traditional protein marker identification to encompass a more detailed view of cellular phenotypes and their functional roles. Recent technological advancements allow the simultaneous measurements of multiple cellular components-transcriptome, proteome, chromatin, epigenetic modifications and metabolites-within single cells, including in spatial contexts within tissues. This has led to the generation of complex multiscale datasets that can include multimodal measurements from the same cells or a mix of paired and unpaired modalities. Modern machine learning (ML) techniques allow for the integration of multiple "omics" data without the need for extensive independent modelling of each modality. This review focuses on recent advancements in ML integrative approaches applied to immunological studies. We highlight the importance of these methods in creating a unified representation of multiscale data collections, particularly for single-cell and spatial profiling technologies. Finally, we discuss the challenges of these holistic approaches and how they will be instrumental in the development of a common coordinate framework for multiscale studies, thereby accelerating research and enabling discoveries in the computational immunology field.

Candida albicans natural diversity: a resource to dissect fungal commensalism and pathogenesis

Candida albicans is a ubiquitous fungus of humans. It is not only a component of the oral and intestinal microbiota of most healthy adults but also a major cause of mucosal disorders and life-threatening disseminated infections. Until recently, research on the biology and pathogenesis of the fungus was largely based on a single clinical isolate. We review investigations that have started to dissect a diverse set of C. albicans strains. Using different approaches to leverage the species' phenotypic and/or genetic diversity, these studies illuminate the wide range of interactions between fungus and host. While connecting genetic variants to phenotypes of interest remains challenging, research on C. albicans' natural diversity is central to understand fungal commensalism and pathogenesis.

Pan-transcriptome reveals a large accessory genome contribution to gene expression variation in yeast

Gene expression is an essential step in the translation of genotypes into phenotypes. However, little is known about the transcriptome architecture and the underlying genetic effects at the species level. Here we generated and analyzed the pan-transcriptome of ~1,000 yeast natural isolates across 4,977 core and 1,468 accessory genes. We found that the accessory genome is an underappreciated driver of transcriptome divergence. Global gene expression patterns combined with population structure showed that variation in heritable expression mainly lies within subpopulation-specific signatures, for which accessory genes are overrepresented. Genome-wide association analyses consistently highlighted that accessory genes are associated with proportionally more variants with larger effect sizes, illustrating the critical role of the accessory genome on the transcriptional landscape within and between populations.

Multiomic foundations of human prefrontal cortex tissue function

The prefrontal cortex (PFC) is a region of the brain that in humans is involved in the production of higher-order functions such as cognition, emotion, perception, and behavior. Neurotransmission in the PFC produces higher-order functions by integrating information from other areas of the brain. At the foundation of neurotransmission, and by extension at the foundation of higher-order brain functions, are an untold number of coordinated molecular processes involving the DNA sequence variants in the genome, RNA transcripts in the transcriptome, and proteins in the proteome. These "multiomic" foundations are poorly understood in humans, perhaps in part because most modern studies that characterize the molecular state of the human PFC use tissue obtained when neurotransmission and higher-order brain functions have ceased (i.e., the postmortem state). Here, analyses are presented on data generated for the Living Brain Project (LBP) to investigate whether PFC tissue from individuals with intact higher-order brain function has characteristic multiomic foundations. Two complementary strategies were employed towards this end. The first strategy was to identify in PFC samples obtained from living study participants a signature of RNA transcript expression associated with neurotransmission measured intracranially at the time of PFC sampling, in some cases while participants performed a task engaging higher-order brain functions. The second strategy was to perform multiomic comparisons between PFC samples obtained from individuals with intact higher-order brain function at the time of sampling (i.e., living study participants) and PFC samples obtained in the postmortem state. RNA transcript expression within multiple PFC cell types was associated with fluctuations of dopaminergic, serotonergic, and/or noradrenergic neurotransmission in the substantia nigra measured while participants played a computer game that engaged higher-order brain functions. A subset of these associations - termed the "transcriptional program associated with neurotransmission" (TPAWN) - were reproduced in analyses of brain RNA transcript expression and intracranial neurotransmission data obtained from a second LBP cohort and from a cohort in an independent study. RNA transcripts involved in TPAWN were found to be (1) enriched for RNA transcripts associated with measures of neurotransmission in rodent and cell models, (2) enriched for RNA transcripts encoded by evolutionarily constrained genes, (3) depleted of RNA transcripts regulated by common DNA sequence variants, and (4) enriched for RNA transcripts implicated in higher-order brain functions by human population genetic studies. In PFC excitatory neurons of living study participants, higher expression of the genes in TPAWN tracked with higher expression of RNA transcripts that in rodent PFC samples are markers of a class of excitatory neurons that connect the PFC to deep brain structures. TPAWN was further reproduced by RNA transcript expression patterns differentiating living PFC samples from postmortem PFC samples, and significant differences between living and postmortem PFC samples were additionally observed with respect to (1) the expression of most primary RNA transcripts, mature RNA transcripts, and proteins, (2) the splicing of most primary RNA transcripts into mature RNA transcripts, (3) the patterns of co-expression between RNA transcripts and proteins, and (4) the effects of some DNA sequence variants on RNA transcript and protein expression. Taken together, this report highlights that studies of brain tissue obtained in a safe and ethical manner from large cohorts of living individuals can help advance understanding of the multiomic foundations of brain function.

Species-wide quantitative transcriptomes and proteomes reveal distinct genetic control of gene expression variation in yeast

Gene expression varies between individuals and corresponds to a key step linking genotypes to phenotypes. However, our knowledge regarding the species-wide genetic control of protein abundance, including its dependency on transcript levels, is very limited. Here, we have determined quantitative proteomes of a large population of 942 diverse natural Saccharomyces cerevisiae yeast isolates. We found that mRNA and protein abundances are weakly correlated at the population gene level. While the protein coexpression network recapitulates major biological functions, differential expression patterns reveal proteomic signatures related to specific populations. Comprehensive genetic association analyses highlight that genetic variants associated with variation in protein (pQTL) and transcript (eQTL) levels poorly overlap (3%). Our results demonstrate that transcriptome and proteome are governed by distinct genetic bases, likely explained by protein turnover. It also highlights the importance of integrating these different levels of gene expression to better understand the genotype-phenotype relationship.

Evolutionary Rate Shifts in Coding and Regulatory Regions Underpin Repeated Adaptation to Sulfidic Streams in Poeciliid Fishes

Adaptation to extreme environments often involves the evolution of dramatic physiological changes. To better understand how organisms evolve these complex phenotypic changes, the repeatability and predictability of evolution, and possible constraints on adapting to an extreme environment, it is important to understand how adaptive variation has evolved. Poeciliid fishes represent a particularly fruitful study system for investigations of adaptation to extreme environments due to their repeated colonization of toxic hydrogen sulfide-rich springs across multiple species within the clade. Previous investigations have highlighted changes in the physiology and gene expression in specific species that are thought to facilitate adaptation to hydrogen sulfide-rich springs. However, the presence of adaptive nucleotide variation in coding and regulatory regions and the degree to which convergent evolution has shaped the genomic regions underpinning sulfide tolerance across taxa are unknown. By sampling across seven independent lineages in which nonsulfidic lineages have colonized and adapted to sulfide springs, we reveal signatures of shared evolutionary rate shifts across the genome. We found evidence of genes, promoters, and putative enhancer regions associated with both increased and decreased convergent evolutionary rate shifts in hydrogen sulfide-adapted lineages. Our analysis highlights convergent evolutionary rate shifts in sulfidic lineages associated with the modulation of endogenous hydrogen sulfide production and hydrogen sulfide detoxification. We also found that regions with shifted evolutionary rates in sulfide spring fishes more often exhibited convergent shifts in either the coding region or the regulatory sequence of a given gene, rather than both.

Genome-wide association study and fine-mapping on Korean biobank to discover renal trait-associated variants

BACKGROUND

Chronic kidney disease is a significant health burden worldwide, with increasing incidence. Although several genome- wide association studies (GWAS) have investigated single nucleotide polymorphisms (SNP) associated with kidney trait, most studies were focused on European ancestry.

METHODS

We utilized clinical and genetic information collected from the Korean Genome and Epidemiology Study (KoGES).

RESULTS

More than five million SNPs from 58,406 participants were analyzed. After meta-GWAS, 1,360 loci associated with estimated glomerular filtration rate (eGFR) at a genome-wide significant level (p = 5 × 10-8) were identified. Among them, 399 loci were validated with at least one other biomarker (blood urea nitrogen [BUN] or eGFRcysC) and 149 loci were validated using both markers. Among them, 18 SNPs (nine known ones and nine novel ones) with 20 putative genes were found. The aggregated effect of genes estimated by MAGMA gene analysis showed that these significant genes were enriched in kidney-associated pathways, with the kidney and liver being the most enriched tissues.

CONCLUSION

In this study, we conducted GWAS for more than 50,000 Korean individuals and identified several variants associated with kidney traits, including eGFR, BUN, and eGFRcysC. We also investigated functions of relevant genes using computational methods to define putative causal variants.

Cis-eQTLs in seven duck tissues identify novel candidate genes for growth and carcass traits

BACKGROUND

Expression quantitative trait loci (eQTL) studies aim to understand the influence of genetic variants on gene expression. The colocalization of eQTL mapping and GWAS strategy could help identify essential candidate genes and causal DNA variants vital to complex traits in human and many farm animals. However, eQTL mapping has not been conducted in ducks. It is desirable to know whether eQTLs within GWAS signals contributed to duck economic traits.

RESULTS

In this study, we conducted an eQTL analysis using publicly available RNA sequencing data from 820 samples, focusing on liver, muscle, blood, adipose, ovary, spleen, and lung tissues. We identified 113,374 cis-eQTLs for 12,266 genes, a substantial fraction 39.1% of which were discovered in at least two tissues. The cis-eQTLs of blood were less conserved across tissues, while cis-eQTLs from any tissue exhibit a strong sharing pattern to liver tissue. Colocalization between cis-eQTLs and genome-wide association studies (GWAS) of 50 traits uncovered new associations between gene expression and potential loci influencing growth and carcass traits. SRSF4, GSS, and IGF2BP1 in liver, NDUFC2 in muscle, ELF3 in adipose, and RUNDC1 in blood could serve as the candidate genes for duck growth and carcass traits.

CONCLUSIONS

Our findings highlight substantial differences in genetic regulation of gene expression across duck primary tissues, shedding light on potential mechanisms through which candidate genes may impact growth and carcass traits. Furthermore, this availability of eQTL data offers a valuable resource for deciphering further genetic association signals that may arise from ongoing extensive endeavors aimed at enhancing duck production traits.

Transcriptional Variabilities in Human hiPSC-derived Cardiomyocytes: All Genes Are Not Equal and Their Robustness May Foretell Donor's Disease Susceptibility

Human induced pluripotent stem cells (hiPSCs) are frequently used to study disease-associated variations. We characterized transcriptional variability from a hiPSC-derived cardiomyocyte (hiPSC-CM) study of left ventricular hypertrophy (LVH) using donor samples from the HyperGEN study. Multiple hiPSC-CM differentiations over reprogramming events (iPSC generation) across 7 donors were used to assess variabilities from reprogramming, differentiation, and donor LVH status. Variability arising from pathological alterations was assessed using a cardiac stimulant applied to the hiPSC-CMs to trigger hypertrophic responses. We found that for most genes (73.3%~85.5%), technical variability was smaller than biological variability. Further, we identified and characterized lists of "noise" genes showing greater technical variability and "signal" genes showing greater biological variability. Together, they support a "genetic robustness" hypothesis of disease-modeling whereby cellular response to relevant stimuli in hiPSC-derived somatic cells from diseased donors tends to show more transcriptional variability. Our findings suggest that hiPSC-CMs can provide a valid model for cardiac hypertrophy and distinguish between technical and disease-relevant transcriptional changes.

Improved multi-ancestry fine-mapping identifies cis-regulatory variants underlying molecular traits and disease risk

Multi-ancestry statistical fine-mapping of cis-molecular quantitative trait loci (cis-molQTL) aims to improve the precision of distinguishing causal cis-molQTLs from tagging variants. However, existing approaches fail to reflect shared genetic architectures. To solve this limitation, we present the Sum of Shared Single Effects (SuShiE) model, which leverages LD heterogeneity to improve fine-mapping precision, infer cross-ancestry effect size correlations, and estimate ancestry-specific expression prediction weights. We apply SuShiE to mRNA expression measured in PBMCs (n=956) and LCLs (n=814) together with plasma protein levels (n=854) from individuals of diverse ancestries in the TOPMed MESA and GENOA studies. We find SuShiE fine-maps cis-molQTLs for 16% more genes compared with baselines while prioritizing fewer variants with greater functional enrichment. SuShiE infers highly consistent cis-molQTL architectures across ancestries on average; however, we also find evidence of heterogeneity at genes with predicted loss-of-function intolerance, suggesting that environmental interactions may partially explain differences in cis-molQTL effect sizes across ancestries. Lastly, we leverage estimated cis-molQTL effect-sizes to perform individual-level TWAS and PWAS on six white blood cell-related traits in AOU Biobank individuals (n=86k), and identify 44 more genes compared with baselines, further highlighting its benefits in identifying genes relevant for complex disease risk. Overall, SuShiE provides new insights into the cis-genetic architecture of molecular traits.

Chromatin accessibility variation provides insights into missing regulation underlying immune-mediated diseases

Most genetic loci associated with complex traits and diseases through genome-wide association studies (GWAS) are noncoding, suggesting that the causal variants likely have gene regulatory effects. However, only a small number of loci have been linked to expression quantitative trait loci (eQTLs) detected currently. To better understand the potential reasons for many trait-associated loci lacking eQTL colocalization, we investigated whether chromatin accessibility QTLs (caQTLs) in lymphoblastoid cell lines (LCLs) explain immune-mediated disease associations that eQTLs in LCLs did not. The power to detect caQTLs was greater than that of eQTLs and was less affected by the distance from the transcription start site of the associated gene. Meta-analyzing LCL eQTL data to increase the sample size to over a thousand led to additional loci with eQTL colocalization, demonstrating that insufficient statistical power is still likely to be a factor. Moreover, further eQTL colocalization loci were uncovered by surveying eQTLs of other immune cell types. Altogether, insufficient power and context-specificity of eQTLs both contribute to the 'missing regulation.'

Recalibrating differential gene expression by genetic dosage variance prioritizes functionally relevant genes

Differential expression (DE) analysis is a widely used method for identifying genes that are functionally relevant for an observed phenotype or biological response. However, typical DE analysis includes selection of genes based on a threshold of fold change in expression under the implicit assumption that all genes are equally sensitive to dosage changes of their transcripts. This tends to favor highly variable genes over more constrained genes where even small changes in expression may be biologically relevant. To address this limitation, we have developed a method to recalibrate each gene's differential expression fold change based on genetic expression variance observed in the human population. The newly established metric ranks statistically differentially expressed genes not by nominal change of expression, but by relative change in comparison to natural dosage variation for each gene. We apply our method to RNA sequencing datasets from rare disease and in-vitro stimulus response experiments. Compared to the standard approach, our method adjusts the bias in discovery towards highly variable genes, and enriches for pathways and biological processes related to metabolic and regulatory activity, indicating a prioritization of functionally relevant driver genes. With that, our method provides a novel view on DE and contributes towards bridging the existing gap between statistical and biological significance. We believe that this approach will simplify the identification of disease causing genes and enhance the discovery of therapeutic targets.

The Impact of Inherited Genetic Variation on DNA Methylation in Prostate Cancer and Benign Tissues of African American and European American Men

BACKGROUND

American men of African ancestry (AA) have higher prostate cancer incidence and mortality rates compared with American men of European ancestry (EA). Differences in genetic susceptibility mechanisms may contribute to this disparity.

METHODS

To gain insights into the regulatory mechanisms of prostate cancer susceptibility variants, we tested the association between SNPs and DNA methylation (DNAm) at nearby CpG sites across the genome in benign and cancer prostate tissue from 74 AA and 74 EA men. Genome-wide SNP data (from benign tissue) and DNAm were generated using Illumina arrays.

RESULTS

Among AA men, we identified 6,298 and 2,641 cis-methylation QTLs (meQTL; FDR of 0.05) in benign and tumor tissue, respectively, with 6,960 and 1,700 detected in EA men. We leveraged genome-wide association study (GWAS) summary statistics to identify previously reported prostate cancer GWAS signals likely to share a common causal variant with a detected meQTL. We identified nine GWAS-meQTL pairs with strong evidence of colocalization (four in EA benign, three in EA tumor, two in AA benign, and three in AA tumor). Among these colocalized GWAS-meQTL pairs, we identified colocalizing expression quantitative trait loci (eQTL) impacting four eGenes with known roles in tumorigenesis.

CONCLUSIONS

These findings highlight epigenetic regulatory mechanisms by which prostate cancer-risk SNPs can modify local DNAm and/or gene expression in prostate tissue.

IMPACT

Overall, our findings showed general consistency in the meQTL landscape of AA and EA men, but meQTLs often differ by tissue type (normal vs. cancer). Ancestry-based linkage disequilibrium differences and lack of AA representation in GWAS decrease statistical power to detect colocalization for some regions.

Diallel panel reveals a significant impact of low-frequency genetic variants on gene expression variation in yeast

Unraveling the genetic sources of gene expression variation is essential to better understand the origins of phenotypic diversity in natural populations. Genome-wide association studies identified thousands of variants involved in gene expression variation, however, variants detected only explain part of the heritability. In fact, variants such as low-frequency and structural variants (SVs) are poorly captured in association studies. To assess the impact of these variants on gene expression variation, we explored a half-diallel panel composed of 323 hybrids originated from pairwise crosses of 26 natural Saccharomyces cerevisiae isolates. Using short- and long-read sequencing strategies, we established an exhaustive catalog of single nucleotide polymorphisms (SNPs) and SVs for this panel. Combining this dataset with the transcriptomes of all hybrids, we comprehensively mapped SNPs and SVs associated with gene expression variation. While SVs impact gene expression variation, SNPs exhibit a higher effect size with an overrepresentation of low-frequency variants compared to common ones. These results reinforce the importance of dissecting the heritability of complex traits with a comprehensive catalog of genetic variants at the population level.

Identification of candidate causal cis -regulatory variants underlying electrocardiographic QT interval GWAS loci

Identifying causal variants among tens or hundreds of associated variants at each locus mapped by genome-wide association studies (GWAS) of complex traits is a challenge. As vast majority of GWAS variants are noncoding, sequence variation at cis -regulatory elements affecting transcriptional expression of specific genes is a widely accepted molecular hypothesis. Following this cis -regulatory hypothesis and combining it with the observation that nucleosome-free open chromatin is a universal hallmark of all types of cis -regulatory elements, we aimed to identify candidate causal regulatory variants underlying electrocardiographic QT interval GWAS loci. At a dozen loci, selected for higher effect sizes and a better understanding of the likely causal gene, we identified and included all common variants in high linkage disequilibrium with the GWAS variants as candidate variants. Using ENCODE DNase-seq and ATAC-seq from multiple human adult cardiac left ventricle tissue samples, we generated genome-wide maps of open chromatin regions marking putative regulatory elements. QT interval associated candidate variants were filtered for overlap with cardiac left ventricle open chromatin regions to identify candidate causal cis -regulatory variants, which were further assessed for colocalizing with a known cardiac GTEx expression quantitative trait locus variant as additional evidence for their causal role. Together, these efforts have generated a comprehensive set of candidate causal variants that are expected to be enriched for cis -regulatory potential and thereby, explaining the observed genetic associations.

Linking environmental factors and gene regulation

A technique called mSTARR-seq sheds light on how DNA methylation may shape responses to external stimuli by altering the activity of sequences that control gene expression.

Rs4862705 in the melatonin receptor 1A gene is associated with renal function decline in type 1 diabetes individuals

Aim

The pathogenesis of chronic diabetes complications has oxidative stress as one of the major elements, and single-nucleotide polymorphisms (SNPs) in genes belonging to antioxidant pathways modulate susceptibility to these complications. Considering that melatonin is a powerful antioxidant compound, our aim was to explore, in a longitudinal cohort study of type 1 diabetes (T1D) individuals, the association of microvascular complications and SNPs in the gene encoding melatonin receptor 1A (MTNR1A).

Methods

Eight SNPs in MTNR1A were genotyped in 489 T1D individuals. Besides cross-sectional analyses of SNPs with each one of the microvascular complications (distal polyneuropathy, cardiovascular autonomic neuropathy, retinopathy, and diabetic kidney disease), a longitudinal analysis evaluated the associations of SNPs with renal function decline in 411 individuals followed up for a median of 8 years. In a subgroup of participants, the association of complications with urinary 6-sulfatoxymelatonin (aMT6s) concentration was investigated.

Results

The group of individuals with a renal function decline ≥ 5 mL min-1 1.73 m-2 year-1 presented a higher frequency of the A allele of rs4862705 in comparison with nondecliners, even after adjustment for confounding variables (OR = 1.84, 95% CI = 1.20-2.82; p = 0.0046). No other significant associations were found.

Conclusions

This is the first study showing an association between a variant in a gene belonging to the melatonin system and renal function decline in the diabetic setting.

Genomic dissection of the correlation between milk yield and various health traits using functional and evolutionary information about imputed sequence variants of 34,497 German Holstein cows

BACKGROUND

Over the last decades, it was subject of many studies to investigate the genomic connection of milk production and health traits in dairy cattle. Thereby, incorporating functional information in genomic analyses has been shown to improve the understanding of biological and molecular mechanisms shaping complex traits and the accuracies of genomic prediction, especially in small populations and across-breed settings. Still, little is known about the contribution of different functional and evolutionary genome partitioning subsets to milk production and dairy health. Thus, we performed a uni- and a bivariate analysis of milk yield (MY) and eight health traits using a set of ~34,497 German Holstein cows with 50K chip genotypes and ~17 million imputed sequence variants divided into 27 subsets depending on their functional and evolutionary annotation. In the bivariate analysis, eight trait-combinations were observed that contrasted MY with each health trait. Two genomic relationship matrices (GRM) were included, one consisting of the 50K chip variants and one consisting of each set of subset variants, to obtain subset heritabilities and genetic correlations. In addition, 50K chip heritabilities and genetic correlations were estimated applying merely the 50K GRM.

RESULTS

In general, 50K chip heritabilities were larger than the subset heritabilities. The largest heritabilities were found for MY, which was 0.4358 for the 50K and 0.2757 for the subset heritabilities. Whereas all 50K genetic correlations were negative, subset genetic correlations were both, positive and negative (ranging from -0.9324 between MY and mastitis to 0.6662 between MY and digital dermatitis). The subsets containing variants which were annotated as noncoding related, splice sites, untranslated regions, metabolic quantitative trait loci, and young variants ranked highest in terms of their contribution to the traits` genetic variance. We were able to show that linkage disequilibrium between subset variants and adjacent variants did not cause these subsets` high effect.

CONCLUSION

Our results confirm the connection of milk production and health traits in dairy cattle via the animals` metabolic state. In addition, they highlight the potential of including functional information in genomic analyses, which helps to dissect the extent and direction of the observed traits` connection in more detail.

The relationship between extreme inter-individual variation in macrophage gene expression and genetic susceptibility to inflammatory bowel disease

The differentiation of resident intestinal macrophages from blood monocytes depends upon signals from the macrophage colony-stimulating factor receptor (CSF1R). Analysis of genome-wide association studies (GWAS) indicates that dysregulation of macrophage differentiation and response to microorganisms contributes to susceptibility to chronic inflammatory bowel disease (IBD). Here, we analyzed transcriptomic variation in monocyte-derived macrophages (MDM) from affected and unaffected sib pairs/trios from 22 IBD families and 6 healthy controls. Transcriptional network analysis of the data revealed no overall or inter-sib distinction between affected and unaffected individuals in basal gene expression or the temporal response to lipopolysaccharide (LPS). However, the basal or LPS-inducible expression of individual genes varied independently by as much as 100-fold between subjects. Extreme independent variation in the expression of pairs of HLA-associated transcripts (HLA-B/C, HLA-A/F and HLA-DRB1/DRB5) in macrophages was associated with HLA genotype. Correlation analysis indicated the downstream impacts of variation in the immediate early response to LPS. For example, variation in early expression of IL1B was significantly associated with local SNV genotype and with subsequent peak expression of target genes including IL23A, CXCL1, CXCL3, CXCL8 and NLRP3. Similarly, variation in early IFNB1 expression was correlated with subsequent expression of IFN target genes. Our results support the view that gene-specific dysregulation in macrophage adaptation to the intestinal milieu is associated with genetic susceptibility to IBD.

Genetic and molecular architecture of complex traits

Human genetics has emerged as one of the most dynamic areas of biology, with a broadening societal impact. In this review, we discuss recent achievements, ongoing efforts, and future challenges in the field. Advances in technology, statistical methods, and the growing scale of research efforts have all provided many insights into the processes that have given rise to the current patterns of genetic variation. Vast maps of genetic associations with human traits and diseases have allowed characterization of their genetic architecture. Finally, studies of molecular and cellular effects of genetic variants have provided insights into biological processes underlying disease. Many outstanding questions remain, but the field is well poised for groundbreaking discoveries as it increases the use of genetic data to understand both the history of our species and its applications to improve human health.

A genetic variant in the immune-related gene ERAP1 affects colorectal cancer prognosis

BACKGROUND

Findings on the association of genetic factors and colorectal cancer (CRC) survival are limited and inconsistent, and revealing the mechanism underlying their prognostic roles is of great importance. This study aimed to explore the relationship between functional genetic variations and the prognosis of CRC and further reveal the possible mechanism.

METHODS

We first systematically performed expression quantitative trait locus (eQTL) analysis using The Cancer Genome Atlas (TCGA) dataset. Then, the Kaplan-Meier analysis was used to filter out the survival-related eQTL target genes of CRC patients in two public datasets (TCGA and GSE39582 dataset from the Gene Expression Omnibus database). The seven most potentially functional eQTL single nucleotide polymorphisms (SNPs) associated with six survival-related eQTL target genes were genotyped in 907 Chinese CRC patients with clinical prognosis data. The regulatory mechanism of the survival-related SNP was further confirmed by functional experiments.

RESULTS

The rs71630754 regulating the expression of endoplasmic reticulum aminopeptidase 1 ( ERAP1 ) was significantly associated with the prognosis of CRC (additive model, hazard ratio [HR]: 1.43, 95% confidence interval [CI]: 1.08-1.88, P = 0.012). The results of dual-luciferase reporter assay and electrophoretic mobility shift assay showed that the A allele of the rs71630754 could increase the binding of transcription factor 3 (TCF3) and subsequently reduce the expression of ERAP1 . The results of bioinformatic analysis showed that lower expression of ERAP1 could affect the tumor immune microenvironment and was significantly associated with severe survival outcomes.

CONCLUSION

The rs71630754 could influence the prognosis of CRC patients by regulating the expression of the immune-related gene ERAP1 .

TRIAL REGISTRATION

No. NCT00454519 ( https://clinicaltrials.gov/ ).

A review and analysis of key biomarkers in Alzheimer's disease

Alzheimer's disease (AD) is a progressive neurodegenerative disorder that affects over 50 million elderly individuals worldwide. Although the pathogenesis of AD is not fully understood, based on current research, researchers are able to identify potential biomarker genes and proteins that may serve as effective targets against AD. This article aims to present a comprehensive overview of recent advances in AD biomarker identification, with highlights on the use of various algorithms, the exploration of relevant biological processes, and the investigation of shared biomarkers with co-occurring diseases. Additionally, this article includes a statistical analysis of key genes reported in the research literature, and identifies the intersection with AD-related gene sets from databases such as AlzGen, GeneCard, and DisGeNet. For these gene sets, besides enrichment analysis, protein-protein interaction (PPI) networks utilized to identify central genes among the overlapping genes. Enrichment analysis, protein interaction network analysis, and tissue-specific connectedness analysis based on GTEx database performed on multiple groups of overlapping genes. Our work has laid the foundation for a better understanding of the molecular mechanisms of AD and more accurate identification of key AD markers.

Mapping Expression Quantitative Trait Loci Targeting Candidate Genes for Pregnancy in Beef Cows

Despite collective efforts to understand the complex regulation of reproductive traits, no causative genes and/or mutations have been reported yet. By integrating genomics and transcriptomics data, potential regulatory mechanisms may be unveiled, providing opportunities to dissect the genetic factors governing fertility. Herein, we identified regulatory variants from RNA-Seq data associated with gene expression regulation in the uterine luminal epithelial cells of beef cows. We identified 4676 cis and 7682 trans eQTLs (expression quantitative trait loci) affecting the expression of 1120 and 2503 genes, respectively (FDR < 0.05). These variants affected the expression of transcription factor coding genes (71 cis and 193 trans eQTLs) and genes previously reported as differentially expressed between pregnant and nonpregnant cows. Functional over-representation analysis highlighted pathways related to metabolism, immune response, and hormone signaling (estrogen and GnRH) affected by eQTL-regulated genes (p-value ≤ 0.01). Furthermore, eQTLs were enriched in QTL regions for 13 reproduction-related traits from the CattleQTLdb (FDR ≤ 0.05). Our study provides novel insights into the genetic basis of reproductive processes in cattle. The underlying causal mechanisms modulating the expression of uterine genes warrant further investigation.

Comparative analysis of models in predicting the effects of SNPs on TF-DNA binding using large-scale in vitro and in vivo data

Non-coding variants associated with complex traits can alter the motifs of transcription factor (TF)-deoxyribonucleic acid binding. Although many computational models have been developed to predict the effects of non-coding variants on TF binding, their predictive power lacks systematic evaluation. Here we have evaluated 14 different models built on position weight matrices (PWMs), support vector machines, ordinary least squares and deep neural networks (DNNs), using large-scale in vitro (i.e. SNP-SELEX) and in vivo (i.e. allele-specific binding, ASB) TF binding data. Our results show that the accuracy of each model in predicting SNP effects in vitro significantly exceeds that achieved in vivo. For in vitro variant impact prediction, kmer/gkm-based machine learning methods (deltaSVM_HT-SELEX, QBiC-Pred) trained on in vitro datasets exhibit the best performance. For in vivo ASB variant prediction, DNN-based multitask models (DeepSEA, Sei, Enformer) trained on the ChIP-seq dataset exhibit relatively superior performance. Among the PWM-based methods, tRap demonstrates better performance in both in vitro and in vivo evaluations. In addition, we find that TF classes such as basic leucine zipper factors could be predicted more accurately, whereas those such as C2H2 zinc finger factors are predicted less accurately, aligning with the evolutionary conservation of these TF classes. We also underscore the significance of non-sequence factors such as cis-regulatory element type, TF expression, interactions and post-translational modifications in influencing the in vivo predictive performance of TFs. Our research provides valuable insights into selecting prioritization methods for non-coding variants and further optimizing such models.

Multi-omics and pathway analyses of genome-wide associations implicate regulation and immunity in verbal declarative memory performance

BACKGROUND

Uncovering the functional relevance underlying verbal declarative memory (VDM) genome-wide association study (GWAS) results may facilitate the development of interventions to reduce age-related memory decline and dementia.

METHODS

We performed multi-omics and pathway enrichment analyses of paragraph (PAR-dr) and word list (WL-dr) delayed recall GWAS from 29,076 older non-demented individuals of European descent. We assessed the relationship between single-variant associations and expression quantitative trait loci (eQTLs) in 44 tissues and methylation quantitative trait loci (meQTLs) in the hippocampus. We determined the relationship between gene associations and transcript levels in 53 tissues, annotation as immune genes, and regulation by transcription factors (TFs) and microRNAs. To identify significant pathways, gene set enrichment was tested in each cohort and meta-analyzed across cohorts. Analyses of differential expression in brain tissues were conducted for pathway component genes.

RESULTS

The single-variant associations of VDM showed significant linkage disequilibrium (LD) with eQTLs across all tissues and meQTLs within the hippocampus. Stronger WL-dr gene associations correlated with reduced expression in four brain tissues, including the hippocampus. More robust PAR-dr and/or WL-dr gene associations were intricately linked with immunity and were influenced by 31 TFs and 2 microRNAs. Six pathways, including type I diabetes, exhibited significant associations with both PAR-dr and WL-dr. These pathways included fifteen MHC genes intricately linked to VDM performance, showing diverse expression patterns based on cognitive status in brain tissues.

CONCLUSIONS

VDM genetic associations influence expression regulation via eQTLs and meQTLs. The involvement of TFs, microRNAs, MHC genes, and immune-related pathways contributes to VDM performance in older individuals.

Haplotype-aware modeling of cis-regulatory effects highlights the gaps remaining in eQTL data

Expression Quantitative Trait Loci (eQTLs) are critical to understanding the mechanisms underlying disease-associated genomic loci. Nearly all protein-coding genes in the human genome have been associated with one or more eQTLs. Here we introduce a multi-variant generalization of allelic Fold Change (aFC), aFC-n, to enable quantification of the cis-regulatory effects in multi-eQTL genes under the assumption that all eQTLs are known and conditionally independent. Applying aFC-n to 458,465 eQTLs in the Genotype-Tissue Expression (GTEx) project data, we demonstrate significant improvements in accuracy over the original model in estimating the eQTL effect sizes and in predicting genetically regulated gene expression over the current tools. We characterize some of the empirical properties of the eQTL data and use this framework to assess the current state of eQTL data in terms of characterizing cis-regulatory landscape in individual genomes. Notably, we show that 77.4% of the genes with an allelic imbalance in a sample show 0.5 log2 fold or more of residual imbalance after accounting for the eQTL data underlining the remaining gap in characterizing regulatory landscape in individual genomes. We further contrast this gap across tissue types, and ancestry backgrounds to identify its correlates and guide future studies.

Non-additive genetic components contribute significantly to population-wide gene expression variation

Gene expression variation, an essential step between genotype and phenotype, is collectively controlled by local (cis) and distant (trans) regulatory changes. Nevertheless, how these regulatory elements differentially influence gene expression variation remains unclear. Here, we bridge this gap by analyzing the transcriptomes of a large diallel panel consisting of 323 unique hybrids originating from genetically divergent Saccharomyces cerevisiae isolates. Our analysis across 5,087 transcript abundance traits showed that non-additive components account for 36% of the gene expression variance on average. By comparing allele-specific read counts in parent-hybrid trios, we found that trans-regulatory changes underlie the majority of gene expression variation in the population. Remarkably, most cis-regulatory variations are also exaggerated or attenuated by additional trans effects. Overall, we showed that the transcriptome is globally buffered at the genetic level mainly due to trans-regulatory variation in the population.

Trans-regulatory variant network contributes to missing heritability

In a recent Cell Genomics article, Tsouris et al.1 analyze the transcriptomes of a large diallel panel of hybrids from Saccharomyces cerevisiae natural isolates to study cis- and trans-regulatory changes underlying gene expression variation. Vanessa Pereira and Elena Kuzmin discuss the authors' findings and the wider context in missing heritability research in this preview.

PigBiobank: a valuable resource for understanding genetic and biological mechanisms of diverse complex traits in pigs

To fully unlock the potential of pigs as both agricultural species for animal-based protein food and biomedical models for human biology and disease, a comprehensive understanding of molecular and cellular mechanisms underlying various complex phenotypes in pigs and how the findings can be translated to other species, especially humans, are urgently needed. Here, within the Farm animal Genotype-Tissue Expression (FarmGTEx) project, we build the PigBiobank (http://pigbiobank.farmgtex.org) to systematically investigate the relationships among genomic variants, regulatory elements, genes, molecular networks, tissues and complex traits in pigs. This first version of the PigBiobank curates 71 885 pigs with both genotypes and phenotypes from over 100 pig breeds worldwide, covering 264 distinct complex traits. The PigBiobank has the following functions: (i) imputed sequence-based genotype-phenotype associations via a standardized and uniform pipeline, (ii) molecular and cellular mechanisms underlying trait-associations via integrating multi-omics data, (iii) cross-species gene mapping of complex traits via transcriptome-wide association studies, and (iv) high-quality results display and visualization. The PigBiobank will be updated timely with the development of the FarmGTEx-PigGTEx project, serving as an open-access and easy-to-use resource for genetically and biologically dissecting complex traits in pigs and translating the findings to other species.

PancanQTLv2.0: a comprehensive resource for expression quantitative trait loci across human cancers

Expression quantitative trait locus (eQTL) analysis is a powerful tool used to investigate genetic variations in complex diseases, including cancer. We previously developed a comprehensive database, PancanQTL, to characterize cancer eQTLs using The Cancer Genome Atlas (TCGA) dataset, and linked eQTLs with patient survival and GWAS risk variants. Here, we present an updated version, PancanQTLv2.0 (https://hanlaboratory.com/PancanQTLv2/), with advancements in fine-mapping causal variants for eQTLs, updating eQTLs overlapping with GWAS linkage disequilibrium regions and identifying eQTLs associated with drug response and immune infiltration. Through fine-mapping analysis, we identified 58 747 fine-mapped eQTLs credible sets, providing mechanic insights of gene regulation in cancer. We further integrated the latest GWAS Catalog and identified a total of 84 592 135 linkage associations between eQTLs and the existing GWAS loci, which represents a remarkable ∼50-fold increase compared to the previous version. Additionally, PancanQTLv2.0 uncovered 659516 associations between eQTLs and drug response and identified 146948 associations between eQTLs and immune cell abundance, providing potentially clinical utility of eQTLs in cancer therapy. PancanQTLv2.0 expanded the resources available for investigating gene expression regulation in human cancers, leading to advancements in cancer research and precision oncology.

Nanomaterials for detection of biomolecules and delivering therapeutic agents in theragnosis: A review

Nanomaterials are emerging facts used to deliver therapeutic agents in living systems. Nanotechnology is used as a compliment by implementing different kinds of nanotechnological applications such as nano-porous structures, functionalized nanomaterials, quantum dots, carbon nanomaterials, and polymeric nanostructures. The applications are in the initial stage, which led to achieving several diagnoses and therapy in clinical practice. This review conveys the importance of nanomaterials in post-genomic employment, which includes the design of immunosensors, immune assays, and drug delivery. In this view, genomics is a molecular tool containing large databases that are useful in choosing an apt molecular inhibitor such as drug, ligand and antibody target in the drug delivery process. This study identifies the expression of genes and proteins in analysis and classification of diseases. Experimentally, the study analyses the design of a disease model. In particular, drug delivery is a boon area to treat cancer. The identified drugs enter different phase trails (Trails I, II, and III). The genomic information conveys more essential entities to the phase I trials and helps to move further for other trails such as trails-II and III. In such cases, the biomarkers play a crucial role by monitoring the unique pathological process. Genetic engineering with recombinant DNA techniques can be employed to develop genetically engineered disease models. Delivering drugs in a specific area is one of the challenging issues achieved using nanoparticles. Therefore, genomics is considered as a vast molecular tool to identify drugs in personalized medicine for cancer therapy.

A compendium of genetic regulatory effects across pig tissues

The Farm Animal Genotype-Tissue Expression (FarmGTEx) project has been established to develop a public resource of genetic regulatory variants in livestock, which is essential for linking genetic polymorphisms to variation in phenotypes, helping fundamental biological discovery and exploitation in animal breeding and human biomedicine. Here we show results from the pilot phase of PigGTEx by processing 5,457 RNA-sequencing and 1,602 whole-genome sequencing samples passing quality control from pigs. We build a pig genotype imputation panel and associate millions of genetic variants with five types of transcriptomic phenotypes in 34 tissues. We evaluate tissue specificity of regulatory effects and elucidate molecular mechanisms of their action using multi-omics data. Leveraging this resource, we decipher regulatory mechanisms underlying 207 pig complex phenotypes and demonstrate the similarity of pigs to humans in gene expression and the genetic regulation behind complex phenotypes, supporting the importance of pigs as a human biomedical model.

Identification of constrained sequence elements across 239 primate genomes

Noncoding DNA is central to our understanding of human gene regulation and complex diseases1,2, and measuring the evolutionary sequence constraint can establish the functional relevance of putative regulatory elements in the human genome3-9. Identifying the genomic elements that have become constrained specifically in primates has been hampered by the faster evolution of noncoding DNA compared to protein-coding DNA10, the relatively short timescales separating primate species11, and the previously limited availability of whole-genome sequences12. Here we construct a whole-genome alignment of 239 species, representing nearly half of all extant species in the primate order. Using this resource, we identified human regulatory elements that are under selective constraint across primates and other mammals at a 5% false discovery rate. We detected 111,318 DNase I hypersensitivity sites and 267,410 transcription factor binding sites that are constrained specifically in primates but not across other placental mammals and validate their cis-regulatory effects on gene expression. These regulatory elements are enriched for human genetic variants that affect gene expression and complex traits and diseases. Our results highlight the important role of recent evolution in regulatory sequence elements differentiating primates, including humans, from other placental mammals.

An Update on the Genetics of IgA Nephropathy

Immunoglobulin A (IgA) nephropathy (IgAN), the most common form of glomerulonephritis, is one of the leading causes of end-stage kidney disease (ESKD). It is widely believed that genetic factors play a significant role in the development of IgAN. Previous studies of IgAN have provided important insights to unravel the genetic architecture of IgAN and its potential pathogenic mechanisms. The genome-wide association studies (GWASs) together have identified over 30 risk loci for IgAN, which emphasizes the importance of IgA production and regulation in the pathogenesis of IgAN. Follow-up fine-mapping studies help to elucidate the candidate causal variant and the potential pathogenic molecular pathway and provide new potential therapeutic targets. With the rapid development of next-generation sequencing technologies, linkage studies based on whole-genome sequencing (WGS)/whole-exome sequencing (WES) also identify rare variants associated with IgAN, accounting for some of the missing heritability. The complexity of pathogenesis and phenotypic variability may be better understood by integrating genetics, epigenetics, and environment. We have compiled a review summarizing the latest advancements in genetic studies on IgAN. We similarly summarized relevant studies examining the involvement of epigenetics in the pathogenesis of IgAN. Future directions and challenges in this field are also proposed.

Single-cell eQTL mapping in yeast reveals a tradeoff between growth and reproduction

Expression quantitative trait loci (eQTLs) provide a key bridge between noncoding DNA sequence variants and organismal traits. The effects of eQTLs can differ among tissues, cell types, and cellular states, but these differences are obscured by gene expression measurements in bulk populations. We developed a one-pot approach to map eQTLs in Saccharomyces cerevisiae by single-cell RNA sequencing (scRNA-seq) and applied it to over 100,000 single cells from three crosses. We used scRNA-seq data to genotype each cell, measure gene expression, and classify the cells by cell-cycle stage. We mapped thousands of local and distant eQTLs and identified interactions between eQTL effects and cell-cycle stages. We took advantage of single-cell expression information to identify hundreds of genes with allele-specific effects on expression noise. We used cell-cycle stage classification to map 20 loci that influence cell-cycle progression. One of these loci influenced the expression of genes involved in the mating response. We showed that the effects of this locus arise from a common variant (W82R) in the gene GPA1, which encodes a signaling protein that negatively regulates the mating pathway. The 82R allele increases mating efficiency at the cost of slower cell-cycle progression and is associated with a higher rate of outcrossing in nature. Our results provide a more granular picture of the effects of genetic variants on gene expression and downstream traits.

Chargaff's second parity rule lies at the origin of additive genetic interactions in quantitative traits to make omnigenic selection possible

Background

Francis Crick's central dogma provides a residue-by-residue mechanistic explanation of the flow of genetic information in living systems. However, this principle may not be sufficient for explaining how random mutations cause continuous variation of quantitative highly polygenic complex traits. Chargaff's second parity rule (CSPR), also referred to as intrastrand DNA symmetry, defined as near-exact equalities G ≈ C and A ≈ T within a single DNA strand, is a statistical property of cellular genomes. The phenomenon of intrastrand DNA symmetry was discovered more than 50 years ago; at present, it remains unclear what its biological role is, what the mechanisms are that force cellular genomes to comply strictly with CSPR, and why genomes of certain noncellular organisms have broken intrastrand DNA symmetry. The present work is aimed at studying a possible link between intrastrand DNA symmetry and the origin of genetic interactions in quantitative traits.

Methods

Computational analysis of single-nucleotide polymorphisms in human and mouse populations and of nucleotide composition biases at different codon positions in bacterial and human proteomes.

Results

The analysis of mutation spectra inferred from single-nucleotide polymorphisms observed in murine and human populations revealed near-exact equalities of numbers of reverse complementary mutations, indicating that random genetic variations obey CSPR. Furthermore, nucleotide compositions of coding sequences proved to be statistically interwoven via CSPR because pyrimidine bias at the 3rd codon position compensates purine bias at the 1st and 2nd positions.

Conclusions

According to Fisher's infinitesimal model, we propose that accumulation of reverse complementary mutations results in a continuous phenotypic variation due to small additive effects of statistically interwoven genetic variations. Therefore, additive genetic interactions can be inferred as a statistical entanglement of nucleotide compositions of separate genetic loci. CSPR challenges the neutral theory of molecular evolution-because all random mutations participate in variation of a trait-and provides an alternative solution to Haldane's dilemma by making a gene function diffuse. We propose that CSPR is symmetry of Fisher's infinitesimal model and that genetic information can be transferred in an implicit contactless manner.

Single-cell methylation analysis of brain tissue prioritizes mutations that alter transcription

Relating genetic variants to behavior remains a fundamental challenge. To assess the utility of DNA methylation marks in discovering causative variants, we examined their relationship to genetic variation by generating single-nucleus methylomes from the hippocampus of eight inbred mouse strains. At CpG sequence densities under 40 CpG/Kb, cells compensate for loss of methylated sites by methylating additional sites to maintain methylation levels. At higher CpG sequence densities, the exact location of a methylated site becomes more important, suggesting that variants affecting methylation will have a greater effect when occurring in higher CpG densities than in lower. We found this to be true for a variant's effect on transcript abundance, indicating that candidate variants can be prioritized based on CpG sequence density. Our findings imply that DNA methylation influences the likelihood that mutations occur at specific sites in the genome, supporting the view that the distribution of mutations is not random.

DNA methylation modulated genetic variant effect on gene transcriptional regulation

BACKGROUND

Expression quantitative trait locus (eQTL) analysis has emerged as an important tool in elucidating the link between genetic variants and gene expression, thereby bridging the gap between risk SNPs and associated diseases. We recently identified and validated a specific case where the methylation of a CpG site influences the relationship between the genetic variant and gene expression.

RESULTS

Here, to systematically evaluate this regulatory mechanism, we develop an extended eQTL mapping method, termed DNA methylation modulated eQTL (memo-eQTL). Applying this memo-eQTL mapping method to 128 normal prostate samples enables identification of 1063 memo-eQTLs, the majority of which are not recognized as conventional eQTLs in the same cohort. We observe that the methylation of the memo-eQTL CpG sites can either enhance or insulate the interaction between SNP and gene expression by altering CTCF-based chromatin 3D structure.

CONCLUSIONS

This study demonstrates the prevalence of memo-eQTLs paving the way to identify novel causal genes for traits or diseases associated with genetic variations.