Gene-Environment Interaction Between CYP1B1 and Oral Contraception on Frontal Fibrosing Alopecia.

Importance: Frontal fibrosing alopecia (FFA) is an increasingly prevalent form of follicular lichen planus, causing irreversible hair loss predominantly in postmenopausal individuals. An earlier genome-wide meta-analysis of female FFA identified risk loci in genes implicated in self-antigen presentation and T-cell homeostasis, including HLA-B*07:02, ST3GAL1, and SEMA4B. However, CYP1B1, which is important for hormone metabolism, was also implicated with the substitution of serine for asparagine at position 453 (c.1358A>G, p.Asn453Ser) exhibiting a protective effect against FFA. Increasing understanding of genetic and environmental variables and their interactions will improve understanding of disease pathogenesis and has the potential to inform risk mitigation strategies. Objective: To investigate whether oral contraceptive pill (OCP) use modulates the protective effect of the common missense variant in CYP1B1 (c.1358A>G, p.Asn453Ser) on FFA risk. Design, Setting, and Participants: This gene-environment interaction study using a case-control design enrolled female patients with FFA from UK-based dermatology clinics. The patients were matched with unrelated age- and ancestry-matched female control individuals derived from UK Biobank in a 1:66 ratio, determined by the first 4 principal components from genome-wide genotypes. Data were collected from July 2015 to September 2017, and analyzed from October 2022 to December 2023. Main Outcome and Measure: The main outcomes were the modulatory effect of OCP use on the contribution of the CYP1B1 missense variant to female FFA risk and a formal gene-environment interaction test evaluated by a logistic regression model with a multiplicative interaction term, under the assumptions of an additive genetic model interaction term, under the assumptions of an additive genetic model. Results: Of the 489 female patients with FFA, the mean (SD) age was 65.8 (9.7) years, and 370 (75.7%) had a history of OCP use. Of the 34 254 age- and ancestry-matched control individuals, the mean (SD) age was 65.0 (8.4) years, and previous OCP use was reported in 31 177 (91.0%). An association between female FFA and the CYP1B1 risk allele was observed in individuals who reported OCP use (odds ratio, 1.90 [95% CI, 1.50-2.40]; P = 8.41 × 10-8) but not in those with no documented exposure to OCPs (odds ratio, 1.16 [95% CI, 0.82-1.64]; P = .39). A full gene-environment interaction model demonstrated a significant additive statistical interaction between c.1358A, p.453Asn, and history of OCP use on FFA risk (OR for interaction, 1.63 [95% CI, 1.07-2.46]; P = .02). Conclusions and Relevance: This gene-environment interaction analysis suggests that the protective effect of the CYP1B1 missense variant on FFA risk might be mediated by exposure to OCPs. The allele that encodes an asparagine at position 453 of CYP1B1 was associated with increased odds of FFA only in participants with OCP history.

Development of antidrug antibodies against adalimumab maps to variation within the HLA-DR peptide-binding groove.

Targeted biologic therapies can elicit an undesirable host immune response characterized by the development of antidrug antibodies (ADA), an important cause of treatment failure. The most widely used biologic across immune-mediated diseases is adalimumab, a tumor necrosis factor inhibitor. This study aimed to identify genetic variants that contribute to the development of ADA against adalimumab, thereby influencing treatment failure. In patients with psoriasis on their first course of adalimumab, in whom serum ADA had been evaluated 6-36 months after starting treatment, we observed a genome-wide association with ADA against adalimumab within the major histocompatibility complex (MHC). The association signal mapped to the presence of tryptophan at position 9 and lysine at position 71 of the HLA-DR peptide-binding groove, with both residues conferring protection against ADA. Underscoring their clinical relevance, these residues were also protective against treatment failure. Our findings highlight antigenic peptide presentation via MHC class II as a critical mechanism in the development of ADA against biologic therapies and downstream treatment response.

Comparative Genetic Analysis of Psoriatic Arthritis and Psoriasis for the Discovery of Genetic Risk Factors and Risk Prediction Modeling.

