Differential Gene Repertoire in Mycobacterium ulcerans Identifies Candidate Genes for Patho-Adaptation

Background

Based on large genomic sequence polymorphisms, several haplotypes belonging to two major lineages of the human pathogen Mycobacterium ulcerans could be distinguished among patient isolates from various geographic origins. However, the biological relevance of insertional/deletional diversity is not understood.

Methodology

Using comparative genomics, we have investigated the genes located in regions of difference recently identified by DNA microarray based hybridisation analysis. The analysed regions of difference comprise ∼7% of the entire M. ulcerans genome.

Principal Findings

Several different mechanisms leading to loss of functional genes were identified, ranging from pseudogenization, caused by frame shift mutations or mobile genetic element interspersing, to large sequence polymorphisms. Four hot spot regions for genetic instability were unveiled. Altogether, 229 coding sequences were found to be differentially inactivated, constituting a repertoire of coding sequence variation in the rather monomorphic M. ulcerans.

Conclusions/Significance

The differential gene inactivation patterns associated with the M. ulcerans haplotypes identified candidate genes that may confer enhanced adaptation upon ablation of expression. A number of gene conversions confined to the classical lineage may contribute to particular virulence of this group comprising isolates from Africa and Australia. Identification of this spectrum of anti-virulence gene candidates expands our understanding of the pathogenicity and ecology of the emerging infectious disease Buruli ulcer.     

Identification of G1-Regulated Genes in Normally Cycling Human Cells

Background

Obtaining synchronous cell populations is essential for cell-cycle studies. Methods such as serum withdrawal or use of drugs which block cells at specific points in the cell cycle alter cellular events upon re-entry into the cell cycle. Regulatory events occurring in early G1 phase of a new cell cycle could have been overlooked.

Methodology and Findings

We used a robotic mitotic shake-off apparatus to select cells in late mitosis for genome-wide gene expression studies. Two separate microarray experiments were conducted, one which involved isolation of RNA hourly for several hours from synchronous cell populations, and one experiment which examined gene activity every 15 minutes from late telophase of mitosis into G1 phase. To verify synchrony of the cell populations under study, we utilized methods including BrdU uptake, FACS, and microarray analyses of histone gene activity. We also examined stress response gene activity. Our analysis enabled identification of 200 early G1-regulated genes, many of which currently have unknown functions. We also confirmed the expression of a set of genes candidates (fos, atf3 and tceb) by qPCR to further validate the newly identified genes.

Conclusion and Significance

Genome-scale expression analyses of the first two hours of G1 in naturally cycling cells enabled the discovery of a unique set of G1-regulated genes, many of which currently have unknown functions, in cells progressing normally through the cell division cycle. This group of genes may contain future targets for drug development and treatment of human disease.     

First Identification of Novel NDM Carbapenemase,NDM-7, in Escherichia coli in France

Background

The NDM-1 carbapenemase has been identified in 2008 in Enterobacteriaceae. Since then, several reports have emphasized its rapid dissemination throughout the world. The spread of NDM carbapenemases involve several bla _NDM gene variants associated with various plasmids among several Gram negative species.

Methodology

A multidrug-resistant E. coli isolate recovered from urine of a patient who had travelled to Burma has been characterized genetically and biochemically.

Principal Findings

E. coli COU was resistant to all antibiotics tested except amikacin, tigecycline, fosfomycin, and chloramphenicol. Analysis of the antibiotic resistance traits identified a metallo-ß-lactamase, a novel NDM variant, NDM-7. It differs from NDM-4 by a single amino acid substitution sharing an identical extended spectrum profile towards carbapenems. The bla _NDM-7 gene was located on an untypeable conjugative plasmid and associated with a close genetic background similar to those described among the bla _NDM-1 genes. The isolate also harbours bla _CTXM-15 and bla _OXA-1 genes and belonged to ST167.

Significance

This study highlights that spread of NDM producers correspond to spread of multiple bla _NDM genes and clones and therefore will be difficult to control.     

Phenotypic and Molecular Characterization of Extended-Spectrum Beta-Lactamase-Producing Escherichia coli in Bangladesh

Background

Resistance to cephalosporins in Enterobacteriaceae is mainly due to the production of extended-spectrum beta-lactamase (ESBL). Little is known about ESBL-producing bacteria in Bangladesh. Therefore, the study presents results of phenotypic and molecular characterization of ESBL-producing Escherichia coli from hospitals in Bangladesh.

