Unravelling genetic architecture of circulatory amino acid levels, and their effect on risk of complex disorders.

Variations in serum amino acid levels are linked to a multitude of complex disorders. We report the largest genome-wide association study (GWAS) on nine serum amino acids in the UK Biobank participants (117 944, European descent). We identified 34 genomic loci for circulatory levels of alanine, 48 loci for glutamine, 44 loci for glycine, 16 loci for histidine, 11 loci for isoleucine, 19 loci for leucine, 9 loci for phenylalanine, 32 loci for tyrosine and 20 loci for valine. Our gene-based analysis mapped 46-293 genes associated with serum amino acids, including MIP, GLS2, SLC gene family, GCKR, LMO1, CPS1 and COBLL1.The gene-property analysis across 30 tissues highlighted enriched expression of the identified genes in liver tissues for all studied amino acids, except for isoleucine and valine, in muscle tissues for serum alanine and glycine, in adrenal gland tissues for serum isoleucine and leucine, and in pancreatic tissues for serum phenylalanine. Mendelian randomization (MR) phenome-wide association study analysis and subsequent two-sample MR analysis provided evidence that every standard deviation increase in valine is associated with 35% higher risk of type 2 diabetes and elevated levels of serum alanine and branched-chain amino acids with higher levels of total cholesterol, triglyceride and low-density lipoprotein, and lower levels of high-density lipoprotein. In contrast to reports by observational studies, MR analysis did not support a causal association between studied amino acids and coronary artery disease, Alzheimer's disease, breast cancer or prostate cancer. In conclusion, we explored the genetic architecture of serum amino acids and provided evidence supporting a causal role of amino acids in cardiometabolic health.

Unravelling morphoea aetiopathogenesis by next-generation sequencing of paired skin biopsies.

BACKGROUND: Morphoea can have a significant disease burden. Aetiopathogenesis remains poorly understood, with very limited existing genetic studies. Linear morphoea (LM) may follow Blascho's lines of epidermal development, providing potential pathogenic clues. OBJECTIVE: The first objective of this study was to identify the presence of primary somatic epidermal mosaicism in LM. The second objective was tTo explore differential gene expression in morphoea epidermis and dermis to identify potential pathogenic molecular pathways and tissue layer cross-talk. METHODOLOGY: Skin biopsies from paired affected and contralateral unaffected skin were taken from 16 patients with LM. Epidermis and dermis were isolated using a 2-step chemical-physical separation protocol. Whole Genome Sequencing (WGS; n = 4 epidermal) and RNA-seq (n = 5-epidermal, n = 5-dermal) with gene expression analysis via GSEA-MSigDBv6.3 and PANTHER-v14.1 pathway analyses, were performed. RTqPCR and immunohistochemistry were used to replicate key results. RESULTS: Sixteen participants (93.8% female, mean age 27.7 yrs disease-onset) were included. Epidermal WGS identified no single affected gene or SNV. However, many potential disease-relevant pathogenic variants were present, including ADAMTSL1 and ADAMTS16. A highly proliferative, inflammatory and profibrotic epidermis was seen, with significantly-overexpressed TNFα-via-NFkB, TGFß, IL6/JAKSTAT and IFN-signaling, apoptosis, p53 and KRAS-responses. Upregulated IFI27 and downregulated LAMA4 potentially represent initiating epidermal 'damage' signals and enhanced epidermal-dermal communication. Morphoea dermis exhibited significant profibrotic, B-cell and IFN-signatures, and upregulated morphogenic patterning pathways such as Wnt. CONCLUSION: This study supports the absence of somatic epidermal mosaicism in LM, and identifies potential disease-driving epidermal mechanisms, epidermal-dermal interactions and disease-specific dermal differential-gene-expression in morphoea. We propose a potential molecular narrative for morphoea aetiopathogenesis which could help guide future targeted studies and therapies.

Role of CD14+ monocyte-derived oxidised mitochondrial DNA in the inflammatory interferon type 1 signature in juvenile dermatomyositis.

