Finding exclusively deleted or amplified genomic areas in lung adenocarcinomas using a novel chromosomal pattern analysis

<p>Abstract</p> <p>Background</p> <p>Genomic copy number alteration (CNA) that are recurrent across multiple samples often harbor critical genes that can drive either the initiation or the progression of cancer disease. Up to now, most researchers investigating recurrent CNAs consider separately the marginal frequencies for copy gain or loss and select the areas of interest based on arbitrary cut-off thresholds of these frequencies. In practice, these analyses ignore the interdependencies between the propensity of being deleted or amplified for a clone. In this context, a joint analysis of the copy number changes across tumor samples may bring new insights about patterns of recurrent CNAs.</p> <p>Methods</p> <p>We propose to identify patterns of recurrent CNAs across tumor samples from high-resolution comparative genomic hybridization microarrays. Clustering is achieved by modeling the copy number state (loss, no-change, gain) as a multinomial distribution with probabilities parameterized through a latent class model leading to nine patterns of recurrent CNAs. This model gives us a powerful tool to identify clones with contrasting propensity of being deleted or amplified across tumor samples. We applied this model to a homogeneous series of 65 lung adenocarcinomas.</p> <p>Results</p> <p>Our latent class model analysis identified interesting patterns of chromosomal aberrations. Our results showed that about thirty percent of the genomic clones were classified either as "exclusively" deleted or amplified recurrent CNAs and could be considered as non random chromosomal events. Most of the known oncogenes or tumor suppressor genes associated with lung adenocarcinoma were located within these areas. We also describe genomic areas of potential interest and show that an increase of the frequency of amplification in these particular areas is significantly associated with poorer survival.</p> <p>Conclusion</p> <p>Analyzing jointly deletions and amplifications through our latent class model analysis allows highlighting specific genomic areas with exclusively amplified or deleted recurrent CNAs which are good candidate for harboring oncogenes or tumor suppressor genes.</p

A constrained polynomial regression procedure for estimating the local False Discovery Rate

<p>Abstract</p> <p>Background</p> <p>In the context of genomic association studies, for which a large number of statistical tests are performed simultaneously, the local False Discovery Rate (<it>lFDR</it>), which quantifies the evidence of a specific gene association with a clinical or biological variable of interest, is a relevant criterion for taking into account the multiple testing problem. The <it>lFDR </it>not only allows an inference to be made for each gene through its specific value, but also an estimate of Benjamini-Hochberg's False Discovery Rate (<it>FDR</it>) for subsets of genes.</p> <p>Results</p> <p>In the framework of estimating procedures without any distributional assumption under the alternative hypothesis, a new and efficient procedure for estimating the <it>lFDR </it>is described. The results of a simulation study indicated good performances for the proposed estimator in comparison to four published ones. The five different procedures were applied to real datasets.</p> <p>Conclusion</p> <p>A novel and efficient procedure for estimating <it>lFDR </it>was developed and evaluated.</p

Monocyte NOTCH2 expression predicts interferon-beta immunogenicity in multiple sclerosis patients

Multiple sclerosis (MS) is an autoimmune disease characterized by CNS inflammation leading to demyelination and axonal damage. IFN-β is an established treatment for MS; however, up to 30% of IFN-β–treated MS patients develop neutralizing antidrug antibodies (nADA), leading to reduced drug bioactivity and efficacy. Mechanisms driving antidrug immunogenicity remain uncertain, and reliable biomarkers to predict immunogenicity development are lacking. Using high-throughput flow cytometry, NOTCH2 expression on CD14+ monocytes and increased frequency of proinflammatory monocyte subsets were identified as baseline predictors of nADA development in MS patients treated with IFN-β. The association of this monocyte profile with nADA development was validated in 2 independent cross-sectional MS patient cohorts and a prospective cohort followed before and after IFN-β administration. Reduced monocyte NOTCH2 expression in nADA+ MS patients was associated with NOTCH2 activation measured by increased expression of Notch-responsive genes, polarization of monocytes toward a nonclassical phenotype, and increased proinflammatory IL-6 production. NOTCH2 activation was T cell dependent and was only triggered in the presence of serum from nADA+ patients. Thus, nADA development was driven by a proinflammatory environment that triggered activation of the NOTCH2 signaling pathway prior to first IFN-β administration

