Projecting genetic associations through gene expression patterns highlights disease etiology and drug mechanisms.

Genes act in concert with each other in specific contexts to perform their functions. Determining how these genes influence complex traits requires a mechanistic understanding of expression regulation across different conditions. It has been shown that this insight is critical for developing new therapies. Transcriptome-wide association studies have helped uncover the role of individual genes in disease-relevant mechanisms. However, modern models of the architecture of complex traits predict that gene-gene interactions play a crucial role in disease origin and progression. Here we introduce PhenoPLIER, a computational approach that maps gene-trait associations and pharmacological perturbation data into a common latent representation for a joint analysis. This representation is based on modules of genes with similar expression patterns across the same conditions. We observe that diseases are significantly associated with gene modules expressed in relevant cell types, and our approach is accurate in predicting known drug-disease pairs and inferring mechanisms of action. Furthermore, using a CRISPR screen to analyze lipid regulation, we find that functionally important players lack associations but are prioritized in trait-associated modules by PhenoPLIER. By incorporating groups of co-expressed genes, PhenoPLIER can contextualize genetic associations and reveal potential targets missed by single-gene strategies.

Frequencies of pharmacogenomic alleles across biogeographic groups in a large-scale biobank.

Pharmacogenomics (PGx) is an integral part of precision medicine and contributes to the maximization of drug efficacy and reduction of adverse drug event risk. Accurate information on PGx allele frequencies improves the implementation of PGx. Nonetheless, curating such information from published allele data is time and resource intensive. The limited number of allelic variants in most studies leads to an underestimation of certain alleles. We applied the Pharmacogenomics Clinical Annotation Tool (PharmCAT) on an integrated 200K UK Biobank genetic dataset (N = 200,044). Based on PharmCAT results, we estimated PGx frequencies (alleles, diplotypes, phenotypes, and activity scores) for 17 pharmacogenes in five biogeographic groups: European, Central/South Asian, East Asian, Afro-Caribbean, and Sub-Saharan African. PGx frequencies were distinct for each biogeographic group. Even biogeographic groups with similar proportions of phenotypes were driven by different sets of dominant PGx alleles. PharmCAT also identified "no-function" alleles that were rare or seldom tested in certain groups by previous studies, e.g., SLCO1B1∗31 in the Afro-Caribbean (3.0%) and Sub-Saharan African (3.9%) groups. Estimated PGx frequencies are disseminated via the PharmGKB (The Pharmacogenomics Knowledgebase: www.pharmgkb.org). We demonstrate that genetic biobanks such as the UK Biobank are a robust resource for estimating PGx frequencies. Improving our understanding of PGx allele and phenotype frequencies provides guidance for future PGx studies and clinical genetic test panel design, and better serves individuals from wider biogeographic backgrounds.

Gene burden analysis identifies genes associated with increased risk and severity of adult-onset hearing loss in a diverse hospital-based cohort.

Loss or absence of hearing is common at both extremes of human lifespan, in the forms of congenital deafness and age-related hearing loss. While these are often studied separately, there is increasing evidence that their genetic basis is at least partially overlapping. In particular, both common and rare variants in genes associated with monogenic forms of hearing loss also contribute to the more polygenic basis of age-related hearing loss. Here, we directly test this model in the Penn Medicine BioBank-a healthcare system cohort of around 40,000 individuals with linked genetic and electronic health record data. We show that increased burden of predicted deleterious variants in Mendelian hearing loss genes is associated with increased risk and severity of adult-onset hearing loss. As a specific example, we identify one gene-TCOF1, responsible for a syndromic form of congenital hearing loss-in which deleterious variants are also associated with adult-onset hearing loss. We also identify four additional novel candidate genes (COL5A1, HMMR, RAPGEF3, and NNT) in which rare variant burden may be associated with hearing loss. Our results confirm that rare variants in Mendelian hearing loss genes contribute to polygenic risk of hearing loss, and emphasize the utility of healthcare system cohorts to study common complex traits and diseases.

How to Run the Pharmacogenomics Clinical Annotation Tool (PharmCAT).

