Analysis of a donor gene region for a variant surface glycoprotein and its expression site in African trypanosomes

African trypanosomes evade the immune response of their mammalian hosts by sequentially expressing genes for different variant surface glycoproteins (VSGs) from telomere-linked VSG expression sites. In the Trypanosoma brucei clone whose genome is being sequenced (GUTat 10.1), we show that the expressed VSG (VSG 10.1) is duplicated from a silent donor VSG located at another telomere-linked site. We have determined two 130 kb sequences representing the VSG 10.1 donor and expression sites. The telomere-linked donor VSG 10.1 resembles metacyclic VSG expression sites, and is preceded by a cluster of 35 or more tandem housekeeping genes, all of which are transcribed away from the telomere. The 45 kb telomere-linked VSG 10.1 expression site contains a promoter followed by seven expression site-associated genes (ESAGs), three pseudo ESAGs, two pseudo VSGs and VSG 10.1. The 80 kb preceding the expression site has few, if any, functional ORFs, but contains 50 bp repeats, INGI retrotransposon-like elements, and novel 4–12 kb repeats found near other telomeres. This analysis provides the first look over a 130 kb range of a telomere-linked donor VSG and its corresponding telomere-linked VSG expression site and forms the basis for studies on antigenic variation in the context of a completely sequenced genome.     

Organization of chromosome ends in the rice blast fungus, Magnaporthe oryzae

Eukaryotic pathogens of humans often evade the immune system by switching the expression of surface proteins encoded by subtelomeric gene families. To determine if plant pathogenic fungi use a similar mechanism to avoid host defenses, we sequenced the 14 chromosome ends of the rice blast pathogen, Magnaporthe oryzae. One telomere is directly joined to ribosomal RNA-encoding genes, at the end of the ~2 Mb rDNA array. Two are attached to chromosome-unique sequences, and the remainder adjoin a distinct subtelomere region, consisting of a telomere-linked RecQ-helicase (TLH) gene flanked by several blocks of tandem repeats. Unlike other microbes, M.oryzae exhibits very little gene amplification in the subtelomere regions—out of 261 predicted genes found within 100 kb of the telomeres, only four were present at more than one chromosome end. Therefore, it seems unlikely that M.oryzae uses switching mechanisms to evade host defenses. Instead, the M.oryzae telomeres have undergone frequent terminal truncation, and there is evidence of extensive ectopic recombination among transposons in these regions. We propose that the M.oryzae chromosome termini play more subtle roles in host adaptation by promoting the loss of terminally-positioned genes that tend to trigger host defenses.     

Sequence and structure of the extrachromosomal palindrome encoding the ribosomal RNA genes in Dictyostelium

Ribosomal RNAs (rRNAs) are encoded by multicopy families of identical genes. In Dictyostelium and other protists, the rDNA is carried on extrachromosomal palindromic elements that comprise up to 20% of the nuclear DNA. We present the sequence of the 88 kb Dictyostelium rDNA element, noting that the rRNA genes are likely to be the only transcribed regions. By interrogating a library of ordered YAC clones, we provide evidence for a chromosomal copy of the rDNA on chromosome 4. This locus may provide master copies for the stable transmission of the extrachromosomal elements. The extrachromosomal elements were also found to form chromosome-sized clusters of DNA within nuclei of nocodazole-treated cells arrested in mitosis. These clusters resemble true chromosomes and may allow the efficient segregation of the rDNA during mitosis. These rDNA clusters may also explain the cytological observations of a seventh chromosome in this organism.     

The complete genomic sequence of Mycoplasma penetrans,an intracellular bacterial pathogen in humans

The complete genomic sequence of an intracellular bacterial pathogen, Mycoplasma penetrans HF-2 strain, was determined. The HF-2 genome consists of a 1 358 633 bp single circular chromosome containing 1038 predicted coding sequences (CDSs), one set of rRNA genes and 30 tRNA genes. Among the 1038 CDSs, 264 predicted proteins are common to the Mycoplasmataceae sequenced thus far and 463 are M.penetrans specific. The genome contains the two-component system but lacks the essential cellular gene, uridine kinase. The relatively large genome of M.penetrans HF-2 among mycoplasma species may be accounted for by both its rich core proteome and the presence of a number of paralog families corresponding to 25.4% of all CDSs. The largest paralog family is the p35 family, which encodes surface lipoproteins including the major antigen, P35. A total of 44 genes for p35 and p35 homologs were identified and 30 of them form one large cluster in the chromosome. The genetic tree of p35 paralogs suggests the occurrence of dynamic chromosomal rearrangement in paralog formation during evolution. Thus, M.penetrans HF-2 may have acquired diverse repertoires of antigenic variation-related genes to allow its persistent infection in humans.     

