Investigation of candidate gene copy number identifies FCGR3B as a potential biomarker for rheumatoid arthritis.

OBJECTIVES: Copy number variants (CNVs) could explain a part of the missing heritability in rheumatoid arthritis (RA). Our goal is to investigate the association of RA with CNVs of three functional candidate genes, Glutathione S-transferase M1 (GSTM1), Glutathione S-transferase T1 (GSTT1) and Fcγ receptor type IIIAB (FCGR3B). METHODS: We quantified the absolute copy number of GSTM1, GSTT1 and FCGR3B genes using droplet digital PCR. Transmission of copy number alleles was investigated in trio families with RA using family-based association tests (Transmission Disequilibrium Test and Genotype Haplotype Relative Risk). Clinical, environmental and biological data on RA patients were also used to stratify patients sample in analysis. RESULTS: Copy numbers from zero to three were identified. Genotype combinations characterised in 182 trios allowed testing the association with RA. Genotypes without null allele of FCGR3B gene were significantly associated with RA (3.41x10-7). Three copy numbers of this gene is observed only in cases of RA (n=14) and a protective effect of null allele was characterised (OR=0.3 (0.17-0.53)). CONCLUSIONS: CNVs in FCGR3B are associated with RA in our set of samples. This gene may play a role in physiopathology of this disease.

MYLK*FLNB and DOCK1*LAMA2 gene-gene interactions associated with rheumatoid arthritis in the focal adhesion pathway.

Rheumatoid arthritis (RA) is a chronic, systemic autoimmune disease caused by a combination of genetic and environmental factors. Rare variants with low predicted effects in genes participating in the same biological function might be involved in developing complex diseases such as RA. From whole-exome sequencing (WES) data, we identified genes containing rare non-neutral variants with complete penetrance and no phenocopy in at least one of nine French multiplex families. Further enrichment analysis highlighted focal adhesion as the most significant pathway. We then tested if interactions between the genes participating in this function would increase or decrease the risk of developing RA disease. The model-based multifactor dimensionality reduction (MB-MDR) approach was used to detect epistasis in a discovery sample (19 RA cases and 11 healthy individuals from 9 families and 98 unrelated CEU controls from the International Genome Sample Resource). We identified 9 significant interactions involving 11 genes (MYLK, FLNB, DOCK1, LAMA2, RELN, PIP5K1C, TNC, PRKCA, VEGFB, ITGB5, and FLT1). One interaction (MYLK*FLNB) increasing RA risk and one interaction decreasing RA risk (DOCK1*LAMA2) were confirmed in a replication sample (200 unrelated RA cases and 91 GBR unrelated controls). Functional and genomic data in RA samples or relevant cell types argue the key role of these genes in RA.

Familial melanoma: clinical factors associated with germline CDKN2A mutations according to the number of patients affected by melanoma in a family.

BACKGROUND: Features associated with an increased frequency of cyclin-dependent kinase inhibitor 2A (CDKN2A) mutations have been identified in families with 3 or more patients with cutaneous melanoma (CM). However, in families with 2 patients with CM, which represent the majority of familial melanoma, these factors have been rarely studied. OBJECTIVE: We investigated association of 3 clinical features with the presence of a CDKN2A mutation in a family by extent of CM family clustering (2 vs ≥3 patients with CM among first-degree relatives in a family). METHODS: We included 483 French families that comprised 387 families with 2 patients with CM (F2 families) and 96 families with 3 or more patients with CM (F3+ families). Three clinical factors were examined individually and in a joint analysis: median age at diagnosis younger than 50 years, and 1 or more patient in a family with multiple primary melanoma or with pancreatic cancer. RESULTS: The frequency of CDKN2A mutations was higher in F3+ families (32%) than in F2 families (13%). Although early age at melanoma diagnosis and occurrence of multiple primary melanoma in 1 or more patient were significantly associated with the risk of a CDKN2A mutation in F2 families, early age at melanoma diagnosis and occurrence of pancreatic cancer in a family were significantly associated with CDKN2A mutations in F3+ families. LIMITATIONS: The study was not population based. CONCLUSIONS: This study shows that factors associated with CDKN2A mutations differ by extent of CM family clustering. It indicates that, in France, families with 2 patients with CM are eligible for genetic testing especially when there is an early age at CM diagnosis and/or 1 or more patients with multiple primary melanoma.

