High-coverage sequencing and annotated assemblies of the budgerigar genome.

BACKGROUND: Parrots belong to a group of behaviorally advanced vertebrates and have an advanced ability of vocal learning relative to other vocal-learning birds. They can imitate human speech, synchronize their body movements to a rhythmic beat, and understand complex concepts of referential meaning to sounds. However, little is known about the genetics of these traits. Elucidating the genetic bases would require whole genome sequencing and a robust assembly of a parrot genome. FINDINGS: We present a genomic resource for the budgerigar, an Australian Parakeet (Melopsittacus undulatus) -- the most widely studied parrot species in neuroscience and behavior. We present genomic sequence data that includes over 300× raw read coverage from multiple sequencing technologies and chromosome optical maps from a single male animal. The reads and optical maps were used to create three hybrid assemblies representing some of the largest genomic scaffolds to date for a bird; two of which were annotated based on similarities to reference sets of non-redundant human, zebra finch and chicken proteins, and budgerigar transcriptome sequence assemblies. The sequence reads for this project were in part generated and used for both the Assemblathon 2 competition and the first de novo assembly of a giga-scale vertebrate genome utilizing PacBio single-molecule sequencing. CONCLUSIONS: Across several quality metrics, these budgerigar assemblies are comparable to or better than the chicken and zebra finch genome assemblies built from traditional Sanger sequencing reads, and are sufficient to analyze regions that are difficult to sequence and assemble, including those not yet assembled in prior bird genomes, and promoter regions of genes differentially regulated in vocal learning brain regions. This work provides valuable data and material for genome technology development and for investigating the genomics of complex behavioral traits.

Whole genome sequencing identifies circulating Beijing-lineage Mycobacterium tuberculosis strains in Guatemala and an associated urban outbreak.

Limited data are available regarding the molecular epidemiology of Mycobacterium tuberculosis (Mtb) strains circulating in Guatemala. Beijing-lineage Mtb strains have gained prevalence worldwide and are associated with increased virulence and drug resistance, but there have been only a few cases reported in Central America. Here we report the first whole genome sequencing of Central American Beijing-lineage strains of Mtb. We find that multiple Beijing-lineage strains, derived from independent founding events, are currently circulating in Guatemala, but overall still represent a relatively small proportion of disease burden. Finally, we identify a specific Beijing-lineage outbreak centered on a poor neighborhood in Guatemala City.

Genomic Sequence and Analysis of EhV-99B1, a New Coccolithovirus from the Norwegian Fjords

Coccolithoviruses are giant dsDNA viruses that infect Emiliania huxleyi, the most ubiquitous marine microalga. Here, we present the genome of the latest coccolithovirus strain to be sequenced, EhV-99B1, and compare it with two other coccolithovirus genomes (EhV-86 and EhV-163). EhV-99B1 shares a pairwise nucleotide identity of 98% with EhV-163 (the two strains were isolated from the same Norwegian fjord but in different years), and just 96.5% with EhV-86 (isolated in the same spring as EhV-99B1 but in the English Channel). We confirmed and extended the list of relevant genomic differences between these EhVs from the Norwegian fjord and EhVs from the English Channel, namely the removal/insertions of: a phosphate permease, an endonuclease, a transposase, and two specific tRNAs. As a whole, this study provided new clues and insights into the diversity and mechanisms driving the evolution of these large oceanic viruses, in particular those processes involving selfish genetic elements.

Meta-analysis of 32 genome-wide linkage studies of schizophrenia

A genome scan meta-analysis (GSMA) was carried out on 32 independent genome-wide linkage scan analyses that included 3255 pedigrees with 7413 genotyped cases affected with schizophrenia (SCZ) or related disorders. The primary GSMA divided the autosomes into 120 bins, rank-ordered the bins within each study according to the most positive linkage result in each bin, summed these ranks (weighted for study size) for each bin across studies and determined the empirical probability of a given summed rank (P-SR) by simulation. Suggestive evidence for linkage was observed in two single bins, on chromosomes 5q (142-168 Mb) and 2q (103-134 Mb). Genome-wide evidence for linkage was detected on chromosome 2q (119-152 Mb) when bin boundaries were shifted to the middle of the previous bins. The primary analysis met empirical criteria for 'aggregate' genome-wide significance, indicating that some or all of 10 bins are likely to contain loci linked to SCZ, including regions of chromosomes 1, 2q, 3q, 4q, 5q, 8p and 10q. In a secondary analysis of 22 studies of European-ancestry samples, suggestive evidence for linkage was observed on chromosome 8p (16-33 Mb). Although the newer genome-wide association methodology has greater power to detect weak associations to single common DNA sequence variants, linkage analysis can detect diverse genetic effects that segregate in families, including multiple rare variants within one locus or several weakly associated loci in the same region. Therefore, the regions supported by this meta-analysis deserve close attention in future studies. Molecular Psychiatry (2009) 14, 774-785; doi:10.1038/mp.2008.135; published online 30 December 2008