OBJECTIVES: Psoriatic arthritis (PsA) has a strong genetic component, and the identification of genetic risk factors could help identify the ~30% of psoriasis patients at high risk of developing PsA. Our objectives were to identify genetic risk factors and pathways that differentiate PsA from cutaneous-only psoriasis (PsC) and to evaluate the performance of PsA risk prediction models. METHODS: Genome-wide meta-analyses were conducted separately for 5,065 patients with PsA and 21,286 healthy controls and separately for 4,340 patients with PsA and 6,431 patients with PsC. The heritability of PsA was calculated as a single-nucleotide polymorphism (SNP)-based heritability estimate (h2 SNP ) and biologic pathways that differentiate PsA from PsC were identified using Priority Index software. The generalizability of previously published PsA risk prediction pipelines was explored, and a risk prediction model was developed with external validation. RESULTS: We identified a novel genome-wide significant susceptibility locus for the development of PsA on chromosome 22q11 (rs5754467; P = 1.61 × 10-9 ), and key pathways that differentiate PsA from PsC, including NF-κB signaling (adjusted P = 1.4 × 10-45 ) and Wnt signaling (adjusted P = 9.5 × 10-58 ). The heritability of PsA in this cohort was found to be moderate (h2 SNP = 0.63), which was similar to the heritability of PsC (h2 SNP = 0.61). We observed modest performance of published classification pipelines (maximum area under the curve 0.61), with similar performance of a risk model derived using the current data. CONCLUSION: Key biologic pathways associated with the development of PsA were identified, but the investigation of risk classification revealed modest utility in the available data sets, possibly because many of the PsC patients included in the present study were receiving treatments that are also effective in PsA. Future predictive models of PsA should be tested in PsC patients recruited from primary care.

Genome-wide association meta-analysis identifies 29 new acne susceptibility loci.

Acne vulgaris is a highly heritable skin disorder that primarily impacts facial skin. Severely inflamed lesions may leave permanent scars that have been associated with long-term psychosocial consequences. Here, we perform a GWAS meta-analysis comprising 20,165 individuals with acne from nine independent European ancestry cohorts. We identify 29 novel genome-wide significant loci and replicate 14 of the 17 previously identified risk loci, bringing the total number of reported acne risk loci to 46. Using fine-mapping and eQTL colocalisation approaches, we identify putative causal genes at several acne susceptibility loci that have previously been implicated in Mendelian hair and skin disorders, including pustular psoriasis. We identify shared genetic aetiology between acne, hormone levels, hormone-sensitive cancers and psychiatric traits. Finally, we show that a polygenic risk score calculated from our results explains up to 5.6% of the variance in acne liability in an independent cohort.

Application of information theoretic feature selection and machine learning methods for the development of genetic risk prediction models.