Inference of an Integrative, Executable Network for Rheumatoid Arthritis Combining Data-Driven Machine Learning Approaches and a State-of-the-Art Mechanistic Disease Map.

Rheumatoid arthritis (RA) is a multifactorial, complex autoimmune disease that involves various genetic, environmental, and epigenetic factors. Systems biology approaches provide the means to study complex diseases by integrating different layers of biological information. Combining multiple data types can help compensate for missing or conflicting information and limit the possibility of false positives. In this work, we aim to unravel mechanisms governing the regulation of key transcription factors in RA and derive patient-specific models to gain more insights into the disease heterogeneity and the response to treatment. We first use publicly available transcriptomic datasets (peripheral blood) relative to RA and machine learning to create an RA-specific transcription factor (TF) co-regulatory network. The TF cooperativity network is subsequently enriched in signalling cascades and upstream regulators using a state-of-the-art, RA-specific molecular map. Then, the integrative network is used as a template to analyse patients' data regarding their response to anti-TNF treatment and identify master regulators and upstream cascades affected by the treatment. Finally, we use the Boolean formalism to simulate in silico subparts of the integrated network and identify combinations and conditions that can switch on or off the identified TFs, mimicking the effects of single and combined perturbations.

A novel nonsense variant in SUPT20H gene associated with Rheumatoid Arthritis identified by Whole Exome Sequencing of multiplex families.

The triggering and development of Rheumatoid Arthritis (RA) is conditioned by environmental and genetic factors. Despite the identification of more than one hundred genetic variants associated with the disease, not all the cases can be explained. Here, we performed Whole Exome Sequencing in 9 multiplex families (N = 30) to identify rare variants susceptible to play a role in the disease pathogenesis. We pre-selected 77 genes which carried rare variants with a complete segregation with RA in the studied families. Follow-up linkage and association analyses with pVAAST highlighted significant RA association of 43 genes (p-value < 0.05 after 106 permutations) and pinpointed their most likely causal variant. We re-sequenced the 10 most significant likely causal variants (p-value ≤ 3.78*10-3 after 106 permutations) in the extended pedigrees and 9 additional multiplex families (N = 110). Only one SNV in SUPT20H: c.73A>T (p.Lys25*), presented a complete segregation with RA in an extended pedigree with early-onset cases. In summary, we identified in this study a new variant associated with RA in SUPT20H gene. This gene belongs to several biological pathways like macro-autophagy and monocyte/macrophage differentiation, which contribute to RA pathogenesis. In addition, these results showed that analyzing rare variants using a family-based approach is a strategy that allows to identify RA risk loci, even with a small dataset.

Features associated with germline CDKN2A mutations: a GenoMEL study of melanoma-prone families from three continents.

BACKGROUND: The major factors individually reported to be associated with an increased frequency of CDKN2A mutations are increased number of patients with melanoma in a family, early age at melanoma diagnosis, and family members with multiple primary melanomas (MPM) or pancreatic cancer. METHODS: These four features were examined in 385 families with > or =3 patients with melanoma pooled by 17 GenoMEL groups, and these attributes were compared across continents. RESULTS: Overall, 39% of families had CDKN2A mutations ranging from 20% (32/162) in Australia to 45% (29/65) in North America to 57% (89/157) in Europe. All four features in each group, except pancreatic cancer in Australia (p = 0.38), individually showed significant associations with CDKN2A mutations, but the effects varied widely across continents. Multivariate examination also showed different predictors of mutation risk across continents. In Australian families, > or =2 patients with MPM, median age at melanoma diagnosis < or =40 years and > or =6 patients with melanoma in a family jointly predicted the mutation risk. In European families, all four factors concurrently predicted the risk, but with less stringent criteria than in Australia. In North American families, only > or =1 patient with MPM and age at diagnosis < or =40 years simultaneously predicted the mutation risk. CONCLUSIONS: The variation in CDKN2A mutations for the four features across continents is consistent with the lower melanoma incidence rates in Europe and higher rates of sporadic melanoma in Australia. The lack of a pancreatic cancer-CDKN2A mutation relationship in Australia probably reflects the divergent spectrum of mutations in families from Australia versus those from North America and Europe. GenoMEL is exploring candidate host, genetic and/or environmental risk factors to better understand the variation observed.

