Large-scale plasma proteomics comparisons through genetics and disease associations.

High-throughput proteomics platforms measuring thousands of proteins in plasma combined with genomic and phenotypic information have the power to bridge the gap between the genome and diseases. Here we performed association studies of Olink Explore 3072 data generated by the UK Biobank Pharma Proteomics Project1 on plasma samples from more than 50,000 UK Biobank participants with phenotypic and genotypic data, stratifying on British or Irish, African and South Asian ancestries. We compared the results with those of a SomaScan v4 study on plasma from 36,000 Icelandic people2, for 1,514 of whom Olink data were also available. We found modest correlation between the two platforms. Although cis protein quantitative trait loci were detected for a similar absolute number of assays on the two platforms (2,101 on Olink versus 2,120 on SomaScan), the proportion of assays with such supporting evidence for assay performance was higher on the Olink platform (72% versus 43%). A considerable number of proteins had genomic associations that differed between the platforms. We provide examples where differences between platforms may influence conclusions drawn from the integration of protein levels with the study of diseases. We demonstrate how leveraging the diverse ancestries of participants in the UK Biobank helps to detect novel associations and refine genomic location. Our results show the value of the information provided by the two most commonly used high-throughput proteomics platforms and demonstrate the differences between them that at times provides useful complementarity.

Actionable Genotypes and Their Association with Life Span in Iceland.

BACKGROUND: In 2021, the American College of Medical Genetics and Genomics (ACMG) recommended reporting actionable genotypes in 73 genes associated with diseases for which preventive or therapeutic measures are available. Evaluations of the association of actionable genotypes in these genes with life span are currently lacking. METHODS: We assessed the prevalence of coding and splice variants in genes on the ACMG Secondary Findings, version 3.0 (ACMG SF v3.0), list in the genomes of 57,933 Icelanders. We assigned pathogenicity to all reviewed variants using reported evidence in the ClinVar database, the frequency of variants, and their associations with disease to create a manually curated set of actionable genotypes (variants). We assessed the relationship between these genotypes and life span and further examined the specific causes of death among carriers. RESULTS: Through manual curation of 4405 sequence variants in the ACMG SF v3.0 genes, we identified 235 actionable genotypes in 53 genes. Of the 57,933 participants, 2306 (4.0%) carried at least one actionable genotype. We found shorter median survival among persons carrying actionable genotypes than among noncarriers. Specifically, we found that carrying an actionable genotype in a cancer gene was associated with survival that was 3 years shorter than that among noncarriers, with causes of death among carriers attributed primarily to cancer-related conditions. Furthermore, we found evidence of association between carrying an actionable genotype in certain genes in the cardiovascular disease group and a reduced life span. CONCLUSIONS: On the basis of the ACMG SF v3.0 guidelines, we found that approximately 1 in 25 Icelanders carried an actionable genotype and that carrying such a genotype was associated with a reduced life span. (Funded by deCODE Genetics-Amgen.).

Reduced anticoagulation variability in patients on warfarin monitored with Fiix-prothrombin time associates with reduced thromboembolism: The Fiix-trial.