Population structure and characterization of viridans group streptococci (VGS) including Streptococcus pneumoniae isolated from adult patients with cystic fibrosis (CF)

A study was undertaken to examine the population structure of viridans group streptococci (VGS) in the sputum of adult patients with cystic fibrosis (CF). Freshly expectorated sputa (n=58) from 45 adult CF patients were examined by selective conventional culture on Mitis-Salivarius agar and yielded 190 isolates of VGS. Sequence analyses of the rpnB and 16-23S rRNA ITS genes identified these isolates to belong to 12 species of VGS and included S. anginosus, S. australis, S. cristatus, S. gordonii, S. infantis, S. mitis, S. mutans, S. oralis, S. parasanguinis, S. pneumoniae, S. salivarius and S. sanguinis. The most frequently VGS organism isolated was S. salivarius (47/190; 24.7%), followed by S. mitts (36/190; 19%), S. sanguinis (25/190; 13.2%), S. oralis (20/190; 11.0%), S. pneumoniae (19/190; 10.0%), S. parasanguinis (16/190; 8.4%), S. infantis (11/190; 5.8%), S. gordonii (7/190; 3.7%), S. anginosus (4/190; 2.1%), S. cristatus (2/190; 1.1%), S. australis (1/190; 0.5%), S. mutans (1/190; 0.5%) and S. agalactiae (1/190; 0.5%). All, but four, patients harboured at least one VGS species, which ranged from one to five streptococcal species, with a mean of 2.85 species per patient. There was no clonality at the subspecies level employing ERIC RAPD PCR. Antibiotic susceptibility was determined by Minimum Inhibitory Concentration (MIC) testing against penicillin, erythromycin and ciprofloxacin. Overall, resistance to penicillin with all VGS was 73/190 (38.4%) and 167/190 (87.9%) for erythromycin. With regard to ciprofloxacin, 27/190 (14.2%) were fully resistant, whilst a further 21/190 (11.1%) showed intermediate resistance, which equated to approximately three quarters (74.7%) of isolates being fully sensitive to this agent. In addition, as a comparator control population, we examined antibiotic susceptibility, as above, in a non-CF population comprising 12 individuals (50 VGS isolates), who were not receiving chronic antibiotics. In comparison, 8% and 38% of VGS isolates from non-CF individuals were resistant by disk susceptibility testing to penicillin and erythromycin, respectively. None of the non-CF VGS organisms were resistant to ciprofloxacin, but 42% showed intermediate resistance. (C) 2010 European Cystic Fibrosis Society. Published by Elsevier B.V. All rights reserved.

Characterisation of Genome-Wide PLZF/RARA Target Genes

The PLZF/RARA fusion protein generated by the t(11;17)(q23;q21) translocation in acute promyelocytic leukaemia (APL) is believed to act as an oncogenic transcriptional regulator recruiting epigenetic factors to genes important for its transforming potential. However, molecular mechanisms associated with PLZF/RARA-dependent leukaemogenesis still remain unclear. We searched for specific PLZF/RARA target genes by ChIP-on-chip in the haematopoietic cell line U937 conditionally expressing PLZF/RARA. By comparing bound regions found in U937 cells expressing endogenous PLZF with PLZF/RARA-induced U937 cells, we isolated specific PLZF/RARA target gene promoters. We next analysed gene expression profiles of our identified target genes in PLZF/RARA APL patients and analysed DNA sequences and epigenetic modification at PLZF/RARA binding sites. We identify 413 specific PLZF/RARA target genes including a number encoding transcription factors involved in the regulation of haematopoiesis. Among these genes, 22 were significantly down regulated in primary PLZF/RARA APL cells. In addition, repressed PLZF/RARA target genes were associated with increased levels of H3K27me3 and decreased levels of H3K9K14ac. Finally, sequence analysis of PLZF/RARA bound sequences reveals the presence of both consensus and degenerated RAREs as well as enrichment for tissue-specific transcription factor motifs, highlighting the complexity of targeting fusion protein to chromatin. Our study suggests that PLZF/RARA directly targets genes important for haematopoietic development and supports the notion that PLZF/RARA acts mainly as an epigenetic regulator of its direct target genes.

Characterization of IS2112, a new insertion sequence from Rhodococcus, and its relationship with mobile elements belonging to the IS110 family

