Genome sequencing with comprehensive variant calling identifies structural variants and repeat expansions in a large fraction of individuals with ataxia and/or neuromuscular disorders.

Introduction: Neuromuscular disorders (NMDs) have a heterogeneous etiology. A genetic diagnosis is key to personalized healthcare and access to targeted treatment for the affected individuals. Methods: In this study, 861 patients with NMDs were analyzed with genome sequencing and comprehensive variant calling including single nucleotide variants, small insertions/deletions (SNVs/INDELs), and structural variants (SVs) in a panel of 895 NMD genes, as well as short tandem repeat expansions (STRs) at 28 loci. In addition, for unsolved cases with an unspecific clinical presentation, the analysis of a panel with OMIM disease genes was added. Results: In the cohort, 27% (232/861) of the patients harbored pathogenic variants, of which STRs and SVs accounted for one-third of the patients (71/232). The variants were found in 107 different NMD genes. Furthermore, 18 pediatric patients harbored pathogenic variants in non-NMD genes. Discussion: Our results highlight that for children with unspecific hypotonia, a genome-wide analysis rather than a disease-based gene panel should be considered as a diagnostic approach. More importantly, our results clearly show that it is crucial to include STR- and SV-analyses in the diagnostics of patients with neuromuscular disorders.

Genome sequencing is a sensitive first-line test to diagnose individuals with intellectual disability.

PURPOSE: Individuals with intellectual disability (ID) and/or neurodevelopment disorders (NDDs) are currently investigated with several different approaches in clinical genetic diagnostics. METHODS: We compared the results from 3 diagnostic pipelines in patients with ID/NDD: genome sequencing (GS) first (N = 100), GS as a secondary test (N = 129), or chromosomal microarray (CMA) with or without FMR1 analysis (N = 421). RESULTS: The diagnostic yield was 35% (GS-first), 26% (GS as a secondary test), and 11% (CMA/FMR1). Notably, the age of diagnosis was delayed by 1 year when GS was performed as a secondary test and the cost per diagnosed individual was 36% lower with GS first than with CMA/FMR1. Furthermore, 91% of those with a negative result after CMA/FMR1 analysis (338 individuals) have not yet been referred for additional genetic testing and remain undiagnosed. CONCLUSION: Our findings strongly suggest that genome analysis outperforms other testing strategies and should replace traditional CMA and FMR1 analysis as a first-line genetic test in individuals with ID/NDD. GS is a sensitive, time- and cost-effective method that results in a confirmed molecular diagnosis in 35% of all referred patients.

Unravelling genetic variants of a swedish family with high risk of prostate cancer.

BACKGROUND: Prostate cancer is the most prevalent cancer in men worldwide. It is a polygenic disease with a substantial proportion of heritability. Identification of novel candidate biomarkers is crucial for clinical cancer prevention and the development of therapeutic strategies. Here, we describe the analysis of rare and common genetic variants that can predispose to the development of prostate cancer. METHODS: Whole-genome sequencing was performed on germline DNA of five Swedish siblings which were diagnosed with prostate cancer. The high-risk variants were identified setting the minor allele frequency < 0.01, CADD > 10 and if tested in PRACTICAL, OR > 1.5, while the low-risk variants were identified minor allele frequency > 0.01, CADD > 10 and if tested in PRACTICAL, OR > 1.1. RESULTS: We identified 38 candidate high-risk gene variants and 332 candidate low-risk gene variants, where 2 and 14 variants were in coding regions, respectively, that were shared by the brothers with prostate cancer. CONCLUSIONS: This study expanded the knowledge of potential risk factor candidates involved in hereditary and familial prostate cancer. Our findings can be beneficial when applying targeted screening in families with a high risk of developing the disease.

Whole-Exome Sequencing of Germline Variants in Non-BRCA Families with Hereditary Breast Cancer.