Characterization of copy number variants for CCL3L1 gene in rheumatoid arthritis for French trio families and Tunisian cases and controls.

Analyses of copy number variants (CNVs) for candidate genes in complex diseases are currently a promising research field. CNVs of C-C chemokine ligand 3-like 1 (CCL3L1) gene are candidate genomic factors in rheumatoid arthritis (RA). We investigated CCL3L1 CNVs association with a case-control study in Tunisians and a transmission analysis in French trio families. Relative copy number (rCN) of CCL3L1 gene was quantified by droplet digital PCR (ddPCR) in 100 French trio families (RA patients and their two parents) and in 166 RA cases and 102 healthy controls from Tunisia. We calculated odds ratio (OR) to investigate association risk for CCL3L1 CNVs in RA. rCN identified varied from 0 to 4 in the French population and from 0 to 7 in the Tunisian population. A significant difference was observed in the distribution of these rCNs between the two populations (p = 2.34 × 10(-10)), as when rCN from French and Tunisian RA patients were compared (p = 2.83 × 10(-5)). CNVs transmission in French RA trios allowed the characterization of genotypes with the presence of tandem duplication and triplication on the same chromosome. RA association tests highlighted a protective effect of rCN = 5 for CCL3L1 gene in the Tunisian population (OR = 0.056; CI 95 % [0.01-0.46]). Characterization of CCL3L1 CNVs with ddPCR methodology highlighted specific CN genotypes in a French family sample. A copy number polymorphism of a RA candidate gene was quantified, and its significant association with RA was revealed in a Tunisian sample.

Melanocortin-1 receptor (MC1R) gene variants and dysplastic nevi modify penetrance of CDKN2A mutations in French melanoma-prone pedigrees.

Germline mutations in CDKN2A gene predispose to melanoma with high but incomplete penetrance. Penetrance of CDKN2A gene was found to be significantly influenced by host factors (nevus phenotypes and sunburn) on one hand and by variants of MC1R gene (RHC variants consistently associated with red hair and fair skin) on the other hand. Our goal was to examine the joint effects of MC1R variants and other potential risk factors [total nevi, dysplastic nevi, pigmentary traits (skin, hair and eye color), skin reactions to sunlight, and degree of sun exposure] on CDKN2A penetrance. Clinical, genetic, and covariate data were recorded in 20 French melanoma-prone families with cosegregating CDKN2A mutations. Analysis of the cotransmission of melanoma and CDKN2A mutations was conducted by likelihood-based methods using the regressive logistic models, which can account for a variation of disease risk with age and can include the aforementioned risk factors as covariates. RHC variants, considered either alone or in the presence of pigmentation and nevus phenotypes, were found to increase significantly CDKN2A penetrance. Multivariate analysis, using a stepwise selection procedure, showed significant effects of two factors on melanoma risk in CDKN2A mutations carriers: RHC variants [odds ratio of hazard function (OR), 2.21; P = 0.03] and dysplastic nevi (OR, 2.93; P < 0.01). Such results may have important consequences to improve the prediction of melanoma risk in families.

Familial correlations and inter-relationships of four asthma-associated quantitative phenotypes in 320 French EGEA families ascertained through asthmatic probands.

Asthma is a complex disease, associated with biological and physiological phenotypes including immunoglobulin E (IgE) levels, sum of positive skin prick tests to allergens (SPTQ), eosinophil counts (EOS) and percent predicted forced expiratory volume in 1 s (%FEV1). We investigated the patterns of familial correlations and the inter-relationships of these four quantitative phenotypes, using the general class D regressive model, in 320 French EGEA nuclear families ascertained through 204 offspring (set A) and 116 parents (set B). Familial correlations of IgE and SPTQ were consistent with a model including no spouse correlation and equal parent-offspring and sib-sib correlations (rhoPO = rhoSS = 0.25 for IgE and 0.15 for SPTQ), this model being compatible with an additive polygenic model in the whole sample and the two family subsets A and B. Different patterns of familial correlations of EOS and %FEV1 were observed in these two sets. In set A, the best fitting model included no spouse correlation and equality of parent-offspring and sib-sib correlations (rhoPO = rhoSS = 0.14 for EOS and 0.23 for %FEV1). In set B, EOS had only a significant rhoSS of 0.28, while %FEV1 had significant rhoMO of 0.28 and rhoSS of 0.16. Analysis of shared familial determinants between these phenotypes indicated an overlap of at most 30% in rhoFO for IgE and SPTQ and in both rhoFO and rhoMO for IgE and EOS, while determinants of %FEV1 and atopy-related phenotypes appear distinct. These results may have implications for further linkage and association studies with genetic markers.