OBJECTIVES: To define the host mechanisms contributing to the pathological interferon (IFN) type 1 signature in Juvenile dermatomyositis (JDM). METHODS: RNA-sequencing was performed on CD4+, CD8+, CD14+ and CD19+ cells sorted from pretreatment and on-treatment JDM (pretreatment n=10, on-treatment n=11) and age/sex-matched child healthy-control (CHC n=4) peripheral blood mononuclear cell (PBMC). Mitochondrial morphology and superoxide were assessed by fluorescence microscopy, cellular metabolism by 13C glucose uptake assays, and oxidised mitochondrial DNA (oxmtDNA) content by dot-blot. Healthy-control PBMC and JDM pretreatment PBMC were cultured with IFN-α, oxmtDNA, cGAS-inhibitor, TLR-9 antagonist and/or n-acetyl cysteine (NAC). IFN-stimulated gene (ISGs) expression was measured by qPCR. Total numbers of patient and controls for functional experiments, JDM n=82, total CHC n=35. RESULTS: Dysregulated mitochondrial-associated gene expression correlated with increased ISG expression in JDM CD14+ monocytes. Altered mitochondrial-associated gene expression was paralleled by altered mitochondrial biology, including 'megamitochondria', cellular metabolism and a decrease in gene expression of superoxide dismutase (SOD)1. This was associated with enhanced production of oxidised mitochondrial (oxmt)DNA. OxmtDNA induced ISG expression in healthy PBMC, which was blocked by targeting oxidative stress and intracellular nucleic acid sensing pathways. Complementary experiments showed that, under in vitro experimental conditions, targeting these pathways via the antioxidant drug NAC, TLR9 antagonist and to a lesser extent cGAS-inhibitor, suppressed ISG expression in pretreatment JDM PBMC. CONCLUSIONS: These results describe a novel pathway where altered mitochondrial biology in JDM CD14+ monocytes lead to oxmtDNA production and stimulates ISG expression. Targeting this pathway has therapeutical potential in JDM and other IFN type 1-driven autoimmune diseases.

Failures of nerve regeneration caused by aging or chronic denervation are rescued by restoring Schwann cell c-Jun.

After nerve injury, myelin and Remak Schwann cells reprogram to repair cells specialized for regeneration. Normally providing strong regenerative support, these cells fail in aging animals, and during chronic denervation that results from slow axon growth. This impairs axonal regeneration and causes significant clinical problems. In mice, we find that repair cells express reduced c-Jun protein as regenerative support provided by these cells declines during aging and chronic denervation. In both cases, genetically restoring Schwann cell c-Jun levels restores regeneration to control levels. We identify potential gene candidates mediating this effect and implicate Shh in the control of Schwann cell c-Jun levels. This establishes that a common mechanism, reduced c-Jun in Schwann cells, regulates success and failure of nerve repair both during aging and chronic denervation. This provides a molecular framework for addressing important clinical problems, suggesting molecular pathways that can be targeted to promote repair in the PNS.

Conserved properties of genetic architecture of renal and fat transcriptomes in rat models of insulin resistance.

To define renal molecular mechanisms that are affected by permanent hyperglycaemia and might promote phenotypes relevant to diabetic nephropathy, we carried out linkage analysis of genome-wide gene transcription in the kidneys of F2 offspring from the Goto-Kakizaki (GK) rat model of type 2 diabetes and normoglycaemic Brown Norway (BN) rats. We mapped 2526 statistically significant expression quantitative trait loci (eQTLs) in the cross. More than 40% of eQTLs mapped in the close vicinity of the linked transcripts, underlying possible cis-regulatory mechanisms of gene expression. We identified eQTL hotspots on chromosomes 5 and 9 regulating the expression of 80-165 genes, sex or cross direction effects, and enriched metabolic and immunological processes by segregating GK alleles. Comparative analysis with adipose tissue eQTLs in the same cross showed that 496 eQTLs, in addition to the top enriched biological pathways, are conserved in the two tissues. Extensive similarities in eQTLs mapped in the GK rat and in the spontaneously hypertensive rat (SHR) suggest a common aetiology of disease phenotypes common to the two strains, including insulin resistance, which is a prominent pathophysiological feature in both GK rats and SHRs. Our data shed light on shared and tissue-specific molecular mechanisms that might underlie aetiological aspects of insulin resistance in the context of spontaneously occurring hyperglycaemia and hypertension.

