Genes with differential expression across ancestries are enriched in ancestry-specific disease effects likely due to gene-by-environment interactions.

Multi-ancestry genome-wide association studies (GWASs) have highlighted the existence of variants with ancestry-specific effect sizes. Understanding where and why these ancestry-specific effects occur is fundamental to understanding the genetic basis of human diseases and complex traits. Here, we characterized genes differentially expressed across ancestries (ancDE genes) at the cell-type level by leveraging single-cell RNA-sequencing data in peripheral blood mononuclear cells for 21 individuals with East Asian (EAS) ancestry and 23 individuals with European (EUR) ancestry (172,385 cells); then, we tested whether variants surrounding those genes were enriched in disease variants with ancestry-specific effect sizes by leveraging ancestry-matched GWASs of 31 diseases and complex traits (average n ∼ 90,000 and ∼ 267,000 in EAS and EUR, respectively). We observed that ancDE genes tended to be cell-type specific and enriched in genes interacting with the environment and in variants with ancestry-specific disease effect sizes, which suggests cell-type-specific, gene-by-environment interactions shared between regulatory and disease architectures. Finally, we illustrated how different environments might have led to ancestry-specific myeloid cell leukemia 1 (MCL1) expression in B cells and ancestry-specific allele effect sizes in lymphocyte count GWASs for variants surrounding MCL1. Our results imply that large single-cell and GWAS datasets from diverse ancestries are required to improve our understanding of human diseases.

RPA3-UMAD1 rs12702634 and rheumatoid arthritis-associated interstitial lung disease in European ancestry.

