BINDER: computationally inferring a gene regulatory network for Mycobacterium abscessus.

BACKGROUND: Although many of the genic features in Mycobacterium abscessus have been fully validated, a comprehensive understanding of the regulatory elements remains lacking. Moreover, there is little understanding of how the organism regulates its transcriptomic profile, enabling cells to survive in hostile environments. Here, to computationally infer the gene regulatory network for Mycobacterium abscessus we propose a novel statistical computational modelling approach: BayesIan gene regulatory Networks inferreD via gene coExpression and compaRative genomics (BINDER). In tandem with derived experimental coexpression data, the property of genomic conservation is exploited to probabilistically infer a gene regulatory network in Mycobacterium abscessus.Inference on regulatory interactions is conducted by combining 'primary' and 'auxiliary' data strata. The data forming the primary and auxiliary strata are derived from RNA-seq experiments and sequence information in the primary organism Mycobacterium abscessus as well as ChIP-seq data extracted from a related proxy organism Mycobacterium tuberculosis. The primary and auxiliary data are combined in a hierarchical Bayesian framework, informing the apposite bivariate likelihood function and prior distributions respectively. The inferred relationships provide insight to regulon groupings in Mycobacterium abscessus. RESULTS: We implement BINDER on data relating to a collection of 167,280 regulator-target pairs resulting in the identification of 54 regulator-target pairs, across 5 transcription factors, for which there is strong probability of regulatory interaction. CONCLUSIONS: The inferred regulatory interactions provide insight to, and a valuable resource for further studies of, transcriptional control in Mycobacterium abscessus, and in the family of Mycobacteriaceae more generally. Further, the developed BINDER framework has broad applicability, useable in settings where computational inference of a gene regulatory network requires integration of data sources derived from both the primary organism of interest and from related proxy organisms.

Functional characterization of the Mycobacterium abscessus genome coupled with condition specific transcriptomics reveals conserved molecular strategies for host adaptation and persistence.

BACKGROUND: Mycobacterium abscessus subsp. abscessus (MAB) is a highly drug resistant mycobacterium and the most common respiratory pathogen among the rapidly growing non-tuberculous mycobacteria. MAB is also one of the most deadly of the emerging cystic fibrosis (CF) pathogens requiring prolonged treatment with multiple antibiotics. In addition to its "mycobacterial" virulence genes, the genome of MAB harbours a large accessory genome, presumably acquired via lateral gene transfer including homologs shared with the CF pathogens Pseudomonas aeruginosa and Burkholderia cepacia. While multiple genome sequences are available there is little functional genomics data available for this important pathogen. RESULTS: We report here the first multi-omics approach to characterize the primary transcriptome, coding potential and potential regulatory regions of the MAB genome utilizing differential RNA sequencing (dRNA-seq), RNA-seq, Ribosome profiling and LC-MS proteomics. In addition we attempt to address the genomes contribution to the molecular systems that underlie MAB's adaptation and persistence in the human host through an examination of MABs transcriptional response to a number of clinically relevant conditions. These include hypoxia, exposure to sub-inhibitory concentrations of antibiotics and growth in an artificial sputum designed to mimic the conditions within the cystic fibrosis lung. CONCLUSIONS: Our integrated map provides the first comprehensive view of the primary transcriptome of MAB and evidence for the translation of over one hundred new short open reading frames (sORFs). Our map will act as a resource for ongoing functional genomics characterization of MAB and our transcriptome data from clinically relevant stresses informs our understanding of MAB's adaptation to life in the CF lung. MAB's adaptation to growth in artificial CF sputum highlights shared metabolic strategies with other CF pathogens including P. aeruginosa and mirrors the transcriptional responses that lead to persistence in mycobacteria. These strategies include an increased reliance on amino acid metabolism, and fatty acid catabolism and highlights the relevance of the glyoxylate shunt to growth in the CF lung. Our data suggests that, similar to what is seen in chronically persisting P. aeruginosa, progression towards a biofilm mode of growth would play a more prominent role in a longer-term MAB infection. Finally, MAB's transcriptional response to antibiotics highlights the role of antibiotic modifications enzymes, active transport and the evolutionarily conserved WhiB7 driven antibiotic resistance regulon.