Pharmacogenomics (PGx) investigates the genetic influence on drug response and is an integral part of precision medicine. While PGx testing is becoming more common in clinical practice and may be reimbursed by Medicare/Medicaid and commercial insurance, interpreting PGx testing results for clinical decision support is still a challenge. The Pharmacogenomics Clinical Annotation Tool (PharmCAT) has been designed to tackle the need for transparent, automatic interpretations of patient genetic data. PharmCAT incorporates a patient's genotypes, annotates PGx information (allele, genotype, and phenotype), and generates a report with PGx guideline recommendations from the Clinical Pharmacogenetics Implementation Consortium (CPIC) and/or the Dutch Pharmacogenetics Working Group (DPWG). PharmCAT has introduced new features in the last 2 years, including a variant call format (VCF) Preprocessor, the inclusion of DPWG guidelines, and functionalities for PGx research. For example, researchers can use the VCF Preprocessor to prepare biobank-scale data for PharmCAT. In addition, PharmCAT enables the assessment of novel partial and combination alleles that are composed of known PGx variants and can call CYP2D6 genotypes based on single and deletions in the input VCF file. This tutorial provides materials and detailed step-by-step instructions for how to use PharmCAT in a versatile way that can be tailored to users' individual needs.

Evaluating the frequency and the impact of pharmacogenetic alleles in an ancestrally diverse Biobank population.

BACKGROUND: Pharmacogenomics (PGx) aims to utilize a patient's genetic data to enable safer and more effective prescribing of medications. The Clinical Pharmacogenetics Implementation Consortium (CPIC) provides guidelines with strong evidence for 24 genes that affect 72 medications. Despite strong evidence linking PGx alleles to drug response, there is a large gap in the implementation and return of actionable pharmacogenetic findings to patients in standard clinical practice. In this study, we evaluated opportunities for genetically guided medication prescribing in a diverse health system and determined the frequencies of actionable PGx alleles in an ancestrally diverse biobank population. METHODS: A retrospective analysis of the Penn Medicine electronic health records (EHRs), which includes ~ 3.3 million patients between 2012 and 2020, provides a snapshot of the trends in prescriptions for drugs with genotype-based prescribing guidelines ('CPIC level A or B') in the Penn Medicine health system. The Penn Medicine BioBank (PMBB) consists of a diverse group of 43,359 participants whose EHRs are linked to genome-wide SNP array and whole exome sequencing (WES) data. We used the Pharmacogenomics Clinical Annotation Tool (PharmCAT), to annotate PGx alleles from PMBB variant call format (VCF) files and identify samples with actionable PGx alleles. RESULTS: We identified ~ 316.000 unique patients that were prescribed at least 2 drugs with CPIC Level A or B guidelines. Genetic analysis in PMBB identified that 98.9% of participants carry one or more PGx actionable alleles where treatment modification would be recommended. After linking the genetic data with prescription data from the EHR, 14.2% of participants (n = 6157) were prescribed medications that could be impacted by their genotype (as indicated by their PharmCAT report). For example, 856 participants received clopidogrel who carried CYP2C19 reduced function alleles, placing them at increased risk for major adverse cardiovascular events. When we stratified by genetic ancestry, we found disparities in PGx allele frequencies and clinical burden. Clopidogrel users of Asian ancestry in PMBB had significantly higher rates of CYP2C19 actionable alleles than European ancestry users of clopidrogrel (p < 0.0001, OR = 3.68). CONCLUSIONS: Clinically actionable PGx alleles are highly prevalent in our health system and many patients were prescribed medications that could be affected by PGx alleles. These results illustrate the potential utility of preemptive genotyping for tailoring of medications and implementation of PGx into routine clinical care.

A Genome-First Approach to Rare Variants in Dominant Postlingual Hearing Loss Genes in a Large Adult Population.

