Filaggrin loss-of-function variants are associated with atopic dermatitis phenotypes in a diverse, early life prospective cohort

Loss-of-Function (LoF) variants in the filaggrin (FLG) gene are the strongest known genetic risk factor for atopic dermatitis (AD), but the impact of these variants on AD outcomes is poorly understood. We comprehensively identified genetic variants through targeted region sequencing of FLG in children (n = 438) participating in the Mechanisms of Progression of Atopic Dermatitis to Asthma in Children (MPAACH) cohort. Twenty FLG LoF variants were identified, including one novel variant and nine variants not previously associated with AD. FLG LoF variants were found in 13.6% of the cohort. Among these children, the presence of one or more FLG LoF variants was associated with moderate/severe AD (odds ratio (OR) = 2.00 (95% CI, 1.23-3.68) compared to those with mild AD. Children with FLG LoF variants had a higher SCORAD (SCORing for Atopic Dermatitis (SCORAD); P = 0.012) and higher likelihood of food allergy within the first 2.5 years of life (OR = 2.81, 1.50-5.26). LoF variants were associated with higher transepidermal Water Loss (TEWL) in both lesional (P = 0.018) and non-lesional skin (P = 0.015). Collectively, our study identifies established and novel AD-associated FLG LoF variants and associates FLG LoF with higher TEWL in lesional and non-lesional skin.

VExD: a curated resource for human gene expression alterations following viral infection

Much of the host antiviral response is mediated through changes to host gene expression levels. Likewise, viruses induce changes to host gene expression levels in order to promote the viral life cycle and evade the host immune system. However, there is no resource that specifically collects human gene expression levels pre- and post-virus infection. Further, public gene expression repositories do not contain enough specialized metadata to easily find relevant experiments. Here, we present the Virus Expression Database (VExD), a freely available website and database, that collects human gene expression datasets in response to viral infection. VExD contains ∼8,000 uniformly processed samples obtained from 289 studies examining 51 distinct human viruses. We show that the VExD processing pipeline captures known antiviral responses in the form of interferon-stimulated genes. We further show that the datasets collected in VExD can be used to quickly identify supporting data for experiments performed in human cells or model organisms. VExD is freely available at https://vexd.cchmc.org/.

Single cell transcriptomic analysis of HPV16-infected epithelium identifies a keratinocyte subpopulation implicated in cancer

Persistent HPV16 infection is a major cause of the global cancer burden. The viral life cycle is dependent on the differentiation program of stratified squamous epithelium, but the landscape of keratinocyte subpopulations which support distinct phases of the viral life cycle has yet to be elucidated. Here, single cell RNA sequencing of HPV16 infected compared to uninfected organoids identifies twelve distinct keratinocyte populations, with a subset mapped to reconstruct their respective 3D geography in stratified squamous epithelium. Instead of conventional terminally differentiated cells, an HPV-reprogrammed keratinocyte subpopulation (HIDDEN cells) forms the surface compartment and requires overexpression of the ELF3/ESE-1 transcription factor. HIDDEN cells are detected throughout stages of human carcinogenesis including primary human cervical intraepithelial neoplasias and HPV positive head and neck cancers, and a possible role in promoting viral carcinogenesis is supported by TCGA analyses. Single cell transcriptome information on HPV-infected versus uninfected epithelium will enable broader studies of the role of individual keratinocyte subpopulations in tumor virus infection and cancer evolution.

Single-cell RNA sequencing of human esophageal epithelium in homeostasis and allergic inflammation

Inflammation of the esophageal epithelium is a hallmark of eosinophilic esophagitis (EoE), an emerging chronic allergic disease. Herein, we probed human esophageal epithelial cells at single-cell resolution during homeostasis and EoE. During allergic inflammation, the epithelial differentiation program was blocked, leading to loss of KRT6high differentiated populations and expansion of TOP2high proliferating and DSPhigh, SERPINB3high transitioning populations; however, there was stability of the stem cell-enriched PDPNhigh basal epithelial compartment. This differentiation program blockade was associated with dysregulation of transcription factors, including nuclear receptor signalers, in the most differentiated epithelial cells and altered NOTCH-related cell-to-cell communication. Each epithelial population expressed genes with allergic disease risk variants, supporting their functional interplay. The esophageal epithelium differed notably between EoE in histologic remission and controls, indicating that remission is a transitory state poised to relapse. Collectively, our data uncover the dynamic nature of the inflamed human esophageal epithelium and provide a framework to better understand esophageal health and disease.

The genetic architecture of pediatric cardiomyopathy

To understand the genetic contribution to primary pediatric cardiomyopathy, we performed exome sequencing in a large cohort of 528 children with cardiomyopathy. Using clinical interpretation guidelines and targeting genes implicated in cardiomyopathy, we identified a genetic cause in 32% of affected individuals. Cardiomyopathy sub-phenotypes differed by ancestry, age at diagnosis, and family history. Infants < 1 year were less likely to have a molecular diagnosis (p < 0.001). Using a discovery set of 1,703 candidate genes and informatic tools, we identified rare and damaging variants in 56% of affected individuals. We see an excess burden of damaging variants in affected individuals as compared to two independent control sets, 1000 Genomes Project (p < 0.001) and SPARK parental controls (p < 1 × 10^-16). Cardiomyopathy variant burden remained enriched when stratified by ancestry, variant type, and sub-phenotype, emphasizing the importance of understanding the contribution of these factors to genetic architecture. Enrichment in this discovery candidate gene set suggests multigenic mechanisms underlie sub-phenotype-specific causes and presentations of cardiomyopathy. These results identify important information about the genetic architecture of pediatric cardiomyopathy and support recommendations for clinical genetic testing in children while illustrating differences in genetic architecture by age, ancestry, and sub-phenotype and providing rationale for larger studies to investigate multigenic contributions.

Mucosal Genomics Implicate Lymphocyte Activation and Lipid Metabolism in Refractory Environmental Enteric Dysfunction

BACKGROUND & AIMS

Environmental enteric dysfunction (EED) limits the Sustainable Development Goals of improved childhood growth and survival. We applied mucosal genomics to advance our understanding of EED.

METHODS

The Study of Environmental Enteropathy and Malnutrition (SEEM) followed 416 children from birth to 24 months in a rural district in Pakistan. Biomarkers were measured at 9 months and tested for association with growth at 24 months. The duodenal methylome and transcriptome were determined in 52 undernourished SEEM participants and 42 North American controls and patients with celiac disease.

RESULTS

After accounting for growth at study entry, circulating insulin-like growth factor-1 (IGF-1) and ferritin predicted linear growth, whereas leptin correlated with future weight gain. The EED transcriptome exhibited suppression of antioxidant, detoxification, and lipid metabolism genes, and induction of anti-microbial response, interferon, and lymphocyte activation genes. Relative to celiac disease, suppression of antioxidant and detoxification genes and induction of antimicrobial response genes were EED-specific. At the epigenetic level, EED showed hyper-methylation of epithelial metabolism and barrier function genes, and hypo-methylation of immune response and cell proliferation genes. Duodenal coexpression modules showed association between lymphocyte proliferation and epithelial metabolic genes and histologic severity, fecal energy loss, and wasting (weight-for-length/height Z < -2.0). Leptin was associated with expression of epithelial carbohydrate metabolism and stem cell renewal genes. Immune response genes were attenuated by giardia colonization.

CONCLUSIONS

Children with reduced circulating IGF-1 are more likely to experience stunting. Leptin and a gene signature for lymphocyte activation and dysregulated lipid metabolism are implicated in wasting, suggesting new approaches for EED refractory to nutritional intervention. ClinicalTrials.gov, Number: NCT03588013. (https://clinicaltrials.gov/ct2/show/NCT03588013).

Genetic Causes of Cardiomyopathy in Children: First Results From the Pediatric Cardiomyopathy Genes Study

Background

Pediatric cardiomyopathy is a genetically heterogeneous disease with substantial morbidity and mortality. Current guidelines recommend genetic testing in children with hypertrophic, dilated, or restrictive cardiomyopathy, but practice variations exist. Robust data on clinical testing practices and diagnostic yield in children are lacking. This study aimed to identify the genetic causes of cardiomyopathy in children and to investigate clinical genetic testing practices.

Methods and Results

Children with familial or idiopathic cardiomyopathy were enrolled from 14 institutions in North America. Probands underwent exome sequencing. Rare sequence variants in 37 known cardiomyopathy genes were assessed for pathogenicity using consensus clinical interpretation guidelines. Of the 152 enrolled probands, 41% had a family history of cardiomyopathy. Of 81 (53%) who had undergone clinical genetic testing for cardiomyopathy before enrollment, 39 (48%) had a positive result. Genetic testing rates varied from 0% to 97% between sites. A positive family history and hypertrophic cardiomyopathy subtype were associated with increased likelihood of genetic testing (P=0.005 and P=0.03, respectively). A molecular cause was identified in an additional 21% of the 63 children who did not undergo clinical testing, with positive results identified in both familial and idiopathic cases and across all phenotypic subtypes.

Conclusions

A definitive molecular genetic diagnosis can be made in a substantial proportion of children for whom the cause and heritable nature of their cardiomyopathy was previously unknown. Practice variations in genetic testing are great and should be reduced. Improvements can be made in comprehensive cardiac screening and predictive genetic testing in first‐degree relatives. Overall, our results support use of routine genetic testing in cases of both familial and idiopathic cardiomyopathy.

Registration

URL: https://www.clinicaltrials.gov; Unique identifier: NCT01873963.

Mucosal Inflammatory and Wound Healing Gene Programs Reveal Targets for Stricturing Behavior in Pediatric Crohn's Disease

BACKGROUND AND AIMS

Ileal strictures are the major indication for resective surgery in Crohn's disease (CD). We aimed to define ileal gene programs present at diagnosis linked with future stricturing behavior during five year follow-up, and to identify potential small molecules to reverse these gene signatures.

METHODS

Antimicrobial serologies and pre-treatment ileal gene expression were assessed in a representative subset of 249 CD patients within the RISK multicenter pediatric CD inception cohort study, including 113 that are unique to this report. These data were used to define genes associated with stricturing behavior and for model testing to predict stricturing behavior. A bioinformatics approach to define small molecules which may reverse the stricturing gene signature was applied.

RESULTS

19 of the 249 patients developed isolated B2 stricturing behavior during follow-up, while 218 remained B1 inflammatory. Using deeper RNA sequencing than in our prior report, we have now defined an inflammatory gene signature including an oncostatin M co-expression signature, tightly associated with extra-cellular matrix (ECM) gene expression in those who developed stricturing complications. We further computationally prioritize small molecules targeting macrophage and fibroblast activation and angiogenesis which may reverse the stricturing gene signature. A model containing ASCA and CBir1 serologies and a refined eight ECM gene set was significantly associated with stricturing development by year five after diagnosis (AUC (95th CI) = 0.82 (0.7-0.94)).

CONCLUSION

An ileal gene program for macrophage and fibroblast activation is linked to stricturing complications in treatment naïve pediatric CD, and may inform novel small molecule therapeutic approaches.

Ulcerative colitis mucosal transcriptomes reveal mitochondriopathy and personalized mechanisms underlying disease severity and treatment response

Molecular mechanisms driving disease course and response to therapy in ulcerative colitis (UC) are not well understood. Here, we use RNAseq to define pre-treatment rectal gene expression, and fecal microbiota profiles, in 206 pediatric UC patients receiving standardised therapy. We validate our key findings in adult and paediatric UC cohorts of 408 participants. We observe a marked suppression of mitochondrial genes and function across cohorts in active UC, and that increasing disease severity is notable for enrichment of adenoma/adenocarcinoma and innate immune genes. A subset of severity genes improves prediction of corticosteroid-induced remission in the discovery cohort; this gene signature is also associated with response to anti-TNFα and anti-α4β7 integrin in adults. The severity and therapeutic response gene signatures were in turn associated with shifts in microbes previously implicated in mucosal homeostasis. Our data provide insights into UC pathogenesis, and may prioritise future therapies for nonresponders to current approaches.

Whole-exome sequencing uncovers oxidoreductases DHTKD1 and OGDHL as linkers between mitochondrial dysfunction and eosinophilic esophagitis

Eosinophilic esophagitis (EoE) is an allergic inflammatory esophageal disorder with a complex underlying genetic etiology often associated with other comorbidities. Using whole-exome sequencing (WES) of 63 patients with EoE and 60 unaffected family members and family-based trio analysis, we sought to uncover rare coding variants. WES analysis identified 5 rare, damaging variants in dehydrogenase E1 and transketolase domain-containing 1 (DHTKD1). Rare variant burden analysis revealed an overabundance of putative, potentially damaging DHTKD1 mutations in EoE (P = 0.01). Interestingly, we also identified 7 variants in the DHTKD1 homolog oxoglutarate dehydrogenase-like (OGDHL). Using shRNA-transduced esophageal epithelial cells and/or patient fibroblasts, we further showed that disruption of normal DHTKD1 or OGDHL expression blunts mitochondrial function. Finally, we demonstrated that the loss of DHTKD1 expression increased ROS production and induced the expression of viperin, a gene previously shown to be involved in production of Th2 cytokines in T cells. Viperin had increased expression in esophageal biopsies of EoE patients compared with control individuals and was upregulated by IL-13 in esophageal epithelial cells. These data identify a series of rare genetic variants implicating DHTKD1 and OGDHL in the genetic etiology of EoE and underscore a potential pathogenic role for mitochondrial dysfunction in EoE.

Long ncRNA Landscape in the Ileum of Treatment-Naive Early-Onset Crohn Disease

BACKGROUND

Long noncoding RNAs (lncRNA) are key regulators of gene transcription and many show tissue-specific expression. We previously defined a novel inflammatory and metabolic ileal gene signature in treatment-naive pediatric Crohn disease (CD). We now extend our analyses to include potential regulatory lncRNA.

METHODS

Using RNAseq, we systematically profiled lncRNAs and protein-coding gene expression in 177 ileal biopsies. Co-expression analysis was used to identify functions and tissue-specific expression. RNA in situ hybridization was used to validate expression. Real-time polymerase chain reaction was used to test lncRNA regulation by IL-1β in Caco-2 enterocytes.

RESULTS

We characterize widespread dysregulation of 459 lncRNAs in the ileum of CD patients. Using only the lncRNA in discovery and independent validation cohorts showed patient classification as accurate as the protein-coding genes, linking lncRNA to CD pathogenesis. Co-expression and functional annotation enrichment analyses across several tissues and cell types 1showed that the upregulated LINC01272 is associated with a myeloid pro-inflammatory signature, whereas the downregulated HNF4A-AS1 exhibits association with an epithelial metabolic signature. We confirmed tissue-specific expression in biopsies using in situ hybridization, and validated regulation of prioritized lncRNA upon IL-1β exposure in differentiated Caco-2 cells. Finally, we identified significant correlations between LINC01272 and HNF4A-AS1 expression and more severe mucosal injury.

CONCLUSIONS

We systematically define differentially expressed lncRNA in the ileum of newly diagnosed pediatric CD. We show lncRNA utility to correctly classify disease or healthy states and demonstrate their regulation in response to an inflammatory signal. These lncRNAs, after mechanistic exploration, may serve as potential new tissue-specific targets for RNA-based interventions.

Integrated Genomic Analysis of Diverse Induced Pluripotent Stem Cells from the Progenitor Cell Biology Consortium

The rigorous characterization of distinct induced pluripotent stem cells (iPSC) derived from multiple reprogramming technologies, somatic sources, and donors is required to understand potential sources of variability and downstream potential. To achieve this goal, the Progenitor Cell Biology Consortium performed comprehensive experimental and genomic analyses of 58 iPSC from ten laboratories generated using a variety of reprogramming genes, vectors, and cells. Associated global molecular characterization studies identified functionally informative correlations in gene expression, DNA methylation, and/or copy-number variation among key developmental and oncogenic regulators as a result of donor, sex, line stability, reprogramming technology, and cell of origin. Furthermore, X-chromosome inactivation in PSC produced highly correlated differences in teratoma-lineage staining and regulator expression upon differentiation. All experimental results, and raw, processed, and metadata from these analyses, including powerful tools, are interactively accessible from a new online portal at https://www.synapse.org to serve as a reusable resource for the stem cell community.

•

Comprehensive portal for diverse iPSC, protocols, metadata, and genomic assays

•

Recurrent CNV occur during reprogramming, impact oncogenes and tumor suppressors

•

DNA methylation is influenced by cell of origin in iPSC

•

PSC X-chromosome inactivation impacts lineage differentiation outcomes

Intrinsic Age-Dependent Changes and Cell-Cell Contacts Regulate Nephron Progenitor Lifespan

During fetal development, nephrons of the metanephric kidney form from a mesenchymal progenitor population that differentiates en masse before or shortly after birth. We explored intrinsic and extrinsic mechanisms controlling progenitor lifespan in a transplantation assay that allowed us to compare engraftment of old and young progenitors into the same young niche. The progenitors displayed an age-dependent decrease in proliferation and concomitant increase in niche exit rates. Single-cell transcriptome profiling revealed progressive age-dependent changes, with heterogeneity increasing in older populations. Age-dependent elevation in mTor and reduction in Fgf20 could contribute to increased exit rates. Importantly, 30% of old progenitors remained in the niche for up to 1 week post engraftment, a net gain of 50% to their lifespan, but only if surrounded by young neighbors. We provide evidence in support of a model in which intrinsic age-dependent changes affect inter-progenitor interactions that drive cessation of nephrogenesis.

IRAK1 is a novel DEK transcriptional target and is essential for head and neck cancer cell survival

The chromatin-binding DEK protein was recently reported to promote the growth of HPV+ and HPV- head and neck squamous cell carcinomas (HNSCCs). Relevant cellular and molecular mechanism(s) controlled by DEK in HNSCC remain poorly understood. While DEK is known to regulate specific transcriptional targets, global DEK-dependent gene networks in HNSCC are unknown. To identify DEK transcriptional signatures we performed RNA-Sequencing (RNA-Seq) in HNSCC cell lines that were either proficient or deficient for DEK. Bioinformatic analyses and subsequent validation revealed that IRAK1, a regulator of inflammatory signaling, and IRAK1-dependent regulatory networks were significantly repressed upon DEK knockdown in HNSCC. According to TCGA data, 14% of HNSCC specimens overexpressed IRAK1, thus supporting possible oncogenic functions. Furthermore, genetic or pharmacologic inhibition of IRAK1 in HNSCC cell lines was sufficient to attenuate downstream signaling such as ERK1/2 and to induce HNSCC cell death by apoptosis. Finally, targeting DEK and IRAK1 simultaneously enhanced cell death as compared to targeting either alone. Our findings reveal that IRAK1 promotes cell survival and is an attractive therapeutic target in HNSCC cells. Thus, we propose a model wherein IRAK1 stimulates tumor signaling and phenotypes both independently and in conjunction with DEK.

Pediatric Crohn disease patients exhibit specific ileal transcriptome and microbiome signature

Interactions between the host and gut microbial community likely contribute to Crohn disease (CD) pathogenesis; however, direct evidence for these interactions at the onset of disease is lacking. Here, we characterized the global pattern of ileal gene expression and the ileal microbial community in 359 treatment-naive pediatric patients with CD, patients with ulcerative colitis (UC), and control individuals. We identified core gene expression profiles and microbial communities in the affected CD ilea that are preserved in the unaffected ilea of patients with colon-only CD but not present in those with UC or control individuals; therefore, this signature is specific to CD and independent of clinical inflammation. An abnormal increase of antimicrobial dual oxidase (DUOX2) expression was detected in association with an expansion of Proteobacteria in both UC and CD, while expression of lipoprotein APOA1 gene was downregulated and associated with CD-specific alterations in Firmicutes. The increased DUOX2 and decreased APOA1 gene expression signature favored oxidative stress and Th1 polarization and was maximally altered in patients with more severe mucosal injury. A regression model that included APOA1 gene expression and microbial abundance more accurately predicted month 6 steroid-free remission than a model using clinical factors alone. These CD-specific host and microbe profiles identify the ileum as the primary inductive site for all forms of CD and may direct prognostic and therapeutic approaches.

Vertical sleeve gastrectomy reduces hepatic steatosis while increasing serum bile acids in a weight-loss-independent manner

OBJECTIVE

Our objective was to investigate the role of bile acids in hepatic steatosis reduction after vertical sleeve gastrectomy (VSG).

DESIGN AND METHODS

High fat diet (HFD)-induced obese C57Bl/6 mice were randomized to VSG, Sham operation (Sham), Sham operation with pair feeding to VSG (Sham-PF), or nonsurgical controls (Naïve). All mice were on HFD until sacrifice. Mice were observed postsurgery and data for body weight, body composition, metabolic parameters, serum bile acid level and composition were collected. Further hepatic gene expression by mRNA-seq and RT-PCR analysis was assessed.

RESULTS

VSG and Sham-PF mice lost equal weight postsurgery while VSG mice had the lowest hepatic triglyceride content at sacrifice. The VSG mice had elevated serum bile acid levels that positively correlated with maximal weight loss. Serum bile composition in the VSG group had increased cholic and tauroursodeoxycholic acid. These bile acid composition changes in VSG mice explained observed downregulation of hepatic lipogenic and bile acid synthetic genes.

CONCLUSION

VSG in obese mice results in greater hepatic steatosis reduction than seen with caloric restriction alone. VSG surgery increases serum bile acids that correlate with weight lost postsurgery and changes serum bile composition that could explain suppression of hepatic genes responsible for lipogenesis.

Increased protein-coding mutations in the mitochondrial genome of African American women with preeclampsia

Preeclampsia occurs more frequently in women of African ancestry. The cause of this hypertensive complication is unclear, but placental oxidative stress may play a role. Because mitochondria are the major sites of oxidative phosphorylation, we hypothesized that placentas of preeclamptic pregnancies harbor mitochondrial DNA (mtDNA) mutations. Next-generation sequencing of placental mtDNA in African American preeclamptics (N = 30) and controls (N = 38) from Chicago revealed significant excesses in preeclamptics of nonsynonymous substitutions in protein-coding genes and mitochondrially encoded nicotinamide adenine dinucleotide dehydrogenase 5 gene and an increase in the substitution rate (P = .0001). Moreover, 88% of preeclamptics and 53% of controls carried at least one nonsynonymous substitution (P = .005; odds ratio [OR] = 6.36, 95% confidence interval [CI]: 1.5-39.1). These results were not replicated in a sample of African American preeclamptics (N = 162) and controls (N = 171) from Detroit. Differences in study design and heterogeneity may account for this lack of replication. Nonsynonymous substitutions in mtDNA may be risk factors for preeclampsia in some African American women, but additional studies are required to establish this relationship.

Identification of genes and pathways involved in retinal neovascularization by microarray analysis of two animal models of retinal angiogenesis

PURPOSE

Comparative retinal gene expression analysis in two rodent models of oxygen-induced retinopathy (OIR) was performed to identify the genes and pathways involved in retinal neovascularization.

METHODS

Three independent experimental runs were conducted for each species, according to standard protocols for induction of OIR. Total retinal RNA was isolated at two time points, corresponding to the early response to relative hypoxia (P13 in mouse, P15 in rat) and to the later phase of maximum retinal neovascularization (P18 in mouse, P20 in rat) and was used to prepare labeled probes for hybridization. Gene expression was compared between normal and experimental conditions for each species at each time point. Probesets with a false-discovery rate of <or=0.05 were considered significantly different and were classified as cellular functions or biological pathways. Changes in expression of selected genes were confirmed by quantitative rtPCR.

RESULTS

At the early time point, there were changes in 43 genes in each species, with two in common. Increased expression of members of the VEGF and ephrin receptor signaling pathways were identified in both models. At the later time point, there were changes in 26 genes in the rat and in 1622 in the mouse, with 13 in common. Four pathways were identified in both models.

CONCLUSIONS

Genes and pathways known to be involved in angiogenesis, as well as other biologically plausible genes and pathways, were identified. This work serves as a comprehensive resource for the study of retinal neovascularization and identification of potential rational targets for antiangiogenic therapy.

Complementary RNA amplification methods enhance microarray identification of transcripts expressed in the C. elegans nervous system

Background

DNA microarrays provide a powerful method for global analysis of gene expression. The application of this technology to specific cell types and tissues, however, is typically limited by small amounts of available mRNA, thereby necessitating amplification. Here we compare microarray results obtained with two different methods of RNA amplification to profile gene expression in the C. elegans larval nervous system.

Results

We used the mRNA-tagging strategy to isolate transcripts specifically from C. elegans larval neurons. The WT-Ovation Pico System (WT-Pico) was used to amplify 2 ng of pan-neural RNA to produce labeled cDNA for microarray analysis. These WT-Pico-derived data were compared to microarray results obtained with a labeled aRNA target generated by two rounds of In Vitro Transcription (IVT) of 25 ng of pan-neural RNA. WT-Pico results in a higher fraction of present calls than IVT, a finding consistent with the proposal that DNA-DNA hybridization results in lower mismatch signals than the RNA-DNA heteroduplexes produced by IVT amplification. Microarray data sets from these samples were compared to a reference profile of all larval cells to identify transcripts with elevated expression in neurons. These results were validated by the high proportion of known neuron-expressed genes detected in these profiles and by promoter-GFP constructs for previously uncharacterized genes in these data sets. Together, the IVT and WT-Pico methods identified 2,173 unique neuron-enriched transcripts. Only about half of these transcripts (1,044), however, are detected as enriched by both IVT and WT-Pico amplification.

Conclusion

We show that two different methods of RNA amplification, IVT and WT-Pico, produce valid microarray profiles of gene expression in the C. elegans larval nervous system with a low rate of false positives. However, our results also show that each method of RNA amplification detects a unique subset of bona fide neural-enriched transcripts and thus a wider array of authentic neural genes are identified by the combination of these data sets than by the microarray profiles obtained with either method of RNA amplification alone. With its relative ease of implementation and greater sensitivity, WT-Pico is the preferred method of amplification for cases in which sample RNA is limiting.