Objective: Recently, a genome-wide association study identified an association between RA-associated interstitial lung disease (ILD) and RPA3-UMAD1 rs12702634 in the Japanese population, especially for patients with a usual interstitial pneumonia (UIP) pattern. We aimed to replicate this association in a European population and test for interaction with MUC5B rs35705950. Methods: In this genetic case-control association study, patients with RA and ILD and controls with RA and no ILD were included from France, the USA and the Netherlands. Only cases and controls from European genetic ancestries determined by principal components analysis were included in the analyses. RA was defined by the 1987 ACR or 2010 ACR/EULAR criteria and ILD by chest high-resolution CT scan, except in the control dataset from the Netherlands, where the absence of ILD was determined by chart review. Patients were genotyped for RPA3-UMAD1 rs12702634 and MUC5B rs35705950. Associations were tested using logistic regression adjusted for sex, age at RA onset, age at ILD onset or at certified absence of ILD, tobacco smoking status and country of origin. Results: Among the 883 patients included, 322 were RA-ILD cases (36.5%). MUC5B rs35705950 was strongly associated with RA-ILD in all datasets {combined adjusted odds ratio [OR] 2.9 [95% CI 2.1, 3.9], P = 1.1 × 10-11. No association between RPA3-UMAD1 rs12702634 and RA-ILD was observed [combined OR 1.2 (95% CI 0.8, 1.6), P = 0.31. No interaction was found between RPA3-UMAD1 rs12702634 and MUC5B rs35705950 (P = 0.70). Conclusion: Our findings did not support a contribution of RPA3-UMAD1 rs12702634 to the overall RA-ILD susceptibility in the European population.

A noncoding regulatory variant in IKZF1 increases acute lymphoblastic leukemia risk in Hispanic/Latino children.

Hispanic/Latino children have the highest risk of acute lymphoblastic leukemia (ALL) in the US compared to other racial/ethnic groups, yet the basis of this remains incompletely understood. Through genetic fine-mapping analyses, we identified a new independent childhood ALL risk signal near IKZF1 in self-reported Hispanic/Latino individuals, but not in non-Hispanic White individuals, with an effect size of ∼1.44 (95% confidence interval = 1.33-1.55) and a risk allele frequency of ∼18% in Hispanic/Latino populations and <0.5% in European populations. This risk allele was positively associated with Indigenous American ancestry, showed evidence of selection in human history, and was associated with reduced IKZF1 expression. We identified a putative causal variant in a downstream enhancer that is most active in pro-B cells and interacts with the IKZF1 promoter. This variant disrupts IKZF1 autoregulation at this enhancer and results in reduced enhancer activity in B cell progenitors. Our study reveals a genetic basis for the increased ALL risk in Hispanic/Latino children.

A scalable approach to characterize pleiotropy across thousands of human diseases and complex traits using GWAS summary statistics.

Genome-wide association studies (GWASs) across thousands of traits have revealed the pervasive pleiotropy of trait-associated genetic variants. While methods have been proposed to characterize pleiotropic components across groups of phenotypes, scaling these approaches to ultra-large-scale biobanks has been challenging. Here, we propose FactorGo, a scalable variational factor analysis model to identify and characterize pleiotropic components using biobank GWAS summary data. In extensive simulations, we observe that FactorGo outperforms the state-of-the-art (model-free) approach tSVD in capturing latent pleiotropic factors across phenotypes while maintaining a similar computational cost. We apply FactorGo to estimate 100 latent pleiotropic factors from GWAS summary data of 2,483 phenotypes measured in European-ancestry Pan-UK BioBank individuals (N = 420,531). Next, we find that factors from FactorGo are more enriched with relevant tissue-specific annotations than those identified by tSVD (p = 2.58E-10) and validate our approach by recapitulating brain-specific enrichment for BMI and the height-related connection between reproductive system and muscular-skeletal growth. Finally, our analyses suggest shared etiologies between rheumatoid arthritis and periodontal condition in addition to alkaline phosphatase as a candidate prognostic biomarker for prostate cancer. Overall, FactorGo improves our biological understanding of shared etiologies across thousands of GWASs.

A scalable variational approach to characterize pleiotropic components across thousands of human diseases and complex traits using GWAS summary statistics.

Genome-wide association studies (GWAS) across thousands of traits have revealed the pervasive pleiotropy of trait-associated genetic variants. While methods have been proposed to characterize pleiotropic components across groups of phenotypes, scaling these approaches to ultra large-scale biobanks has been challenging. Here, we propose FactorGo, a scalable variational factor analysis model to identify and characterize pleiotropic components using biobank GWAS summary data. In extensive simulations, we observe that FactorGo outperforms the state-of-the-art (model-free) approach tSVD in capturing latent pleiotropic factors across phenotypes, while maintaining a similar computational cost. We apply FactorGo to estimate 100 latent pleiotropic factors from GWAS summary data of 2,483 phenotypes measured in European-ancestry Pan-UK BioBank individuals (N=420,531). Next, we find that factors from FactorGo are more enriched with relevant tissue-specific annotations than those identified by tSVD (P=2.58E-10), and validate our approach by recapitulating brain-specific enrichment for BMI and the height-related connection between reproductive system and muscular-skeletal growth. Finally, our analyses suggest novel shared etiologies between rheumatoid arthritis and periodontal condition, in addition to alkaline phosphatase as a candidate prognostic biomarker for prostate cancer. Overall, FactorGo improves our biological understanding of shared etiologies across thousands of GWAS.

Leveraging base-pair mammalian constraint to understand genetic variation and human disease.

Thousands of genomic regions have been associated with heritable human diseases, but attempts to elucidate biological mechanisms are impeded by an inability to discern which genomic positions are functionally important. Evolutionary constraint is a powerful predictor of function, agnostic to cell type or disease mechanism. Single-base phyloP scores from 240 mammals identified 3.3% of the human genome as significantly constrained and likely functional. We compared phyloP scores to genome annotation, association studies, copy-number variation, clinical genetics findings, and cancer data. Constrained positions are enriched for variants that explain common disease heritability more than other functional annotations. Our results improve variant annotation but also highlight that the regulatory landscape of the human genome still needs to be further explored and linked to disease.

Leveraging Base Pair Mammalian Constraint to Understand Genetic Variation and Human Disease.

Although thousands of genomic regions have been associated with heritable human diseases, attempts to elucidate biological mechanisms are impeded by a general inability to discern which genomic positions are functionally important. Evolutionary constraint is a powerful predictor of function that is agnostic to cell type or disease mechanism. Here, single base phyloP scores from the whole genome alignment of 240 placental mammals identified 3.5% of the human genome as significantly constrained, and likely functional. We compared these scores to large-scale genome annotation, genome-wide association studies (GWAS), copy number variation, clinical genetics findings, and cancer data sets. Evolutionarily constrained positions are enriched for variants explaining common disease heritability (more than any other functional annotation). Our results improve variant annotation but also highlight that the regulatory landscape of the human genome still needs to be further explored and linked to disease.

Leveraging fine-mapping and multipopulation training data to improve cross-population polygenic risk scores.

Polygenic risk scores suffer reduced accuracy in non-European populations, exacerbating health disparities. We propose PolyPred, a method that improves cross-population polygenic risk scores by combining two predictors: a new predictor that leverages functionally informed fine-mapping to estimate causal effects (instead of tagging effects), addressing linkage disequilibrium differences, and BOLT-LMM, a published predictor. When a large training sample is available in the non-European target population, we propose PolyPred+, which further incorporates the non-European training data. We applied PolyPred to 49 diseases/traits in four UK Biobank populations using UK Biobank British training data, and observed relative improvements versus BOLT-LMM ranging from +7% in south Asians to +32% in Africans, consistent with simulations. We applied PolyPred+ to 23 diseases/traits in UK Biobank east Asians using both UK Biobank British and Biobank Japan training data, and observed improvements of +24% versus BOLT-LMM and +12% versus PolyPred. Summary statistics-based analogs of PolyPred and PolyPred+ attained similar improvements.

COVID-19 tissue atlases reveal SARS-CoV-2 pathology and cellular targets.

COVID-19, which is caused by SARS-CoV-2, can result in acute respiratory distress syndrome and multiple organ failure1-4, but little is known about its pathophysiology. Here we generated single-cell atlases of 24 lung, 16 kidney, 16 liver and 19 heart autopsy tissue samples and spatial atlases of 14 lung samples from donors who died of COVID-19. Integrated computational analysis uncovered substantial remodelling in the lung epithelial, immune and stromal compartments, with evidence of multiple paths of failed tissue regeneration, including defective alveolar type 2 differentiation and expansion of fibroblasts and putative TP63+ intrapulmonary basal-like progenitor cells. Viral RNAs were enriched in mononuclear phagocytic and endothelial lung cells, which induced specific host programs. Spatial analysis in lung distinguished inflammatory host responses in lung regions with and without viral RNA. Analysis of the other tissue atlases showed transcriptional alterations in multiple cell types in heart tissue from donors with COVID-19, and mapped cell types and genes implicated with disease severity based on COVID-19 genome-wide association studies. Our foundational dataset elucidates the biological effect of severe SARS-CoV-2 infection across the body, a key step towards new treatments.

Large-scale genome-wide association study in a Japanese population identifies novel susceptibility loci across different diseases.

The overwhelming majority of participants in current genetic studies are of European ancestry. To elucidate disease biology in the East Asian population, we conducted a genome-wide association study (GWAS) with 212,453 Japanese individuals across 42 diseases. We detected 320 independent signals in 276 loci for 27 diseases, with 25 novel loci (P < 9.58 × 10-9). East Asian-specific missense variants were identified as candidate causal variants for three novel loci, and we successfully replicated two of them by analyzing independent Japanese cohorts; p.R220W of ATG16L2 (associated with coronary artery disease) and p.V326A of POT1 (associated with lung cancer). We further investigated enrichment of heritability within 2,868 annotations of genome-wide transcription factor occupancy, and identified 378 significant enrichments across nine diseases (false discovery rate < 0.05) (for example, NKX3-1 for prostate cancer). This large-scale GWAS in a Japanese population provides insights into the etiology of complex diseases and highlights the importance of performing GWAS in non-European populations.

MUC5B Promoter Variant and Rheumatoid Arthritis with Interstitial Lung Disease.

BACKGROUND: Given the phenotypic similarities between rheumatoid arthritis (RA)-associated interstitial lung disease (ILD) (hereafter, RA-ILD) and idiopathic pulmonary fibrosis, we hypothesized that the strongest risk factor for the development of idiopathic pulmonary fibrosis, the gain-of-function MUC5B promoter variant rs35705950, would also contribute to the risk of ILD among patients with RA. METHODS: Using a discovery population and multiple validation populations, we tested the association of the MUC5B promoter variant rs35705950 in 620 patients with RA-ILD, 614 patients with RA without ILD, and 5448 unaffected controls. RESULTS: Analysis of the discovery population revealed an association of the minor allele of the MUC5B promoter variant with RA-ILD when patients with RA-ILD were compared with unaffected controls (adjusted odds ratio, 3.8; 95% confidence interval [CI], 2.8 to 5.2; P=9.7×10-17). The MUC5B promoter variant was also significantly overrepresented among patients with RA-ILD, as compared with unaffected controls, in an analysis of the multiethnic case series (adjusted odds ratio, 5.5; 95% CI, 4.2 to 7.3; P=4.7×10-35) and in a combined analysis of the discovery population and the multiethnic case series (adjusted odds ratio, 4.7; 95% CI, 3.9 to 5.8; P=1.3×10-49). In addition, the MUC5B promoter variant was associated with an increased risk of ILD among patients with RA (adjusted odds ratio in combined analysis, 3.1; 95% CI, 1.8 to 5.4; P=7.4×10-5), particularly among those with evidence of usual interstitial pneumonia on high-resolution computed tomography (adjusted odds ratio in combined analysis, 6.1; 95% CI, 2.9 to 13.1; P=2.5×10-6). However, no significant association with the MUC5B promoter variant was observed for the diagnosis of RA alone. CONCLUSIONS: We found that the MUC5B promoter variant was associated with RA-ILD and more specifically associated with evidence of usual interstitial pneumonia on imaging. (Funded by Société Française de Rhumatologie and others.).

Association of modifiers and other genetic factors explain Marfan syndrome clinical variability.

Marfan syndrome (MFS) is a rare autosomal dominant connective tissue disorder related to variants in the FBN1 gene. Prognosis is related to aortic risk of dissection following aneurysm. MFS clinical variability is notable, for age of onset as well as severity and number of clinical manifestations. To identify genetic modifiers, we combined genome-wide approaches in 1070 clinically well-characterized FBN1 disease-causing variant carriers: (1) an FBN1 eQTL analysis in 80 fibroblasts of FBN1 stop variant carriers, (2) a linkage analysis, (3) a kinship matrix association study in 14 clinically concordant and discordant sib-pairs, (4) a genome-wide association study and (5) a whole exome sequencing in 98 extreme phenotype samples.Three genetic mechanisms of variability were found. A new genotype/phenotype correlation with an excess of loss-of-cysteine variants (P = 0.004) in severely affected subjects. A second pathogenic event in another thoracic aortic aneurysm gene or the COL4A1 gene (known to be involved in cerebral aneurysm) was found in nine individuals. A polygenic model involving at least nine modifier loci (named gMod-M1-9) was observed through cross-mapping of results. Notably, gMod-M2 which co-localizes with PRKG1, in which activating variants have already been described in thoracic aortic aneurysm, and gMod-M3 co-localized with a metalloprotease (proteins of extra-cellular matrix regulation) cluster. Our results represent a major advance in understanding the complex genetic architecture of MFS and provide the first steps toward prediction of clinical evolution.

Truncating mutations of TP53AIP1 gene predispose to cutaneous melanoma.

Genetic predisposition to cutaneous malignant melanoma (CMM) involves highly penetrant predisposing genes and low and intermediate penetrant predisposing alleles. However, the missing heritability in (CMM) is still high. For such and in order to identify new genetic factors for CMM, we conducted an exome sequencing study in high-risk CMM patients. Two rounds of exome sequencing were successively performed in 33 and 27 high-risk patients. We focused on genes carrying rare nonsense, frameshift, and splice variants (allelic frequency <1%) that were present in both series of exomes. An extension study was then conducted in a large cohort (1 079 CMM patients and 1 230 Caucasian ethnically matched healthy controls), and the inactivating variants frequency was compared between groups using two-sided Fisher exact test. Two TP53AIP1 truncating mutations were identified in four patients: a frameshift c.63_64insG, p.Q22Afs*81 in two patients from the same family and in the proband of a second family; and a nonsense mutation c.95 C > A, p.Ser32Stop in a patient with multiple CMMs. In all patients, TP53AIP1 truncating variants were strongly associated with CMM risk (two-sided Fisher exact test = 0.004, OR = 3.3[1.3-8.5]). Additionally, we showed that TP53AIP1 mRNA was strongly down-regulated throughout different phases of melanoma progression. TP53AIP1 gene is a TP53 target which plays a key role by inducting apoptosis in response to UV-induced DNA damage. Constitutional mutations of TP53AIP1 had previously been involved in susceptibility to prostate cancer. Our results show that constitutional truncating TP53AIP1 mutations predispose to CMM in the French population. Replication studies in other populations should be performed.

Juvenile myelomonocytic leukemia displays mutations in components of the RAS pathway and the PRC2 network.

Juvenile myelomonocytic leukemia (JMML) is a rare and severe myelodysplastic and myeloproliferative neoplasm of early childhood initiated by germline or somatic RAS-activating mutations. Genetic profiling and whole-exome sequencing of a large JMML cohort (118 and 30 cases, respectively) uncovered additional genetic abnormalities in 56 cases (47%). Somatic events were rare (0.38 events/Mb/case) and restricted to sporadic (49/78; 63%) or neurofibromatosis type 1 (NF1)-associated (8/8; 100%) JMML cases. Multiple concomitant genetic hits targeting the RAS pathway were identified in 13 of 78 cases (17%), disproving the concept of mutually exclusive RAS pathway mutations and defining new pathways activated in JMML involving phosphoinositide 3-kinase (PI3K) and the mTORC2 complex through RAC2 mutation. Furthermore, this study highlights PRC2 loss (26/78; 33% of sporadic JMML cases) that switches the methylation/acetylation status of lysine 27 of histone H3 in JMML cases with altered RAS and PRC2 pathways. Finally, the association between JMML outcome and mutational profile suggests a dose-dependent effect for RAS pathway activation, distinguishing very aggressive JMML rapidly progressing to acute myeloid leukemia.

Heterozygous RTEL1 mutations are associated with familial pulmonary fibrosis.

Pulmonary fibrosis is a fatal disease with progressive loss of respiratory function. Defective telomere maintenance leading to telomere shortening is a cause of pulmonary fibrosis, as mutations in the telomerase component genes TERT (reverse transcriptase) and TERC (RNA component) are found in 15% of familial pulmonary fibrosis (FPF) cases. However, so far, about 85% of FPF remain genetically uncharacterised.Here, in order to identify new genetic causes of FPF, we performed whole-exome sequencing, with a candidate-gene approach, of 47 affected subjects from 35 families with FPF without TERT and TERC mutations.We identified heterozygous mutations in regulator of telomere elongation helicase 1 (RTEL1) in four families. RTEL1 is a DNA helicase with roles in DNA replication, genome stability, DNA repair and telomere maintenance. The heterozygous RTEL1 mutations segregated as an autosomal dominant trait in FPF, and were predicted by structural analyses to severely affect the function and/or stability of RTEL1. In agreement with this, RTEL1-mutated patients exhibited short telomeres in comparison with age-matched controls.Our results provide evidence that heterozygous RTEL1 mutations are responsible for FPF and, thereby, extend the clinical spectrum of RTEL1 deficiency. Thus, RTEL1 enlarges the number of telomere-associated genes implicated in FPF.

Genes involved in the WNT and vesicular trafficking pathways are associated with melanoma predisposition.

Multifactorial predisposition to melanoma includes genes involved in pigmentation, immunity and DNA repair. Nonetheless, missing heritability in melanoma is still important. We studied the role of 335 candidate SNPs in melanoma susceptibility by using a dedicated chip and investigating 110 genes involved in different pathways. A discovery set was comprised of 1069 melanoma patients and 925 controls from France. Data were replicated using validation phases II (1085 cases and 801 controls from Spain) and III (1808 cases and 1894 controls from Germany and a second set of Spanish samples). In addition, an exome sequencing study was performed in three high-risk French melanoma families. Nineteen SNPs in 17 genes were initially associated with melanoma in the French population. Six SNPs were replicated in phase II, including two new SNPs in the WNT3 (rs199524) and VPS41 (rs11773094) genes. The role of VPS41 and WNT3 was confirmed in a meta-analysis (3940 melanoma cases and 3620 controls) with two-side p values of 0.002, (OR = 0.86) and 4.07 × 10(-10) (OR = 0.80), respectively. Exome sequencing revealed a non-synonymous VPS41 variant in one family that was shown to be strongly associated with familial melanoma (OR = 4.46, p = 0.001) in an independent sample of 178 melanoma families. WNT3 belongs to WNT pathway known to play a crucial role in melanoma, whereas VPS41 regulates vesicular trafficking and is thought to play a role in pigmentation. Our work identified two new pathways involved in melanoma predisposition. These results may be useful in the future for identifying individuals highly predisposed to melanoma.

Identification of secreted phosphoprotein 1 gene as a new rheumatoid arthritis susceptibility gene.

OBJECTIVE: To evaluate the contribution of the SPP1 rs11439060 and rs9138 polymorphisms, previously reported as autoimmune risk variants, in the rheumatoid arthritis (RA) genetic background according to anti-citrullinated protein antibodies (ACPAs) status of RA individuals. METHODS: We analysed a total of 11,715 RA cases and 26,493 controls from nine independent cohorts; all individuals were genotyped or had imputed genotypes for SPP1 rs11439060 and rs9138. The effect of the SPP1 rs11439060 and rs9138 risk-allele combination on osteopontin (OPN) expression in macrophages and OPN serum levels was investigated. RESULTS: We provide evidence for a distinct contribution of SPP1 to RA susceptibility according to ACPA status: the combination of ≥3 SPP1 rs11439060 and rs9138 common alleles was associated mainly with ACPA negativity (p=1.29×10(-5), ORACPA-negative 1.257 (1.135 to 1.394)) and less with ACPA positivity (p=0.0148, ORACPA-positive 1.072 (1.014 to 1.134)). The ORs between these subgroups (ie, ACPA-positive and ACPA-negative) significantly differed (p=7.33×10(-3)). Expression quantitative trait locus analysis revealed an association of the SPP1 risk-allele combination with decreased SPP1 expression in peripheral macrophages from 599 individuals. To corroborate these findings, we found an association of the SPP1 risk-allele combination and low serum level of secreted OPN (p=0.0157), as well as serum level of secreted OPN correlated positively with ACPA production (p=0.005; r=0.483). CONCLUSIONS: We demonstrate a significant contribution of the combination of SPP1 rs11439060 and rs9138 frequent alleles to risk of RA, the magnitude of the association being greater in patients negative for ACPAs.

Juvenile myelomonocytic leukaemia and Noonan syndrome.

BACKGROUND: Infants with Noonan syndrome (NS) are predisposed to developing juvenile myelomonocytic leukaemia (JMML) or JMML-like myeloproliferative disorders (MPD). Whereas sporadic JMML is known to be aggressive, JMML occurring in patients with NS is often considered as benign and transitory. However, little information is available regarding the occurrence and characteristics of JMML in NS. METHODS AND RESULTS: Within a large prospective cohort of 641 patients with a germline PTPN11 mutation, we identified MPD features in 36 (5.6%) patients, including 20 patients (3%) who fully met the consensus diagnostic criteria for JMML. Sixty percent of the latter (12/20) had severe neonatal manifestations, and 10/20 died in the first month of life. Almost all (11/12) patients with severe neonatal JMML were males. Two females who survived MPD/JMML subsequently developed another malignancy during childhood. Although the risk of developing MPD/JMML could not be fully predicted by the underlying PTPN11 mutation, some germline PTPN11 mutations were preferentially associated with myeloproliferation: 10/48 patients with NS (20.8%) with a mutation in codon Asp61 developed MPD/JMML in infancy. Patients with a p.Thr73Ile mutation also had more chances of developing MPD/JMML but with a milder clinical course. SNP array and whole exome sequencing in paired tumoral and constitutional samples identified no second acquired somatic mutation to explain the occurrence of myeloproliferation. CONCLUSIONS: JMML represents the first cause of death in PTPN11-associated NS. Few patients have been reported so far, suggesting that JMML may sometimes be overlooked due to early death, comorbidities or lack of confirmatory tests.

A large French case-control study emphasizes the role of rare Mc1R variants in melanoma risk.

BACKGROUND: The MC1R gene implicated in melanogenesis and skin pigmentation is highly polymorphic. Several alleles are associated with red hair and fair skin phenotypes and contribute to melanoma risk. OBJECTIVE: This work aims to assess the effect of different classes of MC1R variants, notably rare variants, on melanoma risk. Methods. MC1R coding region was sequenced in 1131 melanoma patients and 869 healthy controls. MC1R variants were classified as RHC (R) and non-RHC (r). Rare variants (frequency < 1%) were subdivided into two subgroups, predicted to be damaging (D) or not (nD). RESULTS: Both R and r alleles were associated with melanoma (OR = 2.66 [2.20-3.23] and 1.51 [1.32-1.73]) and had similar population attributable risks (15.8% and 16.6%). We also identified 69 rare variants, of which 25 were novel. D variants were strongly associated with melanoma (OR = 2.38 [1.38-4.15]) and clustered in the same MC1R domains as R alleles (intracellular 2, transmembrane 2 and 7). CONCLUSION: This work confirms the role of R and r alleles in melanoma risk in the French population and proposes a novel class of rare D variants as important melanoma risk factors. These findings may improve the definition of high-risk subjects that could be targeted for melanoma prevention and screening.

The MUC5B variant is associated with idiopathic pulmonary fibrosis but not with systemic sclerosis interstitial lung disease in the European Caucasian population.

A polymorphism on the MUC5B promoter (rs35705950) has been associated with idiopathic pulmonary fibrosis (IPF) but not with systemic sclerosis (SSc) with interstitial lung disease (ILD). We genotyped the MUC5B promoter in the first 142 patients of the French national prospective cohort of IPF, in 981 French patients with SSc (346 ILD), 598 Italian patients with SSc (207 ILD), 1383 French controls and 494 Italian controls. A meta-analysis was performed including all American data available. The T risk allele was present in 41.9% of the IPF patients, 10.8% of the controls (P = 2 × 10(-44)), OR 6.3 [4.6-8.7] for heterozygous patients and OR 21.7 [10.4-45.3] for homozygous patients. Prevalence of the T allele was not modified according to age, gender, smoking in IPF patients. However, none of the black patients with IPF presented the T allele. The prevalence of the T risk allele was similar between French (10%) and Italian (12%) cohorts of SSc whatever the presence of an ILD (11.1% and 13.5%, respectively). Meta-analysis confirmed the similarity between French, Italian and American cohorts of IPF or SSc-ILD. This study confirms 1) an association between the T allele risk and IPF, 2) an absence of association with SSc-ILD, suggesting different pathophysiology.