OBJECTIVE: To investigate the importance of rare variants in adult-onset hearing loss. STUDY DESIGN: Genomic association study. SETTING: Large biobank from tertiary care center. METHODS: We investigated rare variants (minor allele frequency <5%) in 42 autosomal dominant (DFNA) postlingual hearing loss (HL) genes in 16,657 unselected individuals in the Penn Medicine Biobank. We determined the prevalence of known pathogenic and predicted deleterious variants in subjects with audiometric-proven sensorineural hearing loss. We scanned across known postlingual DFNA HL genes to determine those most significantly contributing to the phenotype. We replicated findings in an independent cohort (UK Biobank). RESULTS: While rare individually, when considering the accumulation of variants in all postlingual DFNA genes, more than 90% of participants carried at least 1 rare variant. Rare variants predicted to be deleterious were enriched in adults with audiometric-proven hearing loss (pure-tone average >25 dB; P = .015). Patients with a rare predicted deleterious variant had an odds ratio of 1.27 for HL compared with genotypic controls (P = .029). Gene burden in DIABLO, PTPRQ, TJP2, and POU4F3 were independently associated with sensorineural hearing loss. CONCLUSION: Although prior reports have focused on common variants, we find that rare predicted deleterious variants in DFNA postlingual HL genes are enriched in patients with adult-onset HL in a large health care system population. We show the value of investigating rare variants to uncover hearing loss phenotypes related to implicated genes.

How Does the "Cookie-Bite" Audiogram Shape Perform in Discriminating Genetic Hearing Loss in Adults?

"Cookie-bite" or U-shaped audiograms-specifically, those showing midfrequency sensorineural hearing loss (HL)-are traditionally taught to be associated with genetic HL; however, their utility as a screening tool has not been reported. We aim to determine the performance of a cookie-bite audiogram shape in stratifying patients carrying putative loss-of-function variants in known HL genes from wild-type controls. We merged audiometric and exome sequencing data from adults enrolled in a large biobank at a tertiary care center. Of 321 patients, 50 carried a putative loss-of-function variant in an HL gene. The cookie-bite shape was present in 9 of those patients, resulting in low sensitivity (18%) and positive predictive value (15%) in stratifying genetic carrier status; 84% of patients with a cookie-bite audiogram did not carry a genetic variant. A cookie-bite audiogram should not be used to screen adults for possible genetic testing.

An Investigation of the Knowledge Overlap between Pharmacogenomics and Disease Genetics.

Precision medicine faces many challenges, including the gap of knowledge between disease genetics and pharmacogenomics (PGx). Disease genetics interprets the pathogenicity of genetic variants for diagnostic purposes, while PGx investigates the genetic influences on drug responses. Ideally, the quality of health care would be improved from the point of disease diagnosis to drug prescribing if PGx is integrated with disease genetics in clinical care. However, PGx genes or variants are usually not reported as a secondary finding even if they are included in a clinical genetic test for diagnostic purposes. This happens even though the detection of PGx variants can provide valuable drug prescribing recommendations. One underlying reason is the lack of systematic classification of the knowledge overlap between PGx and disease genetics. Here, we address this issue by analyzing gene and genetic variant annotations from multiple expert-curated knowledge databases, including PharmGKB, CPIC, ClinGen and ClinVar. We further classified genes based on the strength of evidence supporting a gene's pathogenic role or PGx effect as well as the level of clinical actionability of a gene. Twenty-six genes were found to have pathogenic variation associated with germline diseases as well as strong evidence for a PGx association. These genes were classified into four sub-categories based on the distinct connection between the gene's pathogenic role and PGx effect. Moreover, we have also found thirteen RYR1 genetic variants that were annotated as pathogenic and at the same time whose PGx effect was supported by a preponderance of evidence and given drug prescribing recommendations. Overall, we identified a nontrivial number of gene and genetic variant overlaps between disease genetics and PGx, which laid out a foundation for combining PGx and disease genetics to improve clinical care from disease diagnoses to drug prescribing and adherence.

From GWAS to Gene: Transcriptome-Wide Association Studies and Other Methods to Functionally Understand GWAS Discoveries.

