An Excess of Risk-Increasing Low-Frequency Variants Can Be a Signal of Polygenic Inheritance in Complex Diseases

In most complex diseases, much of the heritability remains unaccounted for by common variants. It has been postulated that lower-frequency variants contribute to the remaining heritability. Here, we describe a method to test for polygenic inheritance from lower-frequency variants by using GWAS summary association statistics. We explored scenarios with many causal low-frequency variants and showed that there is more power to detect risk variants than to detect protective variants, resulting in an increase in the ratio of detected risk to protective variants (R/P ratio). Such an excess can also occur if risk variants are present and kept at lower frequencies because of negative selection. The R/P ratio can be falsely elevated because of reasons unrelated to polygenic inheritance, such as uneven sample sizes or asymmetric population stratification, so precautions to correct for these confounders are essential. We tested our method on published GWAS results and observed a strong signal in some diseases (schizophrenia and type 2 diabetes) but not others. We also explored the shared genetic component in overlapping phenotypes related to inflammatory bowel disease (Crohn disease [CD] and ulcerative colitis [UC]) and diabetic nephropathy (macroalbuminuria and end-stage renal disease [ESRD]). Although the signal was still present when both CD and UC were jointly analyzed, the signal was lost when macroalbuminuria and ESRD were jointly analyzed, suggesting that these phenotypes should best be studied separately. Thus, our method may also help guide the design of future genetic studies of various traits and diseases.

Annotated Chromosome Maps for Renal Disease

A combination of linkage analyses and association studies are currently employed to promote the identification of genetic factors contributing to inherited renal disease. We have standardized and merged complex genetic data from disparate sources, creating unique chromosomal maps to enhance genetic epidemiological investigations. This database and novel renal maps effectively summarize genomic regions of suggested linkage, association, or chromosomal abnormalities implicated in renal disease. Chromosomal regions associated with potential intermediate clinical phenotypes have been integrated, adding support for particular genomic intervals. More than 500 reports from medical databases, published scientific literature, and the World Wide Web were interrogated for relevant renal-related information. Chromosomal regions highlighted for prioritized investigation of renal complications include 3q13-26, 6q22-27, 10p11-15, 16p11-13, and 18q22. Combined genetic and physical maps are effective tools to organize genetic data for complex diseases. These renal chromosome maps provide insights into renal phenotype-genotype relationships and act as a template for future genetic investigations into complex renal diseases. New data from individual researchers and/or future publications can be readily incorporated to this resource via a user-friendly web-form accessed from the website: www.qub.ac.uk/neph-res/CORGI/index.php.

Genetic association analyses of non-synonymous single nucleotide polymorphisms in diabetic nephropathy

Aims/hypothesis: Diabetic nephropathy, characterised by persistent proteinuria, hypertension and progressive kidney failure, affects a subset of susceptible individuals with diabetes. It is also a leading cause of end-stage renal disease (ESRD). Non-synonymous (ns) single nucleotide polymorphisms (SNPs) have been reported to contribute to genetic susceptibility in both monogenic disorders and common complex diseases. The objective of this study was to investigate whether nsSNPs are involved in susceptibility to diabetic nephropathy using a case-control design.

Methods: White type 1 diabetic patients with (cases) and without (controls) nephropathy from eight centres in the UK and Ireland were genotyped for a selected subset of nsSNPs using Illumina's GoldenGate BeadArray assay. A ? 2 test for trend, stratified by centre, was used to assess differences in genotype distribution between cases and controls. Genomic control was used to adjust for possible inflation of test statistics, and the False Discovery Rate method was used to account for multiple testing.

Results: We assessed 1,111 nsSNPs for association with diabetic nephropathy in 1,711 individuals with type 1 diabetes (894 cases, 817 controls). A number of SNPs demonstrated a significant difference in genotype distribution between groups before but not after correction for multiple testing. Furthermore, neither subgroup analysis (diabetic nephropathy with ESRD or diabetic nephropathy without ESRD) nor stratification by duration of diabetes revealed any significant differences between groups.

Conclusions/interpretation: The nsSNPs investigated in this study do not appear to contribute significantly to the development of diabetic nephropathy in patients with type 1 diabetes.

A Variant in LDLR is Associated with Abdominal Aortic Aneurysm