Fiix-prothrombin time (Fiix-PT) differs from traditional PT in being affected by reduced factor (F) II or FX only. In the randomized controlled Fiix-trial, patients on warfarin monitored with Fiix-PT (Fiix-warfarin patients) had fewer thromboembolisms (TE), similar major bleeding (MB) and more stable anticoagulation than patients monitored with PT (PT-warfarin patients). In the current Fiix-trial report we analyzed how reduced anticoagulation variability during Fiix-PT monitoring was reflected in patients with TE or bleeding. Data from 1143 randomized patients was used. We analyzed the groups for anticoagulation intensity (time within target range; TTR), international normalized ratio (INR) variability (variance growth rate B1; VGR) and dose adjustment frequency. We assessed how these parameters associated with clinically relevant vascular events (CRVE), ie TE or MB or clinically relevant non-MB. TTR was highest in Fiix-warfarin patients without CRVE (median 82%;IQR 72-91) and lowest in PT-warfarin patients with TE (62%;56-81). VGR was lowest in Fiix-warfarin patients without CRVE (median VGR B1 0.17; 95% CI 0.08-0.38) and with TE (0.20;0.07-0.26) and highest in PT-warfarin patients with TE (0.50;0.27-0.90) or MB (0.59;0.07-1.36). The mean annual dose adjustment frequency was lowest in Fiix-warfarin patients with TE (mean 5.4;95% CI 3.9-7.3) and without CRVE (mean 6.0; 5.8-6.2) and highest in PT-warfarin patients with TE (14.2;12.2-16.3). Frequent dose changes predicted MB in both study arms. Compared to patients monitored with PT, high anticoagulation stability in Fiix-warfarin patients coincided with their low TE rate. Those with bleeding had high variability irrespective of monitoring method. Thus, although further improvements are needed to reduce bleeding, stabilization of anticoagulation by Fiix-PT monitoring associates with reduced TE.

Basal cells of the human airways acquire mesenchymal traits in idiopathic pulmonary fibrosis and in culture.

Idiopathic pulmonary fibrosis (IPF) is a progressive interstitial lung disease with high morbidity and mortality. The cellular source of the fibrotic process is currently under debate with one suggested mechanism being epithelial-to-mesenchymal transition (EMT) in the alveolar region. In this study, we show that airway epithelium overlying fibroblastic foci in IPF contains a layer of p63-positive basal cells while lacking ciliated and goblet cells. This basal epithelium shows increased expression of CK14, Vimentin and N-cadherin while retaining E-cadherin. The underlying fibroblastic foci shows both E- and N-cadherin-positive cells. To determine if p63-positive basal cells were able to undergo EMT in culture, we treated VA10, a p63-positive basal cell line, with the serum replacement UltroserG. A sub-population of treated cells acquired a mesenchymal phenotype, including an E- to N-cadherin switch. After isolation, these cells portrayed a phenotype presenting major hallmarks of EMT (loss of epithelial markers, gain of mesenchymal markers, increased migration and anchorage-independent growth). This phenotypic switch was prevented in p63 knockdown (KD) cells. In conclusion, we show that airway epithelium overlying fibroblastic foci in IPF lacks its characteristic functional identity, shows increased reactivity of basal cells and acquisition of a partial EMT phenotype. This study suggests that some p63-positive basal cells are prone to phenotypic changes and could act as EMT progenitors in IPF.

Clinical phenotype in heterozygote and biallelic Bernard-Soulier syndrome--a case control study.

Bernard-Soulier syndrome (BSS) is a rare severe autosomal recessive bleeding disorder. To date heterozygous carriers of BSS mutations have not been shown to have bleeding symptoms. We assessed bleeding using a semi-quantitative questionnaire, platelet parameters, PFA-100 closure times, ristocetin response, GP Ib/IX expression and VWF antigen in 14 BSS patients, 30 heterozygote carriers for related mutations and 29 controls. Eight mutations in GP1BA, GP1BB or GP9 were identified including four previously unknown pathogenic mutations. Subjects with BSS reported markedly more mucocutaneous bleeding than controls. Increased bleeding was also observed in heterozygotes. Compared to controls, patients with BSS had lower optical platelet counts (P < 0.001), CD61-platelet counts (P < 0.001) and higher mean platelet volume (17.7 vs. 7.8 fL, P < 0.001) and ristocetin response and closure times were unmeasurable. Heterozygotes had higher MPV (9.7 fL, P < 0.001) and lower platelet counts (P < 0.001) than controls but response to ristocetin and closure times were normal. The VWF was elevated in both BSS and in heterozygotes (P = 0.005). We conclude that heterozygotes for BSS mutations have lower platelet counts than controls and show a bleeding phenotype albeit much milder than in BSS. Both patients with BSS and heterozygote carriers of pathogenic mutations have raised VWF.

