Whole-genome optical mapping of bone-marrow myeloma cells reveals association of extramedullary multiple myeloma with chromosome 1 abnormalities.

Extramedullary disease (EMM) represents a rare, aggressive and mostly resistant phenotype of multiple myeloma (MM). EMM is frequently associated with high-risk cytogenetics, but their complex genomic architecture is largely unexplored. We used whole-genome optical mapping (Saphyr, Bionano Genomics) to analyse the genomic architecture of CD138+ cells isolated from bone-marrow aspirates from an unselected cohort of newly diagnosed patients with EMM (n = 4) and intramedullary MM (n = 7). Large intrachromosomal rearrangements (> 5 Mbp) within chromosome 1 were detected in all EMM samples. These rearrangements, predominantly deletions with/without inversions, encompassed hundreds of genes and led to changes in the gene copy number on large regions of chromosome 1. Compared with intramedullary MM, EMM was characterised by more deletions (size range of 500 bp-50 kbp) and fewer interchromosomal translocations, and two EMM samples had copy number loss in the 17p13 region. Widespread genomic heterogeneity and novel aberrations in the high-risk IGH/IGK/IGL, 8q24 and 13q14 regions were detected in individual patients but were not specific to EMM/MM. Our pilot study revealed an association of chromosome 1 abnormalities in bone marrow myeloma cells with extramedullary progression. Optical mapping showed the potential for refining the complex genomic architecture in MM and its phenotypes.

Determining optical mapping errors by simulations.

MOTIVATION: Optical mapping is a complementary technology to traditional DNA sequencing technologies, such as next-generation sequencing (NGS). It provides genome-wide, high-resolution restriction maps from single, stained molecules of DNA. It can be used to detect large and small structural variants, copy number variations and complex rearrangements. Optical mapping is affected by different kinds of errors in comparison with traditional DNA sequencing technologies. It is important to understand the source of these errors and how they affect the obtained data. This article proposes a novel approach to modeling errors in the data obtained from the Bionano Genomics Inc. Saphyr system with Direct Label and Stain (DLS) chemistry. Some studies have already addressed this issue for older instruments with nicking enzymes, but we are unaware of a study that addresses this new system. RESULTS: The main result is a framework for studying errors in the data obtained from the Saphyr instrument with DLS chemistry. The framework's main component is a simulation that computes how major sources of errors for this instrument (a false site, a missing site and resolution errors) affect the distribution of fragment lengths in optical maps. The simulation is parametrized by variables describing these errors and we are using a differential evolution algorithm to evaluate parameters that best fit the data from the instrument. Results of the experiments manifest that this approach can be used to study errors in the optical mapping data analysis. AVAILABILITY AND IMPLEMENTATION: Source codes supporting the presented results are available at: https://github.com/mvasinek/olgen-om-error-prediction. The data underlying this article are available on the Bionano Genomics Inc. website, at: https://bionanogenomics.com/library/datasets/. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Excellent option for mass testing during the SARS-CoV-2 pandemic: painless self-collection and direct RT-qPCR.

The early identification of asymptomatic yet infectious cases is vital to curb the 2019 coronavirus (COVID-19) pandemic and to control the disease in the post-pandemic era. In this paper, we propose a fast, inexpensive and high-throughput approach using painless nasal-swab self-collection followed by direct RT-qPCR for the sensitive PCR detection of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). This approach was validated in a large prospective cohort study of 1038 subjects, analysed simultaneously using (1) nasopharyngeal swabs obtained with the assistance of healthcare personnel and analysed by classic two-step RT-qPCR on RNA isolates and (2) nasal swabs obtained by self-collection and analysed with direct RT-qPCR. Of these subjects, 28.6% tested positive for SARS-CoV-2 using nasopharyngeal swab sampling. Our direct RT-qPCR approach for self-collected nasal swabs performed well with results similar to those of the two-step RT-qPCR on RNA isolates, achieving 0.99 positive and 0.98 negative predictive values (cycle threshold [Ct] < 37). Our research also reports on grey-zone viraemia, including samples with near-cut-off Ct values (Ct ≥ 37). In all investigated subjects (n = 20) with grey-zone viraemia, the ultra-small viral load disappeared within hours or days with no symptoms. Overall, this study underscores the importance of painless nasal-swab self-collection and direct RT-qPCR for mass testing during the SARS-CoV-2 pandemic and in the post-pandemic era.

