A Multi-Omics Perspective of Quantitative Trait Loci in Precision Medicine

Prospective of Mitochondrial DNA Variations in Cancer on Genomic Medicine

Mitochondrial DNA (mtDNA) has emerged as a pivotal component in understanding the etiology and susceptibility of cancer. A recent study by Chen and colleagues delineated the germline genetic effect of mtDNA single-nucleotide polymorphisms (SNP) and haplogroups across pan-cancer risk. They identified a subset of mtSNPs and the corresponding risk score, as well as haplogroups A and M7 alongside their genetic interactions, conferring a protective effect against various cancers. These findings underscored the value of mtDNA variations as biomarkers for cancer etiology and as tools for cancer risk stratification. Future investigations are encouraged to integrate comprehensive omics data of genomics, transcriptomics, proteomics, and metabolomics, etc., from nuclear DNA with mtDNA variations, alongside single-cell and spatial technologies, to unravel the tumor mechanism and identify the drug targets. Moreover, the incorporation of polygenic risk score, that included mtDNA variations with both rare and common frequencies, and liquid biopsy-based biomarkers would enhance the predictive performance of cancer risk assessment and refine the risk stratification of population-based cancer screening. This commentary advocates for the validation across diverse populations to harness the full potential of mitochondrial genomics, and ultimately paves the prospective way for advancements in personalized cancer therapeutics and prevention strategies. See related article by Chen and colleagues, Cancer Epidemiol Biomarkers Prev 2024;33:381-8.

Review of Personalized Medicine and Pharmacogenomics of Anti-Cancer Compounds and Natural Products

In recent years, the FDA has approved numerous anti-cancer drugs that are mutation-based for clinical use. These drugs have improved the precision of treatment and reduced adverse effects and side effects. Personalized therapy is a prominent and hot topic of current medicine and also represents the future direction of development. With the continuous advancements in gene sequencing and high-throughput screening, research and development strategies for personalized clinical drugs have developed rapidly. This review elaborates the recent personalized treatment strategies, which include artificial intelligence, multi-omics analysis, chemical proteomics, and computation-aided drug design. These technologies rely on the molecular classification of diseases, the global signaling network within organisms, and new models for all targets, which significantly support the development of personalized medicine. Meanwhile, we summarize chemical drugs, such as lorlatinib, osimertinib, and other natural products, that deliver personalized therapeutic effects based on genetic mutations. This review also highlights potential challenges in interpreting genetic mutations and combining drugs, while providing new ideas for the development of personalized medicine and pharmacogenomics in cancer study.

Combined SNPs sequencing and allele specific proteomics capture reveal functional causality underpinning the 2p25 prostate cancer susceptibility locus

Genome wide association studies (GWASs) have identified numerous risk loci associated with prostate cancer, yet unraveling their functional significance remains elusive. Leveraging our high-throughput SNPs-seq method, we pinpointed rs4519489 within the multi-ancestry GWAS-discovered 2p25 locus as a potential functional SNP due to its significant allelic differences in protein binding. Here, we conduct a comprehensive analysis of rs4519489 and its associated gene, NOL10, employing diverse cohort data and experimental models. Clinical findings reveal a synergistic effect between rs4519489 genotype and NOL10 expression on prostate cancer prognosis and severity. Through unbiased proteomics screening, we reveal that the risk allele A of rs4519489 exhibits enhanced binding to USF1, a novel oncogenic transcription factor (TF) implicated in prostate cancer progression and prognosis, resulting in elevated NOL10 expression. Furthermore, we elucidate that NOL10 regulates cell cycle pathways, fostering prostate cancer progression. The concurrent expression of NOL10 and USF1 correlates with aggressive prostate cancer characteristics and poorer prognosis. Collectively, our study offers a robust strategy for functional SNP screening and TF identification through high-throughput SNPs-seq and unbiased proteomics, highlighting the rs4519489-USF1-NOL10 regulatory axis as a promising biomarker or therapeutic target for clinical diagnosis and treatment of prostate cancer.

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.

The Genetic, Pharmacogenomic, and Immune Landscapes Associated with Protein Expression across Human Cancers

