General Approaches for Combining Multiple Rare Variant Associate Tests Provide Improved Power Across a Wider Range of Genetic Architecture

In the wake of the widespread availability of genome sequencing data made possible by way of nextgeneration technologies, a flood of gene‐based rare variant tests have been proposed. Most methods claim superior power against particular genetic architectures. However, an important practical issue remains for the applied researcher—namely, which test should be used for a particular association study which may consider multiple genes and/or multiple phenotypes. Recently, tests have been proposed which combine individual tests to minimize power loss while improving the robustness to a wide range of genetic architectures. In our analysis, we propose an expansion of these approaches, by providing a general method that works for combining an arbitrarily large number of any gene‐based rare variant test—a flexibility typically not available in other combined testing methods. We provide a theoretical framework for evaluating our combined test to provide direct insights into the relationship between test‐test correlation, test power and the combined test power relative to individual testing approaches and other combined testing approaches. We demonstrate that our flexible combined testing method can provide improved power and robustness against a wide range of genetic architectures. We further demonstrate the performance of our combined test on simulated genotypes, as well as on a dataset of real genotypes with simulated phenotypes. We support the increased use of flexible combined tests in practice to maximize robustness of rare‐variant testing strategies against a wide‐range of genetic architectures

UMP, SCUN seal bilateral collaborations

Universiti Malaysia Pahang (UMP) fortified its international networking when it sealed a Memorandum of Understanding (MoU) with China’s South-Central University for Nationalities (SCUN) in Beijing, China, on June 1

General Approach for Combining Diverse Rare Variant Association Tests Provides Improved Robustness Across a Wider Range of Genetic Architectures

The widespread availability of genome sequencing data made possible by way of next-generation technologies has yielded a flood of different gene-based rare variant association tests. Most of these tests have been published because they have superior power for particular genetic architectures. However, for applied researchers it is challenging to know which test to choose in practice when little is known a priori about genetic architecture. Recently, tests have been proposed which combine two particular individual tests (one burden and one variance components) to minimize power loss while improving robustness to a wider range of genetic architectures. In our analysis we propose an expansion of these approaches, yielding a general method that works for combining any number of individual tests. We demonstrate that running multiple different tests on the same dataset and using a Bonferroni correction for multiple testing is never better than combining tests using our general method. We also find that using a test statistic that is highly robust to the inclusion of non-causal variants (Joint-infinity) together with a previously published combined test (SKAT-O) provides improved robustness to a wide range of genetic architectures and should be considered for use in practice. Software for this approach is supplied. We support the increased use of combined tests in practice-- as well as further exploration of novel combined testing approaches using the general framework provided here--to maximize robustness of rare-variant testing strategies against a wide range of genetic architectures

Admixture mapping implicates 13q33.3 as ancestry-of-origin locus for Alzheimer disease in Hispanic and Latino populations

