A hierarchical and modular approach to the discovery of robust associations in genome-wide association studies from pooled DNA samples

Genetic Variation and Sickle Cell Disease Severity: A Systematic Review and Meta-Analysis

Importance

Sickle cell disease (SCD) is a monogenic disorder, yet clinical outcomes are influenced by additional genetic factors. Despite decades of research, the genetics of SCD remain poorly understood.

Objective

To assess all reported genetic modifiers of SCD, evaluate the design of associated studies, and provide guidelines for future analyses according to modern genetic study recommendations.

Data Sources

PubMed, Web of Science, and Scopus were searched through May 16, 2023, identifying 5290 publications.

Study Selection

At least 2 reviewers identified 571 original, peer-reviewed English-language publications reporting genetic modifiers of human SCD phenotypes, wherein the outcome was not treatment response, and the comparison was not between SCD subtypes or including healthy controls.

Data Extraction and Synthesis

Data relevant to all genetic modifiers of SCD were extracted, evaluated, and presented following STREGA and PRISMA guidelines. Weighted z score meta-analyses and pathway analyses were conducted.

Main Outcomes and Measures

Outcomes were aggregated into 25 categories, grouped as acute complications, chronic conditions, hematologic parameters or biomarkers, and general or mixed measures of SCD severity.

Results

The 571 included studies reported on 29 670 unique individuals (50% ≤ 18 years of age) from 43 countries. Of the 17 757 extracted results (4890 significant) in 1552 genes, 3675 results met the study criteria for meta-analysis: reported phenotype and genotype, association size and direction, variability measure, sample size, and statistical test. Only 173 results for 62 associations could be cross-study combined. The remaining associations could not be aggregated because they were only reported once or methods (eg, study design, reporting practice) and genotype or phenotype definitions were insufficiently harmonized. Gene variants regulating fetal hemoglobin and α-thalassemia (important markers for SCD severity) were frequently identified: 19 single-nucleotide variants in BCL11A, HBS1L-MYB, and HBG2 were significantly associated with fetal hemoglobin (absolute value of Z = 4.00 to 20.66; P = 8.63 × 10-95 to 6.19 × 10-5), and α-thalassemia deletions were significantly associated with increased hemoglobin level and reduced risk of albuminuria, abnormal transcranial Doppler velocity, and stroke (absolute value of Z = 3.43 to 5.16; P = 2.42 × 10-7 to 6.00 × 10-4). However, other associations remain unconfirmed. Pathway analyses of significant genes highlighted the importance of cellular adhesion, inflammation, oxidative and toxic stress, and blood vessel regulation in SCD (23 of the top 25 Gene Ontology pathways involve these processes) and suggested future research areas.

Conclusions and Relevance

The findings of this comprehensive systematic review and meta-analysis of all published genetic modifiers of SCD indicated that implementation of standardized phenotypes, statistical methods, and reporting practices should accelerate discovery and validation of genetic modifiers and development of clinically actionable genetic profiles.

Association of the ATG9B gene polymorphisms with coronary artery disease susceptibility: A case-control study

Introduction: Endothelial nitric oxide synthase (eNOS), the main regulator of cardiac cell functioning, is regulated post-transcriptionally by autophagy-related 9B (ATG9B) gene. The proper function of the heart is partly determined by the intact interaction of these molecules. The present study aimed to investigate the effects of ATG9B rs2373929 and rs7830 gene polymorphisms on the predisposition to coronary artery disease (CAD).

Methods: In this hospital-based case-control study, 150 patients with CAD compared with 150 healthy subjects for the genotype distributions of rs2373929 and rs7830 polymorphisms using T-ARMS PCR and ARMS PCR, respectively.

Results: Considering rs2373929 polymorphism, increased risk of CAD observed in the presence of TT genotype (OR: 3.65; 95% CI: 1.77-7.53; P < 0.001) and also in the recessive model for T allele (OR: 3.41; 95% CI: 1.76- 6.60; P < 0.001). The frequency of the T allele was higher in cases compared to controls (OR: 1.71; 95% CI: 1.24-2.28; P = 0.001). The genotype and allele frequencies of the rs7830 polymorphism did not differ between the two study groups.

