Adaptive sequence evolution in a color gene involved in the formation of the characteristic egg-dummies of male haplochromine cichlid fishes.

BACKGROUND: The exceptionally diverse species flocks of cichlid fishes in East Africa are prime examples of parallel adaptive radiations. About 80% of East Africa's more than 1 800 endemic cichlid species, and all species of the flocks of Lakes Victoria and Malawi, belong to a particularly rapidly evolving lineage, the haplochromines. One characteristic feature of the haplochromines is their possession of egg-dummies on the males' anal fins. These egg-spots mimic real eggs and play an important role in the mating system of these maternal mouthbrooding fish. RESULTS: Here, we show that the egg-spots of haplochromines are made up of yellow pigment cells, xanthophores, and that a gene coding for a type III receptor tyrosine kinase, colony-stimulating factor 1 receptor a (csf1ra), is expressed in egg-spot tissue. Molecular evolutionary analyses reveal that the extracellular ligand-binding and receptor-interacting domain of csf1ra underwent adaptive sequence evolution in the ancestral lineage of the haplochromines, coinciding with the emergence of egg-dummies. We also find that csf1ra is expressed in the egg-dummies of a distantly related cichlid species, the ectodine cichlid Ophthalmotilapia ventralis, in which markings with similar functions evolved on the pelvic fin in convergence to those of the haplochromines. CONCLUSION: We conclude that modifications of existing signal transduction mechanisms might have evolved in the haplochromine lineage in association with the origination of anal fin egg-dummies. That positive selection has acted during the evolution of a color gene that seems to be involved in the morphogenesis of a sexually selected trait, the egg-dummies, highlights the importance of further investigations of the comparative genomic basis of the phenotypic diversification of cichlid fishes.

Genetic diversity and ecological success of Staphylococcus aureus strains colonizing humans.

The genetic determinants and phenotypic traits which make a Staphylococcus aureus strain a successful colonizer are largely unknown. The genetic diversity and population structure of 133 S. aureus isolates from healthy, generally risk-free adult carriers were investigated using four different typing methods: multilocus sequence typing (MLST), amplified fragment length polymorphism analysis (AFLP), double-locus sequence typing (DLST), and spa typing were compared. Carriage isolates displayed great genetic diversity which could only be revealed fully by DLST. Results of AFLP and MLST were highly concordant in the delineation of genotypic clusters of closely related isolates, roughly equivalent to clonal complexes. spa typing and DLST provided considerably less phylogenetic information. The resolution of spa typing was similar to that of AFLP and inferior to that of DLST. AFLP proved to be the most universal method, combining a phylogeny-building capacity similar to that of MLST with a much higher resolution. However, it had a lower reproducibility than sequencing-based MLST, DLST, and spa typing. We found two cases of methicillin-resistant S. aureus colonization, both of which were most likely associated with employment at a health service. Of 21 genotypic clusters detected, 2 were most prevalent: cluster 45 and cluster 30 each colonized 24% of the carrier population. The number of bacteria found in nasal samples varied significantly among the clusters, but the most prevalent clusters were not particularly numerous in the nasal samples. We did not find much evidence that genotypic clusters were associated with different carrier characteristics, such as age, sex, medical conditions, or antibiotic use. This may provide empirical support for the idea that genetic clusters in bacteria are maintained in the absence of adaptation to different niches. Alternatively, carrier characteristics other than those evaluated here or factors other than human hosts may exert selective pressure maintaining genotypic clusters.

Analysis of potential transcriptomic biomarkers for Huntington's disease in peripheral blood.