Since their inception, genome-wide association studies (GWAS) have identified more than a hundred thousand single nucleotide polymorphism (SNP) loci that are associated with various complex human diseases or traits. The majority of GWAS discoveries are located in non-coding regions of the human genome and have unknown functions. The valley between non-coding GWAS discoveries and downstream affected genes hinders the investigation of complex disease mechanism and the utilization of human genetics for the improvement of clinical care. Meanwhile, advances in high-throughput sequencing technologies reveal important genomic regulatory roles that non-coding regions play in the transcriptional activities of genes. In this review, we focus on data integrative bioinformatics methods that combine GWAS with functional genomics knowledge to identify genetically regulated genes. We categorize and describe two types of data integrative methods. First, we describe fine-mapping methods. Fine-mapping is an exploratory approach that calibrates likely causal variants underneath GWAS signals. Fine-mapping methods connect GWAS signals to potentially causal genes through statistical methods and/or functional annotations. Second, we discuss gene-prioritization methods. These are hypothesis generating approaches that evaluate whether genetic variants regulate genes via certain genetic regulatory mechanisms to influence complex traits, including colocalization, mendelian randomization, and the transcriptome-wide association study (TWAS). TWAS is a gene-based association approach that investigates associations between genetically regulated gene expression and complex diseases or traits. TWAS has gained popularity over the years due to its ability to reduce multiple testing burden in comparison to other variant-based analytic approaches. Multiple types of TWAS methods have been developed with varied methodological designs and biological hypotheses over the past 5 years. We dive into discussions of how TWAS methods differ in many aspects and the challenges that different TWAS methods face. Overall, TWAS is a powerful tool for identifying complex trait-associated genes. With the advent of single-cell sequencing, chromosome conformation capture, gene editing technologies, and multiplexing reporter assays, we are expecting a more comprehensive understanding of genomic regulation and genetically regulated genes underlying complex human diseases and traits in the future.

Tissue specificity-aware TWAS (TSA-TWAS) framework identifies novel associations with metabolic, immunologic, and virologic traits in HIV-positive adults.

As a type of relatively new methodology, the transcriptome-wide association study (TWAS) has gained interest due to capacity for gene-level association testing. However, the development of TWAS has outpaced statistical evaluation of TWAS gene prioritization performance. Current TWAS methods vary in underlying biological assumptions about tissue specificity of transcriptional regulatory mechanisms. In a previous study from our group, this may have affected whether TWAS methods better identified associations in single tissues versus multiple tissues. We therefore designed simulation analyses to examine how the interplay between particular TWAS methods and tissue specificity of gene expression affects power and type I error rates for gene prioritization. We found that cross-tissue identification of expression quantitative trait loci (eQTLs) improved TWAS power. Single-tissue TWAS (i.e., PrediXcan) had robust power to identify genes expressed in single tissues, but, often found significant associations in the wrong tissues as well (therefore had high false positive rates). Cross-tissue TWAS (i.e., UTMOST) had overall equal or greater power and controlled type I error rates for genes expressed in multiple tissues. Based on these simulation results, we applied a tissue specificity-aware TWAS (TSA-TWAS) analytic framework to look for gene-based associations with pre-treatment laboratory values from AIDS Clinical Trial Group (ACTG) studies. We replicated several proof-of-concept transcriptionally regulated gene-trait associations, including UGT1A1 (encoding bilirubin uridine diphosphate glucuronosyltransferase enzyme) and total bilirubin levels (p = 3.59×10-12), and CETP (cholesteryl ester transfer protein) with high-density lipoprotein cholesterol (p = 4.49×10-12). We also identified several novel genes associated with metabolic and virologic traits, as well as pleiotropic genes that linked plasma viral load, absolute basophil count, and/or triglyceride levels. By highlighting the advantages of different TWAS methods, our simulation study promotes a tissue specificity-aware TWAS analytic framework that revealed novel aspects of HIV-related traits.

Genome-first approach to rare EYA4 variants and cardio-auditory phenotypes in adults.

