Statistical analysis strategies for association studies involving rare variants

Genomic analysis reveals the association of KIT and MITF variants with the white spotting in swamp buffaloes

BACKGROUND

Swamp-type buffaloes with varying degrees of white spotting are found exclusively in Tana Toraja, South Sulawesi, Indonesia, where spotted buffalo bulls are highly valued in accordance with the Torajan customs. The white spotting depigmentation is caused by the absence of melanocytes. However, the genetic variants that cause this phenotype have not been fully characterized. The objective of this study was to identify the genomic regions and variants responsible for this unique coat-color pattern.

RESULTS

Genome-wide association study (GWAS) and selection signature analysis identified MITF as a key gene based on the whole-genome sequencing data of 28 solid and 39 spotted buffaloes, while KIT was also found to be involved in the development of this phenotype by a candidate gene approach. Alternative candidate mutations included, in addition to the previously reported nonsense mutation c.649 C > T (p.Arg217*) and splice donor mutation c.1179 + 2T > A in MITF, a nonsense mutation c.2028T > A (p.Tyr676*) in KIT. All these three mutations were located in the genomic regions that were highly conserved exclusively in Indonesian swamp buffaloes and they accounted largely (95%) for the manifestation of white spotting. Last but not the least, ADAMTS20 and TWIST2 may also contribute to the diversification of this coat-color pattern.

CONCLUSIONS

The alternative mutations identified in this study affect, at least partially and independently, the development of melanocytes. The presence and persistence of such mutations may be explained by significant financial and social value of spotted buffaloes used in historical Rambu Solo ceremony in Tana Toraja, Indonesia. Several de novo spontaneous mutations have therefore been favored by traditional breeding for the spotted buffaloes.

Biostatistical Aspects of Whole Genome Sequencing Studies: Preprocessing and Quality Control

Rapid advances in high-throughput DNA sequencing technologies have enabled large-scale whole genome sequencing (WGS) studies. Before performing association analysis between phenotypes and genotypes, preprocessing and quality control (QC) of the raw sequence data need to be performed. Because many biostatisticians have not been working with WGS data so far, we first sketch Illumina's short-read sequencing technology. Second, we explain the general preprocessing pipeline for WGS studies. Third, we provide an overview of important QC metrics, which are applied to WGS data: on the raw data, after mapping and alignment, after variant calling, and after multisample variant calling. Fourth, we illustrate the QC with the data from the GENEtic SequencIng Study Hamburg-Davos (GENESIS-HD), a study involving more than 9000 human whole genomes. All samples were sequenced on an Illumina NovaSeq 6000 with an average coverage of 35× using a PCR-free protocol. For QC, one genome in a bottle (GIAB) trio was sequenced in four replicates, and one GIAB sample was successfully sequenced 70 times in different runs. Fifth, we provide empirical data on the compression of raw data using the DRAGEN original read archive (ORA). The most important quality metrics in the application were genetic similarity, sample cross-contamination, deviations from the expected Het/Hom ratio, relatedness, and coverage. The compression ratio of the raw files using DRAGEN ORA was 5.6:1, and compression time was linear by genome coverage. In summary, the preprocessing, joint calling, and QC of large WGS studies are feasible within a reasonable time, and efficient QC procedures are readily available.

Genetic modifiers of cognitive decline in PSEN1 E280A Alzheimer's disease

INTRODUCTION

Rate of cognitive decline (RCD) in Alzheimer's disease (AD) determines the degree of impairment for patients and of burden for caretakers. We studied the association of RCD with genetic variants in AD.

METHODS

RCD was evaluated in 62 familial AD (FAD) and 53 sporadic AD (SAD) cases, and analyzed by whole-exome sequencing for association with common exonic functional variants. Findings were validated in post mortem brain tissue.

RESULTS

One hundred seventy-two gene variants in FAD, and 227 gene variants in SAD associated with RCD. In FAD, performance decline of the immediate recall of the Rey-Osterrieth figure test associated with 122 genetic variants. Olfactory receptor OR51B6 showed the highest number of associated variants. Its expression was detected in temporal cortex neurons.

DISCUSSION

Impaired olfactory function has been associated with cognitive impairment in AD. Genetic variants in these or other genes could help to identify risk of faster memory decline in FAD and SAD patients.

DYNATE: Localizing rare-variant association regions via multiple testing embedded in an aggregation tree

Rare-variants (RVs) genetic association studies enable researchers to uncover the variation in phenotypic traits left unexplained by common variation. Traditional single-variant analysis lacks power; thus, researchers have developed various methods to aggregate the effects of RVs across genomic regions to study their collective impact. Some existing methods utilize a static delineation of genomic regions, often resulting in suboptimal effect aggregation, as neutral subregions within the test region will result in an attenuation of signal. Other methods use varying windows to search for signals but often result in long regions containing many neutral RVs. To pinpoint short genomic regions enriched for disease-associated RVs, we developed a novel method, DYNamic Aggregation TEsting (DYNATE). DYNATE dynamically and hierarchically aggregates smaller genomic regions into larger ones and performs multiple testing for disease associations with a controlled weighted false discovery rate. DYNATE's main advantage lies in its strong ability to identify short genomic regions highly enriched for disease-associated RVs. Extensive numerical simulations demonstrate the superior performance of DYNATE under various scenarios compared with existing methods. We applied DYNATE to an amyotrophic lateral sclerosis study and identified a new gene, EPG5, harboring possibly pathogenic mutations.

Rare variant association on unrelated individuals in case-control studies using aggregation tests: existing methods and current limitations

Over the past years, progress made in next-generation sequencing technologies and bioinformatics have sparked a surge in association studies. Especially, genome-wide association studies (GWASs) have demonstrated their effectiveness in identifying disease associations with common genetic variants. Yet, rare variants can contribute to additional disease risk or trait heterogeneity. Because GWASs are underpowered for detecting association with such variants, numerous statistical methods have been recently proposed. Aggregation tests collapse multiple rare variants within a genetic region (e.g. gene, gene set, genomic loci) to test for association. An increasing number of studies using such methods successfully identified trait-associated rare variants and led to a better understanding of the underlying disease mechanism. In this review, we compare existing aggregation tests, their statistical features and scope of application, splitting them into the five classical classes: burden, adaptive burden, variance-component, omnibus and other. Finally, we describe some limitations of current aggregation tests, highlighting potential direction for further investigations.

A genotype-to-phenotype approach suggests under-reporting of single nucleotide variants in nephrocystin-1 (NPHP1) related disease (UK 100,000 Genomes Project)

Autosomal recessive whole gene deletions of nephrocystin-1 (NPHP1) result in abnormal structure and function of the primary cilia. These deletions can result in a tubulointerstitial kidney disease known as nephronophthisis and retinal (Senior-Løken syndrome) and neurological (Joubert syndrome) diseases. Nephronophthisis is a common cause of end-stage kidney disease (ESKD) in children and up to 1% of adult onset ESKD. Single nucleotide variants (SNVs) and small insertions and deletions (Indels) have been less well characterised. We used a gene pathogenicity scoring system (GenePy) and a genotype-to-phenotype approach on individuals recruited to the UK Genomics England (GEL) 100,000 Genomes Project (100kGP) (n = 78,050). This approach identified all participants with NPHP1-related diseases reported by NHS Genomics Medical Centres and an additional eight participants. Extreme NPHP1 gene scores, often underpinned by clear recessive inheritance, were observed in patients from diverse recruitment categories, including cancer, suggesting the possibility of a more widespread disease than previously appreciated. In total, ten participants had homozygous CNV deletions with eight homozygous or compound heterozygous with SNVs. Our data also reveals strong in-silico evidence that approximately 44% of NPHP1 related disease may be due to SNVs with AlphaFold structural modelling evidence for a significant impact on protein structure. This study suggests historical under-reporting of SNVS in NPHP1 related diseases compared with CNVs.

Relationship between the TGFBR1 Gene and Molar Incisor Hypomineralization

BACKGROUND

