Feasibility and Impact of Embedding an Extended DNA and RNA Tissue-Based Sequencing Panel for the Routine Care of Patients with Advanced Melanoma in Spain.

Targeted NGS allows a fast and efficient multi-gene analysis and the detection of key gene aberrations in melanoma. In this study, we aim to describe the genetic alterations in a series of 87 melanoma cases using the oncomine focus assay (OFA), relate these results with the clinicopathological features of the patients, and compare them with our previous study results in which we used a smaller panel, the oncomine solid tumor (OST) DNA kit. Patients diagnosed with advanced melanoma at our center from 2020 to 2022 were included and DNA and RNA were extracted for sequencing. Common mutated genes were BRAF (29%), NRAS (28%), ALK, KIT, and MAP2K1 (5% each). Co-occurring mutations were detected in 29% of the samples, including BRAF with KIT, CTNNB1, EGFR, ALK, HRAS, or MAP2K1. Amplifications and rearrangements were detected in 5% of cases. Only BRAF mutation showed a significant statistical association with sun exposure. For patients with a given genetic profile, the melanoma survival and recurrence-free survival rates were equivalent, but not for stage and LDH values. This expanded knowledge of molecular alterations has helped to more comprehensively characterize our patients and has provided relevant information for deciding the best treatment strategy.

Genomic Chaos (Multiple Copy Number Variations and Structural Reorganization) Detected in Two Prenatal Cases.

The use of new technologies in the routine diagnosis of constitutional abnormalities, such as high-resolution chromosomal microarray and next-generation sequencing, has unmasked new mechanisms for generating structural variation of the human genome. For example, complex chromosome rearrangements can originate by a chromosome catastrophe phenomenon in which numerous genomic rearrangements are apparently acquired in a single catastrophic event. This phenomenon is named chromoanagenesis (from the Greek "chromo" for chromosome and "anagenesis" for rebirth). Herein, we report 2 cases of genomic chaos detected at prenatal diagnosis. The terms "chromothripsis" and "chromoanasynthesis" and the challenge of genetic counseling are discussed.

Interactions between rheumatoid arthritis antibodies are associated with the response to anti-tumor necrosis factor therapy.

BACKGROUND: Blocking of the Tumor Necrosis Factor (TNF) activity is a successful therapeutic approach for 50-60% of rheumatoid arthritis (RA) patients. However, there are yet no biomarkers to stratify patients for anti-TNF therapy. Rheumatoid factor (RF) and anti-cyclic-citrullinated antibodies (anti-CCP) have been evaluated as biomarkers of response but the results have shown limited consistency. Anti-carbamylated protein (anti-CarP) and anti-peptidylarginine deiminase type 4 (anti-PAD4) antibodies have been much less studied. Despite being linked to common immune processes, the interaction between these markers has not been evaluated yet. Our aim was to analyze the interaction between these four antibodies in relation to the response to anti-TNF therapy. METHODS: For this objective, a prospective cohort of n = 80 RA patients starting anti-TNF therapy was recruited. Serum determinations at baseline were performed for RF, anti-CCP, anti-CarP and anti-PAD4 antibodies using enzyme-linked immunosorbent assays (ELISA). The clinical response to anti-TNF therapy was determined at week 12 using the change in DAS28 score. Association was performed using multivariate linear regression adjusting for baseline DAS28, sex and age. RESULTS: The interaction between pairs of antibodies was tested by the addition of an interaction term. We found two highly significant antibody interactions associated with treatment response: anti-CarP with anti-PAD4 (p = 0.0062), and anti-CCP with RF (p = 0.00068). The latter antibody interaction was replicated in an independent retrospective cohort of RA patients (n = 199, p = 0.04). CONCLUSIONS: The results of this study suggest that antibody interaction effects are important factors in the response to anti-TNF therapy in RA.

Epigenome-wide association study of rheumatoid arthritis identifies differentially methylated loci in B cells.