Association of genetic variants in CDK6 and XRCC1 with the risk of dysplastic nevi in melanoma-prone families.

Dysplastic nevi (DN) is a strong risk factor for cutaneous malignant melanoma (CMM), and it frequently occurs in melanoma-prone families. To identify genetic variants for DN, we genotyped 677 tagSNPs in 38 melanoma candidate genes that are involved in pigmentation, DNA repair, cell cycle control, and melanocyte proliferation pathways in a total of 504 individuals (310 with DN, 194 without DN) from 53 melanoma-prone families (23 CDKN2A mutation positive and 30 negative). Conditional logistic regression, conditioning on families, was used to estimate the association between DN and each single-nucleotide polymorphism (SNP) separately, adjusted for age, sex, CMM, and CDKN2A status. P-values for SNPs in the same gene were combined to yield gene-specific P-values. Two genes, CDK6 (cyclin-dependent kinase 6) and XRCC1, were significantly associated with DN after Bonferroni correction for multiple testing (P=0.0001 and 0.00025, respectively), whereas neither gene was significantly associated with CMM. Associations for CDK6 SNPs were stronger in CDKN2A mutation-positive families (rs2079147, Pinteraction=0.0033), whereas XRCC1 SNPs had similar effects in mutation-positive and -negative families. The association for one of the associated SNPs in XRCC1 (rs25487) was replicated in two independent data sets (random-effect meta-analysis: P<0.0001). Our findings suggest that some genetic variants may contribute to DN risk independently of their association with CMM in melanoma-prone families.

A variant in FTO shows association with melanoma risk not due to BMI.

We report the results of an association study of melanoma that is based on the genome-wide imputation of the genotypes of 1,353 cases and 3,566 controls of European origin conducted by the GenoMEL consortium. This revealed an association between several SNPs in intron 8 of the FTO gene, including rs16953002, which replicated using 12,313 cases and 55,667 controls of European ancestry from Europe, the USA and Australia (combined P = 3.6 × 10(-12), per-allele odds ratio for allele A = 1.16). In addition to identifying a new melanoma-susceptibility locus, this is to our knowledge the first study to identify and replicate an association with SNPs in FTO not related to body mass index (BMI). These SNPs are not in intron 1 (the BMI-related region) and exhibit no association with BMI. This suggests FTO's function may be broader than the existing paradigm that FTO variants influence multiple traits only through their associations with BMI and obesity.

Characteristics of the coexistence of melanoma and renal cell carcinoma.

BACKGROUND: Patients with melanoma (MM) have an increased risk of kidney cancer, and there is an excess risk of MM among patients with renal cell carcinoma (RCC). The objective of the current study was to analyze a series of 42 patients with both MM and RCC to identify clinical and pathologic features as well as risk factors of this association. METHODS: Clinical and pathologic characteristics of 42 patients who developed both MM and RCC (the MM + RCC series) were compared with 2 published series in each cancer alone: a series of 293 patients with MM (MM series) and a series of 1527 patients with RCC (RCC series). RESULTS: RCC was diagnosed concomitantly or after MM in 83% of patients in the MM + RCC series. Those patients displayed a high proportion of asymptomatic RCC at diagnosis (70%) and a higher frequency of stage I tumors (61%) than patients in the RCC series. Compared with the MM series, patients in the MM + RCC series more often were men, had a higher frequency of blond/red hair, had poor tanning ability, and had a higher number of nevi. In addition, patients in the MM + RCC series had a high aggregation of other malignancies (mainly skin cancers) and a significantly higher frequency of family history of MM (P = .005). Only 2 cyclin-dependent kinase 2A gene (CDKN2A) germline mutations were identified among patients in the MM + RCC series who also were members of MM-prone families. CONCLUSIONS: The high aggregation of cancers among patients in the MM + RCC series and the familial clustering of MM argued for a genetic predisposition that may be partly independent of CDKN2A.