Proteomics is a powerful approach that can rapidly enhance our understanding of cancer development. Detailed characterization of the genetic, pharmacogenomic, and immune landscape in relation to protein expression in patients with cancer could provide new insights into the functional roles of proteins in cancer. By taking advantage of the genotype data from The Cancer Genome Atlas and protein expression data from The Cancer Proteome Atlas, we characterized the effects of genetic variants on protein expression across 31 cancer types and identified approximately 100,000 protein quantitative trait loci (pQTL). Among these, over 8000 pQTLs were associated with patient overall survival. Furthermore, characterization of the impact of protein expression on more than 350 imputed anticancer drug responses in patients revealed nearly 230,000 significant associations. In addition, approximately 21,000 significant associations were identified between protein expression and immune cell abundance. Finally, a user-friendly data portal, GPIP (https://hanlaboratory.com/GPIP), was developed featuring multiple modules that enable researchers to explore, visualize, and browse multidimensional data. This detailed analysis reveals the associations between the proteomic landscape and genetic variation, patient outcome, the immune microenvironment, and drug response across cancer types, providing a resource that may offer valuable clinical insights and encourage further functional investigations of proteins in cancer.

SIGNIFICANCE

Comprehensive characterization of the relationship between protein expression and the genetic, pharmacogenomic, and immune landscape of tumors across cancer types provides a foundation for investigating the role of protein expression in cancer development and treatment.

Leveraging molecular quantitative trait loci to comprehend complex diseases/traits from the omics perspective

Comprehending the molecular basis of quantitative genetic variation is a principal goal for complex diseases or traits. Molecular quantitative trait loci (molQTLs) have made it possible to investigate the effects of genetic variants hiding behind large-scale omics data. A deeper understanding of molQTL is urgently required in light of the multi-dimensionalization of omics data to more fully elucidate the pertinent biological mechanisms. Herein, we reviewed molQTLs with the corresponding resource from the omics perspective and further discussed the integrative strategy of GWAS-molQTL to infer their causal effects. Subsequently, we described the opportunities and challenges encountered by molQTL. The case studies showed that molQTL is essential for complex diseases and traits, whether single- or multi-omics QTLs. Overall, we highlighted the functional significance of genetic variants to employ the discovery of molQTL in complex diseases and traits.

ClipperQTL: ultrafast and powerful eGene identification method

A central task in expression quantitative trait locus (eQTL) analysis is to identify cis-eGenes (henceforth "eGenes"), i.e., genes whose expression levels are regulated by at least one local genetic variant. Among the existing eGene identification methods, FastQTL is considered the gold standard but is computationally expensive as it requires thousands of permutations for each gene. Alternative methods such as eigenMT and TreeQTL have lower power than FastQTL. In this work, we propose ClipperQTL, which reduces the number of permutations needed from thousands to 20 for data sets with large sample sizes (> 450) by using the contrastive strategy developed in Clipper; for data sets with smaller sample sizes, it uses the same permutation-based approach as FastQTL. We show that ClipperQTL performs as well as FastQTL and runs about 500 times faster if the contrastive strategy is used and 50 times faster if the conventional permutation-based approach is used. The R package ClipperQTL is available at https://github.com/heatherjzhou/ClipperQTL.

Genetic analysis of blood molecular phenotypes reveals common properties in the regulatory networks affecting complex traits

We evaluate the shared genetic regulation of mRNA molecules, proteins and metabolites derived from whole blood from 3029 human donors. We find abundant allelic heterogeneity, where multiple variants regulate a particular molecular phenotype, and pleiotropy, where a single variant associates with multiple molecular phenotypes over multiple genomic regions. The highest proportion of share genetic regulation is detected between gene expression and proteins (66.6%), with a further median shared genetic associations across 49 different tissues of 78.3% and 62.4% between plasma proteins and gene expression. We represent the genetic and molecular associations in networks including 2828 known GWAS variants, showing that GWAS variants are more often connected to gene expression in trans than other molecular phenotypes in the network. Our work provides a roadmap to understanding molecular networks and deriving the underlying mechanism of action of GWAS variants using different molecular phenotypes in an accessible tissue.

mGWAS-Explorer 2.0: Causal Analysis and Interpretation of Metabolite-Phenotype Associations

Metabolomics-based genome-wide association studies (mGWAS) are key to understanding the genetic regulations of metabolites in complex phenotypes. We previously developed mGWAS-Explorer 1.0 to link single-nucleotide polymorphisms (SNPs), metabolites, genes and phenotypes for hypothesis generation. It has become clear that identifying potential causal relationships between metabolites and phenotypes, as well as providing deep functional insights, are crucial for further downstream applications. Here, we introduce mGWAS-Explorer 2.0 to support the causal analysis between >4000 metabolites and various phenotypes. The results can be interpreted within the context of semantic triples and molecular quantitative trait loci (QTL) data. The underlying R package is released for reproducible analysis. Using two case studies, we demonstrate that mGWAS-Explorer 2.0 is able to detect potential causal relationships between arachidonic acid and Crohn's disease, as well as between glycine and coronary heart disease.

Human Pangenomics: Promises and Challenges of a Distributed Genomic Reference

A pangenome is a collection of the common and unique genomes that are present in a given species. It combines the genetic information of all the genomes sampled, resulting in a large and diverse range of genetic material. Pangenomic analysis offers several advantages compared to traditional genomic research. For example, a pangenome is not bound by the physical constraints of a single genome, so it can capture more genetic variability. Thanks to the introduction of the concept of pangenome, it is possible to use exceedingly detailed sequence data to study the evolutionary history of two different species, or how populations within a species differ genetically. In the wake of the Human Pangenome Project, this review aims at discussing the advantages of the pangenome around human genetic variation, which are then framed around how pangenomic data can inform population genetics, phylogenetics, and public health policy by providing insights into the genetic basis of diseases or determining personalized treatments, targeting the specific genetic profile of an individual. Moreover, technical limitations, ethical concerns, and legal considerations are discussed.

An autoimmune pleiotropic SNP modulates IRF5 alternative promoter usage through ZBTB3-mediated chromatin looping

Genetic sharing is extensively observed for autoimmune diseases, but the causal variants and their underlying molecular mechanisms remain largely unknown. Through systematic investigation of autoimmune disease pleiotropic loci, we found most of these shared genetic effects are transmitted from regulatory code. We used an evidence-based strategy to functionally prioritize causal pleiotropic variants and identify their target genes. A top-ranked pleiotropic variant, rs4728142, yielded many lines of evidence as being causal. Mechanistically, the rs4728142-containing region interacts with the IRF5 alternative promoter in an allele-specific manner and orchestrates its upstream enhancer to regulate IRF5 alternative promoter usage through chromatin looping. A putative structural regulator, ZBTB3, mediates the allele-specific loop to promote IRF5-short transcript expression at the rs4728142 risk allele, resulting in IRF5 overactivation and M1 macrophage polarization. Together, our findings establish a causal mechanism between the regulatory variant and fine-scale molecular phenotype underlying the dysfunction of pleiotropic genes in human autoimmunity.

Gene-environment interactions and their impact on human health

The molecular processes underlying human health and disease are highly complex. Often, genetic and environmental factors contribute to a given disease or phenotype in a non-additive manner, yielding a gene-environment (G × E) interaction. In this work, we broadly review current knowledge on the impact of gene-environment interactions on human health. We first explain the independent impact of genetic variation and the environment. We next detail well-established G × E interactions that impact human health involving environmental toxicants, pollution, viruses, and sex chromosome composition. We conclude with possibilities and challenges for studying G × E interactions.

Chinese Medicine Phenomics (Chinmedphenomics): Personalized, Precise and Promising

The systematicness of phenomics and Traditional Chinese Medicine (TCM) enable these two disciplines to interlink with each other. This article discussed the similarity in theory and application between TCM and phenomics and illustrates their respective advantages in diagnosis and treatment of diseases, forming a new discipline eventually. Chinese medicine phenomics (Chinmedphenomics) is built on classic TCM, combined with phenomics technology, and the development of which needs the mega cohort with TCM syndrome and the characteristics of precision medicine as well as multi-disciplinary cooperation, which is personalized, precise and promising, providing unique scientific insights into understanding human health.

QTLbase2: an enhanced catalog of human quantitative trait loci on extensive molecular phenotypes

Deciphering the fine-scale molecular mechanisms that shape the genetic effects at disease-associated loci from genome-wide association studies (GWAS) remains challenging. The key avenue is to identify the essential molecular phenotypes that mediate the causal variant and disease under particular biological conditions. Therefore, integrating GWAS signals with context-specific quantitative trait loci (QTLs) (such as different tissue/cell types, disease states, and perturbations) from extensive molecular phenotypes would present important strategies for full understanding of disease genetics. Via persistent curation and systematic data processing of large-scale human molecular trait QTLs (xQTLs), we updated our previous QTLbase database (now QTLbase2, http://mulinlab.org/qtlbase) to comprehensively analyze and visualize context-specific QTLs across 22 molecular phenotypes and over 95 tissue/cell types. Overall, the resource features the following major updates and novel functions: (i) 960 more genome-wide QTL summary statistics from 146 independent studies; (ii) new data for 10 previously uncompiled QTL types; (iii) variant query scope expanded to fit 195 QTL datasets based on whole-genome sequencing; (iv) supports filtering and comparison of QTLs for different biological conditions, such as stimulation types and disease states; (v) a new linkage disequilibrium viewer to facilitate variant prioritization across tissue/cell types and QTL types.

PCA outperforms popular hidden variable inference methods for molecular QTL mapping

BACKGROUND

Estimating and accounting for hidden variables is widely practiced as an important step in molecular quantitative trait locus (molecular QTL, henceforth "QTL") analysis for improving the power of QTL identification. However, few benchmark studies have been performed to evaluate the efficacy of the various methods developed for this purpose.

RESULTS

Here we benchmark popular hidden variable inference methods including surrogate variable analysis (SVA), probabilistic estimation of expression residuals (PEER), and hidden covariates with prior (HCP) against principal component analysis (PCA)-a well-established dimension reduction and factor discovery method-via 362 synthetic and 110 real data sets. We show that PCA not only underlies the statistical methodology behind the popular methods but is also orders of magnitude faster, better-performing, and much easier to interpret and use.

CONCLUSIONS

To help researchers use PCA in their QTL analysis, we provide an R package PCAForQTL along with a detailed guide, both of which are freely available at https://github.com/heatherjzhou/PCAForQTL . We believe that using PCA rather than SVA, PEER, or HCP will substantially improve and simplify hidden variable inference in QTL mapping as well as increase the transparency and reproducibility of QTL research.

SRTdb: an omnibus for human tissue and cancer-specific RNA transcripts

The production of functional mature RNA transcripts from genes undergoes various pre-transcriptional regulation and post-transcriptional modifications. Accumulating studies demonstrated that gene transcription carries out in tissue and cancer type-dependent ways. However, RNA transcript-level specificity analysis in large-scale transcriptomics data across different normal tissue and cancer types is lacking. We applied reference-based de novo transcript assembly and quantification of 27,741 samples across 33 cancer types, 29 tissue types, and 25 cancer cell line types. We totally identified 231,836 specific RNA transcripts (SRTs) across various tissue and cancer types, most of which are found independent of specific genes. Almost half of tumor SRTs are also tissue-specific but in different tissues. Furthermore, we found that 10 ~ 20% of tumor SRTs in most tumor types were testis-specific. The SRT database (SRTdb) was constructed based on these resources. Taking liver cancer as an example, we showed how SRTdb resource is utilized to optimize the identification of RNA transcripts for more precision diagnosis of particular cancers. Our results provide a useful resource for exploring transcript specificity across various cancer and tissue types, and boost the precision medicine for tumor patients.

Interpreting Functional Impact of Genetic Variations by Network QTL for Genotype–Phenotype Association Study

An enormous challenge in the post-genome era is to annotate and resolve the consequences of genetic variation on diverse phenotypes. The genome-wide association study (GWAS) is a well-known method to identify potential genetic loci for complex traits from huge genetic variations, following which it is crucial to identify expression quantitative trait loci (eQTL). However, the conventional eQTL methods usually disregard the systematical role of single-nucleotide polymorphisms (SNPs) or genes, thereby overlooking many network-associated phenotypic determinates. Such a problem motivates us to recognize the network-based quantitative trait loci (QTL), i.e., network QTL (nQTL), which is to detect the cascade association as genotype → network → phenotype rather than conventional genotype → expression → phenotype in eQTL. Specifically, we develop the nQTL framework on the theory and approach of single-sample networks, which can identify not only network traits (e.g., the gene subnetwork associated with genotype) for analyzing complex biological processes but also network signatures (e.g., the interactive gene biomarker candidates screened from network traits) for characterizing targeted phenotype and corresponding subtypes. Our results show that the nQTL framework can efficiently capture associations between SNPs and network traits (i.e., edge traits) in various simulated data scenarios, compared with traditional eQTL methods. Furthermore, we have carried out nQTL analysis on diverse biological and biomedical datasets. Our analysis is effective in detecting network traits for various biological problems and can discover many network signatures for discriminating phenotypes, which can help interpret the influence of nQTL on disease subtyping, disease prognosis, drug response, and pathogen factor association. Particularly, in contrast to the conventional approaches, the nQTL framework could also identify many network traits from human bulk expression data, validated by matched single-cell RNA-seq data in an independent or unsupervised manner. All these results strongly support that nQTL and its detection framework can simultaneously explore the global genotype–network–phenotype associations and the underlying network traits or network signatures with functional impact and importance.

Genetic, Pharmacogenomic and Immune landscapes of enhancer RNAs across human cancers

Enhancer RNAs (eRNA) regulate gene expression and play critical roles in cancer. Using large-scale omics data from The Cancer Genome Atlas (TCGA), we systematically investigated the impact of genetic variants on eRNA expression and identified ~1 million eRNA quantitative trait loci (eRNA-QTL) as cis- and trans-acting. Over 16,000 eRNA-QTLs were associated with patient overall survival. Assessing the impact of eRNAs on >1,000 imputed anti-cancer drug responses across ~10,000 cancer patients revealed > 7 million significant associations. Furthermore, ~240,000 significant associations were identified between eRNA expression and immune cell abundance deconvoluted by TIMER, CIBERSORT, ImmuCellAI, and ImmuneCellGSVA. Finally, a user-friendly data portal was generated: Genetic, Pharmacogenomic, and Immune landscapes of eRNAs (GPIeR, https://hanlab.tamhsc.edu/GPIeR/). GPIeR is a large-scale multi-dimensional data portal that can be used to explore eRNA-associated genetic variants, drug responses, and immune infiltration with the purpose of facilitating functional and clinical investigations of eRNAs in cancer.

Genetic approaches for increasing fitness in endangered species

The global rate of wildlife extinctions is accelerating, and the persistence of many species requires conservation breeding programs. A central paradigm of these programs is to preserve the genetic diversity of the founder populations. However, this may preserve original characteristics that make them vulnerable to extinction. We introduce targeted genetic intervention (TGI) as an alternative approach that promotes traits that enable species to persist in the face of threats by changing the incidence of alleles that impact on fitness. The TGI toolkit includes methods with established efficacy in model organisms and agriculture but are largely untried for conservation, such as synthetic biology and artificial selection. We explore TGI approaches as a species-restoration tool for intractable threats including infectious disease and climate change.

Maintenance of Complex Trait Variation: Classic Theory and Modern Data

Numerous studies have found evidence that GWAS loci experience negative selection, which increases in intensity with the effect size of identified variants. However, there is also accumulating evidence that this selection is not entirely mediated by the focal trait and contains a substantial pleiotropic component. Understanding how selective constraint shapes phenotypic variation requires advancing models capable of balancing these and other components of selection, as well as empirical analyses capable of inferring this balance and how it is generated by the underlying biology. We first review the classic theory connecting phenotypic selection to selection at individual loci as well as approaches and findings from recent analyses of negative selection in GWAS data. We then discuss geometric theories of pleiotropic selection with the potential to guide future modeling efforts. Recent findings revealing the nature of pleiotropic genetic variation provide clues to which genetic relationships are important and should be incorporated into analyses of selection, while findings that effect sizes vary between populations indicate that GWAS measurements could be misleading if effect sizes have also changed throughout human history.

Using Genomic Techniques in Sports and Exercise Science: Current Status and Future Opportunities

ABSTRACT

The past two decades have built on the successes of the Human Genome Project identifying the impact of genetics and genomics on human traits. Given the importance of exercise in the physical and psychological health of individuals across the lifespan, using genomics to understand the impact of genes in the sports medicine field is an emerging field. Given the complexity of the systems involved, high-throughput genomics is required to understand genetic variants, their functions, and ultimately their effect on the body. Consequently, genomic studies have been performed across several domains of sports medicine with varying degrees of success. While the breadth of these is great, they focus largely on the following three areas: 1) performance; 2) injury susceptibility; and 3) sports associated chronic conditions, such as osteoarthritis. Herein, we review literature on genetics and genomics in sports medicine, offer suggestions to bolster existing studies, and suggest ways to ideally impact clinical care.