Highly quantitative biomarkers of neurodegenerative disease remain an important need in the urgent quest for disease-modifying therapies. For Huntington's disease (HD), a genetic test is available (trait marker), but necessary state markers are still in development. In this report, we describe a large battery of transcriptomic tests explored as state biomarker candidates. In an attempt to exploit the known neuroinflammatory and transcriptional perturbations of disease, we measured relevant mRNAs in peripheral blood cells. The performance of these potential markers was weak overall, with only one mRNA, immediate early response 3 (IER3), showing a modest but significant increase of 32% in HD samples compared with controls. No statistically significant differences were found for any other mRNAs tested, including a panel of 12 RNA biomarkers identified in a previous report [Borovecki F, Lovrecic L, Zhou J, Jeong H, Then F, Rosas HD, Hersch SM, Hogarth P, Bouzou B, Jensen RV, et al. (2005) Proc Natl Acad Sci USA 102:11023-11028]. The present results may nonetheless inform the future design and testing of HD biomarker strategies.

In vitro susceptibility of Actinobaculum schaalii to 12 antimicrobial agents and molecular analysis of fluoroquinolone resistance.

OBJECTIVES: To assess the in vitro susceptibility of Actinobaculum schaalii to 12 antimicrobial agents as well as to dissect the genetic basis of fluoroquinolone resistance. METHODS: Forty-eight human clinical isolates of A. schaalii collected in Switzerland and France were studied. Each isolate was identified by 16S rRNA sequencing. MICs of amoxicillin, ceftriaxone, gentamicin, vancomycin, clindamycin, linezolid, ciprofloxacin, levofloxacin, moxifloxacin, co-trimoxazole, nitrofurantoin and metronidazole were determined using the Etest method. Interpretation of results was made according to EUCAST clinical breakpoints. The quinolone-resistance-determining regions (QRDRs) of gyrA and parC genes were also identified and sequence analysis was performed for all 48 strains. RESULTS: All isolates were susceptible to amoxicillin, ceftriaxone, gentamicin, clindamycin (except three), vancomycin, linezolid and nitrofurantoin, whereas 100% and 85% were resistant to ciprofloxacin/metronidazole and co-trimoxazole, respectively. Greater than or equal to 90% of isolates were susceptible to the other tested fluoroquinolones, and only one strain was highly resistant to levofloxacin (MIC ?32 mg/L) and moxifloxacin (MIC 8 mg/L). All isolates that were susceptible or low-level resistant to levofloxacin/moxifloxacin (n?=?47) showed identical GyrA and ParC amino acid QRDR sequences. In contrast, the isolate exhibiting high-level resistance to levofloxacin and moxifloxacin possessed a unique mutation in GyrA, Ala83Val (Escherichia coli numbering), whereas no mutation was present in ParC. CONCLUSIONS: When an infection caused by A. schaalii is suspected, there is a risk of clinical failure by treating with ciprofloxacin or co-trimoxazole, and ?-lactams should be preferred. In addition, acquired resistance to fluoroquinolones more active against Gram-positive bacteria is possible.

Staphylococcus aureus host range and human-bovine host shift.

Staphylococcus aureus is a major agent of bovine mastitis. The concomitant emergence of pig-associated methicillin-resistant S. aureus (MRSA) in human carriage and infection requires a reexamination of the host range and specificity of human- and cow-associated S. aureus strains, something which has not been systematically studied previously. The genetic relatedness of 500 S. aureus isolates from bovine mastitis cases, 57 isolates from nasal carriage of farmers, and 133 isolates from nonfarmers was determined by amplified fragment length polymorphism (AFLP) analysis and spa typing. Multilocus sequence typing (MLST) was conducted on a subset of isolates to match AFLP clusters with MLST clonal complexes (CCs). This data set allowed us to study host range and host specificity and to estimate the extent of bovine-to-human transmission. The genotype compositions of S. aureus isolates from farmers and nonfarmers were very similar, while the mastitis isolates were quite distinct. Overall, transmission was low, but specific genotypes did show increased cow-to-human transmission. Unexpectedly, more than one-third of mastitis isolates belonged to CC8, a lineage which has not been considered to be bovine mastitis associated, but it is well known from human carriage and infection (i.e., USA300). Despite the fact that we did detect some transmission of other genotypes from cows to farmers, no transmission of CC8 isolates to farmers was detected, except for one tentative case. This was despite the close genetic relatedness of mastitis CC8 strains to nonfarmer carriage strains. These results suggest that the emergence of the new bovine-adapted genotype was due to a recent host shift from humans to cows concurrent with a loss of the ability to colonize humans. More broadly, our results indicate that host specificity is a lineage-specific trait that can rapidly evolve.