Epigenetic regulation of immune cell types could be critical for the development and maintenance of autoimmune diseases like rheumatoid arthritis (RA). B cells are highly relevant in RA, since patients express autoantibodies and depleting this cell type is a successful therapeutic approach. Epigenetic variation, such as DNA methylation, may mediate the pathogenic activity of B cells. In this study, we performed an epigenome-wide association study (EWAS) for RA with three different replication cohorts, to identify disease-specific alterations in DNA methylation in B cells. CpG methylation in isolated B lymphocytes was assayed on the Illumina HumanMethylation450 BeadChip in a discovery cohort of RA patients (N = 50) and controls (N = 75). Differential methylation was observed in 64 CpG sites (q < 0.05). Six biological pathways were also differentially methylated in RA B cells. Analysis in an independent cohort of patients (N = 15) and controls (N = 15) validated the association of 10 CpG sites located on 8 genes CD1C, TNFSF10, PARVG, NID1, DHRS12, ITPK1, ACSF3 and TNFRSF13C, and 2 intergenic regions. Differential methylation at the CBL signaling pathway was replicated. Using an additional case-control cohort (N = 24), the association between RA risk and CpGs cg18972751 at CD1C (P = 2.26 × 10-9) and cg03055671 at TNFSF10 (P = 1.67 × 10-8) genes was further validated. Differential methylation at genes CD1C, TNFSF10, PARVG, NID1, DHRS12, ITPK1, ACSF3, TNFRSF13C and intergenic region chr10p12.31 was replicated in a cohort of systemic lupus erythematosus (SLE) patients (N = 47) and controls (N = 56). Our results highlight genes that may drive the pathogenic activity of B cells in RA and suggest shared methylation patterns with SLE.

A New Case with Corpus Callosum Abnormalities, Microcephaly and Seizures Associated with a 2.3-Mb 1q43-q44 Deletion.

1q44 deletion is a rare syndrome associated with facial dysmorphism and developmental delay, in particular related with expressive speech, seizures, and hypotonia (ORPHA:238769). Until today, the distinct genetic causes for the different symptoms remain not entirely clear. We present a patient with a 2.3-Mb 1q44 deletion, including AKT3, ZBTB18, and HNRNPU, who shows microcephaly, developmental delay, abnormal corpus callosum, and seizures. The genetic findings in this case and a review of the literature spotlight a region between 243 Mb and 245 Mb on chromosome 1q related to the genesis of the typical symptoms of 1q44 deletion.

Genetic association between CD96 locus and immunogenicity to anti-TNF therapy in Crohn's disease.

The production of antibodies to anti-tumor necrosis factor alpha (TNF) agents is one of the main causes of treatment failure in Crohn's disease (CD). To date, however, the contribution of genetics to anti-TNF immunogenicity in CD is still unknown. The objective of the present study was to identify genetic variation associated with anti-TNF immunogenicity in CD. We performed a two-stage genome-wide association study in a cohort of 96 and 123 adalimumab-treated patients, respectively. In the discovery stage, we identified a genome-wide significant association between the CD96 locus and the production of antibodies to anti-TNF treatment (P = 1.88e-09). This association was validated in the replication stage (P < 0.05). The risk allele for anti-TNF immunogenicity was found to be also associated with a lack of response to anti-TNF therapy (P = 0.019). These findings represent an important step toward the understanding of the immunogenicity-based mechanisms that underlie anti-TNF response in CD.

A deletion at ADAMTS9-MAGI1 locus is associated with psoriatic arthritis risk.

OBJECTIVE: Copy number variants (CNVs) have been associated with the risk to develop multiple autoimmune diseases. Our objective was to identify CNVs associated with the risk to develop psoriatic arthritis (PsA) using a genome-wide analysis approach. METHODS: A total of 835 patients with PsA and 1498 healthy controls were genotyped for CNVs using the Illumina HumanHap610 BeadChip genotyping platform. Genomic CNVs were characterised using CNstream analysis software and analysed for association using the χ(2) test. The most significant genomic CNV associations with PsA risk were independently tested in a validation sample of 1133 patients with PsA and 1831 healthy controls. In order to test for the specificity of the variants with PsA aetiology, we also analysed the association to a cohort of 822 patients with purely cutaneous psoriasis (PsC). RESULTS: A total of 165 common CNVs were identified in the genome-wide analysis. We found a highly significant association of an intergenic deletion between ADAMTS9 and MAGI1 genes on chromosome 3p14.1 (p=0.00014). Using the independent patient and control cohort, we validated the association between ADAMTS9-MAGI1 deletion and PsA risk (p=0.032). Using next-generation sequencing, we characterised the 26 kb associated deletion. Finally, analysing the PsC cohort we found a lower frequency of the deletion compared with the PsA cohort (p=0.0088) and a similar frequency to that of healthy controls (p>0.3). CONCLUSIONS: The present genome-wide scan for CNVs associated with PsA risk has identified a new deletion associated with disease risk and which is also differential from PsC risk.

