Boosting GWAS using biological networks: A study on susceptibility to familial breast cancer.

Genome-wide association studies (GWAS) explore the genetic causes of complex diseases. However, classical approaches ignore the biological context of the genetic variants and genes under study. To address this shortcoming, one can use biological networks, which model functional relationships, to search for functionally related susceptibility loci. Many such network methods exist, each arising from different mathematical frameworks, pre-processing steps, and assumptions about the network properties of the susceptibility mechanism. Unsurprisingly, this results in disparate solutions. To explore how to exploit these heterogeneous approaches, we selected six network methods and applied them to GENESIS, a nationwide French study on familial breast cancer. First, we verified that network methods recovered more interpretable results than a standard GWAS. We addressed the heterogeneity of their solutions by studying their overlap, computing what we called the consensus. The key gene in this consensus solution was COPS5, a gene related to multiple cancer hallmarks. Another issue we observed was that network methods were unstable, selecting very different genes on different subsamples of GENESIS. Therefore, we proposed a stable consensus solution formed by the 68 genes most consistently selected across multiple subsamples. This solution was also enriched in genes known to be associated with breast cancer susceptibility (BLM, CASP8, CASP10, DNAJC1, FGFR2, MRPS30, and SLC4A7, P-value = 3 × 10-4). The most connected gene was CUL3, a regulator of several genes linked to cancer progression. Lastly, we evaluated the biases of each method and the impact of their parameters on the outcome. In general, network methods preferred highly connected genes, even after random rewirings that stripped the connections of any biological meaning. In conclusion, we present the advantages of network-guided GWAS, characterize their shortcomings, and provide strategies to address them. To compute the consensus networks, implementations of all six methods are available at https://github.com/hclimente/gwas-tools.

Gene- and pathway-level analyses of iCOGS variants highlight novel signaling pathways underlying familial breast cancer susceptibility.

Single-nucleotide polymorphisms (SNPs) in over 180 loci have been associated with breast cancer (BC) through genome-wide association studies involving mostly unselected population-based case-control series. Some of them modify BC risk of women carrying a BRCA1 or BRCA2 (BRCA1/2) mutation and may also explain BC risk variability in BC-prone families with no BRCA1/2 mutation. Here, we assessed the contribution of SNPs of the iCOGS array in GENESIS consisting of BC cases with no BRCA1/2 mutation and a sister with BC, and population controls. Genotyping data were available for 1281 index cases, 731 sisters with BC, 457 unaffected sisters and 1272 controls. In addition to the standard SNP-level analysis using index cases and controls, we performed pedigree-based association tests to capture transmission information in the sibships. We also performed gene- and pathway-level analyses to maximize the power to detect associations with lower-frequency SNPs or those with modest effect sizes. While SNP-level analyses identified 18 loci, gene-level analyses identified 112 genes. Furthermore, 31 Kyoto Encyclopedia of Genes and Genomes and 7 Atlas of Cancer Signaling Network pathways were highlighted (false discovery rate of 5%). Using results from the "index case-control" analysis, we built pathway-derived polygenic risk scores (PRS) and assessed their performance in the population-based CECILE study and in a data set composed of GENESIS-affected sisters and CECILE controls. Although these PRS had poor predictive value in the general population, they performed better than a PRS built using our SNP-level findings, and we found that the joint effect of family history and PRS needs to be considered in risk prediction models.

WDR34, a candidate gene for non-syndromic rod-cone dystrophy.

Rod-cone dystrophy (RCD), also called retinitis pigmentosa, is characterized by rod followed by cone photoreceptor degeneration, leading to gradual visual loss. Mutations in over 65 genes have been associated with non-syndromic RCD explaining 60% to 70% of cases, with novel gene defects possibly accounting for the unsolved cases. Homozygosity mapping and whole-exome sequencing applied to a case of autosomal recessive non-syndromic RCD from a consanguineous union identified a homozygous variant in WDR34. Mutations in WDR34 have been previously associated with severe ciliopathy syndromes possibly associated with a retinal dystrophy. This is the first report of a homozygous mutation in WDR34 associated with non-syndromic RCD.

Familial breast cancer and DNA repair genes: Insights into known and novel susceptibility genes from the GENESIS study, and implications for multigene panel testing.