GPEdit: the genetic and pharmacogenomic landscape of A-to-I RNA editing in cancers

Altered A-to-I RNA editing has been widely observed in many human cancers and some editing sites are associated with drug sensitivity, implicating its therapeutic potential. Increasing evidence has demonstrated that a quantitative trait loci mapping approach is effective to understanding the genetic basis of RNA editing. We systematically performed RNA editing quantitative trait loci (edQTL) analysis in 33 human cancer types for >10 000 cancer samples and identified 320 029 edQTLs. We also identified 1688 ed-QTLs associated with patient overall survival and 4672 ed-QTLs associated with GWAS risk loci. Furthermore, we demonstrated the associations between RNA editing and >1000 anti-cancer drug response with ∼3.5 million significant associations. We developed GPEdit (https://hanlab.uth.edu/GPEdit/) to facilitate a global map of the genetic and pharmacogenomic landscape of RNA editing. GPEdit is a user-friendly and comprehensive database that provides an opportunity for a better understanding of the genetic impact and the effects on drug response of RNA editing in cancers.

Cap analysis of gene expression (CAGE) and noncoding regulatory elements

Cap analysis of gene expression (CAGE) was developed to detect the 5' end of RNA. Trapping of the RNA 5'-cap structure enables the enrichment and selective sequencing of complete transcripts. Upscaled high-throughput versions of CAGE have enabled the genome-wide identification of transcription start sites, including transcriptionally active promoters and enhancers. CAGE sequencing can be exploited to draw comprehensive maps of active genomic regulatory elements in a cell type- and activation-specific manner. The cells of the immune system are among the best candidates to be analyzed in humans, since they are easily accessible. In this review, we discuss how CAGE data are instrumental for integrative analyses with quantitative trait loci and omics data, and their usefulness in the mechanistic interpretation of the effects of genetic variations over the entire human genome. Integrating CAGE data with the currently available omics information will contribute to better understanding of the genome-wide association study variants that lie outside of annotated genes, deepening our knowledge on human diseases, and enabling the targeted design of more specific therapeutic interventions.

A panoramic view of proteomics and multiomics in precision health

Health is often qualitatively defined as a status free from disease and its quantitative definition requires finding the boundary separating health from pathological conditions. Since many complex diseases have a strong genetic component, substantial efforts have been made to sequence large-scale personal genomes; however, we are not yet able to effectively quantify health status from personal genomes. Since mutational impacts are ultimately manifested at the protein level, we envision that introducing a panoramic proteomic view of complex diseases will allow us to mechanistically understand the molecular etiologies of human diseases. In this perspective article, we will highlight key proteomic approaches to identify pathogenic mutations and map their convergent pathways underlying disease pathogenesis and the integration of omics data at multiple levels to define the borderline between health and disease.

CoffeeProt: an online tool for correlation and functional enrichment of systems genetics data

The integration of genomics, transcriptomics, proteomics and phenotypic traits across genetically diverse populations is a powerful approach to discover novel biological regulators. The increasing volume of complex data require new and easy-to-use tools accessible to a variety of scientists for the discovery and visualization of functionally relevant associations. To meet this requirement, we developed CoffeeProt, an open-source tool that analyses genetic variants associated to protein networks, other omics datatypes and phenotypic traits. CoffeeProt uses transcriptomics or proteomics data to perform correlation network analyses and annotates results with protein-protein interactions, subcellular localisations and drug associations. It then integrates genetic variants associated with gene expression (eQTLs) or protein abundance (pQTLs) and includes predictions of the potential consequences of variants on gene function. Finally, genetic variants are co-mapped to molecular or phenotypic traits either provided by the user or retrieved directly from publicly available GWAS results. We demonstrate its utility with the analysis of mouse and human population data enabling the rapid identification of genetic variants associated with druggable proteins and clinical traits. We expect that CoffeeProt will serve the systems genetics and basic science research communities, leading to the discovery of novel biologically relevant associations. CoffeeProt is available at www.coffeeprot.com.

Unveiling the Pathogenesis of Psychiatric Disorders Using Network Models

Psychiatric disorders are complex brain disorders with a high degree of genetic heterogeneity, affecting millions of people worldwide. Despite advances in psychiatric genetics, the underlying pathogenic mechanisms of psychiatric disorders are still largely elusive, which impedes the development of novel rational therapies. There has been accumulating evidence suggesting that the genetics of complex disorders can be viewed through an omnigenic lens, which involves contextualizing genes in highly interconnected networks. Thus, applying network-based multi-omics integration methods could cast new light on the pathophysiology of psychiatric disorders. In this review, we first provide an overview of the recent advances in psychiatric genetics and highlight gaps in translating molecular associations into mechanistic insights. We then present an overview of network methodologies and review previous applications of network methods in the study of psychiatric disorders. Lastly, we describe the potential of such methodologies within a multi-tissue, multi-omics approach, and summarize the future directions in adopting diverse network approaches.

Alternative polyadenylation-associated loci interpret human traits and diseases

Alternative polyadenylation (APA)-associated genetic variants have been proposed to impact diverse human phenotypes and disorders. In a recent study, Li et al. established a landscape of 3'-untranslated region (UTR) APA quantitative trait loci (3'aQTLs) across multiple human tissues, revealing substantial 3'aQTLs that contribute to complex human traits and diseases.

