Systematic discovery of gene fusions in pediatric cancer by integrating RNA-seq and WGS.

BACKGROUND: Gene fusions are important cancer drivers in pediatric cancer and their accurate detection is essential for diagnosis and treatment. Clinical decision-making requires high confidence and precision of detection. Recent developments show RNA sequencing (RNA-seq) is promising for genome-wide detection of fusion products but hindered by many false positives that require extensive manual curation and impede discovery of pathogenic fusions. METHODS: We developed Fusion-sq to overcome existing disadvantages of detecting gene fusions. Fusion-sq integrates and "fuses" evidence from RNA-seq and whole genome sequencing (WGS) using intron-exon gene structure to identify tumor-specific protein coding gene fusions. Fusion-sq was then applied to the data generated from a pediatric pan-cancer cohort of 128 patients by WGS and RNA sequencing. RESULTS: In a pediatric pan-cancer cohort of 128 patients, we identified 155 high confidence tumor-specific gene fusions and their underlying structural variants (SVs). This includes all clinically relevant fusions known to be present in this cohort (30 patients). Fusion-sq distinguishes healthy-occurring from tumor-specific fusions and resolves fusions in amplified regions and copy number unstable genomes. A high gene fusion burden is associated with copy number instability. We identified 27 potentially pathogenic fusions involving oncogenes or tumor-suppressor genes characterized by underlying SVs, in some cases leading to expression changes indicative of activating or disruptive effects. CONCLUSIONS: Our results indicate how clinically relevant and potentially pathogenic gene fusions can be identified and their functional effects investigated by combining WGS and RNA-seq. Integrating RNA fusion predictions with underlying SVs advances fusion detection beyond extensive manual filtering. Taken together, we developed a method for identifying candidate gene fusions that is suitable for precision oncology applications. Our method provides multi-omics evidence for assessing the pathogenicity of tumor-specific gene fusions for future clinical decision making.

Mutational analysis using Sanger and next generation sequencing in sporadic spindle cell hemangiomas: A study of 19 cases.

Spindle cell hemangioma (SCH) is a distinct vascular soft-tissue lesion characterized by cavernous blood vessels and a spindle cell component mainly occurring in the distal extremities of young adults. The majority of cases harbor heterozygous mutations in IDH1/2 sporadically or rarely in association with Maffucci syndrome. However, based on mosaicism and accordingly a low percentage of lesional cells harboring a mutant allele, detection can be challenging. We tested 19 sporadic SCHs by Sanger sequencing, multiplex ligation-dependent probe amplification (MLPA), conventional next generation sequencing (NGS), and NGS using a single molecule molecular inversion probes (smMIP)-based library preparation to compare their diagnostic value. Out of 10 cases tested by Sanger sequencing and 2 analyzed using MLPA, 4 and 1, respectively, revealed a mutation in IDH1 (p.R132C). The 7 remaining negative cases and additional 6 cases were investigated using smMIP/NGS, showing hot spot mutations in IDH1 (p.R132C) (8 cases) and IDH2 (3 cases; twice p.R172S and once p.R172G, respectively). One case was negative. Owing to insufficient DNA quality and insufficient coverage, 2 cases were excluded. In total, in 16 out of 17 cases successfully tested, an IDH1/2 mutation was found. Given that IDH1/2 mutations were absent in 161 other vascular lesions tested by smMIP/NGS, the mutation can be considered as highly specific for SCH.

Seasonal changes in gene expression represent cell-type composition in whole blood.

Seasonal patterns in behavior and biological parameters are widespread. Here, we examined seasonal changes in whole blood gene expression profiles of 233 healthy subjects. Using weighted gene co-expression network analysis, we identified three co-expression modules showing circannual patterns. Enrichment analysis suggested that this signal stems primarily from red blood cells and blood platelets. Indeed, a large clinical database with 51 142 observations of blood cell counts over 3 years confirmed a corresponding seasonal pattern of counts of red blood cells, reticulocytes and platelets. We found no direct evidence that these changes are linked to genes known to be key players in regulating immune function or circadian rhythm. It is likely, however, that these seasonal changes in cell counts and gene expression profiles in whole blood represent biological and clinical relevant phenomena. Moreover, our findings highlight possible confounding factors relevant to the study of gene expression profiles in subjects collected at geographical locations with disparaging seasonality patterns.