Methods

A total of 339 E. coli isolated from patients with urinary tract and wound infections attending three different medical hospitals in urban and rural areas of Bangladesh between 2003–2007 were screened for ESBL-production by the double disk diffusion test. Isolates with ESBL-phenotype were further characterized by antibiotic susceptibility testing, PCR and sequencing of different β-lactamase and virulence genes, serotyping, and XbaI-macrorestriction followed by pulsed-field gel electrophoresis (PFGE).

Results

We identified 40 E. coli with ESBL phenotype. These isolates were resistant to ceftriaxone, ceftazidime, cefotaxime, aztreonam, cefepime, and nalidixic acid but remained susceptible to imipenem. All but one isolate were additionally resistant to ciprofloxacin, and 3 isolates were resistant to cefoxitin. ESBL genes of blaCTX-M-1-group were detected in all isolates; blaTEM-type and blaOXA-1-type genes were detected in 33 (82.5%) and 19 (47.5%) isolates, respectively. Virulence genes that are present in diarrhoeagenic E. coli were not found. Class-1 integron was present in 20 (50%) isolates. All the ESBL-producing E. coli isolates harbored plasmids ranging between 1.1 and 120 MDa. PFGE-typing revealed 26 different pulsotypes, but identical pulsotype showed 6 isolates of serotype O25:H4.

Conclusion

The prevalence of multidrug-resistant ESBL-producing E. coli isolates appears to be high and the majority of the isolates were positive for bla _CTX-M. Although there was genetic heterogeneity among isolates, presence of a cluster of isolates belonging to serotype O25:H4 indicates dissemination of the pandemic uropathogenic E. coli clone in Bangladesh.     

Nitrogen stress response of a hybrid species: a gene expression study

Background and Aims

Low soil fertility limits growth and productivity in many natural and agricultural systems, where the ability to sense and respond to nutrient limitation is important for success. Helianthus anomalus is an annual sunflower of hybrid origin that is adapted to desert sand-dune substrates with lower fertility than its parental species, H. annuus and H. petiolaris. Previous studies have shown that H. anomalus has traits generally associated with adaptation to low-fertility habitats, including a lower inherent relative growth rate and longer leaf lifetime.

Methods

Here, a cDNA microarray is used to identify gene expression differences that potentially contribute to increased tolerance of low fertility of the hybrid species by comparing the nitrogen stress response of all three species with high- and low-nutrient treatments.

Key Results

Relative to the set of genes on the microarray, the genes showing differential expression in the hybrid species compared with its parents are enriched in stress-response genes, developmental genes, and genes involved in responses to biotic or abiotic stimuli. After a correction for multiple comparisons, five unique genes show a significantly different response to nitrogen limitation in H. anomalus compared with H. petiolaris and H. annuus. The Arabidopsis thaliana homologue of one of the five genes, catalase 1, has been shown to affect the timing of leaf senescence, and thus leaf lifespan.

Conclusions

The five genes identified in this analysis will be examined further as candidate genes for the adaptive stress response in H. anomalus. Genes that improve growth and productivity under nutrient stress could be used to improve crops for lower soil fertility which is common in marginal agricultural settings.     

Clustered gene expression changes flank targeted gene loci in knockout mice

Background

Gene expression profiling using microarrays is a powerful technology widely used to study regulatory networks. Profiling of mRNA levels in mutant organisms has the potential to identify genes regulated by the mutated protein.

Methodology/Principle Findings

Using tissues from multiple lines of knockout mice we have examined genome-wide changes in gene expression. We report that a significant proportion of changed genes were found near the targeted gene.

Conclusions/Significance

The apparent clustering of these genes was explained by the presence of flanking DNA from the parental ES cell. We provide recommendations for the analysis and reporting of microarray data from knockout mice     

Closely Related NDM-1-Encoding Plasmids from Escherichia coli and Klebsiella pneumoniae in Taiwan

Objective

Two plasmids carrying bla _NDM-1 isolated from carbapenem-resistant Klebsiella pneumoniae (CR-KP) and carbapenem-resistant Escherichia coli (CR-EC) were sequenced. CR-KP and CR-EC were isolated from two Taiwanese patients without travel histories.

Methods

Complete sequencing of the plasmids (pLK75 and pLK78) was conducted using a shotgun approach. Annotation of the contigs was performed using the RAST Server, followed by manual inspection and correction.