Beyond association: successes and challenges in linking non-coding genetic variation to functional consequences that modulate Alzheimer's disease risk

Alzheimer's disease (AD) is the most common type of dementia, affecting millions of people worldwide; however, no disease-modifying treatments are currently available. Genome-wide association studies (GWASs) have identified more than 40 loci associated with AD risk. However, most of the disease-associated variants reside in non-coding regions of the genome, making it difficult to elucidate how they affect disease susceptibility. Nonetheless, identification of the regulatory elements, genes, pathways and cell type/tissue(s) impacted by these variants to modulate AD risk is critical to our understanding of disease pathogenesis and ability to develop effective therapeutics. In this review, we provide an overview of the methods and approaches used in the field to identify the functional effects of AD risk variants in the causal path to disease risk modification as well as describe the most recent findings. We first discuss efforts in cell type/tissue prioritization followed by recent progress in candidate causal variant and gene nomination. We discuss statistical methods for fine-mapping as well as approaches that integrate multiple levels of evidence, such as epigenomic and transcriptomic data, to identify causal variants and risk mechanisms of AD-associated loci. Additionally, we discuss experimental approaches and data resources that will be needed to validate and further elucidate the effects of these variants and genes on biological pathways, cellular phenotypes and disease risk. Finally, we discuss future steps that need to be taken to ensure that AD GWAS functional mapping efforts lead to novel findings and bring us closer to finding effective treatments for this devastating disease.