Common variants conferring risk of schizophrenia.

Schizophrenia is a complex disorder, caused by both genetic and environmental factors and their interactions. Research on pathogenesis has traditionally focused on neurotransmitter systems in the brain, particularly those involving dopamine. Schizophrenia has been considered a separate disease for over a century, but in the absence of clear biological markers, diagnosis has historically been based on signs and symptoms. A fundamental message emerging from genome-wide association studies of copy number variations (CNVs) associated with the disease is that its genetic basis does not necessarily conform to classical nosological disease boundaries. Certain CNVs confer not only high relative risk of schizophrenia but also of other psychiatric disorders. The structural variations associated with schizophrenia can involve several genes and the phenotypic syndromes, or the 'genomic disorders', have not yet been characterized. Single nucleotide polymorphism (SNP)-based genome-wide association studies with the potential to implicate individual genes in complex diseases may reveal underlying biological pathways. Here we combined SNP data from several large genome-wide scans and followed up the most significant association signals. We found significant association with several markers spanning the major histocompatibility complex (MHC) region on chromosome 6p21.3-22.1, a marker located upstream of the neurogranin gene (NRGN) on 11q24.2 and a marker in intron four of transcription factor 4 (TCF4) on 18q21.2. Our findings implicating the MHC region are consistent with an immune component to schizophrenia risk, whereas the association with NRGN and TCF4 points to perturbation of pathways involved in brain development, memory and cognition.

Investigation of the involvement of MIR185 and its target genes in the development of schizophrenia.

BACKGROUND: Schizophrenia is a complex neuropsychiatric disorder of unclear etiology. The strongest known genetic risk factor is the 22q11.2 microdeletion. Research has yet to confirm which genes within the deletion region are implicated in schizophrenia. The minimal 1.5 megabase deletion contains MIR185, which encodes microRNA 185. METHODS: We determined miR-185 expression in embryonic and adult mouse brains. Common and rare variants at this locus were then investigated using a human genetics approach. First, we performed gene-based analyses for MIR185 common variants and target genes using Psychiatric Genomics Consortium genome-wide association data. Second, MIR185 was resequenced in German patients (n = 1000) and controls (n = 500). We followed up promising variants by genotyping an additional European sample (patients, n = 3598; controls, n = 4082). RESULTS: In situ hybridization in mice revealed miR-185 expression in brain regions implicated in schizophrenia. Gene-based tests revealed association between common variants in 3 MIR185 target genes (ATAT1, SH3PXD2A, NTRK3) and schizophrenia. Further analyses in mice revealed overlapping expression patterns for these target genes and miR-185. Resequencing identified 2 rare patient-specific novel variants flanking MIR185. However, follow-up genotyping provided no further evidence of their involvement in schizophrenia. LIMITATIONS: Power to detect rare variant associations was limited. CONCLUSION: Human genetic analyses generated no evidence of the involvement of MIR185 in schizophrenia. However, the expression patterns of miR-185 and its target genes in mice, and the genetic association results for the 3 target genes, suggest that further research into the involvement of miR-185 and its downstream pathways in schizophrenia is warranted.

Myosin IXB variant increases the risk of celiac disease and points toward a primary intestinal barrier defect.