In view of the growth of clinical risk prediction models using genetic data, there is an increasing need for studies that use appropriate methods to select the optimum number of features from a large number of genetic variants with a high degree of redundancy between features due to linkage disequilibrium (LD). Filter feature selection methods based on information theoretic criteria, are well suited to this challenge and will identify a subset of the original variables that should result in more accurate prediction. However, data collected from cohort studies are often high-dimensional genetic data with potential confounders presenting challenges to feature selection and risk prediction machine learning models. Patients with psoriasis are at high risk of developing a chronic arthritis known as psoriatic arthritis (PsA). The prevalence of PsA in this patient group can be up to 30% and the identification of high risk patients represents an important clinical research which would allow early intervention and a reduction of disability. This also provides us with an ideal scenario for the development of clinical risk prediction models and an opportunity to explore the application of information theoretic criteria methods. In this study, we developed the feature selection and psoriatic arthritis (PsA) risk prediction models that were applied to a cross-sectional genetic dataset of 1462 PsA cases and 1132 cutaneous-only psoriasis (PsC) cases using 2-digit HLA alleles imputed using the SNP2HLA algorithm. We also developed stratification method to mitigate the impact of potential confounder features and illustrate that confounding features impact the feature selection. The mitigated dataset was used in training of seven supervised algorithms. 80% of data was randomly used for training of seven supervised machine learning methods using stratified nested cross validation and 20% was selected randomly as a holdout set for internal validation. The risk prediction models were then further validated in UK Biobank dataset containing data on 1187 participants and a set of features overlapping with the training dataset.Performance of these methods has been evaluated using the area under the curve (AUC), accuracy, precision, recall, F1 score and decision curve analysis(net benefit). The best model is selected based on three criteria: the 'lowest number of feature subset' with the 'maximal average AUC over the nested cross validation' and good generalisability to the UK Biobank dataset. In the original dataset, with over 100 different bootstraps and seven feature selection (FS) methods, HLA_C_*06 was selected as the most informative genetic variant. When the dataset is mitigated the single most important genetic features based on rank was identified as HLA_B_*27 by the seven different feature selection methods, consistent with previous analyses of this data using regression based methods. However, the predictive accuracy of these single features in post mitigation was found to be moderate (AUC= 0.54 (internal cross validation), AUC=0.53 (internal hold out set), AUC=0.55(external data set)). Sequentially adding additional HLA features based on rank improved the performance of the Random Forest classification model where 20 2-digit features selected by Interaction Capping (ICAP) demonstrated (AUC= 0.61 (internal cross validation), AUC=0.57 (internal hold out set), AUC=0.58 (external dataset)). The stratification method for mitigation of confounding features and filter information theoretic feature selection can be applied to a high dimensional dataset with the potential confounders.

Coagulation Factor XIII-A Subunit Missense Mutation in the Pathobiology of Autosomal Dominant Multiple Dermatofibromas.

Dermatofibromas are common benign skin lesions, the etiology of which is poorly understood. We identified two unrelated pedigrees in which there was autosomal dominant transmission of multiple dermatofibromas. Whole exome sequencing revealed a rare shared heterozygous missense variant in the F13A1 gene encoding factor XIII subunit A (FXIII-A), a transglutaminase involved in hemostasis, wound healing, tumor growth, and apoptosis. The variant (p.Lys679Met) has an allele frequency of 0.0002 and is predicted to be a damaging mutation. Recombinant human Lys679Met FXIII-A demonstrated reduced fibrin crosslinking activity in vitro. Of note, the treatment of fibroblasts with media containing Lys679Met FXIII-A led to enhanced adhesion, proliferation, and type I collagen synthesis. Immunostaining revealed co-localization between FXIII-A and α4ß1 integrins, more prominently for Lys679Met FXIII-A than the wild type. In addition, both the α4ß1 inhibitors and the mutation of the FXIII-A Isoleucine-Leucine-Aspartate-Threonine (ILDT) motif prevented Lys679Met FXIII-A-dependent proliferation and collagen synthesis of fibroblasts. Our data suggest that the Lys679Met mutation may lead to a conformational change in the FXIII-A protein that enhances α4-integrin binding and provides insight into an unexpected role for FXIII-A in the pathobiology of familial dermatofibroma.

Frequency of Pathogenic Germline Variants in CDH1, BRCA2, CHEK2, PALB2, BRCA1, and TP53 in Sporadic Lobular Breast Cancer.