High-throughput transcriptomics reveals common and strain-specific responses of human macrophages to infection with Mycobacterium abscessus Smooth and Rough variants.

BACKGROUND: Mycobacterium abscessus (MAB) is an emerging pathogen causing pulmonary infections in those with inflammatory lung disorders, such as Cystic Fibrosis (CF), and is associated with the highest fatality rate among rapidly growing mycobacteria (RGM). Phenotypically, MAB manifests as either a Smooth (MAB-S) or a Rough (MAB-R) morphotype, which differ in their levels of cell wall glycopeptidolipids (GPLs) and in their pathogenicity in vivo. As one of the primary immune cells encountered by MAB, we sought to examine the early transcriptional events within macrophages, following infection with both MAB-S or MAB-R. RESULTS: We sampled the transcriptomes (mRNA and miRNA) of THP-1-derived macrophages infected with both MAB-R and MAB-S at 1, 4 and 24 h post-infection (hpi) using RNA-seq. A core set of 606 genes showed consistent expression profiles in response to both morphotypes, corresponding to the early transcriptional response to MAB. The core response is type I Interferon (IFN)-driven, involving the NF-κB and MAPK signaling pathways with concomitant pro-inflammatory cytokine production, and network analysis identified STAT1, EGR1, and SRC as key hub and bottleneck genes. MAB-S elicited a more robust transcriptional response at both the mRNA and miRNA levels, which was reflected in higher cytokine levels in culture supernatants. The transcriptional profiles of macrophages infected with both morphotypes were highly correlated, however, and a direct comparison identified few genes to distinguish them. Most of the induced miRNAs have previously been associated with mycobacterial infection and overall miRNA expression patterns were similarly highly correlated between the morphotypes. CONCLUSIONS: The report here details the first whole transcriptome analysis of the early macrophage response to MAB infection. The overall picture at the early stages of macrophage infection is similar to that of other mycobacteria, reflected in a core type I IFN and pro-inflammatory cytokine response. Large-scale transcriptional differences in the host response to the different MAB morphotypes are not evident in the early stages of infection, however the subset of genes with distinct expression profiles suggest potentially interesting differences in internal trafficking of MAB within macrophages.

Genome sequencing of the extinct Eurasian wild aurochs, Bos primigenius, illuminates the phylogeography and evolution of cattle.

BACKGROUND: Domestication of the now-extinct wild aurochs, Bos primigenius, gave rise to the two major domestic extant cattle taxa, B. taurus and B. indicus. While previous genetic studies have shed some light on the evolutionary relationships between European aurochs and modern cattle, important questions remain unanswered, including the phylogenetic status of aurochs, whether gene flow from aurochs into early domestic populations occurred, and which genomic regions were subject to selection processes during and after domestication. Here, we address these questions using whole-genome sequencing data generated from an approximately 6,750-year-old British aurochs bone and genome sequence data from 81 additional cattle plus genome-wide single nucleotide polymorphism data from a diverse panel of 1,225 modern animals. RESULTS: Phylogenomic analyses place the aurochs as a distinct outgroup to the domestic B. taurus lineage, supporting the predominant Near Eastern origin of European cattle. Conversely, traditional British and Irish breeds share more genetic variants with this aurochs specimen than other European populations, supporting localized gene flow from aurochs into the ancestors of modern British and Irish cattle, perhaps through purposeful restocking by early herders in Britain. Finally, the functions of genes showing evidence for positive selection in B. taurus are enriched for neurobiology, growth, metabolism and immunobiology, suggesting that these biological processes have been important in the domestication of cattle. CONCLUSIONS: This work provides important new information regarding the origins and functional evolution of modern cattle, revealing that the interface between early European domestic populations and wild aurochs was significantly more complex than previously thought.

Relaxed selection drives a noisy noncoding transcriptome in members of the Mycobacterium tuberculosis complex.