While newborns and children with hearing loss are routinely offered genetic testing, adults are rarely clinically tested for a genetic etiology. One clinically actionable result from genetic testing in children is the discovery of variants in syndromic hearing loss genes. EYA4 is a known hearing loss gene which is also involved in important pathways in cardiac tissue. The pleiotropic effects of rare EYA4 variants are poorly understood and their prevalence in a large cohort has not been previously reported. We investigated cardio-auditory phenotypes in 11,451 individuals in a large biobank using a rare variant, genome-first approach to EYA4. We filtered 256 EYA4 variants carried by 6737 participants to 26 rare and predicted deleterious variants carried by 42 heterozygotes. We aggregated predicted deleterious EYA4 gene variants into a combined variable (i.e. "gene burden") and performed association studies across phenotypes compared to wildtype controls. We validated findings with replication in three independent cohorts and human tissue expression data. EYA4 gene burden was significantly associated with audiometric-proven HL (p = [Formula: see text], Mobitz Type II AV block (p = [Formula: see text]) and the syndromic presentation of HL and primary cardiomyopathy (p = 0.0194). Analyses on audiogram, echocardiogram, and electrocardiogram data validated these associations. Prior reports have focused on identifying variants in families with severe or syndromic phenotypes. In contrast, we found, using a genotype-first approach, that gene burden in EYA4 is associated with more subtle cardio-auditory phenotypes in an adult medical biobank population, including cardiac conduction disorders which have not been previously reported. We show the value of using a focused approach to uncover human disease related to pleiotropic gene variants and suggest a role for genetic testing in adults presenting with hearing loss.

Influence of tissue context on gene prioritization for predicted transcriptome-wide association studies.

Transcriptome-wide association studies (TWAS) have recently gained great attention due to their ability to prioritize complex trait-associated genes and promote potential therapeutics development for complex human diseases. TWAS integrates genotypic data with expression quantitative trait loci (eQTLs) to predict genetically regulated gene expression components and associates predictions with a trait of interest. As such, TWAS can prioritize genes whose differential expressions contribute to the trait of interest and provide mechanistic explanation of complex trait(s). Tissue-specific eQTL information grants TWAS the ability to perform association analysis on tissues whose gene expression profiles are otherwise hard to obtain, such as liver and heart. However, as eQTLs are tissue context-dependent, whether and how the tissue-specificity of eQTLs influences TWAS gene prioritization has not been fully investigated. In this study, we addressed this question by adopting two distinct TWAS methods, PrediXcan and UTMOST, which assume single tissue and integrative tissue effects of eQTLs, respectively. Thirty-eight baseline laboratory traits in 4,360 antiretroviral treatment-naïve individuals from the AIDS Clinical Trials Group (ACTG) studies comprised the input dataset for TWAS. We performed TWAS in a tissue-specific manner and obtained a total of 430 significant gene-trait associations (q-value < 0.05) across multiple tissues. Single tissue-based analysis by PrediXcan contributed 116 of the 430 associations including 64 unique gene-trait pairs in 28 tissues. Integrative tissue-based analysis by UTMOST found the other 314 significant associations that include 50 unique gene-trait pairs across all 44 tissues. Both analyses were able to replicate some associations identified in past variant-based genome-wide association studies (GWAS), such as high-density lipoprotein (HDL) and CETP (PrediXcan, q-value = 3.2e-16). Both analyses also identified novel associations. Moreover, single tissue-based and integrative tissuebased analysis shared 11 of 103 unique gene-trait pairs, for example, PSRC1-low-density lipoprotein (PrediXcan's lowest q-value = 8.5e-06; UTMOST's lowest q-value = 1.8e-05). This study suggests that single tissue-based analysis may have performed better at discovering gene-trait associations when combining results from all tissues. Integrative tissue-based analysis was better at prioritizing genes in multiple tissues and in trait-related tissue. Additional exploration is needed to confirm this conclusion. Finally, although single tissue-based and integrative tissue-based analysis shared significant novel discoveries, tissue context-dependency of eQTLs impacted TWAS gene prioritization. This study provides preliminary data to support continued work on tissue contextdependency of eQTL studies and TWAS.

Collective feature selection to identify crucial epistatic variants.

BACKGROUND: Machine learning methods have gained popularity and practicality in identifying linear and non-linear effects of variants associated with complex disease/traits. Detection of epistatic interactions still remains a challenge due to the large number of features and relatively small sample size as input, thus leading to the so-called "short fat data" problem. The efficiency of machine learning methods can be increased by limiting the number of input features. Thus, it is very important to perform variable selection before searching for epistasis. Many methods have been evaluated and proposed to perform feature selection, but no single method works best in all scenarios. We demonstrate this by conducting two separate simulation analyses to evaluate the proposed collective feature selection approach. RESULTS: Through our simulation study we propose a collective feature selection approach to select features that are in the "union" of the best performing methods. We explored various parametric, non-parametric, and data mining approaches to perform feature selection. We choose our top performing methods to select the union of the resulting variables based on a user-defined percentage of variants selected from each method to take to downstream analysis. Our simulation analysis shows that non-parametric data mining approaches, such as MDR, may work best under one simulation criteria for the high effect size (penetrance) datasets, while non-parametric methods designed for feature selection, such as Ranger and Gradient boosting, work best under other simulation criteria. Thus, using a collective approach proves to be more beneficial for selecting variables with epistatic effects also in low effect size datasets and different genetic architectures. Following this, we applied our proposed collective feature selection approach to select the top 1% of variables to identify potential interacting variables associated with Body Mass Index (BMI) in ~ 44,000 samples obtained from Geisinger's MyCode Community Health Initiative (on behalf of DiscovEHR collaboration). CONCLUSIONS: In this study, we were able to show that selecting variables using a collective feature selection approach could help in selecting true positive epistatic variables more frequently than applying any single method for feature selection via simulation studies. We were able to demonstrate the effectiveness of collective feature selection along with a comparison of many methods in our simulation analysis. We also applied our method to identify non-linear networks associated with obesity.

Evaluation of PrediXcan for prioritizing GWAS associations and predicting gene expression.

Genome-wide association studies (GWAS) have been successful in facilitating the understanding of genetic architecture behind human diseases, but this approach faces many challenges. To identify disease-related loci with modest to weak effect size, GWAS requires very large sample sizes, which can be computational burdensome. In addition, the interpretation of discovered associations remains difficult. PrediXcan was developed to help address these issues. With built in SNP-expression models, PrediXcan is able to predict the expression of genes that are regulated by putative expression quantitative trait loci (eQTLs), and these predicted expression levels can then be used to perform gene-based association studies. This approach reduces the multiple testing burden from millions of variants down to several thousand genes. But most importantly, the identified associations can reveal the genes that are under regulation of eQTLs and consequently involved in disease pathogenesis. In this study, two of the most practical functions of PrediXcan were tested: 1) predicting gene expression, and 2) prioritizing GWAS results. We tested the prediction accuracy of PrediXcan by comparing the predicted and observed gene expression levels, and also looked into some potential influential factors and a filter criterion with the aim of improving PrediXcan performance. As for GWAS prioritization, predicted gene expression levels were used to obtain gene-trait associations, and background regions of significant associations were examined to decrease the likelihood of false positives. Our results showed that 1) PrediXcan predicted gene expression levels accurately for some but not all genes; 2) including more putative eQTLs into prediction did not improve the prediction accuracy; and 3) integrating predicted gene expression levels from the two PrediXcan whole blood models did not eliminate false positives. Still, PrediXcan was able to prioritize GWAS associations that were below the genome-wide significance threshold in GWAS, while retaining GWAS significant results. This study suggests several ways to consider PrediXcan's performance that will be of value to eQTL and complex human disease research.

