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.

Histopathology and levels of proteins in plasma associate with survival after colorectal cancer diagnosis.

BACKGROUND: The TNM system is used to assess prognosis after colorectal cancer (CRC) diagnosis. Other prognostic factors reported include histopathological assessments of the tumour, tumour mutations and proteins in the blood. As some of these factors are strongly correlated, it is important to evaluate the independent effects they may have on survival. METHODS: Tumour samples from 2162 CRC patients were visually assessed for amount of tumour stroma, severity of lymphocytic infiltrate at the tumour margins and the presence of lymphoid follicles. Somatic mutations in the tumour were assessed for 2134 individuals. Pre-surgical levels of 4963 plasma proteins were measured in 128 individuals. The associations between these features and prognosis were inspected by a Cox Proportional Hazards Model (CPH). RESULTS: Levels of stroma, lymphocytic infiltration and presence of lymphoid follicles all associate with prognosis, along with high tumour mutation burden, high microsatellite instability and TP53 and BRAF mutations. The somatic mutations are correlated with the histopathology and none of the somatic mutations associate with survival in a multivariate analysis. Amount of stroma and lymphocytic infiltration associate with local invasion of tumours. Elevated levels of two plasma proteins, CA-125 and PPP1R1A, associate with a worse prognosis. CONCLUSIONS: Tumour stroma and lymphocytic infiltration variables are strongly associated with prognosis of CRC and capture the prognostic effects of tumour mutation status. CA-125 and PPP1R1A may be useful prognostic biomarkers in CRC.

Physical and cognitive impact following SARS-CoV-2 infection in a large population-based case-control study.

BACKGROUND: Persistent symptoms are common after SARS-CoV-2 infection but correlation with objective measures is unclear. METHODS: We invited all 3098 adults who tested SARS-CoV-2 positive in Iceland before October 2020 to the deCODE Health Study. We compared multiple symptoms and physical measures between 1706 Icelanders with confirmed prior infection (cases) who participated, and 619 contemporary and 13,779 historical controls. Cases participated in the study 5-18 months after infection. RESULTS: Here we report that 41 of 88 symptoms are associated with prior infection, most significantly disturbed smell and taste, memory disturbance, and dyspnea. Measured objectively, cases had poorer smell and taste results, less grip strength, and poorer memory recall. Differences in grip strength and memory recall were small. No other objective measure associated with prior infection including heart rate, blood pressure, postural orthostatic tachycardia, oxygen saturation, exercise tolerance, hearing, and traditional inflammatory, cardiac, liver, and kidney blood biomarkers. There was no evidence of more anxiety or depression among cases. We estimate the prevalence of long Covid to be 7% at a median of 8 months after infection. CONCLUSIONS: We confirm that diverse symptoms are common months after SARS-CoV-2 infection but find few differences between cases and controls in objective parameters measured. These discrepancies between symptoms and physical measures suggest a more complicated contribution to symptoms related to prior infection than is captured with conventional tests. Traditional clinical assessment is not expected to be particularly informative in relating symptoms to a past SARS-CoV-2 infection.

Persistent symptoms are commonly reported after SARS-CoV-2 infection, and this is often described as long Covid. We compared different symptoms reported following SARS-CoV- 2 infection with the results obtained during various medical evaluations that are often used to assess health, such as blood tests, smell tests, taste tests, hearing tests, etc. We compared symptoms and test results between 1,706 Icelanders who had been infected previously with SARS-CoV-2 infection (cases) and 14,398 individuals who had not been infected (controls). Out of 88 assessed symptoms, 41 were more common in cases than controls. However, relatively few differences were seen in the results obtained from the various medical evaluations (cases had poorer smell and taste test results, slightly less grip strength, and slightly poorer memory recall than controls). The differences seen between symptoms and results of medical evaluations suggests that conventional clinical tests may not be informative in relating symptoms to a past SARS-CoV-2 infection.

HLA alleles, disease severity, and age associate with T-cell responses following infection with SARS-CoV-2.