Addressing trypsin bias in large scale (phospho)proteome analysis by size exclusion chromatography and secondary digestion of large post-trypsin peptides.

In the vast majority of bottom-up proteomics studies, protein digestion is performed using only mammalian trypsin. Although it is clearly the best enzyme available, the sole use of trypsin rarely leads to complete sequence coverage, even for abundant proteins. It is commonly assumed that this is because many tryptic peptides are either too short or too long to be identified by RPLC-MS/MS. We show through in silico analysis that 20-30% of the total sequence of three proteomes (Schizosaccharomyces pombe, Saccharomyces cerevisiae, and Homo sapiens) is expected to be covered by Large post-Trypsin Peptides (LpTPs) with M(r) above 3000 Da. We then established size exclusion chromatography to fractionate complex yeast tryptic digests into pools of peptides based on size. We found that secondary digestion of LpTPs followed by LC-MS/MS analysis leads to a significant increase in identified proteins and a 32-50% relative increase in average sequence coverage compared to trypsin digestion alone. Application of the developed strategy to analyze the phosphoproteomes of S. pombe and of a human cell line identified a significant fraction of novel phosphosites. Overall our data indicate that specific targeting of LpTPs can complement standard bottom-up workflows to reveal a largely neglected portion of the proteome.

Pseudogenes in the ENCODE regions: consensus annotation, analysis of transcription, and evolution.

Arising from either retrotransposition or genomic duplication of functional genes, pseudogenes are "genomic fossils" valuable for exploring the dynamics and evolution of genes and genomes. Pseudogene identification is an important problem in computational genomics, and is also critical for obtaining an accurate picture of a genome's structure and function. However, no consensus computational scheme for defining and detecting pseudogenes has been developed thus far. As part of the ENCyclopedia Of DNA Elements (ENCODE) project, we have compared several distinct pseudogene annotation strategies and found that different approaches and parameters often resulted in rather distinct sets of pseudogenes. We subsequently developed a consensus approach for annotating pseudogenes (derived from protein coding genes) in the ENCODE regions, resulting in 201 pseudogenes, two-thirds of which originated from retrotransposition. A survey of orthologs for these pseudogenes in 28 vertebrate genomes showed that a significant fraction ( approximately 80%) of the processed pseudogenes are primate-specific sequences, highlighting the increasing retrotransposition activity in primates. Analysis of sequence conservation and variation also demonstrated that most pseudogenes evolve neutrally, and processed pseudogenes appear to have lost their coding potential immediately or soon after their emergence. In order to explore the functional implication of pseudogene prevalence, we have extensively examined the transcriptional activity of the ENCODE pseudogenes. We performed systematic series of pseudogene-specific RACE analyses. These, together with complementary evidence derived from tiling microarrays and high throughput sequencing, demonstrated that at least a fifth of the 201 pseudogenes are transcribed in one or more cell lines or tissues.

Improved multiple-locus variable-number tandem-repeat assay for Staphylococcus aureus genotyping, providing a highly informative technique together with strong phylogenetic value.

We describe an improved multiple-locus variable-number tandem-repeat (VNTR) analysis (MLVA) scheme for genotyping Staphylococcus aureus. We compare its performance to those of multilocus sequence typing (MLST) and spa typing in a survey of 309 strains. This collection includes 87 epidemic methicillin-resistant S. aureus (MRSA) strains of the Harmony collection, 75 clinical strains representing the major MLST clonal complexes (CCs) (50 methicillin-sensitive S. aureus [MSSA] and 25 MRSA), 135 nasal carriage strains (133 MSSA and 2 MRSA), and 13 published S. aureus genome sequences. The results show excellent concordance between the techniques' results and demonstrate that the discriminatory power of MLVA is higher than those of both MLST and spa typing. Two hundred forty-two genotypes are discriminated with 14 VNTR loci (diversity index, 0.9965; 95% confidence interval, 0.9947 to 0.9984). Using a cutoff value of 45%, 21 clusters are observed, corresponding to the CCs previously defined by MLST. The variability of the different tandem repeats allows epidemiological studies, as well as follow-up of the evolution of CCs and the identification of potential ancestors. The 14 loci can conveniently be analyzed in two steps, based upon a first-line simplified assay comprising a subset of 10 loci (panel 1) and a second subset of 4 loci (panel 2) that provides higher resolution when needed. In conclusion, the MLVA scheme proposed here, in combination with available on-line genotyping databases (including http://mlva.u-psud.fr/), multiplexing, and automatic sizing, can provide a basis for almost-real-time large-scale population monitoring of S. aureus.