[Studied on the identification of Fructus amomi from its confusable varieties by Fourier transform infrared spectroscopy].

This paper uses Fourier transform infrared spectrometer with OMNI sampler to distinguish Fructus amomi from their confusable varieties, i. e. Amomum aurantiacum H. T. Tsai et S. W. Zhao, Amomum chinense Chun ex T. L. Wu, Alpinia chinensis (Ketz.) Rosc and Alpinia japonica (Thunb.) Miq. IRs of Amomum villosum Lour., Amomum longiligulare T. L. Wu and Amomum villosum Lour. Var xanthioides T. L. Wu et Senjen are resemble, and they are markedly different from the FTIR of the confusable varieties. Repeat experiments were processed with different samples of the same set, and the probability is 1.000. The result shows that FTIR can be directly used to distinguish Fructus amomi from their confusable varieties.

[Direct identification of Myristica fragrans and Myristica sp. by FTIR].

OBJECTIVE: To directly and accurately identify Myristica fragrans and Myristica sp. METHODS: Fourier transform infrare(FTIR) spectrum method was used. RESULT: There were obvious differences between the FTIR spectrums of above-mentioned plants. CONCLUSION: Myristica fragrans and Myristica sp. were identificated by FTIR directly, fast and accurately.

[Direct identification of Celosia argentea from its confused species by FTIR].

OBJECTIVE: To directly and accurately identify Celosia argentea L. from its confused species. METHODS: Fourier Transform Infrared Spectrum. RESULT: Obvious characteristics for the identification in FTIR were revealed, which can be used to identify Celosia argentea L. and its confused species, such as Celosia cristata L., Amoranthus retroflexus. A. tricolor L. and A. patulus Bertol. CONCLUSION: Celosia argentea L. and its confused species were identified by FTIR directly, fastly and accurately.