Systems Genetics of Hepatic Metabolome Reveals Octopamine as a Target for Non-Alcoholic Fatty Liver Disease Treatment.

Non-alcoholic fatty liver disease (NAFLD) is often associated with obesity and type 2 diabetes. To disentangle etiological relationships between these conditions and identify genetically-determined metabolites involved in NAFLD processes, we mapped 1H nuclear magnetic resonance (NMR) metabolomic and disease-related phenotypes in a mouse F2 cross derived from strains showing resistance (BALB/c) and increased susceptibility (129S6) to these diseases. Quantitative trait locus (QTL) analysis based on single nucleotide polymorphism (SNP) genotypes identified diet responsive QTLs in F2 mice fed control or high fat diet (HFD). In HFD fed F2 mice we mapped on chromosome 18 a QTL regulating liver micro- and macrovesicular steatosis and inflammation, independently from glucose intolerance and adiposity, which was linked to chromosome 4. Linkage analysis of liver metabolomic profiling data identified a QTL for octopamine, which co-localised with the QTL for liver histopathology in the cross. Functional relationship between these two QTLs was validated in vivo in mice chronically treated with octopamine, which exhibited reduction in liver histopathology and metabolic benefits, underlining its role as a mechanistic biomarker of fatty liver with potential therapeutic applications.

CD19+CD24hiCD38hi B Cells Are Expanded in Juvenile Dermatomyositis and Exhibit a Pro-Inflammatory Phenotype After Activation Through Toll-Like Receptor 7 and Interferon-α.

Juvenile dermatomyositis (JDM) is a rare form of childhood autoimmune myositis that presents with proximal muscle weakness and skin rash. B cells are strongly implicated in the pathogenesis of the disease, but the underlying mechanisms are unknown. Therefore, the main objective of our study was to investigate mechanisms driving B cell lymphocytosis and define pathological features of B cells in JDM patients. Patients were recruited through the UK JDM Cohort and Biomarker study. Peripheral blood B cell subpopulations were immunophenotyped by flow cytometry. The results identified that immature transitional B cells were significantly expanded in active JDM, actively dividing, and correlated positively with disease activity. Protein and RNAseq analysis revealed high interferon alpha (IFNα) and TLR7-pathway signatures pre-treatment. Stimulation of B cells through TLR7/8 promoted both IL-10 and IL-6 production in controls but failed to induce IL-10 in JDM patient cells. Interrogation of the CD40-CD40L pathway (known to induce B cell IL-10 and IL-6) revealed similar expression of IL-10 and IL-6 in B cells cultured with CD40L from both JDM patients and controls. In conclusion, JDM patients with active disease have a significantly expanded immature transitional B cell population which correlated with the type I IFN signature. Activation through TLR7 and IFNα may drive the expansion of immature transitional B cells in JDM and skew the cells toward a pro-inflammatory phenotype.

Topological analysis of metabolic networks integrating co-segregating transcriptomes and metabolomes in type 2 diabetic rat congenic series.