Phylogeography and Pleistocene refugia of the adder (Vipera berus) as inferred from mitochondrial DNA sequence data.

In order to contribute to the debate about southern glacial refugia used by temperate species and more northern refugia used by boreal or cold-temperate species, we examined the phylogeography of a widespread snake species (Vipera berus) inhabiting Europe up to the Arctic Circle. The analysis of the mitochondrial DNA (mtDNA) sequence variation in 1043 bp of the cytochrome b gene and in 918 bp of the noncoding control region was performed with phylogenetic approaches. Our results suggest that both the duplicated control region and cytochrome b evolve at a similar rate in this species. Phylogenetic analysis showed that V. berus is divided into three major mitochondrial lineages, probably resulting from an Italian, a Balkan and a Northern (from France to Russia) refugial area in Eastern Europe, near the Carpathian Mountains. In addition, the Northern clade presents an important substructure, suggesting two sequential colonization events in Europe. First, the continent was colonized from the three main refugial areas mentioned above during the Lower-Mid Pleistocene. Second, recolonization of most of Europe most likely originated from several refugia located outside of the Mediterranean peninsulas (Carpathian region, east of the Carpathians, France and possibly Hungary) during the Mid-Late Pleistocene, while populations within the Italian and Balkan Peninsulas fluctuated only slightly in distribution range, with larger lowland populations during glacial times and with refugial mountain populations during interglacials, as in the present time. The phylogeographical structure revealed in our study suggests complex recolonization dynamics of the European continent by V. berus, characterized by latitudinal as well as altitudinal range shifts, driven by both climatic changes and competition with related species.

GSVA: gene set variation analysis for microarray and RNA-seq data

Gene set enrichment (GSE) analysis is a popular framework for condensing information from gene expression profiles into a pathway or signature summary. The strengths of this approach over single gene analysis include noise and dimension reduction, as well as greater biological interpretability. As molecular profiling experiments move beyond simple case-control studies, robust and flexible GSE methodologies are needed that can model pathway activity within highly heterogeneous data sets. To address this challenge, we introduce Gene Set Variation Analysis (GSVA), a GSE method that estimates variation of pathway activity over a sample population in an unsupervised manner. We demonstrate the robustness of GSVA in a comparison with current state of the art sample-wise enrichment methods. Further, we provide examples of its utility in differential pathway activity and survival analysis. Lastly, we show how GSVA works analogously with data from both microarray and RNA-seq experiments. GSVA provides increased power to detect subtle pathway activity changes over a sample population in comparison to corresponding methods. While GSE methods are generally regarded as end points of a bioinformatic analysis, GSVA constitutes a starting point to build pathway-centric models of biology. Moreover, GSVA contributes to the current need of GSE methods for RNA-seq data. GSVA is an open source software package for R which forms part of the Bioconductor project and can be downloaded at http://www.bioconductor.org.

Study of simple sequence repeat markers from coffee expressed sequences associated to leaf miner resistance