Risk of recurrence of gastrointestinal stromal tumour after surgery: an analysis of pooled population-based cohorts.

BACKGROUND: The risk of recurrence of gastrointestinal stromal tumour (GIST) after surgery needs to be estimated when considering adjuvant systemic therapy. We assessed prognostic factors of patients with operable GIST, to compare widely used risk-stratification schemes and to develop a new method for risk estimation. METHODS: Population-based cohorts of patients diagnosed with operable GIST, who were not given adjuvant therapy, were identified from the literature. Data from ten series and 2560 patients were pooled. Risk of tumour recurrence was stratified using the National Institute of Health (NIH) consensus criteria, the modified consensus criteria, and the Armed Forces Institute of Pathology (AFIP) criteria. Prognostic factors were examined using proportional hazards and non-linear models. The results were validated in an independent centre-based cohort consisting of 920 patients with GIST. FINDINGS: Estimated 15-year recurrence-free survival (RFS) after surgery was 59·9% (95% CI 56·2-63·6); few recurrences occurred after the first 10 years of follow-up. Large tumour size, high mitosis count, non-gastric location, presence of rupture, and male sex were independent adverse prognostic factors. In receiver operating characteristics curve analysis of 10-year RFS, the NIH consensus criteria, modified consensus criteria, and AFIP criteria resulted in an area under the curve (AUC) of 0·79 (95% CI 0·76-0·81), 0·78 (0·75-0·80), and 0·82 (0·80-0·85), respectively. The modified consensus criteria identified a single high-risk group. Since tumour size and mitosis count had a non-linear association with the risk of GIST recurrence, novel prognostic contour maps were generated using non-linear modelling of tumour size and mitosis count, and taking into account tumour site and rupture. The non-linear model accurately predicted the risk of recurrence (AUC 0·88, 0·86-0·90). INTERPRETATION: The risk-stratification schemes assessed identify patients who are likely to be cured by surgery alone. Although the modified NIH classification is the best criteria to identify a single high-risk group for consideration of adjuvant therapy, the prognostic contour maps resulting from non-linear modelling are appropriate for estimation of individualised outcomes. FUNDING: Academy of Finland, Cancer Society of Finland, Sigrid Juselius Foundation and Helsinki University Research Funds.

Genetics and epidemiology of mutational barcode-defined clonal hematopoiesis.

Clonal hematopoiesis (CH) arises when a substantial proportion of mature blood cells is derived from a single hematopoietic stem cell lineage. Using whole-genome sequencing of 45,510 Icelandic and 130,709 UK Biobank participants combined with a mutational barcode method, we identified 16,306 people with CH. Prevalence approaches 50% in elderly participants. Smoking demonstrates a dosage-dependent impact on risk of CH. CH associates with several smoking-related diseases. Contrary to published claims, we find no evidence that CH is associated with cardiovascular disease. We provide evidence that CH is driven by genes that are commonly mutated in myeloid neoplasia and implicate several new driver genes. The presence and nature of a driver mutation alters the risk profile for hematological disorders. Nevertheless, most CH cases have no known driver mutations. A CH genome-wide association study identified 25 loci, including 19 not implicated previously in CH. Splicing, protein and expression quantitative trait loci were identified for CD164 and TCL1A.

Genetic architecture of band neutrophil fraction in Iceland.

The characteristic lobulated nuclear morphology of granulocytes is partially determined by composition of nuclear envelope proteins. Abnormal nuclear morphology is primarily observed as an increased number of hypolobulated immature neutrophils, called band cells, during infection or in rare envelopathies like Pelger-Huët anomaly. To search for sequence variants affecting nuclear morphology of granulocytes, we performed a genome-wide association study using band neutrophil fraction from 88,101 Icelanders. We describe 13 sequence variants affecting band neutrophil fraction at nine loci. Five of the variants are at the Lamin B receptor (LBR) locus, encoding an inner nuclear membrane protein. Mutations in LBR are linked to Pelger-Huët anomaly. In addition, we identify cosegregation of a rare stop-gain sequence variant in LBR and Pelger Huët anomaly in an Icelandic eight generation pedigree, initially reported in 1963. Two of the other loci include genes which, like LBR, play a role in the nuclear membrane function and integrity. These GWAS results highlight the role proteins of the inner nuclear membrane have as important for neutrophil nuclear morphology.