Direct-RT-qPCR Detection of SARS-CoV-2 without RNA Extraction as Part of a COVID-19 Testing Strategy: From Sample to Result in One Hour.

Due to the lack of protective immunity in the general population and the absence of effective antivirals and vaccines, the Coronavirus disease 2019 (COVID-19) pandemic continues in some countries, with local epicentres emerging in others. Due to the great demand for effective COVID-19 testing programmes to control the spread of the disease, we have suggested such a testing programme that includes a rapid RT-qPCR approach without RNA extraction. The Direct-One-Step-RT-qPCR (DIOS-RT-qPCR) assay detects severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in less than one hour while maintaining the high sensitivity and specificity required of diagnostic tools. This optimised protocol allows for the direct use of swab transfer media (14 µL) without the need for RNA extraction, achieving comparable sensitivity to the standard method that requires the time-consuming and costly step of RNA isolation. The limit of detection for DIOS-RT-qPCR was lower than seven copies/reaction, which translates to 550 virus copies/mL of swab. The speed, ease of use and low price of this assay make it suitable for high-throughput screening programmes. The use of fast enzymes allows RT-qPCR to be performed under standard laboratory conditions within one hour, making it a potential point-of-care solution on high-speed cycling instruments. This protocol also implements the heat inactivation of SARS-CoV-2 (75 °C for 10 min), which renders samples non-infectious, enabling testing in BSL-2 facilities. Moreover, we discuss the critical steps involved in developing tests for the rapid detection of COVID-19. Implementing rapid, easy, cost-effective methods can help control the worldwide spread of the COVID-19 infection.

Hypomethylation of IL1RN and NFKB1 genes is linked to the dysbalance in IL1ß/IL-1Ra axis in female patients with type 2 diabetes mellitus.

Inflammation has received considerable attention in the pathogenesis of type 2 diabetes mellitus (T2DM). Supporting this concept, enhanced expression of interleukin (IL)-1ß and increased infiltration of macrophages are observed in pancreatic islets of patients with T2DM. Although IL-1 receptor antagonist (IL-1Ra) plays a major role in controlling of IL-1ß-mediated inflammation, its counteraction effects and epigenetic alterations in T2DM are less studied. Thus, we aimed to analyze the DNA methylation status in IL1RN, RELA (p65) and NFKB1 (p50) genes in peripheral blood mononuclear cells (PBMCs) from treated T2DM patients (n = 35) and age-/sex- matched healthy controls (n = 31). Production of IL-1ß and IL-1Ra was analyzed in plasma and supernatants from LPS-induced PBMCs. Immunomodulatory effects of IL-1ß and IL-1Ra were studied on THP-1 cells. Average DNA methylation level of IL1RN and NFKB1 gene promoters was significantly decreased in T2DM patients in comparison with healthy controls (P< 0.05), which was associated with the increased IL-1Ra (P< 0.001) and IL-1ß (P = 0.039) plasma levels in T2DM patients. Negative association between average methylation of IL1RN gene and IL-1Ra plasma levels were observed in female T2DM patients. Methylation of NFKB1 gene was negatively correlated with IL-1Ra levels in the patients and positively with IL-1ß levels in female patients. LPS-stimulated PBMCs from female patients failed to raise IL-1ß production, while the cells from healthy females increased IL-1ß production in comparison with unstimulated cells (P< 0.001). Taken together, the findings suggest that hypomethylation of IL1RN and NFKB1 gene promoters may promote the increased IL-1ß/IL-1Ra production and regulate chronic inflammation in T2DM. Further studies are necessary to elucidate the causal direction of these associations and potential role of IL-1Ra in anti-inflammatory processes in treated patients with T2DM.