BACKGROUND: Invasive lobular breast cancer (ILC) accounts for approximately 15% of invasive breast carcinomas and is commonly associated with lobular carcinoma in situ (LCIS). Both have been shown to have higher familial risks than the more common ductal cancers. However, there are little data on the prevalence of the known high and moderate penetrance breast cancer predisposition genes in ILC. The aim of this study was to assess the frequency of germline variants in CDH1, BRCA2, BRCA1, CHEK2, PALB2, and TP53 in sporadic ILC and LCIS diagnosed in women ages ≤60 years. METHODS: Access Array technology (Fluidigm) was used to amplify all exons of CDH1, BRCA2, BRCA1, TP53, CHEK2, and PALB2 using a custom-made targeted sequencing panel in 1,434 cases of ILC and 368 cases of pure LCIS together with 1,611 controls. RESULTS: Case-control analysis revealed an excess of pathogenic variants in BRCA2, CHEK2, PALB2, and CDH1 in women with ILC. CHEK2 was the only gene that showed an association with pure LCIS [OR = 9.90; 95% confidence interval (CI), 3.42-28.66, P = 1.4 × 10-5] with a larger effect size seen in LCIS compared with ILC (OR = 4.31; 95% CI, 1.61-11.58, P = 1.7 × 10-3). CONCLUSIONS: Eleven percent of patients with ILC ages ≤40 years carried germline variants in known breast cancer susceptibility genes. IMPACT: Women with ILC ages ≤40 years should be offered genetic screening using a panel of genes that includes BRCA2, CHEK2, PALB2, and CDH1.

Frequency of pathogenic germline variants in BRCA1, BRCA2, PALB2, CHEK2 and TP53 in ductal carcinoma in situ diagnosed in women under the age of 50 years.

INTRODUCTION: Ductal carcinoma in situ (DCIS) is a non-obligate precursor of invasive ductal breast cancer, and approximately 20% of screen-detected tumours are pure DCIS. Most risk factors for breast cancer have similar associations with DCIS and IDC; however, there is limited data on the prevalence of the known high and moderate penetrance breast cancer predisposition genes in DCIS and which women with DCIS should be referred for genetic screening. The aim of this study was to assess the frequency of germline variants in BRCA2, BRCA1, CHEK2, PALB2 and TP53 in DCIS in women aged less than 50 years of age. METHODS: After DNA extraction from the peripheral blood, Access Array technology (Fluidigm) was used to amplify all exons of these five known breast cancer predisposition genes using a custom made targeted sequencing panel in 655 cases of pure DCIS presenting in women under the age of 50 years together with 1611 controls. RESULTS: Case-control analysis revealed an excess of pathogenic variants in BRCA2 (OR = 27.96, 95%CI 6.56-119.26, P = 2.0 × 10-10) and CHEK2 (OR = 8.04, 95%CI 2.93-22.05, P = 9.0 × 10-6), with weaker associations with PALB2 (P = 0.003), BRCA1 (P = 0.007) and TP53 (P = 0.02). For oestrogen receptor (ER)-positive DCIS the frequency of pathogenic variants was 9% under the age of 50 (14% with a family history of breast cancer) and 29% under the age of 40 (42% with a family history of breast cancer). For ER-negative DCIS, the frequency was 9% (16% with a family history of breast cancer) and 8% (11% with a family history of breast cancer) under the ages of 50 and 40, respectively. CONCLUSIONS: This study has shown that breast tumourigenesis in women with pathogenic variants in BRCA2, CHEK2, PALB2, BRCA1 and TP53 can involve a DCIS precursor stage and that the focus of genetic testing in DCIS should be on women under the age of 40 with ER-positive DCIS.

Author Correction: Novel mutations in PIEZO1 cause an autosomal recessive generalized lymphatic dysplasia with non-immune hydrops fetalis.

This Article contains an error in the last sentence of the 'Variant analysis suggests they are pathogenic' section of the Results, which incorrectly reads 'No truncated PIEZO1 protein products were identified in western blot analysis in GLD1:II.3 and GLD2:II.2 (Fig. 2, Supplementary Fig. 6), suggesting that the truncated protein is not stable and therefore degraded.' This should read 'No full-size PIEZO1 protein products were identified in western blot analysis in GLD1:II.3 and GLD2:II.2 (Fig. 2, Supplementary Fig. 6); the three nonsense mutations are predicted to lead to premature termination of the protein, hence it is possible that those truncated proteins will be non-functional or even unstable and degraded.' The error has not been fixed in the PDF or HTML versions of the Article.

Sequencing of human genomes with nanopore technology.

