Meta-analysis of Gene-Level Associations for Rare Variants Based on Single-Variant Statistics

Meta-analysis of genome-wide association studies (GWASs) has led to the discoveries of many common variants associated with complex human diseases. There is a growing recognition that identifying “causal” rare variants also requires large-scale meta-analysis. The fact that association tests with rare variants are performed at the gene level rather than at the variant level poses unprecedented challenges in the meta-analysis. First, different studies may adopt different gene-level tests, so the results are not compatible. Second, gene-level tests require multivariate statistics (i.e., components of the test statistic and their covariance matrix), which are difficult to obtain. To overcome these challenges, we propose to perform gene-level tests for rare variants by combining the results of single-variant analysis (i.e., p values of association tests and effect estimates) from participating studies. This simple strategy is possible because of an insight that multivariate statistics can be recovered from single-variant statistics, together with the correlation matrix of the single-variant test statistics, which can be estimated from one of the participating studies or from a publicly available database. We show both theoretically and numerically that the proposed meta-analysis approach provides accurate control of the type I error and is as powerful as joint analysis of individual participant data. This approach accommodates any disease phenotype and any study design and produces all commonly used gene-level tests. An application to the GWAS summary results of the Genetic Investigation of ANthropometric Traits (GIANT) consortium reveals rare and low-frequency variants associated with human height. The relevant software is freely available.     

Improved Detection of Common Variants Associated with Schizophrenia by Leveraging Pleiotropy with Cardiovascular-Disease Risk Factors

Several lines of evidence suggest that genome-wide association studies (GWASs) have the potential to explain more of the “missing heritability” of common complex phenotypes. However, reliable methods for identifying a larger proportion of SNPs are currently lacking. Here, we present a genetic-pleiotropy-informed method for improving gene discovery with the use of GWAS summary-statistics data. We applied this methodology to identify additional loci associated with schizophrenia (SCZ), a highly heritable disorder with significant missing heritability. Epidemiological and clinical studies suggest comorbidity between SCZ and cardiovascular-disease (CVD) risk factors, including systolic blood pressure, triglycerides, low- and high-density lipoprotein, body mass index, waist-to-hip ratio, and type 2 diabetes. Using stratified quantile-quantile plots, we show enrichment of SNPs associated with SCZ as a function of the association with several CVD risk factors and a corresponding reduction in false discovery rate (FDR). We validate this “pleiotropic enrichment” by demonstrating increased replication rate across independent SCZ substudies. Applying the stratified FDR method, we identified 25 loci associated with SCZ at a conditional FDR level of 0.01. Of these, ten loci are associated with both SCZ and CVD risk factors, mainly triglycerides and low- and high-density lipoproteins but also waist-to-hip ratio, systolic blood pressure, and body mass index. Together, these findings suggest the feasibility of using genetic-pleiotropy-informed methods for improving gene discovery in SCZ and identifying potential mechanistic relationships with various CVD risk factors.     

Epistatic Gene-Based Interaction Analyses for Glaucoma in eMERGE and NEIGHBOR Consortium

Primary open angle glaucoma (POAG) is a complex disease and is one of the major leading causes of blindness worldwide. Genome-wide association studies have successfully identified several common variants associated with glaucoma; however, most of these variants only explain a small proportion of the genetic risk. Apart from the standard approach to identify main effects of variants across the genome, it is believed that gene-gene interactions can help elucidate part of the missing heritability by allowing for the test of interactions between genetic variants to mimic the complex nature of biology. To explain the etiology of glaucoma, we first performed a genome-wide association study (GWAS) on glaucoma case-control samples obtained from electronic medical records (EMR) to establish the utility of EMR data in detecting non-spurious and relevant associations; this analysis was aimed at confirming already known associations with glaucoma and validating the EMR derived glaucoma phenotype. Our findings from GWAS suggest consistent evidence of several known associations in POAG. We then performed an interaction analysis for variants found to be marginally associated with glaucoma (SNPs with main effect p-value <0.01) and observed interesting findings in the electronic MEdical Records and GEnomics Network (eMERGE) network dataset. Genes from the top epistatic interactions from eMERGE data (Likelihood Ratio Test i.e. LRT p-value <1e-05) were then tested for replication in the NEIGHBOR consortium dataset. To replicate our findings, we performed a gene-based SNP-SNP interaction analysis in NEIGHBOR and observed significant gene-gene interactions (p-value <0.001) among the top 17 gene-gene models identified in the discovery phase. Variants from gene-gene interaction analysis that we found to be associated with POAG explain 3.5% of additional genetic variance in eMERGE dataset above what is explained by the SNPs in genes that are replicated from previous GWAS studies (which was only 2.1% variance explained in eMERGE dataset); in the NEIGHBOR dataset, adding replicated SNPs from gene-gene interaction analysis explain 3.4% of total variance whereas GWAS SNPs alone explain only 2.8% of variance. Exploring gene-gene interactions may provide additional insights into many complex traits when explored in properly designed and powered association studies.     

Structure of the scientific community modelling the evolution of resistance