Celiac disease is probably the best-understood immune-related disorder. The disease presents in the small intestine and results from the interplay between multiple genes and gluten, the triggering environmental factor. Although HLA class II genes explain 40% of the heritable risk, non-HLA genes accounting for most of the familial clustering have not yet been identified. Here we report significant and replicable association (P = 2.1 x 10(-6)) to a common variant located in intron 28 of the gene myosin IXB (MYO9B), which encodes an unconventional myosin molecule that has a role in actin remodeling of epithelial enterocytes. Individuals homozygous with respect to the at-risk allele have a 2.3-times higher risk of celiac disease (P = 1.55 x 10(-5)). This result is suggestive of a primary impairment of the intestinal barrier in the etiology of celiac disease, which may explain why immunogenic gluten peptides are able to pass through the epithelial barrier.

Common variants at VRK2 and TCF4 conferring risk of schizophrenia.

Common sequence variants have recently joined rare structural polymorphisms as genetic factors with strong evidence for association with schizophrenia. Here we extend our previous genome-wide association study and meta-analysis (totalling 7 946 cases and 19 036 controls) by examining an expanded set of variants using an enlarged follow-up sample (up to 10 260 cases and 23 500 controls). In addition to previously reported alleles in the major histocompatibility complex region, near neurogranin (NRGN) and in an intron of transcription factor 4 (TCF4), we find two novel variants showing genome-wide significant association: rs2312147[C], upstream of vaccinia-related kinase 2 (VRK2) [odds ratio (OR) = 1.09, P = 1.9 × 10(-9)] and rs4309482[A], between coiled-coiled domain containing 68 (CCDC68) and TCF4, about 400 kb from the previously described risk allele, but not accounted for by its association (OR = 1.09, P = 7.8 × 10(-9)).

Genetic variation and the risk of haloperidol-related parkinsonism in elderly patients: a candidate gene approach.

Factors that influence the variation in occurrence of antipsychotic-related parkinsonism in elderly have not been well elucidated. The aim of this study was to investigate whether previous identified and studied genetic polymorphisms at DRD2, ANKK1, DRD3, HTR2A, HTR2C, RGS2, COMT, and BDNF genes are associated with antipsychotic-related parkinsonism in elderly patients.This cross-sectional study included 150 inpatients aged 65 years and older who were treated with haloperidol. Parkinsonism assessed by the Simpson Angus Scale was present in 46% of the included patients. The investigated predictors were polymorphisms in DRD2 (141CIns/Del and C957T), ANNK1 (TaqIA), DRD3 (Ser9Gly), HTR2A (-1438G>A and His452Tyr), HTR2C (Cys23Ser and -759C/T), RGS2 (+2971C>G), COMT (G158A), and BDNF (Val66Met). Frequencies of the -759 T allele of the HTR2C gene and the 158A allele of the COMT gene were significantly higher in patients without antipsychotic-induced parkinsonism (AIP) (nominal P = 0.03 and P = 0.02, respectively). -759 T allele carriership in females was associated with a lower risk of AIP (adjusted odds ratio, 0.31; 95% confidence interval, 0.11-0.85). The decrease in risk of AIP in carriers of the COMT 158A allele did not reach statistical significance. No significant associations were found between AIP and the remaining selected polymorphisms.Although validation is needed, this study suggests that carriership of the -759 T allele of the HTR2C gene in females may be protective against development of parkinsonism in elderly patients during treatment with haloperidol.

Common inversion polymorphism at 17q21.31 affects expression of multiple genes in tissue-specific manner.

BACKGROUND: Chromosome 17q21.31 contains a common inversion polymorphism of approximately 900 kb in populations with European ancestry. Two divergent MAPT haplotypes, H1 and H2 are described with distinct linkage disequilibrium patterns across the region reflecting the inversion status at this locus. The MAPT H1 haplotype has been associated with progressive supranuclear palsy, corticobasal degeneration, Parkinson's disease and Alzheimer's disease, while the H2 is linked to recurrent deletion events associated with the 17q21.31 microdeletion syndrome, a disease characterized by developmental delay and learning disability. RESULTS: In this study, we investigate the effect of the inversion on the expression of genes in the 17q21.31 region. We find the expression of several genes in and at the borders of the inversion to be affected; specific either to whole blood or different regions of the human brain. The H1 haplotype was found to be associated with an increased expression of LRRC37A4, PLEKH1M and MAPT. In contrast, a decreased expression of MGC57346, LRRC37A and CRHR1 was associated with H1. CONCLUSIONS: Studies thus far have focused on the expression of MAPT in the inversion region. However, our results show that the inversion status affects expression of other genes in the 17q21.31 region as well. Given the link between the inversion status and different neurological diseases, these genes may also be involved in disease pathology, possibly in a tissue-specific manner.

