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.

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.).

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.

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.

Genome-wide meta-analysis identifies eight new susceptibility loci for cutaneous squamous cell carcinoma.

Cutaneous squamous cell carcinoma (SCC) is one of the most common cancers in the United States. Previous genome-wide association studies (GWAS) have identified 14 single nucleotide polymorphisms (SNPs) associated with cutaneous SCC. Here, we report the largest cutaneous SCC meta-analysis to date, representing six international cohorts and totaling 19,149 SCC cases and 680,049 controls. We discover eight novel loci associated with SCC, confirm all previously associated loci, and perform fine mapping of causal variants. The novel SNPs occur within skin-specific regulatory elements and implicate loci involved in cancer development, immune regulation, and keratinocyte differentiation in SCC susceptibility.

Common and Rare Sequence Variants Influencing Tumor Biomarkers in Blood.

BACKGROUND: Alpha-fetoprotein (AFP), cancer antigens 15.3, 19.9, and 125, carcinoembryonic antigen, and alkaline phosphatase (ALP) are widely measured in attempts to detect cancer and to monitor treatment response. However, due to lack of sensitivity and specificity, their utility is debated. The serum levels of these markers are affected by a number of nonmalignant factors, including genotype. Thus, it may be possible to improve both sensitivity and specificity by adjusting test results for genetic effects. METHODS: We performed genome-wide association studies of serum levels of AFP (N = 22,686), carcinoembryonic antigen (N = 22,309), cancer antigens 15.3 (N = 7,107), 19.9 (N = 9,945), and 125 (N = 9,824), and ALP (N = 162,774). We also examined the correlations between levels of these biomarkers and the presence of cancer, using data from a nationwide cancer registry. RESULTS: We report a total of 84 associations of 79 sequence variants with levels of the six biomarkers, explaining between 2.3% and 42.3% of the phenotypic variance. Among the 79 variants, 22 are cis (in- or near the gene encoding the biomarker), 18 have minor allele frequency less than 1%, 31 are coding variants, and 7 are associated with gene expression in whole blood. We also find multiple conditions associated with higher biomarker levels. CONCLUSIONS: Our results provide insights into the genetic contribution to diversity in concentration of tumor biomarkers in blood. IMPACT: Genetic correction of biomarker values could improve prediction algorithms and decision-making based on these biomarkers.

Characterizing mutagenic effects of recombination through a sequence-level genetic map.

Genetic diversity arises from recombination and de novo mutation (DNM). Using a combination of microarray genotype and whole-genome sequence data on parent-child pairs, we identified 4,531,535 crossover recombinations and 200,435 DNMs. The resulting genetic map has a resolution of 682 base pairs. Crossovers exhibit a mutagenic effect, with overrepresentation of DNMs within 1 kilobase of crossovers in males and females. In females, a higher mutation rate is observed up to 40 kilobases from crossovers, particularly for complex crossovers, which increase with maternal age. We identified 35 loci associated with the recombination rate or the location of crossovers, demonstrating extensive genetic control of meiotic recombination, and our results highlight genes linked to the formation of the synaptonemal complex as determinants of crossovers.

Multiple transmissions of de novo mutations in families.

De novo mutations (DNMs) cause a large proportion of severe rare diseases of childhood. DNMs that occur early may result in mosaicism of both somatic and germ cells. Such early mutations can cause recurrence of disease. We scanned 1,007 sibling pairs from 251 families and identified 878 DNMs shared by siblings (ssDNMs) at 448 genomic sites. We estimated DNM recurrence probability based on parental mosaicism, sharing of DNMs among siblings, parent-of-origin, mutation type and genomic position. We detected 57.2% of ssDNMs in the parental blood. The recurrence probability of a DNM decreases by 2.27% per year for paternal DNMs and 1.78% per year for maternal DNMs. Maternal ssDNMs are more likely to be T>C mutations than paternal ssDNMs, and less likely to be C>T mutations. Depending on the properties of the DNM, the recurrence probability ranges from 0.011% to 28.5%. We have launched an online calculator to allow estimation of DNM recurrence probability for research purposes.

Genome-wide associations for benign prostatic hyperplasia reveal a genetic correlation with serum levels of PSA.

