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.

Rare variants with large effects provide functional insights into the pathology of migraine subtypes, with and without aura

Migraine is a complex neurovascular disease with a range of severity and symptoms, yet mostly studied as one phenotype in genome-wide association studies (GWAS). Here we combine large GWAS datasets from six European populations to study the main migraine subtypes, migraine with aura (MA) and migraine without aura (MO). We identified four new MA-associated variants (in PRRT2, PALMD, ABO and LRRK2) and classified 13 MO-associated variants. Rare variants with large effects highlight three genes. A rare frameshift variant in brain-expressed PRRT2 confers large risk of MA and epilepsy, but not MO. A burden test of rare loss-of-function variants in SCN11A, encoding a neuron-expressed sodium channel with a key role in pain sensation, shows strong protection against migraine. Finally, a rare variant with cis-regulatory effects on KCNK5 confers large protection against migraine and brain aneurysms. Our findings offer new insights with therapeutic potential into the complex biology of migraine and its subtypes.

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.

Sequence variant affects GCSAML splicing, mast cell specific proteins, and risk of urticaria

Urticaria is a skin disorder characterized by outbreaks of raised pruritic wheals. In order to identify sequence variants associated with urticaria, we performed a meta-analysis of genome-wide association studies for urticaria with a total of 40,694 cases and 1,230,001 controls from Iceland, the UK, Finland, and Japan. We also performed transcriptome- and proteome-wide analyses in Iceland and the UK. We found nine sequence variants at nine loci associating with urticaria. The variants are at genes participating in type 2 immune responses and/or mast cell biology (CBLB, FCER1A, GCSAML, STAT6, TPSD1, ZFPM1), the innate immunity (C4), and NF-κB signaling. The most significant association was observed for the splice-donor variant rs56043070[A] (hg38: chr1:247556467) in GCSAML (MAF = 6.6%, OR = 1.24 (95%CI: 1.20-1.28), P-value = 3.6 × 10^-44). We assessed the effects of the variants on transcripts, and levels of proteins relevant to urticaria pathophysiology. Our results emphasize the role of type 2 immune response and mast cell activation in the pathogenesis of urticaria. Our findings may point to an IgE-independent urticaria pathway that could help address unmet clinical need.

Deficit of homozygosity among 1.52 million individuals and genetic causes of recessive lethality

Genotypes causing pregnancy loss and perinatal mortality are depleted among living individuals and are therefore difficult to find. To explore genetic causes of recessive lethality, we searched for sequence variants with deficit of homozygosity among 1.52 million individuals from six European populations. In this study, we identified 25 genes harboring protein-altering sequence variants with a strong deficit of homozygosity (10% or less of predicted homozygotes). Sequence variants in 12 of the genes cause Mendelian disease under a recessive mode of inheritance, two under a dominant mode, but variants in the remaining 11 have not been reported to cause disease. Sequence variants with a strong deficit of homozygosity are over-represented among genes essential for growth of human cell lines and genes orthologous to mouse genes known to affect viability. The function of these genes gives insight into the genetics of intrauterine lethality. We also identified 1077 genes with homozygous predicted loss-of-function genotypes not previously described, bringing the total set of genes completely knocked out in humans to 4785.

Multiomics study of nonalcoholic fatty liver disease

Nonalcoholic fatty liver (NAFL) and its sequelae are growing health problems. We performed a genome-wide association study of NAFL, cirrhosis and hepatocellular carcinoma, and integrated the findings with expression and proteomic data. For NAFL, we utilized 9,491 clinical cases and proton density fat fraction extracted from 36,116 liver magnetic resonance images. We identified 18 sequence variants associated with NAFL and 4 with cirrhosis, and found rare, protective, predicted loss-of-function variants in MTARC1 and GPAM, underscoring them as potential drug targets. We leveraged messenger RNA expression, splicing and predicted coding effects to identify 16 putative causal genes, of which many are implicated in lipid metabolism. We analyzed levels of 4,907 plasma proteins in 35,559 Icelanders and 1,459 proteins in 47,151 UK Biobank participants, identifying multiple proteins involved in disease pathogenesis. We show that proteomics can discriminate between NAFL and cirrhosis. The present study provides insights into the development of noninvasive evaluation of NAFL and new therapeutic options.

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.