Conclusion: The ATG9B gene rs2373929 polymorphism might involve in the pathogenesis of the CAD and can be considered as a screening molecular marker in the subjects prone to CAD.

A genome-wide association study of essential hypertension in an Australian population using a DNA pooling approach

Despite the success of genome-wide association studies (GWAS) in detecting genetic loci involved in complex traits, few susceptibility genes have been detected for essential hypertension (EH). We aimed to use pooled DNA GWAS approach to identify and validate novel genomic loci underlying EH susceptibility in an Australian case-control population. Blood samples and questionnaires detailing medical history, blood pressure, and prescribed medications were collected for 409 hypertensives and 409 age-, sex- and ethnicity-matched normotensive controls. Case and control DNA were pooled in quadruplicate and hybridized to Illumina 1 M-Duo arrays. Allele frequencies agreed with those reported in reference data and known EH association signals were represented in the top-ranked SNPs more frequently than expected by chance. Validation showed that pooled DNA GWAS gave reliable estimates of case and control allele frequencies. Although no markers reached Bonferroni-corrected genome-wide significance levels (5.0 × 10^-8), the top marker rs34870220 near ASGR1 approached significance (p = 4.32 × 10^-7), as did several candidate loci (p < 1 × 10^-6) on chromosomes 2, 4, 6, 9, 12, and 17. Four markers (located in or near genes NHSL1, NKFB1, GLI2, and LRRC10) from the top ten ranked SNPs were individually genotyped in pool samples and were tested for association between cases and controls using the χ ² test. Of these, rs1599961 (NFKB1) and rs12711538 (GLI2) showed significant difference between cases and controls (p < 0.01). Additionally, four top-ranking markers within NFKB1 were found to be in LD, suggesting a single strong association signal for this gene.

Increased circulating PEDF and low sICAM-1 are associated with sickle cell retinopathy

Sickle cell retinopathy (SCR) develops in up to 30% of sickle cell disease patients (SCD) during the second decade of life. Treatment for this affection remains palliative, so studies on its pathophysiology may contribute to the future development of novel therapies. SCR is more frequently observed in hemoglobin SC disease and derives from vaso-occlusion in the microvasculature of the retina leading to neovascularization and, eventually, to blindness. Circulating inflammatory cytokines, angiogenic factors, and their interaction may contribute to the pathophysiology of this complication. Angiopoietin (Ang)-1, Ang-2, soluble vascular cell adhesion molecule-1, intercellular adhesion molecule (ICAM)-1, E-selectin, P-selectin, IL1-β, TNF-α, pigment epithelium derived factor (PEDF) and vascular endothelial growth factor plasmatic levels were determined in 37 SCD patients with retinopathy, 34 without retinopathy, and healthy controls. We observed that sICAM-1 is significantly decreased, whereas PEDF is elevated in HbSC patients with retinopathy (P=0.012 and P=0.031, respectively). Ang-1, Ang-2 and IL1-β levels were elevated in SCD patients (P=0.001, P<0.001 and P=0.001, respectively), compared to controls, and HbSS patients presented higher levels of Ang-2 compared to HbSC (P<0.001). Our study supports the possible influence of sICAM-1 and PEDF on the pathophysiology of retinal neovascularization in SCD patients.

Bayesian methods for multivariate modeling of pleiotropic SNP associations and genetic risk prediction