Benign prostatic hyperplasia and associated lower urinary tract symptoms (BPH/LUTS) are common conditions affecting the majority of elderly males. Here we report the results of a genome-wide association study of symptomatic BPH/LUTS in 20,621 patients and 280,541 controls of European ancestry, from Iceland and the UK. We discovered 23 genome-wide significant variants, located at 14 loci. There is little or no overlap between the BPH/LUTS variants and published prostate cancer risk variants. However, 15 of the variants reported here also associate with serum levels of prostate specific antigen (PSA) (at a Bonferroni corrected P < 0.0022). Furthermore, there is a strong genetic correlation, rg = 0.77 (P = 2.6 × 10-11), between PSA and BPH/LUTS, and one standard deviation increase in a polygenic risk score (PRS) for BPH/LUTS increases PSA levels by 12.9% (P = 1.6×10-55). These results shed a light on the genetic background of BPH/LUTS and its substantial influence on PSA levels.

Insights into imprinting from parent-of-origin phased methylomes and transcriptomes.

Imprinting is the preferential expression of one parental allele over the other. It is controlled primarily through differential methylation of cytosine at CpG dinucleotides. Here we combine 285 methylomes and 11,617 transcriptomes from peripheral blood samples with parent-of-origin phased haplotypes, to produce a new map of imprinted methylation and gene expression patterns across the human genome. We demonstrate how imprinted methylation is a continuous rather than a binary characteristic. We describe at high resolution the parent-of-origin methylation pattern at the 15q11.2 Prader-Willi/Angelman syndrome locus, with nearly confluent stochastic paternal methylation punctuated by 'spikes' of maternal methylation. We find examples of polymorphic imprinted methylation unrelated (at VTRNA2-1 and PARD6G) or related (at CHRNE) to nearby SNP genotypes. We observe RNA isoform-specific imprinted expression patterns suggestive of a methylation-sensitive transcriptional elongation block. Finally, we gain new insights into parent-of-origin-specific effects on phenotypes at the DLK1/MEG3 and GNAS loci.

Variants associating with uterine leiomyoma highlight genetic background shared by various cancers and hormone-related traits.

Uterine leiomyomas are common benign tumors of the myometrium. We performed a meta-analysis of two genome-wide association studies of leiomyoma in European women (16,595 cases and 523,330 controls), uncovering 21 variants at 16 loci that associate with the disease. Five variants were previously reported to confer risk of various malignant or benign tumors (rs78378222 in TP53, rs10069690 in TERT, rs1800057 and rs1801516 in ATM, and rs7907606 at OBFC1) and four signals are located at established risk loci for hormone-related traits (endometriosis and breast cancer) at 1q36.12 (CDC42/WNT4), 2p25.1 (GREB1), 20p12.3 (MCM8), and 6q26.2 (SYNE1/ESR1). Polygenic score for leiomyoma, computed using UKB data, is significantly correlated with risk of cancer in the Icelandic population. Functional annotation suggests that the non-coding risk variants affect multiple genes, including ESR1. Our results provide insights into the genetic background of leiomyoma that are shared by other benign and malignant tumors and highlight the role of hormones in leiomyoma growth.

Association of BRCA2 K3326* With Small Cell Lung Cancer and Squamous Cell Cancer of the Skin.

Background: Most pathogenic mutations in the BRCA2 gene carry a high risk of hereditary breast and ovarian cancer (HBOC). However, a stop-gain mutation, K3326* (rs11571833), confers risk of lung cancer and cancers of the upper-aero-digestive tract but only a modest risk of breast or ovarian cancer. The Icelandic population provides an opportunity for comprehensive characterization of the cancer risk profiles of K3326* and HBOC mutations because a single mutation, BRCA2 999del5, is responsible for almost all BRCA2-related HBOC in the population. Methods: Genotype information on 43 641 cancer patients and 370 971 control subjects from Iceland, the Netherlands, and the United States was used to assess the cancer risk profiles of K3326* and BRCA2 999del5. BRCA2 expression was assessed using RNAseq data from blood (n = 2233), as well as 52 tissues reported in the GTEx database. Results: The cancer risks associated with K3326* are fundamentally different from those associated with 999del5. We report for the first time an association between K3326* and small cell lung cancer (odds ratio [OR] = 2.06, 95% confidence interval [CI] = 1.35 to 3.16) and squamous cell carcinoma of the skin (OR = 1.69, 95% CI = 1.26 to 2.26). Individuals homozygous for K3326* reach old age and have children. Unlike BRCA2 999del5, the K3326* allele does not affect the level of BRCA2 transcripts, and the allele is expressed to the same extent as the wild-type allele. Conclusions: K3326* associates primarily with cancers that have strong environmental genotoxic risk factors. Expression of the K3326* allele suggests that a variant protein may be made that retains the DNA repair capabilities important to hormone-responsive tissues but may be less efficient in responding to genotoxic stress.

Whole genome characterization of sequence diversity of 15,220 Icelanders.