Natural variation of the Y chromosome suppresses sex ratio distortion and modulates testis-specific gene expression in Drosophila simulans

X-linked sex-ratio distorters that disrupt spermatogenesis can cause a deficiency in functional Y-bearing sperm and a female-biased sex ratio. Y-linked modifiers that restore a normal sex ratio might be abundant and favored when a X-linked distorter is present. Here we investigated natural variation of Y-linked suppressors of sex-ratio in the Winters systems and the ability of these chromosomes to modulate gene expression in Drosophila simulans. Seventy-eight Y chromosomes of worldwide origin were assayed for their resistance to the X-linked sex-ratio distorter gene Dox. Y chromosome diversity caused males to sire ∼63% to ∼98% female progeny. Genome-wide gene expression analysis revealed hundreds of genes differentially expressed between isogenic males with sensitive (high sex ratio) and resistant (low sex ratio) Y chromosomes from the same population. Although the expression of about 75% of all testis-specific genes remained unchanged across Y chromosomes, a subset of post-meiotic genes was upregulated by resistant Y chromosomes. Conversely, a set of accessory gland-specific genes and mitochondrial genes were downregulated in males with resistant Y chromosomes. The D. simulans Y chromosome also modulated gene expression in XXY females in which the Y-linked protein-coding genes are not transcribed. The data suggest that the Y chromosome might exert its regulatory functions through epigenetic mechanisms that do not require the expression of protein-coding genes. The gene network that modulates sex ratio distortion by the Y chromosome is poorly understood, other than that it might include interactions with mitochondria and enriched for genes expressed in post-meiotic stages of spermatogenesis.     

New Assembly,Reannotation and Analysis of the Entamoeba histolytica Genome Reveal New Genomic Features and Protein Content Information

Background

In order to maintain genome information accurately and relevantly, original genome annotations need to be updated and evaluated regularly. Manual reannotation of genomes is important as it can significantly reduce the propagation of errors and consequently diminishes the time spent on mistaken research. For this reason, after five years from the initial submission of the Entamoeba histolytica draft genome publication, we have re-examined the original 23 Mb assembly and the annotation of the predicted genes.

Principal Findings

The evaluation of the genomic sequence led to the identification of more than one hundred artifactual tandem duplications that were eliminated by re-assembling the genome. The reannotation was done using a combination of manual and automated genome analysis. The new 20 Mb assembly contains 1,496 scaffolds and 8,201 predicted genes, of which 60% are identical to the initial annotation and the remaining 40% underwent structural changes. Functional classification of 60% of the genes was modified based on recent sequence comparisons and new experimental data. We have assigned putative function to 3,788 proteins (46% of the predicted proteome) based on the annotation of predicted gene families, and have identified 58 protein families of five or more members that share no homology with known proteins and thus could be entamoeba specific. Genome analysis also revealed new features such as the presence of segmental duplications of up to 16 kb flanked by inverted repeats, and the tight association of some gene families with transposable elements.

Significance

This new genome annotation and analysis represents a more refined and accurate blueprint of the pathogen genome, and provides an upgraded tool as reference for the study of many important aspects of E. histolytica biology, such as genome evolution and pathogenesis.     

Comparative genomic analyses of seventeen Streptococcus pneumoniae strains: insights into the pneumococcal supragenome

The distributed-genome hypothesis (DGH) states that pathogenic bacteria possess a supragenome that is much larger than the genome of any single bacterium and that these pathogens utilize genetic recombination and a large, noncore set of genes as a means of diversity generation. We sequenced the genomes of eight nasopharyngeal strains of Streptococcus pneumoniae isolated from pediatric patients with upper respiratory symptoms and performed quantitative genomic analyses among these and nine publicly available pneumococcal strains. Coding sequences from all strains were grouped into 3,170 orthologous gene clusters, of which 1,454 (46%) were conserved among all 17 strains. The majority of the gene clusters, 1,716 (54%), were not found in all strains. Genic differences per strain pair ranged from 35 to 629 orthologous clusters, with each strain's genome containing between 21 and 32% noncore genes. The distribution of the orthologous clusters per genome for the 17 strains was entered into the finite-supragenome model, which predicted that (i) the S. pneumoniae supragenome contains more than 5,000 orthologous clusters and (ii) 99% of the orthologous clusters (~3,000) that are represented in the S. pneumoniae population at frequencies of ≥0.1 can be identified if 33 representative genomes are sequenced. These extensive genic diversity data support the DGH and provide a basis for understanding the great differences in clinical phenotype associated with various pneumococcal strains. When these findings are taken together with previous studies that demonstrated the presence of a supragenome for Streptococcus agalactiae and Haemophilus influenzae, it appears that the possession of a distributed genome is a common host interaction strategy.     