BACKGROUND: The genetic regulation of metabolic phenotypes (i.e., metabotypes) in type 2 diabetes mellitus occurs through complex organ-specific cellular mechanisms and networks contributing to impaired insulin secretion and insulin resistance. Genome-wide gene expression profiling systems can dissect the genetic contributions to metabolome and transcriptome regulations. The integrative analysis of multiple gene expression traits and metabolic phenotypes (i.e., metabotypes) together with their underlying genetic regulation remains a challenge. Here, we introduce a systems genetics approach based on the topological analysis of a combined molecular network made of genes and metabolites identified through expression and metabotype quantitative trait locus mapping (i.e., eQTL and mQTL) to prioritise biological characterisation of candidate genes and traits. METHODS: We used systematic metabotyping by 1H NMR spectroscopy and genome-wide gene expression in white adipose tissue to map molecular phenotypes to genomic blocks associated with obesity and insulin secretion in a series of rat congenic strains derived from spontaneously diabetic Goto-Kakizaki (GK) and normoglycemic Brown-Norway (BN) rats. We implemented a network biology strategy approach to visualize the shortest paths between metabolites and genes significantly associated with each genomic block. RESULTS: Despite strong genomic similarities (95-99 %) among congenics, each strain exhibited specific patterns of gene expression and metabotypes, reflecting the metabolic consequences of series of linked genetic polymorphisms in the congenic intervals. We subsequently used the congenic panel to map quantitative trait loci underlying specific mQTLs and genome-wide eQTLs. Variation in key metabolites like glucose, succinate, lactate, or 3-hydroxybutyrate and second messenger precursors like inositol was associated with several independent genomic intervals, indicating functional redundancy in these regions. To navigate through the complexity of these association networks we mapped candidate genes and metabolites onto metabolic pathways and implemented a shortest path strategy to highlight potential mechanistic links between metabolites and transcripts at colocalized mQTLs and eQTLs. Minimizing the shortest path length drove prioritization of biological validations by gene silencing. CONCLUSIONS: These results underline the importance of network-based integration of multilevel systems genetics datasets to improve understanding of the genetic architecture of metabotype and transcriptomic regulation and to characterize novel functional roles for genes determining tissue-specific metabolism.

Transcriptome Profiling in Rat Inbred Strains and Experimental Cross Reveals Discrepant Genetic Architecture of Genome-Wide Gene Expression.

To test the impact of genetic heterogeneity on cis- and trans-mediated mechanisms of gene expression regulation, we profiled the transcriptome of adipose tissue in 20 inbred congenic strains derived from diabetic Goto-Kakizaki (GK) rats and Brown-Norway (BN) controls, which contain well-defined blocks (1-183 Mb) of genetic polymorphisms, and in 123 genetically heterogeneous rats of an (GK × BN)F2 offspring. Within each congenic we identified 73-1351 differentially expressed genes (DEGs), only 7.7% of which mapped within the congenic blocks, and which may be regulated in cis The remainder localized outside the blocks, and therefore must be regulated in trans Most trans-regulated genes exhibited approximately twofold expression changes, consistent with monoallelic expression. Altered biological pathways were replicated between congenic strains sharing blocks of genetic polymorphisms, but polymorphisms at different loci also had redundant effects on transcription of common distant genes and pathways. We mapped 2735 expression quantitative trait loci (eQTL) in the F2 cross, including 26% predominantly cis-regulated genes, which validated DEGs in congenic strains. A hotspot of >300 eQTL in a 10 cM region of chromosome 1 was enriched in DEGs in a congenic strain. However, many DEGs among GK, BN and congenic strains did not replicate as eQTL in F2 hybrids, demonstrating distinct mechanisms of gene expression when alleles segregate in an outbred population or are fixed homozygous across the entire genome or in short genomic regions. Our analysis provides conceptual advances in our understanding of the complex architecture of genome expression and pathway regulation, and suggests a prominent impact of epistasis and monoallelic expression on gene transcription.

Genetically distinct leukemic stem cells in human CD34- acute myeloid leukemia are arrested at a hemopoietic precursor-like stage.

Our understanding of the perturbation of normal cellular differentiation hierarchies to create tumor-propagating stem cell populations is incomplete. In human acute myeloid leukemia (AML), current models suggest transformation creates leukemic stem cell (LSC) populations arrested at a progenitor-like stage expressing cell surface CD34. We show that in ∼25% of AML, with a distinct genetic mutation pattern where >98% of cells are CD34(-), there are multiple, nonhierarchically arranged CD34(+) and CD34(-) LSC populations. Within CD34(-) and CD34(+) LSC-containing populations, LSC frequencies are similar; there are shared clonal structures and near-identical transcriptional signatures. CD34(-) LSCs have disordered global transcription profiles, but these profiles are enriched for transcriptional signatures of normal CD34(-) mature granulocyte-macrophage precursors, downstream of progenitors. But unlike mature precursors, LSCs express multiple normal stem cell transcriptional regulators previously implicated in LSC function. This suggests a new refined model of the relationship between LSCs and normal hemopoiesis in which the nature of genetic/epigenetic changes determines the disordered transcriptional program, resulting in LSC differentiation arrest at stages that are most like either progenitor or precursor stages of hemopoiesis.