Understanding of sequence diversity is the cornerstone of analysis of genetic disorders, population genetics, and evolutionary biology. Here, we present an update of our sequencing set to 15,220 Icelanders who we sequenced to an average genome-wide coverage of 34X. We identified 39,020,168 autosomal variants passing GATK filters: 31,079,378 SNPs and 7,940,790 indels. Calling de novo mutations (DNMs) is a formidable challenge given the high false positive rate in sequencing datasets relative to the mutation rate. Here we addressed this issue by using segregation of alleles in three-generation families. Using this transmission assay, we controlled the false positive rate and identified 108,778 high quality DNMs. Furthermore, we used our extended family structure and read pair tracing of DNMs to a panel of phased SNPs, to determine the parent of origin of 42,961 DNMs.

Parental influence on human germline de novo mutations in 1,548 trios from Iceland.

The characterization of mutational processes that generate sequence diversity in the human genome is of paramount importance both to medical genetics and to evolutionary studies. To understand how the age and sex of transmitting parents affect de novo mutations, here we sequence 1,548 Icelanders, their parents, and, for a subset of 225, at least one child, to 35× genome-wide coverage. We find 108,778 de novo mutations, both single nucleotide polymorphisms and indels, and determine the parent of origin of 42,961. The number of de novo mutations from mothers increases by 0.37 per year of age (95% CI 0.32-0.43), a quarter of the 1.51 per year from fathers (95% CI 1.45-1.57). The number of clustered mutations increases faster with the mother's age than with the father's, and the genomic span of maternal de novo mutation clusters is greater than that of paternal ones. The types of de novo mutation from mothers change substantially with age, with a 0.26% (95% CI 0.19-0.33%) decrease in cytosine-phosphate-guanine to thymine-phosphate-guanine (CpG>TpG) de novo mutations and a 0.33% (95% CI 0.28-0.38%) increase in C>G de novo mutations per year, respectively. Remarkably, these age-related changes are not distributed uniformly across the genome. A striking example is a 20 megabase region on chromosome 8p, with a maternal C>G mutation rate that is up to 50-fold greater than the rest of the genome. The age-related accumulation of maternal non-crossover gene conversions also mostly occurs within these regions. Increased sequence diversity and linkage disequilibrium of C>G variants within regions affected by excess maternal mutations indicate that the underlying mutational process has persisted in humans for thousands of years. Moreover, the regional excess of C>G variation in humans is largely shared by chimpanzees, less by gorillas, and is almost absent from orangutans. This demonstrates that sequence diversity in humans results from evolving interactions between age, sex, mutation type, and genomic location.

Clonal hematopoiesis, with and without candidate driver mutations, is common in the elderly.

Clonal hematopoiesis (CH) arises when a substantial proportion of mature blood cells is derived from a single dominant hematopoietic stem cell lineage. Somatic mutations in candidate driver (CD) genes are thought to be responsible for at least some cases of CH. Using whole-genome sequencing of 11 262 Icelanders, we found 1403 cases of CH by using barcodes of mosaic somatic mutations in peripheral blood, whether or not they have a mutation in a CD gene. We find that CH is very common in the elderly, trending toward inevitability. We show that somatic mutations in TET2, DNMT3A, ASXL1, and PPM1D are associated with CH at high significance. However, known CD mutations were evident in only a fraction of CH cases. Nevertheless, the highly prevalent CH we detect associates with increased mortality rates, risk for hematological malignancy, smoking behavior, telomere length, Y-chromosome loss, and other phenotypic characteristics. Modeling suggests some CH cases could arise in the absence of CD mutations as a result of neutral drift acting on a small population of active hematopoietic stem cells. Finally, we find a germline deletion in intron 3 of the telomerase reverse transcriptase (TERT) gene that predisposes to CH (rs34002450; P = 7.4 × 10-12; odds ratio, 1.37).

A genome-wide association study yields five novel thyroid cancer risk loci.

The great majority of thyroid cancers are of the non-medullary type. Here we report findings from a genome-wide association study of non-medullary thyroid cancer, including in total 3,001 patients and 287,550 controls from five study groups of European descent. Our results yield five novel loci (all with Pcombined<3 × 10-8): 1q42.2 (rs12129938 in PCNXL2), 3q26.2 (rs6793295 a missense mutation in LRCC34 near TERC), 5q22.1 (rs73227498 between NREP and EPB41L4A), 10q24.33 (rs7902587 near OBFC1), and two independently associated variants at 15q22.33 (rs2289261 and rs56062135; both in SMAD3). We also confirm recently published association results from a Chinese study of a variant on 5p15.33 (rs2736100 near the TERT gene) and present a stronger association result for a moderately correlated variant (rs10069690; OR=1.20, P=3.2 × 10-7) based on our study of individuals of European ancestry. In combination, these results raise several opportunities for future studies of the pathogenesis of thyroid cancer.

