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 Lynch syndrome risk variants in the Romanian population.

Two familial forms of colorectal cancer (CRC), Lynch syndrome (LS) and familial adenomatous polyposis (FAP), are caused by rare mutations in DNA mismatch repair genes (MLH1, MSH2, MSH6, PMS2) and the genes APC and MUTYH, respectively. No information is available on the presence of high-risk CRC mutations in the Romanian population. We performed whole-genome sequencing of 61 Romanian CRC cases with a family history of cancer and/or early onset of disease, focusing the analysis on candidate variants in the LS and FAP genes. The frequencies of all candidate variants were assessed in a cohort of 688 CRC cases and 4567 controls. Immunohistochemical (IHC) staining for MLH1, MSH2, MSH6, and PMS2 was performed on tumour tissue. We identified 11 candidate variants in 11 cases; six variants in MLH1, one in MSH6, one in PMS2, and three in APC. Combining information on the predicted impact of the variants on the proteins, IHC results and previous reports, we found three novel pathogenic variants (MLH1:p.Lys84ThrfsTer4, MLH1:p.Ala586CysfsTer7, PMS2:p.Arg211ThrfsTer38), and two novel variants that are unlikely to be pathogenic. Also, we confirmed three previously published pathogenic LS variants and suggest to reclassify a previously reported variant of uncertain significance to pathogenic (MLH1:c.1559-1G>C).

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

A frameshift deletion in the sarcomere gene MYL4 causes early-onset familial atrial fibrillation.

AIMS: Atrial fibrillation (AF) is the most common sustained cardiac arrhythmia in man, causing substantial morbidity and mortality with a major worldwide public health impact. It is increasingly recognized as a highly heritable condition. This study aimed to determine genetic risk factors for early-onset AF. METHODS AND RESULTS: We sequenced the whole genomes of 8453 Icelanders and imputed genotypes of the 25.5 million sequence variants we discovered into 1799 Icelanders with early-onset AF (diagnosed before 60 years of age) and 337 453 controls. Each sequence variant was tested for association based on multiplicative and recessive inheritance models. We discovered a rare frameshift deletion in the myosin MYL4 gene (c.234delC) that associates with early-onset AF under a recessive mode of inheritance (allelic frequency = 0.58%). We found eight homozygous carriers of the mutation, all of whom had early-onset AF. Six of the homozygotes were diagnosed by the age of 30 and the remaining two in their 50s. Three of the homozygotes had received pacemaker implantations due to sick sinus syndrome, three had suffered an ischemic stroke, and one suffered sudden cardiac death. CONCLUSIONS: Through a population approach we found a loss of function mutation in the myosin gene MYL4 that, in the homozygous state, is completely penetrant for early-onset AF. The finding may provide novel mechanistic insight into the pathophysiology of this complex arrhythmia.

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.

Genome-wide association study yields variants at 20p12.2 that associate with urinary bladder cancer.

Genome-wide association studies (GWAS) of urinary bladder cancer (UBC) have yielded common variants at 12 loci that associate with risk of the disease. We report here the results of a GWAS of UBC including 1670 UBC cases and 90 180 controls, followed by replication analysis in additional 5266 UBC cases and 10 456 controls. We tested a dataset containing 34.2 million variants, generated by imputation based on whole-genome sequencing of 2230 Icelanders. Several correlated variants at 20p12, represented by rs62185668, show genome-wide significant association with UBC after combining discovery and replication results (OR = 1.19, P = 1.5 × 10(-11) for rs62185668-A, minor allele frequency = 23.6%). The variants are located in a non-coding region approximately 300 kb upstream from the JAG1 gene, an important component of the Notch signaling pathways that may be oncogenic or tumor suppressive in several forms of cancer. Our results add to the growing number of UBC risk variants discovered through GWAS.

A genome-wide copy number association study of osteoporotic fractures points to the 6p25.1 locus.