Revealed heterogeneity in rheumatoid arthritis based on multivariate innate signature analysis.

OBJECTIVES: A growing body of evidence highlights the persistent activation of the innate immune system and type I interferon (IFN) signature in the pathogenesis of rheumatoid arthritis (RA) and its association with disease activity. Since the recent study revealed heterogeneity in the IFN signature in RA, we investigated for the first time the heterogeneity in innate signature in RA. METHODS: The innate gene expression signature (10 TLRs, 7 IL1/IL1R family members, and CXCL8/IL8) was assessed in peripheral blood mononuclear cells from RA patients (n=67), both with active (DAS28≥3.2, n=32) and inactive disease (DAS28<3.2, n=35), and in healthy control subjects (n=55). RESULTS: Of the 13 deregulated innate genes (TLR2, TLR3, TLR4, TLR5, TLR8, TLR10, IL1B, IL1RN, IL18, IL18R1, IL1RAP, and SIGIRR/IL1R8) associated with RA, TLR10 and IL1RAP are being reported for the first time. Multivariate analysis based on utilising patient similarity networks revealed the existence of four patient's subsets (clusters) based on different TLR8 and IL1RN expression profiles, two in active and two in inactive RA. Moreover, neural network analysis identified two main gene sets describing active RA within an activity-related innate signature (TLR1, TLR2, TLR3, TLR7, TLR8, CXCL8/IL8, IL1RN, IL18R1). When comparing active and inactive RA, upregulated TLR2, TLR4, TLR6, and TLR8 and downregulated TLR10 (P<0.04) expression was associated with the disease activity. CONCLUSIONS: Our study on the comprehensive innate gene profiling together with multivariate analysis revealed a certain heterogeneity in innate signature within RA patients. Whether the heterogeneity of RA elucidated from diversity in innate signatures may impact the disease course and treatment response deserves future investigations.

Cross-Disease Innate Gene Signature: Emerging Diversity and Abundance in RA Comparing to SLE and SSc.

Overactivation of the innate immune system together with the impaired downstream pathway of type I interferon-responding genes is a hallmark of rheumatoid arthritis (RA), systemic lupus erythematosus (SLE), and systemic sclerosis (SSc). To date, limited data on the cross-disease innate gene signature exists among those diseases. We compared therefore an innate gene signature of Toll-like receptors (TLRs), seven key members of the interleukin (IL)1/IL1R family, and CXCL8/IL8 in peripheral blood mononuclear cells from well-defined patients with active stages of RA (n = 36, DAS28 ≥ 3.2), SLE (n = 28, SLEDAI > 6), and SSc (n = 22, revised EUSTAR index > 2.25). Emerging diversity and abundance of the innate signature in RA patients were detected: RA was characterized by the upregulation of TLR3, TLR5, IL1RAP/IL1R3, IL18R1, and SIGIRR/IL1R8 when compared to SSc (P corr < 0.02) and of TLR2, TLR5, and SIGIRR/IL1R8 when compared to SLE (P corr < 0.02). Applying the association rule analysis, six rules (combinations and expression of genes describing disease) were identified for RA (most frequently included high TLR3 and/or IL1RAP/IL1R3) and three rules for SLE (low IL1RN and IL18R1) and SSc (low TLR5 and IL18R1). This first cross-disease study identified emerging heterogeneity in the innate signature of RA patients with many upregulated innate genes compared to that of SLE and SSc.

Inflammation time-axis in aseptic loosening of total knee arthroplasty: A preliminary study.