Genetic analysis of DNA methylation and gene expression levels in whole blood of healthy human subjects.

BACKGROUND: The predominant model for regulation of gene expression through DNA methylation is an inverse association in which increased methylation results in decreased gene expression levels. However, recent studies suggest that the relationship between genetic variation, DNA methylation and expression is more complex. RESULTS: Systems genetic approaches for examining relationships between gene expression and methylation array data were used to find both negative and positive associations between these levels. A weighted correlation network analysis revealed that i) both transcriptome and methylome are organized in modules, ii) co-expression modules are generally not preserved in the methylation data and vice-versa, and iii) highly significant correlations exist between co-expression and co-methylation modules, suggesting the existence of factors that affect expression and methylation of different modules (i.e., trans effects at the level of modules). We observed that methylation probes associated with expression in cis were more likely to be located outside CpG islands, whereas specificity for CpG island shores was present when methylation, associated with expression, was under local genetic control. A structural equation model based analysis found strong support in particular for a traditional causal model in which gene expression is regulated by genetic variation via DNA methylation instead of gene expression affecting DNA methylation levels. CONCLUSIONS: Our results provide new insights into the complex mechanisms between genetic markers, epigenetic mechanisms and gene expression. We find strong support for the classical model of genetic variants regulating methylation, which in turn regulates gene expression. Moreover we show that, although the methylation and expression modules differ, they are highly correlated.

A large genome scan for rare CNVs in amyotrophic lateral sclerosis.

Amyotrophic lateral sclerosis (ALS) is a fatal neurodegenerative disease selectively affecting motor neurons in the brain and spinal cord. Recent genome-wide association studies (GWASs) have identified several common variants which increase disease susceptibility. In contrast, rare copy-number variants (CNVs), which have been associated with several neuropsychiatric traits, have not been studied for ALS in well-powered study populations. To examine the role of rare CNVs in ALS susceptibility, we conducted a CNV association study including over 19,000 individuals. In a genome-wide screen of 1875 cases and 8731 controls, we did not find evidence for a difference in global CNV burden between cases and controls. In our association analyses, we identified two loci that met our criteria for follow-up: the DPP6 locus (OR = 3.59, P = 6.6 × 10(-3)), which has already been implicated in ALS pathogenesis, and the 15q11.2 locus, containing NIPA1 (OR = 12.46, P = 9.3 × 10(-5)), the gene causing hereditary spastic paraparesis type 6 (HSP 6). We tested these loci in a replication cohort of 2559 cases and 5887 controls. Again, results were suggestive of association, but did not meet our criteria for independent replication: DPP6 locus: OR = 1.92, P = 0.097, pooled results: OR = 2.64, P = 1.4 × 10(-3); NIPA1: OR = 3.23, P = 0.041, pooled results: OR = 6.20, P = 2.2 × 10(-5)). Our results highlight DPP6 and NIPA1 as candidates for more in-depth studies. Unlike other complex neurological and psychiatric traits, rare CNVs with high effect size do not play a major role in ALS pathogenesis.

A large replication study and meta-analysis in European samples provides further support for association of AHI1 markers with schizophrenia.

