Development and application of in vivo expression technology (IVET) for analysing microbial gene expression in complex environments

Establishing the mechanisms by which microbes interact with their environment, including eukaryotic hosts, is a major challenge that is essential for the economic utilisation of microbes and their products. Techniques for determining global gene expression profiles of microbes, such as microarray analyses, are often hampered by methodological restraints, particularly the recovery of bacterial transcripts (RNA) from complex mixtures and rapid degradation of RNA. A pioneering technology that avoids this problem is In Vivo Expression Technology (IVET). IVET is a 'promoter-trapping' methodology that can be used to capture nearly all bacterial promoters (genes) upregulated during a microbe-environment interaction. IVET is especially useful because there is virtually no limit to the type of environment used (examples to date include soil, oomycete, a host plant or animal) to select for active microbial promoters. Furthermore, IVET provides a powerful method to identify genes that are often overlooked during genomic annotation, and has proven to be a flexible technology that can provide even more information than identification of gene expression profiles. A derivative of IVET, termed resolvase-IVET (RIVET), can be used to provide spatio-temporal information about environment-specific gene expression. More recently, niche-specific genes captured during an IVET screen have been exploited to identify the regulatory mechanisms controlling their expression. Overall, IVET and its various spin-offs have proven to be a valuable and robust set of tools for analysing microbial gene expression in complex environments and providing new targets for biotechnological development.

Molecular comparisons of the Culex pipiens (L.) complex esterase gene amplicons

The amplification of carboxylesterase genes is a mechanism of organophosphate resistance in Culex mosquitoes. Amplified carboxylesterase genes from an insecticide resistant Culex pipiens strain collected in Cyprus were analysed and compared to other Culex amplified carboxylesterase alleles. A 12 kb section of genomic DNA containing two gene loci coding for carboxylesterase alleles A5 and B5 was cloned and sequenced. A comparison between this amplicon and one from a strain with co-amplified carboxylesterase alleles A2 and B2 revealed a number of differences. The intergenic spacer was 3.7 kb in length in the A5-B5 amplicon (2.7 kb in A2-B2) and contained putative Juan and transposable elements upstream of B5. A fragment of a gene with high homology to aldehyde oxidase was also present immediately downstream of A5. The comparison revealed no differences that would explain the successful spread of the A2-B2 amplicon worldwide whilst the A5-B5 amplicon is restricted to the Mediterranean. (C) 2004 Elsevier Ltd. All rights reserved.

Use of global gene expression patterns in mechanistic studies of oestrogen action in MCF7 human breast cancer cells

Over the years, the MCF7 human breast cancer cell line has provided a model system for the study of cellular and molecular mechanisms in oestrogen regulation of cell proliferation and in progression to oestrogen and antioestrogen independent growth. Global gene expression profiling has shown that oestrogen action in MCF7 cells involves the coordinated regulation of hundreds of genes across a wide range of functional groupings and that more genes are down regulated than upregulated. Adaptation to long-term oestrogen deprivation, which results in loss of oestrogen-responsive growth, involves alterations to gene patterns not only at early time points (0-4 weeks) but continuing through to later times (20-55 weeks), and even involves alterations to patterns of oestrogen-regulated gene expression. Only 48% of the genes which were regulated >= 2-fold by oestradiol in oestrogen-responsive cells retained this responsiveness after long-term oestrogen deprivation but other genes developed de novo oestrogen regulation. Long-term exposure to fulvestrant, which resulted in loss of growth inhibition by the antioestrogen, resulted in some very large fold changes in gene expression up to 10,000-fold. Comparison of gene profiles produced by environmental chemicals with oestrogenic properties showed that each ligand gave its own unique expression profile which suggests that environmental oestrogens entering the human breast may give rise to a more complex web of interference in cell function than simply mimicking oestrogen action at inappropriate times. (C) 2009 Elsevier Ltd. All rights reserved.

Evidence for complex multigenic inheritance of radiation AML susceptibility in mice revealed using a surrogate phenotypic assay

