Genetics of osteopontin in patients with chronic kidney disease: The German Chronic Kidney Disease study.

Osteopontin (OPN), encoded by SPP1, is a phosphorylated glycoprotein predominantly synthesized in kidney tissue. Increased OPN mRNA and protein expression correlates with proteinuria, reduced creatinine clearance, and kidney fibrosis in animal models of kidney disease. But its genetic underpinnings are incompletely understood. We therefore conducted a genome-wide association study (GWAS) of OPN in a European chronic kidney disease (CKD) population. Using data from participants of the German Chronic Kidney Disease (GCKD) study (N = 4,897), a GWAS (minor allele frequency [MAF]≥1%) and aggregated variant testing (AVT, MAF<1%) of ELISA-quantified serum OPN, adjusted for age, sex, estimated glomerular filtration rate (eGFR), and urinary albumin-to-creatinine ratio (UACR) was conducted. In the project, GCKD participants had a mean age of 60 years (SD 12), median eGFR of 46 mL/min/1.73m2 (p25: 37, p75: 57) and median UACR of 50 mg/g (p25: 9, p75: 383). GWAS revealed 3 loci (p<5.0E-08), two of which replicated in the population-based Young Finns Study (YFS) cohort (p<1.67E-03): rs10011284, upstream of SPP1 encoding the OPN protein and related to OPN production, and rs4253311, mapping into KLKB1 encoding prekallikrein (PK), which is processed to kallikrein (KAL) implicated through the kinin-kallikrein system (KKS) in blood pressure control, inflammation, blood coagulation, cancer, and cardiovascular disease. The SPP1 gene was also identified by AVT (p = 2.5E-8), comprising 7 splice-site and missense variants. Among others, downstream analyses revealed colocalization of the OPN association signal at SPP1 with expression in pancreas tissue, and at KLKB1 with various plasma proteins in trans, and with phenotypes (bone disorder, deep venous thrombosis) in human tissue. In summary, this GWAS of OPN levels revealed two replicated associations. The KLKB1 locus connects the function of OPN with PK, suggestive of possible further post-translation processing of OPN. Further studies are needed to elucidate the complex role of OPN within human (patho)physiology.

The membrane transporter SLC25A48 enables transport of choline into human mitochondria.

Choline has important physiological functions as a precursor for essential cell components, signaling molecules, phospholipids, and the neurotransmitter acetylcholine. Choline is a water-soluble charged molecule requiring transport proteins to cross biological membranes. Although transporters continue to be identified, membrane transport of choline is incompletely understood and knowledge about choline transport into intracellular organelles such as mitochondria remains limited. Here we show that SLC25A48 imports choline into human mitochondria. Human loss-of-function mutations in SLC25A48 show impaired choline transport into mitochondria and are associated with elevated urine and plasma choline levels. Thus, our studies may have implications for understanding and treating conditions related to choline metabolism.

pgainsim: an R-package to assess the mode of inheritance for quantitative trait loci in GWAS.

MOTIVATION: When performing genome-wide association studies conventionally the additive genetic model is used to explore whether a single nucleotide polymorphism (SNP) is associated with a quantitative trait. But for variants, which do not follow an intermediate mode of inheritance (MOI), the recessive or the dominant genetic model can have more power to detect associations and furthermore the MOI is important for downstream analyses and clinical interpretation. When multiple MOIs are modelled the question arises, which describes the true underlying MOI best. RESULTS: We developed an R-package allowing for the first time to determine study specific critical values when one of the three models is more informative than the other ones for a quantitative trait locus. The package allows for user-friendly simulations to determine these critical values with predefined minor allele frequencies and study sizes. For application scenarios with extensive multiple testing we integrated an interpolation functionality to determine critical values already based on a moderate number of random draws. AVAILABILITY AND IMPLEMENTATION: The R-package pgainsim is freely available for download on Github at https://github.com/genepi-freiburg/pgainsim. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

New insights into the hypothalamic-pituitary-thyroid axis: a transcriptome- and proteome-wide association study.