Prospective evaluation of quality of life in total versus supracervical laparoscopic hysterectomy.

STUDY OBJECTIVE: To evaluate and compare recovery times and quality of life (QOL) in patients undergoing a total laparoscopic hysterectomy (TLH) and laparoscopic supracervical hysterectomy (LSH). DESIGN: Patients underwent either a TLH or LSH. After surgery, patients maintained a daily log documenting pain, nausea, use of pain medications, and return to daily activities. They also completed a QOL questionnaire (SF-36) before and after surgery. DESIGN CLASSIFICATION: Prospective cohort study (Canadian Task Force Classification II-1). SETTING: University teaching hospital. PATIENTS: A total of 122 women undergoing laparoscopic hysterectomy. MEASUREMENTS AND MAIN RESULTS: A total of 122 women underwent TLH (n = 71) or LSH (n = 51) for benign indications from February 2008 to January 2010. There was a significantly higher postoperative improvement of QOL scores in the LSH group in 6 of 10 questionnaire categories and summary scores, including physical functioning (p =.03), role physical (p =.002), and bodily pain (p =.03). There were no significant differences in use of pain medications, level of pain, level of nausea, or return to normal activities. CONCLUSION: LSH appears to provide greater improvement in short-term postoperative QOL compared with TLH. No significant differences were found in postoperative pain or return to daily activities.

Genetic variants associated with platelet count are predictive of human disease and physiological markers.

Platelets play an important role in hemostasis and other aspects of vascular biology. We conducted a meta-analysis of platelet count GWAS using data on 536,974 Europeans and identified 577 independent associations. To search for mechanisms through which these variants affect platelets, we applied cis-expression quantitative trait locus, DEPICT and IPA analyses and assessed genetic sharing between platelet count and various traits using polygenic risk scoring. We found genetic sharing between platelet count and counts of other blood cells (except red blood cells), in addition to several other quantitative traits, including markers of cardiovascular, liver and kidney functions, height, and weight. Platelet count polygenic risk score was predictive of myeloproliferative neoplasms, rheumatoid arthritis, ankylosing spondylitis, hypertension, and benign prostate hyperplasia. Taken together, these results advance understanding of diverse aspects of platelet biology and how they affect biological processes in health and disease.

Eleven genomic loci affect plasma levels of chronic inflammation marker soluble urokinase-type plasminogen activator receptor.

Soluble urokinase-type plasminogen activator receptor (suPAR) is a chronic inflammation marker associated with the development of a range of diseases, including cancer and cardiovascular disease. The genetics of suPAR remain unexplored but may shed light on the biology of the marker and its connection to outcomes. We report a heritability estimate of 60% for the variation in suPAR and performed a genome-wide association meta-analysis on suPAR levels measured in Iceland (N = 35,559) and in Denmark (N = 12,177). We identified 13 independently genome-wide significant sequence variants associated with suPAR across 11 distinct loci. Associated variants were found in and around genes encoding uPAR (PLAUR), its ligand uPA (PLAU), the kidney-disease-associated gene PLA2R1 as well as genes with relations to glycosylation, glycoprotein biosynthesis, and the immune response. These findings provide new insight into the causes of variation in suPAR plasma levels, which may clarify suPAR's potential role in associated diseases, as well as the underlying mechanisms that give suPAR its prognostic value as a unique marker of chronic inflammation.

Large-scale integration of the plasma proteome with genetics and disease.