Response to Biologic Drugs in Patients with Rheumatoid Arthritis and Antidrug Antibodies

Importance: There are conflicting data on the association of antidrug antibodies with response to biologic disease-modifying antirheumatic drugs (bDMARDs) in rheumatoid arthritis (RA). Objective: To analyze the association of antidrug antibodies with response to treatment for RA. Design, Setting, and Participants: This cohort study analyzed data from the ABI-RA (Anti-Biopharmaceutical Immunization: Prediction and Analysis of Clinical Relevance to Minimize the Risk of Immunization in Rheumatoid Arthritis Patients) multicentric, open, prospective study of patients with RA from 27 recruiting centers in 4 European countries (France, Italy, the Netherlands, and the UK). Eligible patients were 18 years or older, had RA diagnosis, and were initiating a new bDMARD. Recruitment spanned from March 3, 2014, to June 21, 2016. The study was completed in June 2018, and data were analyzed in June 2022. Exposures: Patients were treated with a new bDMARD: adalimumab, infliximab (grouped as anti-tumor necrosis factor [TNF] monoclonal antibodies [mAbs]), etanercept, tocilizumab, and rituximab according to the choice of the treating physician. Main Outcomes and Measures: The primary outcome was the association of antidrug antibody positivity with EULAR (European Alliance of Associations for Rheumatology; formerly, European League Against Rheumatism) response to treatment at month 12 assessed through univariate logistic regression. The secondary end points were the EULAR response at month 6 and at visits from month 6 to months 15 to 18 using generalized estimating equation models. Detection of antidrug antibody serum levels was performed at months 1, 3, 6, 12, and 15 to 18 using electrochemiluminescence (Meso Scale Discovery) and drug concentration for anti-TNF mAbs, and etanercept in the serum was measured using enzyme-linked immunosorbent assay. Results: Of the 254 patients recruited, 230 (mean [SD] age, 54.3 [13.7] years; 177 females [77.0%]) were analyzed. At month 12, antidrug antibody positivity was 38.2% in patients who were treated with anti-TNF mAbs, 6.1% with etanercept, 50.0% with rituximab, and 20.0% with tocilizumab. There was an inverse association between antidrug antibody positivity (odds ratio [OR], 0.19; 95% CI, 0.09-0.38; P <.001) directed against all biologic drugs and EULAR response at month 12. Analyzing all the visits starting at month 6 using generalized estimating equation models confirmed the inverse association between antidrug antibody positivity and EULAR response (OR, 0.35; 95% CI, 0.18-0.65; P <.001). A similar association was found for tocilizumab alone (OR, 0.18; 95% CI, 0.04-0.83; P =.03). In the multivariable analysis, antidrug antibodies, body mass index, and rheumatoid factor were independently inversely associated with response to treatment. There was a significantly higher drug concentration of anti-TNF mAbs in patients with antidrug antibody-negative vs antidrug antibody-positive status (mean difference, -9.6 [95% CI, -12.4 to -6.9] mg/L; P < 001). Drug concentrations of etanercept (mean difference, 0.70 [95% CI, 0.2-1.2] mg/L; P =.005) and adalimumab (mean difference, 1.8 [95% CI, 0.4-3.2] mg/L; P =.01) were lower in nonresponders vs responders. Methotrexate comedication at baseline was inversely associated with antidrug antibodies (OR, 0.50; 95% CI, 0.25-1.00; P =.05). Conclusions and Relevance: Results of this prospective cohort study suggest an association between antidrug antibodies and nonresponse to bDMARDs in patients with RA. Monitoring antidrug antibodies could be considered in the treatment of these patients, particularly nonresponders to biologic RA drugs

Differential co-expression framework to quantify goodness of biclusters and compare biclustering algorithms