Results

These similar plasmids were obtained from two patients with overlapping stays at the same hospital. The pLK75 and pLK78 plasmids were 56,489-bp and 56,072-bp in length, respectively. Plasmid annotation revealed a common backbone similar to the IncN plasmid pR46. The regions flanking the bla _NDM-1 genes in these plasmids were very similar to plasmid pNDM-HU01 in Japan, which contains a complex class 1 integron located next to an ISCR1 element. The ISCR1 element has been suggested to provide a powerful mechanism for mobilising antibiotic resistance genes.

Conclusion

Two indigenous NDM-1-producing Enterobacteriaceae cases were identified for the first time in Taiwan, highlighting the alarming introduction of NDM-1-producing Enterobacteriaceae in this region.     

Design Parameters to Control Synthetic Gene Expression in Escherichia coli

Background

Production of proteins as therapeutic agents, research reagents and molecular tools frequently depends on expression in heterologous hosts. Synthetic genes are increasingly used for protein production because sequence information is easier to obtain than the corresponding physical DNA. Protein-coding sequences are commonly re-designed to enhance expression, but there are no experimentally supported design principles.

Principal Findings

To identify sequence features that affect protein expression we synthesized and expressed in E. coli two sets of 40 genes encoding two commercially valuable proteins, a DNA polymerase and a single chain antibody. Genes differing only in synonymous codon usage expressed protein at levels ranging from undetectable to 30% of cellular protein. Using partial least squares regression we tested the correlation of protein production levels with parameters that have been reported to affect expression. We found that the amount of protein produced in E. coli was strongly dependent on the codons used to encode a subset of amino acids. Favorable codons were predominantly those read by tRNAs that are most highly charged during amino acid starvation, not codons that are most abundant in highly expressed E. coli proteins. Finally we confirmed the validity of our models by designing, synthesizing and testing new genes using codon biases predicted to perform well.

Conclusion

The systematic analysis of gene design parameters shown in this study has allowed us to identify codon usage within a gene as a critical determinant of achievable protein expression levels in E. coli. We propose a biochemical basis for this, as well as design algorithms to ensure high protein production from synthetic genes. Replication of this methodology should allow similar design algorithms to be empirically derived for any expression system.     

Immunomodulatory Effects of Escherichia coli ATCC 25922 on Allergic Airway Inflammation in a Mouse Model

Background

Hygiene hypothesis demonstrates that the lack of microbial exposure would promote the development of allergic airway disease (AAD). Therefore, the gut microbiota, including Escherichia coli (E. coli), would probably offer a potential strategy for AAD.

Objective

To investigate whether E. coli infection is able to suppress the induction of AAD and to elucidate the underlying mechanisms.

Methods

Nonpathogenic E. coli ATCC 25922 was infected by gavage before AAD phase in three patterns: 10⁸ or 10⁶ CFU in neonates or 10⁸ CFU in adults. Then mice were sensitized and challenged with ovalbumin (OVA) to induce allergic inflammation in both the upper and lower airways. Hallmarks of AAD, in terms of eosinophil infiltration and goblet cell metaplasia in subepithelial mucosa, Th2 skewing of the immune response, and levels of T regulate cells (Tregs), were examined by histological analysis, ELISA, and flow cytometry, respectively.

Results

E. coli, especially neonatally infected with an optimal dose, attenuated allergic responses, including a decrease in nasal rubbing and sneezing, a reduction in eosinophil inflammation and goblet cell metaplasia in subepithelial mucosa, decreased serum levels of OVA-specific IgE, and reduced Th2 (IL-4) cytokines. In contrast, this effect came with an increase of Th1 (IFN-r and IL-2) cytokines, and an enhancement of IL-10-secreting Tregs in paratracheal lymph nodes (PTLN).

Conclusion

E. coli suppresses allergic responses in mice, probably via a shift from Th1 to Th2 and/or induction of Tregs. Moreover, this infection is age- and dose-dependent, which may open up novel possibilities for new therapeutic interventions.     

Genome-wide association analysis of oxidative stress resistance in Drosophila melanogaster

Background

Aerobic organisms are susceptible to damage by reactive oxygen species. Oxidative stress resistance is a quantitative trait with population variation attributable to the interplay between genetic and environmental factors. Drosophila melanogaster provides an ideal system to study the genetics of variation for resistance to oxidative stress.