Molar Incisor Hypomineralization Syndrome (MIH) is a problem of increasing incidence that represents a new challenge in the dental treatment of many of the children we see in our dental offices. Understanding the etiology of this syndrome (still unknown) will help us to prevent the appearance of this process. Lately a certain genetic relationship has been suggested in the syndrome. The aim of the present study was to explore the relationship between activation of the TGFBR1 gene and the development of MIH, as recent studies suggest that there may be an association in this regard.

MATERIALS AND METHODS

The study sample consisted of 50 children between 6-17 years of age with MIH, each with at least one parent and a sibling with or without MIH, and a group control of 100 children without MIH. The condition of the permanent molars and incisors was evaluated and recorded based on the criteria of Mathu-Muju and Wright. Saliva samples were collected after washing and rinsing of the oral cavity. Genotyping was performed with the saliva samples for the selection of a target polymorphism of the studied gene (TGFBR1).

RESULTS

The mean age was 9.7 years (SD 2.36). Of the 50 children with MIH, 56% were boys and 44% girls. The degree of MIH was predominantly severe (58%), with moderate and mild involvement in 22% and 20% of the cases, respectively, according to the classification of Mathu-Muju. The allelic frequencies were seen to behave as expected. The logistic regression analysis aimed to relate each polymorphism to the presence or absence of the factors. These results were inconclusive, with no evidence suggesting an alteration of the TGFBR1 gene to be related to the appearance of MIH.

CONCLUSIONS

Within the limitations posed by a study of these characteristics, it can be affirmed that no relationship has been found between the TGFBR1 gene and the appearance of molar incisor hypomineralization.

Common and rare variant associations with latent traits underlying depression, bipolar disorder, and schizophrenia

Genetic studies in psychiatry have primarily focused on the effects of common genetic variants, but few have investigated the role of rare genetic variants, particularly for major depression. In order to explore the role of rare variants in the gap between estimates of single nucleotide polymorphism (SNP) heritability and twin study heritability, we examined the contribution of common and rare genetic variants to latent traits underlying psychiatric disorders using high-quality imputed genotype data from the UK Biobank. Using a pre-registered analysis, we used items from the UK Biobank Mental Health Questionnaire relevant to three psychiatric disorders: major depression (N = 134,463), bipolar disorder (N = 117,376) and schizophrenia (N = 130,013) and identified a general hierarchical factor for each that described participants' responses. We calculated participants' scores on these latent traits and conducted single-variant genetic association testing (MAF > 0.05%), gene-based burden testing and pathway association testing associations with these latent traits. We tested for enrichment of rare variants (MAF 0.05-1%) in genes that had been previously identified by common variant genome-wide association studies, and genes previously associated with Mendelian disorders having relevant symptoms. We found moderate genetic correlations between the latent traits in our study and case-control phenotypes in previous genome-wide association studies, and identified one common genetic variant (rs72657988, minor allele frequency = 8.23%, p = 1.01 × 10^-9) associated with the general factor of schizophrenia, but no other single variants, genes or pathways passed significance thresholds in this analysis, and we did not find enrichment in previously identified genes.

A fast and efficient approach for gene-based association studies of ordinal phenotypes

Many human disease conditions need to be measured by ordinal phenotypes, so analysis of ordinal phenotypes is valuable in genome-wide association studies (GWAS). However, existing association methods for dichotomous or quantitative phenotypes are not appropriate to ordinal phenotypes. Therefore, based on an aggregated Cauchy association test, we propose a fast and efficient association method to test the association between genetic variants and an ordinal phenotype. To enrich association signals of rare variants, we first use the burden method to aggregate rare variants. Then we respectively test the significance of the aggregated rare variants and other common variants. Finally, the combination of transformed variant-level P values is taken as test statistic, that approximately follows Cauchy distribution under the null hypothesis. Extensive simulation studies and analysis of GAW19 show that our proposed method is powerful and computationally fast as a gene-based method. Especially, in the presence of an extremely low proportion of causal variants in a gene, our method has better performance.

A framework for detecting noncoding rare-variant associations of large-scale whole-genome sequencing studies

Large-scale whole-genome sequencing studies have enabled analysis of noncoding rare-variant (RV) associations with complex human diseases and traits. Variant-set analysis is a powerful approach to study RV association. However, existing methods have limited ability in analyzing the noncoding genome. We propose a computationally efficient and robust noncoding RV association detection framework, STAARpipeline, to automatically annotate a whole-genome sequencing study and perform flexible noncoding RV association analysis, including gene-centric analysis and fixed window-based and dynamic window-based non-gene-centric analysis by incorporating variant functional annotations. In gene-centric analysis, STAARpipeline uses STAAR to group noncoding variants based on functional categories of genes and incorporate multiple functional annotations. In non-gene-centric analysis, STAARpipeline uses SCANG-STAAR to incorporate dynamic window sizes and multiple functional annotations. We apply STAARpipeline to identify noncoding RV sets associated with four lipid traits in 21,015 discovery samples from the Trans-Omics for Precision Medicine (TOPMed) program and replicate several of them in an additional 9,123 TOPMed samples. We also analyze five non-lipid TOPMed traits.

Genetic and functional analyses implicate microRNA 499A in bipolar disorder development

Bipolar disorder (BD) is a complex mood disorder with a strong genetic component. Recent studies suggest that microRNAs contribute to psychiatric disorder development. In BD, specific candidate microRNAs have been implicated, in particular miR-137, miR-499a, miR-708, miR-1908 and miR-2113. The aim of the present study was to determine the contribution of these five microRNAs to BD development. For this purpose, we performed: (i) gene-based tests of the five microRNA coding genes, using data from a large genome-wide association study of BD; (ii) gene-set analyses of predicted, brain-expressed target genes of the five microRNAs; (iii) resequencing of the five microRNA coding genes in 960 BD patients and 960 controls and (iv) in silico and functional studies for selected variants. Gene-based tests revealed a significant association with BD for MIR499A, MIR708, MIR1908 and MIR2113. Gene-set analyses revealed a significant enrichment of BD associations in the brain-expressed target genes of miR-137 and miR-499a-5p. Resequencing identified 32 distinct rare variants (minor allele frequency < 1%), all of which showed a non-significant numerical overrepresentation in BD patients compared to controls (p = 0.214). Seven rare variants were identified in the predicted stem-loop sequences of MIR499A and MIR2113. These included rs142927919 in MIR2113 (p_nom = 0.331) and rs140486571 in MIR499A (p_nom = 0.297). In silico analyses predicted that rs140486571 might alter the miR-499a secondary structure. Functional analyses showed that rs140486571 significantly affects miR-499a processing and expression. Our results suggest that MIR499A dysregulation might contribute to BD development. Further research is warranted to elucidate the contribution of the MIR499A regulated network to BD susceptibility.

The Genetic Architecture of Non-Syndromic Rhegmatogenous Retinal Detachment

Rhegmatogenous retinal detachment (RRD) is the most common form of retinal detachment (RD), affecting 1 in 10,000 patients per year. The condition has significant ocular morbidity, with a sizeable proportion of patients obtaining poor visual outcomes. Despite this, the genetics underpinning Idiopathic Retinal Detachment (IRD) remain poorly understood; this is likely due to small sample sizes in relevant studies. The majority of research pertains to the well-characterised Mende lian syndromes, such as Sticklers and Wagners, associated with RRD. Nevertheless, in recent years, there has been an increasing body of literature identifying the common genetic mutations and mechanisms associated with IRD. Several recent Genomic Wide Association Studies (GWAS) studies have identified a number of genetic loci related to the development of IRD. Our review aims to provide an up-to-date summary of the significant genetic mechanisms and associations of Idiopathic RRD.

Solid Waste Image Classification Using Deep Convolutional Neural Network

