Variation in transposable element sequence and activity in the nematode Caenorhabditis elegans

Eukaryotic genomes are replete with transposable elements. The nematode C. elegans will be the first multicellular organism to have its genome completely sequenced. This sequence will allow identification of all the transposon and transposon-related sequences from a single genome. In anticipation of the complete genome sequence I have initiated a series of analyses of sequences from the C. elegans genome database that share significant similarity to known families of transposons. Several members of known transposon families were observed along with a plethora of sequences related to these known transposons. Cladistic analyses were used to describe the relationships among transposons and transposon families. These analyses suggest that transposons in C. elegans may be found in both autonomous and nonautonomous forms. The differences between related element families lies mostly in the length of the inverted repeats and the presence of open reading frames. Differences between sequences within an element family suggest several mechanisms for generating length variation in inverted repeats. Characterization of the consequences of Tc1 insertion requires a means of detecting insertions. I describe reverse genetic methodology for identifying new transposon insertions. To study the regulation of transposon activity I focused on the tissue-specific and developmental regulation of Tc1. I identified sites that are frequent targets for Tc1 insertion. In the most dramatic example, insertion of Tc1 was detected at the same site in the unc-54 gene in nearly every animal screened. This site was previously shown to be a hotspot for germ-line insertion, although at a frequency several orders of magnitude less than the levels now detected. I believe these insertions are somatic events because they increase in frequency during development but are not transmitted to progeny based on both genetic and molecular evidence and because I detect them in animals lacking a germline. Additional sites in unc-54 and src-1, another C. elegans gene, were identified as frequent targets for insertion of Tc1; however, none are hit as frequently as the unc-54 hotspot . Somatic insertion of Tc1 depends on genetic background and may be suppressed early in development

Reordering contigs of draft genomes using the Mauve Aligner

Summary: Mauve Contig Mover provides a new method for proposing the relative order of contigs that make up a draft genome based on comparison to a complete or draft reference genome. A novel application of the Mauve aligner and viewer provides an automated reordering algorithm coupled with a powerful drill-down display allowing detailed exploration of results

Analyzing the dose-dependence of the Saccharomyces cerevisiae global transcriptional response to methyl methanesulfonate and ionizing radiation

BACKGROUND: One of the most crucial tasks for a cell to ensure its long term survival is preserving the integrity of its genetic heritage via maintenance of DNA structure and sequence. While the DNA damage response in the yeast Saccharomyces cerevisiae, a model eukaryotic organism, has been extensively studied, much remains to be elucidated about how the organism senses and responds to different types and doses of DNA damage. We have measured the global transcriptional response of S. cerevisiae to multiple doses of two representative DNA damaging agents, methyl methanesulfonate (MMS) and gamma radiation. RESULTS: Hierarchical clustering of genes with a statistically significant change in transcription illustrated the differences in the cellular responses to MMS and gamma radiation. Overall, MMS produced a larger transcriptional response than gamma radiation, and many of the genes modulated in response to MMS are involved in protein and translational regulation. Several clusters of coregulated genes whose responses varied with DNA damaging agent dose were identified. Perhaps the most interesting cluster contained four genes exhibiting biphasic induction in response to MMS dose. All of the genes (DUN1, RNR2, RNR4, and HUG1) are involved in the Mec1p kinase pathway known to respond to MMS, presumably due to stalled DNA replication forks. The biphasic responses of these genes suggest that the pathway is induced at lower levels as MMS dose increases. The genes in this cluster with a threefold or greater transcriptional response to gamma radiation all showed an increased induction with increasing gamma radiation dosage. CONCLUSION: Analyzing genome-wide transcriptional changes to multiple doses of external stresses enabled the identification of cellular responses that are modulated by magnitude of the stress, providing insights into how a cell deals with genotoxicity

Genome-wide detection and analysis of homologous recombination among sequenced strains of Escherichia coli

BACKGROUND: Comparisons of complete bacterial genomes reveal evidence of lateral transfer of DNA across otherwise clonally diverging lineages. Some lateral transfer events result in acquisition of novel genomic segments and are easily detected through genome comparison. Other more subtle lateral transfers involve homologous recombination events that result in substitution of alleles within conserved genomic regions. This type of event is observed infrequently among distantly related organisms. It is reported to be more common within species, but the frequency has been difficult to quantify since the sequences under comparison tend to have relatively few polymorphic sites. RESULTS: Here we report a genome-wide assessment of homologous recombination among a collection of six complete Escherichia coli and Shigella flexneri genome sequences. We construct a whole-genome multiple alignment and identify clusters of polymorphic sites that exhibit atypical patterns of nucleotide substitution using a random walk-based method. The analysis reveals one large segment (approximately 100 kb) and 186 smaller clusters of single base pair differences that suggest lateral exchange between lineages. These clusters include portions of 10% of the 3,100 genes conserved in six genomes. Statistical analysis of the functional roles of these genes reveals that several classes of genes are over-represented, including those involved in recombination, transport and motility. CONCLUSION: We demonstrate that intraspecific recombination in E. coli is much more common than previously appreciated and may show a bias for certain types of genes. The described method provides high-specificity, conservative inference of past recombination events