Mouse Genetic Reference Populations: Cellular Platforms for Integrative Systems Genetics

Interrogation of disease-relevant cellular and molecular traits exhibited by genetically diverse cell populations enables in vitro systems genetics approaches for uncovering the basic properties of cellular function and identity. Primary cells, stem cells, and organoids derived from genetically diverse mouse strains, such as Collaborative Cross and Diversity Outbred populations, offer the opportunity for parallel in vitro/in vivo screening. These panels provide genetic resolution for variant discovery and functional characterization, as well as disease modeling and in vivo validation capabilities. Here we review mouse cellular systems genetics approaches for characterizing the influence of genetic variation on signaling networks and phenotypic diversity, and we discuss approaches for data integration and cross-species validation.

Deep Learning in Head and Neck Tumor Multiomics Diagnosis and Analysis: Review of the Literature

Head and neck tumors are the sixth most common neoplasms. Multiomics integrates multiple dimensions of clinical, pathologic, radiological, and biological data and has the potential for tumor diagnosis and analysis. Deep learning (DL), a type of artificial intelligence (AI), is applied in medical image analysis. Among the DL techniques, the convolution neural network (CNN) is used for image segmentation, detection, and classification and in computer-aided diagnosis. Here, we reviewed multiomics image analysis of head and neck tumors using CNN and other DL neural networks. We also evaluated its application in early tumor detection, classification, prognosis/metastasis prediction, and the signing out of the reports. Finally, we highlighted the challenges and potential of these techniques.