Separating household waste into categories such as organic and recyclable is a critical part of waste management systems to make sure that valuable materials are recycled and utilised. This is beneficial to human health and the environment because less risky treatments are used at landfill and/or incineration, ultimately leading to improved circular economy. Conventional waste separation relies heavily on manual separation of objects by humans, which is inefficient, expensive, time consuming, and prone to subjective errors caused by limited knowledge of waste classification. However, advances in artificial intelligence research has led to the adoption of machine learning algorithms to improve the accuracy of waste classification from images. In this paper, we used a waste classification dataset to evaluate the performance of a bespoke five-layer convolutional neural network when trained with two different image resolutions. The dataset is publicly available and contains 25,077 images categorised into 13,966 organic and 11,111 recyclable waste. Many researchers have used the same dataset to evaluate their proposed methods with varying accuracy results. However, these results are not directly comparable to our approach due to fundamental issues observed in their method and validation approach, including the lack of transparency in the experimental setup, which makes it impossible to replicate results. Another common issue associated with image classification is high computational cost which often results to high development time and prediction model size. Therefore, a lightweight model with high accuracy and a high level of methodology transparency is of particular importance in this domain. To investigate the computational cost issue, we used two image resolution sizes (i.e., 225×264 and 80×45) to explore the performance of our bespoke five-layer convolutional neural network in terms of development time, model size, predictive accuracy, and cross-entropy loss. Our intuition is that smaller image resolution will lead to a lightweight model with relatively high and/or comparable accuracy than the model trained with higher image resolution. In the absence of reliable baseline studies to compare our bespoke convolutional network in terms of accuracy and loss, we trained a random guess classifier to compare our results. The results show that small image resolution leads to a lighter model with less training time and the accuracy produced (80.88%) is better than the 76.19% yielded by the larger model. Both the small and large models performed better than the baseline which produced 50.05% accuracy. To encourage reproducibility of our results, all the experimental artifacts including preprocessed dataset and source code used in our experiments are made available in a public repository.

Pedigree investigation, clinical characteristics, and prognosis analysis of haematological disease patients with germline TET2 mutation

Background

Increasing germline gene mutations have been discovered in haematological malignancies with the development of next-generation sequencing (NGS), which is critical for proper clinical management and long-term follow-up of affected individuals. Tet methylcytosine dioxygenase 2 (TET2) is one of the most common mutations in haematological neoplasms. We aimed to compare the clinical characteristics of patients with germline and somatic TET2 mutations in haematological diseases and to analyse whether germline TET2 mutations have a family aggregation and tumour predisposition.

Methods

Out of 612 patients who underwent NGS of 34 recurrently mutated genes in haematological diseases, 100 haematological patients with TET2 mutations were selected for further study. Somatic mutations were detected by NGS in bone marrow/peripheral blood genomic DNA (gDNA). Germline TET2 mutations were validated in nail/hair gDNA by Sanger sequencing. Digital data were extracted from the haematology department of the West China Hospital of Sichuan University. TET2 mutation results were analysed by referencing online public databases (COSMIC and ClinVar).

Results

One hundred patients were studied, including 33 patients with germline and 67 patients with somatic TET2 mutations. For germline TET2 mutations, the variant allele frequency (VAF) was more stable (50.58% [40.5–55], P < 0.0001), and mutation sites recurrently occurred in three sites, unlike somatic TET2 mutations. Patients with germline TET2 mutations were younger (median age 48, 16–82 years) (P = 0.0058) and mainly suffered from myelodysplastic syndromes (MDS) (n = 13, 39.4%), while patients with somatic TET2 mutations were mainly affected by acute myeloid leukemia (AML) (n = 26, 38.8%) (P = 0.0004). Germline TET2 mutation affected the distribution of cell counts in the peripheral blood and bone marrow (P < 0.05); it was a poor prognostic factor for MDS patients via univariate analysis (HR = 5.3, 95% CI: 0.89–32.2, P = 0.0209) but not in multivariate analysis using the Cox regression model (P = 0.062).

Conclusions

Germline TET2 mutation might have a family aggregation, and TET2 may be a predisposition gene for haematological malignancy under the other gene mutations as the second hit. Germline TET2 mutation may play a role in the proportion of blood and bone marrow cells and, most importantly, may be an adverse factor for MDS patients.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12885-022-09347-0.

Simulation Research on the Methods of Multi-Gene Region Association Analysis Based on a Functional Linear Model

Genome-wide association analysis is an important approach to identify genetic variants associated with complex traits. Complex traits are not only affected by single gene loci, but also by the interaction of multiple gene loci. Studies of association between gene regions and quantitative traits are of great significance in revealing the genetic mechanism of biological development. There have been a lot of studies on single-gene region association analysis, but the application of functional linear models in multi-gene region association analysis is still less. In this paper, a functional multi-gene region association analysis test method is proposed based on the functional linear model. From the three directions of common multi-gene region method, multi-gene region weighted method and multi-gene region loci weighted method, that test method is studied combined with computer simulation. The following conclusions are obtained through computer simulation: (a) The functional multi-gene region association analysis test method has higher power than the functional single gene region association analysis test method; (b) The functional multi-gene region weighted method performs better than the common functional multi-gene region method; (c) the functional multi-gene region loci weighted method is the best method for association analysis on three directions of the common multi-gene region method; (d) the performance of the Step method and Multi-gene region loci weighted Step for multi-gene regions is the best in general. Functional multi-gene region association analysis test method can theoretically provide a feasible method for the study of complex traits affected by multiple genes.

Performing Genome-Wide Association Studies Using rMVP

Genome wide association study (GWAS), which is a powerful tool to detect the relationship between the traits of interest and high-density markers, has provided unprecedented insights into the genetic basis of quantitative variation for complex traits. Along with the development of high-throughput sequencing technology, both sample sizes and marker sizes are increasing rapidly, which make computations more challenging than ever. Therefore, to efficiently process big data with limited computing resources in a reasonable time and to use state-of-the-art statistical models to reduce false positive and false negative rates have always been hot topics in the domain of GWAS. In this chapter, we describe how to perform GWAS using an R package, rMVP, which includes data preparation, evaluation of population structure, association tests by different models, and high-quality visualization of GWAS results.

Unraveling the Genomic Architecture of the CYP3A Locus and ADME Genes for Personalized Tacrolimus Dosing

BACKGROUND

Tacrolimus (TAC) is an immunosuppressant widely prescribed following an allogenic organ transplant. Due to wide interindividual pharmacokinetic (PK) variability, optimizing TAC dosing based on genetic factors is required to minimize nephrotoxicity and acute rejections.

METHODS

We enrolled 1133 participants receiving TAC from 4 cohorts, consisting of 3 with kidney transplant recipients and 1 with healthy males from clinical trials. The effects of clinical factors were estimated to appropriately control confounding variables. A genome-wide association study, haplotype analysis, and a gene-based association test were conducted using the Korea Biobank Array or targeted sequencing for 114 pharmacogenes.

RESULTS

Genome-wide association study verified that CYP3A5*3 is the only common variant associated with TAC PK variability in Koreans. We detected several CYP3A5 and CYP3A4 rare variants that could potentially affect TAC metabolism. The haplotype structure of CYP3A5 stratified by CYP3A5*3 was a significant factor for CYP3A5 rare variant interpretation. CYP3A4 rare variant carriers among CYP3A5 intermediate metabolizers displayed higher TAC trough levels. Gene-based association tests in the 61 absorption, distribution, metabolism, and excretion genes revealed that CYP1A1 are associated with additional TAC PK variability: CYP1A1 rare variant carriers among CYP3A5 poor metabolizers showed lower TAC trough levels than the noncarrier controls.

CONCLUSIONS

Our study demonstrates that rare variant profiling of CYP3A5 and CYP3A4, combined with the haplotype structures of CYP3A locus, provide additive value for personalized TAC dosing. We also identified a novel association between CYP1A1 rare variants and TAC PK variability in the CYP3A5 nonexpressers that needs to be further investigated.

A two-stage testing strategy for detecting genes×environment interactions in association studies

Identifying gene×environment (G×E) interactions, especially when rare variants are included in genome-wide association studies, is a major challenge in statistical genetics. However, the detection of G×E interactions is very important for understanding the etiology of complex diseases. Although currently some statistical methods have been developed to detect the interactions between genes and environment, the detection of the interactions for the case of rare variants is still limited. Therefore, it is particularly important to develop a new method to detect the interactions between genes and environment for rare variants. In this study, we extend an existing method of adaptive combination of P-values (ADA) and design a novel strategy (called iSADA) for testing the effects of G×E interactions for rare variants. We propose a new two-stage test to detect the interactions between genes and environment in a certain region of a chromosome or even for the whole genome. First, the score statistic is used to test the associations between trait value and the interaction terms of genes and environment and obtain the original P-values. Then, based on the idea of the ADA method, we further construct a full test statistic via the P-values of the preliminary tests in the first stage, so that we can comprehensively test the interactions between genes and environment in the considered genome region. Simulation studies are conducted to compare our proposed method with other existing methods. The results show that the iSADA has higher power than other methods in each case. A GAW17 data set is also applied to illustrate the applicability of the new method.