Genetic control of differential acetylation in diabetic rats.

Post-translational protein modifications such as acetylation have significant regulatory roles in metabolic processes, but their relationship to both variation in gene expression and DNA sequence is unclear. We address this question in the Goto-Kakizaki (GK) rat inbred strain, a model of polygenic type 2 diabetes. Expression of the NAD-dependent deacetylase Sirtuin-3 is down-regulated in GK rats compared to normoglycemic Brown Norway (BN) rats. We show first that a promoter SNP causes down-regulation of Sirtuin-3 expression in GK rats. We then use mass-spectrometry to identify proteome-wide differential lysine acetylation of putative Sirtuin-3 protein targets in livers of GK and BN rats. These include many proteins in pathways connected to diabetes and metabolic syndrome. We finally sequence GK and BN liver transcriptomes and find that mRNA expression of these targets does not differ significantly between GK and BN rats, in contrast to other components of the same pathways. We conclude that physiological differences between GK and BN rats are mediated by a combination of differential protein acetylation and gene transcription and that genetic variation can modulate acetylation independently of expression.

A bi-modal function of Wnt signalling directs an FGF activity gradient to spatially regulate neuronal differentiation in the midbrain.

FGFs and Wnts are important morphogens during midbrain development, but their importance and potential interactions during neurogenesis are poorly understood. We have employed a combination of genetic and pharmacological manipulations in zebrafish to show that during neurogenesis FGF activity occurs as a gradient along the anterior-posterior axis of the dorsal midbrain and directs spatially dynamic expression of the Hairy gene her5. As FGF activity diminishes during development, Her5 is lost and differentiation of neuronal progenitors occurs in an anterior-posterior manner. We generated mathematical models to explain how Wnt and FGFs direct the spatial differentiation of neurons in the midbrain through Wnt regulation of FGF signalling. These models suggested that a negative-feedback loop controlled by Wnt is crucial for regulating FGF activity. We tested Sprouty genes as mediators of this regulatory loop using conditional mouse knockouts and pharmacological manipulations in zebrafish. These reveal that Sprouty genes direct the positioning of early midbrain neurons and are Wnt responsive in the midbrain. We propose a model in which Wnt regulates FGF activity at the isthmus by driving both FGF and Sprouty gene expression. This controls a dynamic, posteriorly retracting expression of her5 that directs neuronal differentiation in a precise spatiotemporal manner in the midbrain.

Genome sequencing reveals loci under artificial selection that underlie disease phenotypes in the laboratory rat.

Large numbers of inbred laboratory rat strains have been developed for a range of complex disease phenotypes. To gain insights into the evolutionary pressures underlying selection for these phenotypes, we sequenced the genomes of 27 rat strains, including 11 models of hypertension, diabetes, and insulin resistance, along with their respective control strains. Altogether, we identified more than 13 million single-nucleotide variants, indels, and structural variants across these rat strains. Analysis of strain-specific selective sweeps and gene clusters implicated genes and pathways involved in cation transport, angiotensin production, and regulators of oxidative stress in the development of cardiovascular disease phenotypes in rats. Many of the rat loci that we identified overlap with previously mapped loci for related traits in humans, indicating the presence of shared pathways underlying these phenotypes in rats and humans. These data represent a step change in resources available for evolutionary analysis of complex traits in disease models.

Causes and consequences of chromatin variation between inbred mice.

Variation at regulatory elements, identified through hypersensitivity to digestion by DNase I, is believed to contribute to variation in complex traits, but the extent and consequences of this variation are poorly characterized. Analysis of terminally differentiated erythroblasts in eight inbred strains of mice identified reproducible variation at approximately 6% of DNase I hypersensitive sites (DHS). Only 30% of such variable DHS contain a sequence variant predictive of site variation. Nevertheless, sequence variants within variable DHS are more likely to be associated with complex traits than those in non-variant DHS, and variants associated with complex traits preferentially occur in variable DHS. Changes at a small proportion (less than 10%) of variable DHS are associated with changes in nearby transcriptional activity. Our results show that whilst DNA sequence variation is not the major determinant of variation in open chromatin, where such variants exist they are likely to be causal for complex traits.