Related species are often used to understand the molecular underpinning of virulence through examination of a shared set of biological features attributable to a core genome of orthologous genes. An important but insufficiently studied issue, however, is the extent to which the regulatory architectures are similarly conserved. A small number of studies have compared the primary transcriptomes of different bacterial species, but few have compared closely related species with clearly divergent evolutionary histories. We addressed the impact of differing modes of evolution within the genus Mycobacterium through comparison of the primary transcriptome of M. marinum with that of a closely related lineage, M. bovis. Both are thought to have evolved from an ancestral generalist species, with M. bovis and other members of the M. tuberculosis complex having subsequently undergone downsizing of their genomes during the transition to obligate pathogenicity. M. marinum, in contrast, has retained a large genome, appropriate for an environmental organism, and is a broad-host-range pathogen. We also examined changes over a shorter evolutionary time period through comparison of the primary transcriptome of M. bovis with that of another member of the M. tuberculosis complex (M. tuberculosis) which possesses an almost identical genome but maintains a distinct host preference. Importance: Our comparison of the transcriptional start site (TSS) maps of M. marinum and M. bovis uncovers a pillar of conserved promoters, noncoding RNA (NCRNA), and a genome-wide signal in the -35 promoter regions of both species. We identify evolutionarily conserved transcriptional attenuation and highlight its potential contribution to multidrug resistance mediated through the transcriptional regulator whiB7. We show that a species population history is reflected in its transcriptome and posit relaxed selection as the main driver of an abundance of canonical -10 promoter sites in M. bovis relative to M. marinum. It appears that transcriptome composition in mycobacteria is driven primarily by the availability of such sites and that their frequencies diverge significantly across the mycobacterial clade. Finally, through comparison of M. bovis and M. tuberculosis, we illustrate that single nucleotide polymorphism (SNP)-driven promoter differences likely underpin many of the transcriptional differences between M. tuberculosis complex lineages.

(Non-)translational medicine: targeting bacterial RNA.

The rise and spread of antibiotic resistance is among the most severe challenges facing modern medicine. Despite this fact, attempts to develop novel classes of antibiotic have been largely unsuccessful. The traditional mechanisms by which antibiotics work are subject to relatively rapid bacterial resistance via mutation, and hence have a limited period of efficacy. One promising strategy to ameliorate this problem is to shift from the use of chemical compounds targeting protein structures and processes to a new era of RNA-based therapeutics. RNA-mediated regulation (riboregulation) has evolved naturally in bacteria and is therefore a highly efficient means by which gene expression can be manipulated. Here, we describe recent advances toward the development of effective anti-bacterial therapies, which operate through various strategies centered on RNA.

Whole-transcriptome, high-throughput RNA sequence analysis of the bovine macrophage response to Mycobacterium bovis infection in vitro.

BACKGROUND: Mycobacterium bovis, the causative agent of bovine tuberculosis, is an intracellular pathogen that can persist inside host macrophages during infection via a diverse range of mechanisms that subvert the host immune response. In the current study, we have analysed and compared the transcriptomes of M. bovis-infected monocyte-derived macrophages (MDM) purified from six Holstein-Friesian females with the transcriptomes of non-infected control MDM from the same animals over a 24 h period using strand-specific RNA sequencing (RNA-seq). In addition, we compare gene expression profiles generated using RNA-seq with those previously generated by us using the high-density Affymetrix® GeneChip® Bovine Genome Array platform from the same MDM-extracted RNA. RESULTS: A mean of 7.2 million reads from each MDM sample mapped uniquely and unambiguously to single Bos taurus reference genome locations. Analysis of these mapped reads showed 2,584 genes (1,392 upregulated; 1,192 downregulated) and 757 putative natural antisense transcripts (558 upregulated; 119 downregulated) that were differentially expressed based on sense and antisense strand data, respectively (adjusted P-value ≤ 0.05). Of the differentially expressed genes, 694 were common to both the sense and antisense data sets, with the direction of expression (i.e. up- or downregulation) positively correlated for 693 genes and negatively correlated for the remaining gene. Gene ontology analysis of the differentially expressed genes revealed an enrichment of immune, apoptotic and cell signalling genes. Notably, the number of differentially expressed genes identified from RNA-seq sense strand analysis was greater than the number of differentially expressed genes detected from microarray analysis (2,584 genes versus 2,015 genes). Furthermore, our data reveal a greater dynamic range in the detection and quantification of gene transcripts for RNA-seq compared to microarray technology. CONCLUSIONS: This study highlights the value of RNA-seq in identifying novel immunomodulatory mechanisms that underlie host-mycobacterial pathogen interactions during infection, including possible complex post-transcriptional regulation of host gene expression involving antisense RNA.