Iceland has used four different SARS-CoV-2 vaccines in various combinations. Here, the authors describe differences in the immune responses elicited by different initial/booster vaccine combinations, and then use population-level data to assess the effects of booster doses against Delta and Omicron 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.

A study of 768 convalescent SARS CoV-2-infected and 500 uninfected Icelanders reveals broad and stable T-cell responses 3-8 months from infection. HLA alleles, disease severity, and age contribute to the heterogeneity of cellular immunity.

Multiomics analysis of rheumatoid arthritis yields sequence variants that have large effects on risk of the seropositive subset

OBJECTIVES

To find causal genes for rheumatoid arthritis (RA) and its seropositive (RF and/or ACPA positive) and seronegative subsets.

METHODS

We performed a genome-wide association study (GWAS) of 31 313 RA cases (68% seropositive) and ~1 million controls from Northwestern Europe. We searched for causal genes outside the HLA-locus through effect on coding, mRNA expression in several tissues and/or levels of plasma proteins (SomaScan) and did network analysis (Qiagen).

RESULTS

We found 25 sequence variants for RA overall, 33 for seropositive and 2 for seronegative RA, altogether 37 sequence variants at 34 non-HLA loci, of which 15 are novel. Genomic, transcriptomic and proteomic analysis of these yielded 25 causal genes in seropositive RA and additional two overall. Most encode proteins in the network of interferon-alpha/beta and IL-12/23 that signal through the JAK/STAT-pathway. Highlighting those with largest effect on seropositive RA, a rare missense variant in STAT4 (rs140675301-A) that is independent of reported non-coding STAT4-variants, increases the risk of seropositive RA 2.27-fold (p=2.1×10^-9), more than the rs2476601-A missense variant in PTPN22 (OR=1.59, p=1.3×10^-160). STAT4 rs140675301-A replaces hydrophilic glutamic acid with hydrophobic valine (Glu128Val) in a conserved, surface-exposed loop. A stop-mutation (rs76428106-C) in FLT3 increases seropositive RA risk (OR=1.35, p=6.6×10^-11). Independent missense variants in TYK2 (rs34536443-C, rs12720356-C, rs35018800-A, latter two novel) associate with decreased risk of seropositive RA (ORs=0.63-0.87, p=10^-9-10^-27) and decreased plasma levels of interferon-alpha/beta receptor 1 that signals through TYK2/JAK1/STAT4.

CONCLUSION

Sequence variants pointing to causal genes in the JAK/STAT pathway have largest effect on seropositive RA, while associations with seronegative RA remain scarce.

The sequences of 150,119 genomes in the UK Biobank

Detailed knowledge of how diversity in the sequence of the human genome affects phenotypic diversity depends on a comprehensive and reliable characterization of both sequences and phenotypic variation. Over the past decade, insights into this relationship have been obtained from whole-exome sequencing or whole-genome sequencing of large cohorts with rich phenotypic data^1,2. Here we describe the analysis of whole-genome sequencing of 150,119 individuals from the UK Biobank³. This constitutes a set of high-quality variants, including 585,040,410 single-nucleotide polymorphisms, representing 7.0% of all possible human single-nucleotide polymorphisms, and 58,707,036 indels. This large set of variants allows us to characterize selection based on sequence variation within a population through a depletion rank score of windows along the genome. Depletion rank analysis shows that coding exons represent a small fraction of regions in the genome subject to strong sequence conservation. We define three cohorts within the UK Biobank: a large British Irish cohort, a smaller African cohort and a South Asian cohort. A haplotype reference panel is provided that allows reliable imputation of most variants carried by three or more sequenced individuals. We identified 895,055 structural variants and 2,536,688 microsatellites, groups of variants typically excluded from large-scale whole-genome sequencing studies. Using this formidable new resource, we provide several examples of trait associations for rare variants with large effects not found previously through studies based on whole-exome sequencing and/or imputation.