The light responsive transcriptome of the zebrafish: function and regulation.

Most organisms possess circadian clocks that are able to anticipate the day/night cycle and are reset or "entrained" by the ambient light. In the zebrafish, many organs and even cultured cell lines are directly light responsive, allowing for direct entrainment of the clock by light. Here, we have characterized light induced gene transcription in the zebrafish at several organizational levels. Larvae, heart organ cultures and cell cultures were exposed to 1- or 3-hour light pulses, and changes in gene expression were compared with controls kept in the dark. We identified 117 light regulated genes, with the majority being induced and some repressed by light. Cluster analysis groups the genes into five major classes that show regulation at all levels of organization or in different subset combinations. The regulated genes cover a variety of functions, and the analysis of gene ontology categories reveals an enrichment of genes involved in circadian rhythms, stress response and DNA repair, consistent with the exposure to visible wavelengths of light priming cells for UV-induced damage repair. Promoter analysis of the induced genes shows an enrichment of various short sequence motifs, including E- and D-box enhancers that have previously been implicated in light regulation of the zebrafish period2 gene. Heterologous reporter constructs with sequences matching these motifs reveal light regulation of D-box elements in both cells and larvae. Morpholino-mediated knock-down studies of two homologues of the D-box binding factor Tef indicate that these are differentially involved in the cell autonomous light induction in a gene-specific manner. These findings suggest that the mechanisms involved in period2 regulation might represent a more general pathway leading to light induced gene expression.

Aplexone targets the HMG-CoA reductase pathway and differentially regulates arteriovenous angiogenesis.

Arterial and venous endothelial cells exhibit distinct molecular characteristics at early developmental stages. These lineage-specific molecular programs are instructive to the development of distinct vascular architectures and physiological conditions of arteries and veins, but their roles in angiogenesis remain unexplored. Here, we show that the caudal vein plexus in zebrafish forms by endothelial cell sprouting, migration and anastomosis, providing a venous-specific angiogenesis model. Using this model, we have identified a novel compound, aplexone, which effectively suppresses venous, but not arterial, angiogenesis. Multiple lines of evidence indicate that aplexone differentially regulates arteriovenous angiogenesis by targeting the HMG-CoA reductase (HMGCR) pathway. Treatment with aplexone affects the transcription of enzymes in the HMGCR pathway and reduces cellular cholesterol levels. Injecting mevalonate, a metabolic product of HMGCR, reverses the inhibitory effect of aplexone on venous angiogenesis. In addition, aplexone treatment inhibits protein prenylation and blocking the activity of geranylgeranyl transferase induces a venous angiogenesis phenotype resembling that observed in aplexone-treated embryos. Furthermore, endothelial cells of venous origin have higher levels of proteins requiring geranylgeranylation than arterial endothelial cells and inhibiting the activity of Rac or Rho kinase effectively reduces the migration of venous, but not arterial, endothelial cells. Taken together, our findings indicate that angiogenesis is differentially regulated by the HMGCR pathway via an arteriovenous-dependent requirement for protein prenylation in zebrafish and human endothelial cells.

Simplet controls cell proliferation and gene transcription during zebrafish caudal fin regeneration.

Two hallmarks of vertebrate epimorphic regeneration are a significant increase in the proliferation of normally quiescent cells and a re-activation of genes that are active during embryonic development. It is unclear what the molecular determinants are that regulate these events and how they are coordinated. Zebrafish have the ability to regenerate several compound structures by regulating cell proliferation and gene transcription. We report that fam53b/simplet (smp) regulates both cell proliferation and the transcription of specific genes. In situ hybridization and quantitative RT-PCR experiments showed that amputation of zebrafish hearts and fins resulted in strong up-regulation of the smp gene. In regenerating adult fin, smp expression remained strong in the distal mesenchyme which later expanded to the basal layers of the distal epidermis and distal tip epithelium. Morpholino knockdown of smp reduced regenerative outgrowth by decreasing cell proliferation as measured by BrdU incorporation and histone H3 phosphorylation. In addition, smp knockdown increased the expression of msxb, msxc, and shh, as well as the later formation of ectopic bone. Taken together, these data indicate a requirement for smp in fin regeneration through control of cell proliferation, the regulation of specific genes and proper bone patterning.