OBJECTIVE: Aseptic loosening (AL) is the most frequent long-term reason for revision of total knee arthroplasty (TKA) affecting about 15-20% patients within 20 years after the surgery. Although there is a solid body of evidence about the crucial role of inflammation in the AL pathogenesis, scared information on inflammation signature and its time-axis in tissues around TKA exists. DESIGN: The inflammation protein signatures in pseudosynovial tissues collected at revision surgery from patients with AL (AL, n = 12) and those with no clinical/radiographic signs of AL (non-AL, n = 9) were investigated by Proximity Extension Assay (PEA)-Immunoassay and immunohistochemistry. RESULTS: AL tissues had elevated levels of TNF-family members sTNFR2, TNFSF14, sFasL, sBAFF, cytokines/chemokines IL8, CCL2, IL1RA/IL36, sIL6R, and growth factors sAREG, CSF1, comparing to non-AL. High interindividual variability in protein levels was evident particularly in non-AL. Levels of sTNFR2, sBAFF, IL8, sIL6R, and MPO discriminated between AL and non-AL and were associated with the time from index surgery, suggesting the cumulative character of inflammatory osteolytic response to prosthetic byproducts. The source of elevated inflammatory molecules was macrophages and multinucleated osteoclast-like cells in AL and histiocytes and osteoclast-like cells in non-AL tissues, respectively. All proteins were present in higher levels in osteoclast-like cells than in macrophages. CONCLUSIONS: Our study revealed a differential inflammation signature between AL and non-AL stages of TKA. It also highlighted the unique patient's response to TKA in non-AL stages. Further confirmation of our preliminary results on a larger cohort is needed. Analysis of the time-axis of processes ongoing around TKA implantation may help to understand the mechanisms driving periprosthetic bone resorption needed for diagnostic/preventative strategies.

Serum protein fingerprinting by PEA immunoassay coupled with a pattern-recognition algorithms distinguishes MGUS and multiple myeloma.

Serum protein fingerprints associated with MGUS and MM and their changes in MM after autologous stem cell transplantation (MM-ASCT, day 100) remain unexplored. Using highly-sensitive Proximity Extension ImmunoAssay on 92 cancer biomarkers (Proseek Multiplex, Olink), enhanced serum levels of Adrenomedullin (ADM, Pcorr= .0004), Growth differentiation factor 15 (GDF15, Pcorr= .003), and soluble Major histocompatibility complex class I-related chain A (sMICA, Pcorr= .023), all prosurvival and chemoprotective factors for myeloma cells, were detected in MM comparing to MGUS. Comparison of MGUS and healthy subjects revealed elevation of angiogenic and antia-poptotic midkine (Pcorr= .0007) and downregulation of Transforming growth factor beta 1 (TGFB1, Pcorr= .005) in MGUS. Importantly, altered serum pattern was associated with MM-ASCT compared to paired MM at the diagnosis as well as to healthy controls, namely by upregulated B-Cell Activating Factor (sBAFF) (Pcorr< .006) and sustained elevation of other pro-tumorigenic factors. In conclusion, the serum fingerprints of MM and MM-ASCT were characteristic by elevated levels of prosurvival and chemoprotective factors for myeloma cells.

Excellent Diagnostic Characteristics for Ultrafast Gene Profiling of DEFA1-IL1B-LTF in Detection of Prosthetic Joint Infections.

