Ozren Polasek

SLC2A9 is a newly identified urate transporter influencing serum urate concentration, urate excretion and gout

Uric acid is the end product of purine metabolism in humans and great apes, which have lost hepatic uricase activity, leading to uniquely high serum uric acid concentrations (200–500 μM) compared with other mammals (3–120 μM). About 70% of daily urate disposal occurs via the kidneys, and in 5–25% of the human population, impaired renal excretion leads to hyperuricemia. About 10% of people with hyperuricemia develop gout, an inflammatory arthritis that results from deposition of monosodium urate crystals in the joint. We have identified genetic variants within a transporter gene, SLC2A9, that explain 1.7–5.3% of the variance in serum uric acid concentrations, following a genome-wide association scan in a Croatian population sample. SLC2A9 variants were also associated with low fractional excretion of uric acid and/or gout in UK, Croatian and German population samples. SLC2A9 is a known fructose transporter, and we now show that it has strong uric acid transport activity in Xenopus laevis oocytes.

A Comprehensive Evaluation of Potential Lung Function Associated Genes in the SpiroMeta General Population Sample

Rationale Lung function measures are heritable traits that predict population morbidity and mortality and are essential for the diagnosis of chronic obstructive pulmonary disease (COPD). Variations in many genes have been reported to affect these traits, but attempts at replication have provided conflicting results. Recently, we undertook a meta-analysis of Genome Wide Association Study (GWAS) results for lung function measures in 20,288 individuals from the general population (the SpiroMeta consortium). Objectives To comprehensively analyse previously reported genetic associations with lung function measures, and to investigate whether single nucleotide polymorphisms (SNPs) in these genomic regions are associated with lung function in a large population sample. Methods We analysed association for SNPs tagging 130 genes and 48 intergenic regions (+/−10 kb), after conducting a systematic review of the literature in the PubMed database for genetic association studies reporting lung function associations. Results The analysis included 16,936 genotyped and imputed SNPs. No loci showed overall significant association for FEV1 or FEV1/FVC traits using a carefully defined significance threshold of 1.3×10−5. The most significant loci associated with FEV1 include SNPs tagging MACROD2 (P = 6.81×10−5), CNTN5 (P = 4.37×10−4), and TRPV4 (P = 1.58×10−3). Among ever-smokers, SERPINA1 showed the most significant association with FEV1 (P = 8.41×10−5), followed by PDE4D (P = 1.22×10−4). The strongest association with FEV1/FVC ratio was observed with ABCC1 (P = 4.38×10−4), and ESR1 (P = 5.42×10−4) among ever-smokers. Conclusions Polymorphisms spanning previously associated lung function genes did not show strong evidence for association with lung function measures in the SpiroMeta consortium population. Common SERPINA1 polymorphisms may affect FEV1 among smokers in the general population.

High throughput isolation and glycosylation analysis of IgG - variability and heritability of the IgG glycome in three isolated human populations

All immunoglobulin G molecules carry N-glycans, which modulate their biological activity. Changes in N-glycosylation of IgG associate with various diseases and affect the activity of therapeutic antibodies and intravenous immunoglobulins. We have developed a novel 96-well protein G monolithic plate and used it to rapidly isolate IgG from plasma of 2298 individuals from three isolated human populations. N-glycans were released by PNGase F, labeled with 2-aminobenzamide and analyzed by hydrophilic interaction chromatography with fluorescence detection. The majority of the structural features of the IgG glycome were consistent with previous studies, but sialylation was somewhat higher than reported previously. Sialylation was particularly prominent in core fucosylated glycans containing two galactose residues and bisecting GlcNAc where median sialylation level was nearly 80%. Very high variability between individuals was observed, approximately three times higher than in the total plasma glycome. For example, neutral IgG glycans without core fucose varied between 1.3 and 19%, a difference that significantly affects the effector functions of natural antibodies, predisposing or protecting individuals from particular diseases. Heritability of IgG glycans was generally between 30 and 50%. The individuals age was associated with a significant decrease in galactose and increase of bisecting GlcNAc, whereas other functional elements of IgG glycosylation did not change much with age. Gender was not an important predictor for any IgG glycan. An important observation is that competition between glycosyltransferases, which occurs in vitro, did not appear to be relevant in vivo, indicating that the final glycan structures are not a simple result of competing enzymatic activities, but a carefully regulated outcome designed to meet the prevailing physiological needs.

A General Approach for Haplotype Phasing across the Full Spectrum of Relatedness