Physical and Genetic Map of the Obligate Intracellular Bacterium Coxiella burnetii

Pulsed-field gel electrophoresis and PCR techniques have been used to construct a NotI macrorestriction map of the obligate intracellular bacterium Coxiella burnetii Nine Mile. The size of the chromosome has been determined to be 2,103 kb comprising 29 NotI restriction fragments. The average resolution is 72.5 kb, or about 3.5% of the genome. Experimental data support the presence of a linear chromosome. Published genes were localized on the physical map by Southern hybridization. One gene, recognized as transposable element, was found to be present in at least nine sites evenly distributed over the whole chromosome. There is only one copy of a 16S rRNA gene. The putative oriC has been located on a 27.5-kb NotI fragment. Gene organization upstream the oriC is almost identical to that of Pseudomonas putida and Bacillus subtilis, whereas gene organization downstream the oriC seems to be unique among bacteria. The physical map will be helpful in investigations of the great heterogeneity in restriction fragment length polymorphism patterns of different isolates and the great variation in genome size. The genetic map will help to determine whether gene order in different isolates is conserved.     

Genome Sequence of Cronobacter sakazakii BAA-894 and Comparative Genomic Hybridization Analysis with Other Cronobacter Species

Background

The genus Cronobacter (formerly called Enterobacter sakazakii) is composed of five species; C. sakazakii, C. malonaticus, C. turicensis, C. muytjensii, and C. dublinensis. The genus includes opportunistic human pathogens, and the first three species have been associated with neonatal infections. The most severe diseases are caused in neonates and include fatal necrotizing enterocolitis and meningitis. The genetic basis of the diversity within the genus is unknown, and few virulence traits have been identified.

Methodology/Principal Findings

We report here the first sequence of a member of this genus, C. sakazakii strain BAA-894. The genome of Cronobacter sakazakii strain BAA-894 comprises a 4.4 Mb chromosome (57% GC content) and two plasmids; 31 kb (51% GC) and 131 kb (56% GC). The genome was used to construct a 387,000 probe oligonucleotide tiling DNA microarray covering the whole genome. Comparative genomic hybridization (CGH) was undertaken on five other C. sakazakii strains, and representatives of the four other Cronobacter species. Among 4,382 annotated genes inspected in this study, about 55% of genes were common to all C. sakazakii strains and 43% were common to all Cronobacter strains, with 10–17% absence of genes.

Conclusions/Significance

CGH highlighted 15 clusters of genes in C. sakazakii BAA-894 that were divergent or absent in more than half of the tested strains; six of these are of probable prophage origin. Putative virulence factors were identified in these prophage and in other variable regions. A number of genes unique to Cronobacter species associated with neonatal infections (C. sakazakii, C. malonaticus and C. turicensis) were identified. These included a copper and silver resistance system known to be linked to invasion of the blood-brain barrier by neonatal meningitic strains of Escherichia coli. In addition, genes encoding for multidrug efflux pumps and adhesins were identified that were unique to C. sakazakii strains from outbreaks in neonatal intensive care units.     

Complete Genome Sequence of Halalkalicoccus jeotgali B3T,an Extremely Halophilic Archaeon