The Abelson helper integration site 1 (AHI1) gene locus on chromosome 6q23 is among a group of candidate loci for schizophrenia susceptibility that were initially identified by linkage followed by linkage disequilibrium mapping, and subsequent replication of the association in an independent sample. Here, we present results of a replication study of AHI1 locus markers, previously implicated in schizophrenia, in a large European sample (in total 3907 affected and 7429 controls). Furthermore, we perform a meta-analysis of the implicated markers in 4496 affected and 18,920 controls. Both the replication study of new samples and the meta-analysis show evidence for significant overrepresentation of all tested alleles in patients compared with controls (meta-analysis; P = 8.2 x 10(-5)-1.7 x 10(-3), common OR = 1.09-1.11). The region contains two genes, AHI1 and C6orf217, and both genes-as well as the neighbouring phosphodiesterase 7B (PDE7B)-may be considered candidates for involvement in the genetic aetiology of schizophrenia.

Molecular Characterization Reveals Subclasses of 1q Gain in Intermediate Risk Wilms Tumors.

Chromosomal alterations have recurrently been identified in Wilms tumors (WTs) and some are associated with poor prognosis. Gain of 1q (1q+) is of special interest given its high prevalence and is currently actively studied for its prognostic value. However, the underlying mutational mechanisms and functional effects remain unknown. In a national unbiased cohort of 30 primary WTs, we integrated somatic SNVs, CNs and SVs with expression data and distinguished four clusters characterized by affected biological processes: muscle differentiation, immune system, kidney development and proliferation. Combined genome-wide CN and SV profiles showed that tumors profoundly differ in both their types of 1q+ and genomic stability and can be grouped into WTs with co-occurring 1p-/1q+, multiple chromosomal gains or CN neutral tumors. We identified 1q+ in eight tumors that differ in mutational mechanisms, subsequent rearrangements and genomic contexts. Moreover, 1q+ tumors were present in all four expression clusters reflecting activation of various biological processes, and individual tumors overexpress different genes on 1q. In conclusion, by integrating CNs, SVs and gene expression, we identified subgroups of 1q+ tumors reflecting differences in the functional effect of 1q gain, indicating that expression data is likely needed for further risk stratification of 1q+ WTs.

Improved Gene Fusion Detection in Childhood Cancer Diagnostics Using RNA Sequencing.

PURPOSE: Gene fusions play a significant role in cancer etiology, making their detection crucial for accurate diagnosis, prognosis, and determining therapeutic targets. Current diagnostic methods largely focus on either targeted or low-resolution genome-wide techniques, which may be unable to capture rare events or both fusion partners. We investigate if RNA sequencing can overcome current limitations with traditional diagnostic techniques to identify gene fusion events. METHODS: We first performed RNA sequencing on a validation cohort of 24 samples with a known gene fusion event, after which a prospective pan-pediatric cancer cohort (n = 244) was tested by RNA sequencing in parallel to existing diagnostic procedures. This cohort included hematologic malignancies, tumors of the CNS, solid tumors, and suspected neoplastic samples. All samples were processed in the routine diagnostic workflow and analyzed for gene fusions using standard-of-care methods and RNA sequencing. RESULTS: We identified a clinically relevant gene fusion in 83 of 244 cases in the prospective cohort. Sixty fusions were detected by both routine diagnostic techniques and RNA sequencing, and one fusion was detected only in routine diagnostics, but an additional 24 fusions were detected solely by RNA sequencing. RNA sequencing, therefore, increased the diagnostic yield by 38%-39%. In addition, RNA sequencing identified both gene partners involved in the gene fusion, in contrast to most routine techniques. For two patients, the newly identified fusion by RNA sequencing resulted in treatment with targeted agents. CONCLUSION: We show that RNA sequencing is sufficiently robust for gene fusion detection in routine diagnostics of childhood cancers and can make a difference in treatment decisions.

Recurrent CNVs disrupt three candidate genes in schizophrenia patients.