MC1R variant alleles and malignant melanoma risk in Israel.

To evaluate the contribution of MC1R variants to malignant melanoma risk in Israeli Jews, sequencing of the MC1R gene was performed in 132 melanoma patients and 184 ethnically matched controls. Overall, 22 MC1R variants were detected, two were novel (M73I and 496_497insG). Using age and sex-adjusted logistic regression, one specific variant, R151C, conferred significantly increased melanoma risk among Ashkenazim (OR=2.6, 95% CI: 1.3-5.3; p=0.05 after Bonferroni correction). A gene dosage effect was noted, with significantly increased melanoma risk being observed in subjects with at least two variants whether when all variants are pooled (OR=4.8, 95% CI: 2.0-11.2; p=0.002 after Bonferroni correction) or when red hair colour (RHC) variants and non-RHC variants are distinguished (OR=7.6, 95% CI: 2.8-20.3; p=0.0004 after Bonferroni correction). If further studies support these findings, the assessment of MC1R status may be useful in identifying Jewish Israeli individuals at high risk for melanoma.

Protective effect of copy number polymorphism of glutathione S-transferase T1 gene on melanoma risk in presence of CDKN2A mutations, MC1R variants and host-related phenotypes.

The effect of CDKN2A, the major high-risk melanoma susceptibility gene, has been shown to be modified by host-related phenotypes and variants of MC1R gene. The glutathione S-transferase (GSTs) genes, implicated in detoxification of metabolites after UV exposure, are candidates for modulating CDKN2A penetrance. Few case-control studies have investigated the effect of GSTs on melanoma risk, and have led to controversial results while these genes have not yet been studied in CDKN2A melanoma-prone families. We examined the effect of GSTP1, GSTM1 and GSTT1 genotypes on melanoma risk in 25 multi-generational melanoma-prone families with CDKN2A mutations, in presence of MC1R gene variants, sun exposure, and host-related phenotypes. These data included 195 genotyped subjects for all studied genes. We applied the GEE (Generalized Estimating Equations) approach to test for the effect of GSTs while adjusting for age, sex and CDKN2A mutation status and including successively MC1R, sun exposure and host factors in the model. No significant effect of null GSTM1 allele and GSTP1 variants (p.I105V, p.A114V) on melanoma risk was found. However, a significant protective effect of carrying >or=1 null GSTT1 allele was shown: OR(adjusted for age,sex,CDKN2A ) = 0.41 (0.18-0.94) and OR(adjusted for age,sex,CDKN2A,MC1R ) = 0.24 (0.15-0.58). Altogether, the factors modifying significantly the melanoma risk associated with CDKN2A mutations (stepwise procedure) were: MC1R and dysplastic nevi (increasing the risk) and GSTT1 (decreasing the risk). This study shows that even when a high-risk gene (CDKN2A) has been identified, multiple genetic modifiers influence melanoma risk.

Cutaneous phenotype and MC1R variants as modifying factors for the development of melanoma in CDKN2A G101W mutation carriers from 4 countries.

The G101W founder mutation is the most common CDKN2A mutation in Italy, Spain, and France. As the background of modifying genes, environmental exposures, and sun behavior vary across countries, studying G101W carriers from distinct countries offers a unique opportunity to evaluate possible modifying factors in melanoma development. We evaluated 76 G101W cases and 59 carrier controls from France, Italy, Spain, and the United States. Hair color and dysplastic nevi distributions differed significantly in cases and controls across the 4 study groups. Cases also varied significantly for eye color, freckling, and nevi. The distribution of MC1R variants in cases differed significantly across study groups because 12% of Italian melanoma patients had > or =2 MC1R variants vs. >50% for the other case groups. Several MC1R covariates showed significant associations with melanoma risk in all groups combined and in the American, French, and Spanish samples; no significant findings were observed in the Italian sample. In multiple-case families, the number and type of MC1R variants varied significantly between multiple-primary-melanoma and single-primary-melanoma patients from the 4 groups; there was also a significant decrease in median age at melanoma diagnosis as the number or type of MC1R variants increased. The variation in the effects of the cutaneous phenotypic and MC1R factors across the study sample suggests that these factors differentially contribute to development of melanoma even on a common genetic background of a germline CDKN2A mutation. Differences in melanoma risk across geographic regions justify the need for individual studies in each country before counseling should be considered.