The plasma proteome can help bridge the gap between the genome and diseases. Here we describe genome-wide association studies (GWASs) of plasma protein levels measured with 4,907 aptamers in 35,559 Icelanders. We found 18,084 associations between sequence variants and levels of proteins in plasma (protein quantitative trait loci; pQTL), of which 19% were with rare variants (minor allele frequency (MAF) < 1%). We tested plasma protein levels for association with 373 diseases and other traits and identified 257,490 associations. We integrated pQTL and genetic associations with diseases and other traits and found that 12% of 45,334 lead associations in the GWAS Catalog are with variants in high linkage disequilibrium with pQTL. We identified 938 genes encoding potential drug targets with variants that influence levels of possible biomarkers. Combining proteomics, genomics and transcriptomics, we provide a valuable resource that can be used to improve understanding of disease pathogenesis and to assist with drug discovery and development.

A genome-wide meta-analysis yields 46 new loci associating with biomarkers of iron homeostasis.

Iron is essential for many biological functions and iron deficiency and overload have major health implications. We performed a meta-analysis of three genome-wide association studies from Iceland, the UK and Denmark of blood levels of ferritin (N = 246,139), total iron binding capacity (N = 135,430), iron (N = 163,511) and transferrin saturation (N = 131,471). We found 62 independent sequence variants associating with iron homeostasis parameters at 56 loci, including 46 novel loci. Variants at DUOX2, F5, SLC11A2 and TMPRSS6 associate with iron deficiency anemia, while variants at TF, HFE, TFR2 and TMPRSS6 associate with iron overload. A HBS1L-MYB intergenic region variant associates both with increased risk of iron overload and reduced risk of iron deficiency anemia. The DUOX2 missense variant is present in 14% of the population, associates with all iron homeostasis biomarkers, and increases the risk of iron deficiency anemia by 29%. The associations implicate proteins contributing to the main physiological processes involved in iron homeostasis: iron sensing and storage, inflammation, absorption of iron from the gut, iron recycling, erythropoiesis and bleeding/menstruation.

Airway branching morphogenesis in three dimensional culture.

BACKGROUND: Lungs develop from the fetal digestive tract where epithelium invades the vascular rich stroma in a process called branching morphogenesis. In organogenesis, endothelial cells have been shown to be important for morphogenesis and the maintenance of organ structure. The aim of this study was to recapitulate human lung morphogenesis in vitro by establishing a three dimensional (3D) co-culture model where lung epithelial cells were cultured in endothelial-rich stroma. METHODS: We used a human bronchial epithelial cell line (VA10) recently developed in our laboratory. This cell line cell line maintains a predominant basal cell phenotype, expressing p63 and other basal markers such as cytokeratin-5 and -14. Here, we cultured VA10 with human umbilical vein endothelial cells (HUVECs), to mimic the close interaction between these cell types during lung development. Morphogenesis and differentiation was monitored by phase contrast microscopy, immunostainings and confocal imaging. RESULTS: We found that in co-culture with endothelial cells, the VA10 cells generated bronchioalveolar like structures, suggesting that lung epithelial branching is facilitated by the presence of endothelial cells. The VA10 derived epithelial structures display various complex patterns of branching and show partial alveolar type-II differentiation with pro-Surfactant-C expression. The epithelial origin of the branching VA10 colonies was confirmed by immunostaining. These bronchioalveolar-like structures were polarized with respect to integrin expression at the cell-matrix interface. The endothelial-induced branching was mediated by soluble factors. Furthermore, fibroblast growth factor receptor-2 (FGFR-2) and sprouty-2 were expressed at the growing tips of the branching structures and the branching was inhibited by the FGFR-small molecule inhibitor SU5402. DISCUSSION: In this study we show that a human lung epithelial cell line can be induced by endothelial cells to form branching bronchioalveolar-like structures in 3-D culture. This novel model of human airway morphogenesis can be used to study critical events in human lung development and suggests a supportive role for the endothelium in promoting branching of airway epithelium.

Predicted loss and gain of function mutations in ACO1 are associated with erythropoiesis.