Introduction: Thyroid hormones have systemic effects on the human body and play a key role in the development and function of virtually all tissues. They are regulated via the hypothalamic-pituitary-thyroid (HPT) axis and have a heritable component. Using genetic information, we applied tissue-specific transcriptome-wide association studies (TWAS) and plasma proteome-wide association studies (PWAS) to elucidate gene products related to thyrotropin (TSH) and free thyroxine (FT4) levels. Results: TWAS identified 297 and 113 transcripts associated with TSH and FT4 levels, respectively (25 shared), including transcripts not identified by genome-wide association studies (GWAS) of these traits, demonstrating the increased power of this approach. Testing for genetic colocalization revealed a shared genetic basis of 158 transcripts with TSH and 45 transcripts with FT4, including independent, FT4-associated genetic signals within the CAPZB locus that were differentially associated with CAPZB expression in different tissues. PWAS identified 18 and ten proteins associated with TSH and FT4, respectively (HEXIM1 and QSOX2 with both). Among these, the cognate genes of five TSH- and 7 FT4-associated proteins mapped outside significant GWAS loci. Colocalization was observed for five plasma proteins each with TSH and FT4. There were ten TSH and one FT4-related gene(s) significant in both TWAS and PWAS. Of these, ANXA5 expression and plasma annexin A5 levels were inversely associated with TSH (PWAS: P = 1.18 × 10-13, TWAS: P = 7.61 × 10-12 (whole blood), P = 6.40 × 10-13 (hypothalamus), P = 1.57 × 10-15 (pituitary), P = 4.27 × 10-15 (thyroid)), supported by colocalizations. Conclusion: Our analyses revealed new thyroid function-associated genes and prioritized candidates in known GWAS loci, contributing to a better understanding of transcriptional regulation and protein levels relevant to thyroid function.

SLC17A1/3 transporters mediate renal excretion of Lac-Phe in mice and humans.

N-lactoyl-phenylalanine (Lac-Phe) is a lactate-derived metabolite that suppresses food intake and body weight. Little is known about the mechanisms that mediate Lac-Phe transport across cell membranes. Here we identify SLC17A1 and SLC17A3, two kidney-restricted plasma membrane-localized solute carriers, as physiologic urine Lac-Phe transporters. In cell culture, SLC17A1/3 exhibit high Lac-Phe efflux activity. In humans, levels of Lac-Phe in urine exhibit a strong genetic association with the SLC17A1-4 locus. Urine Lac-Phe levels are increased following a Wingate sprint test. In mice, genetic ablation of either SLC17A1 or SLC17A3 reduces urine Lac-Phe levels. Despite these differences, both knockout strains have normal blood Lac-Phe and body weights, demonstrating SLC17A1/3-dependent de-coupling of urine and plasma Lac-Phe pools. Together, these data establish SLC17A1/3 family members as the physiologic urine Lac-Phe transporters and uncover a biochemical pathway for the renal excretion of this signaling metabolite.

SLC17 transporters mediate renal excretion of Lac-Phe in mice and humans.

N-lactoyl-phenylalanine (Lac-Phe) is a lactate-derived metabolite that suppresses food intake and body weight. Little is known about the mechanisms that mediate Lac-Phe transport across cell membranes. Here we identify SLC17A1 and SLC17A3, two kidney-restricted plasma membrane-localized solute carriers, as physiologic urine Lac-Phe transporters. In cell culture, SLC17A1/3 exhibit high Lac-Phe efflux activity. In humans, levels of Lac-Phe in urine exhibit a strong genetic association with the SLC17A1-4 locus. Urine Lac-Phe levels are also increased following a Wingate sprint test. In mice, genetic ablation of either SLC17A1 or SLC17A3 reduces urine Lac-Phe levels. Despite these differences, both knockout strains have normal blood Lac-Phe and body weights, demonstrating that urine and plasma Lac-Phe pools are functionally and biochemically de-coupled. Together, these data establish SLC17 family members as the physiologic urine transporters for Lac-Phe and uncover a biochemical pathway for the renal excretion of this signaling metabolite.