Whole-genome sequencing (WGS) is becoming widely used in clinical medicine in diagnostic contexts and to inform treatment choice. Here we evaluate the potential of the Oxford Nanopore Technologies (ONT) MinION long-read sequencer for routine WGS by sequencing the reference sample NA12878 and the genome of an individual with ataxia-pancytopenia syndrome and severe immune dysregulation. We develop and apply a novel reference panel-free analytical method to infer and then exploit phase information which improves single-nucleotide variant (SNV) calling performance from otherwise modest levels. In the clinical sample, we identify and directly phase two non-synonymous de novo variants in SAMD9L, (OMIM #159550) inferring that they lie on the same paternal haplotype. Whilst consensus SNV-calling error rates from ONT data remain substantially higher than those from short-read methods, we demonstrate the substantial benefits of analytical innovation. Ongoing improvements to base-calling and SNV-calling methodology must continue for nanopore sequencing to establish itself as a primary method for clinical WGS.

The ADAMTS13-VWF axis is dysregulated in chronic thromboembolic pulmonary hypertension.

Chronic thromboembolic pulmonary hypertension (CTEPH) is an important consequence of pulmonary embolism that is associated with abnormalities in haemostasis. We investigated the ADAMTS13-von Willebrand factor (VWF) axis in CTEPH, including its relationship with disease severity, inflammation, ABO groups and ADAMTS13 genetic variants.ADAMTS13 and VWF plasma antigen levels were measured in patients with CTEPH (n=208), chronic thromboembolic disease without pulmonary hypertension (CTED) (n=35), resolved pulmonary embolism (n=28), idiopathic pulmonary arterial hypertension (n=30) and healthy controls (n=68). CTEPH genetic ABO associations and protein quantitative trait loci were investigated. ADAMTS13-VWF axis abnormalities were assessed in CTEPH and healthy control subsets by measuring ADAMTS13 activity, D-dimers and VWF multimeric size.Patients with CTEPH had decreased ADAMTS13 (adjusted ß -23.4%, 95% CI -30.9- -15.1%, p<0.001) and increased VWF levels (ß +75.5%, 95% CI 44.8-113%, p<0.001) compared to healthy controls. ADAMTS13 levels remained low after reversal of pulmonary hypertension by pulmonary endarterectomy surgery and were equally reduced in CTED. We identified a genetic variant near the ADAMTS13 gene associated with ADAMTS13 protein that accounted for ∼8% of the variation in levels.The ADAMTS13-VWF axis is dysregulated in CTEPH. This is unrelated to pulmonary hypertension, disease severity or markers of systemic inflammation and implicates the ADAMTS13-VWF axis in CTEPH pathobiology.

MED12, TERT promoter and RBM15 mutations in primary and recurrent phyllodes tumours.

BACKGROUND: MED12 and TERT promoter mutations have been shown to be the most common somatic mutations in phyllodes tumours (PTs). The aims of this study were to determine the frequency of these mutations in recurrent PTs, assess whether TERT promoter mutations could be helpful in distinguishing fibroadenomas (FAs) from PTs and identify novel mutations that may be driving malignant progression. METHODS: MED12 and the TERT promoter were Sanger sequenced in 75 primary PTs, 21 recurrences, 19 single FAs and 2 cases of multiple FAs with benign PTs. Whole-exome sequencing was performed on one borderline PT. RESULTS: Recurrent PTs and multiple FAs showed temporal discordance in MED12 but not TERT. Recurrent samples did acquire TERT mutations, with recurrent benign PTs more likely to have mutations in both genes. TERT mutations were not helpful in differentiating between benign PTs and FAs in cases of multiple FAs/PTs. Exome sequencing revealed a nonsense mutation in RBM15 and Sanger sequencing revealed another three RBM15 mutations in malignant/borderline PTs. CONCLUSIONS: This study has shown that MED12 mutations can be heterogeneous in both synchronous and recurrent PTs unlike TERT mutations. We have also shown that RBM15 mutations may be important in the pathogenesis of borderline/malignant PTs.

Biallelic Mutations in KDSR Disrupt Ceramide Synthesis and Result in a Spectrum of Keratinization Disorders Associated with Thrombocytopenia.

Mutations in ceramide biosynthesis pathways have been implicated in a few Mendelian disorders of keratinization, although ceramides are known to have key roles in several biological processes in skin and other tissues. Using whole-exome sequencing in four probands with undiagnosed skin hyperkeratosis/ichthyosis, we identified compound heterozygosity for mutations in KDSR, encoding an enzyme in the de novo synthesis pathway of ceramides. Two individuals had hyperkeratosis confined to palms, soles, and anogenital skin, whereas the other two had more severe, generalized harlequin ichthyosis-like skin. Thrombocytopenia was present in all patients. The mutations in KDSR were associated with reduced ceramide levels in skin and impaired platelet function. KDSR enzymatic activity was variably reduced in all patients, resulting in defective acylceramide synthesis. Mutations in KDSR have recently been reported in inherited recessive forms of progressive symmetric erythrokeratoderma, but our study shows that biallelic mutations in KDSR are implicated in an extended spectrum of disorders of keratinization in which thrombocytopenia is also part of the phenotype. Mutations in KDSR cause defective ceramide biosynthesis, underscoring the importance of ceramide and sphingosine synthesis pathways in skin and platelet biology.

Cross-phenotype association mapping of the MHC identifies genetic variants that differentiate psoriatic arthritis from psoriasis.

OBJECTIVES: Psoriatic arthritis (PsA) is a chronic inflammatory arthritis, with a strong heritable component, affecting patients with psoriasis. Here we attempt to identify genetic variants within the major histocompatibility complex (MHC) that differentiate patients with PsA from patients with cutaneous psoriasis alone (PsC). METHODS: 2808 patients with PsC, 1945 patients with PsA and 8920 population controls were genotyped. We imputed SNPs, amino acids and classical HLA alleles across the MHC and tested for association with PsA compared to population controls and the PsC patient group. In addition we investigated the impact of the age of disease onset on associations. RESULTS: HLA-C*06:02 was protective of PsA compared to PsC (p=9.57×10-66, OR 0.37). The HLA-C*06:02 risk allele was associated with a younger age of psoriasis onset in all patients (p=1.01×10-59). After controlling for the age of psoriasis onset no association of PsA to HLA-C*06:02 (p=0.07) was observed; instead, the most significant association was to amino acid at position 97 of HLA-B (p=1.54×10-9) where the presence of asparagine or serine residue increased PsA risk. Asparagine at position 97 of HLA-B defines the HLA-B*27 alleles. CONCLUSIONS: By controlling for the age of psoriasis onset, we show, for the first time, that HLA-C*06:02 is not associated with PsA and that amino acid position 97 of HLA-B differentiates PsA from PsC. This amino acid also represents the largest genetic effect for ankylosing spondylitis, thereby refining the genetic overlap of these two spondyloarthropathies. Correcting for bias has important implications for cross-phenotype genetic studies.

PRUNE is crucial for normal brain development and mutated in microcephaly with neurodevelopmental impairment.

PRUNE is a member of the DHH (Asp-His-His) phosphoesterase protein superfamily of molecules important for cell motility, and implicated in cancer progression. Here we investigated multiple families from Oman, India, Iran and Italy with individuals affected by a new autosomal recessive neurodevelopmental and degenerative disorder in which the cardinal features include primary microcephaly and profound global developmental delay. Our genetic studies identified biallelic mutations of PRUNE1 as responsible. Our functional assays of disease-associated variant alleles revealed impaired microtubule polymerization, as well as cell migration and proliferation properties, of mutant PRUNE. Additionally, our studies also highlight a potential new role for PRUNE during microtubule polymerization, which is essential for the cytoskeletal rearrangements that occur during cellular division and proliferation. Together these studies define PRUNE as a molecule fundamental for normal human cortical development and define cellular and clinical consequences associated with PRUNE mutation.