Multi-parametric characterization of drug effects on cells

We present here a novel multi-parametric approach for the characterization of multiple cellular features, using images acquired by high-throughput and high-definition light microscopy. We specifically used this approach for deep and unbiased analysis of the effects of a drug library on five cultured cell lines. The presented method enables the acquisition and analysis of millions of images, of treated and control cells, followed by an automated identification of drugs inducing strong responses, evaluating the median effect concentrations and those cellular properties that are most highly affected by the drug. The tools described here provide standardized quantification of multiple attributes for systems level dissection of complex functions in normal and diseased cells, using multiple perturbations. Such analysis of cells, derived from pathological samples, may help in the diagnosis and follow-up of treatment in patients.

Where Are the Disease-Associated eQTLs?

Most disease-associated variants, although located in putatively regulatory regions, do not have detectable effects on gene expression. One explanation could be that we have not examined gene expression in the cell types or conditions that are most relevant for disease. Even large-scale efforts to study gene expression across tissues are limited to human samples obtained opportunistically or postmortem, mostly from adults. In this review we evaluate recent findings and suggest an alternative strategy, drawing on the dynamic and highly context-specific nature of gene regulation. We discuss new technologies that can extend the standard regulatory mapping framework to more diverse, disease-relevant cell types and states.

aScan: A Novel Method for the Study of Allele Specific Expression in Single Individuals