Genome of Acanthamoeba castellanii highlights extensive lateral gene transfer and early evolution of tyrosine kinase signaling.

BACKGROUND: The Amoebozoa constitute one of the primary divisions of eukaryotes, encompassing taxa of both biomedical and evolutionary importance, yet its genomic diversity remains largely unsampled. Here we present an analysis of a whole genome assembly of Acanthamoeba castellanii (Ac) the first representative from a solitary free-living amoebozoan. RESULTS: Ac encodes 15,455 compact intron-rich genes, a significant number of which are predicted to have arisen through inter-kingdom lateral gene transfer (LGT). A majority of the LGT candidates have undergone a substantial degree of intronization and Ac appears to have incorporated them into established transcriptional programs. Ac manifests a complex signaling and cell communication repertoire, including a complete tyrosine kinase signaling toolkit and a comparable diversity of predicted extracellular receptors to that found in the facultatively multicellular dictyostelids. An important environmental host of a diverse range of bacteria and viruses, Ac utilizes a diverse repertoire of predicted pattern recognition receptors, many with predicted orthologous functions in the innate immune systems of higher organisms. CONCLUSIONS: Our analysis highlights the important role of LGT in the biology of Ac and in the diversification of microbial eukaryotes. The early evolution of a key signaling facility implicated in the evolution of metazoan multicellularity strongly argues for its emergence early in the Unikont lineage. Overall, the availability of an Ac genome should aid in deciphering the biology of the Amoebozoa and facilitate functional genomic studies in this important model organism and environmental host.

Polycomb PHF19 binds H3K36me3 and recruits PRC2 and demethylase NO66 to embryonic stem cell genes during differentiation.

Polycomb group proteins are repressive chromatin modifiers with essential roles in metazoan development, cellular differentiation and cell fate maintenance. How Polycomb proteins access active chromatin to confer transcriptional silencing during lineage transitions remains unclear. Here we show that the Polycomb repressive complex 2 (PRC2) component PHF19 binds trimethylated histone H3 Lys36 (H3K36me3), a mark of active chromatin, via its Tudor domain. PHF19 associates with the H3K36me3 demethylase NO66, and it is required to recruit the PRC2 complex and NO66 to stem cell genes during differentiation, leading to PRC2-mediated trimethylation of histone H3 Lys27 (H3K27), loss of H3K36me3 and transcriptional silencing. We propose a model whereby PHF19 functions during mouse embryonic stem cell differentiation to transiently bind the H3K36me3 mark via its Tudor domain, forming essential contact points that allow recruitment of PRC2 and H3K36me3 demethylase activity to active gene loci during their transition to a Polycomb-repressed state.

Molecular aspects of mucin biosynthesis and mucus formation in the bovine cervix during the periestrous period.

Mucus within the cervical canal represents a hormonally regulated barrier that reconciles the need to exclude the vaginal microflora from the uterus during progesterone dominance, while permitting sperm transport at estrus. Its characteristics change during the estrous cycle to facilitate these competing functional requirements. Hydrated mucin glycoproteins synthesized by the endocervical epithelium form the molecular scaffold of this mucus. This study uses the bovine cervix as a model to examine functional groups of genes related to mucin biosynthesis and mucus production over the periestrous period when functional changes in cervical barrier function are most prominent. Cervical tissue samples were collected from 30 estrus synchronized beef heifers. Animals were slaughtered in groups starting 12 h after the withdrawal of intravaginal progesterone releasing devices (controlled internal drug releases) until 7 days postonset of estrus (luteal phase). Subsequent groupings represented proestrus, early estrus, late estrus, metestrus, and finally the early luteal phase. Tissues were submitted to next generation RNA-seq transcriptome analysis. We identified 114 genes associated with biosynthesis and intracellular transport of mucins, and postsecretory modifications of cervical; 53 of these genes showed at least a twofold change in one or more experimental group in relation to onset of estrus, and the differences between groups were significant (P < 0.05). The majority of these genes showed the greatest alteration in their expression in the 48 h postestrus and luteal phase groups.

Methicillin resistance alters the biofilm phenotype and attenuates virulence in Staphylococcus aureus device-associated infections.