Genome-wide association studies (GWAS) have identified numerous associations between genetic loci and individual phenotypes; however, relatively few GWAS have attempted to detect pleiotropic associations, in which loci are simultaneously associated with multiple distinct phenotypes. We show that pleiotropic associations can be directly modeled via the construction of simple Bayesian networks, and that these models can be applied to produce single or ensembles of Bayesian classifiers that leverage pleiotropy to improve genetic risk prediction. The proposed method includes two phases: (1) Bayesian model comparison, to identify Single-Nucleotide Polymorphisms (SNPs) associated with one or more traits; and (2) cross-validation feature selection, in which a final set of SNPs is selected to optimize prediction. To demonstrate the capabilities and limitations of the method, a total of 1600 case-control GWAS datasets with two dichotomous phenotypes were simulated under 16 scenarios, varying the association strengths of causal SNPs, the size of the discovery sets, the balance between cases and controls, and the number of pleiotropic causal SNPs. Across the 16 scenarios, prediction accuracy varied from 90 to 50%. In the 14 scenarios that included pleiotropically associated SNPs, the pleiotropic model search and prediction methods consistently outperformed the naive model search and prediction. In the two scenarios in which there were no true pleiotropic SNPs, the differences between the pleiotropic and naive model searches were minimal. To further evaluate the method on real data, a discovery set of 1071 sickle cell disease (SCD) patients was used to search for pleiotropic associations between cerebral vascular accidents and fetal hemoglobin level. Classification was performed on a smaller validation set of 352 SCD patients, and showed that the inclusion of pleiotropic SNPs may slightly improve prediction, although the difference was not statistically significant. The proposed method is robust, computationally efficient, and provides a powerful new approach for detecting and modeling pleiotropic disease loci.

A single nucleotide polymorphism of the adenosine deaminase, RNA-specific gene is associated with the serum triglyceride level, abdominal circumference, and serum adiponectin concentration

BACKGROUND

Single nucleotide polymorphisms (SNPs) of the adenosine deaminase, RNA-specific (ADAR) gene were reported to be associated with human longevity. There are possibilities that ADAR is associated with major risk factors of atherosclerotic cardiovascular diseases (CVD), such as hypertension, diabetes, dyslipidemia, and obesity.

OBJECTIVE

To investigate the association between SNPs of the ADAR gene and clinical data associated with major risk factors of atherosclerotic CVD.

SUBJECTS

A total of 1504 general population residents (586 males and 918 females) of two towns, Tanno-cho and Sobestu-cho, in Hokkaido, Japan.

METHODS

Clinical data associated with risk factors of atherosclerotic CVD were collected from these study subjects. DNA from peripheral blood and written informed consent were obtained. Three single nucleotide polymorphisms of ADARB1 and ADARB2, which were previously reported to be associated with longevity, were genotyped employing the TaqMan PCR method. The associations between SNPs in ADARB1 and ADARB2 and clinical parameters related to risk factors of atherosclerosis were analyzed.

RESULTS

On uni- and multivariate analyses, rs2805533 in ADARB2 was significantly associated with the abdominal circumference, body mass index, serum triglyceride level, and serum adiponectin level. The subjects with the AA genotype of rs2805533 had a greater abdominal circumference, higher body mass index, higher triglyceride level, and lower adiponectin level than those with AG and GG genotypes.

CONCLUSION

The SNP in ADARB2 related to longevity is associated with metabolic disorders. This finding suggests that genetic factors modulate human longevity via the regulation of metabolic factors such as abdominal obesity and lipid profiles.

Scanning for genomic regions subject to selective sweeps using SNP-MaP strategy

Population genomic approaches, which take advantages of high-throughput genotyping, are powerful yet costly methods to scan for selective sweeps. DNA-pooling strategies have been widely used for association studies because it is a cost-effective alternative to large-scale individual genotyping. Here, we performed an SNP-MaP (single nucleotide polymorphism microarrays and pooling) analysis using samples from Eurasia to evaluate the efficiency of pooling strategy in genome-wide scans for selection. By conducting simulations of allelotype data, we first demonstrated that the boxplot with average heterozygosity (HET) is a promising method to detect strong selective sweeps with a moderate level of pooling error. Based on this, we used a sliding window analysis of HET to detect the large contiguous regions (LCRs) putatively under selective sweeps from Eurasia datasets. This survey identified 63 LCRs in a European population. These signals were further supported by the integrated haplotype score (iHS) test using HapMap II data. We also confirmed the European-specific signatures of positive selection from several previously identified genes(KEL, TRPV5, TRPV6, EPHB6). In summary, our results not only revealed the high credibility of SNP-MaP strategy in scanning for selective sweeps, but also provided an insight into the population differentiation.