The timely and exact diagnosis of prosthetic joint infection (PJI) is crucial for surgical decision-making. Intraoperatively, delivery of the result within an hour is required. Alpha-defensin lateral immunoassay of joint fluid (JF) is precise for the intraoperative exclusion of PJI; however, for patients with a limited amount of JF and/or in cases where the JF is bloody, this test is unhelpful. Important information is hidden in periprosthetic tissues that may much better reflect the current status of implant pathology. We therefore investigated the utility of the gene expression patterns of 12 candidate genes (TLR1, -2, -4, -6, and 10, DEFA1, LTF, IL1B, BPI, CRP, IFNG, and DEFB4A) previously associated with infection for detection of PJI in periprosthetic tissues of patients with total joint arthroplasty (TJA) (n = 76) reoperated for PJI (n = 38) or aseptic failure (n = 38), using the ultrafast quantitative reverse transcription-PCR (RT-PCR) Xxpress system (BJS Biotechnologies Ltd.). Advanced data-mining algorithms were applied for data analysis. For PJI, we detected elevated mRNA expression levels of DEFA1 (P < 0.0001), IL1B (P < 0.0001), LTF (P < 0.0001), TLR1 (P = 0.02), and BPI (P = 0.01) in comparison to those in tissues from aseptic cases. A feature selection algorithm revealed that the DEFA1-IL1B-LTF pattern was the most appropriate for detection/exclusion of PJI, achieving 94.5% sensitivity and 95.7% specificity, with likelihood ratios (LRs) for positive and negative results of 16.3 and 0.06, respectively. Taken together, the results show that DEFA1-IL1B-LTF gene expression detection by use of ultrafast qRT-PCR linked to an electronic calculator allows detection of patients with a high probability of PJI within 45 min after sampling. Further testing on a larger cohort of patients is needed.

A Novel Non-Immunoglobulin (non-Ig)/BCL6 Translocation in Diffuse Large B-Cell Lymphoma Involving Chromosome 10q11.21 Loci and Review on Clinical Consequences of BCL6 Rearrangements.

BCL6 rearrangements (3q27) are the most common chromosomal abnormalities in diffuse large B-cell lymphoma (DLBCL), with numerous immunoglobulin (Ig) and non-Ig genes as partners. In DLBCL, the translocations occur predominantly in the "major breakpoint region" encompassing the first noncoding exon and a part of the first intron of BCL6; few cases with "alternative breakpoint cluster" located 245-285 kb 5' BCL6 were also described. The regulatory sequences of known Ig and non-Ig partners replace the 5' untranslated region of the BCL6 in the same transcriptional orientation. Contrary to Ig/BCL6 fusions typical by high BCL6 gene expression, in non-Ig/BCL6 translocations were observed unexpectedly low BCL6 mRNA levels. From the clinical point of view, the survival rate of DLBCL patients with non-Ig partners is inferior to those with Ig/BCL6 translocations, suggesting that non-Ig/BCL6 fusion is a poor prognostic indicator. Hereby we provide comprehensive information about known non-Ig translocation partners and clinical consequences of BCL6 rearrangements in DLBCL. Moreover, we describe a novel reciprocal translocation t(3;10) in refractory patient with DLBCL with the breaking points at 5' untranslated region of BCL6 and 5' untranslated region of the RASGEF1A gene on chromosome 10q11.21 loci; this rearrangement was associated with low BCL6 and RASGEF1A gene expressions. Our patient harbouring dual chromosomal rearrangement involving BCL2 and BCL6 genes relapsed three-times and died soon; thus, further supporting the notion that non-Ig/BCL6 fusion is a poor prognostic indicator of DLBCL. There is evidence of prognostic value of BCL6 rearrangements also in rituximab era.

Correlation Network Analysis Reveals Relationships between MicroRNAs, Transcription Factor T-bet, and Deregulated Cytokine/Chemokine-Receptor Network in Pulmonary Sarcoidosis.

Sarcoidosis is an inflammatory granulomatous disease with unknown etiology driven by cytokines and chemokines. There is limited information regarding the regulation of cytokine/chemokine-receptor network in bronchoalveolar lavage (BAL) cells in pulmonary sarcoidosis, suggesting contribution of miRNAs and transcription factors. We therefore investigated gene expression of 25 inflammation-related miRNAs, 27 cytokines/chemokines/receptors, and a Th1-transcription factor T-bet in unseparated BAL cells obtained from 48 sarcoidosis patients and 14 control subjects using quantitative RT-PCR. We then examined both miRNA-mRNA expressions to enrich relevant relationships. This first study on miRNAs in sarcoid BAL cells detected deregulation of miR-146a, miR-150, miR-202, miR-204, and miR-222 expression comparing to controls. Subanalysis revealed higher number of miR-155, let-7c transcripts in progressing (n = 20) comparing to regressing (n = 28) disease as assessed by 2-year follow-up. Correlation network analysis revealed relationships between microRNAs, transcription factor T-bet, and deregulated cytokine/chemokine-receptor network in sarcoid BAL cells. Furthermore, T-bet showed more pronounced regulatory capability to sarcoidosis-associated cytokines/chemokines/receptors than miRNAs, which may function rather as "fine-tuners" of cytokine/chemokine expression. Our correlation network study implies contribution of both microRNAs and Th1-transcription factor T-bet to the regulation of cytokine/chemokine-receptor network in BAL cells in sarcoidosis. Functional studies are needed to confirm biological relevance of the obtained relationships.