Clinical isolates of Staphylococcus aureus can express biofilm phenotypes promoted by the major cell wall autolysin and the fibronectin-binding proteins or the icaADBC-encoded polysaccharide intercellular adhesin/poly-N-acetylglucosamine (PIA/PNAG). Biofilm production in methicillin-susceptible S. aureus (MSSA) strains is typically dependent on PIA/PNAG whereas methicillin-resistant isolates express an Atl/FnBP-mediated biofilm phenotype suggesting a relationship between susceptibility to ß-lactam antibiotics and biofilm. By introducing the methicillin resistance gene mecA into the PNAG-producing laboratory strain 8325-4 we generated a heterogeneously resistant (HeR) strain, from which a homogeneous, high-level resistant (HoR) derivative was isolated following exposure to oxacillin. The HoR phenotype was associated with a R602H substitution in the DHHA1 domain of GdpP, a recently identified c-di-AMP phosphodiesterase with roles in resistance/tolerance to ß-lactam antibiotics and cell envelope stress. Transcription of icaADBC and PNAG production were impaired in the 8325-4 HoR derivative, which instead produced a proteinaceous biofilm that was significantly inhibited by antibodies against the mecA-encoded penicillin binding protein 2a (PBP2a). Conversely excision of the SCCmec element in the MRSA strain BH1CC resulted in oxacillin susceptibility and reduced biofilm production, both of which were complemented by mecA alone. Transcriptional activity of the accessory gene regulator locus was also repressed in the 8325-4 HoR strain, which in turn was accompanied by reduced protease production and significantly reduced virulence in a mouse model of device infection. Thus, homogeneous methicillin resistance has the potential to affect agr- and icaADBC-mediated phenotypes, including altered biofilm expression and virulence, which together are consistent with the adaptation of healthcare-associated MRSA strains to the antibiotic-rich hospital environment in which they are frequently responsible for device-related infections in immuno-compromised patients.

Sequencing illustrates the transcriptional response of Legionella pneumophila during infection and identifies seventy novel small non-coding RNAs.

Second generation sequencing has prompted a number of groups to re-interrogate the transcriptomes of several bacterial and archaeal species. One of the central findings has been the identification of complex networks of small non-coding RNAs that play central roles in transcriptional regulation in all growth conditions and for the pathogen's interaction with and survival within host cells. Legionella pneumophila is a gram-negative facultative intracellular human pathogen with a distinct biphasic lifestyle. One of its primary environmental hosts in the free-living amoeba Acanthamoeba castellanii and its infection by L. pneumophila mimics that seen in human macrophages. Here we present analysis of strand specific sequencing of the transcriptional response of L. pneumophila during exponential and post-exponential broth growth and during the replicative and transmissive phase of infection inside A. castellanii. We extend previous microarray based studies as well as uncovering evidence of a complex regulatory architecture underpinned by numerous non-coding RNAs. Over seventy new non-coding RNAs could be identified; many of them appear to be strain specific and in configurations not previously reported. We discover a family of non-coding RNAs preferentially expressed during infection conditions and identify a second copy of 6S RNA in L. pneumophila. We show that the newly discovered putative 6S RNA as well as a number of other non-coding RNAs show evidence for antisense transcription. The nature and extent of the non-coding RNAs and their expression patterns suggests that these may well play central roles in the regulation of Legionella spp. specific traits and offer clues as to how L. pneumophila adapts to its intracellular niche. The expression profiles outlined in the study have been deposited into Genbank's Gene Expression Omnibus (GEO) database under the series accession GSE27232.

Characterization of the equine skeletal muscle transcriptome identifies novel functional responses to exercise training.