Memory T-cell responses following SARS-CoV-2 infection have been extensively investigated but many studies have been small with a limited range of disease severity. Here we analyze SARS-CoV-2 reactive T-cell responses in 768 convalescent SARS-CoV-2-infected (cases) and 500 uninfected (controls) Icelanders. The T-cell responses are stable three to eight months after SARS-CoV-2 infection, irrespective of disease severity and even those with the mildest symptoms induce broad and persistent T-cell responses. Robust CD4+ T-cell responses are detected against all measured proteins (M, N, S and S1) while the N protein induces strongest CD8+ T-cell responses. CD4+ T-cell responses correlate with disease severity, humoral responses and age, whereas CD8+ T-cell responses correlate with age and functional antibodies. Further, CD8+ T-cell responses associate with several class I HLA alleles. Our results, provide new insight into HLA restriction of CD8+ T-cell immunity and other factors contributing to heterogeneity of T-cell responses following SARS-CoV-2 infection.

Effect of booster vaccination against Delta and Omicron SARS-CoV-2 variants in Iceland.

By the end of July 2021, the majority of the Icelandic population had received vaccination against COVID-19. In mid-July a wave of SARS-CoV-2 infections, dominated by the Delta variant, spread through the population, followed by an Omicron wave in December. A booster vaccination campaign was initiated to curb the spread of the virus. We estimate the risk of infection for different vaccine combinations using vaccination data from 276,028 persons and 963,557 qPCR tests for 277,687 persons. We measure anti-Spike-RBD antibody levels and ACE2-Spike binding inhibitory activity in 371 persons who received one of four recommended vaccination schedules with or without an mRNA vaccine booster. Overall, we find different antibody levels and inhibitory activity in recommended vaccination schedules, reflected in the observed risk of SARS-CoV-2 infections. We observe an increased protection following mRNA boosters, against both Omicron and Delta variant infections, although BNT162b2 boosters provide greater protection against Omicron than mRNA-1273 boosters.

Functional dissection of inherited non-coding variation influencing multiple myeloma risk.

Thousands of non-coding variants have been associated with increased risk of human diseases, yet the causal variants and their mechanisms-of-action remain obscure. In an integrative study combining massively parallel reporter assays (MPRA), expression analyses (eQTL, meQTL, PCHiC) and chromatin accessibility analyses in primary cells (caQTL), we investigate 1,039 variants associated with multiple myeloma (MM). We demonstrate that MM susceptibility is mediated by gene-regulatory changes in plasma cells and B-cells, and identify putative causal variants at six risk loci (SMARCD3, WAC, ELL2, CDCA7L, CEP120, and PREX1). Notably, three of these variants co-localize with significant plasma cell caQTLs, signaling the presence of causal activity at these precise genomic positions in an endogenous chromosomal context in vivo. Our results provide a systematic functional dissection of risk loci for a hematologic malignancy.

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.

The CRTAC1 Protein in Plasma Is Associated With Osteoarthritis and Predicts Progression to Joint Replacement: A Large-Scale Proteomics Scan in Iceland.

OBJECTIVE: Biomarkers for diagnosis and progression of osteoarthritis (OA) are lacking. This study was undertaken to identify circulating biomarkers for OA that could predict disease occurrence and/or progression to joint replacement. METHODS: Using the SomaScan platform, we measured 4,792 proteins in plasma from 37,278 individuals, of whom 12,178 individuals had OA and 2,524 had undergone joint replacement. We performed a case-control study for identification of potential protein biomarkers for hip, knee, and/or hand OA, and a prospective study for identification of biomarkers for joint replacement. RESULTS: Among the large panel of plasma proteins assessed, cartilage acidic protein 1 (CRTAC1) was the most strongly associated with both OA diagnosis (odds ratio 1.46 [95% confidence interval 1.41-1.52] for knee OA, odds ratio 1.36 [95% confidence interval 1.29-1.43] for hip OA, and odds ratio 1.33 [95% confidence interval 1.26-1.40] for hand OA) and progression to joint replacement (hazard ratio 1.40 [95% confidence interval 1.30-1.51] for knee replacement and hazard ratio 1.31 [95% confidence interval 1.19-1.45] for hip replacement). Patients with OA who were in the highest quintile of risk of joint replacement, based on known risk factors (i.e., age, sex, and body mass index) and plasma CRTAC1 level, were 16 times more likely to undergo knee replacement within 5 years of plasma sample collection than those in the lowest quintile, and 6.5 times more likely to undergo hip replacement. CRTAC1 was not associated with other types of inflammatory arthritis. A specific protein profile was identified in those patients who had undergone joint replacement prior to plasma sample collection. CONCLUSION: Through a hypothesis-free approach, we identified CRTAC1 in plasma as a novel promising candidate biomarker for OA that is both associated with occurrence of OA and predictive of progression to joint replacement. This biomarker might also be useful in the selection of suitable patients for clinical trial enrollment.