Breast cancer is the most prevalent malignancy among women worldwide and hereditary breast cancer (HBC) accounts for about 5−10% of the cases. Today, the most recurrent genes known are BRCA1 and BRCA2, accounting for around 25% of familial cases. Although thousands of loss-of-function variants in more than twenty predisposing genes have been found, the majority of familial cases of HBC remain unexplained. The aim of this study was to identify new predisposing genes for HBC in three non-BRCA families with autosomal dominant inheritance pattern using whole-exome sequencing and functional prediction tools. No pathogenic variants in known hereditary cancer-related genes could explain the breast cancer susceptibility in these families. Among 2122 exonic variants with maximum minor allele frequency (MMAF) < 0.1%, between 17−35 variants with combined annotation-dependent depletion (CADD) > 20 segregated with disease in the three analyzed families. Selected candidate genes, i.e., UBASH3A, MYH13, UTP11L, and PAX7, were further evaluated using protein expression analysis but no alterations of cancer-related pathways were observed. In conclusion, identification of new high-risk cancer genes using whole-exome sequencing has been more challenging than initially anticipated, in spite of selected families with pronounced family history of breast cancer. A combination of low- and intermediate-genetic-risk variants may instead contribute the breast cancer susceptibility in these families.

A Swedish Genome-Wide Haplotype Association Analysis Identifies a Novel Breast Cancer Susceptibility Locus in 8p21.2 and Characterizes Three Loci on Chromosomes 10, 11 and 16.

(1) Background: The heritability of breast cancer is partly explained but much of the genetic contribution remains to be identified. Haplotypes are often used as markers of ethnicity as they are preserved through generations. We have previously demonstrated that haplotype analysis, in addition to standard SNP association studies, could give novel and more detailed information on genetic cancer susceptibility. (2) Methods: In order to examine the association of a SNP or a haplotype to breast cancer risk, we performed a genome wide haplotype association study, using sliding window analysis of window sizes 1−25 and 50 SNPs, in 3200 Swedish breast cancer cases and 5021 controls. (3) Results: We identified a novel breast cancer susceptibility locus in 8p21.1 (OR 2.08; p 3.92 × 10−8), confirmed three known loci in 10q26.13, 11q13.3, 16q12.1-2 and further identified novel subloci within these three loci. Altogether 76 risk SNPs, 3302 risk haplotypes of window size 2−25 and 113 risk haplotypes of window size 50 at p < 5 × 10−8 on chromosomes 8, 10, 11 and 16 were identified. In the known loci haplotype analysis reached an OR of 1.48 in overall breast cancer and in familial cases OR 1.68. (4) Conclusions: Analyzing haplotypes, rather than single variants, could detect novel susceptibility loci even in small study populations but the method requires a fairly homogenous study population.

Sequencing for germline mutations in Swedish breast cancer families reveals novel breast cancer risk genes.

Identifying genetic cancer risk factors will lead to improved genetic counseling, cancer prevention and cancer care. Analyzing families with a strong history of breast cancer (BC) has been a successful method to identify genes that contribute to the disease. This has led to discoveries of high-risk genes like the BRCA-genes. Nevertheless, many BC incidences are of unknown causes. In this study, exome sequencing on 59 BC patients from 24 Swedish families with a strong history of BC was performed to identify variants in known and novel BC predisposing genes. First, we screened known BC genes and identified two pathogenic variants in the BRIP1 and PALB2 genes. Secondly, to identify novel BC genes, rare and high impact variants and segregating in families were analyzed to identify 544 variants in novel BC candidate genes. Of those, 22 variants were defined as high-risk variants. Several interesting genes, either previously linked with BC or in pathways that when flawed could contribute to BC, were among the detected genes. The strongest candidates identified are the FANCM gene, involved in DNA double-strand break repair, and the RAD54L gene, involved in DNA recombination. Our study shows identifying pathogenic variants is challenging despite a strong family history of BC. Several interesting candidates were observed here that need to be further studied.

Identification of known and novel familial cancer genes in Swedish colorectal cancer families.