BACKGROUND: Digital gene expression profiling was used to characterize the assembly of genes expressed in equine skeletal muscle and to identify the subset of genes that were differentially expressed following a ten-month period of exercise training. The study cohort comprised seven Thoroughbred racehorses from a single training yard. Skeletal muscle biopsies were collected at rest from the gluteus medius at two time points: T(1) - untrained, (9 +/- 0.5 months old) and T(2) - trained (20 +/- 0.7 months old). RESULTS: The most abundant mRNA transcripts in the muscle transcriptome were those involved in muscle contraction, aerobic respiration and mitochondrial function. A previously unreported over-representation of genes related to RNA processing, the stress response and proteolysis was observed. Following training 92 tags were differentially expressed of which 74 were annotated. Sixteen genes showed increased expression, including the mitochondrial genes ACADVL, MRPS21 and SLC25A29 encoded by the nuclear genome. Among the 58 genes with decreased expression, MSTN, a negative regulator of muscle growth, had the greatest decrease.Functional analysis of all expressed genes using FatiScan revealed an asymmetric distribution of 482 Gene Ontology (GO) groups and 18 KEGG pathways. Functional groups displaying highly significant (P < 0.0001) increased expression included mitochondrion, oxidative phosphorylation and fatty acid metabolism while functional groups with decreased expression were mainly associated with structural genes and included the sarcoplasm, laminin complex and cytoskeleton. CONCLUSION: Exercise training in Thoroughbred racehorses results in coordinate changes in the gene expression of functional groups of genes related to metabolism, oxidative phosphorylation and muscle structure.

A complete mitochondrial genome sequence from a mesolithic wild aurochs (Bos primigenius).

BACKGROUND: The derivation of domestic cattle from the extinct wild aurochs (Bos primigenius) has been well-documented by archaeological and genetic studies. Genetic studies point towards the Neolithic Near East as the centre of origin for Bos taurus, with some lines of evidence suggesting possible, albeit rare, genetic contributions from locally domesticated wild aurochsen across Eurasia. Inferences from these investigations have been based largely on the analysis of partial mitochondrial DNA sequences generated from modern animals, with limited sequence data from ancient aurochsen samples. Recent developments in DNA sequencing technologies, however, are affording new opportunities for the examination of genetic material retrieved from extinct species, providing new insight into their evolutionary history. Here we present DNA sequence analysis of the first complete mitochondrial genome (16,338 base pairs) from an archaeologically-verified and exceptionally-well preserved aurochs bone sample. METHODOLOGY: DNA extracts were generated from an aurochs humerus bone sample recovered from a cave site located in Derbyshire, England and radiocarbon-dated to 6,738+/-68 calibrated years before present. These extracts were prepared for both Sanger and next generation DNA sequencing technologies (Illumina Genome Analyzer). In total, 289.9 megabases (22.48%) of the post-filtered DNA sequences generated using the Illumina Genome Analyzer from this sample mapped with confidence to the bovine genome. A consensus B. primigenius mitochondrial genome sequence was constructed and was analysed alongside all available complete bovine mitochondrial genome sequences. CONCLUSIONS: For all nucleotide positions where both Sanger and Illumina Genome Analyzer sequencing methods gave high-confidence calls, no discrepancies were observed. Sequence analysis reveals evidence of heteroplasmy in this sample and places this mitochondrial genome sequence securely within a previously identified aurochsen haplogroup (haplogroup P), thus providing novel insights into pre-domestic patterns of variation. The high proportion of authentic, endogenous aurochs DNA preserved in this sample bodes well for future efforts to determine the complete genome sequence of a wild ancestor of domestic cattle.

The Aedes aegypti glutathione transferase family.

In this report, we describe the glutathione transferase (GST) gene family in the dengue vector Aedes aegypti and suggest a novel role for a new class of mosquito GSTs. Twenty-six GST genes are present in Ae. aegypti, two of which are alternatively spliced to give a total of 29 transcripts for cytosolic GSTs. The six classes identified in other insect species are all represented and, as in Anopheles gambiae, the majority of the mosquito GSTs belong to the insect-specific Delta and Epsilon classes with eight members each. Sixteen secure 1:1 orthologs were identified between GSTs in Ae. aegypti and An. gambiae, but only four of these have recognisable orthologs in Drosophila melanogaster. Three mosquito-specific GSTs were identified which did not belong to any previously recognised GST classes. One of these, GSTx2, has been previously implicated in conferring 1,1,1-trichloro-2,2-bis-(p-chlorophenyl)ethane (DDT) resistance in Ae. aegypti from South America. However, we found no evidence for increased levels of this GST protein in DDT/pyrethroid-resistant populations from Thailand. Furthermore, we show that the recombinant GSTX2-2 protein is unable to metabolise DDT. Interestingly, GSTX2-2 showed an affinity for hematin, and this, together with the restricted distribution of this class to haematophagous insects, may indicate a role for these enzymes in protecting mosquitoes against heme toxicity during blood feeding.