Lifelong Reduction in LDL (Low-Density Lipoprotein) Cholesterol due to a Gain-of-Function Mutation in LDLR.

BACKGROUND: Loss-of-function mutations in the LDL (low-density lipoprotein) receptor gene (LDLR) cause elevated levels of LDL cholesterol and premature cardiovascular disease. To date, a gain-of-function mutation in LDLR with a large effect on LDL cholesterol levels has not been described. Here, we searched for sequence variants in LDLR that have a large effect on LDL cholesterol levels. METHODS: We analyzed whole-genome sequencing data from 43 202 Icelanders. Single-nucleotide polymorphisms and structural variants including deletions, insertions, and duplications were genotyped using whole-genome sequencing-based data. LDL cholesterol associations were carried out in a sample of >100 000 Icelanders with genetic information (imputed or whole-genome sequencing). Molecular analyses were performed using RNA sequencing and protein expression assays in Epstein-Barr virus-transformed lymphocytes. RESULTS: We discovered a 2.5-kb deletion (del2.5) overlapping the 3' untranslated region of LDLR in 7 heterozygous carriers from a single family. Mean level of LDL cholesterol was 74% lower in del2.5 carriers than in 101 851 noncarriers, a difference of 2.48 mmol/L (96 mg/dL; P=8.4×10-8). Del2.5 results in production of an alternative mRNA isoform with a truncated 3' untranslated region. The truncation leads to a loss of target sites for microRNAs known to repress translation of LDLR. In Epstein-Barr virus-transformed lymphocytes derived from del2.5 carriers, expression of alternative mRNA isoform was 1.84-fold higher than the wild-type isoform (P=0.0013), and there was 1.79-fold higher surface expression of the LDL receptor than in noncarriers (P=0.0086). We did not find a highly penetrant detrimental impact of lifelong very low levels of LDL cholesterol due to del2.5 on health of the carriers. CONCLUSIONS: Del2.5 is the first reported gain-of-function mutation in LDLR causing a large reduction in LDL cholesterol. These data point to a role for alternative polyadenylation of LDLR mRNA as a potent regulator of LDL receptor expression in humans.

Humoral Immune Response to SARS-CoV-2 in Iceland.

BACKGROUND: Little is known about the nature and durability of the humoral immune response to infection with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). METHODS: We measured antibodies in serum samples from 30,576 persons in Iceland, using six assays (including two pan-immunoglobulin [pan-Ig] assays), and we determined that the appropriate measure of seropositivity was a positive result with both pan-Ig assays. We tested 2102 samples collected from 1237 persons up to 4 months after diagnosis by a quantitative polymerase-chain-reaction (qPCR) assay. We measured antibodies in 4222 quarantined persons who had been exposed to SARS-CoV-2 and in 23,452 persons not known to have been exposed. RESULTS: Of the 1797 persons who had recovered from SARS-CoV-2 infection, 1107 of the 1215 who were tested (91.1%) were seropositive; antiviral antibody titers assayed by two pan-Ig assays increased during 2 months after diagnosis by qPCR and remained on a plateau for the remainder of the study. Of quarantined persons, 2.3% were seropositive; of those with unknown exposure, 0.3% were positive. We estimate that 0.9% of Icelanders were infected with SARS-CoV-2 and that the infection was fatal in 0.3%. We also estimate that 56% of all SARS-CoV-2 infections in Iceland had been diagnosed with qPCR, 14% had occurred in quarantined persons who had not been tested with qPCR (or who had not received a positive result, if tested), and 30% had occurred in persons outside quarantine and not tested with qPCR. CONCLUSIONS: Our results indicate that antiviral antibodies against SARS-CoV-2 did not decline within 4 months after diagnosis. We estimate that the risk of death from infection was 0.3% and that 44% of persons infected with SARS-CoV-2 in Iceland were not diagnosed by qPCR.

Human Embryonic-Derived Mesenchymal Progenitor Cells (hES-MP Cells) are Fully Supported in Culture with Human Platelet Lysates.