The mapping of genes which affect individual cancer risk is an important but complex challenge. A surrogate assay of susceptibility to radiation-induced acute myeloid leukaemia (AML) in the mouse based on chromosomal radiosensitivity has been developed and validated. This assay was applied to the mapping of radiation-induced AML risk modifier loci by association with microsatellite markers. A region on chromosome (chr) 18 with strong association is identified and confirmed by backcross analysis. Additional loci on chrs 8 and 13 show significant association. A key candidate gene Rbbp8 on chr18 is identified. Rbbp8 is shown to be upregulated in response to X-irradiation in the AML sensitive CBA strain but not AML resistant C57BL/6 strain. This study demonstrates the strength of utilizing surrogate endpoints of cancer susceptibility in the mapping of mouse loci and identifies additional loci that may affect radiation cancer risk.

Restricted expression of NADPH oxidase/peroxidase gene (Duox) in zone VII of the ascidian endostyle

The ascidian Ciona intestinalis, a marine invertebrate chordate, is an emerging model system for developmental and evolutionary studies. The endostyle, one of the characteristic organs of ascidians, is a pharyngeal structure with iodine-concentrating and peroxidase activities and is therefore considered to be homologous to the follicular thyroid of higher vertebrates. We have previously reported that a limited part of the endostyle (zone VII) is marked by the expression of orthologs of the thyroid peroxidase (TPO) and thyroid transcription factor-2 (TTF-2/FoxE) genes. In this study, we have identified the Ciona homolog of NADPH oxidase/peroxidase (Duox), which provides hydrogen peroxide (H2O2) for iodine metabolism by TPO in the vertebrate thyroid. Expression patterns assessed by in situ hybridization have revealed that Ciona Duox (Ci-Duox) is predominantly expressed in the dorsal part of zone VII of the endostyle. Furthermore, two-color fluorescent in situ hybridization with Ci-Duox and Ciona TPO (CiTPO) has revealed that the ventral boundary of the Ci-Duox domain of expression is more dorsal than that of CiTPO. We have also characterized several genes, such as Ci-Fgf8/17/18, 5HT7, and Ci-NK4, which are predominantly expressed in the ventral part of zone VII, in a region complementary to the Ci-Duox expression domain. These observations suggest that, at the molecular level, zone VII has a complex organization that might have some impact on the specification of cell types and functions in this thyroid-equivalent element of the ascidian endostyle.

Global iron-dependent gene regulation in Escherichia coli - A new mechanism for iron homeostasis

Organisms generally respond to iron deficiency by increasing their capacity to take up iron and by consuming intracellular iron stores. Escherichia coli, in which iron metabolism is particularly well understood, contains at least 7 iron-acquisition systems encoded by 35 iron-repressed genes. This Fe-dependent repression is mediated by a transcriptional repressor, Fur ( ferric uptake regulation), which also controls genes involved in other processes such as iron storage, the Tricarboxylic Acid Cycle, pathogenicity, and redox-stress resistance. Our macroarray-based global analysis of iron- and Fur-dependent gene expression in E. coli has revealed several novel Fur-repressed genes likely to specify at least three additional iron- transport pathways. Interestingly, a large group of energy metabolism genes was found to be iron and Fur induced. Many of these genes encode iron- rich respiratory complexes. This iron- and Fur-dependent regulation appears to represent a novel iron-homeostatic mechanism whereby the synthesis of many iron- containing proteins is repressed under iron- restricted conditions. This mechanism thus accounts for the low iron contents of fur mutants and explains how E. coli can modulate its iron requirements. Analysis of Fe-55-labeled E. coli proteins revealed a marked decrease in iron- protein composition for the fur mutant, and visible and EPR spectroscopy showed major reductions in cytochrome b and d levels, and in iron- sulfur cluster contents for the chelator-treated wild-type and/or fur mutant, correlating well with the array and quantitative RT-PCR data. In combination, the results provide compelling evidence for the regulation of intracellular iron consumption by the Fe2+-Fur complex.

Effect of Ginkgo biloba (EGb 761®) on differential gene expression

Supplementation of diets with plant extracts such as ginkgo biloba extract (EGb 761®) (definition see editorial) for health and prevention of degenerative diseases is popular. However, it is often difficult to analyse the biological activities of plant extracts due to their complex nature and the possible synergistic and/or antagonistic effects of their components. Genome-wide expression monitoring with high-density oligonucleotide arrays provides one way to examine the molecular targets of plant extracts and may prove a useful tool in evaluating their therapeutic claims. Here, we will briefly describe some of our work on the effect of EGb 761® on differential gene expression in relation to its potential anti-carcinogenic, photoprotective and neuromodulatory properties.

Complex history of a chromosomal paralogy region: insights from amphioxus aromatic amino acid hydroxylase genes and insulin-related genes