Controlling for human population stratification in rare variant association studies

Population stratification is a confounder of genetic association studies. In analyses of rare variants, corrections based on principal components (PCs) and linear mixed models (LMMs) yield conflicting conclusions. Studies evaluating these approaches generally focused on limited types of structure and large sample sizes. We investigated the properties of several correction methods through a large simulation study using real exome data, and several within- and between-continent stratification scenarios. We considered different sample sizes, with situations including as few as 50 cases, to account for the analysis of rare disorders. Large samples showed that accounting for stratification was more difficult with a continental than with a worldwide structure. When considering a sample of 50 cases, an inflation of type-I-errors was observed with PCs for small numbers of controls (≤ 100), and with LMMs for large numbers of controls (≥ 1000). We also tested a novel local permutation method (LocPerm), which maintained a correct type-I-error in all situations. Powers were equivalent for all approaches pointing out that the key issue is to properly control type-I-errors. Finally, we found that power of analyses including small numbers of cases can be increased, by adding a large panel of external controls, provided an appropriate stratification correction was used.

Integrative analysis of multi-omics data for discovering low-frequency variants associated with low-density lipoprotein cholesterol levels

MOTIVATION

The abundance of omics data has facilitated integrative analyses of single and multiple molecular layers with genome-wide association studies focusing on common variants. Built on its successes, we propose a general analysis framework to leverage multi-omics data with sequencing data to improve the statistical power of discovering new associations and understanding of the disease susceptibility due to low-frequency variants. The proposed test features its robustness to model misspecification, high power across a wide range of scenarios and the potential of offering insights into the underlying genetic architecture and disease mechanisms.

RESULTS

Using the Framingham Heart Study data, we show that low-frequency variants are predictive of DNA methylation, even after conditioning on the nearby common variants. In addition, DNA methylation and gene expression provide complementary information to functional genomics. In the Avon Longitudinal Study of Parents and Children with a sample size of 1497, one gene CLPTM1 is identified to be associated with low-density lipoprotein cholesterol levels by the proposed powerful adaptive gene-based test integrating information from gene expression, methylation and enhancer-promoter interactions. It is further replicated in the TwinsUK study with 1706 samples. The signal is driven by both low-frequency and common variants.

AVAILABILITY AND IMPLEMENTATION

Models are available at https://github.com/ytzhong/DNAm.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Genome-wide association study and its applications in the non-model crop Sesamum indicum

BACKGROUND

Sesame is a rare example of non-model and minor crop for which numerous genetic loci and candidate genes underlying features of interest have been disclosed at relatively high resolution. These progresses have been achieved thanks to the applications of the genome-wide association study (GWAS) approach. GWAS has benefited from the availability of high-quality genomes, re-sequencing data from thousands of genotypes, extensive transcriptome sequencing, development of haplotype map and web-based functional databases in sesame.

RESULTS