C-MYC and C-FOS expression changes and cellular aspects of the photodynamic reaction with photosensitizers TMPyP and ClAlPcS2.

Photodynamic therapy (PDT) is based on the tumor-selective accumulation of photosensitizer followed by irradiation with light of an appropriate wavelength. After irradiation and in the presence of oxygen, photosensitizer induces cellular damage. The aim of this study was to evaluate effects of two photosensitizers TMPyP and ClAlPcS2 on cell lines to obtain better insight into their mechanisms of action. We determined cell viability, reactive oxygen species (ROS) generation and changes in expression levels of two important early response genes, C-MYC and C-FOS, on tumor MCF7 (human breast adenocarcinoma) and G361 (human melanoma) cell lines and non-tumor BJ cell line (human fibroblast) after photodynamic reaction with TMPyP and ClAlPcS2 as photosensitizers. In addition TMPyP and ClAlPcS2 cellular uptake and clearance and antioxidant capacity of the mentioned cell lines were investigated. We found appropriate therapeutic doses and confirmed that both tested photosensitizers are photodynamically efficient in treatment used cells in vitro. TMPyP is more efficient; it had higher ROS production and toxicity after irradiation by intermediate therapeutic doses than ClAlPcS2. We revealed that both TMPyP and ClAlPcS2-PDT increased C-FOS expression on tumor cell lines (G361 and MCF7), but not on non-tumor BJ cell line. Conversely, both TMPyP and ClAlPcS2-PDT decreased C-MYC expression on non-tumor BJ cell line but not on tumor cell lines. As first we tested these photosensitizers in such extent and we believe that it can help to better understand mechanisms of PDT and increase its efficiency and applicability.

A novel frameshift mutation in the AFG3L2 gene in a patient with spinocerebellar ataxia.

Spinocerebellar ataxia type 28 (SCA28) is an autosomal dominant neurodegenerative disorder caused by missense AFG3L2 mutations. To examine the occurrence of SCA28 in the Czech Republic, we screened 288 unrelated ataxic patients with hereditary (N = 49) and sporadic or unknown (N = 239) form of ataxia for mutations in exons 15 and 16, the AFG3L2 mutation hotspots. A single significant variant, frameshift mutation c.1958dupT leading to a premature termination codon, was identified in a patient with slowly progressive speech and gait problems starting at the age of 68 years. Neurological examination showed cerebellar ataxia, mild Parkinsonian features with predominant bradykinesia, polyneuropathy of the lower limbs, and cognitive decline. However, other common SCA28 features like pyramidal tract signs (lower limb hyperreflexia, positive Babinski sign), ophthalmoparesis or ptosis were absent. The mutation was also found in a patient's unaffected daughter in whom a targeted examination at 53 years of age revealed mild imbalance signs. RNA analysis showed a decreased ratio of the transcript from the mutated AFG3L2 allele relative to the normal transcript in the peripheral lymphocytes of both patients. The ratio was increased by puromycin treatment, indicating that the mutated transcript can be degraded via nonsense-mediated RNA decay. The causal link between the mutation and the phenotype of the patient is currently unclear but a pathogenic mechanism based on AFG3L2 haploinsufficiency rather than the usual dominant-negative effect of missense AFG3L2 mutations reported in SCA28, cannot be excluded.