Aromatic amino acid hydroxylase (AAAH) genes and insulin-like genes form part of an extensive paralogy region shared by human chromosomes 11 and 12, thought to have arisen by tetraploidy in early vertebrate evolution. Cloning of a complementary DNA (cDNA) for an amphioxus (Branchiostoma floridae) hydroxylase gene (AmphiPAH) allowed us to investigate the ancestry of the human chromosome 11/12 paralogy region. Molecular phylogenetic evidence reveals that AmphiPAH is orthologous to vertebrate phenylalanine (PAH) genes; the implication is that all three vertebrate AAAH genes arose early in metazoan evolution, predating vertebrates. In contrast, our phylogenetic analysis of amphioxus and vertebrate insulin-related gene sequences is consistent with duplication of these genes during early chordate ancestry. The conclusion is that two tightly linked gene families on human chromosomes 11 and 12 were not duplicated coincidentally. We rationalize this paradox by invoking gene loss in the AAAH gene family and conclude that paralogous genes shared by paralogous chromosomes need not have identical evolutionary histories.

Endosomal proteolysis regulates calcitonin gene-related peptide responses in mesenteric arteries

Background and Purpose: Calcitonin gene‐related peptide (CGRP) is a potent vasodilator, implicated in the pathogenesis of migraine. CGRP activates a receptor complex comprising, calcitonin receptor‐like receptor (CLR) and receptor activity‐modifying protein 1 (RAMP1). In vitro studies indicate recycling of CLR•RAMP1 is regulated by degradation of CGRP in early endosomes by endothelin‐converting enzyme‐1 (ECE‐1). However, it is not known if ECE‐1 regulates the resensitization of CGRP‐induced responses in functional arterial tissue. Experimental Approach: CLR, ECE‐1a‐d and RAMP1 expression in rat mesenteric artery smooth muscle cells (RMA‐SMCs) and mesenteric arteries was analyzed by RT‐PCR and by immunofluorescence and confocal microscopy. CGRP‐induced signaling in cells was examined by measuring cAMP production and ERK activation. CGRP‐induced relaxation of arteries was measured by isometric wire myography. ECE‐1 was inhibited using the specific inhibitor, SM‐19712. Key Results: RMA‐SMCs and arteries contained mRNA for CLR, ECE‐1a‐d and RAMP1. ECE‐1 was present in early endosomes of RMA‐SMCs and in the smooth muscle layer of arteries. CGRP induced endothelium‐independent relaxation of arteries. ECE‐1 inhibition had no effect on initial CGRP‐induced responses but reduced cAMP generation in RMA‐SMCs and vasodilation in mesenteric arteries responses to subsequent CGRP challenges. Conclusions and Implications: ECE‐1 regulates the resensitization of responses to CGRP in RMA‐SMCs and mesenteric arteries. CGRP‐induced relaxation does not involve endothelium‐derived pathways. This is the first report of ECE‐1 regulating CGRP responses in SMCs and arteries. ECE‐1 inhibitors may attenuate an important vasodilatory pathway, implicated in primary headaches and may represent a new therapeutic approach for the treatment of migraine.

Modulation of stretch-induced myocyte remodeling and gene expression by nitric oxide: a novel role for lipoma preferred partner in myofibrillogenesis

Prolonged hemodynamic load as a result of hypertension eventually leads to maladaptive cardiac adaptation and heart failure. The signalling pathways that underlie these changes are still poorly understood. The adaptive response to mechanical load is mediated by mechanosensors which convert the mechanical stimuli into a biological response. We examined the effect of cyclic mechanical stretch on myocyte adaptation using neonatal rat ventricular myocytes with 10% (adaptive) or 20% (maladaptive) maximum strain, 1Hz for 48 hours to mimic in vivo mechanical stress. Cells were also treated with and without L-NAME, a general nitric oxide synthase (NOS) inhibitor to suppress NO production. Maladaptive 20% mechanical stretch led to a significant loss of intact sarcomeres which was rescued by LNAME (P<0.05, n≥5 cultures). We hypothesized that the mechanism was through NOinduced alteration of myocyte gene expression. L-NAME up-regulated the mechanosensing proteins Muscle LIM protein (MLP (by 100%, p<0.05, n=4 cultures)) and lipoma preferred partner, a novel cardiac protein (LPP (by 80%, p<0.05, n=4 cultures)). L-NAME also significantly altered the subcellular localisation of LPP and MLP in a manner that favoured growth and adaptation. These findings suggest that NO participates in stretch-mediated adaptation. The use of isoform selective NOS inhibitors indicated a complex interaction between iNOS and nNOS isoforms regulate gene expression. LPP knockdown by siRNA led to formation of α-actinin aggregates and Z-bodies showing that myofibrillogenesis was impaired. There was an up-regulation of E3 ubiquitin ligase (MUL1) by 75% (P<0.05, n=5 cultures). This indicates that NO contributes to stretch-mediated adaptation via the upregulation of proteins associated mechansensing and myofibrillogenesis, thereby presenting potential therapeutic targets during the progression of heart failure. Keywords: Mechanotransduction, heart failure, stretch, heart, hypertrophy