Analysis of 14 BAC sequences from the Aedes aegypti genome: a benchmark for genome annotation and assembly.

BACKGROUND: Aedes aegypti is the principal vector of yellow fever and dengue viruses throughout the tropical world. To provide a set of manually curated and annotated sequences from the Ae. aegypti genome, 14 mapped bacterial artificial chromosome (BAC) clones encompassing 1.57 Mb were sequenced, assembled and manually annotated using a combination of computational gene-finding, expressed sequence tag (EST) matches and comparative protein homology. PCR and sequencing were used to experimentally confirm expression and sequence of a subset of these transcripts. RESULTS: Of the 51 manual annotations, 50 and 43 demonstrated a high level of similarity to Anopheles gambiae and Drosophila melanogaster genes, respectively. Ten of the 12 BAC sequences with more than one annotated gene exhibited synteny with the A. gambiae genome. Putative transcripts from eight BAC clones were found in multiple copies (two copies in most cases) in the Aedes genome assembly, which point to the probable presence of haplotype polymorphisms and/or misassemblies. CONCLUSION: This study not only provides a benchmark set of manually annotated transcripts for this genome that can be used to assess the quality of the auto-annotation pipeline and the assembly, but it also looks at the effect of a high repeat content on the genome assembly and annotation pipeline.

Unique organisation of tRNA genes in Entamoeba histolytica.

The genome sequence of the protistan parasite Entamoeba histolytica HM-1:IMSS has been completed recently. Among the findings has been a unique organisation for the tRNA genes in this organism. Forty-two of the tRNA isoacceptor types are encoded in tandem arrays that vary in unit length from 490 to 1775 basepairs and contain from 1 to 5 tRNA genes. In three cases a 5S RNA gene is also present in the unit. An estimated 10% of the genome is made up of these arrays. Interspersed between RNA-encoding sequences are short tandem repeats that are polymorphic between isolates and, in some cases, within isolates. The number and organisation of tRNA genes in E. histolytica is unprecedented. In addition to encoding the tRNAs of the organism we propose that the arrays may fulfil a structural role in the genome.

Gene discovery in the Acanthamoeba castellanii genome.

Acanthamoeba castellanii is a free-living amoeba found in soil, freshwater, and marine environments and an important predator of bacteria. Acanthamoeba castellanii is also an opportunistic pathogen of clinical interest, responsible for several distinct diseases in humans. In order to provide a genomic platform for the study of this ubiquitous and important protist, we generated a sequence survey of approximately 0.5 x coverage of the genome. The data predict that A. castellanii exhibits a greater biosynthetic capacity than the free-living Dictyostelium discoideum and the parasite Entamoeba histolytica, providing an explanation for the ability of A. castellanii to inhabit a diversity of environments. Alginate lyase may provide access to bacteria within biofilms by breaking down the biofilm matrix, and polyhydroxybutyrate depolymerase may facilitate utilization of the bacterial storage compound polyhydroxybutyrate as a food source. Enzymes for the synthesis and breakdown of cellulose were identified, and they likely participate in encystation and excystation as in D. discoideum. Trehalose-6-phosphate synthase is present, suggesting that trehalose plays a role in stress adaptation. Detection and response to a number of stress conditions is likely accomplished with a large set of signal transduction histidine kinases and a set of putative receptor serine/threonine kinases similar to those found in E. histolytica. Serine, cysteine and metalloproteases were identified, some of which are likely involved in pathogenicity.

Entamoeba histolytica: observations on metabolism based on the genome sequence.

The sequencing of the genome of Entamoeba histolytica has allowed a reconstruction of its metabolic pathways, many of which are unusual for a eukaryote. Based on the genome sequence, it appears that amino acids may play a larger role than previously thought in energy metabolism, with roles in both ATP synthesis and NAD regeneration. Arginine decarboxylase may be involved in survival of E. histolytica during its passage through the stomach. The usual pyrimidine synthesis pathway is absent, but a partial pyrimidine degradation pathway could be part of a novel pyrimidine synthesis pathway. Ribonucleotide reductase was not found in the E. histolytica genome, but it was found in the close relatives Entamoeba invadens and Entamoeba moshkovskii, suggesting a recent loss from E. histolytica. The usual eukaryotic glucose transporters are not present, but members of a prokaryotic monosaccharide transporter family are present.