Schizophrenia is a severe psychiatric disease with complex etiology, affecting approximately 1% of the general population. Most genetics studies so far have focused on disease association with common genetic variation, such as single-nucleotide polymorphisms (SNPs), but it has recently become apparent that large-scale genomic copy-number variants (CNVs) are involved in disease development as well. To assess the role of rare CNVs in schizophrenia, we screened 54 patients with deficit schizophrenia using Affymetrix's GeneChip 250K SNP arrays. We identified 90 CNVs in total, 77 of which have been reported previously in unaffected control cohorts. Among the genes disrupted by the remaining rare CNVs are MYT1L, CTNND2, NRXN1, and ASTN2, genes that play an important role in neuronal functioning but--except for NRXN1--have not been associated with schizophrenia before. We studied the occurrence of CNVs at these four loci in an additional cohort of 752 patients and 706 normal controls from The Netherlands. We identified eight additional CNVs, of which the four that affect coding sequences were found only in the patient cohort. Our study supports a role for rare CNVs in schizophrenia susceptibility and identifies at least three candidate genes for this complex disorder.

Gene expression profiling in C57BL/6J and A/J mouse inbred strains reveals gene networks specific for brain regions independent of genetic background.

BACKGROUND: We performed gene expression profiling of the amygdala and hippocampus taken from inbred mouse strains C57BL/6J and A/J. The selected brain areas are implicated in neurobehavioral traits while these mouse strains are known to differ widely in behavior. Consequently, we hypothesized that comparing gene expression profiles for specific brain regions in these strains might provide insight into the molecular mechanisms of human neuropsychiatric traits. We performed a whole-genome gene expression experiment and applied a systems biology approach using weighted gene co-expression network analysis. RESULTS: We were able to identify modules of co-expressed genes that distinguish a strain or brain region. Analysis of the networks that are most informative for hippocampus and amygdala revealed enrichment in neurologically, genetically and psychologically related pathways. Close examination of the strain-specific gene expression profiles, however, revealed no functional relevance but a significant enrichment of single nucleotide polymorphisms in the probe sequences used for array hybridization. This artifact was not observed for the modules of co-expressed genes that distinguish amygdala and hippocampus. CONCLUSIONS: The brain-region specific modules were found to be independent of genetic background and are therefore likely to represent biologically relevant molecular networks that can be studied to complement our knowledge about pathways in neuropsychiatric disease.

Association study of NRG1, DTNBP1, RGS4, G72/G30, and PIP5K2A with schizophrenia and symptom severity in a Hungarian sample.

Genetic association studies have yielded extensive but frequently inconclusive data about genetic risk factors for schizophrenia. Clinical and genetic heterogeneity are possible factors explaining the inconsistent findings. The objective of this study was to test the association of commonly incriminated candidate genes with two clinically divergent subgroups, non-deficit (SZ-ND) and deficit-schizophrenia (SZ-D), and symptom severity, in order to test for replication of previously reported results. A homogeneous sample of 280 schizophrenia patients and 230 healthy controls of Hungarian, Caucasian descent were genotyped for polymorphisms in schizophrenia candidate genes NRG1, DTNBP1, RGS4, G72/G30, and PIP5K2A. Patients were divided into the diagnostic subgroups of SZ-ND and SZ-D using the Schedule for Deficit Syndrome (SDS), and assessed clinically by the Positive and Negative Symptom Scale (PANSS). SNP8NRG241930 in NRG1 and rs1011313 in DTNBP1 were associated with SZ-ND (P = 0.04 and 0.03, respectively). Polymorphisms in RGS4, G72/G30, and PIP5K2A were neither associated with SZ-ND nor with SZ-D. SNP8NRG241930 showed association with the PANSS cognitive and hostility/excitability factors, rs1011313 with the negative factor and SDS total score, and rs10917670 in RGS4 was associated with the depression factor. Although these results replicate earlier findings about the genetic background of SZ-ND and SZ-D only partially, our data seem to confirm previously reported association of NRG1 with schizophrenia without prominent negative symptoms. It was possible to detect associations of small-to-medium effect size between the investigated candidate genes and symptom severity. Such studies have the potential to unravel the possible connection between genetic and clinical heterogeneity in schizophrenia.