BACKGROUND: Osteoporosis is a systemic skeletal disease characterised by reduced bone mineral density and increased susceptibility to fracture; these traits are highly heritable. Both common and rare copy number variants (CNVs) potentially affect the function of genes and may influence disease risk. AIM: To identify CNVs associated with osteoporotic bone fracture risk. METHOD: We performed a genome-wide CNV association study in 5178 individuals from a prospective cohort in the Netherlands, including 809 osteoporotic fracture cases, and performed in silico lookups and de novo genotyping to replicate in several independent studies. RESULTS: A rare (population prevalence 0.14%, 95% CI 0.03% to 0.24%) 210 kb deletion located on chromosome 6p25.1 was associated with the risk of fracture (OR 32.58, 95% CI 3.95 to 1488.89; p = 8.69 × 10(-5)). We performed an in silico meta-analysis in four studies with CNV microarray data and the association with fracture risk was replicated (OR 3.11, 95% CI 1.01 to 8.22; p = 0.02). The prevalence of this deletion showed geographic diversity, being absent in additional samples from Australia, Canada, Poland, Iceland, Denmark, and Sweden, but present in the Netherlands (0.34%), Spain (0.33%), USA (0.23%), England (0.15%), Scotland (0.10%), and Ireland (0.06%), with insufficient evidence for association with fracture risk. CONCLUSIONS: These results suggest that deletions in the 6p25.1 locus may predispose to higher risk of fracture in a subset of populations of European origin; larger and geographically restricted studies will be needed to confirm this regional association. This is a first step towards the evaluation of the role of rare CNVs in osteoporosis.

European genome-wide association study identifies SLC14A1 as a new urinary bladder cancer susceptibility gene.

Three genome-wide association studies in Europe and the USA have reported eight urinary bladder cancer (UBC) susceptibility loci. Using extended case and control series and 1000 Genomes imputations of 5 340 737 single-nucleotide polymorphisms (SNPs), we searched for additional loci in the European GWAS. The discovery sample set consisted of 1631 cases and 3822 controls from the Netherlands and 603 cases and 37 781 controls from Iceland. For follow-up, we used 3790 cases and 7507 controls from 13 sample sets of European and Iranian ancestry. Based on the discovery analysis, we followed up signals in the urea transporter (UT) gene SLC14A. The strongest signal at this locus was represented by a SNP in intron 3, rs17674580, that reached genome-wide significance in the overall analysis of the discovery and follow-up groups: odds ratio = 1.17, P = 7.6 × 10(-11). SLC14A1 codes for UTs that define the Kidd blood group and are crucial for the maintenance of a constant urea concentration gradient in the renal medulla and, through this, the kidney's ability to concentrate urine. It is speculated that rs17674580, or other sequence variants in LD with it, indirectly modifies UBC risk by affecting urine production. If confirmed, this would support the 'urogenous contact hypothesis' that urine production and voiding frequency modify the risk of UBC.

Efficacy and safety of erlotinib versus chemotherapy in second-line treatment of patients with advanced, non-small-cell lung cancer with poor prognosis (TITAN): a randomised multicentre, open-label, phase 3 study.