Regulatory hotspots are associated with plant gene expression under varying soil phosphorus supply in Brassica rapa

Gene expression is a quantitative trait that can be mapped genetically in structured populations to identify expression quantitative trait loci (eQTL). Genes and regulatory networks underlying complex traits can subsequently be inferred. Using a recently released genome sequence, we have defined cis- and trans-eQTL and their environmental response to low phosphorus (P) availability within a complex plant genome and found hotspots of trans-eQTL within the genome. Interval mapping, using P supply as a covariate, revealed 18,876 eQTL. trans-eQTL hotspots occurred on chromosomes A06 and A01 within Brassica rapa; these were enriched with P metabolism-related Gene Ontology terms (A06) as well as chloroplast-and photosynthesis-related terms (A01). We have also attributed heritability components to measures of gene expression across environments, allowing the identification of novel gene expression markers and gene expression changes associated with low P availability. Informative gene expression markers were used to map eQTL and P use efficiency-related QTL. Genes responsive to P supply had large environmental and heritable variance components. Regulatory loci and genes associated with P use efficiency identified through eQTL analysis are potential targets for further characterization and may have potential for crop improvement.

A Brassica exon array for whole-transcript gene expression profiling

Affymetrix GeneChip (R) arrays are used widely to study transcriptional changes in response to developmental and environmental stimuli. GeneChip (R) arrays comprise multiple 25-mer oligonucleotide probes per gene and retain certain advantages over direct sequencing. For plants, there are several public GeneChip (R) arrays whose probes are localised primarily in 39 exons. Plant whole-transcript (WT) GeneChip (R) arrays are not yet publicly available, although WT resolution is needed to study complex crop genomes such as Brassica, which are typified by segmental duplications containing paralogous genes and/or allopolyploidy. Available sequence data were sampled from the Brassica A and C genomes, and 142,997 gene models identified. The assembled gene models were then used to establish a comprehensive public WT exon array for transcriptomics studies. The Affymetrix GeneChip (R) Brassica Exon 1.0 ST Array is a 5 mu M feature size array, containing 2.4 million 25-base oligonucleotide probes representing 135,201 gene models, with 15 probes per gene distributed among exons. Discrimination of the gene models was based on an E-value cut-off of 1E(-5), with <= 98 sequence identity. The 135 k Brassica Exon Array was validated by quantifying transcriptome differences between leaf and root tissue from a reference Brassica rapa line (R-o-18), and categorisation by Gene Ontologies (GO) based on gene orthology with Arabidopsis thaliana. Technical validation involved comparison of the exon array with a 60-mer array platform using the same starting RNA samples. The 135 k Brassica Exon Array is a robust platform. All data relating to the array design and probe identities are available in the public domain and are curated within the BrassEnsembl genome viewer at http://www.brassica.info/BrassEnsembl/index.html.

Modelling negative feedback networks for activating transcription factor 3 predicts a dominant role for miRNAs in immediate early gene regulation