RNA surveillance-an emerging role for RNA regulatory networks in aging

In this review, we describe recent advances in the field of RNA regulatory biology and relate these advances to aging science. We introduce a new term, RNA surveillance, an RNA regulatory process that is conserved in metazoans, and describe how RNA surveillance represents molecular cross-talk between two emerging RNA regulatory systems-RNA interference and RNA editing. We discuss how RNA surveillance mechanisms influence mRNA and microRNA expression and activity during lifespan. Additionally, we summarize recent data from our own laboratory linking the RNA editor, ADAR, with exceptional longevity in humans and lifespan in Caenorhabditis elegans. We present data showing that transcriptional knockdown of RNA interference restores lifespan losses in the context of RNA editing defects, further suggesting that interaction between these two systems influences lifespan. Finally, we discuss the implications of RNA surveillance for sarcopenia and muscle maintenance, as frailty is a universal feature of aging. We end with a discussion of RNA surveillance as a robust regulatory system that can change in response to environmental stressors and represents a novel axis in aging science.

Pooled genome-wide analysis to identify novel risk loci for pediatric allergic asthma

Background

Genome-wide association studies of pooled DNA samples were shown to be a valuable tool to identify candidate SNPs associated to a phenotype. No such study was up to now applied to childhood allergic asthma, even if the very high complexity of asthma genetics is an appropriate field to explore the potential of pooled GWAS approach.

Methodology/Principal Findings

We performed a pooled GWAS and individual genotyping in 269 children with allergic respiratory diseases comparing allergic children with and without asthma. We used a modular approach to identify the most significant loci associated with asthma by combining silhouette statistics and physical distance method with cluster-adapted thresholding. We found 97% concordance between pooled GWAS and individual genotyping, with 36 out of 37 top-scoring SNPs significant at individual genotyping level. The most significant SNP is located inside the coding sequence of C5, an already identified asthma susceptibility gene, while the other loci regulate functions that are relevant to bronchial physiopathology, as immune- or inflammation-mediated mechanisms and airway smooth muscle contraction. Integration with gene expression data showed that almost half of the putative susceptibility genes are differentially expressed in experimental asthma mouse models.

Conclusion/Significance

Combined silhouette statistics and cluster-adapted physical distance threshold analysis of pooled GWAS data is an efficient method to identify candidate SNP associated to asthma development in an allergic pediatric population.

Severe sickle cell anemia is associated with increased plasma levels of TNF-R1 and VCAM-1

Sickle cell anemia (SCA, HBB glu6val) is characterized by multiple complications and a high degree of phenotypic variability: some subjects have only sporadic pain crises and few acute hospitalizations, while others experience multiple serious complications, high levels of morbidity, and accelerated mortality [1]. The tumor necrosis factor-α (TNF-α) signaling pathway plays important roles in inflammation and the immune response; variation in this pathway might be expected to modify the overall severity of SCA through the pathway's effects on the vascular endothelium [2,3]. We examined plasma biomarkers of TNF-α activity and endothelial cell activation for associations with SCA severity in 24 adults (12 mild, 12 severe). Two biomarkers, tumor necrosis factor-α receptor-1 (TNF-R1) and vascular cell adhesion molecule-1 (VCAM-1) were significantly higher in subjects with severe SCA. Along with these biomarker differences, we also examined data from a genome-wide association study (GWAS) using SCA severity as a disease phenotype, and found evidence of genetic association between disease severity and a single nucleotide polymorphism (SNP) in VCAM1, which codes for VCAM-1, and several SNPs in ARFGEF2, a gene involved in TNF-R1 release [4].

Identification of a common variant affecting human episodic memory performance using a pooled genome-wide association approach: a case study of disease gene identification