In diploid organisms, two copies of each allele are normally inherited from parents. Paternal and maternal alleles can be regulated and expressed unequally, which is referred to as allele-specific expression (ASE). In this work, we present aScan, a novel method for the identification of ASE from the analysis of matched individual genomic and RNA sequencing data. By performing extensive analyses of both real and simulated data, we demonstrate that aScan can correctly identify ASE with high accuracy and sensitivity in different experimental settings. Additionally, by applying our method to a small cohort of individuals that are not included in publicly available databases of human genetic variation, we outline the value of possible applications of ASE analysis in single individuals for deriving a more accurate annotation of "private" low-frequency genetic variants associated with regulatory effects on transcription. All in all, we believe that aScan will represent a beneficial addition to the set of bioinformatics tools for the analysis of ASE. Finally, while our method was initially conceived for the analysis of RNA-seq data, it can in principle be applied to any quantitative NGS assay for which matched genotypic and expression data are available. AVAILABILITY: aScan is currently available in the form of an open source standalone software package at: https://github.com/Federico77z/aScan/. aScan version 1.0.3, available at https://github.com/Federico77z/aScan/releases/tag/1.0.3, has been used for all the analyses included in this manuscript. A Docker image of the tool has also been made available at https://github.com/pmandreoli/aScanDocker.