Reconstruction of a large-scale outbreak of SARS-CoV-2 infection in Iceland informs vaccination strategies

OBJECTIVES

The spread of SARS-CoV-2 is dependent on several factors, both biological and behavioural. The effectiveness of nonpharmaceutical interventions can be attributed largely to changes in human behaviour, but quantifying this effect remains challenging. Reconstructing the transmission tree of the third wave of SARS-CoV-2 infections in Iceland using contact tracing and viral sequence data from 2522 cases enables us to directly compare the infectiousness of distinct groups of persons.

METHODS

The transmission tree enables us to model the effect that a given population prevalence of vaccination would have had on the third wave had one of three different vaccination strategies been implemented before that time. This allows us to compare the effectiveness of the strategies in terms of minimizing the number of cases, deaths, critical cases, and severe cases.

RESULTS

We found that people diagnosed outside of quarantine (Rˆ=1.31) were 89% more infectious than those diagnosed while in quarantine (Rˆ=0.70) and that infectiousness decreased as a function of time spent in quarantine before diagnosis, with people diagnosed outside of quarantine being 144% more infectious than those diagnosed after ≥3 days in quarantine (Rˆ=0.54). People of working age, 16 to 66 years (Rˆ=1.08), were 46% more infectious than those outside of that age range (Rˆ=0.74).

DISCUSSION

We found that vaccinating the population in order of ascending age or uniformly at random would have prevented more infections per vaccination than vaccinating in order of descending age, without significantly affecting the expected number of deaths, critical cases, or severe cases.

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.

Evgenia Mikaelsdottir et al. report a study of variants associated with platelet count among European individuals where they identify 577 associations. They also report a genetic overlap between platelet count and human diseases, including myeloproliferative neoplasms, rheumatoid arthritis, and hypertension, as well as a genetic overlap between platelet count and various physiological markers.

The genetic architecture of age-related hearing impairment revealed by genome-wide association analysis

Age-related hearing impairment (ARHI) is the most common sensory disorder in older adults. We conducted a genome-wide association meta-analysis of 121,934 ARHI cases and 591,699 controls from Iceland and the UK. We identified 21 novel sequence variants, of which 13 are rare, under either additive or recessive models. Of special interest are a missense variant in LOXHD1 (MAF = 1.96%) and a tandem duplication in FBF1 covering 4 exons (MAF = 0.22%) associating with ARHI (OR = 3.7 for homozygotes, P = 1.7 × 10-22 and OR = 4.2 for heterozygotes, P = 5.7 × 10-27, respectively). We constructed an ARHI genetic risk score (GRS) using common variants and showed that a common variant GRS can identify individuals at risk comparable to carriers of rare high penetrance variants. Furthermore, we found that ARHI and tinnitus share genetic causes. This study sheds a new light on the genetic architecture of ARHI, through several rare variants in both Mendelian deafness genes and genes not previously linked to hearing.

Loss-of-Function Variants in the Tumor-Suppressor Gene PTPN14 Confer Increased Cancer Risk