The evolution of metabolic networks of E. coli

<p>Abstract</p> <p>Background</p> <p>Despite the availability of numerous complete genome sequences from <it>E. coli </it>strains, published genome-scale metabolic models exist only for two commensal <it>E. coli </it>strains. These models have proven useful for many applications, such as engineering strains for desired product formation, and we sought to explore how constructing and evaluating additional metabolic models for <it>E. coli </it>strains could enhance these efforts.</p> <p>Results</p> <p>We used the genomic information from 16 <it>E. coli </it>strains to generate an <it>E. coli </it>pangenome metabolic network by evaluating their collective 76,990 ORFs. Each of these ORFs was assigned to one of 17,647 ortholog groups including ORFs associated with reactions in the most recent metabolic model for <it>E. coli </it>K-12. For orthologous groups that contain an ORF already represented in the MG1655 model, the gene to protein to reaction associations represented in this model could then be easily propagated to other <it>E. coli </it>strain models. All remaining orthologous groups were evaluated to see if new metabolic reactions could be added to generate a pangenome-scale metabolic model (iEco1712_pan). The pangenome model included reactions from a metabolic model update for <it>E. coli </it>K-12 MG1655 (iEco1339_MG1655) and enabled development of five additional strain-specific genome-scale metabolic models. These additional models include a second K-12 strain (iEco1335_W3110) and four pathogenic strains (two enterohemorrhagic <it>E. coli </it>O157:H7 and two uropathogens). When compared to the <it>E. coli </it>K-12 models, the metabolic models for the enterohemorrhagic (iEco1344_EDL933 and iEco1345_Sakai) and uropathogenic strains (iEco1288_CFT073 and iEco1301_UTI89) contained numerous lineage-specific gene and reaction differences. All six <it>E. coli </it>models were evaluated by comparing model predictions to carbon source utilization measurements under aerobic and anaerobic conditions, and to batch growth profiles in minimal media with 0.2% (w/v) glucose. An ancestral genome-scale metabolic model based on conserved ortholog groups in all 16 <it>E. coli </it>genomes was also constructed, reflecting the conserved ancestral core of <it>E. coli </it>metabolism (iEco1053_core). Comparative analysis of all six strain-specific <it>E. coli </it>models revealed that some of the pathogenic <it>E. coli </it>strains possess reactions in their metabolic networks enabling higher biomass yields on glucose. Finally the lineage-specific metabolic traits were compared to the ancestral core model predictions to derive new insight into the evolution of metabolism within this species.</p> <p>Conclusion</p> <p>Our findings demonstrate that a pangenome-scale metabolic model can be used to rapidly construct additional <it>E. coli </it>strain-specific models, and that quantitative models of different strains of <it>E. coli </it>can accurately predict strain-specific phenotypes. Such pangenome and strain-specific models can be further used to engineer metabolic phenotypes of interest, such as designing new industrial <it>E. coli </it>strains.</p

Gene Ontology annotation highlights shared and divergent pathogenic strategies of type III effector proteins deployed by the plant pathogen Pseudomonas syringae pv tomato DC3000 and animal pathogenic Escherichia coli strains

Genome-informed identification and characterization of Type III effector repertoires in various bacterial strains and species is revealing important insights into the critical roles that these proteins play in the pathogenic strategies of diverse bacteria. However, non-systematic discipline-specific approaches to their annotation impede analysis of the accumulating wealth of data and inhibit easy communication of findings among researchers working on different experimental systems. The development of Gene Ontology (GO) terms to capture biological processes occurring during the interaction between organisms creates a common language that facilitates cross-genome analyses. The application of these terms to annotate type III effector genes in different bacterial species – the plant pathogen Pseudomonas syringae pv tomato DC3000 and animal pathogenic strains of Escherichia coli – illustrates how GO can effectively describe fundamental similarities and differences among different gene products deployed as part of diverse pathogenic strategies. In depth descriptions of the GO annotations for P. syringae pv tomato DC3000 effector AvrPtoB and the E. coli effector Tir are described, with special emphasis given to GO capability for capturing information about interacting proteins and taxa. GO-highlighted similarities in biological process and molecular function for effectors from additional pathosystems are also discussed

Genomic Mining for Aspergillus Natural Products