Background—Abdominal aortic aneurysm (AAA) is a common cardiovascular disease among older people and demonstrates significant heritability. In contrast to similar complex diseases, relatively few genetic associations with AAA have been confirmed. We reanalysed our genome-wide study and carried through to replication suggestive discovery associations at a lower level of significance.

Methods and Results—A genome-wide association study was conducted using 1,830 cases from the UK, New Zealand and Australia with infra-renal aorta diameter =30mm or ruptured AAA and 5,435 unscreened controls from the 1958 Birth Cohort and National Blood Service cohort from the Wellcome Trust Case Control Consortium. Eight suggestive associations with P<1x10-4 were carried through to in silico replication in 1,292 AAA cases and 30,503 controls. One SNP associated with P<0.05 after Bonferroni correction in the in silico study underwent further replication (706 AAA cases and 1,063 controls from the UK, 507 AAA cases and 199 controls from Denmark and 885 AAA cases and 1,000 controls from New Zealand). Low density lipoprotein receptor (LDLR) rs6511720 A, was significantly associated overall and in three of five individual replication studies. The full study showed an association that reached genome-wide significance (odds ratio 0.76; 95% confidence interval 0.70 to 0.83; P=2.08x10-10).

Conclusions—LDLR rs6511720 is associated with abdominal aortic aneurysm. This finding is consistent with established effects of this variant on coronary artery disease. Shared aetiological pathways with other cardiovascular diseases may present novel opportunities for preventative and therapeutic strategies for AAA.

Harnessing the complexity of gene expression data from cancer: from single gene to structural pathway methods

High-dimensional gene expression data provide a rich source of information because they capture the expression level of genes in dynamic states that reflect the biological functioning of a cell. For this reason, such data are suitable to reveal systems related properties inside a cell, e.g., in order to elucidate molecular mechanisms of complex diseases like breast or prostate cancer. However, this is not only strongly dependent on the sample size and the correlation structure of a data set, but also on the statistical hypotheses tested. Many different approaches have been developed over the years to analyze gene expression data to (I) identify changes in single genes, (II) identify changes in gene sets or pathways, and (III) identify changes in the correlation structure in pathways. In this paper, we review statistical methods for all three types of approaches, including subtypes, in the context of cancer data and provide links to software implementations and tools and address also the general problem of multiple hypotheses testing. Further, we provide recommendations for the selection of such analysis methods.

DNA hypermethylation and DNA hypomethylation is present at different loci in chronic kidney disease

Genetic risk factors for chronic kidney disease (CKD) are being identified through international collaborations. By comparison, epigenetic risk factors for CKD have only recently been considered using population-based approaches. DNA methylation is a major epigenetic modification that is associated with complex diseases, so we investigated methylome-wide loci for association with CKD. A total of 485,577 unique features were evaluated in 255 individuals with CKD (cases) and 152 individuals without evidence of renal disease (controls). Following stringent quality control, raw data were quantile normalized and β values calculated to reflect the methylation status at each site. The difference in methylation status was evaluated between cases and controls with resultant P values adjusted for multiple testing. Genes with significantly increased and decreased levels of DNA methylation were considered for biological relevance by functional enrichment analysis using KEGG pathways in Partek Genomics Suite. Twenty-three genes, where more than one CpG per loci was identified with Padjusted < 10−8, demonstrated significant methylation changes associated with CKD and additional support for these associated loci was sought from published literature. Strong biological candidates for CKD that showed statistically significant differential methylation include CUX1, ELMO1, FKBP5, INHBA-AS1, PTPRN2, and PRKAG2 genes; several genes are differentially methylated in kidney tissue and RNA-seq supports a functional role for differential methylation in ELMO1 and PRKAG2 genes. This study reports the largest, most comprehensive, genome-wide quantitative evaluation of DNA methylation for association with CKD. Evidence confirming methylation sites influence development of CKD would stimulate research to identify epigenetic therapies that might be clinically useful for CKD.

Novel genetic susceptibility loci for diabetic end-stage renal disease identified through robust naive Bayes classification