Pathogenic variants in BRCA1 and BRCA2 only explain the underlying genetic cause of about 10% of hereditary breast and ovarian cancer families. Because of cost-effectiveness, multigene panel testing is often performed even if the clinical utility of testing most of the genes remains questionable. The purpose of our study was to assess the contribution of rare, deleterious-predicted variants in DNA repair genes in familial breast cancer (BC) in a well-characterized and homogeneous population. We analyzed 113 DNA repair genes selected from either an exome sequencing or a candidate gene approach in the GENESIS study, which includes familial BC cases with no BRCA1 or BRCA2 mutation and having a sister with BC (N = 1,207), and general population controls (N = 1,199). Sequencing data were filtered for rare loss-of-function variants (LoF) and likely deleterious missense variants (MV). We confirmed associations between LoF and MV in PALB2, ATM and CHEK2 and BC occurrence. We also identified for the first time associations between FANCI, MAST1, POLH and RTEL1 and BC susceptibility. Unlike other associated genes, carriers of an ATM LoF had a significantly higher risk of developing BC than carriers of an ATM MV (ORLoF = 17.4 vs. ORMV = 1.6; p Het = 0.002). Hence, our approach allowed us to specify BC relative risks associated with deleterious-predicted variants in PALB2, ATM and CHEK2 and to add MAST1, POLH, RTEL1 and FANCI to the list of DNA repair genes possibly involved in BC susceptibility. We also highlight that different types of variants within the same gene can lead to different risk estimates.

Mutations in IFT172 cause isolated retinal degeneration and Bardet-Biedl syndrome.

Primary cilia are sensory organelles present on most mammalian cells. The assembly and maintenance of primary cilia are facilitated by intraflagellar transport (IFT), a bidirectional protein trafficking along the cilium. Mutations in genes coding for IFT components have been associated with a group of diseases called ciliopathies. These genetic disorders can affect a variety of organs including the retina. Using whole exome sequencing in three families, we identified mutations in Intraflagellar Transport 172 Homolog [IFT172 (Chlamydomonas)] that underlie an isolated retinal degeneration and Bardet-Biedl syndrome. Extensive functional analyses of the identified mutations in cell culture, rat retina and in zebrafish demonstrated their hypomorphic or null nature. It has recently been reported that mutations in IFT172 cause a severe ciliopathy syndrome involving skeletal, renal, hepatic and retinal abnormalities (Jeune and Mainzer-Saldino syndromes). Here, we report for the first time that mutations in this gene can also lead to an isolated form of retinal degeneration. The functional data for the mutations can partially explain milder phenotypes; however, the involvement of modifying alleles in the IFT172-associated phenotypes cannot be excluded. These findings expand the spectrum of disease associated with mutations in IFT172 and suggest that mutations in genes originally reported to be associated with syndromic ciliopathies should also be considered in subjects with non-syndromic retinal dystrophy.

Investigation of DNA repair-related SNPs underlying susceptibility to papillary thyroid carcinoma reveals MGMT as a novel candidate gene in Belarusian children exposed to radiation.

BACKGROUND: Genetic factors may influence an individual's sensitivity to ionising radiation and therefore modify his/her risk of developing papillary thyroid carcinoma (PTC). Previously, we reported that common single nucleotide polymorphisms (SNPs) within the DNA damage recognition gene ATM contribute to PTC risk in Belarusian children exposed to fallout from the Chernobyl power plant accident. Here we explored in the same population the contribution of a panel of DNA repair-related SNPs in genes acting downstream of ATM. METHODS: The association of 141 SNPs located in 43 DNA repair genes was examined in 75 PTC cases and 254 controls from the Gomel region in Belarus. All subjects were younger than 15 years at the time of the Chernobyl accident. Conditional logistic regressions accounting for radiation dose were performed with PLINK using the additive allelic inheritance model, and a linkage disequilibrium (LD)-based Bonferroni correction was used for correction for multiple testing. RESULTS: The intronic SNP rs2296675 in MGMT was associated with an increased PTC risk [per minor allele odds ratio (OR) 2.54 95% CI 1.50, 4.30, P per allele = 0.0006, P corr.= 0.05], and gene-wide association testing highlighted a possible role for ERCC5 (P Gene = 0.01) and PCNA (P Gene = 0.05) in addition to MGMT (P Gene = 0.008). CONCLUSIONS: These findings indicate that several genes acting in distinct DNA repair mechanisms contribute to PTC risk. Further investigation is needed to decipher the functional properties of the methyltransferase encoded by MGMT and to understand how alteration of such functions may lead to the development of the most common type of thyroid cancer.

CRB1 mutations in inherited retinal dystrophies.

Mutations in the CRB1 gene are associated with variable phenotypes of severe retinal dystrophies, ranging from leber congenital amaurosis (LCA) to rod-cone dystrophy, also called retinitis pigmentosa (RP). Moreover, retinal dystrophies resulting from CRB1 mutations may be accompanied by specific fundus features: preservation of the para-arteriolar retinal pigment epithelium (PPRPE) and retinal telangiectasia with exudation (also referred to as Coats-like vasculopathy). In this publication, we report seven novel mutations and classify over 150 reported CRB1 sequence variants that were found in more that 240 patients. The data from previous reports were used to analyze a potential correlation between CRB1 variants and the clinical features of respective patients. This meta-analysis suggests that the differential phenotype of patients with CRB1 mutations is due to additional modifying factors rather than particular mutant allele combination.