Activated phenotype of circulating neutrophils in familial Mediterranean fever.

Familial Mediterranean fever (FMF) is autoinflammatory disorder, characterized by MEFV gene mutations and recurrent episodes of fever and serosal or synovial inflammation. Neutrophils are the predominant effector cells of acute inflammatory attacks in FMF; however pathogenic role and molecular phenotype of these cells remain largely unknown. To gain insight into the processes that contribute to the self-directed autoinflammation we characterized expression of a spectrum of genes involved in regulation of inflammation in unstimulated and LPS-activated neutrophils from FMF patients. Expression of 12 candidate immune genes encoding for inflammation-related molecules was assessed by quantitative RT-PCR in freshly isolated and LPS-stimulated peripheral polymorphonuclear neutrophils from fifteen FMF patients in attack-free period and ten healthy volunteers as controls. The relative expression was calculated using the second derivative method; the target gene expression was normalized to the expression of RPL32 gene. FMF neutrophils were characterized by up-regulated baseline gene expression of c-FOS (9.5-fold, p < 0.05), IL-8 (12-fold, p < 0.05), MMP9 (8-fold, p < 0.01), TLR2 (7-fold, p < 0.05) compared to the neutrophils from control subjects, a trend was also evident towards increased caspase-1 expression (3-fold, p = 0.09). Discriminant analysis clustered the patient and control subjects into two distinct groups (Wilks's lambda = 0.165, p = 0.042). Further, LPS-induced alterations of expression profiles were shared between FMF and healthy neutrophils, the profile consisting namely of up-regulated IL-1ß, TLR4, IL-8, and TNFAIP6 transcripts. Present study demonstrates distinct expression patterns of pre-activated neutrophils during attack-free period of FMF when compared to neutrophils from healthy controls. Furthermore, our data emphasize the importance of host-derived ligands in activation of FMF neutrophils.

PSMB2 and RPL32 are suitable denominators to normalize gene expression profiles in bronchoalveolar cells.

BACKGROUND: For accuracy of quantitative reverse transcriptase-polymerase chain reaction (qRT-PCR), normalisation with suitable reference genes is required. To date, no reference genes have been validated for expression studies of bronchoalveolar (BAL) cells. The aims of this study were to identify gene(s) with stable mRNA expression in BAL cells irrespective of gender, smoking, BAL cellular composition, lung pathology, treatment; and to assess the influence of reference genes on target gene expression data. RESULTS: The mRNA expression of ten housekeeping genes (ACTB, ARF1, CANX, G6PD, GAPDH, GPS1, GNB2L1, PSMB2, PSMD2, RPL32) was investigated by qRT-PCR in BAL cells from 71 subjects across a spectrum of lung diseases. The analyses were validated in an independent BAL cohort from 63 sarcoidosis patients and 17 control subjects. A second derivative method was used to calculate expression values (CTt); an equivalence test, applets BestKeeper, geNorm and NormFinder were applied to investigate gene expression stability. Of the investigated genes, PSMB2 (CTt +/- SD, 23.66 +/- 0.86) and RPL32 (18.65 +/- 0.92) were the most stable; both were constantly expressed in BAL samples from parallel investigated cohorts irrespective of evaluated variables. Finally, to demonstrate effect of traditional (ACTB/GAPDH) and novel (PSMB2/RPL32) reference genes as denominators, expression of two cytokines known associated with sarcoidosis was investigated in sarcoid BAL cells. While normalization with PSMB2/RPL32 resulted in elevated IFNG mRNA expression (p = 0.004); no change was observed using GAPDH/ACTB (p > 0.05). CCL2 mRNA up-regulation was observed only when PSMB2/RPL32 were used as denominators (p < 0.03). CONCLUSION: PSMB2 and RPL32 are, therefore, suitable reference genes to normalize qRT-PCR in BAL cells in sarcoidosis, and other interstitial lung disease.