Functional genomics of autoimmune diseases

For more than a decade, genome-wide association studies have been applied to autoimmune diseases and have expanded our understanding on the pathogeneses. Genetic risk factors associated with diseases and traits are essentially causative. However, elucidation of the biological mechanism of disease from genetic factors is challenging. In fact, it is difficult to identify the causal variant among multiple variants located on the same haplotype or linkage disequilibrium block and thus the responsible biological genes remain elusive. Recently, multiple studies have revealed that the majority of risk variants locate in the non-coding region of the genome and they are the most likely to regulate gene expression such as quantitative trait loci. Enhancer, promoter and long non-coding RNA appear to be the main target mechanisms of the risk variants. In this review, we discuss functional genetics to challenge these puzzles.

Unique roles of rare variants in the genetics of complex diseases in humans

Genome-wide association studies have identified >10,000 genetic variants associated with various phenotypes and diseases. Although the majority are common variants, rare variants with >0.1% of minor allele frequency have been investigated by imputation and using disease-specific custom SNP arrays. Rare variants sequencing analysis mainly revealed have played unique roles in the genetics of complex diseases in humans due to their distinctive features, in contrast to common variants. Unique roles are hypothesis-free evidence for gene causality, a precise target of functional analysis for understanding disease mechanisms, a new favorable target for drug development, and a genetic marker with high disease risk for personalized medicine. As whole-genome sequencing continues to identify more rare variants, the roles associated with rare variants will also increase. However, a better estimation of the functional impact of rare variants across whole genome is needed to enhance their contribution to improvements in human health.

The genetic and pharmacogenomic landscape of snoRNAs in human cancer

Emerging evidence has revealed significant roles for small nucleolar RNAs (snoRNAs) in tumorigenesis. However, the genetic and pharmacogenomic landscape of snoRNAs has not been characterized. Using the genotype and snoRNA expression data from The Cancer Genome Atlas, we characterized the effects of genetic variants on snoRNAs across 29 cancer types and further linked related alleles with patient survival as well as genome-wide association study risk loci. Furthermore, we characterized the impact of snoRNA expression on drug response in patients to facilitate the clinical utility of snoRNAs in cancer. We also developed a user-friendly data resource, GPSno (http://hanlab.uth.edu/GPSno), with multiple modules for researchers to visualize, browse, and download multi-dimensional data. Our study provides a comprehensive genetic and pharmacogenomic landscape of snoRNAs, which will shed light on future clinical considerations for the development of snoRNA-based targeted therapies.