Many existing cohorts contain a range of relatedness between genotyped individuals, either by design or by chance. Haplotype estimation in such cohorts is a central step in many downstream analyses. Using genotypes from six cohorts from isolated populations and two cohorts from non-isolated populations, we have investigated the performance of different phasing methods designed for nominally ‘unrelated’ individuals. We find that SHAPEIT2 produces much lower switch error rates in all cohorts compared to other methods, including those designed specifically for isolated populations. In particular, when large amounts of IBD sharing is present, SHAPEIT2 infers close to perfect haplotypes. Based on these results we have developed a general strategy for phasing cohorts with any level of implicit or explicit relatedness between individuals. First SHAPEIT2 is run ignoring all explicit family information. We then apply a novel HMM method (duoHMM) to combine the SHAPEIT2 haplotypes with any family information to infer the inheritance pattern of each meiosis at all sites across each chromosome. This allows the correction of switch errors, detection of recombination events and genotyping errors. We show that the method detects numbers of recombination events that align very well with expectations based on genetic maps, and that it infers far fewer spurious recombination events than Merlin. The method can also detect genotyping errors and infer recombination events in otherwise uninformative families, such as trios and duos. The detected recombination events can be used in association scans for recombination phenotypes. The method provides a simple and unified approach to haplotype estimation, that will be of interest to researchers in the fields of human, animal and plant genetics.

Genome-wide association analyses identify multiple loci associated with central corneal thickness and keratoconus

Central corneal thickness (CCT) is associated with eye conditions including keratoconus and glaucoma. We performed a meta-analysis on >20,000 individuals in European and Asian populations that identified 16 new loci associated with CCT at genome-wide significance (P < 5 × 10−8). We further showed that 2 CCT-associated loci, FOXO1 and FNDC3B, conferred relatively large risks for keratoconus in 2 cohorts with 874 cases and 6,085 controls (rs2721051 near FOXO1 had odds ratio (OR) = 1.62, 95% confidence interval (CI) = 1.4–1.88, P = 2.7 × 10−10, and rs4894535 in FNDC3B had OR = 1.47, 95% CI = 1.29–1.68, P = 4.9 × 10−9). FNDC3B was also associated with primary open-angle glaucoma (P = 5.6 × 10−4; tested in 3 cohorts with 2,979 cases and 7,399 controls). Further analyses implicate the collagen and extracellular matrix pathways in the regulation of CCT.

Genome-Wide Association Study Identifies Novel Loci Associated with Circulating Phospho- and Sphingolipid Concentrations

Phospho- and sphingolipids are crucial cellular and intracellular compounds. These lipids are required for active transport, a number of enzymatic processes, membrane formation, and cell signalling. Disruption of their metabolism leads to several diseases, with diverse neurological, psychiatric, and metabolic consequences. A large number of phospholipid and sphingolipid species can be detected and measured in human plasma. We conducted a meta-analysis of five European family-based genome-wide association studies (N = 4034) on plasma levels of 24 sphingomyelins (SPM), 9 ceramides (CER), 57 phosphatidylcholines (PC), 20 lysophosphatidylcholines (LPC), 27 phosphatidylethanolamines (PE), and 16 PE-based plasmalogens (PLPE), as well as their proportions in each major class. This effort yielded 25 genome-wide significant loci for phospholipids (smallest P-value = 9.88×10−204) and 10 loci for sphingolipids (smallest P-value = 3.10×10−57). After a correction for multiple comparisons (P-value<2.2×10−9), we observed four novel loci significantly associated with phospholipids (PAQR9, AGPAT1, PKD2L1, PDXDC1) and two with sphingolipids (PLD2 and APOE) explaining up to 3.1% of the variance. Further analysis of the top findings with respect to within class molar proportions uncovered three additional loci for phospholipids (PNLIPRP2, PCDH20, and ABDH3) suggesting their involvement in either fatty acid elongation/saturation processes or fatty acid specific turnover mechanisms. Among those, 14 loci (KCNH7, AGPAT1, PNLIPRP2, SYT9, FADS1-2-3, DLG2, APOA1, ELOVL2, CDK17, LIPC, PDXDC1, PLD2, LASS4, and APOE) mapped into the glycerophospholipid and 12 loci (ILKAP, ITGA9, AGPAT1, FADS1-2-3, APOA1, PCDH20, LIPC, PDXDC1, SGPP1, APOE, LASS4, and PLD2) to the sphingolipid pathways. In large meta-analyses, associations between FADS1-2-3 and carotid intima media thickness, AGPAT1 and type 2 diabetes, and APOA1 and coronary artery disease were observed. In conclusion, our study identified nine novel phospho- and sphingolipid loci, substantially increasing our knowledge of the genetic basis for these traits.

Glycans Are a Novel Biomarker of Chronological and Biological Ages