Identifying new candidate colorectal cancer (CRC) genes and mutations are important for clinical cancer prevention as well as in cancer care. Genetic counseling is already implemented for known high-risk variants; however, the majority of CRC are of unknown causes. In our study, 110 CRC patients in 55 Swedish families with a strong history of CRC but unknown genetic causes were analyzed with the aim of identifying novel candidate CRC predisposing genes. Exome sequencing was used to identify rare and high-impact variants enriched in the families. No clear pathogenic variants were found in known CRC predisposing genes; however, potential pathogenic variants in novel CRC predisposing genes were identified. Over 3000 variants with minor allele frequency (MAF) <0.01 and Combined Annotation Dependent Depletion (CADD) > 20 were seen aggregating in the CRC families. Of those, 27 variants with MAF < 0.001 and CADD>25 were considered high-risk mutations. Interestingly, more than half of the high-risk variants were detected in three families, suggesting cumulating contribution of several variants to CRC. In summary, our study shows that despite a strong history of CRC within families, identifying pathogenic variants is challenging. In a small number of families, few rare mutations were shared by affected family members. This could indicate that in the absence of known CRC predisposing genes, a cumulating contribution of mutations leads to CRC observed in these families.

Whole genome DNA sequencing provides an atlas of somatic mutagenesis in healthy human cells and identifies a tumor-prone cell type.

BACKGROUND: The lifelong accumulation of somatic mutations underlies age-related phenotypes and cancer. Mutagenic forces are thought to shape the genome of aging cells in a tissue-specific way. Whole genome analyses of somatic mutation patterns, based on both types and genomic distribution of variants, can shed light on specific processes active in different human tissues and their effect on the transition to cancer. RESULTS: To analyze somatic mutation patterns, we compile a comprehensive genetic atlas of somatic mutations in healthy human cells. High-confidence variants are obtained from newly generated and publicly available whole genome DNA sequencing data from single non-cancer cells, clonally expanded in vitro. To enable a well-controlled comparison of different cell types, we obtain single genome data (92% mean coverage) from multi-organ biopsies from the same donors. These data show multiple cell types that are protected from mutagens and display a stereotyped mutation profile, despite their origin from different tissues. Conversely, the same tissue harbors cells with distinct mutation profiles associated to different differentiation states. Analyses of mutation rate in the coding and non-coding portions of the genome identify a cell type bearing a unique mutation pattern characterized by mutation enrichment in active chromatin, regulatory, and transcribed regions. CONCLUSIONS: Our analysis of normal cells from healthy donors identifies a somatic mutation landscape that enhances the risk of tumor transformation in a specific cell population from the kidney proximal tubule. This unique pattern is characterized by high rate of mutation accumulation during adult life and specific targeting of expressed genes and regulatory regions.

Somatic mutation that affects transcription factor binding upstream of CD55 in the temporal cortex of a late-onset Alzheimer disease patient.

Alzheimer's disease (AD) is the most common neurodegenerative disease worldwide. Familial cases suggest genetic components; however, monogenetic causes are few, and the vast majority of incidences have unknown cause. Sequencing efforts have focused on germline mutations, but improved technology has opened up for studies on somatic mutations in affected brain tissue samples. Here we use ultra-deep sequencing on brain and blood from early-onset AD (EOAD) and late-onset AD (LOAD) patients and non-AD individuals (n = 16). In total, 2.86 Mb of genomic regions, previously associated with AD, were targeted included 28 genes and upstream and downstream regulatory regions. Tailored downstream bioinformatics filtering identified 11 somatic single nucleotide variants in the temporal cortex in AD patients and none in the controls. One variant was validated to be present at 0.4% allele frequency in temporal cortex of a LOAD patient. This variant was predicted to affect transcription factor binding sites upstream of the CD55 gene, contributing to AD pathogenesis by affecting the complement system. Our results suggest that future studies targeting larger portions of the genome for somatic mutation analysis are important to obtain an increased understanding for the molecular basis of both EOAD and LOAD.

Accumulation of Progerin Affects the Symmetry of Cell Division and Is Associated with Impaired Wnt Signaling and the Mislocalization of Nuclear Envelope Proteins.

Hutchinson-Gilford progeria syndrome (HGPS) is the result of a defective form of the lamin A protein called progerin. While progerin is known to disrupt the properties of the nuclear lamina, the underlying mechanisms responsible for the pathophysiology of HGPS remain less clear. Previous studies in our laboratory have shown that progerin expression in murine epidermal basal cells results in impaired stratification and halted development of the skin. Stratification and differentiation of the epidermis is regulated by asymmetric stem cell division. Here, we show that expression of progerin impairs the ability of stem cells to maintain tissue homeostasis as a result of altered cell division. Quantification of basal skin cells showed an increase in symmetric cell division that correlated with progerin accumulation in HGPS mice. Investigation of the mechanisms underlying this phenomenon revealed a putative role of Wnt/ß-catenin signaling. Further analysis suggested an alteration in the nuclear translocation of ß-catenin involving the inner and outer nuclear membrane proteins, emerin and nesprin-2. Taken together, our results suggest a direct involvement of progerin in the transmission of Wnt signaling and normal stem cell division. These insights into the molecular mechanisms of progerin may help develop new treatment strategies for HGPS.