Insertion of an SVA-E retrotransposon into the CASP8 gene is associated with protection against prostate cancer.

Transcriptional and splicing anomalies have been observed in intron 8 of the CASP8 gene (encoding procaspase-8) in association with cutaneous basal-cell carcinoma (BCC) and linked to a germline SNP rs700635. Here, we show that the rs700635[C] allele, which is associated with increased risk of BCC and breast cancer, is protective against prostate cancer [odds ratio (OR) = 0.91, P = 1.0 × 10(-6)]. rs700635[C] is also associated with failures to correctly splice out CASP8 intron 8 in breast and prostate tumours and in corresponding normal tissues. Investigation of rs700635[C] carriers revealed that they have a human-specific short interspersed element-variable number of tandem repeat-Alu (SINE-VNTR-Alu), subfamily-E retrotransposon (SVA-E) inserted into CASP8 intron 8. The SVA-E shows evidence of prior activity, because it has transduced some CASP8 sequences during subsequent retrotransposition events. Whole-genome sequence (WGS) data were used to tag the SVA-E with a surrogate SNP rs1035142[T] (r(2) = 0.999), which showed associations with both the splicing anomalies (P = 6.5 × 10(-32)) and with protection against prostate cancer (OR = 0.91, P = 3.8 × 10(-7)).

Loss-of-function variants in ATM confer risk of gastric cancer.

Gastric cancer is a serious health problem worldwide, with particularly high prevalence in eastern Asia. Genome-wide association studies (GWAS) in Asian populations have identified several loci that associate with gastric cancer risk. Here we report a GWAS of gastric cancer in a European population, using information on 2,500 population-based gastric cancer cases and 205,652 controls. We found a new gastric cancer association with loss-of-function mutations in ATM (gene test, P = 8.0 × 10(-12); odds ratio (OR) = 4.74). The combination of the loss-of-function variants p.Gln852*, p.Ser644* and p.Tyr103* (combined minor allele frequency (MAF) = 0.3%) also associates with pancreatic and prostate cancers (OR = 3.81 and 2.18, respectively) and gives an indication of risk of breast and colorectal cancers (OR = 1.82 and 1.97, respectively). Cancers in those carrying loss-of-function ATM mutations are diagnosed at a significantly earlier age than in non-carriers. Our results confirm an association between gastric cancer in Europeans and three loci previously reported in Asians, MUC1, PRKAA1 and PSCA, refine the association signal at PRKAA1 and support a pathogenic role for the tandem repeat identified in MUC1.

New basal cell carcinoma susceptibility loci.

In an ongoing screen for DNA sequence variants that confer risk of cutaneous basal cell carcinoma (BCC), we conduct a genome-wide association study (GWAS) of 24,988,228 SNPs and small indels detected through whole-genome sequencing of 2,636 Icelanders and imputed into 4,572 BCC patients and 266,358 controls. Here we show the discovery of four new BCC susceptibility loci: 2p24 MYCN (rs57244888[C], OR=0.76, P=4.7 × 10(-12)), 2q33 CASP8-ALS2CR12 (rs13014235[C], OR=1.15, P=1.5 × 10(-9)), 8q21 ZFHX4 (rs28727938[G], OR=0.70, P=3.5 × 10(-12)) and 10p14 GATA3 (rs73635312[A], OR=0.74, P=2.4 × 10(-16)). Fine mapping reveals that two variants correlated with rs73635312[A] occur in conserved binding sites for the GATA3 transcription factor. In addition, expression microarrays and RNA-seq show that rs13014235[C] and a related SNP rs700635[C] are associated with expression of CASP8 splice variants in which sequences from intron 8 are retained.

Large-scale whole-genome sequencing of the Icelandic population.

Here we describe the insights gained from sequencing the whole genomes of 2,636 Icelanders to a median depth of 20×. We found 20 million SNPs and 1.5 million insertions-deletions (indels). We describe the density and frequency spectra of sequence variants in relation to their functional annotation, gene position, pathway and conservation score. We demonstrate an excess of homozygosity and rare protein-coding variants in Iceland. We imputed these variants into 104,220 individuals down to a minor allele frequency of 0.1% and found a recessive frameshift mutation in MYL4 that causes early-onset atrial fibrillation, several mutations in ABCB4 that increase risk of liver diseases and an intronic variant in GNAS associating with increased thyroid-stimulating hormone levels when maternally inherited. These data provide a study design that can be used to determine how variation in the sequence of the human genome gives rise to human diversity.