Aims/hypothesis: Diabetic nephropathy is a major diabetic complication, and diabetes is the leading cause of end-stage renal disease (ESRD). Family studies suggest a hereditary component for diabetic nephropathy. However, only a few genes have been associated with diabetic nephropathy or ESRD in diabetic patients. Our aim was to detect novel genetic variants associated with diabetic nephropathy and ESRD. Methods: We exploited a novel algorithm, ‘Bag of Naive Bayes’, whose marker selection strategy is complementary to that of conventional genome-wide association models based on univariate association tests. The analysis was performed on a genome-wide association study of 3,464 patients with type 1 diabetes from the Finnish Diabetic Nephropathy (FinnDiane) Study and subsequently replicated with 4,263 type 1 diabetes patients from the Steno Diabetes Centre, the All Ireland-Warren 3-Genetics of Kidneys in Diabetes UK collection (UK–Republic of Ireland) and the Genetics of Kidneys in Diabetes US Study (GoKinD US). Results: Five genetic loci (WNT4/ZBTB40-rs12137135, RGMA/MCTP2-rs17709344, MAPRE1P2-rs1670754, SEMA6D/SLC24A5-rs12917114 and SIK1-rs2838302) were associated with ESRD in the FinnDiane study. An association between ESRD and rs17709344, tagging the previously identified rs12437854 and located between the RGMA and MCTP2 genes, was replicated in independent case–control cohorts. rs12917114 near SEMA6D was associated with ESRD in the replication cohorts under the genotypic model (p < 0.05), and rs12137135 upstream of WNT4 was associated with ESRD in Steno. Conclusions/interpretation: This study supports the previously identified findings on the RGMA/MCTP2 region and suggests novel susceptibility loci for ESRD. This highlights the importance of applying complementary statistical methods to detect novel genetic variants in diabetic nephropathy and, in general, in complex diseases.

Do persistent organic pollutants interact with the stress response? Individual compounds, and their mixtures, interaction with the glucocorticoid receptor

Persistent organic pollutants (POPs) are toxic substances, highly resistant to environmental degradation, which can bio-accumulate and have long-range atmospheric transport potential (UNEP 2001). The majority of studies on endocrine disruption have focused on interferences on the sexual steroid hormones and so have overlooked disruption to glucocorticoid hormones. Here the endocrine disrupting potential of individual POPs and their mixtures has been investigated in vitro to identify any disruption to glucocorticoid nuclear receptor transcriptional activity. POP mixtures were screened for glucocorticoid receptor (GR) translocation using a GR redistribution assay (RA) on a CellInsight(TM) NXT High Content Screening (HCS) platform. A mammalian reporter gene assay (RGA) was then used to assess the individual POPs, and their mixtures, for effects on glucocorticoid nuclear receptor transactivation. POP mixtures did not induce GR translocation in the GR RA or produce an agonist response in the GR RGA. However, in the antagonist test, in the presence of cortisol, an individual POP, p,p'-dichlorodiphenyldichloroethylene (DDE), was found to decrease glucocorticoid nuclear receptor transcriptional activity to 72.5% (in comparison to the positive cortisol control). Enhanced nuclear transcriptional activity, in the presence of cortisol, was evident for the two lowest concentrations of perfluorodecanoic acid (PFOS) potassium salt (0.0147mg/ml and 0.0294mg/ml), the two highest concentrations of perfluorodecanoic acid (PFDA) (0.0025mg/ml and 0.005mg/ml) and the highest concentration of 2,2',4,4'-tetrabromodiphenyl ether (BDE 47) (0.0000858mg/ml). It is important to gain a better understanding of how POPs can interact with GRs as the disruption of glucocorticoid action is thought to contribute to complex diseases.

Inflammation markers are associated with Cardiovascular diseases risk in adolescents The Young Hearts Project 2000

Purose: The traditional approach for identifying subjects at risk from cardiovascular diseases (CVD) is to determine the extent of clustering of biological risk factors adjusted for lifestyle. Recently, markers of endothelial dysfunction and low grade inflammation, including high sensitivity C-reactive protein (hsCRP), soluble intercellular adhesion molecules (sICAM), and soluble vascular adhesion molecules (sVCAM), have been included in the detection for high risk individuals. However, the relationship of these novel biomarkers with CVD risk in adolescents remains unclear. The purpose of this study, therefore, was to establish the association of hsCRP, sICAM, and sVCAM with CVD risk in an adolescent population.
Methods: Data from the Young Hearts 2000 cross-sectional cohort study, carried out in 1999-2001, were used. From a total of 2,017 male and female participants, 95 obese subjects were identified and matched according to age, sex, and cigarette smoking, with 95 overweight and 95 normal-weight adolescents. Clustered CVD risk was computed using a sum of Z-scores of biological risk factors. The relationship was described using multiple linear regression analyses.
Results: hsCRP, sICAM, and sVCAM showed significant associations with CVD risk. hsCRP and sICAM had a positive relation with CVD risk, whereas sVCAM showed an inverse relationship. In this study, lifestyle factors showed no relation with CVD risk.
Conclusion: The results fit the hypothesized role of low grade inflammation and endothelial dysfunction in CVD risk in asymptomatic adolescents. The inverse relationship of VCAM, however, is hard to explain and indicates the complex mechanisms underlying CVD. Further research is needed to draw firm conclusions on the biomarkers used.