Fine structural details of glycans attached to the conserved N-glycosylation site significantly not only affect function of individual immunoglobulin G (IgG) molecules but also mediate inflammation at the systemic level. By analyzing IgG glycosylation in 5,117 individuals from four European populations, we have revealed very complex patterns of changes in IgG glycosylation with age. Several IgG glycans (including FA2B, FA2G2, and FA2BG2) changed considerably with age and the combination of these three glycans can explain up to 58% of variance in chronological age, significantly more than other markers of biological age like telomere lengths. The remaining variance in these glycans strongly correlated with physiological parameters associated with biological age. Thus, IgG glycosylation appears to be closely linked with both chronological and biological ages. Considering the important role of IgG glycans in inflammation, and because the observed changes with age promote inflammation, changes in IgG glycosylation also seem to represent a factor contributing to aging. Significance Statement Glycosylation is the key posttranslational mechanism that regulates function of immunoglobulins, with multiple systemic repercussions to the immune system. Our study of IgG glycosylation in 5,117 individuals from four European populations has revealed very extensive and complex changes in IgG glycosylation with age. The combined index composed of only three glycans explained up to 58% of variance in age, considerably more than other biomarkers of age like telomere lengths. The remaining variance in these glycans strongly correlated with physiological parameters associated with biological age; thus, IgG glycosylation appears to be closely linked with both chronological and biological ages. The ability to measure human biological aging using molecular profiling has practical applications for diverse fields such as disease prevention and treatment, or forensics.

Comparative Performance of Four Methods for High-throughput Glycosylation Analysis of Immunoglobulin G in Genetic and Epidemiological Research

The biological and clinical relevance of glycosylation is becoming increasingly recognized, leading to a growing interest in large-scale clinical and population-based studies. In the past few years, several methods for high-throughput analysis of glycans have been developed, but thorough validation and standardization of these methods is required before significant resources are invested in large-scale studies. In this study, we compared liquid chromatography, capillary gel electrophoresis, and two MS methods for quantitative profiling of N-glycosylation of IgG in the same data set of 1201 individuals. To evaluate the accuracy of the four methods we then performed analysis of association with genetic polymorphisms and age. Chromatographic methods with either fluorescent or MS-detection yielded slightly stronger associations than MS-only and multiplexed capillary gel electrophoresis, but at the expense of lower levels of throughput. Advantages and disadvantages of each method were identified, which should inform the selection of the most appropriate method in future studies.

Estimating Human Inbreeding Coefficients: Comparison of Genealogical and Marker Heterozygosity Approaches

We have used genealogies and genomic polymorphisms to estimate individual inbreeding coefficients (F) in 50 subjects with an expected range (based on recent genealogies) of F from 0.0 to 0.0625. The estimates were based on two approaches, using genotypes respectively from 410 microsatellite markers (410‐STR panel) and from 10,000 SNPs (10K‐SNP panel). The latter was performed in a sub‐sample of 15 individuals. We concluded that for both marker panels measures of inbreeding based on the excess of homozygosity over Hardy‐Weinberg expectation were not closely correlated with 4‐5 generation genealogical F‐values. For the 10K‐SNP panel we found two alternative measures which correlated more closely with F, based respectively on standard errors and on paired homozygosity of nearby SNPs over distances of 2‐4 cM. We propose an empirical method for estimating standard errors and hence individual F‐values, based on the variation between individual autosomes. This method could provide useful estimates of average F‐values for groups of individuals in population‐based studies of the effects of inbreeding/homozygosity on quantitative traits.

Genome-wide meta-analysis of myopia and hyperopia provides evidence for replication of 11 loci

Refractive error (RE) is a complex, multifactorial disorder characterized by a mismatch between the optical power of the eye and its axial length that causes object images to be focused off the retina. The two major subtypes of RE are myopia (nearsightedness) and hyperopia (farsightedness), which represent opposite ends of the distribution of the quantitative measure of spherical refraction. We performed a fixed effects meta-analysis of genome-wide association results of myopia and hyperopia from 9 studies of European-derived populations: AREDS, KORA, FES, OGP-Talana, MESA, RSI, RSII, RSIII and ERF. One genome-wide significant region was observed for myopia, corresponding to a previously identified myopia locus on 8q12 (p = 1.25×10−8), which has been reported by Kiefer et al. as significantly associated with myopia age at onset and Verhoeven et al. as significantly associated to mean spherical-equivalent (MSE) refractive error. We observed two genome-wide significant associations with hyperopia. These regions overlapped with loci on 15q14 (minimum p value = 9.11×10−11) and 8q12 (minimum p value 1.82×10−11) previously reported for MSE and myopia age at onset. We also used an intermarker linkage- disequilibrium-based method for calculating the effective number of tests in targeted regional replication analyses. We analyzed myopia (which represents the closest phenotype in our data to the one used by Kiefer et al.) and showed replication of 10 additional loci associated with myopia previously reported by Kiefer et al. This is the first replication of these loci using myopia as the trait under analysis. “Replication-level” association was also seen between hyperopia and 12 of Kiefer et al.s published loci. For the loci that show evidence of association to both myopia and hyperopia, the estimated effect of the risk alleles were in opposite directions for the two traits. This suggests that these loci are important contributors to variation of refractive error across the distribution.