Amplicon-Based Targeted Next-Generation Sequencing of Formalin-Fixed, Paraffin-Embedded Tissue.

Next-generation sequencing (NGS) is rapidly becoming the method of choice for mutation analysis in both research and diagnostics. The benefit of targeted NGS compared to whole-genome and whole-exome sequencing is that smaller amounts of input material can be used as well as qualitatively suboptimal tissue samples, like formalin-fixed, paraffin-embedded archival tissue.Here, we describe the protocol for targeted next-generation sequencing using the Ion Torrent PGM platform in combination with Ion Ampliseq NGS gene panels for formalin-fixed, paraffin-embedded tissues. Both the manual and the automated workflow are described as well as the bioinformatics for data analysis.

E-cadherin loss induces targetable autocrine activation of growth factor signalling in lobular breast cancer.

Despite the fact that loss of E-cadherin is causal to the development and progression of invasive lobular carcinoma (ILC), options to treat this major breast cancer subtype are limited if tumours develop resistance to anti-oestrogen treatment regimens. This study aimed to identify clinically targetable pathways that are aberrantly active downstream of E-cadherin loss in ILC. Using a combination of reverse-phase protein array (RPPA) analyses, mRNA sequencing, conditioned medium growth assays and CRISPR/Cas9-based knock-out experiments, we demonstrate that E-cadherin loss causes increased responsiveness to autocrine growth factor receptor (GFR)-dependent activation of phosphatidylinositol-4,5-bisphosphate 3-kinase (PI3K)/Akt signalling. Autocrine activation of GFR signalling and its downstream PI3K/Akt hub was independent of oncogenic mutations in PIK3CA, AKT1 or PTEN. Analyses of human ILC samples confirmed growth factor production and pathway activity. Pharmacological inhibition of Akt using AZD5363 or MK2206 resulted in robust inhibition of cell growth and survival of ILC cells, and impeded tumour growth in a mouse ILC model. Because E-cadherin loss evokes hypersensitisation of PI3K/Akt activation independent of oncogenic mutations in this pathway, we propose clinical intervention of PI3K/Akt in ILC based on functional E-cadherin inactivation, irrespective of activating pathway mutations.

Neutrophil recruitment and barrier impairment in celiac disease: a genomic study.

BACKGROUND & AIMS: Celiac disease is an enteropathy featuring villous atrophy, crypt hyperplasia, and lymphocytosis. Tissue remodeling is driven by an inflammatory reaction to gluten in genetically susceptible individuals. The adaptive pathway is considered the major immune response but recent evidence has indicated the involvement of innate immunity as well. To assess the contribution of either immune response we performed global gene expression profiling of the regenerating mucosa. METHODS: Microarray hybridizations were performed with biopsy samples from 13 untreated patients, 31 patients on a gluten-free diet in various stages of remission, and 21 controls. Additional data were generated using low-density array and conventional quantitative reverse-transcription polymerase chain reaction, and immunohistochemistry. RESULTS: A total of 108 differentially expressed immune-related genes were identified (50 innate, 43 adaptive, 9 both innate/adaptive, and 6 immunoregulatory). Expression levels showed a gradual change as opposed to the discrete histological transitions. In addition to details provided on the adaptive and innate immune pathways used, we observed a chronic recruitment of activated neutrophils. Neutrophil involvement was unabated in otherwise completely normalized remission patients. CONCLUSIONS: We observed a contribution of both the innate and adaptive immune response in celiac disease pathogenesis. The discrepancy between the histological classification and the observed incremental change in immune-gene expression may have consequences for current diagnostic inclusion criteria. Enhanced neutrophil infiltration in both active and remission patients points to a genetic impairment of the intestinal barrier that may contribute to the cause rather than the consequence of celiac disease.