<p>Abstract</p> <p>Background</p> <p>Biclustering is an important analysis procedure to understand the biological mechanisms from microarray gene expression data. Several algorithms have been proposed to identify biclusters, but very little effort was made to compare the performance of different algorithms on real datasets and combine the resultant biclusters into one unified ranking.</p> <p>Results</p> <p>In this paper we propose differential co-expression framework and a differential co-expression scoring function to objectively quantify quality or goodness of a bicluster of genes based on the observation that genes in a bicluster are co-expressed in the conditions belonged to the bicluster and not co-expressed in the other conditions. Furthermore, we propose a scoring function to stratify biclusters into three types of co-expression. We used the proposed scoring functions to understand the performance and behavior of the four well established biclustering algorithms on six real datasets from different domains by combining their output into one unified ranking.</p> <p>Conclusions</p> <p>Differential co-expression framework is useful to provide quantitative and objective assessment of the goodness of biclusters of co-expressed genes and performance of biclustering algorithms in identifying co-expression biclusters. It also helps to combine the biclusters output by different algorithms into one unified ranking i.e. meta-biclustering.</p

Estimation of Parent Specific DNA Copy Number in Tumors using High-Density Genotyping Arrays

Chromosomal gains and losses comprise an important type of genetic change in tumors, and can now be assayed using microarray hybridization-based experiments. Most current statistical models for DNA copy number estimate total copy number, which do not distinguish between the underlying quantities of the two inherited chromosomes. This latter information, sometimes called parent specific copy number, is important for identifying allele-specific amplifications and deletions, for quantifying normal cell contamination, and for giving a more complete molecular portrait of the tumor. We propose a stochastic segmentation model for parent-specific DNA copy number in tumor samples, and give an estimation procedure that is computationally efficient and can be applied to data from the current high density genotyping platforms. The proposed method does not require matched normal samples, and can estimate the unknown genotypes simultaneously with the parent specific copy number. The new method is used to analyze 223 glioblastoma samples from the Cancer Genome Atlas (TCGA) project, giving a more comprehensive summary of the copy number events in these samples. Detailed case studies on these samples reveal the additional insights that can be gained from an allele-specific copy number analysis, such as the quantification of fractional gains and losses, the identification of copy neutral loss of heterozygosity, and the characterization of regions of simultaneous changes of both inherited chromosomes

Clinicogenomic factors of biotherapy immunogenicity in autoimmune disease: A prospective multicohort study of the ABIRISK consortium