In this paper, we reviewed the GWAS methods, the underlying statistical models and the applications for genetic discovery of important traits in sesame. A novel online database SiGeDiD ( http://sigedid.ucad.sn/ ) has been developed to provide access to all genetic and genomic discoveries through GWAS in sesame. We also tested for the first time, applications of various new GWAS multi-locus models in sesame.

CONCLUSIONS

Collectively, this work portrays steps and provides guidelines for efficient GWAS implementation in sesame, a non-model crop.

Identification of Novel Genomic-Variant Patterns of OR56A5, OR52L1, and CTSD in Retinitis Pigmentosa Patients by Whole-Exome Sequencing

Inherited retinal dystrophies (IRDs) are rare but highly heterogeneous genetic disorders that affect individuals and families worldwide. However, given its wide variability, its analysis of the driver genes for over 50% of the cases remains unexplored. The present study aims to identify novel driver genes, disease-causing variants, and retinitis pigmentosa (RP)-associated pathways. Using family-based whole-exome sequencing (WES) to identify putative RP-causing rare variants, we identified a total of five potentially pathogenic variants located in genes OR56A5, OR52L1, CTSD, PRF1, KBTBD13, and ATP2B4. Of the variants present in all affected individuals, genes OR56A5, OR52L1, CTSD, KBTBD13, and ATP2B4 present as missense mutations, while PRF1 and CTSD present as frameshift variants. Sanger sequencing confirmed the presence of the novel pathogenic variant PRF1 (c.124_128del) that has not been reported previously. More causal-effect or evidence-based studies will be required to elucidate the precise roles of these SNPs in the RP pathogenesis. Taken together, our findings may allow us to explore the risk variants based on the sequencing data and upgrade the existing variant annotation database in Taiwan. It may help detect specific eye diseases such as retinitis pigmentosa in East Asia.

Interplay between IL6 and CRIM1 in thiopurine intolerance due to hematological toxicity in leukemic patients with wild-type NUDT15 and TPMT

NUDT15 and TPMT variants are strong genetic determinants of thiopurine-induced hematological toxicity. Despite the impact of homozygous CRIM1 on thiopurine toxicity, several patients with wild-type NUDT15, TPMT, and CRIM1 experience thiopurine toxicity, therapeutic failure, and relapse of acute lymphoblastic leukemia (ALL). Novel pharmacogenetic interactions associated with thiopurine intolerance from hematological toxicities were investigated using whole-exome sequencing for last-cycle 6-mercaptopurine dose intensity percentages (DIP) tolerated by pediatric ALL patients (N = 320). IL6 rs13306435 carriers (N = 19) exhibited significantly lower DIP (48.0 ± 27.3%) than non-carriers (N = 209, 69.9 ± 29.0%; p = 0.0016 and 0.0028 by t test and multiple linear regression, respectively). Among 19 carriers, 7 with both heterozygous IL6 rs13306435 and CRIM1 rs3821169 showed significantly decreased DIP (24.7 ± 8.9%) than those with IL6 (N = 12, 61.6 ± 25.1%) or CRIM1 (N = 94, 68.1 ± 28.4%) variants. IL6 and CRIM1 variants showed marked inter-ethnic variability. Four-gene-interplay models revealed the best odds ratio (8.06) and potential population impact [relative risk (5.73), population attributable fraction (58%), number needed to treat (3.67), and number needed to genotype (12.50)]. Interplay between IL6 rs13306435 and CRIM1 rs3821169 was suggested as an independent and/or additive genetic determinant of thiopurine intolerance beyond NUDT15 and TPMT in pediatric ALL.

Detecting Clustered Independent Rare Variant Associations Using Genetic Algorithms

The availability of an increasing collection of sequencing data provides the opportunity to study genetic variation with an unprecedented level of detail. There is much interest in uncovering the role of rare variants and their contribution to disease. However, detecting associations of rare variants with small minor allele frequencies (MAF) and modest effects remains a challenge for rare variant association methods. Due to this low signal-to-noise ratio, most methods are underpowered to detect associations even when conducting rare variant association tests at the gene level. We present a new method for detecting rare variant associations. The algorithm consists of two steps. In the first step, a genetic algorithm searches for a promising genomic region containing a collection of genes with causal rare variants. In the second step, a genetic algorithm aims at removing false positives from the located genomic region. We tested the proposed method with a collection of datasets obtained from real exome data. The proposed method possesses sufficient power for detecting associations of rare variants with complex phenotypes. This method can be used for studying the contribution of rare variants with complex disease, particularly in cases where single-variant or gene-based tests are underpowered.

Whole-genome sequencing reveals KRTAP1-1 as a novel genetic variant associated with antidepressant treatment outcomes

Achieving remission following initial antidepressant therapy in patients with major depressive disorder (MDD) is an important clinical result. Making predictions based on genetic markers holds promise for improving the remission rate. However, genetic variants found in previous genetic studies do not provide robust evidence to aid pharmacogenetic decision-making in clinical settings. Thus, the objective of this study was to perform whole-genome sequencing (WGS) using genomic DNA to identify genetic variants associated with the treatment outcomes of selective serotonin reuptake inhibitors (SSRIs). We performed WGS on 100 patients with MDD who were treated with escitalopram (discovery set: 36 remitted and 64 non-remitted). The findings were applied to an additional 553 patients with MDD who were treated with SSRIs (replication set: 185 remitted and 368 non-remitted). A novel loss-of-function variant (rs3213755) in keratin-associated protein 1-1 (KRTAP1-1) was identified in this study. This rs3213755 variant was significantly associated with remission following antidepressant treatment (p = 0.0184, OR 3.09, 95% confidence interval [CI] 1.22-7.80 in the discovery set; p = 0.00269, OR 1.75, 95% CI 1.22-2.53 in the replication set). Moreover, the expression level of KRTAP1-1 in surgically resected human temporal lobe samples was significantly associated with the rs3213755 genotype. WGS studies on a larger sample size in various ethnic groups are needed to investigate genetic markers useful in the pharmacogenetic prediction of remission following antidepressant treatment.

Identification of Influential Variants in Significant Aggregate Rare Variant Tests

INTRODUCTION

Studies that examine the role of rare variants in both simple and complex disease are increasingly common. Though the usual approach of testing rare variants in aggregate sets is more powerful than testing individual variants, it is of interest to identify the variants that are plausible drivers of the association. We present a novel method for prioritization of rare variants after a significant aggregate test by quantifying the influence of the variant on the aggregate test of association.

METHODS

In addition to providing a measure used to rank variants, we use outlier detection methods to present the computationally efficient Rare Variant Influential Filtering Tool (RIFT) to identify a subset of variants that influence the disease association. We evaluated several outlier detection methods that vary based on the underlying variance measure: interquartile range (Tukey fences), median absolute deviation, and SD. We performed 1,000 simulations for 50 regions of size 3 kb and compared the true and false positive rates. We compared RIFT using the Inner Tukey to 2 existing methods: adaptive combination of p values (ADA) and a Bayesian hierarchical model (BeviMed). Finally, we applied this method to data from our targeted resequencing study in idiopathic pulmonary fibrosis (IPF).

RESULTS

All outlier detection methods observed higher sensitivity to detect uncommon variants (0.001 < minor allele frequency, MAF > 0.03) compared to very rare variants (MAF <0.001). For uncommon variants, RIFT had a lower median false positive rate compared to the ADA. ADA and RIFT had significantly higher true positive rates than that observed for BeviMed. When applied to 2 regions found previously associated with IPF including 100 rare variants, we identified 6 polymorphisms with the greatest evidence for influencing the association with IPF.

DISCUSSION

In summary, RIFT has a high true positive rate while maintaining a low false positive rate for identifying polymorphisms influencing rare variant association tests. This work provides an approach to obtain greater resolution of the rare variant signals within significant aggregate sets; this information can provide an objective measure to prioritize variants for follow-up experimental studies and insight into the biological pathways involved.

A fast and powerful aggregated Cauchy association test for joint analysis of multiple phenotypes

BACKGROUND

Pleiotropy is a widespread phenomenon in complex human diseases. Jointly analyzing multiple phenotypes can improve power performance of detecting genetic variants and uncover the underlying genetic mechanism.

OBJECTIVE

This study aims to detect the association between genetic variants in a genomic region and multiple phenotypes.

METHODS

We develop the aggregated Cauchy association test to detect the association between rare variants in a genomic region and multiple phenotypes (abbreviated as "Multi-ACAT"). Multi-ACAT first detects the association between each rare variant and multiple phenotypes based on reverse regression and obtains variant-level p-values, then takes linear combination of transformed p-values as the test statistic which approximately follows Cauchy distribution under the null hypothesis.

RESULTS

Extensive simulation studies show that when the proportion of causal variants in a genomic region is extremely small, Multi-ACAT is more powerful than the other several methods and is robust to bi-directional effects of causal variants. Finally, we illustrate our proposed method by analyzing two phenotypes [systolic blood pressure (SBP) and diastolic blood pressure (DBP)] from Genetic Analysis Workshop 19 (GAW19).

CONCLUSION

The Multi-ACAT computes extremely fast, does not consider complex distributions of multiple correlated phenotypes, and can be applied to the case with noise phenotypes.

Identifying rare variants for quantitative traits in extreme samples of population via Kullback-Leibler distance

Background

The rapid development of sequencing technology and simultaneously the availability of large quantities of sequence data has facilitated the identification of rare variant associated with quantitative traits. However, existing statistical methods depend on certain assumptions and thus lacking uniform power. The present study focuses on mapping rare variant associated with quantitative traits.

Results

In the present study, we proposed a two-stage strategy to identify rare variant of quantitative traits using phenotype extreme selection design and Kullback-Leibler distance, where the first stage was association analysis and the second stage was fine mapping. We presented a statistic and a linkage disequilibrium measure for the first stage and the second stage, respectively. Theory analysis and simulation study showed that (1) the power of the proposed statistic for association analysis increased with the stringency of the sample selection and was affected slightly by non-causal variants and opposite effect variants, (2) the statistic here achieved higher power than three commonly used methods, and (3) the linkage disequilibrium measure for fine mapping was independent of the frequencies of non-causal variants and simply dependent on the frequencies of causal variants.

Conclusions

We conclude that the two-stage strategy here can be used effectively to mapping rare variant associated with quantitative traits.

Simultaneous Detection of Signal Regions Using Quadratic Scan Statistics With Applications to Whole Genome Association Studies

We consider in this paper detection of signal regions associated with disease outcomes in whole genome association studies. Gene- or region-based methods have become increasingly popular in whole genome association analysis as a complementary approach to traditional individual variant analysis. However, these methods test for the association between an outcome and the genetic variants in a pre-specified region, e.g., a gene. In view of massive intergenic regions in whole genome sequencing (WGS) studies, we propose a computationally efficient quadratic scan (Q-SCAN) statistic based method to detect the existence and the locations of signal regions by scanning the genome continuously. The proposed method accounts for the correlation (linkage disequilibrium) among genetic variants, and allows for signal regions to have both causal and neutral variants, and the effects of signal variants to be in different directions. We study the asymptotic properties of the proposed Q-SCAN statistics. We derive an empirical threshold that controls for the family-wise error rate, and show that under regularity conditions the proposed method consistently selects the true signal regions. We perform simulation studies to evaluate the finite sample performance of the proposed method. Our simulation results show that the proposed procedure outperforms the existing methods, especially when signal regions have causal variants whose effects are in different directions, or are contaminated with neutral variants. We illustrate Q-SCAN by analyzing the WGS data from the Atherosclerosis Risk in Communities study.

Dynamic incorporation of multiple in silico functional annotations empowers rare variant association analysis of large whole-genome sequencing studies at scale

Large-scale whole-genome sequencing studies have enabled the analysis of rare variants (RVs) associated with complex phenotypes. Commonly used RV association tests have limited scope to leverage variant functions. We propose STAAR (variant-set test for association using annotation information), a scalable and powerful RV association test method that effectively incorporates both variant categories and multiple complementary annotations using a dynamic weighting scheme. For the latter, we introduce 'annotation principal components', multidimensional summaries of in silico variant annotations. STAAR accounts for population structure and relatedness and is scalable for analyzing very large cohort and biobank whole-genome sequencing studies of continuous and dichotomous traits. We applied STAAR to identify RVs associated with four lipid traits in 12,316 discovery and 17,822 replication samples from the Trans-Omics for Precision Medicine Program. We discovered and replicated new RV associations, including disruptive missense RVs of NPC1L1 and an intergenic region near APOC1P1 associated with low-density lipoprotein cholesterol.

Detecting X-linked common and rare variant effects in family-based sequencing studies

The breakthroughs in next generation sequencing have allowed us to access data consisting of both common and rare variants, and in particular to investigate the impact of rare genetic variation on complex diseases. Although rare genetic variants are thought to be important components in explaining genetic mechanisms of many diseases, discovering these variants remains challenging, and most studies are restricted to population-based designs. Further, despite the shift in the field of genome-wide association studies (GWAS) towards studying rare variants due to the "missing heritability" phenomenon, little is known about rare X-linked variants associated with complex diseases. For instance, there is evidence that X-linked genes are highly involved in brain development and cognition when compared with autosomal genes; however, like most GWAS for other complex traits, previous GWAS for mental diseases have provided poor resources to deal with identification of rare variant associations on X-chromosome. In this paper, we address the two issues described above by proposing a method that can be used to test X-linked variants using sequencing data on families. Our method is much more general than existing methods, as it can be applied to detect both common and rare variants, and is applicable to autosomes as well. Our simulation study shows that the method is efficient, and exhibits good operational characteristics. An application to the University of Miami Study on Genetics of Autism and Related Disorders also yielded encouraging results.

Gene-Based Association Testing of Dichotomous Traits With Generalized Functional Linear Mixed Models Using Extended Pedigrees: Applications to Age-Related Macular Degeneration

Genetics plays a role in age-related macular degeneration (AMD), a common cause of blindness in the elderly. There is a need for powerful methods for carrying out region-based association tests between a dichotomous trait like AMD and genetic variants on family data. Here, we apply our new generalized functional linear mixed models (GFLMM) developed to test for gene-based association in a set of AMD families. Using common and rare variants, we observe significant association with two known AMD genes: CFH and ARMS2. Using rare variants, we find suggestive signals in four genes: ASAH1, CLEC6A, TMEM63C, and SGSM1. Intriguingly, ASAH1 is down-regulated in AMD aqueous humor, and ASAH1 deficiency leads to retinal inflammation and increased vulnerability to oxidative stress. These findings were made possible by our GFLMM which model the effect of a major gene as a fixed mean, the polygenic contributions as a random variation, and the correlation of pedigree members by kinship coefficients. Simulations indicate that the GFLMM likelihood ratio tests (LRTs) accurately control the Type I error rates. The LRTs have similar or higher power than existing retrospective kernel and burden statistics. Our GFLMM-based statistics provide a new tool for conducting family-based genetic studies of complex diseases. Supplementary materials for this article, including a standardized description of the materials available for reproducing the work, are available as an online supplement.

Mutations in sphingolipid metabolism genes are associated with ADHD

Attention deficit hyperactivity disorder (ADHD) is the most prevalent neurodevelopmental disorder in children, with genetic factors accounting for 75-80% of the phenotypic variance. Recent studies have suggested that ADHD patients might present with atypical central myelination that can persist into adulthood. Given the essential role of sphingolipids in myelin formation and maintenance, we explored genetic variation in sphingolipid metabolism genes for association with ADHD risk. Whole-exome genotyping was performed in three independent cohorts from disparate regions of the world, for a total of 1520 genotyped subjects. Cohort 1 (MTA (Multimodal Treatment study of children with ADHD) sample, 371 subjects) was analyzed as the discovery cohort, while cohorts 2 (Paisa sample, 298 subjects) and 3 (US sample, 851 subjects) were used for replication. A set of 58 genes was manually curated based on their roles in sphingolipid metabolism. A targeted exploration for association between ADHD and 137 markers encoding for common and rare potentially functional allelic variants in this set of genes was performed in the screening cohort. Single- and multi-locus additive, dominant and recessive linear mixed-effect models were used. During discovery, we found statistically significant associations between ADHD and variants in eight genes (GALC, CERS6, SMPD1, SMPDL3B, CERS2, FADS3, ELOVL5, and CERK). Successful local replication for associations with variants in GALC, SMPD1, and CERS6 was demonstrated in both replication cohorts. Variants rs35785620, rs143078230, rs398607, and rs1805078, associated with ADHD in the discovery or replication cohorts, correspond to missense mutations with predicted deleterious effects. Expression quantitative trait loci analysis revealed an association between rs398607 and increased GALC expression in the cerebellum.

Views on GWAS statistical analysis

Genome-wide association study (GWAS) is a popular approach to investigate relationships between genetic information and diseases. A number of associations are tested in a study and the results are often corrected using multiple adjustment methods. It is observed that GWAS studies suffer adequate statistical power for reliability. Hence, we document known models for reliability assessment using improved statistical power in GWAS analysis.

Novel exomic rare variants associated with venous thrombosis

Exomic rare variant polymorphisms (c. 300 000) were analysed in the Scripps Venous Thrombosis (VTE) registry (subjects aged <55 years). Besides coagulation factor V (F5) single nucleotide polymorphisms (SNPs), family with sequence similarity 134, member B (FAM134B; rs78314670, Arg127Cys) and myosin heavy chain 8 (MYH8; rs111567318, Glu1838Ala) SNPs were associated with recurrent VTE (n = 34 cases) (false discovery rate-adjusted P < 0·05). FAM134B (rs78314670) was associated with low plasma levels of anticoagulant glucosylceramide. Analysis of 50 chr17p13.1 MYH rare SNPs (clustered skeletal myosin heavy chain genes) using collapsing methods was associated with recurrent VTE (P = 2·70 ×10^-16 ). When intravenously injected, skeletal muscle myosin was pro-coagulant in a haemophilia mouse tail bleeding model. Thus, FAM134B and MYH genetic variants are plausibly linked to VTE risk.

Determining population stratification and subgroup effects in association studies of rare genetic variants for nicotine dependence

Supplemental Digital Content is available in the text.

Background

Rare variants (minor allele frequency < 1% or 5 %) can help researchers to deal with the confounding issue of ‘missing heritability’ and have a proven role in dissecting the etiology for human diseases and complex traits.

Methods

We extended the combined multivariate and collapsing (CMC) and weighted sum statistic (WSS) methods and accounted for the effects of population stratification and subgroup effects using stratified analyses by the principal component analysis, named here as ‘str-CMC’ and ‘str-WSS’. To evaluate the validity of the extended methods, we analyzed the Genetic Architecture of Smoking and Smoking Cessation database, which includes African Americans and European Americans genotyped on Illumina Human Omni2.5, and we compared the results with those obtained with the sequence kernel association test (SKAT) and its modification, SKAT-O that included population stratification and subgroup effect as covariates. We utilized the Cochran–Mantel–Haenszel test to check for possible differences in single nucleotide polymorphism allele frequency between subgroups within a gene. We aimed to detect rare variants and considered population stratification and subgroup effects in the genomic region containing 39 acetylcholine receptor-related genes.

Results

The Cochran–Mantel–Haenszel test as applied to GABRG2 (P = 0.001) was significant. However, GABRG2 was detected both by str-CMC (P= 8.04E-06) and str-WSS (P= 0.046) in African Americans but not by SKAT or SKAT-O.

Conclusions

Our results imply that if associated rare variants are only specific to a subgroup, a stratified analysis might be a better approach than a combined analysis.

Whole-exome sequencing of 81 individuals from 27 multiply affected bipolar disorder families

Bipolar disorder (BD) is a highly heritable neuropsychiatric disease characterized by recurrent episodes of depression and mania. Research suggests that the cumulative impact of common alleles explains 25-38% of phenotypic variance, and that rare variants may contribute to BD susceptibility. To identify rare, high-penetrance susceptibility variants for BD, whole-exome sequencing (WES) was performed in three affected individuals from each of 27 multiply affected families from Spain and Germany. WES identified 378 rare, non-synonymous, and potentially functional variants. These spanned 368 genes, and were carried by all three affected members in at least one family. Eight of the 368 genes harbored rare variants that were implicated in at least two independent families. In an extended segregation analysis involving additional family members, five of these eight genes harbored variants showing full or nearly full cosegregation with BD. These included the brain-expressed genes RGS12 and NCKAP5, which were considered the most promising BD candidates on the basis of independent evidence. Gene enrichment analysis for all 368 genes revealed significant enrichment for four pathways, including genes reported in de novo studies of autism (p_adj < 0.006) and schizophrenia (p_adj = 0.015). These results suggest a possible genetic overlap with BD for autism and schizophrenia at the rare-sequence-variant level. The present study implicates novel candidate genes for BD development, and may contribute to an improved understanding of the biological basis of this common and often devastating disease.

Regulation of glucose and lipid metabolism in health and disease

Glucose and fatty acids are the major sources of energy for human body. Cholesterol, the most abundant sterol in mammals, is a key component of cell membranes although it does not generate ATP. The metabolisms of glucose, fatty acids and cholesterol are often intertwined and regulated. For example, glucose can be converted to fatty acids and cholesterol through de novo lipid biosynthesis pathways. Excessive lipids are secreted in lipoproteins or stored in lipid droplets. The metabolites of glucose and lipids are dynamically transported intercellularly and intracellularly, and then converted to other molecules in specific compartments. The disorders of glucose and lipid metabolism result in severe diseases including cardiovascular disease, diabetes and fatty liver. This review summarizes the major metabolic aspects of glucose and lipid, and their regulations in the context of physiology and diseases.

Familial Alzheimer's Disease and Recessive Modifiers

Alzheimer’s disease (AD) is progressive brain disorder that affects ~ 50 million people worldwide and has no current effective treatment. AD age of onset (ADAOO) has shown to be critical for the identification of genes that modify the appearance of AD signs and symptoms in a specific population. We clinically characterized and whole-exome genotyped 71 individuals with AD from the Paisa genetic isolate, segregating the (PSEN1) E280A dominant fully penetrant mutation, and analyzed the potential recessive effects of ~ 50,000 common functional genomic variants to the ADAOO. Standard quality control and filtering procedures were applied, and recessive single- and multi-locus linear mixed-effects models were used. We identified genetic variants in the SLC9C1, CSN1S1, and LOXL4 acting recessively to delay ADAOO up to ~ 11, ~ 6, and ~ 9 years on average, respectively. In contrast, the CC recessive genotype in marker DHRS4L2-rs2273946 accelerates ADAOO by ~ 8 years. This study, reports new recessive variants modifying ADAOO in PSEN1 E280A mutation carriers. This set of genes are implicated in important biological processes and molecular functions commonly affected by genes associated with the etiology of AD such as APP, APOE, and CLU. Future functional studies using modern techniques such as induced pluripotent stem cells will allow a better understanding of the over expression and down regulation of these recessive modifier variants and hence the pathogenesis of AD. These results are important for prediction of AD and ultimately, substantial to develop new therapeutic strategies for individuals at risk or affected by AD.

A general statistic to test an optimally weighted combination of common and/or rare variants

Both genome-wide association study and next-generation sequencing data analyses are widely employed to identify disease susceptible common and/or rare genetic variants. Rare variants generally have large effects though they are hard to detect due to their low frequencies. Currently, many existing statistical methods for rare variants association studies employ a weighted combination scheme, which usually puts subjective weights or suboptimal weights based on some adhoc assumptions (e.g., ignoring dependence between rare variants). In this study, we analytically derived optimal weights for both common and rare variants and proposed a general and novel approach to test association between an optimally weighted combination of variants (G-TOW) in a gene or pathway for a continuous or dichotomous trait while easily adjusting for covariates. Results of the simulation studies show that G-TOW has properly controlled type I error rates and it is the most powerful test among the methods we compared when testing effects of either both rare and common variants or rare variants only. We also illustrate the effectiveness of G-TOW using the Genetic Analysis Workshop 17 (GAW17) data. Additionally, we applied G-TOW and other competitive methods to test disease-associated genes in real data of schizophrenia. The G-TOW has successfully verified genes FYN and VPS39 which are associated with schizophrenia reported in existing publications. Both of these genes are missed by the weighted sum statistic and the sequence kernel association test. Simulation study and real data analysis indicate that G-TOW is a powerful test.

Genetic Variation Underpinning ADHD Risk in a Caribbean Community

Attention Deficit Hyperactivity Disorder (ADHD) is a highly heritable and prevalent neurodevelopmental disorder that frequently persists into adulthood. Strong evidence from genetic studies indicates that single nucleotide polymorphisms (SNPs) harboured in the ADGRL3 (LPHN3), SNAP25, FGF1, DRD4, and SLC6A2 genes are associated with ADHD. We genotyped 26 SNPs harboured in genes previously reported to be associated with ADHD and evaluated their potential association in 386 individuals belonging to 113 nuclear families from a Caribbean community in Barranquilla, Colombia, using family-based association tests. SNPs rs362990-SNAP25 (T allele; p = 2.46 × 10⁻⁴), rs2282794-FGF1 (A allele; p = 1.33 × 10⁻²), rs2122642-ADGRL3 (C allele, p = 3.5 × 10⁻²), and ADGRL3 haplotype CCC (markers rs1565902-rs10001410-rs2122642, OR = 1.74, P_permuted = 0.021) were significantly associated with ADHD. Our results confirm the susceptibility to ADHD conferred by SNAP25, FGF1, and ADGRL3 variants in a community with a significant African American component, and provide evidence supporting the existence of specific patterns of genetic stratification underpinning the susceptibility to ADHD. Knowledge of population genetics is crucial to define risk and predict susceptibility to disease.

Whole genome sequencing and rare variant analysis in essential tremor families

Essential tremor (ET) is one of the most common movement disorders. The etiology of ET remains largely unexplained. Whole genome sequencing (WGS) is likely to be of value in understanding a large proportion of ET with Mendelian and complex disease inheritance patterns. In ET families with Mendelian inheritance patterns, WGS may lead to gene identification where WES analysis failed to identify the causative single nucleotide variant (SNV) or indel due to incomplete coverage of the entire coding region of the genome, in addition to accurate detection of larger structural variants (SVs) and copy number variants (CNVs). Alternatively, in ET families with complex disease inheritance patterns with gene x gene and gene x environment interactions enrichment of functional rare coding and non-coding variants may explain the heritability of ET. We performed WGS in eight ET families (n = 40 individuals) enrolled in the Family Study of Essential Tremor. The analysis included filtering WGS data based on allele frequency in population databases, rare SNV and indel classification and association testing using the Mixed-Model Kernel Based Adaptive Cluster (MM-KBAC) test. A separate analysis of rare SV and CNVs segregating within ET families was also performed. Prioritization of candidate genes identified within families was performed using phenolyzer. WGS analysis identified candidate genes for ET in 5/8 (62.5%) of the families analyzed. WES analysis in a subset of these families in our previously published study failed to identify candidate genes. In one family, we identified a deleterious and damaging variant (c.1367G>A, p.(Arg456Gln)) in the candidate gene, CACNA1G, which encodes the pore forming subunit of T-type Ca(2+) channels, CaV3.1, and is expressed in various motor pathways and has been previously implicated in neuronal autorhythmicity and ET. Other candidate genes identified include SLIT3 which encodes an axon guidance molecule and in three families, phenolyzer prioritized genes that are associated with hereditary neuropathies (family A, KARS, family B, KIF5A and family F, NTRK1). Functional studies of CACNA1G and SLIT3 suggest a role for these genes in ET disease pathogenesis.

Star Allele-Based Haplotyping versus Gene-Wise Variant Burden Scoring for Predicting 6-Mercaptopurine Intolerance in Pediatric Acute Lymphoblastic Leukemia Patients

Nudix Hydrolase 15 (NUDT15) and Thiopurine S-Methyltransferase (TPMT) are strong genetic determinants of thiopurine toxicity in pediatric acute lymphoblastic leukemia (ALL) patients. Since patients with NUDT15 or TPMT deficiency suffer severe adverse drug reactions, star (*) allele-based haplotypes have been used to predict an optimal 6-mercaptopurine (6-MP) dosing. However, star allele haplotyping suffers from insufficient, inconsistent, and even conflicting designations with uncertain and/or unknown functional alleles. Gene-wise variant burden (GVB) scoring enables us to utilize next-generation sequencing (NGS) data to predict 6-MP intolerance in children with ALL. Whole exome sequencing was performed for 244 pediatric ALL patients under 6-MP treatments. We assigned star alleles with PharmGKB haplotype set translational table. GVB for NUDT15 and TPMT was computed by aggregating in silico deleteriousness scores of multiple coding variants for each gene. Poor last-cycle dose intensity percent (DIP < 25%) was considered as 6-MP intolerance, resulting therapeutic failure of ALL. DIPs showed significant differences ( p < 0.05) among NUDT15 poor (PM, n = 1), intermediate (IM, n = 48), and normal (NM, n = 195) metabolizers. TPMT exhibited no PM and only seven IMs. GVB showed significant differences among the different haplotype groups of both NUDT15 and TPMT ( p < 0.05). Kruskal–Wallis test for DIP values showed statistical significances for the seven different GVB score bins of NUDT15. GVB^NUDT15 outperformed the star allele-based haplotypes in predicting patients with reduced last-cycle DIPs at all DIP threshold levels (i.e., 5%, 10%, 15%, and 25%). In NUDT15-and-TPMT combined interaction analyses, GVB^NUDT15,TPMT outperformed star alleles [area under the receiver operating curve (AUROC) = 0.677 vs. 0.645] in specificity (0.813 vs. 0.796), sensitivity (0.526 vs. 0.474), and positive (0.192 vs. 0.164) and negative (0.953 vs. 0.947) predictive values. Overall, GVB correctly classified five more patients (i.e., one into below and four into above 25% DIP groups) than did star allele haplotypes. GVB analysis demonstrated that 6-MP intolerance in pediatric ALL can be reliably predicted by aggregating NGS-based common, rare, and novel variants together without hampering the predictive power of the conventional haplotype analysis.

Telomerase Variants in Patients with Cirrhosis Awaiting Liver Transplantation

Telomeres are repetitive DNA sequences that protect the ends of linear chromosomes, and they are maintained by a ribonucleoprotein complex called telomerase. Variants in genes encoding for telomerase components have been associated with a spectrum of disease in the lung, skin, bone marrow, and liver. Mutations in the telomerase reverse transcriptase and telomerase RNA component genes have been observed at a higher prevalence in patients with liver disease compared with the general population; however, the presence of variants in other components of the telomerase complex and their impact on clinical outcomes has not been explored. We evaluated 86 patients with end-stage liver disease for variants in an expanded panel of eight genes, and found that 17 patients (20%) had likely deleterious variants by in silico analysis. Seven unique likely deleterious variants were identified in the regulator of telomere elongation helicase 1 (RTEL1) gene that encodes for a DNA helicase important in telomere maintenance and genomic stability. In gene burden association analysis of their clinical data, the presence of any RTEL1 variant was associated with a 29% lower baseline white blood cell count (95% confidence interval [CI], -7% to -46%; P Value = 0.01) compared with patients without RTEL1 variants, and the presence of any exonic missense RTEL1 variant was associated with a 42% lower baseline platelet count (95% CI, -5% to -65%: P Value = 0.03). The presence of any telomerase variant was associated with an increased number of readmissions within 1 year after transplantation demonstrated by an incident rate ratio (IRR) of 3.15 (95% CI, 1.22 to 8.57). No association with survival was observed. Conclusion: Among patients who underwent liver transplantation, the presence of any exonic missense variant was associated with a longer postoperative length of stay with an IRR of 2.16 (95% CI, 1.31 to 3.68).

GIGI-Quick: a fast approach to impute missing genotypes in genome-wide association family data

Summary

Genome-wide association studies have become common over the last ten years, with a shift towards targeting rare variants, especially in pedigree-data. Despite lower costs, sequencing for rare variants still remains expensive. To have a relatively large sample with acceptable cost, imputation approaches may be used, such as GIGI for pedigree data. GIGI is an imputation method that handles large pedigrees and is particularly good for rare variant imputation. GIGI requires a subset of individuals in a pedigree to be fully sequenced, while other individuals are sequenced only at relevant markers. The imputation will infer the missing genotypes at untyped markers. Running GIGI on large pedigrees for large numbers of markers can be very time consuming. We present GIGI-Quick as a method to efficiently split GIGI's input, run GIGI in parallel and efficiently merge the output to reduce the runtime with the number of cores. This allows obtaining imputation results faster, and therefore all subsequent association analyses.

Availability and and implementation

GIGI-Quick is open source and publicly available via: https://cse-git.qcri.org/Imputation/GIGI-Quick.

Contact

msaad@hbku.edu.qa.

Supplementary information

Supplementary data are available at Bioinformatics online.

Association study of SLC6A2 gene Thr99Ile variant (rs1805065) with athletic status in the Brazilian population

Genetic variants in monoamine neurotransmitter genes have been recurrently associated with panic disorder, addiction and mood disorders. Recent evidence also indicates that norepinephrine neurotransmission can influence a series of psychophysical and psychobiological parameters related to athletic performance, and the presence of variants in the SLC6A2 (solute carrier family 6 member 2) gene, which encodes the norepinephrine transporter, can be detrimental to an adequate noradrenergic signaling. Accordingly, the objective of the present study was to explore the SLC6A2 Thr99Ile variant (rs1805065) in a cohort composed of highly-trained individuals and non-trained individuals. A total of 1556 Brazilians: 926 non-athletes and 630 athletes (322 endurance athletes and 308 power athletes) were compared in this case-control association study. The Thr99Ile variant showed only two genotypes (C/C or C/T), and a low minor allele frequency of ≈1%. However, none of the power athletes had the mutant T-allele (i.e., the C/T genotype), which may be related to decreased norepinephrine transporter activity. The genotype distribution and allele frequency observed in power athletes were significantly different when compared to non-athletes or endurance athletes. Therefore, the presence of the T-allele may decrease the chance of belonging to the group of athletes involved in explosive physical tasks. These results still need to be replicated in independent cohorts. However, it appears reasonable to assume that there is an association between the SLC6A2 gene variant and power athletic status.

Dynamic Scan Procedure for Detecting Rare-Variant Association Regions in Whole-Genome Sequencing Studies

Whole-genome sequencing (WGS) studies are being widely conducted in order to identify rare variants associated with human diseases and disease-related traits. Classical single-marker association analyses for rare variants have limited power, and variant-set-based analyses are commonly used by researchers for analyzing rare variants. However, existing variant-set-based approaches need to pre-specify genetic regions for analysis; hence, they are not directly applicable to WGS data because of the large number of intergenic and intron regions that consist of a massive number of non-coding variants. The commonly used sliding-window method requires the pre-specification of fixed window sizes, which are often unknown as a priori, are difficult to specify in practice, and are subject to limitations given that the sizes of genetic-association regions are likely to vary across the genome and phenotypes. We propose a computationally efficient and dynamic scan-statistic method (Scan the Genome [SCANG]) for analyzing WGS data; this method flexibly detects the sizes and the locations of rare-variant association regions without the need to specify a prior, fixed window size. The proposed method controls for the genome-wise type I error rate and accounts for the linkage disequilibrium among genetic variants. It allows the detected sizes of rare-variant association regions to vary across the genome. Through extensive simulated studies that consider a wide variety of scenarios, we show that SCANG substantially outperforms several alternative methods for detecting rare-variant-associations while controlling for the genome-wise type I error rates. We illustrate SCANG by analyzing the WGS lipids data from the Atherosclerosis Risk in Communities (ARIC) study.

Rediscovering the value of families for psychiatric genetics research

As it is likely that both common and rare genetic variation are important for complex disease risk, studies that examine the full range of the allelic frequency distribution should be utilized to dissect the genetic influences on mental illness. The rate limiting factor for inferring an association between a variant and a phenotype is inevitably the total number of copies of the minor allele captured in the studied sample. For rare variation, with minor allele frequencies of 0.5% or less, very large samples of unrelated individuals are necessary to unambiguously associate a locus with an illness. Unfortunately, such large samples are often cost prohibitive. However, by using alternative analytic strategies and studying related individuals, particularly those from large multiplex families, it is possible to reduce the required sample size while maintaining statistical power. We contend that using whole genome sequence (WGS) in extended pedigrees provides a cost-effective strategy for psychiatric gene mapping that complements common variant approaches and WGS in unrelated individuals. This was our impetus for forming the "Pedigree-Based Whole Genome Sequencing of Affective and Psychotic Disorders" consortium. In this review, we provide a rationale for the use of WGS with pedigrees in modern psychiatric genetics research. We begin with a focused review of the current literature, followed by a short history of family-based research in psychiatry. Next, we describe several advantages of pedigrees for WGS research, including power estimates, methods for studying the environment, and endophenotypes. We conclude with a brief description of our consortium and its goals.