Seven-chain adaptive immune receptor repertoire analysis in rheumatoid arthritis reveals novel features associated with disease and clinically relevant phenotypes.

BACKGROUND: In rheumatoid arthritis (RA), the activation of T and B cell clones specific for self-antigens leads to the chronic inflammation of the synovium. Here, we perform an in-depth quantitative analysis of the seven chains that comprise the adaptive immune receptor repertoire (AIRR) in RA. RESULTS: In comparison to controls, we show that RA patients have multiple and strong differences in the B cell receptor repertoire including reduced diversity as well as altered isotype, chain, and segment frequencies. We demonstrate that therapeutic tumor necrosis factor inhibition partially restores this alteration but find a profound difference in the underlying biochemical reactivities between responders and non-responders. Combining the AIRR with HLA typing, we identify the specific T cell receptor repertoire associated with disease risk variants. Integrating these features, we further develop a molecular classifier that shows the utility of the AIRR as a diagnostic tool. CONCLUSIONS: Simultaneous sequencing of the seven chains of the human AIRR reveals novel features associated with the disease and clinically relevant phenotypes, including response to therapy. These findings show the unique potential of AIRR to address precision medicine in immune-related diseases.

Downregulation of Sfrp5 promotes beta cell proliferation during obesity in the rat.

AIMS/HYPOTHESIS: During obesity, the increment in beta cell mass in response to the rising demand for insulin is essential to maintain normal glucose homeostasis. However, the precise cellular and molecular mechanisms involved in beta cell mass plasticity remain poorly understood. The Wnt signalling pathway has been suggested as one possible modulator of beta cell proliferation, which represents the principal process involved in beta cell mass expansion. Here, we sought to determine the mechanisms involved in beta cell mass proliferation using diet-induced obese rats. METHODS: Wistar rats aged 8 weeks old were fed a standard or cafeteria diet. Global transcriptomic analysis of pancreatic rat islets was performed using microarray analysis. Genetic loss-of-function approaches were performed in dispersed primary rat islets and the beta cell line INS1E. Gene expression was measured by real-time PCR, protein levels by immunoblot analysis, proliferation rates by ELISA and apoptosis by flow cytometry. RESULTS: Sfrp5, coding for secreted frizzled-related protein 5, is downregulated in the pancreatic islets of cafeteria-diet-fed rats as well as in the pancreatic islets of human obese patients. We demonstrate that silencing Sfrp5 increases beta cell proliferation, which correlates with activation of Wnt signalling and enhanced levels of proliferation markers. In addition, we show that expression of Sfrp5 in beta cells is modulated by IGF binding protein 3 (IGFBP3) secreted from visceral adipose tissue. CONCLUSIONS/INTERPRETATION: Together, these findings reveal an important role for SFRP5 and Wnt signalling in the regulation of beta cell proliferation in obesity.

Genome-wide transcriptional analysis of T cell activation reveals differential gene expression associated with psoriasis.

BACKGROUND: Psoriasis is a chronic autoimmune disease in which T cells have a predominant role in initiating and perpetuating the chronic inflammation in skin. However, the mechanisms that regulate T cell activation in psoriasis are still incompletely understood. The objective of the present study was to characterize the main genetic pathways associated with T cell activation in psoriasis. RESULTS: Gene expression profiles from in vitro activated T cells were obtained from 17 psoriasis patients and 7 healthy controls using Illumina HT-12 v4 microarrays. From a total of 47,321 analyzed transcripts, 42 genes were found to be differentially expressed between psoriasis and controls (FDR p-value < 0.1, absolute fold-change > 1.2). Using an independent cohort of 8 patients and 8 healthy controls we validated the overexpression of SPATS2L (p-value =0.0009) and KLF6 (p-value =0.0012) genes in activated T cells from psoriasis patients. Using weighted correlation analysis we identified SPATS2L and KLF6 coexpression networks, which were also significantly associated with psoriasis (p-value < 0.05). Gene Ontology analysis allowed the identification of several biological processes associated with each coexpression network. Finally, using Gene Set Enrichment Analysis over the global T cell transcriptome we also found additional genetic pathways strongly associated with psoriasis (p-value < 0.0001). CONCLUSIONS: This study has identified two new genes, SPATS2L and KLF6, strongly associated with T cell activation in psoriasis. Functional analyses of the gene expression profiles also revealed new biological processes and genetic pathways associated with psoriasis. The results of this study provide an important insight into the biology of this common chronic inflammatory disease.