A novel DFNB31 mutation associated with Usher type 2 syndrome showing variable degrees of auditory loss in a consanguineous Portuguese family.

PURPOSE: To identify the genetic defect of a consanguineous Portuguese family with rod-cone dystrophy and varying degrees of decreased audition. METHODS: A detailed ophthalmic and auditory examination was performed on a Portuguese patient with severe autosomal recessive rod-cone dystrophy. Known genetic defects were excluded by performing autosomal recessive retinitis pigmentosa (arRP) genotyping microarray analysis and by Sanger sequencing of the coding exons and flanking intronic regions of eyes shut homolog-drosophila (EYS) and chromosome 2 open reading frame 71 (C2orf71). Subsequently, genome-wide homozygosity mapping was performed in DNA samples from available family members using a 700K single nucleotide polymorphism (SNP) microarray. Candidate genes present in the significantly large homozygous regions were screened for mutations using Sanger sequencing. RESULTS: The largest homozygous region (~11 Mb) in the affected family members was mapped to chromosome 9, which harbors deafness, autosomal recessive 31 (DFNB31; a gene previously associated with Usher syndrome). Mutation analysis of DFNB31 in the index patient identified a novel one-base-pair deletion (c.737delC), which is predicted to lead to a truncated protein (p.Pro246HisfsX13) and co-segregated with the disease in the family. Ophthalmic examination of the index patient and the affected siblings showed severe rod-cone dystrophy. Pure tone audiometry revealed a moderate hearing loss in the index patient, whereas the affected siblings were reported with more profound and early onset hearing impairment. CONCLUSIONS: We report a novel truncating mutation in DFNB31 associated with severe rod-cone dystrophy and varying degrees of hearing impairment in a consanguineous family of Portuguese origin. This is the second report of DFNB31 implication in Usher type 2.

Gene network and biological pathways associated with susceptibility to differentiated thyroid carcinoma.

Variants identified in earlier genome-wide association studies (GWAS) on differentiated thyroid carcinoma (DTC) explain about 10% of the overall estimated genetic contribution and could not provide complete insights into biological mechanisms involved in DTC susceptibility. Integrating systems biology information from model organisms, genome-wide expression data from tumor and matched normal tissue and GWAS data could help identifying DTC-associated genes, and pathways or functional networks in which they are involved. We performed data mining of GWAS data of the EPITHYR consortium (1551 cases and 1957 controls) using various pathways and protein-protein interaction (PPI) annotation databases and gene expression data from The Cancer Genome Atlas. We identified eight DTC-associated genes at known loci 2q35 (DIRC3), 8p12 (NRG1), 9q22 (FOXE1, TRMO, HEMGN, ANP32B, NANS) and 14q13 (MBIP). Using the EW_dmGWAS approach we found that gene networks related to glycogenolysis, glycogen metabolism, insulin metabolism and signal transduction pathways associated with muscle contraction were overrepresented with association signals (false discovery rate adjusted p-value < 0.05). Additionally, suggestive association of 21 KEGG and 75 REACTOME pathways with DTC indicate a link between DTC susceptibility and functions related to metabolism of cholesterol, amino sugar and nucleotide sugar metabolism, steroid biosynthesis, and downregulation of ERBB2 signaling pathways. Together, our results provide novel insights into biological mechanisms contributing to DTC risk.

EIF2AK4 mutations cause pulmonary veno-occlusive disease, a recessive form of pulmonary hypertension.

Pulmonary veno-occlusive disease (PVOD) is a rare and devastating cause of pulmonary hypertension that is characterized histologically by widespread fibrous intimal proliferation of septal veins and preseptal venules and is frequently associated with pulmonary capillary dilatation and proliferation. PVOD is categorized into a separate pulmonary arterial hypertension-related group in the current classification of pulmonary hypertension. PVOD presents either sporadically or as familial cases with a seemingly recessive mode of transmission. Using whole-exome sequencing, we detected recessive mutations in EIF2AK4 (also called GCN2) that cosegregated with PVOD in all 13 families studied. We also found biallelic EIF2AK4 mutations in 5 of 20 histologically confirmed sporadic cases of PVOD. All mutations, either in a homozygous or compound-heterozygous state, disrupted the function of the gene. These findings point to EIF2AK4 as the major gene that is linked to PVOD development and contribute toward an understanding of the complex genetic architecture of pulmonary hypertension.

Positive selection of protective variants for type 2 diabetes from the Neolithic onward: a case study in Central Asia.