Halalkalicoccus jeotgali B3^T, isolated from salt-fermented seafood from South Korea, is an extremely halophilic archaeon belonging to the family Halobacteriaceae. Here, we present the complete genome sequence of the type strain H. jeotgali B3^T (3,698,650 bp, with a G+C content of 62.5%), which consists of one chromosome and six plasmids. This is the first complete genome sequence of the Halalkalicoccus species.Extremely halophilic archaea (haloarchaea) are adapted to hypersaline environments and grow optimally in NaCl solutions of 2.6 M or higher (12). These haloarchaea are classified within the family Halobacteriaceae in the order Halobacteriales; currently, this family comprises 28 genera (3), and only 11 complete genome sequences in Halobacteriaceae have been reported. In a study of archaeal diversity in salt-fermented small shrimp or shellfish from South Korea, our laboratory isolated and characterized 5 novel, extremely halophilic archaeal strains of Halobacteriaceae. These strains included Natronococcus jeotgali (9), Halalkalicoccus jeotgali (11), Halorubrum cibi (7), Haloterrigena jeotgali (10) and Haladaptatus cibarius (8). We have now sequenced the genome of Halalkalicoccus jeotgali B3^T; genome sequencing had not been completed or initiated for any strain in this genus when our sequencing project was begun. The genus Halalkalicoccus currently contains only two species, Halalkalicoccus tibetensis (13) and H. jeotgali, and these species exhibit 98.6% gene sequence similarity in their 16S rRNA. The genome of H. jeotgali B3^T is the first of this genus to be sequenced.The complete genome sequence of H. jeotgali B3^T was determined by a whole-genome shotgun strategy using Roche 454 GS (FLX Titanium) pyrosequencing (898,168 reads totaling ∼348 Mb; ∼94-fold coverage of the genome) and a fosmid library (514 reads totaling ∼680 kb) at the Genome Resource Center, KRIBB (Korea Research Institute of Bioscience and Biotechnology). Genome sequences from pyrosequencing were processed by Roche''s software according to the manufacturer''s instructions, and sequences from the fosmid library were processed by PESTAS (6). A total of 898,196 reads were assembled using Newbler Assembler 2.3 (454 Life Science), which generated 54 large contigs (>100 bp in size) with bases having quality scores of 40 and above. The gaps between contigs were closed by primer walking and sequencing of PCR products across the gaps. The annotation was done by merging results obtained from the RAST (Rapid Annotation using Subsystem Technology) pipeline (1), Glimmer 3.02 (2), tRNAscan-SE 1.21 (5), and RNAmmer 1.2 (4).The H. jeotgali B3^T genome is 3,698,650 bases long with a 62.5% G+C content. The chromosome consists of a single circular chromosome (2,809,118 bp, with a G+C content of 65.0%) and six plasmids (406,285 bp, 55.3%; 363,534 bp, 54.2%; 44,576 bp, 58.9%; 44,459 bp, 54.9%; 23,727 bp, 47.6%; 6,951 bp, 60.6%). The genome contains 3,860 predicted coding sequences and 52 RNA genes (determined using RAST). The chromosome is predicted to contain 3,101 coding sequences with a coding intensity of 90.0%, including 47 tRNA genes, 1 5S rRNA gene, 1 16S rRNA gene, and 1 23S rRNA gene. The largest plasmid contains 466 coding sequences with a coding intensity of 81.2% and 2 tRNA genes, while the other five plasmids contain 425, 44, 48, 29, and 5 coding sequences with coding intensities of 80.2%, 84.2%, 83.0%, 69.6%, and 22.8%, respectively (determined using Glimmer3). More detailed analysis of this genome and comparative analysis with other haloarchaea will provide further insight into the genomic differences and metabolism of the extremely halophilic archaea.     

MCM-BP Is Required for Repression of Life-Cycle Specific Genes Transcribed by RNA Polymerase I in the Mammalian Infectious Form of Trypanosoma brucei

Trypanosoma brucei variant surface glycoprotein (VSG) expression is a classic example of allelic exclusion. While the genome of T. brucei contains >2,000 VSG genes and VSG pseudogenes, only one allele is expressed at the surface of each infectious trypanosome and the others are repressed. Along with recombinatorial VSG switching, allelic exclusion provides a major host evasion mechanism for trypanosomes, a phenomenon known as antigenic variation. To extend our understanding of how trypanosomes escape host immunity by differential expression of VSGs, we attempted to identify genes that contribute to VSG silencing, by performing a loss-of-silencing screen in T. brucei using a transposon-mediated random insertional mutagenesis. One identified gene, which we initially named LOS1, encodes a T. brucei MCM-Binding Protein (TbMCM-BP). Here we show that TbMCM-BP is essential for viability of infectious bloodstream-form (BF) trypanosome and is required for proper cell-cycle progression. Tandem affinity purification of TbMCM-BP followed by mass spectrometry identified four subunits (MCM4-MCM7) of the T. brucei MCM complex, a replicative helicase, and MCM8, a subunit that is uniquely co-purified with TbMCM-BP. TbMCM-BP is required not only for repression of subtelomeric VSGs but also for silencing of life-cycle specific, insect-stage genes, procyclin and procyclin-associated genes (PAGs), that are normally repressed in BF trypanosomes and are transcribed by RNA polymerase I. Our study uncovers a functional link between chromosome maintenance and RNA pol I-mediated gene silencing in T. brucei.     

A Fine-Structure Map of Spontaneous Mitotic Crossovers in the Yeast Saccharomyces cerevisiae