The objective of this work was to identify expressed simple sequence repeats (SSR) markers associated to leaf miner resistance in coffee progenies. Identification of SSR markers was accomplished by directed searches on the Brazilian Coffee Expressed Sequence Tags (EST) database. Sequence analysis of 32 selected SSR loci showed that 65% repeats are of tetra-, 21% of tri- and 14% of dinucleotides. Also, expressed SSR are localized frequently in the 5'-UTR of gene transcript. Moreover, most of the genes containing SSR are associated with defense mechanisms. Polymorphisms were analyzed in progenies segregating for resistance to the leaf miner and corresponding to advanced generations of a Coffea arabica x Coffea racemosa hybrid. Frequency of SSR alleles was 2.1 per locus. However, no polymorphism associated with leaf miner resistance was identified. These results suggest that marker-assisted selection in coffee breeding should be performed on the initial cross, in which genetic variability is still significant.

No longitudinal mitochondrial DNA sequence changes in HIV-infected individuals with and without lipoatrophy.

The potential for mitochondrial (mt) DNA mutation accumulation during antiretroviral therapy (ART), and preferential accumulation in patients with lipoatrophy compared with control participants, remains controversial. We sequenced the entire mitochondrial genome, both before ART and after ART exposure, in 29 human immunodeficiency virus (HIV)-infected Swiss HIV Cohort Study participants initiating a first-line thymidine analogue-containing ART regimen. No accumulation of mtDNA mutations or deletions was detected in 13 participants who developed lipoatrophy or in 16 control participants after significant and comparable ART exposure (median duration, 3.3 and 3.7 years, respectively). In HIV-infected persons, the development of lipoatrophy is unlikely to be associated with accumulation of mtDNA mutations detectable in peripheral blood.

The genome sequence of taurine cattle: a window to ruminant biology and evolution.

To understand the biology and evolution of ruminants, the cattle genome was sequenced to about sevenfold coverage. The cattle genome contains a minimum of 22,000 genes, with a core set of 14,345 orthologs shared among seven mammalian species of which 1217 are absent or undetected in noneutherian (marsupial or monotreme) genomes. Cattle-specific evolutionary breakpoint regions in chromosomes have a higher density of segmental duplications, enrichment of repetitive elements, and species-specific variations in genes associated with lactation and immune responsiveness. Genes involved in metabolism are generally highly conserved, although five metabolic genes are deleted or extensively diverged from their human orthologs. The cattle genome sequence thus provides a resource for understanding mammalian evolution and accelerating livestock genetic improvement for milk and meat production.

Molecular analysis of a novel gene cluster encoding an insect toxin in plant-associated strains of Pseudomonas fluorescens

Pseudomonas fluorescens CHA0 and the related strain Pf-5 are well-characterized representatives of rhizosphere bacteria that have the capacity to protect crop plants from fungal root diseases, mainly by releasing a variety of exoproducts that are toxic to plant pathogenic fungi. Here, we report that the two plant-beneficial pseudomonads also exhibit potent insecticidal activity. Anti-insect activity is linked to a novel genomic locus encoding a large protein toxin termed Fit (for P. fluorescensinsecticidal toxin) that is related to the insect toxin Mcf (Makes caterpillars floppy) of the entomopathogen Photorhabdus luminescens, a mutualist of insect-invading nematodes. When injected into the haemocoel, even low doses of P. fluorescens CHA0 or Pf-5 killed larvae of the tobacco hornworm Manduca sexta and the greater wax moth Galleria mellonella. In contrast, mutants of CHA0 or Pf-5 with deletions in the Fit toxin gene were significantly less virulent to the larvae. When expressed from an inducible promoter in a non-toxic Escherichia coli host, the Fit toxin gene was sufficient to render the bacterium toxic to both insect hosts. Our findings establish the Fit gene products of P. fluorescens CHA0 and Pf-5 as potent insect toxins that define previously unappreciated anti-insect properties of these plant-colonizing bacteria

Genome sequence of Plesiomonas shigelloides strain 302-73 (serotype O1)

Plesiomonas shigelloides, the only species of the genus, is an emergent pathogenic bacterium associated with human diarrheal and extraintestinal disease. We present the whole-genome sequence analysis of the representative strain for the O1 serotype (strain 302-73), providing a tool for studying bacterial outbreaks, virulence factors, and accurate diagnostic methods.