SummaryThe genus Aspergillus is renowned for its ability to produce a myriad of bioactive secondary metabolites. Although the propensity of biosynthetic genes to form contiguous clusters greatly facilitates assignment of putative secondary metabolite genes in the completed Aspergillus genomes, such analysis cannot predict gene expression and, ultimately, product formation. To circumvent this deficiency, we have examined Aspergillus nidulans microarrays for expressed secondary metabolite gene clusters by using the transcriptional regulator LaeA. Deletion or overexpression of laeA clearly identified numerous secondary metabolite clusters. A gene deletion in one of the clusters eliminated the production of the antitumor compound terrequinone A, a metabolite not described, from A. nidulans. In this paper, we highlight that LaeA-based genome mining helps decipher the secondary metabolome of Aspergilli and provides an unparalleled view to assess secondary metabolism gene regulation

Evolution of the metabolic and regulatory networks associated with oxygen availability in two phytopathogenic enterobacteria

<p>Abstract</p> <p>Background</p> <p><it>Dickeya dadantii </it>and <it>Pectobacterium atrosepticum </it>are phytopathogenic enterobacteria capable of facultative anaerobic growth in a wide range of O₂concentrations found in plant and natural environments. The transcriptional response to O₂remains under-explored for these and other phytopathogenic enterobacteria although it has been well characterized for animal-associated genera including <it>Escherichia coli </it>and <it>Salmonella enterica</it>. Knowledge of the extent of conservation of the transcriptional response across orthologous genes in more distantly related species is useful to identify rates and patterns of regulon evolution. Evolutionary events such as loss and acquisition of genes by lateral transfer events along each evolutionary branch results in lineage-specific genes, some of which may have been subsequently incorporated into the O₂-responsive stimulon. Here we present a comparison of transcriptional profiles measured using densely tiled oligonucleotide arrays for two phytopathogens, <it>Dickeya dadantii </it>3937 and <it>Pectobacterium atrosepticum </it>SCRI1043, grown to mid-log phase in MOPS minimal medium (0.1% glucose) with and without O₂.</p> <p>Results</p> <p>More than 7% of the genes of each phytopathogen are differentially expressed with greater than 3-fold changes under anaerobic conditions. In addition to anaerobic metabolism genes, the O₂responsive stimulon includes a variety of virulence and pathogenicity-genes. Few of these genes overlap with orthologous genes in the anaerobic stimulon of <it>E. coli</it>. We define these as the conserved core, in which the transcriptional pattern as well as genetic architecture are well preserved. This conserved core includes previously described anaerobic metabolic pathways such as fermentation. Other components of the anaerobic stimulon show variation in genetic content, genome architecture and regulation. Notably formate metabolism, nitrate/nitrite metabolism, and fermentative butanediol production, differ between <it>E. coli </it>and the phytopathogens. Surprisingly, the overlap of the anaerobic stimulon between the phytopathogens is also relatively small considering that they are closely related, occupy similar niches and employ similar strategies to cause disease. There are cases of interesting divergences in the pattern of transcription of genes between <it>Dickeya </it>and <it>Pectobacterium </it>for virulence-associated subsystems including the type VI secretion system (T6SS), suggesting that fine-tuning of the stimulon impacts interaction with plants or competing microbes.</p> <p>Conclusions</p> <p>The small number of genes (an even smaller number if we consider operons) comprising the conserved core transcriptional response to O₂limitation demonstrates the extent of regulatory divergence prevalent in the Enterobacteriaceae. Our orthology-driven comparative transcriptomics approach indicates that the adaptive response in the eneterobacteria is a result of interaction of core (regulators) and lineage-specific (structural and regulatory) genes. Our subsystems based approach reveals that similar phenotypic outcomes are sometimes achieved by each organism using different genes and regulatory strategies.</p

Using Comparative Genomics for Inquiry-Based Learning to Dissect Virulence of Escherichia coli O157:H7 and Yersinia pestis

Genomics and bioinformatics are topics of increasing interest in undergraduate biological science curricula. Many existing exercises focus on gene annotation and analysis of a single genome. In this paper, we present two educational modules designed to enable students to learn and apply fundamental concepts in comparative genomics using examples related to bacterial pathogenesis. Students first examine alignments of genomes of Escherichia coli O157:H7 strains isolated from three food-poisoning outbreaks using the multiple-genome alignment tool Mauve. Students investigate conservation of virulence factors using the Mauve viewer and by browsing annotations available at the A Systematic Annotation Package for Community Analysis of Genomes database. In the second module, students use an alignment of five Yersinia pestis genomes to analyze single-nucleotide polymorphisms of three genes to classify strains into biovar groups. Students are then given sequences of bacterial DNA amplified from the teeth of corpses from the first and second pandemics of the bubonic plague and asked to classify these new samples. Learning-assessment results reveal student improvement in self-efficacy and content knowledge, as well as students’ ability to use BLAST to identify genomic islands and conduct analyses of virulence factors from E. coli O157:H7 or Y. pestis. Each of these educational modules offers educators new ready-to-implement resources for integrating comparative genomic topics into their curricula