Genome-wide association studies (GWAS) are an important tool for discovering novel genes associated with disease or traits. Careful design of case-control groups greatly facilitates the efficacy of these studies. Here we describe a pooled GWAS study undertaken to find novel genes associated with human episodic memory performance. A genomic locus for the WW and C2 domain-containing 1 protein, KIBRA (also known as WWC1), was found to be associated with memory performance in three cognitively normal cohorts from Switzerland and the USA. This result was further supported by correlation of KIBRA genotype and differences in hippocampal activation as measured by functional magnetic resonance imaging (fMRI). These findings provide an excellent example of the application of GWAS using a pooled genomic DNA approach to successfully identify a locus with strong effects on human memory.

Use of pathway information in molecular epidemiology

Candidate gene studies are generally motivated by some form of pathway reasoning in the selection of genes to be studied, but seldom has the logic of the approach been carried through to the analysis. Marginal effects of polymorphisms in the selected genes, and occasionally pairwise gene-gene or gene-environment interactions, are often presented, but a unified approach to modelling the entire pathway has been lacking. In this review, a variety of approaches to this problem is considered, focusing on hypothesis-driven rather than purely exploratory methods. Empirical modelling strategies are based on hierarchical models that allow prior knowledge about the structure of the pathway and the various reactions to be included as 'prior covariates'. By contrast, mechanistic models aim to describe the reactions through a system of differential equations with rate parameters that can vary between individuals, based on their genotypes. Some ways of combining the two approaches are suggested and Bayesian model averaging methods for dealing with uncertainty about the true model form in either framework is discussed. Biomarker measurements can be incorporated into such analyses, and two-phase sampling designs stratified on some combination of disease, genes and exposures can be an efficient way of obtaining data that would be too expensive or difficult to obtain on a full candidate gene sample. The review concludes with some thoughts about potential uses of pathways in genome-wide association studies.

Genetic modifiers of the severity of sickle cell anemia identified through a genome-wide association study

We conducted a genome-wide association study (GWAS) to discover single nucleotide polymorphisms (SNPs) associated with the severity of sickle cell anemia in 1,265 patients with either "severe" or "mild" disease based on a network model of disease severity. We analyzed data using single SNP analysis and a novel SNP set enrichment analysis (SSEA) developed to discover clusters of associated SNPs. Single SNP analysis discovered 40 SNPs that were strongly associated with sickle cell severity (odds for association >1,000); of the 32 that we could analyze in an independent set of 163 patients, five replicated, eight showed consistent effects although failed to reach statistical significance, whereas 19 did not show any convincing association. Among the replicated associations are SNPs in KCNK6 a K(+) channel gene. SSEA identified 27 genes with a strong enrichment of significant SNPs (P < 10(-6)); 20 were replicated with varying degrees of confidence. Among the novel findings identified by SSEA is the telomere length regulator gene TNKS. These studies are the first to use GWAS to understand the genetic diversity that accounts the phenotypic heterogeneity sickle cell anemia as estimated by an integrated model of severity. Additional validation, resequencing, and functional studies to understand the biology and reveal mechanisms by which candidate genes might have their effects are the future goals of this work.

RNA editing genes associated with extreme old age in humans and with lifespan in C. elegans

BACKGROUND

The strong familiality of living to extreme ages suggests that human longevity is genetically regulated. The majority of genes found thus far to be associated with longevity primarily function in lipoprotein metabolism and insulin/IGF-1 signaling. There are likely many more genetic modifiers of human longevity that remain to be discovered.

METHODOLOGY/PRINCIPAL FINDINGS

Here, we first show that 18 single nucleotide polymorphisms (SNPs) in the RNA editing genes ADARB1 and ADARB2 are associated with extreme old age in a U.S. based study of centenarians, the New England Centenarian Study. We describe replications of these findings in three independently conducted centenarian studies with different genetic backgrounds (Italian, Ashkenazi Jewish and Japanese) that collectively support an association of ADARB1 and ADARB2 with longevity. Some SNPs in ADARB2 replicate consistently in the four populations and suggest a strong effect that is independent of the different genetic backgrounds and environments. To evaluate the functional association of these genes with lifespan, we demonstrate that inactivation of their orthologues adr-1 and adr-2 in C. elegans reduces median survival by 50%. We further demonstrate that inactivation of the argonaute gene, rde-1, a critical regulator of RNA interference, completely restores lifespan to normal levels in the context of adr-1 and adr-2 loss of function.