Human embryonic stem cell-derived mesenchymal progenitor (hES-MP) cells are mesenchymal-like cells, derived from human embryonic stem cells without the aid of feeder cells. They have been suggested as a potential alternative to mesenchymal stromal cells (MSCs) in regenerative medicine due to their mesenchymal-like proliferation and differentiation characteristics. Cells and cell products intended for regenerative medicine in humans should be derived, expanded and differentiated using conditions free of animal-derived products to minimize risk of animal-transmitted disease and immune reactions to foreign proteins. Human platelets are rich in growth factors needed for cell culture and have been used successfully as an animal serum replacement for MSC expansion and differentiation. In this study, we compared the proliferation of hES-MP cells and MSCs; the hES-MP cell growth was sustained for longer than that of MSCs. Growth factors, gene expression, and surface marker expression in hES-MP cells cultured with either human platelet lysate (hPL) or fetal bovine serum (FBS) supplementation were compared, along with differentiation to osteogenic and chondrogenic lineages. Despite some differences between hES-MP cells grown in hPL- and FBS-supplemented media, hPL was found to be a suitable replacement for FBS. In this paper, we demonstrate for the first time that hES-MP cells can be grown using platelet lysates from expired platelet concentrates (hPL).

FLT3 stop mutation increases FLT3 ligand level and risk of autoimmune thyroid disease.

Autoimmune thyroid disease is the most common autoimmune disease and is highly heritable1. Here, by using a genome-wide association study of 30,234 cases and 725,172 controls from Iceland and the UK Biobank, we find 99 sequence variants at 93 loci, of which 84 variants are previously unreported2-7. A low-frequency (1.36%) intronic variant in FLT3 (rs76428106-C) has the largest effect on risk of autoimmune thyroid disease (odds ratio (OR) = 1.46, P = 2.37 × 10-24). rs76428106-C is also associated with systemic lupus erythematosus (OR = 1.90, P = 6.46 × 10-4), rheumatoid factor and/or anti-CCP-positive rheumatoid arthritis (OR = 1.41, P = 4.31 × 10-4) and coeliac disease (OR = 1.62, P = 1.20 × 10-4). FLT3 encodes fms-related tyrosine kinase 3, a receptor that regulates haematopoietic progenitor and dendritic cells. RNA sequencing revealed that rs76428106-C generates a cryptic splice site, which introduces a stop codon in 30% of transcripts that are predicted to encode a truncated protein, which lacks its tyrosine kinase domains. Each copy of rs76428106-C doubles the plasma levels of the FTL3 ligand. Activating somatic mutations in FLT3 are associated with acute myeloid leukaemia8 with a poor prognosis and rs76428106-C also predisposes individuals to acute myeloid leukaemia (OR = 1.90, P = 5.40 × 10-3). Thus, a predicted loss-of-function germline mutation in FLT3 causes a reduction in full-length FLT3, with a compensatory increase in the levels of its ligand and an increased disease risk, similar to that of a gain-of-function mutation.

Publisher Correction: GWAS of bone size yields twelve loci that also affect height, BMD, osteoarthritis or fractures.

The original HTML version of this Article was updated shortly after publication to add links to the Peer Review file.In addition, affiliations 16 and 17 incorrectly read 'School of Medicine Sydney, University of Notre Dame Australia, Sydney, WA, 6160, Australia' and 'St Vincent's Clinical School, University of New South Wales Medicine, University of New South Wales, Sydney, NSW, 2052, Australia.' This has now been corrected in both the PDF and HTML versions of the Article.

GWAS of bone size yields twelve loci that also affect height, BMD, osteoarthritis or fractures.

Bone area is one measure of bone size that is easily derived from dual-energy X-ray absorptiometry (DXA) scans. In a GWA study of DXA bone area of the hip and lumbar spine (N ≥ 28,954), we find thirteen independent association signals at twelve loci that replicate in samples of European and East Asian descent (N = 13,608 - 21,277). Eight DXA area loci associate with osteoarthritis, including rs143384 in GDF5 and a missense variant in COL11A1 (rs3753841). The strongest DXA area association is with rs11614913[T] in the microRNA MIR196A2 gene that associates with lumbar spine area (P = 2.3 × 10-42, ß = -0.090) and confers risk of hip fracture (P = 1.0 × 10-8, OR = 1.11). We demonstrate that the risk allele is less efficient in repressing miR-196a-5p target genes. We also show that the DXA area measure contributes to the risk of hip fracture independent of bone density.