Somatic mutagenesis in satellite cells associates with human skeletal muscle aging.

Human aging is associated with a decline in skeletal muscle (SkM) function and a reduction in the number and activity of satellite cells (SCs), the resident stem cells. To study the connection between SC aging and muscle impairment, we analyze the whole genome of single SC clones of the leg muscle vastus lateralis from healthy individuals of different ages (21-78 years). We find an accumulation rate of 13 somatic mutations per genome per year, consistent with proliferation of SCs in the healthy adult muscle. SkM-expressed genes are protected from mutations, but aging results in an increase in mutations in exons and promoters, targeting genes involved in SC activity and muscle function. In agreement with SC mutations affecting the whole tissue, we detect a missense mutation in a SC propagating to the muscle. Our results suggest somatic mutagenesis in SCs as a driving force in the age-related decline of SkM function.

HLA class II sequence variants influence tuberculosis risk in populations of European ancestry.

Mycobacterium tuberculosis infections cause 9 million new tuberculosis cases and 1.5 million deaths annually. To identify variants conferring risk of tuberculosis, we tested 28.3 million variants identified through whole-genome sequencing of 2,636 Icelanders for association with tuberculosis (8,162 cases and 277,643 controls), pulmonary tuberculosis (PTB) and M. tuberculosis infection. We found association of three variants in the region harboring genes encoding the class II human leukocyte antigens (HLAs): rs557011[T] (minor allele frequency (MAF) = 40.2%), associated with M. tuberculosis infection (odds ratio (OR) = 1.14, P = 3.1 × 10(-13)) and PTB (OR = 1.25, P = 5.8 × 10(-12)), and rs9271378[G] (MAF = 32.5%), associated with PTB (OR = 0.78, P = 2.5 × 10(-12))--both located between HLA-DQA1 and HLA-DRB1--and a missense variant encoding p.Ala210Thr in HLA-DQA1 (MAF = 19.1%, rs9272785), associated with M. tuberculosis infection (P = 9.3 × 10(-9), OR = 1.14). We replicated association of these variants with PTB in samples of European ancestry from Russia and Croatia (P < 5.9 × 10(-4)). These findings show that the HLA class II region contributes to genetic risk of tuberculosis, possibly through reduced presentation of protective M. tuberculosis antigens to T cells.

Heterozygote carriers for CNVs in PARK2 are at increased risk of Parkinson's disease.

Together with point mutations, homozygous deletions or duplications in PARK2 are responsible for the majority of autosomal recessive juvenile Parkinsonism. It is debated, however, whether heterozygous carriers of these mutations are at increased risk of Parkinson's disease (PD). Our goal was to determine whether heterozygous carriers of copy number variants (CNVs) affecting exons of the PARK2 gene are at risk of PD that is greater than that of non-carriers. We searched for CNVs affecting exons of PARK2 in a sample of 105 749 genotyped Icelanders. In total, 989 carriers, including 24 diagnosed with PD, were identified. The heterozygous carriers were tested for association in a sample of 1415 PD patients and 40 474 controls ≥65 years of age. PD patients were more often heterozygous carriers of PARK2 CNVs than controls [odds ratio (OR) = 1.69, P = 0.03] and compound heterozygous PD patients for a CNV and a missense mutation were not found. Furthermore, we conducted a meta-analysis of studies reporting on case-control samples screened for heterozygous PARK2 CNVs. Ten studies were included in the final analysis, with 4538 cases and 4213 controls. The pooled OR and P-value for the published and Icelandic results showed significant association between PARK2 CNVs and risk of PD (OR = 2.11, P = 2.54 × 10(-6)). Our analysis shows that heterozygous carriers of CNVs affecting exons of PARK2 have greater risk of PD than non-carriers.