A new insertion sequence (IS2112) was identified in the genome of the 1-haloalkane-utilizing bacterium Rhodococcus rhodochrous NCIMB 13064. The insertion element is 1415 bp long, does not contain terminal inverted repeats, and is not flanked by directly repeated sequences. IS2112 belongs to the IS110 family of transposable elements, and forms a separate subfamily, along with IS116, Two copies of IS2112 were found in R, rhodochrous NCIMB 13064 and one, two or three copies of a similar sequence were detected in five other 1-haloalkane-degrading Rhodococcus strains. There were no sequences homologous to IS2112 found in the l-haloalkane-degrading 'Pseudomonas pavonaceae' 170 and Rhodococcus sp, HA1 or in several Rhodococcus strains which do not utilize haloalkanes, IS2112 was originally found in plasmid pRTL1 of R. rhodochrous NCIMB 13064 which harbours genes encoding utilization of l-haloalkanes, and was located 5 kbp upstream of the haloalkane dehalogenase gene (dhaA), Although the second copy of IS2112 in strain NCIMB 13064 was also present on the pRTL1 plasmid, these sequences do not apparently comprise a single composite transposon encoding haloalkane utilization. An analysis of derivatives of NCIMB 13064 revealed that IS2112 was involved in genome rearrangements. IS2112 appeared to change its location as a result of transposition and as a result of other rearrangements of the NCIMB 13064 genome.

A nine-base nucleotide sequence in the porcine circovirus type 2 (PCV2) nucleocapsid gene determines viral replication and virulence

Porcine circovirus type 2 (PCV2) was retrospectively identified by serology in swine populations as an asymptomatic infection at least 25 years prior to the first reported case of PCV2-associated postweaning multisystemit wasting syndrome (PMWS). To investigate the sudden emergence of PMWS, viral sequences were amplified from frozen archived (1970-1971) porcine tissues and the complete genome of archival PCV2 was determined. The ORF1 gene product (viral DNA replicase) was homologous to contemporary PCV2 ORF1. In ORF2 (viral nucleocapsid gene) archival PCV2, a consistent linear nine-base sequence difference at base positions 1331 through 1339 was observed.The deduced amino acid sequence from these base changes alters the nucleocapsid conformation within the second immunogenic epitope from a hydrophobic (contemporary PCV2) to a hydrophilic (archival PCV2) configuration.

Detecting sequence dependent transcriptional pauses from RNA and protein number time series

Background: Evidence suggests that in prokaryotes sequence-dependent transcriptional pauses a?ect the dynamics of transcription and translation, as well as of small genetic circuits. So far, a few pause-prone sequences have been identi?ed from in vitro measurements of transcription elongation kinetics.

Results: Using a stochastic model of gene expression at the nucleotide and codon levels with realistic parameter values, we investigate three di?erent but related questions and present statistical methods for their analysis. First, we show that information from in vivo RNA and protein temporal numbers is su?cient to discriminate between models with and without a pause site in their coding sequence. Second, we demonstrate that it is possible to separate a large variety of models from each other with pauses of various durations and locations in the template by means of a hierarchical clustering and a random forest classi?er. Third, we introduce an approximate likelihood function that allows to estimate the location of a pause site.

Conclusions: This method can aid in detecting unknown pause-prone sequences from temporal measurements of RNA and protein numbers at a genome-wide scale and thus elucidate possible roles that these sequences play in the dynamics of genetic networks and phenotype.

An Expanded Multilocus Sequence Typing Scheme for Propionibacterium acnes: Investigation of ‘Pathogenic’, ‘Commensal’ and Antibiotic Resistant Strains