Hemoglobin is the essential oxygen-carrying molecule in humans and is regulated by cellular iron and oxygen sensing mechanisms. To search for novel variants associated with hemoglobin concentration, we performed genome-wide association studies of hemoglobin concentration using a combined set of 684,122 individuals from Iceland and the UK. Notably, we found seven novel variants, six rare coding and one common, at the ACO1 locus associating with either decreased or increased hemoglobin concentration. Of these variants, the missense Cys506Ser and the stop-gained Lys334Ter mutations are specific to eight and ten generation pedigrees, respectively, and have the two largest effects in the study (EffectCys506Ser = -1.61 SD, CI95 = [-1.98, -1.35]; EffectLys334Ter = 0.63 SD, CI95 = [0.36, 0.91]). We also find Cys506Ser to associate with increased risk of persistent anemia (OR = 17.1, P = 2 × 10-14). The strong bidirectional effects seen in this study implicate ACO1, a known iron sensing molecule, as a major homeostatic regulator of hemoglobin concentration.

A truncating mutation in EPOR leads to hypo-responsiveness to erythropoietin with normal haemoglobin.

The cytokine erythropoietin (EPO), signalling through the EPO receptor (EPO-R), is essential for the formation of red blood cells. We performed a genome-wide association study (GWAS) testing 32.5 million sequence variants for association with serum EPO levels in a set of 4187 individuals. We detect an association between a rare and well imputed stop-gained variant rs370865377[A] (p.Gln82Ter) in EPOR, carried by 1 in 550 Icelanders, and increased serum EPO levels (MAF = 0.09%, Effect = 1.47 SD, P = 3.3 × 10-7). We validated these findings by measuring serum EPO levels in 34 additional pairs of carriers and matched controls and found carriers to have 3.23-fold higher EPO levels than controls (P = 1.7 × 10-6; P combined = 1.6 × 10-11). In contrast to previously reported EPOR mutations, p.Gln82Ter does not associate with haemoglobin levels (Effect = -0.045 SD, P = 0.32, N = 273,160), probably due to a compensatory EPO upregulation in response to EPO-R hypo-responsiveness.

Inhibition of PTP1B disrupts cell-cell adhesion and induces anoikis in breast epithelial cells.

Protein tyrosine phosphatase 1B (PTP1B) is a well-known inhibitor of insulin signaling pathways and inhibitors against PTP1B are being developed as promising drug candidates for treatment of obesity. PTP1B has also been linked to breast cancer both as a tumor suppressor and as an oncogene. Furthermore, PTP1B has been shown to be a regulator of cell adhesion and migration in normal and cancer cells. In this study, we analyzed the PTP1B expression in normal breast tissue, primary breast cells and the breast epithelial cell line D492. In normal breast tissue and primary breast cells, PTP1B is widely expressed in both epithelial and stromal cells, with highest expression in myoepithelial cells and fibroblasts. PTP1B is widely expressed in branching structures generated by D492 when cultured in 3D reconstituted basement membrane (3D rBM). Inhibition of PTP1B in D492 and another mammary epithelial cell line HMLE resulted in reduced cell proliferation and induction of anoikis. These changes were seen when cells were cultured both in monolayer and in 3D rBM. PTP1B inhibition affected cell attachment, expression of cell adhesion proteins and actin polymerization. Moreover, epithelial to mesenchymal transition (EMT) sensitized cells to PTP1B inhibition. A mesenchymal sublines of D492 and HMLE (D492M and HMLEmes) were more sensitive to PTP1B inhibition than D492 and HMLE. Reversion of D492M to an epithelial state using miR-200c-141 restored resistance to detachment induced by PTP1B inhibition. In conclusion, we have shown that PTP1B is widely expressed in the human breast gland with highest expression in myoepithelial cells and fibroblasts. Inhibition of PTP1B in D492 and HMLE affects cell-cell adhesion and induces anoikis-like effects. Finally, cells with an EMT phenotype are more sensitive to PTP1B inhibitors making PTP1B a potential candidate for further studies as a target for drug development in cancer involving the EMT phenotype.