Methods and Findings

We used 167 wild-derived inbred lines of the Drosophila Genetic Reference Panel for a genome-wide association study of acute oxidative stress resistance to two oxidizing agents, paraquat and menadione sodium bisulfite. We found significant genetic variation for both stressors. Single nucleotide polymorphisms (SNPs) associated with variation in oxidative stress resistance were often sex-specific and agent-dependent, with a small subset common for both sexes or treatments. Associated SNPs had moderately large effects, with an inverse relationship between effect size and allele frequency. Linear models with up to 12 SNPs explained 67–79% and 56–66% of the phenotypic variance for resistance to paraquat and menadione sodium bisulfite, respectively. Many genes implicated were novel with no known role in oxidative stress resistance. Bioinformatics analyses revealed a cellular network comprising DNA metabolism and neuronal development, consistent with targets of oxidative stress-inducing agents. We confirmed associations of seven candidate genes associated with natural variation in oxidative stress resistance through mutational analysis.

Conclusions

We identified novel candidate genes associated with variation in resistance to oxidative stress that have context-dependent effects. These results form the basis for future translational studies to identify oxidative stress susceptibility/resistance genes that are evolutionary conserved and might play a role in human disease.     

Using effective subnetworks to predict selected properties of gene networks

Background

Difficulties associated with implementing gene therapy are caused by the complexity of the underlying regulatory networks. The forms of interactions between the hundreds of genes, proteins, and metabolites in these networks are not known very accurately. An alternative approach is to limit consideration to genes on the network. Steady state measurements of these influence networks can be obtained from DNA microarray experiments. However, since they contain a large number of nodes, the computation of influence networks requires a prohibitively large set of microarray experiments. Furthermore, error estimates of the network make verifiable predictions impossible.

Methodology/Principal Findings

Here, we propose an alternative approach. Rather than attempting to derive an accurate model of the network, we ask what questions can be addressed using lower dimensional, highly simplified models. More importantly, is it possible to use such robust features in applications? We first identify a small group of genes that can be used to affect changes in other nodes of the network. The reduced effective empirical subnetwork (EES) can be computed using steady state measurements on a small number of genetically perturbed systems. We show that the EES can be used to make predictions on expression profiles of other mutants, and to compute how to implement pre-specified changes in the steady state of the underlying biological process. These assertions are verified in a synthetic influence network. We also use previously published experimental data to compute the EES associated with an oxygen deprivation network of E.coli, and use it to predict gene expression levels on a double mutant. The predictions are significantly different from the experimental results for less than of genes.

Conclusions/Significance

The constraints imposed by gene expression levels of mutants can be used to address a selected set of questions about a gene network.     

Genetic diversity of selected genes that are potentially economically important in feral sheep of New Zealand

Background

Feral sheep are considered to be a source of genetic variation that has been lost from their domestic counterparts through selection.

Methods

This study investigates variation in the genes KRTAP1-1, KRT33, ADRB3 and DQA2 in Merino-like feral sheep populations from New Zealand and its offshore islands. These genes have previously been shown to influence wool, lamb survival and animal health.

Results

All the genes were polymorphic, but no new allele was identified in the feral populations. In some of these populations, allele frequencies differed from those observed in commercial Merino sheep and other breeds found in New Zealand. Heterozygosity levels were comparable to those observed in other studies on feral sheep. Our results suggest that some of the feral populations may have been either inbred or outbred over the duration of their apparent isolation.

Conclusion

The variation described here allows us to draw some conclusions about the likely genetic origin of the populations and selective pressures that may have acted upon them, but they do not appear to be a source of new genetic material, at least for these four genes.     

Cross-platform microarray data normalisation for regulatory network inference

Background

Inferring Gene Regulatory Networks (GRNs) from time course microarray data suffers from the dimensionality problem created by the short length of available time series compared to the large number of genes in the network. To overcome this, data integration from diverse sources is mandatory. Microarray data from different sources and platforms are publicly available, but integration is not straightforward, due to platform and experimental differences.

Methods

We analyse here different normalisation approaches for microarray data integration, in the context of reverse engineering of GRN quantitative models. We introduce two preprocessing approaches based on existing normalisation techniques and provide a comprehensive comparison of normalised datasets.

Conclusions

Results identify a method based on a combination of Loess normalisation and iterative K-means as best for time series normalisation for this problem.