The Gram-positive bacterium Propionibacterium acnes is a member of the normal human skin microbiota and is associated with various infections and clinical conditions. There is tentative evidence to suggest that certain lineages may be associated with disease and others with health. We recently described a multilocus sequence typing scheme (MLST) for P. acnes based on seven housekeeping genes (http://pubmlst.org/pacnes). We now describe an expanded eight gene version based on six housekeeping genes and two ‘putative virulence’ genes (eMLST) that provides improved high resolution
typing (91eSTs from 285 isolates), and generates phylogenies congruent with those based on whole genome analysis. When compared with the nine gene MLST scheme developed at the University of Bath, UK, and utilised by researchers at Aarhus University, Denmark, the eMLST method offers greater resolution. Using the scheme, we examined 208 isolates from disparate clinical sources, and 77 isolates from healthy skin. Acne was predominately associated with type IA₁ clonal complexes CC1, CC3 and CC4; with eST1 and eST3 lineages being highly represented. In contrast, type IA₂ strains were recovered at a rate similar to type IB and II organisms. Ophthalmic infections were predominately associated with type IA₁ and IA₂ strains, while type IB and II were more frequently recovered from soft tissue and retrieved medical devices. Strains with rRNA mutations conferring resistance to antibiotics used in acne treatment were dominated by eST3, with some evidence for intercontinental spread. In contrast, despite its high association with acne, only a small number of resistant CC1 eSTs were identified. A number of eSTs were only recovered from healthy skin, particularly eSTs representing CC72 (type II) and CC77 (type III). Collectively our data lends support to the view that pathogenic versus truly commensal lineages of P. acnes may exist. This is likely to have important therapeutic and diagnostic implications.

Identification of a new mumps virus lineage by nucleotide sequence analysis of the SH gene of ten different strains

The SH gene and its flanking sequences have been analysed for 10 strains of mumps virus (MuV) and compared to 5 others. A new lineage has been identified among UK isolates. Changes in the transcription pattern could not be correlated with differences in the sequences of the F-SH and SH-HN intergenic regions of the genome.

Pan-genome analysis of the emerging foodborne pathogen Cronobacter spp. suggests a species-level bidirectional divergence driven by niche adaptation

Background: Members of the genus Cronobacter are causes of rare but severe illness in neonates and preterm infants following the ingestion of contaminated infant formula. Seven species have been described and two of the species genomes were subsequently published. In this study, we performed comparative genomics on eight strains of Cronobacter, including six that we sequenced (representing six of the seven species) and two previously published, closed genomes.

Results: We identified and characterized the features associated with the core and pan genome of the genus Cronobacter in an attempt to understand the evolution of these bacteria and the genetic content of each species. We identified 84 genomic regions that are present in two or more Cronobacter genomes, along with 45 unique genomic regions. Many potentially horizontally transferred genes, such as lysogenic prophages, were also identified. Most notable among these were several type six secretion system gene clusters, transposons that carried tellurium, copper and/or silver resistance genes, and a novel integrative conjugative element.

Conclusions: Cronobacter have diverged into two clusters, one consisting of C. dublinensis and C. muytjensii (Cdub-Cmuy) and the other comprised of C. sakazakii, C. malonaticus, C. universalis, and C. turicensis, (Csak-Cmal-Cuni-Ctur) from the most recent common ancestral species. While several genetic determinants for plant-association and human virulence could be found in the core genome of Cronobacter, the four Cdub-Cmuy clade genomes contained several accessory genomic regions important for survival in a plant-associated environmental niche, while the Csak-Cmal-Cuni-Ctur clade genomes harbored numerous virulence-related genetic traits.

Whole-mitochondrial genome sequencing in primary open-angle glaucoma using massively parallel sequencing identifies novel and known pathogenic variants

Purpose:The aim of this study was to determine whether mutations in mitochondrial DNA play a role in high-pressure primary open-angle glaucoma (OMIM 137760) by analyzing new data from massively parallel sequencing of mitochondrial DNA.
Methods:Glaucoma patients with high-tension primary open-angle glaucoma and ethnically matched and age-matched control subjects without glaucoma were recruited. The entire human mitochondrial genome was amplified in two overlapping fragments by long-range polymerase chain reaction and used as a template for massively parallel sequencing on an Ion Torrent Personal Genome Machine. All variants were confirmed by conventional Sanger sequencing.
Results:Whole-mitochondrial genome sequencing was performed in 32 patients with primary open-angle glaucoma from India (n = 16) and Ireland (n = 16). In 16 of the 32 patients with primary open-angle glaucoma (50% of cases), there were 22 mitochondrial DNA mutations consisting of 7 novel mutations and 8 previously reported disease-associated sequence variants. Eight of 22 (36.4%) of the mitochondrial DNA mutations were in complex I mitochondrial genes.
Conclusion:Massively parallel sequencing using the Ion Torrent Personal Genome Machine with confirmation by Sanger sequencing detected a pathogenic mitochondrial DNA mutation in 50% of the primary open-angle glaucoma cohort. Our findings support the emerging concept that mitochondrial dysfunction results in the development of glaucoma and, more specifically, that complex I defects play a significant role in primary open-angle glaucoma pathogenesis.

Complete nucleotide sequence of φCHU: A Luz24likevirus infecting Pseudomonas aeruginosa and displaying a unique host range

A complete nucleotide sequence of the new Pseudomonas aeruginosa Luz24likevirus phiCHU was obtained. This virus was shown to have a unique host range whereby it grew poorly on the standard laboratory strain PAO1, but infected 26 of 46 clinical isolates screened, and strains harboring IncP2 plasmid pMG53. It was demonstrated that phiCHU has single strand interruptions in its genome. Analysis of the phiCHU genome also suggested that recombination event(s) participated in the evolution of the leftmost portion of the genome, presumably encoding early genes.