Alzheimer disease (AD) is the most common form of senile dementia, with high incidence late in life in many populations including Caribbean Hispanic (CH) populations. Such admixed populations, descended from more than one ancestral population, can present challenges for genetic studies, including limited sample sizes and unique analytical constraints. Therefore, CH populations and other admixed populations have not been well represented in studies of AD, and much of the genetic variation contributing to AD risk in these populations remains unknown. Here, we conduct genome-wide analysis of AD in multiplex CH families from the Alzheimer Disease Sequencing Project (ADSP). We developed, validated, and applied an implementation of a logistic mixed model for admixture mapping with binary traits that leverages genetic ancestry to identify ancestry-of-origin loci contributing to AD. We identified three loci on chromosome 13q33.3 associated with reduced risk of AD, where associations were driven by Native American (NAM) ancestry. This AD admixture mapping signal spans the FAM155A, ABHD13, TNFSF13B, LIG4, and MYO16 genes and was supported by evidence for association in an independent sample from the Alzheimer's Genetics in Argentina—Alzheimer Argentina consortium (AGA-ALZAR) study with considerable NAM ancestry. We also provide evidence of NAM haplotypes and key variants within 13q33.3 that segregate with AD in the ADSP whole-genome sequencing data. Interestingly, the widely used genome-wide association study approach failed to identify associations in this region. Our findings underscore the potential of leveraging genetic ancestry diversity in recently admixed populations to improve genetic mapping, in this case for AD-relevant loci.Fil: Horimoto, Andrea R.V.R.. University of Washington; Estados UnidosFil: Boyken, Lisa A.. University of Washington; Estados UnidosFil: Blue, Elizabeth E.. University of Washington; Estados Unidos. Brotman Baty Institute for Precision Medicine; Estados UnidosFil: Grinde, Kelsey E.. University of Washington; Estados Unidos. Macalester College; Estados UnidosFil: Nafikov, Rafael A.. University of Washington; Estados UnidosFil: Sohi, Harkirat K.. University of Washington; Estados UnidosFil: Nato, Alejandro Q.. University of Washington; Estados Unidos. Marshall University; Estados UnidosFil: Bis, Joshua C.. University of Washington; Estados UnidosFil: Brusco, Luis Ignacio. Universidad de Buenos Aires. Facultad de Medicina; Argentina. Consejo Nacional de Investigaciones Científicas y Técnicas; ArgentinaFil: Morelli, Laura. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Parque Centenario. Instituto de Investigaciones Bioquímicas de Buenos Aires. Fundación Instituto Leloir. Instituto de Investigaciones Bioquímicas de Buenos Aires; ArgentinaFil: Ramirez, Alfredo Jose. University Of Cologne; Alemania. Universitat Bonn; Alemania. German Center for Neurodegenerative Diseases; Alemania. University Of Texas Health Science Center At San Antonio (ut Health San Antonio) ; University Of Texas At San Antonio; . Universidad Nacional Arturo Jauretche. Unidad Ejecutora de Estudios en Neurociencias y Sistemas Complejos. Provincia de Buenos Aires. Ministerio de Salud. Hospital Alta Complejidad en Red El Cruce Dr. Néstor Carlos Kirchner Samic. Unidad Ejecutora de Estudios en Neurociencias y Sistemas Complejos. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - La Plata. Unidad Ejecutora de Estudios en Neurociencias y Sistemas Complejos; ArgentinaFil: Dalmasso, Maria Carolina. Universidad Nacional Arturo Jauretche. Unidad Ejecutora de Estudios en Neurociencias y Sistemas Complejos. Provincia de Buenos Aires. Ministerio de Salud. Hospital Alta Complejidad en Red El Cruce Dr. Néstor Carlos Kirchner Samic. Unidad Ejecutora de Estudios en Neurociencias y Sistemas Complejos. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - La Plata. Unidad Ejecutora de Estudios en Neurociencias y Sistemas Complejos; Argentina. University Of Cologne; AlemaniaFil: Temple, Seth. University of Washington; Estados UnidosFil: Satizabal, Claudia. University Of Texas Health Science Center At San Antonio (ut Health San Antonio) ; University Of Texas At San Antonio; . University of Texas at San Antonio; Estados UnidosFil: Browning, Sharon R.. University of Washington; Estados UnidosFil: Seshadri, Sudha. University Of Texas Health Science Center At San Antonio (ut Health San Antonio) ; University Of Texas At San Antonio; . University of Texas at San Antonio; Estados UnidosFil: Wijsman, Ellen M.. University of Washington; Estados UnidosFil: Thornton, Timothy A.. University of Washington; Estados Unido

Genome-wide analysis of dental caries and periodontitis combining clinical and self-reported data

Dental caries and periodontitis account for a vast burden of morbidity and healthcare spending, yet their genetic basis remains largely uncharacterized. Here, we identify self-reported dental disease proxies which have similar underlying genetic contributions to clinical disease measures and then combine these in a genome-wide association study meta-analysis, identifying 47 novel and conditionally-independent risk loci for dental caries. We show that the heritability of dental caries is enriched for conserved genomic regions and partially overlapping with a range of complex traits including smoking, education, personality traits and metabolic measures. Using cardio-metabolic traits as an example in Mendelian randomization analysis, we estimate causal relationships and provide evidence suggesting that the processes contributing to dental caries may have undesirable downstream effects on health

Statistical Inference in Admixed Populations

Thesis (Ph.D.)--University of Washington, 2019Understanding the genetic causes of human diseases and traits has long been of interest in the scientific community. However, the large majority of research in this area has been conducted in European populations. This dissertation focuses on developing statistical methods for genetic studies in admixed populations, such as African Americans and Hispanics/Latinos, that have been historically underrepresented in genetics research. The diverse, mixed ancestry of admixed populations presents unique opportunities for statistical inference, many of which are explored in this work. Here, we focus in particular on two important tasks: inferring genetic ancestry from genotype and sequence data, and identifying genetic variants associated with complex traits and diseases. We propose and evaluate methods for inferring local ancestry on chromosome X, correcting for multiple testing in genome-wide admixture mapping studies, and controlling for confounding by global ancestry in admixture mapping and genome-wide association studies in admixed populations. We motivate our proposed methods with theoretical results, simulation studies, and applications to genotype and whole genome sequence data from large studies of African American and Hispanic/Latino individuals. Our work provides solutions to a number of the statistical challenges posed by genetic studies in admixed populations, and we hope that our results will help guide future studies in these populations