Trisomy 18p caused by a supernumerary marker with a chromosome 13/21 centromere: a possible recurrent chromosome aberration.

We present a clinical and molecular cytogenetic characterization of two new patients with a complex supernumerary marker consisting of the entire short arm of chromosome 18 with a chromosome 13/21 centromere. One patient is a girl with a nonsyndromic intellectual disability and the second is a prenatally diagnosed fetus. To our knowledge, these are the fourth and fifth such cases to be described in the literature, suggesting the existence of a possible recurring constitutional structural chromosome abnormality.

Sampling time-dependent artifacts in single-cell genomics studies.

Robust protocols and automation now enable large-scale single-cell RNA and ATAC sequencing experiments and their application on biobank and clinical cohorts. However, technical biases introduced during sample acquisition can hinder solid, reproducible results, and a systematic benchmarking is required before entering large-scale data production. Here, we report the existence and extent of gene expression and chromatin accessibility artifacts introduced during sampling and identify experimental and computational solutions for their prevention.

A Combined Transcriptomic and Genomic Analysis Identifies a Gene Signature Associated With the Response to Anti-TNF Therapy in Rheumatoid Arthritis.

Background: Rheumatoid arthritis (RA) is the most frequent autoimmune disease involving the joints. Although anti-TNF therapies have proven effective in the management of RA, approximately one third of patients do not show a significant clinical response. The objective of this study was to identify new genetic variation associated with the clinical response to anti-TNF therapy in RA. Methods: We performed a sequential multi-omic analysis integrating different sources of molecular information. First, we extracted the RNA from synovial biopsies of 11 RA patients starting anti-TNF therapy to identify gene coexpression modules (GCMs) in the RA synovium. Second, we analyzed the transcriptomic association between each GCM and the clinical response to anti-TNF therapy. The clinical response was determined at week 14 using the EULAR criteria. Third, we analyzed the association between the GCMs and anti-TNF response at the genetic level. For this objective, we used genome-wide data from a cohort of 348 anti-TNF treated patients from Spain. The GCMs that were significantly associated with the anti-TNF response were then tested for validation in an independent cohort of 2,706 anti-TNF treated patients. Finally, the functional implication of the validated GCMs was evaluated via pathway and cell type epigenetic enrichment analyses. Results: A total of 149 GCMs were identified in the RA synovium. From these, 13 GCMs were found to be significantly associated with anti-TNF response (P < 0.05). At the genetic level, we detected two of the 13 GCMs to be significantly associated with the response to adalimumab (P = 0.0015) and infliximab (P = 0.021) in the Spain cohort. Using the independent cohort of RA patients, we replicated the association of the GCM associated with the response to adalimumab (P = 0.0019). The validated module was found to be significantly enriched for genes involved in the nucleotide metabolism (P = 2.41e-5) and epigenetic marks from immune cells, including CD4+ regulatory T cells (P = 0.041). Conclusions: These findings show the existence of a drug-specific genetic basis for anti-TNF response, thereby supporting treatment stratification in the search for response biomarkers in RA.

Genome-wide pathway analysis identifies VEGF pathway association with oral ulceration in systemic lupus erythematosus.