Influence of genes, nevi, and sun sensitivity on melanoma risk in a family sample unselected by family history and in melanoma-prone families.

BACKGROUND: Few family studies have investigated the effects of genetic, environmental, and host factors on melanoma risk, and most have been restricted to high-risk families. We assessed the role of these factors on melanoma risk in two types of families: families ascertained through melanoma probands but unselected by family history and melanoma-prone families. METHODS: Data on pigmentary traits, nevus phenotypes, exposure to sun, and reactions to sunlight were collected from 295 families unselected by family history and 53 melanoma-prone families. We modeled melanoma risk using a logistic regressive model incorporating the effect of a melanoma-predisposing gene, familial dependence, and potential risk factors (e.g., pigmentary traits, nevus phenotypes, history of sun exposure, skin reactions to sunlight). Maximum-likelihood estimates of the parameters of the regressive model were used to compute odds ratios associated with each risk factor and age-specific melanoma risk depending on the genotype at the melanoma-predisposing gene and the effects of risk factors. All statistical tests were two-sided. RESULTS: In the families unselected by family history, there was statistically significant evidence (P<.001) for a dominant gene, with melanoma risk reaching 0.49 and 0.67 by age 80 years in male and female gene carriers, respectively. Melanoma risk was statistically significantly influenced by total nevi (odds ratio of hazard function [OR] = 5.81, 95% confidence interval [CI] = 3.47 to 8.99), sun exposure (OR = 5.37, 95% CI = 4.44 to 6.36), and sunburn interacting with the gene (OR = 26.31, 95% CI = 7.56 to 99.22 in gene carriers and OR = 1.67, 95% CI = 1.36 to 2.03 in noncarriers). Twenty of the 53 melanoma-prone families had cosegregating mutations in CDKN2A, a gene known to be associated with melanoma. In these 53 families, three risk factors in addition to CDKN2A mutations increased melanoma risk: dysplastic nevi (OR = 2.32, 95% CI = 2.08 to 2.58), total nevi (OR = 1.99, 95% CI = 1.61 to 2.20) and sunburn (OR = 5.16, 95% CI = 4.82 to 5.52). CONCLUSIONS: Together, a melanoma-predisposing gene (identified as being CDKN2A in melanoma-prone families), number of nevi and/or dysplastic nevi, and sun-related covariates influence melanoma risk in both families unselected by family history and melanoma-prone families.

Interactions between genetic and reproductive factors in breast cancer risk in a population-based sample of African-American families.

Incidence of breast cancer (BC) varies among ethnic groups, with higher rates in white than in African-American women. Until now, most epidemiological and genetic studies have been carried out in white women. To investigate whether interactions between genetic and reproductive risk factors may explain part of the ethnic disparity in BC incidence, a genetic epidemiology study was conducted, between 1989 and 1994, at the Howard University Cancer Center (Washington, DC), which led to the recruitment of 245 African-American families. Segregation analysis of BC was performed by use of the class D regressive logistic model that allows for censored data to account for a variable age of onset of disease, as implemented in the REGRESS program. Segregation analysis of BC was consistent with a putative dominant gene effect (P < 0.000001) and residual sister-dependence (P < 0.0001). This putative gene was found to interact significantly with age at menarche (P = 0.048), and an interaction with a history of spontaneous abortions was suggested (P = 0.08). A late age at menarche increased BC risk in gene carriers but had a protective effect in non-gene carriers. A history of spontaneous abortions had a protective effect in gene carriers and increased BC risk in non-gene carriers. Our findings agree partially with a similar analysis of French families showing a significant gene x parity interaction and a suggestive gene x age at menarche interaction. Investigating gene x risk factor interactions in different populations may have important implications for further biological investigations and for BC risk assessment.