Severe osteoarthritis of the hand associates with common variants within the ALDH1A2 gene and with rare variants at 1p31.

Osteoarthritis is the most common form of arthritis and is a major cause of pain and disability in the elderly. To search for sequence variants that confer risk of osteoarthritis of the hand, we carried out a genome-wide association study (GWAS) in subjects with severe hand osteoarthritis, using variants identified through the whole-genome sequencing of 2,230 Icelanders. We found two significantly associated loci in the Icelandic discovery set: at 15q22 (frequency of 50.7%, odds ratio (OR) = 1.51, P = 3.99 × 10(-10)) in the ALDH1A2 gene and at 1p31 (frequency of 0.02%, OR = 50.6, P = 9.8 × 10(-10)). Among the carriers of the variant at 1p31 is a family with several members in whom the risk allele segregates with osteoarthritis. The variants within the ALDH1A2 gene were confirmed in replication sets from The Netherlands and the UK, yielding an overall association of OR = 1.46 and P = 1.1 × 10(-11) (rs3204689).

Identification of low-frequency and rare sequence variants associated with elevated or reduced risk of type 2 diabetes.

Through whole-genome sequencing of 2,630 Icelanders and imputation into 11,114 Icelandic cases and 267,140 controls followed by testing in Danish and Iranian samples, we discovered 4 previously unreported variants affecting risk of type 2 diabetes (T2D). A low-frequency (1.47%) variant in intron 1 of CCND2, rs76895963[G], reduces risk of T2D by half (odds ratio (OR) = 0.53, P = 5.0 × 10(-21)) and is correlated with increased CCND2 expression. Notably, this variant is also associated with both greater height and higher body mass index (1.17 cm per allele, P = 5.5 × 10(-12) and 0.56 kg/m(2) per allele, P = 6.5 × 10(-7), respectively). In addition, two missense variants in PAM, encoding p.Asp563Gly (frequency of 4.98%) and p.Ser539Trp (frequency of 0.65%), confer moderately higher risk of T2D (OR = 1.23, P = 3.9 × 10(-10) and OR = 1.47, P = 1.7 × 10(-5), respectively), and a rare (0.20%) frameshift variant in PDX1, encoding p.Gly218Alafs*12, associates with high risk of T2D (OR = 2.27, P = 7.3 × 10(-7)).

Germline sequence variants in TGM3 and RGS22 confer risk of basal cell carcinoma.

To search for new sequence variants that confer risk of cutaneous basal cell carcinoma (BCC), we conducted a genome-wide association study of 38.5 million single nucleotide polymorphisms (SNPs) and small indels identified through whole-genome sequencing of 2230 Icelanders. We imputed genotypes for 4208 BCC patients and 109 408 controls using Illumina SNP chip typing data, carried out association tests and replicated the findings in independent population samples. We found new BCC susceptibility loci at TGM3 (rs214782[G], P = 5.5 × 10(-17), OR = 1.29) and RGS22 (rs7006527[C], P = 8.7 × 10(-13), OR = 0.77). TGM3 encodes transglutaminase type 3, which plays a key role in production of the cornified envelope during epidermal differentiation.

A meta-analysis of genome-wide association studies identifies novel variants associated with osteoarthritis of the hip.

OBJECTIVES: Osteoarthritis (OA) is the most common form of arthritis with a clear genetic component. To identify novel loci associated with hip OA we performed a meta-analysis of genome-wide association studies (GWAS) on European subjects. METHODS: We performed a two-stage meta-analysis on more than 78,000 participants. In stage 1, we synthesised data from eight GWAS whereas data from 10 centres were used for 'in silico' or 'de novo' replication. Besides the main analysis, a stratified by sex analysis was performed to detect possible sex-specific signals. Meta-analysis was performed using inverse-variance fixed effects models. A random effects approach was also used. RESULTS: We accumulated 11,277 cases of radiographic and symptomatic hip OA. We prioritised eight single nucleotide polymorphism (SNPs) for follow-up in the discovery stage (4349 OA cases); five from the combined analysis, two male specific and one female specific. One locus, at 20q13, represented by rs6094710 (minor allele frequency (MAF) 4%) near the NCOA3 (nuclear receptor coactivator 3) gene, reached genome-wide significance level with p=7.9×10(-9) and OR=1.28 (95% CI 1.18 to 1.39) in the combined analysis of discovery (p=5.6×10(-8)) and follow-up studies (p=7.3×10(-4)). We showed that this gene is expressed in articular cartilage and its expression was significantly reduced in OA-affected cartilage. Moreover, two loci remained suggestive associated; rs5009270 at 7q31 (MAF 30%, p=9.9×10(-7), OR=1.10) and rs3757837 at 7p13 (MAF 6%, p=2.2×10(-6), OR=1.27 in male specific analysis). CONCLUSIONS: Novel genetic loci for hip OA were found in this meta-analysis of GWAS.