BACKGROUND: Systemic lupus erythematosus (SLE) is a genetically complex rheumatic disease characterized by heterogeneous clinical manifestations of unknown etiology. Recent studies have suggested the existence of a genetic basis for SLE heterogeneity. The objective of the present study was to identify new genetic variation associated with the clinically relevant phenotypes in SLE. METHODS: A two-stage pathway-based approach was used to identify the genetic variation associated with the main clinical phenotypes in SLE. In the discovery stage, 482 SLE patients were genotyped using Illumina Human Quad610 microarrays. Association between 798 reference genetic pathways from the Molecular Signatures Database and 11 SLE phenotypes was tested using the set-based method implemented in PLINK software. Pathways significantly associated after multiple test correction were subsequently tested for replication in an independent cohort of 425 SLE patients. Using an in silico approach, we analyzed the functional effects of common SLE therapies on the replicated genetic pathways. The association of known SLE risk variants with the development of the clinical phenotypes was also analyzed. RESULTS: In the discovery stage, we found a significant association between the vascular endothelial growth factor (VEGF) pathway and oral ulceration (P value for false discovery rate (P FDR) < 0.05), and between the negative regulation signaling pathway of retinoic acid inducible gene-I/melanoma differentiation associated gene 5 and the production of antinuclear antibodies (P FDR < 0.05). In the replication stage, we validated the association between the VEGF pathway and oral ulceration. Therapies commonly used to treat mucocutaneous phenotypes in SLE were found to strongly influence VEGF pathway gene expression (P = 4.60e-4 to 5.38e-14). Analysis of known SLE risk loci identified a strong association between PTPN22 and the risk of hematologic disorder and with the development of antinuclear antibodies. CONCLUSIONS: The present study has identified VEGF genetic pathway association with the risk of oral ulceration in SLE. New therapies targeting the VEGF pathway could be more effective in reducing the severity of this phenotype. These findings represent a first step towards the understanding of the genetic basis of phenotype heterogeneity in SLE.

Genetic variation associated with cardiovascular risk in autoimmune diseases.

Autoimmune diseases have a higher prevalence of cardiovascular events compared to the general population. The objective of this study was to investigate the genetic basis of cardiovascular disease (CVD) risk in autoimmunity. We analyzed genome-wide genotyping data from 6,485 patients from six autoimmune diseases that are associated with a high socio-economic impact. First, for each disease, we tested the association of established CVD risk loci. Second, we analyzed the association of autoimmune disease susceptibility loci with CVD. Finally, to identify genetic patterns associated with CVD risk, we applied the cross-phenotype meta-analysis approach (CPMA) on the genome-wide data. A total of 17 established CVD risk loci were significantly associated with CVD in the autoimmune patient cohorts. From these, four loci were found to have significantly different genetic effects across autoimmune diseases. Six autoimmune susceptibility loci were also found to be associated with CVD risk. Genome-wide CPMA analysis identified 10 genetic clusters strongly associated with CVD risk across all autoimmune diseases. Two of these clusters are highly enriched in pathways previously associated with autoimmune disease etiology (TNFα and IFNγ cytokine pathways). The results of this study support the presence of specific genetic variation associated with the increase of CVD risk observed in autoimmunity.

Genome-Wide Pathway Analysis Identifies Genetic Pathways Associated with Psoriasis.

Psoriasis is a chronic inflammatory disease with a complex genetic architecture. To date, the psoriasis heritability is only partially explained. However, there is increasing evidence that the missing heritability in psoriasis could be explained by multiple genetic variants of low effect size from common genetic pathways. The objective of this study was to identify new genetic variation associated with psoriasis risk at the pathway level. We genotyped 598,258 single nucleotide polymorphisms in a discovery cohort of 2,281 case-control individuals from Spain. We performed a genome-wide pathway analysis using 1,053 reference biological pathways. A total of 14 genetic pathways (PFDR ≤ 2.55 × 10(-2)) were found to be significantly associated with psoriasis risk. Using an independent validation cohort of 7,353 individuals from the UK, a total of 6 genetic pathways were significantly replicated (PFDR ≤ 3.46 × 10(-2)). We found genetic pathways that had not been previously associated with psoriasis risk such as retinol metabolism (Pcombined = 1.84 × 10(-4)), the transport of inorganic ions and amino acids (Pcombined = 1.57 × 10(-7)), and post-translational protein modification (Pcombined = 1.57 × 10(-7)). In the latter pathway, MGAT5 showed a strong network centrality, and its association with psoriasis risk was further validated in an additional case-control cohort of 3,429 individuals (P < 0.05). These findings provide insights into the biological mechanisms associated with psoriasis susceptibility.