Homologous recombination is an important mechanism for the repair of DNA damage in mitotically dividing cells. Mitotic crossovers between homologues with heterozygous alleles can produce two homozygous daughter cells (loss of heterozygosity), whereas crossovers between repeated genes on non-homologous chromosomes can result in translocations. Using a genetic system that allows selection of daughter cells that contain the reciprocal products of mitotic crossing over, we mapped crossovers and gene conversion events at a resolution of about 4 kb in a 120-kb region of chromosome V of Saccharomyces cerevisiae. The gene conversion tracts associated with mitotic crossovers are much longer (averaging about 12 kb) than the conversion tracts associated with meiotic recombination and are non-randomly distributed along the chromosome. In addition, about 40% of the conversion events have patterns of marker segregation that are most simply explained as reflecting the repair of a chromosome that was broken in G1 of the cell cycle.     

TOPO3α Influences Antigenic Variation by Monitoring Expression-Site-Associated VSG Switching in Trypanosoma brucei

Homologous recombination (HR) mediates one of the major mechanisms of trypanosome antigenic variation by placing a different variant surface glycoprotein (VSG) gene under the control of the active expression site (ES). It is believed that the majority of VSG switching events occur by duplicative gene conversion, but only a few DNA repair genes that are central to HR have been assigned a role in this process. Gene conversion events that are associated with crossover are rarely seen in VSG switching, similar to mitotic HR. In other organisms, TOPO3α (Top3 in yeasts), a type IA topoisomerase, is part of a complex that is involved in the suppression of crossovers. We therefore asked whether a related mechanism might suppress VSG recombination. Using a set of reliable recombination and switching assays that could score individual switching mechanisms, we discovered that TOPO3α function is conserved in Trypanosoma brucei and that TOPO3α plays a critical role in antigenic switching. Switching frequency increased 10–40-fold in the absence of TOPO3α and this hyper-switching phenotype required RAD51. Moreover, the preference of 70-bp repeats for VSG recombination was mitigated, while homology regions elsewhere in ES were highly favored, in the absence of TOPO3α. Our data suggest that TOPO3α may remove undesirable recombination intermediates constantly arising between active and silent ESs, thereby balancing ES integrity against VSG recombination.     

Identification and replication of three novel myopia common susceptibility gene loci on chromosome 3q26 using linkage and linkage disequilibrium mapping

Refractive error is a highly heritable quantitative trait responsible for considerable morbidity. Following an initial genome-wide linkage study using microsatellite markers, we confirmed evidence for linkage to chromosome 3q26 and then conducted fine-scale association mapping using high-resolution linkage disequilibrium unit (LDU) maps. We used a preliminary discovery marker set across the 30-Mb region with an average SNP density of 1 SNP/15 kb (Map 1). Map 1 was divided into 51 LDU windows and additional SNPs were genotyped for six regions (Map 2) that showed preliminary evidence of multi-marker association using composite likelihood. A total of 575 cases and controls selected from the tails of the trait distribution were genotyped for the discovery sample. Malecot model estimates indicate three loci with putative common functional variants centred on MFN1 (180,566 kb; 95% confidence interval 180,505–180, 655 kb), approximately 156 kb upstream from alternate-splicing SOX2OT (182,595 kb; 95% CI 182,533–182,688 kb) and PSARL (184,386 kb; 95% CI 184,356–184,411 kb), with the loci showing modest to strong evidence of association for the Map 2 discovery samples (p<10⁻⁷, p<10⁻¹⁰, and p=0.01, respectively). Using an unselected independent sample of 1,430 individuals, results replicated for the MFN1 (p=0.006), SOX2OT (p=0.0002), and PSARL (p=0.0005) gene regions. MFN1 and PSARL both interact with OPA1 to regulate mitochondrial fusion and the inhibition of mitochondrial-led apoptosis, respectively. That two mitochondrial regulatory processes in the retina are implicated in the aetiology of myopia is surprising and is likely to provide novel insight into the molecular genetic basis of common myopia.     

Mosaic VSGs and the Scale of Trypanosoma brucei Antigenic Variation

A main determinant of prolonged Trypanosoma brucei infection and transmission and success of the parasite is the interplay between host acquired immunity and antigenic variation of the parasite variant surface glycoprotein (VSG) coat. About 0.1% of trypanosome divisions produce a switch to a different VSG through differential expression of an archive of hundreds of silent VSG genes and pseudogenes, but the patterns and extent of the trypanosome diversity phenotype, particularly in chronic infection, are unclear. We applied longitudinal VSG cDNA sequencing to estimate variant richness and test whether pseudogenes contribute to antigenic variation. We show that individual growth peaks can contain at least 15 distinct variants, are estimated computationally to comprise many more, and that antigenically distinct ‘mosaic’ VSGs arise from segmental gene conversion between donor VSG genes or pseudogenes. The potential for trypanosome antigenic variation is probably much greater than VSG archive size; mosaic VSGs are core to antigenic variation and chronic infection.