The high prevalence of type 2 diabetes and its uneven distribution among human populations is both a major public health concern and a puzzle in evolutionary biology. Why is this deleterious disease so common, while the associated genetic variants should be removed by natural selection? The 'thrifty genotype' hypothesis proposed that the causal genetic variants were advantageous and selected for during the majority of human evolution. It remains, however, unclear whether genetic data support this scenario. In this study, we characterized patterns of selection at 10 variants associated with type 2 diabetes, contrasting one herder and one farmer population from Central Asia. We aimed at identifying which alleles (risk or protective) are under selection, dating the timing of selective events, and investigating the effect of lifestyle on selective patterns. We did not find any evidence of selection on risk variants, as predicted by the thrifty genotype hypothesis. Instead, we identified clear signatures of selection on protective variants, in both populations, dating from the beginning of the Neolithic, which suggests that this major transition was accompanied by a selective advantage for non-thrifty variants. Combining our results with worldwide data further suggests that East Asia was particularly prone to such recent selection of protective haplotypes. As much effort has been devoted so far to searching for thrifty variants, we argue that more attention should be paid to the evolution of non-thrifty variants.

HLA-DRB1*01 associated with cutaneous hypersensitivity induced by nevirapine and efavirenz.

HLA typing, demographic and immunological risk factors for nevirapine and efavirenz reactions were studied in a French cohort of HIV patients. Cases with isolated rash were significantly associated with HLA-DRB101 allele. No liver toxicity was observed and no association was detected with the percentage of CD4 T-cells. This study suggests that HLA-DRB101 allele plays an important role in susceptibility to cutaneous reactions associated with nevirapine and efavirenz in HIV patients.

A European study of HLA-B in Stevens-Johnson syndrome and toxic epidermal necrolysis related to five high-risk drugs.

BACKGROUND: Stevens-Johnson syndrome (SJS) and its severe form, toxic epidermal necrolysis (TEN), are rare but life-threatening cutaneous adverse reactions to drugs, especially to allopurinol, carbamazepine, lamotrigine, phenobarbital, phenytoine, sulfamethoxazole, oxicam and nevirapine. Recently, a strong association between carbamazepine and allopurinol induced SJS or TEN has been described with respectively, HLA-B*1502 and HLA-B*5801 in a Han Chinese population from Taiwan and other Asian countries. OBJECTIVE: The objective is to further investigate the relationship between SJS/TEN and HLA-B in a large number of patients in a European population. METHODS: HLA-B genotyping was performed on 150 patients included in a European study (RegiSCAR) of SJS and TEN. We focused on patients related to 'high-risk' drugs including: 31 cases related to allopurinol, 28 to sulfamethoxazole, 19 to lamotrigine and 14 to oxicam. RESULTS: Sixty-one percent of 31 allopurinol-induced SJS/TEN patients carried the HLA-B*5801 allele and the figure was 55% for 27 patients of European ancestry [odds ratio=80 (34-187)], (P<10(-6)) as previously observed in Han Chinese. For other drugs, two rare alleles showed a weaker association with SJS/TEN in a limited number of patients: B*38 for sulfamethoxazole or lamotrigine-related patients, and B*73 for oxicam. CONCLUSION: At variance with prior results in Asia, this study shows that even when HLA-B alleles behave as strong risk factors, as for allopurinol, they are neither sufficient nor necessary to explain the disease. Further investigations are necessary to delineate the exact role of the HLA region in SJS/TEN, and to look for other associations in other regions of the genome.

Linkage disequilibrium in human populations.

Whereas the human linkage map appears on limited evidence to be constant over populations, maps of linkage disequilibrium (LD) vary among populations that differ in gene history. The greatest difference is between populations of sub-Saharan origin and populations remotely derived from Africa after a major bottleneck that reduced their heterozygosity and altered their Malecot parameters, increasing the intercept M that reflects association in founders and decreasing the exponential decline epsilon. Variation among populations within this ethnic dichotomy is much smaller. These observations validate use of a cosmopolitan LD map based on a sizeable sample representing a large population reliably typed for markers at high density. Then an LD map for a region or isolate within an ethnic group may be created by fitting the sample LD to the cosmopolitan map, estimating Malecot parameters simultaneously. The cosmopolitan map scaled by epsilon recovers 95% of the information that a local map at the same density gives and therefore more than the information in a low-resolution local map. Relative to a Eurasian cosmopolitan map the scaling factors are estimated to be 0.82 for isolates of European descent, 1.53 for Yorubans, and 1.74 for African Americans. These observations are consistent with a common bottleneck (perhaps but not necessarily speciation) approximately 173,500 years ago, if the bottleneck associated with migration out of Africa was 100,000 years ago. Eurasian populations (especially isolates with numerous cases) are efficient for genome scans, and populations of recent African origin (such as African Americans) are efficient for identification of causal polymorphisms within a candidate sequence.