Variation at interleukin-6 receptor gene is associated to joint damage in rheumatoid arthritis.

INTRODUCTION: Interleukin-6 (IL-6) cytokine signaling is key in Rheumatoid Arthritis (RA) pathophysiology. Blocking IL-6 receptor (IL6R) has proven to be a highly effective treatment to prevent joint damage. This study was performed to investigate the association between the genetic variation at IL6R gene and the severity of joint damage in RA. METHODS: IL6R gene tagging SNPs (n = 5) were genotyped in a discovery group of 527 RA patients from 5 different university hospitals from Spain. For each marker, a linear regression analysis was performed using an additive model and adjusting for the years of evolution of the disease, autoantibody status, gender and age. Haplotypes combining the SNPs were also estimated and tested for association with the level of joint destruction. Using an independent cohort of 705 RA patients from 6 university hospitals we performed a validation study of the SNPs associated in the discovery phase. RESULTS: In the discovery group we found a highly significant association between IL6R SNP rs4845618 and the level of joint destruction in RA (P = 0.0058, P corrected = 0.026), and a moderate association with SNP rs4453032 (P = 0.02, P corrected = 0.05). The resulting haplotype from both SNPs was more significantly associated with joint damage (P = 0.0037, P corrected = 0.011). Using the validation cohort, we replicated the association between the two IL-6R SNPs with the degree of joint destruction in RA (P = 0.007 and P = 0.04, meta-analysis P = 0.00011 and P = 0.0021, respectively), and the haplotype association (P = 0.0058, meta-analysis P = 6.64 e-5). CONCLUSIONS: Genetic variation at IL6R gene is associated with joint damage in RA.

GWAS replication study confirms the association of PDE3A-SLCO1C1 with anti-TNF therapy response in rheumatoid arthritis.

AIM: The present study was undertaken to replicate the association of candidate genes for anti-TNF response in rheumatoid arthritis. Candidate genes were selected from a recent genome-wide association study on anti-TNF response performed in a population from Denmark. MATERIALS & METHODS: Genomic DNA was obtained from 315 Spanish rheumatoid arthritis patients having received an anti-TNF agent as their first biological therapy. SNPs from NR2FR2, MAP2K6, CBLN2 and PDE3A-SLCO1C1 candidate loci were genotyped. RESULTS: The PDE3A-SLCO1C1 locus rs3794271 SNP showed a highly significant association with anti-TNF treatment response (p = 1.74 × 10⁻5). Combining the statistical evidence from the Spanish and Danish rheumatoid arthritis cohorts, the associated rs3794271 SNP reached a genome-wide significance level of association (p = 3.3 × 10⁻¹°). CONCLUSION: The present findings establish the PDE3A-SLCO1C1 locus as a strong genetic marker of anti-TNF therapy response.

Role of IGFBP-3 in the regulation of ß-cell mass during obesity: adipose tissue/ß-cell cross talk.

In obesity an increase in ß-cell mass occurs to cope with the rise in insulin demand. This ß-cell plasticity is essential to avoid the onset of hyperglycemia, although the molecular mechanisms that regulate this process remain unclear. This study analyzed the role of adipose tissue in the control of ß-cell replication. Using a diet-induced model of obesity, we obtained conditioned media from three different white adipose tissue depots. Only in the adipose tissue depot surrounding the pancreas did the diet induce changes that led to an increase in INS1E cells and the islet replication rate. To identify the factors responsible for this proliferative effect, adipose tissue gene expression analysis was conducted by microarrays and quantitative RT-PCR. Of all the differentially expressed proteins, only the secreted ones were studied. IGF binding protein 3 (Igfbp3) was identified as the candidate for this effect. Furthermore, in the conditioned media, although the blockage of IGFBP3 led to an increase in the proliferation rate, the blockage of IGF-I receptor decreased it. Taken together, these data show that obesity induces specific changes in the expression profile of the adipose tissue depot surrounding the pancreas, leading to a decrease in IGFBP3 secretion. This decrease acts in a paracrine manner, stimulating the ß-cell proliferation rate, probably through an IGF-I-dependent mechanism. This cross talk between the visceral-pancreatic adipose tissue and ß-cells is a novel mechanism that participates in the control of ß-cell plasticity.