Differential gene expression as a toxicant-sensitive endpoint in zebrafish embryos and larvae.

The zebrafish (Danio rerio) embryo toxicity test (DarT) is under consideration as an alternative to the acute fish toxicity test. Microscopically visible developmental disorders or death are the endpoints used to report on toxicity in DarT. These endpoints are easily observed. They, however, rarely reveal mechanisms leading to a toxic effect and are relatively insensitive compared to chronic toxic effects. We hypothesized that, by using gene expression profiles as an additional endpoint, it may be possible to increase the sensitivity and predictive value of DarT. Therefore, as a proof of principle, we exposed zebrafish embryos to the reference compound 3,4-dichloroaniline (3,4-DCA) and analyzed gene expression patterns with a 14k oligonucleotide array. Important stress response genes not included in the microarray were additionally quantified by reverse transcriptase polymerase chain reaction. Six genes involved in biotransformation (cyp1a, ahr2), stress response (nfe212, maft, hmox1) and cell cycle control (fzr1) were significantly regulated. With the exception of fzr1, these genes proved to be differentially expressed in post hatch life stages as well. The identified genes point toward an aryl hydrocarbon receptor-mediated response. Differential gene expression in embryos exposed for 48 h was observed at 3,4-DCA concentrations as low as 0.78 microM, which is more than 10-fold below the concentrations that elicited visible toxic effects. Upon exposure for 5 days, differential expression was detected at concentrations as low as 0.22 microM of 3,4-DCA, which was close to the lowest observed effect concentration (0.11 microM) in the 30-day early life stage test. This study therefore indicates that gene expression analysis in DarT is able to reveal mechanistic information and may also be exploited for the development of replacement methods for chronic fish tests.

Transcriptome profiling of adult zebrafish at the late stage of chronic tuberculosis due to Mycobacterium marinum infection.

The Mycobacterium marinum-zebrafish infection model was used in this study for analysis of a host transcriptome response to mycobacterium infection at the organismal level. RNA isolated from adult zebrafish that showed typical signs of fish tuberculosis due to a chronic progressive infection with M. marinum was compared with RNA from healthy fish in microarray analyses. Spotted oligonucleotide sets (designed by Sigma-Compugen and MWG) and Affymetrix GeneChips were used, in total comprising 45,465 zebrafish transcript annotations. Based on a detailed comparative analysis and quantitative reverse transcriptase-PCR analysis, we present a validated reference set of 159 genes whose regulation is strongly affected by mycobacterial infection in the three types of microarrays analyzed. Furthermore, we analyzed the separate datasets of the microarrays with special emphasis on the expression profiles of immune-related genes. Upregulated genes include many known components of the inflammatory response and several genes that have previously been implicated in the response to mycobacterial infections in cell cultures of other organisms. Different marker genes of the myeloid lineage that have been characterized in zebrafish also showed increased expression. Furthermore, the zebrafish homologs of many signal transduction genes with relationship to the immune response were induced by M. marinum infection. Future functional analysis of these genes may contribute to understanding the mechanisms of mycobacterial pathogenesis. Since a large group of genes linked to immune responses did not show altered expression in the infected animals, these results suggest specific responses in mycobacterium-induced disease.

AppaDB: an AcedB database for the nematode satellite organism Pristionchus pacificus.

Pristionchus pacificus is a free-living nematode of the Diplogastridae family and was recently developed as a satellite system in evolutionary developmental biology. AppaDB, a P.pacificus database, was created (http://appadb.eb.tuebingen. mpg.de) to integrate the genomic data of P.pacificus, comprising the physical map, genetic linkage map, EST and BAC end sequence and hybridization data. This developing database serves as a repository to search and find any information regarding physical contigs or genetic markers required for mapping of mutants. Additionally, it provides a platform for the Caenorhabditis elegans community to compare nematode genetic data in an evolutionary perspective.