BACKGROUND: Erlotinib, docetaxel, and pemetrexed are approved for the second-line treatment of non-small-cell lung cancer (NSCLC), but no head-to-head data from large clinical trials are available. We undertook the Tarceva In Treatment of Advanced NSCLC (TITAN) study to assess the efficacy and tolerability of second-line erlotinib versus chemotherapy in patients with refractory NSCLC. METHODS: TITAN was an international, randomised multicentre, open-label, phase 3 study that was done at 77 sites in 24 countries. Chemotherapy-naive patients with locally advanced, recurrent, or metastatic NSCLC received up to four cycles of first-line platinum doublet chemotherapy, after which patients with disease progression during or immediately after chemotherapy were offered enrolment into TITAN. Enrolled patients were randomly assigned (1:1) by a minimisation method to ensure balanced stratification, to receive erlotinib 150 mg/day or chemotherapy (standard docetaxel or pemetrexed regimens, at the treating investigators' discretion), until unacceptable toxicity, disease progression, or death. Patients were stratified by disease stage, Eastern Cooperative Oncology Group performance status, smoking history, and region of residence. The primary endpoint was overall survival in the intention-to-treat population. TITAN was halted prematurely because of slow recruitment. This study is registered with ClinicalTrials.gov, number NCT00556322. FINDINGS: Between April 10, 2006, and Feb 24, 2010, 2590 chemotherapy-naive patients were treated with first-line platinum doublet chemotherapy, of whom 424 had disease progression and were enrolled into TITAN. 203 patients were randomly assigned to receive erlotinib and 221 were assigned to receive chemotherapy. Median follow-up was 27·9 months (IQR 11·0-36·0) in the erlotinib group and 24·8 months (12·1-41·6) in the chemotherapy group. Median overall survival was 5·3 months (95% CI 4·0-6·0) with erlotinib and 5·5 months (4·4-7·1) with chemotherapy (hazard ratio [HR] 0·96, 95% CI 0·78-1·19; log-rank p=0·73). The adverse-event profile of each group was in line with previous studies. Rash (98/196 [50%] in the erlotinib group vs 10/213 [5%] in the chemotherapy group for all grades; nine [5%] vs none for grade 3 or 4) and diarrhoea (36 [18%] vs four [2%] for all grades; five [3%] vs none for grade 3 or 4) were the most common treatment-related adverse events with erlotinib, whereas alopecia (none vs 23 [11%] for all grades; none vs one [<1%] for grade 3/4) was the most common treatment-related adverse event with chemotherapy. INTERPRETATION: No significant differences in efficacy were noted between patients treated with erlotinib and those treated with docetaxel or pemetrexed. Since the toxicity profiles of erlotinib and chemotherapy differ, second-line treatment decisions should take into account patient preference and specific toxicity risk profiles. FUNDING: F Hoffmann-La Roche.

Host cell factors and functions involved in vesicular stomatitis virus entry.

Vesicular stomatitis virus (VSV) is an animal virus that based on electron microscopy and its dependence on acidic cellular compartments for infection is thought to enter its host cells in a clathrin-dependent manner. The exact cellular mechanism, however, is largely unknown. In this study, we characterized the entry kinetics of VSV and elucidated viral requirements for host cell factors during infection in HeLa cells. We found that endocytosis of VSV was a fast process with a half time of 2.5 to 3 min and that acid activation occurred within 1 to 2 min after internalization in early endosomes. The majority of viral particles were endocytosed in a clathrin-based, dynamin-2-dependent manner. Although associated with some of the surface-bound viruses, the classical adaptor protein complex AP-2 was not required for infection. Time-lapse microscopy revealed that the virus either entered preformed clathrin-coated pits or induced de novo formation of pits. Dynamin-2 was recruited to plasma membrane-confined virus particles. Thus, VSV can induce productive internalization by exploiting a specific combination of the clathrin-associated proteins and cellular functions.

Genetic predisposition to hypertension is associated with preeclampsia in European and Central Asian women.

Preeclampsia is a serious complication of pregnancy, affecting both maternal and fetal health. In genome-wide association meta-analysis of European and Central Asian mothers, we identify sequence variants that associate with preeclampsia in the maternal genome at ZNF831/20q13 and FTO/16q12. These are previously established variants for blood pressure (BP) and the FTO variant has also been associated with body mass index (BMI). Further analysis of BP variants establishes that variants at MECOM/3q26, FGF5/4q21 and SH2B3/12q24 also associate with preeclampsia through the maternal genome. We further show that a polygenic risk score for hypertension associates with preeclampsia. However, comparison with gestational hypertension indicates that additional factors modify the risk of preeclampsia.

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.

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.

Whole-genome sequencing identifies rare genotypes in COMP and CHADL associated with high risk of hip osteoarthritis.

We performed a genome-wide association study of total hip replacements, based on variants identified through whole-genome sequencing, which included 4,657 Icelandic patients and 207,514 population controls. We discovered two rare signals that strongly associate with osteoarthritis total hip replacement: a missense variant, c.1141G>C (p.Asp369His), in the COMP gene (allelic frequency = 0.026%, P = 4.0 × 10-12, odds ratio (OR) = 16.7) and a frameshift mutation, rs532464664 (p.Val330Glyfs*106), in the CHADL gene that associates through a recessive mode of inheritance (homozygote frequency = 0.15%, P = 4.5 × 10-18, OR = 7.71). On average, c.1141G>C heterozygotes and individuals homozygous for rs532464664 had their hip replacement operation 13.5 years and 4.9 years earlier than others (P = 0.0020 and P = 0.0026), respectively. We show that the full-length CHADL transcript is expressed in cartilage. Furthermore, the premature stop codon introduced by the CHADL frameshift mutation results in nonsense-mediated decay of the mutant transcripts.

Lack of HIN-1 methylation in BRCA1-linked and "BRCA1-like" breast tumors.

We recently identified a candidate tumor suppressor gene, HIN-1, that is silenced due to methylation in the majority of sporadic breast carcinomas and is localized to 5q33-qter, an area frequently lost in BRCA1 tumors and thought to harbor a BRCA1 modifier gene. To establish whether germ-line mutations in HIN-1 may influence breast cancer risk, we sequenced the HIN-1 coding region in 10 familial breast cancer patients with positive logarithm of the odds scores of at least one of the markers flanking HIN-1. We also sequenced the HIN-1 coding region in 15 BRCA1 and 35 sporadic breast tumors to determine whether HIN-1 is the target of the frequent 5q loss in BRCA1 tumors. No sequence alterations were found in any of the cases analyzed. However, analysis of HIN-1 promoter methylation status revealed that in striking contrast to sporadic cases, there is a nearly complete lack of HIN-1 methylation in BRCA1 tumors (P < 0.0001). Sporadic breast tumors with a "BRCA1-like" histopathological phenotype also demonstrated significantly lower frequency of HIN-1 promoter methylation (P = 0.01) compared with other cancer types, and there was also a difference among tumors based on their estrogen receptor and HER2 status (P = 0.006), suggesting that HIN-1 methylation patterns are associated with specific breast cancer subtypes.

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.

Loss-of-function variants in ABCA7 confer risk of Alzheimer's disease.

We conducted a search for rare, functional variants altering susceptibility to Alzheimer's disease that exploited knowledge of common variants associated with the same disease. We found that loss-of-function variants in ABCA7 confer risk of Alzheimer's disease in Icelanders (odds ratio (OR) = 2.12, P = 2.2 × 10(-13)) and discovered that the association replicated in study groups from Europe and the United States (combined OR = 2.03, P = 6.8 × 10(-15)).

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.

Deletions on chromosome 4 in sporadic and BRCA mutated tumors and association with pathological variables.

BACKGROUND: Chromosomal aberrations in breast tumors from BRCA1 and BRCA2 germ-line mutation carriers are considerably more frequent than what is seen in sporadic breast tumors. According to Comparative Genomic Hybridisation analysis (CGH), deletions on chromosome 4 are one of the most frequent events in BRCA1-associated tumors, suggesting inactivation of specific tumor suppressor genes. MATERIALS AND METHODS: In the present study, 16 microsatellite markers covering chromosome 4 were used to map loss of heterozygosity (LOH) in tumors from BRCA1 (n=41) as well as in tumors from BRCA2 (n=66) mutation carriers and in tumors from unselected cases of breast cancer (n =68). RESULTS: The frequency of LOH in these groups ranged from 16-73% in BRCA1-associated tumors, 13-42% in BRCA2-associated tumors and 8-33% in unselected tumors. LOH was significantly more frequent in BRCA1-associated tumors as compared to BRCA2-associated tumors and unselected tumors, and particularly high (over 70%) at 4q35.2. Pathological variables that were found significantly associated (p< or =0.05) with LOH at specific markers were: high percentage of cells in S-phase, negative estrogen receptor status, young age at diagnosis and large tumors. Deletion mapping indicates the existence of seven non-overlapping regions at chromosome 4, which were identified in all three groups of tumors. Three of these seven regions, 4p16.3-p16.1, 4q27-q32.1 and 4q35.1-4qter, have not been reported in breast cancer previously. CONCLUSION: The results manifest the frequent alterations of chromosome 4 in BRCA1-associated breast tumors and indicate the location of several genes of potential importance in breast cancer development.