Circular RNA Profiles in Viremia and ART Suppression Predict Competing circRNA–miRNA–mRNA Networks Exclusive to HIV-1 Viremic Patients

Since the onset of the HIV-1/AIDS epidemic in 1981, 75 million people have been infected with the virus, and the disease remains a public health crisis worldwide. Circular RNAs (circRNAs) are derived from excised exons and introns during backsplicing, a form of alternative splicing. The relevance of unconventional, non-capped, and non-poly(A) transcripts to transcriptomics studies remains to be routinely investigated. Knowledge gaps to be filled are the interface between host-encoded circRNAs and viral replication in chronically progressed patients and upon treatment with antiviral drugs. We implemented a bioinformatic pipeline and repurpose publicly archived RNA sequence reads from the blood of 19 HIV-1-positive patients that previously compared transcriptomes during viremia and viremia suppression by antiretroviral therapy (ART). The in silico analysis identified viremic patients’ circRNA that became undetectable after ART. The circRNAs originated from a subset of host genes enriched in the HDAC biological pathway. These circRNAs and parental mRNAs held in common a small collection of miRNA response elements (MREs), some of which were present in HIV-1 mRNAs. The function of the MRE-containing target mRNA enriched the RNA polymerase II GO pathway. To visualize the interplay between individual circRNA–miRNA–target mRNA, important for HIV-1 and potentially other diseases, an Interactive Circos tool was developed to efficiently parse the intricately competing endogenous network of circRNA–miRNA–mRNA interactions originating from seven circRNA singled out in viremic versus non-viremic patients. The combined downregulation of the identified circRNAs warrants investigation as a novel antiviral targeting strategy

A multistep approach to single nucleotide polymorphism–set analysis: an evaluation of power and type I error of gene-based tests of association after pathway-based association tests

The aggregation of functionally associated variants given a priori biological information can aid in the discovery of rare variants associated with complex diseases. Many methods exist that aggregate rare variants into a set and compute a single p value summarizing association between the set of rare variants and a phenotype of interest. These methods are often called gene-based, rare variant tests of association because the variants in the set are often all contained within the same gene. A reasonable extension of these approaches involves aggregating variants across an even larger set of variants (eg, all variants contained in genes within a pathway). Testing sets of variants such as pathways for association with a disease phenotype reduces multiple testing penalties, may increase power, and allows for straightforward biological interpretation. However, a significant variant-set association test does not indicate precisely which variants contained within that set are causal. Because pathways often contain many variants, it may be helpful to follow-up significant pathway tests by conducting gene-based tests on each gene in that pathway to narrow in on the region of causal variants. In this paper, we propose such a multistep approach for variant-set analysis that can also account for covariates and complex pedigree structure. We demonstrate this approach on simulated phenotypes from Genetic Analysis Workshop 19. We find generally better power for the multistep approach when compared to a more conventional, single-step approach that simply runs gene-based tests of association on each gene across the genome. Further work is necessary to evaluate the multistep approach on different data sets with different characteristics

Circular RNA Profiles in Viremia and ART Suppression Predict Competing circRNA–miRNA–mRNA Networks Exclusive to HIV-1 Viremic Patients

Since the onset of the HIV-1/AIDS epidemic in 1981, 75 million people have been infected with the virus, and the disease remains a public health crisis worldwide. Circular RNAs (circRNAs) are derived from excised exons and introns during backsplicing, a form of alternative splicing. The relevance of unconventional, non-capped, and non-poly(A) transcripts to transcriptomics studies remains to be routinely investigated. Knowledge gaps to be filled are the interface between host-encoded circRNAs and viral replication in chronically progressed patients and upon treatment with antiviral drugs. We implemented a bioinformatic pipeline and repurpose publicly archived RNA sequence reads from the blood of 19 HIV-1-positive patients that previously compared transcriptomes during viremia and viremia suppression by antiretroviral therapy (ART). The in silico analysis identified viremic patients’ circRNA that became undetectable after ART. The circRNAs originated from a subset of host genes enriched in the HDAC biological pathway. These circRNAs and parental mRNAs held in common a small collection of miRNA response elements (MREs), some of which were present in HIV-1 mRNAs. The function of the MRE-containing target mRNA enriched the RNA polymerase II GO pathway. To visualize the interplay between individual circRNA–miRNA–target mRNA, important for HIV-1 and potentially other diseases, an Interactive Circos tool was developed to efficiently parse the intricately competing endogenous network of circRNA–miRNA–mRNA interactions originating from seven circRNA singled out in viremic versus non-viremic patients. The combined downregulation of the identified circRNAs warrants investigation as a novel antiviral targeting strategy