The success of genome-wide association studies (GWAS) in identifying common, low-penetrance variant-cancer associations for the past decade is undisputed. However, discovering additional high-penetrance cancer mutations in unknown cancer predisposing genes requires detection of variant-cancer association of ultra-rare coding variants. Consequently, large-scale next-generation sequence data with associated phenotype information are needed. Here, we used genotype data on 166,281 Icelanders, of which, 49,708 were whole-genome sequenced and 408,595 individuals from the UK Biobank, of which, 41,147 were whole-exome sequenced, to test for association between loss-of-function burden in autosomal genes and basal cell carcinoma (BCC), the most common cancer in Caucasians. A total of 25,205 BCC cases and 683,058 controls were tested. Rare germline loss-of-function variants in PTPN14 conferred substantial risks of BCC (OR, 8.0; P = 1.9 × 10^-12), with a quarter of carriers getting BCC before age 70 and over half in their lifetime. Furthermore, common variants at the PTPN14 locus were associated with BCC, suggesting PTPN14 as a new, high-impact BCC predisposition gene. A follow-up investigation of 24 cancers and three benign tumor types showed that PTPN14 loss-of-function variants are associated with high risk of cervical cancer (OR, 12.7, P = 1.6 × 10^-4) and low age at diagnosis. Our findings, using power-increasing methods with high-quality rare variant genotypes, highlight future prospects for new discoveries on carcinogenesis. SIGNIFICANCE: This study identifies the tumor-suppressor gene PTPN14 as a high-impact BCC predisposition gene and indicates that inactivation of PTPN14 by germline sequence variants may also lead to increased risk of cervical cancer.

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.

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.

Spread of SARS-CoV-2 in the Icelandic Population

BACKGROUND

During the current worldwide pandemic, coronavirus disease 2019 (Covid-19) was first diagnosed in Iceland at the end of February. However, data are limited on how SARS-CoV-2, the virus that causes Covid-19, enters and spreads in a population.

METHODS

We targeted testing to persons living in Iceland who were at high risk for infection (mainly those who were symptomatic, had recently traveled to high-risk countries, or had contact with infected persons). We also carried out population screening using two strategies: issuing an open invitation to 10,797 persons and sending random invitations to 2283 persons. We sequenced SARS-CoV-2 from 643 samples.

RESULTS

As of April 4, a total of 1221 of 9199 persons (13.3%) who were recruited for targeted testing had positive results for infection with SARS-CoV-2. Of those tested in the general population, 87 (0.8%) in the open-invitation screening and 13 (0.6%) in the random-population screening tested positive for the virus. In total, 6% of the population was screened. Most persons in the targeted-testing group who received positive tests early in the study had recently traveled internationally, in contrast to those who tested positive later in the study. Children under 10 years of age were less likely to receive a positive result than were persons 10 years of age or older, with percentages of 6.7% and 13.7%, respectively, for targeted testing; in the population screening, no child under 10 years of age had a positive result, as compared with 0.8% of those 10 years of age or older. Fewer females than males received positive results both in targeted testing (11.0% vs. 16.7%) and in population screening (0.6% vs. 0.9%). The haplotypes of the sequenced SARS-CoV-2 viruses were diverse and changed over time. The percentage of infected participants that was determined through population screening remained stable for the 20-day duration of screening.

CONCLUSIONS

In a population-based study in Iceland, children under 10 years of age and females had a lower incidence of SARS-CoV-2 infection than adolescents or adults and males. The proportion of infected persons identified through population screening did not change substantially during the screening period, which was consistent with a beneficial effect of containment efforts. (Funded by deCODE Genetics-Amgen.).

Genome-wide association identifies seven loci for pelvic organ prolapse in Iceland and the UK Biobank

Pelvic organ prolapse (POP) is a downward descent of one or more of the pelvic organs, resulting in a protrusion of the vaginal wall and/or uterus. We performed a genome-wide association study of POP using data from Iceland and the UK Biobank, a total of 15,010 cases with hospital-based diagnosis code and 340,734 female controls, and found eight sequence variants at seven loci associating with POP (P < 5 × 10-8); seven common (minor allele frequency >5%) and one with minor allele frequency of 4.87%. Some of the variants associating with POP also associated with traits of similar pathophysiology. Of these, rs3820282, which may alter the estrogen-based regulation of WNT4, also associates with leiomyoma of uterus, gestational duration and endometriosis. Rs3791675 at EFEMP1, a gene involved in connective tissue homeostasis, also associates with hernias and carpal tunnel syndrome. Our results highlight the role of connective tissue metabolism and estrogen exposure in the etiology of POP.

GraphTyper2 enables population-scale genotyping of structural variation using pangenome graphs