Identification of a general secretory pathway in a human isolate of Burkholderia vietnamiensis (formerly B. cepacia complex genomovar V) that is required for the secretion of hemolysin and phospholipase C activities

Members of the Burkholderia cepacia complex can secrete proteases, lipases, and hemolysins. We report in this study the identification of a general secretory pathway present in a B. vietnamiensis (formerly genomovar V) clinical isolate, which is required for the efficient secretion of phospholipase C and hemolysin activities. Southern blot hybridization experiments revealed that this general secretion pathway is highly conserved among the different genomovars of the B. cepacia complex and is homologous to a similar system described in B. pseudomallei. We also show that this pathway appears not to be necessary for intracellular survival of B. vietnamiensis within Acanthamoeba polyphaga.

A PCR-based method to distinguish fungi of the rice sheath-blight complex, Rhizoctonia solani, R. oryzae and R. oryzae-sativae

Identification of Rhizoctonia solani, R. oryzae and R. oryzae-sativae, components of the rice sheath disease complex, is extremely difficult and often inaccurate and as a result may hinder the success of extensive breeding programmes throughout Asia. In this study, primers designed from unique regions within the rDNA internal transcribed spacers have been used to develop a rapid PCR-based diagnostic test to provide an accurate identification of the species on rice. Tests on the specificity of the primers concerned showed that they provide the means for accurate identification of the Rhizoctonia species responsible for sheath diseases in rice.

Partitioning Heritability of Regulatory and Cell-Type-Specific Variants across 11 Common Diseases

Regulatory and coding variants are known to be enriched with associations identified by genome-wide association studies (GWASs) of complex disease, but their contributions to trait heritability are currently unknown. We applied variance-component methods to imputed genotype data for 11 common diseases to partition the heritability explained by genotyped SNPs () across functional categories (while accounting for shared variance due to linkage disequilibrium). Extensive simulations showed that in contrast to current estimates from GWAS summary statistics, the variance-component approach partitions heritability accurately under a wide range of complex-disease architectures. Across the 11 diseases DNaseI hypersensitivity sites (DHSs) from 217 cell types spanned 16% of imputed SNPs (and 24% of genotyped SNPs) but explained an average of 79% (SE = 8%) of  from imputed SNPs (5.1× enrichment; p = 3.7 × 10⁻¹⁷) and 38% (SE = 4%) of  from genotyped SNPs (1.6× enrichment, p = 1.0 × 10⁻⁴). Further enrichment was observed at enhancer DHSs and cell-type-specific DHSs. In contrast, coding variants, which span 1% of the genome, explained <10% of  despite having the highest enrichment. We replicated these findings but found no significant contribution from rare coding variants in independent schizophrenia cohorts genotyped on GWAS and exome chips. Our results highlight the value of analyzing components of heritability to unravel the functional architecture of common disease.

Clonal complex Pseudomonas aeruginosa in horses

Pseudomonas aeruginosa is associated with infectious endometritis in horses. Although infectious endometritis is often considered a venereal infection, there is relatively limited genotypic-based evidence to support this mode of transmission. The study sought to determine the relatedness between genital P. aeruginosa isolates collected from a limited geographical region using molecular strain typing. Enterobacterial repetitive intergenic consensus PCR typing was performed on 93 isolates collected between 2005 and 2009 from 2058 thoroughbred horses (including 18 stallions) at 66 studs. While P. aeruginosa was not detected in the stallions, 53/93 (57%) mares harbouring P. aeruginosa had clonally related strains, which included a single dominant genotype detected in 42 (45%) mares from 13 different studs. These novel findings suggest that most equine genital P. aeruginosa infections in this region may have been acquired from mechanisms other than direct horse to horse transmission. Instead, other potential acquisition pathways, as well as strain specific adaptation to the equine genital tract, should be investigated.