Nonsense mutation in the LGR4 gene is associated with several human diseases and other traits.

Low bone mineral density (BMD) is used as a parameter of osteoporosis. Genome-wide association studies of BMD have hitherto focused on BMD as a quantitative trait, yielding common variants of small effects that contribute to the population diversity in BMD. Here we use BMD as a dichotomous trait, searching for variants that may have a direct effect on the risk of pathologically low BMD rather than on the regulation of BMD in the healthy population. Through whole-genome sequencing of Icelandic individuals, we found a rare nonsense mutation within the leucine-rich-repeat-containing G-protein-coupled receptor 4 (LGR4) gene (c.376C>T) that is strongly associated with low BMD, and with osteoporotic fractures. This mutation leads to termination of LGR4 at position 126 and fully disrupts its function. The c.376C>T mutation is also associated with electrolyte imbalance, late onset of menarche and reduced testosterone levels, as well as an increased risk of squamous cell carcinoma of the skin and biliary tract cancer. Interestingly, the phenotype of carriers of the c.376C>T mutation overlaps that of Lgr4 mutant mice.

Identification of low-frequency variants associated with gout and serum uric acid levels.

We tested 16 million SNPs, identified through whole-genome sequencing of 457 Icelanders, for association with gout and serum uric acid levels. Genotypes were imputed into 41,675 chip-genotyped Icelanders and their relatives, for effective sample sizes of 968 individuals with gout and 15,506 individuals for whom serum uric acid measurements were available. We identified a low-frequency missense variant (c.1580C>G) in ALDH16A1 associated with gout (OR = 3.12, P = 1.5 × 10(-16), at-risk allele frequency = 0.019) and serum uric acid levels (effect = 0.36 s.d., P = 4.5 × 10(-21)). We confirmed the association with gout by performing Sanger sequencing on 6,017 Icelanders. The association with gout was stronger in males relative to females. We also found a second variant on chromosome 1 associated with gout (OR = 1.92, P = 0.046, at-risk allele frequency = 0.986) and serum uric acid levels (effect = 0.48 s.d., P = 4.5 × 10(-16)). This variant is close to a common variant previously associated with serum uric acid levels. This work illustrates how whole-genome sequencing data allow the detection of associations between low-frequency variants and complex traits.

Mutations in BRIP1 confer high risk of ovarian cancer.

Ovarian cancer causes more deaths than any other gynecologic malignancy in developed countries. Sixteen million sequence variants, identified through whole-genome sequencing of 457 Icelanders, were imputed to 41,675 Icelanders genotyped using SNP chips, as well as to their relatives. Sequence variants were tested for association with ovarian cancer (N of affected individuals = 656). We discovered a rare (0.41% allelic frequency) frameshift mutation, c.2040_2041insTT, in the BRIP1 (FANCJ) gene that confers an increase in ovarian cancer risk (odds ratio (OR) = 8.13, P = 2.8 × 10(-14)). The mutation was also associated with increased risk of cancer in general and reduced lifespan by 3.6 years. In a Spanish population, another frameshift mutation in BRIP1, c.1702_1703del, was seen in 2 out of 144 subjects with ovarian cancer and 1 out of 1,780 control subjects (P = 0.016). This allele was also associated with breast cancer (seen in 6/927 cases; P = 0.0079). Ovarian tumors from heterozygous carriers of the Icelandic mutation show loss of the wild-type allele, indicating that BRIP1 behaves like a classical tumor suppressor gene in ovarian cancer.