Analysis of sequence diversity in the human genome is fundamental for genetic studies. Structural variants (SVs) are frequently omitted in sequence analysis studies, although each has a relatively large impact on the genome. Here, we present GraphTyper2, which uses pangenome graphs to genotype SVs and small variants using short-reads. Comparison to the syndip benchmark dataset shows that our SV genotyping is sensitive and variant segregation in families demonstrates the accuracy of our approach. We demonstrate that incorporating public assembly data into our pipeline greatly improves sensitivity, particularly for large insertions. We validate 6,812 SVs on average per genome using long-read data of 41 Icelanders. We show that GraphTyper2 can simultaneously genotype tens of thousands of whole-genomes by characterizing 60 million small variants and half a million SVs in 49,962 Icelanders, including 80 thousand SVs with high-confidence.

A PRPH splice-donor variant associates with reduced sural nerve amplitude and risk of peripheral neuropathy

Nerve conduction (NC) studies generate measures of peripheral nerve function that can reveal underlying pathology due to axonal loss, demyelination or both. We perform a genome-wide association study of sural NC amplitude and velocity in 7045 Icelanders and find a low-frequency splice-donor variant in PRPH (c.996+1G>A; MAF = 1.32%) associating with decreased NC amplitude but not velocity. PRPH encodes peripherin, an intermediate filament (IF) protein involved in cytoskeletal development and maintenance of neurons. Through RNA and protein studies, we show that the variant leads to loss-of-function (LoF), as when over-expressed in a cell line devoid of other IFs, it does not allow formation of the normal filamentous structure of peripherin, yielding instead punctate protein inclusions. Recall of carriers for neurological assessment confirms that from an early age, homozygotes have significantly lower sural NC amplitude than non-carriers and are at risk of a mild, early-onset, sensory-negative, axonal polyneuropathy.

Diagnosis and classification of peripheral neuropathy (PN) is facilitated by nerve conduction (NC) studies. Here, Bjornsdottir et al. find a low-frequency PRPH splice-donor variant that associates with NC amplitude and neurological assessment of recalled PRPH variant carriers reveals increased risk of a mild sensory-negative PN.

Coding variants in RPL3L and MYZAP increase risk of atrial fibrillation

Most sequence variants identified hitherto in genome-wide association studies (GWAS) of atrial fibrillation are common, non-coding variants associated with risk through unknown mechanisms. We performed a meta-analysis of GWAS of atrial fibrillation among 29,502 cases and 767,760 controls from Iceland and the UK Biobank with follow-up in samples from Norway and the US, focusing on low-frequency coding and splice variants aiming to identify causal genes. We observe associations with one missense (OR = 1.20) and one splice-donor variant (OR = 1.50) in RPL3L, the first ribosomal gene implicated in atrial fibrillation to our knowledge. Analysis of 167 RNA samples from the right atrium reveals that the splice-donor variant in RPL3L results in exon skipping. We also observe an association with a missense variant in MYZAP (OR = 1.38), encoding a component of the intercalated discs of cardiomyocytes. Both discoveries emphasize the close relationship between the mechanical and electrical function of the heart.

Pseudoalignment for metagenomic read assignment

Motivation

Read assignment is an important first step in many metagenomic analysis workflows, providing the basis for identification and quantification of species. However ambiguity among the sequences of many strains makes it difficult to assign reads at the lowest level of taxonomy, and reads are typically assigned to taxonomic levels where they are unambiguous. We explore connections between metagenomic read assignment and the quantification of transcripts from RNA-Seq data in order to develop novel methods for rapid and accurate quantification of metagenomic strains.

Results

We find that the recent idea of pseudoalignment introduced in the RNA-Seq context is highly applicable in the metagenomics setting. When coupled with the Expectation-Maximization (EM) algorithm, reads can be assigned far more accurately and quickly than is currently possible with state of the art software, making it possible and practical for the first time to analyze abundances of individual genomes in metagenomics projects.

Availability and Implementation

Pipeline and analysis code can be downloaded from http://github.com/pachterlab/metakallisto.

Contact

lpachter@math.berkeley.edu.