Faced with the recurrent evolution of resistance to pesticides and drugs, the scientific community has developed theoretical models aimed at identifying the main factors of this evolution and predicting the efficiency of resistance management strategies. The evolutionary forces considered by these models are generally similar for viruses, bacteria, fungi, plants or arthropods facing drugs or pesticides, so interaction between scientists working on different biological organisms would be expected. We tested this by analysing co-authorship and co-citation networks using a database of 187 articles published from 1977 to 2006 concerning models of resistance evolution to all major classes of pesticides and drugs. These analyses identified two main groups. One group, led by ecologists or agronomists, is interested in agricultural crop or stock pests and diseases. It mainly uses a population genetics approach to model the evolution of resistance to insecticidal proteins, insecticides, herbicides, antihelminthic drugs and miticides. By contrast, the other group, led by medical scientists, is interested in human parasites and mostly uses epidemiological models to study the evolution of resistance to antibiotic and antiviral drugs. Our analyses suggested that there is also a small scientific group focusing on resistance to antimalaria drugs, and which is only poorly connected with the two larger groups. The analysis of cited references indicates that each of the two large communities publishes its research in a different set of literature and has its own keystone references: citations with a large impact in one group are almost never cited by the other. We fear the lack of exchange between the two communities might slow progress concerning resistance evolution which is currently a major issue for society.     

The Genographic Project public participation mitochondrial DNA database

The Genographic Project is studying the genetic signatures of ancient human migrations and creating an open-source research database. It allows members of the public to participate in a real-time anthropological genetics study by submitting personal samples for analysis and donating the genetic results to the database. We report our experience from the first 18 months of public participation in the Genographic Project, during which we have created the largest standardized human mitochondrial DNA (mtDNA) database ever collected, comprising 78,590 genotypes. Here, we detail our genotyping and quality assurance protocols including direct sequencing of the mtDNA HVS-I, genotyping of 22 coding-region SNPs, and a series of computational quality checks based on phylogenetic principles. This database is very informative with respect to mtDNA phylogeny and mutational dynamics, and its size allows us to develop a nearest neighbor-based methodology for mtDNA haplogroup prediction based on HVS-I motifs that is superior to classic rule-based approaches. We make available to the scientific community and general public two new resources: a periodically updated database comprising all data donated by participants, and the nearest neighbor haplogroup prediction tool.     

In situ proteolysis for protein crystallization and structure determination

We tested the general applicability of in situ proteolysis to form protein crystals suitable for structure determination by adding a protease (chymotrypsin or trypsin) digestion step to crystallization trials of 55 bacterial and 14 human proteins that had proven recalcitrant to our best efforts at crystallization or structure determination. This is a work in progress; so far we determined structures of 9 bacterial proteins and the human aminoimidazole ribonucleotide synthetase (AIRS) domain.     

Proteome-wide Mendelian randomization identifies causal links between blood proteins and severe COVID-19

In November 2021, the COVID-19 pandemic death toll surpassed five million individuals. We applied Mendelian randomization including >3,000 blood proteins as exposures to identify potential biomarkers that may indicate risk for hospitalization or need for respiratory support or death due to COVID-19, respectively. After multiple testing correction, using genetic instruments and under the assumptions of Mendelian Randomization, our results were consistent with higher blood levels of five proteins GCNT4, CD207, RAB14, C1GALT1C1, and ABO being causally associated with an increased risk of hospitalization or respiratory support/death due to COVID-19 (ORs = 1.12–1.35). Higher levels of FAAH2 were solely associated with an increased risk of hospitalization (OR = 1.19). On the contrary, higher levels of SELL, SELE, and PECAM-1 decrease risk of hospitalization or need for respiratory support/death (ORs = 0.80–0.91). Higher levels of LCTL, SFTPD, KEL, and ATP2A3 were solely associated with a decreased risk of hospitalization (ORs = 0.86–0.93), whilst higher levels of ICAM-1 were solely associated with a decreased risk of respiratory support/death of COVID-19 (OR = 0.84). Our findings implicate blood group markers and binding proteins in both hospitalization and need for respiratory support/death. They, additionally, suggest that higher levels of endocannabinoid enzymes may increase the risk of hospitalization. Our research replicates findings of blood markers previously associated with COVID-19 and prioritises additional blood markers for risk prediction of severe forms of COVID-19. Furthermore, we pinpoint druggable targets potentially implicated in disease pathology.     

An ancient haplotype containing antimicrobial peptide gene variants is associated with severe fungal skin disease in Persian cats

Dermatophytosis, also known as ringworm, is a contagious fungal skin disease affecting humans and animals worldwide. Persian cats exhibit severe forms of the disease more commonly than other breeds of cat, including other long-haired breeds. Certain types of severe dermatophytosis in humans are reportedly caused by monogenic inborn errors of immunity. The goal of this study was to identify genetic variants in Persian cats contributing to the phenotype of severe dermatophytosis. Whole-genome sequencing of case and control Persian cats followed by a genome-wide association study identified a highly divergent, disease-associated haplotype on chromosome F1 containing the S100 family of genes. S100 calcium binding protein A9 (S100A9), which encodes a subunit of the antimicrobial heterodimer known as calprotectin, contained 13 nonsynonymous variants between cases and controls. Evolutionary analysis of S100A9 haplotypes comparing cases, controls, and wild felids suggested the divergent disease-associated haplotype was likely introgressed into the domestic cat lineage and maintained via balancing selection. We demonstrated marked upregulation of calprotectin expression in the feline epidermis during dermatophytosis, suggesting involvement in disease pathogenesis. Given this divergent allele has been maintained in domestic cat and wildcat populations, this haplotype may have beneficial effects against other pathogens. The pathogen specificity of this altered protein should be investigated before attempting to reduce the allele frequency in the Persian cat breed. Further work is needed to clarify if severe Persian dermatophytosis is a monogenic disease or if hidden disease-susceptibility loci remain to be discovered. Consideration should be given to engineering antimicrobial peptides such as calprotectin for topical treatment of dermatophytosis in humans and animals.