SLC26A1 is a major determinant of sulfate homeostasis in humans.

Sulfate plays a pivotal role in numerous physiological processes in the human body, including bone and cartilage health. A role of the anion transporter SLC26A1 (Sat1) for sulfate reabsorption in the kidney is supported by the observation of hyposulfatemia and hypersulfaturia in Slc26a1-knockout mice. The impact of SLC26A1 on sulfate homeostasis in humans remains to be defined. By combining clinical genetics, functional expression assays, and population exome analysis, we identify SLC26A1 as a sulfate transporter in humans and experimentally validate several loss-of-function alleles. Whole-exome sequencing from a patient presenting with painful perichondritis, hyposulfatemia, and renal sulfate wasting revealed a homozygous mutation in SLC26A1, which has not been previously described to the best of our knowledge. Whole-exome data analysis of more than 5,000 individuals confirmed that rare, putatively damaging SCL26A1 variants were significantly associated with lower plasma sulfate at the population level. Functional expression assays confirmed a substantial reduction in sulfate transport for the SLC26A1 mutation of our patient, which we consider to be novel, as well as for the additional variants detected in the population study. In conclusion, combined evidence from 3 complementary approaches supports SLC26A1 activity as a major determinant of sulfate homeostasis in humans. In view of recent evidence linking sulfate homeostasis with back pain and intervertebral disc disorder, our study identifies SLC26A1 as a potential target for modulation of musculoskeletal health.

Genetic studies of paired metabolomes reveal enzymatic and transport processes at the interface of plasma and urine.

The kidneys operate at the interface of plasma and urine by clearing molecular waste products while retaining valuable solutes. Genetic studies of paired plasma and urine metabolomes may identify underlying processes. We conducted genome-wide studies of 1,916 plasma and urine metabolites and detected 1,299 significant associations. Associations with 40% of implicated metabolites would have been missed by studying plasma alone. We detected urine-specific findings that provide information about metabolite reabsorption in the kidney, such as aquaporin (AQP)-7-mediated glycerol transport, and different metabolomic footprints of kidney-expressed proteins in plasma and urine that are consistent with their localization and function, including the transporters NaDC3 (SLC13A3) and ASBT (SLC10A2). Shared genetic determinants of 7,073 metabolite-disease combinations represent a resource to better understand metabolic diseases and revealed connections of dipeptidase 1 with circulating digestive enzymes and with hypertension. Extending genetic studies of the metabolome beyond plasma yields unique insights into processes at the interface of body compartments.

Genome-wide studies reveal factors associated with circulating uromodulin and its relationships to complex diseases.

Uromodulin (UMOD) is a major risk gene for monogenic and complex forms of kidney disease. The encoded kidney-specific protein uromodulin is highly abundant in urine and related to chronic kidney disease, hypertension, and pathogen defense. To gain insights into potential systemic roles, we performed genome-wide screens of circulating uromodulin using complementary antibody-based and aptamer-based assays. We detected 3 and 10 distinct significant loci, respectively. Integration of antibody-based results at the UMOD locus with functional genomics data (RNA-Seq, ATAC-Seq, Hi-C) of primary human kidney tissue highlighted an upstream variant with differential accessibility and transcription in uromodulin-synthesizing kidney cells as underlying the observed cis effect. Shared association patterns with complex traits, including chronic kidney disease and blood pressure, placed the PRKAG2 locus in the same pathway as UMOD. Experimental validation of the third antibody-based locus, B4GALNT2, showed that the p.Cys466Arg variant of the encoded N-acetylgalactosaminyltransferase had a loss-of-function effect leading to higher serum uromodulin levels. Aptamer-based results pointed to enzymes writing glycan marks present on uromodulin and to their receptors in the circulation, suggesting that this assay permits investigating uromodulin's complex glycosylation rather than its quantitative levels. Overall, our study provides insights into circulating uromodulin and its emerging functions.