Activating transcription factor 3 (Atf3) is rapidly and transiently upregulated in numerous systems, and is associated with various disease states. Atf3 is required for negative feedback regulation of other genes, but is itself subject to negative feedback regulation possibly by autorepression. In cardiomyocytes, Atf3 and Egr1 mRNAs are upregulated via ERK1/2 signalling and Atf3 suppresses Egr1 expression. We previously developed a mathematical model for the Atf3-Egr1 system. Here, we adjusted and extended the model to explore mechanisms of Atf3 feedback regulation. Introduction of an autorepressive loop for Atf3 tuned down its expression and inhibition of Egr1 was lost, demonstrating that negative feedback regulation of Atf3 by Atf3 itself is implausible in this context. Experimentally, signals downstream from ERK1/2 suppress Atf3 expression. Mathematical modelling indicated that this cannot occur by phosphorylation of pre-existing inhibitory transcriptional regulators because the time delay is too short. De novo synthesis of an inhibitory transcription factor (ITF) with a high affinity for the Atf3 promoter could suppress Atf3 expression, but (as with the Atf3 autorepression loop) inhibition of Egr1 was lost. Developing the model to include newly-synthesised miRNAs very efficiently terminated Atf3 protein expression and, with a 4-fold increase in the rate of degradation of mRNA from the mRNA/miRNA complex, profiles for Atf3 mRNA, Atf3 protein and Egr1 mRNA approximated to the experimental data. Combining the ITF model with that of the miRNA did not improve the profiles suggesting that miRNAs are likely to play a dominant role in switching off Atf3 expression post-induction.

Molecular species delimitation of a symbiotic fig-pollinating wasp species complex reveals extreme deviation from reciprocal partner specificity

Background: Symbiotic relationships have contributed to major evolutionary innovations, the maintenance of fundamental ecosystem functions, and the generation and maintenance of biodiversity. However, the exact nature of host/symbiont associations, which has important consequences for their dynamics, is often poorly known due to limited understanding of symbiont taxonomy and species diversity. Among classical symbioses, figs and their pollinating wasps constitute a highly diverse keystone resource in tropical forest and savannah environments. Historically, they were considered to exemplify extreme reciprocal partner specificity (one-to-one host-symbiont species relationships), but recent work has revealed several more complex cases. However, there is a striking lack of studies with the specific aims of assessing symbiont diversity and how this varies across the geographic range of the host. Results: Here, we use molecular methods to investigate cryptic diversity in the pollinating wasps of a widespread Australian fig species. Standard barcoding genes and methods were not conclusive, but incorporation of phylogenetic analyses and a recently developed nuclear barcoding gene (ITS2), gave strong support for five pollinator species. Each pollinator species was most common in a different geographic region, emphasising the importance of wide geographic sampling to uncover diversity, and the scope for divergence in coevolutionary trajectories across the host plant range. In addition, most regions had multiple coexisting pollinators, raising the question of how they coexist in apparently similar or identical resource niches. Conclusion: Our study offers a striking example of extreme deviation from reciprocal partner specificity over the full geographical range of a fig-wasp system. It also suggests that superficially identical species may be able to co-exist in a mutualistic setting albeit at different frequencies in relation to their fig host’s range. We show that comprehensive sampling and molecular taxonomic techniques may be required to uncover the true structure of cryptic biodiversity underpinning intimate ecological interactions.

Genetic variation in the oxytocin receptor (OXTR) gene is associated with Asperger Syndrome

Background Autism Spectrum Conditions (ASC) are a group of neurodevelopmental conditions characterized by impairments in communication and social interaction, alongside unusually repetitive behaviors and narrow interests. ASC are highly heritable and have complex patterns of inheritance where multiple genes are involved, alongside environmental and epigenetic factors. Asperger Syndrome (AS) is a subgroup of these conditions, where there is no history of language or cognitive delay. Animal models suggest a role for oxytocin (OXT) and oxytocin receptor (OXTR) genes in social-emotional behaviors, and several studies indicate that the oxytocin/oxytocin receptor system is altered in individuals with ASC. Previous studies have reported associations between genetic variations in the OXTR gene and ASC. Methods The present study tested for an association between nine single nucleotide polymorphisms (SNPs) in the OXTR gene and AS in 530 individuals of Caucasian origin, using SNP association test and haplotype analysis. Results There was a significant association between rs2268493 in OXTR and AS. Multiple haplotypes that include this SNP (rs2268493-rs2254298, rs2268490-rs2268493-rs2254298, rs2268493-rs2254298-rs53576, rs237885-rs2268490-rs2268493-rs2254298, rs2268490-rs2268493-rs2254298-rs53576) were also associated with AS. rs2268493 has been previously associated with ASC and putatively alters several transcription factor-binding sites and regulates chromatin states, either directly or through other variants in linkage disequilibrium (LD). Conclusions This study reports a significant association of the sequence variant rs2268493 in the OXTR gene and associated haplotypes with AS.