CONCLUSIONS/SIGNIFICANCE

Our results suggest that RNA editors may be an important regulator of aging in humans and that, when evaluated in C. elegans, this pathway may interact with the RNA interference machinery to regulate lifespan.

Methodological Issues in Multistage Genome-wide Association Studies

Because of the high cost of commercial genotyping chip technologies, many investigations have used a two-stage design for genome-wide association studies, using part of the sample for an initial discovery of "promising" SNPs at a less stringent significance level and the remainder in a joint analysis of just these SNPs using custom genotyping. Typical cost savings of about 50% are possible with this design to obtain comparable levels of overall type I error and power by using about half the sample for stage I and carrying about 0.1% of SNPs forward to the second stage, the optimal design depending primarily upon the ratio of costs per genotype for stages I and II. However, with the rapidly declining costs of the commercial panels, the generally low observed ORs of current studies, and many studies aiming to test multiple hypotheses and multiple endpoints, many investigators are abandoning the two-stage design in favor of simply genotyping all available subjects using a standard high-density panel. Concern is sometimes raised about the absence of a "replication" panel in this approach, as required by some high-profile journals, but it must be appreciated that the two-stage design is not a discovery/replication design but simply a more efficient design for discovery using a joint analysis of the data from both stages. Once a subset of highly-significant associations has been discovered, a truly independent "exact replication" study is needed in a similar population of the same promising SNPs using similar methods. This can then be followed by (1) "generalizability" studies to assess the full scope of replicated associations across different races, different endpoints, different interactions, etc.; (2) fine-mapping or re-sequencing to try to identify the causal variant; and (3) experimental studies of the biological function of these genes. Multistage sampling designs may be more useful at this stage, say for selecting subsets of subjects for deep re-sequencing of regions identified in the GWAS.

Genome-wide association studies and the genetic dissection of complex traits

The availability of affordable high throughput technology for parallel genotyping has opened the field of genetics to genome-wide association studies (GWAS), and in the last few years hundreds of articles reporting results of GWAS for a variety of heritable traits have been published. What do these results tell us? Although GWAS have discovered a few hundred reproducible associations, this number is underwhelming in relation to the huge amount of data produced, and challenges the conjecture that common variants may be the genetic causes of common diseases. We argue that the massive amount of genetic data that result from these studies remains largely unexplored and unexploited because of the challenge of mining and modeling enormous data sets, the difficulty of using nontraditional computational techniques and the focus of accepted statistical analyses on controlling the false positive rate rather than limiting the false negative rate. In this article, we will review the common approach to analysis of GWAS data and then discuss options to learn more from these data. We will use examples from our ongoing studies of sickle cell anemia and also GWAS in multigenic traits.

Molecular genetics of adult ADHD: converging evidence from genome-wide association and extended pedigree linkage studies

A genome-wide association (GWA) study with pooled DNA in adult attention-deficit/hyperactivity disorder (ADHD) employing approximately 500K SNP markers identifies novel risk genes and reveals remarkable overlap with findings from recent GWA scans in substance use disorders. Comparison with results from our previously reported high-resolution linkage scan in extended pedigrees confirms several chromosomal loci, including 16q23.1-24.3 which also reached genome-wide significance in a recent meta-analysis of seven linkage studies (Zhou et al. in Am J Med Genet Part B, 2008). The findings provide additional support for a common effect of genes coding for cell adhesion molecules (e.g., CDH13, ASTN2) and regulators of synaptic plasticity (e.g., CTNNA2, KALRN) despite the complex multifactorial etiologies of adult ADHD and addiction vulnerability.