BACKGROUND: Biopharmaceutical products (BPs) are widely used to treat autoimmune diseases, but immunogenicity limits their efficacy for an important proportion of patients. Our knowledge of patient-related factors influencing the occurrence of antidrug antibodies (ADAs) is still limited. METHODS AND FINDINGS: The European consortium ABIRISK (Anti-Biopharmaceutical Immunization: prediction and analysis of clinical relevance to minimize the RISK) conducted a clinical and genomic multicohort prospective study of 560 patients with multiple sclerosis (MS, n = 147), rheumatoid arthritis (RA, n = 229), Crohn's disease (n = 148), or ulcerative colitis (n = 36) treated with 8 different biopharmaceuticals (etanercept, n = 84; infliximab, n = 101; adalimumab, n = 153; interferon [IFN]-beta-1a intramuscularly [IM], n = 38; IFN-beta-1a subcutaneously [SC], n = 68; IFN-beta-1b SC, n = 41; rituximab, n = 31; tocilizumab, n = 44) and followed during the first 12 months of therapy for time to ADA development. From the bioclinical data collected, we explored the relationships between patient-related factors and the occurrence of ADAs. Both baseline and time-dependent factors such as concomitant medications were analyzed using Cox proportional hazard regression models. Mean age and disease duration were 35.1 and 0.85 years, respectively, for MS; 54.2 and 3.17 years for RA; and 36.9 and 3.69 years for inflammatory bowel diseases (IBDs). In a multivariate Cox regression model including each of the clinical and genetic factors mentioned hereafter, among the clinical factors, immunosuppressants (adjusted hazard ratio [aHR] = 0.408 [95% confidence interval (CI) 0.253-0.657], p < 0.001) and antibiotics (aHR = 0.121 [0.0437-0.333], p < 0.0001) were independently negatively associated with time to ADA development, whereas infections during the study (aHR = 2.757 [1.616-4.704], p < 0.001) and tobacco smoking (aHR = 2.150 [1.319-3.503], p < 0.01) were positively associated. 351,824 Single-Nucleotide Polymorphisms (SNPs) and 38 imputed Human Leukocyte Antigen (HLA) alleles were analyzed through a genome-wide association study. We found that the HLA-DQA1*05 allele significantly increased the rate of immunogenicity (aHR = 3.9 [1.923-5.976], p < 0.0001 for the homozygotes). Among the 6 genetic variants selected at a 20% false discovery rate (FDR) threshold, the minor allele of rs10508884, which is situated in an intron of the CXCL12 gene, increased the rate of immunogenicity (aHR = 3.804 [2.139-6.764], p < 1 × 10-5 for patients homozygous for the minor allele) and was chosen for validation through a CXCL12 protein enzyme-linked immunosorbent assay (ELISA) on patient serum at baseline before therapy start. CXCL12 protein levels were higher for patients homozygous for the minor allele carrying higher ADA risk (mean: 2,693 pg/ml) than for the other genotypes (mean: 2,317 pg/ml; p = 0.014), and patients with CXCL12 levels above the median in serum were more prone to develop ADAs (aHR = 2.329 [1.106-4.90], p = 0.026). A limitation of the study is the lack of replication; therefore, other studies are required to confirm our findings. CONCLUSION: In our study, we found that immunosuppressants and antibiotics were associated with decreased risk of ADA development, whereas tobacco smoking and infections during the study were associated with increased risk. We found that the HLA-DQA1*05 allele was associated with an increased rate of immunogenicity. Moreover, our results suggest a relationship between CXCL12 production and ADA development independent of the disease, which is consistent with its known function in affinity maturation of antibodies and plasma cell survival. Our findings may help physicians in the management of patients receiving biotherapies

Pooled Analysis of Prognostic Impact of Urokinase-Type Plasminogen Activator and Its Inhibitor PAI-1 in 8377 Breast Cancer Patients

Background: Urokinase-type plasminogen activator (uPA) and its inhibitor (PAI-1) play essential roles in tumor invasion and metastasis. High levels of both uPA and PAI-1 are associated with poor prognosis in breast cancer patients. To confirm the prognostic value of uPA and PAI-1 in primary breast cancer, we reanalyzed individual patient data provided by members of the European Organization for Research and Treatment of Cancer-Receptor and Biomarker Group (EORTC-RBG). Methods: The study included 18 datasets involving 8377 breast cancer patients. During follow-up (median 79 months), 35% of the patients relapsed and 27% died. Levels of uPA and PAI-1 in tumor tissue extracts were determined by different immunoassays; values were ranked within each dataset and divided by the number of patients in that dataset to produce fractional ranks that could be compared directly across datasets. Associations of ranks of uPA and PAI-1 levels with relapse-free survival (RFS) and overall survival (OS) were analyzed by Cox multivariable regression analysis stratified by dataset, including the following traditional prognostic variables: age, menopausal status, lymph node status, tumor size, histologic grade, and steroid hormone-receptor status. All P values were two-sided. Results: Apart from lymph node status, high levels of uPA and PAI-1 were the strongest predictors of both poor RFS and poor OS in the analyses of all patients. Moreover, in both lymph node-positive and lymph node-negative patients, higher uPA and PAI-1 values were independently associated with poor RFS and poor OS. For (untreated) lymph node-negative patients in particular, uPA and PAI-1 included together showed strong prognostic ability (all P<.001). Conclusions: This pooled analysis of the EORTC-RBG datasets confirmed the strong and independent prognostic value of uPA and PAI-1 in primary breast cancer. For patients with lymph node-negative breast cancer, uPA and PAI-1 measurements in primary tumors may be especially useful for designing individualized treatment strategie