Prediction of whole-genome DNA-DNA similarity, determination of G+C content and phylogenetic analysis within the family Pasteurellaceae by multilocus sequence analysis (MLSA)

Genome predictions based on selected genes would be a very welcome approach for taxonomic studies, including DNA-DNA similarity, G+C content and representative phylogeny of bacteria. At present, DNA-DNA hybridizations are still considered the gold standard in species descriptions. However, this method is time-consuming and troublesome, and datasets can vary significantly between experiments as well as between laboratories. For the same reasons, full matrix hybridizations are rarely performed, weakening the significance of the results obtained. The authors established a universal sequencing approach for the three genes recN, rpoA and thdF for the Pasteurellaceae, and determined if the sequences could be used for predicting DNA-DNA relatedness within the family. The sequence-based similarity values calculated using a previously published formula proved most useful for species and genus separation, indicating that this method provides better resolution and no experimental variation compared to hybridization. By this method, cross-comparisons within the family over species and genus borders easily become possible. The three genes also serve as an indicator of the genome G+C content of a species. A mean divergence of around 1 % was observed from the classical method, which in itself has poor reproducibility. Finally, the three genes can be used alone or in combination with already-established 16S rRNA, rpoB and infB gene-sequencing strategies in a multisequence-based phylogeny for the family Pasteurellaceae. It is proposed to use the three sequences as a taxonomic tool, replacing DNA-DNA hybridization.

Molecular characterization and phylogenetic analysis of small ruminant lentiviruses isolated from Canadian sheep and goats

Background: Small Ruminant Lentiviruses (SRLV) are widespread in Canadian sheep and goats and represent an important health issue in these animals. There is however no data about the genetic diversity of Caprine Arthritis Encephalitis Virus (CAEV) or Maedi Visna Virus (MVV) in this country. Findings: We performed a molecular and phylogenetic analysis of sheep and goat lentiviruses from a small geographic area in Canada using long sequences from the gag region of 30 infected sheep and 36 infected goats originating from 14 different flocks. Pairwise DNA distance and phylogenetic analyses revealed that all SRLV sequences obtained from sheep clustered tightly with prototypical Maedi visna sequences from America. Similarly, all SRLV strains obtained from goats clustered tightly with prototypical US CAEV-Cork strain. Conclusions: The data reported in this study suggests that Canadian and US SRLV strains share common origins. In addition, the molecular data failed to bring to light any evidence of past cross species transmission between sheep and goats, which is consistent with the type of farming practiced in this part of the country where single species flocks predominate and where opportunities of cross species transmissions are proportionately low.

Phylogenetic analysis of the gag region encoding the matrix protein of small ruminant lentiviruses: comparative analysis and molecular epidemiological applications

Little sequence information exists on the matrix-protein (MA) encoding region of small ruminant lentiviruses (SRLV). Fifty-two novel sequences were established and permitted a first phylogenetic analysis of this region of the SRLV genome. The variability of the MA encoding region is higher compared to the gag region encoding the capsid protein and surprisingly close to that reported for the env gene. In contrast to primate lentiviruses, the deduced amino acid sequences of the N- and C-terminal domains of MA are variable. This permitted to pinpoint a basic domain in the N-terminal domain that is conserved in all lentiviruses and likely to play an important functional role. Additionally, a seven amino acid insertion was detected in all MVV strains, which may be used to differentiate CAEV and MVV isolates. A molecular epidemiology analysis based on these sequences indicates that the Italian lentivirus strains are closely related to each other and to the CAEV-CO strain, a prototypic strain isolated three decades ago in the US. This suggests a common origin of the SRLV circulating in the monitored flocks, possibly related to the introduction of infected goats in a negative population. Finally, this study shows that the MA region is suitable for phylogenetic studies and may be applied to monitor SRLV eradication programs.

Phylogenetic analysis of "Candidatus Mycoplasma turicensis" isolates from pet cats in the United Kingdom, Australia, and South Africa, with analysis of risk factors for infection

Two hemotropic mycoplasmas have been recognized in cats, Mycoplasma haemofelis and "Candidatus Mycoplasma haemominutum." We recently described a third feline hemoplasma species, designated "Candidatus Mycoplasma turicensis," in a Swiss cat with hemolytic anemia. This isolate induced anemia after experimental transmission to two specific-pathogen-free cats and analysis of the 16S rRNA gene revealed its close relationship to rodent hemotropic mycoplasmas. The agent was recently shown to be prevalent in Swiss pet cats. We sought to investigate the presence and clinical importance of "Candidatus Mycoplasma turicensis" infection in pet cats outside of Switzerland and to perform the molecular characterization of isolates from different countries. A "Candidatus Mycoplasma turicensis"-specific real-time PCR assay was applied to blood samples from 426 United Kingdom (UK), 147 Australian, and 69 South African pet cats. The 16S rRNA genes of isolates from different countries were sequenced and signalment and laboratory data for the cats were evaluated for associations with "Candidatus Mycoplasma turicensis" infection. Infections were detected in samples from UK, Australian, and South African pet cats. Infection was associated with the male gender, and "Candidatus Mycoplasma haemominutum" and M. haemofelis coinfection. Coinfected cats exhibited significantly lower packed cell volume (PCV) values than uninfected cats. Phylogenetic analyses revealed that some Australian and South African "Candidatus Mycoplasma turicensis" isolates branched away from the remaining isolates. In summary, "Candidatus Mycoplasma turicensis" infection in pet cats exists over a wide geographical area and significantly decreased PCV values are observed in cats coinfected with other feline hemoplasmas.

Combining spatial-temporal and phylogenetic analysis approaches for improved understanding on global H5N1 transmission

Background Since late 2003, the highly pathogenic influenza A H5N1 had initiated several outbreak waves that swept across the Eurasia and Africa continents. Getting prepared for reassortment or mutation of H5N1 viruses has become a global priority. Although the spreading mechanism of H5N1 has been studied from different perspectives, its main transmission agents and spread route problems remain unsolved. Methodology/Principal Findings Based on a compilation of the time and location of global H5N1 outbreaks from November 2003 to December 2006, we report an interdisciplinary effort that combines the geospatial informatics approach with a bioinformatics approach to form an improved understanding on the transmission mechanisms of H5N1 virus. Through a spherical coordinate based analysis, which is not conventionally done in geographical analyses, we reveal obvious spatial and temporal clusters of global H5N1 cases on different scales, which we consider to be associated with two different transmission modes of H5N1 viruses. Then through an interdisciplinary study of both geographic and phylogenetic analysis, we obtain a H5N1 spreading route map. Our results provide insight on competing hypotheses as to which avian hosts are responsible for the spread of H5N1. Conclusions/Significance We found that although South China and Southeast Asia may be the virus pool of avian flu, East Siberia may be the source of the H5N1 epidemic. The concentration of migratory birds from different places increases the possibility of gene mutation. Special attention should be paid to East Siberia, Middle Siberia and South China for improved surveillance of H5N1 viruses and monitoring of migratory birds.

Genome-wide analysis of rare copy number variations reveals PARK2 as a candidate gene for attention-deficit/hyperactivity disorder

Attention-deficit/hyperactivity disorder (ADHD) is a common, highly heritable neurodevelopmental disorder. Genetic loci have not yet been identified by genome-wide association studies. Rare copy number variations (CNVs), such as chromosomal deletions or duplications, have been implicated in ADHD and other neurodevelopmental disorders. To identify rare (frequency 1%) CNVs that increase the risk of ADHD, we performed a whole-genome CNV analysis based on 489 young ADHD patients and 1285 adult population-based controls and identified one significantly associated CNV region. In tests for a global burden of large (>500 kb) rare CNVs, we observed a nonsignificant (P=0.271) 1.126-fold enriched rate of subjects carrying at least one such CNV in the group of ADHD cases. Locus-specific tests of association were used to assess if there were more rare CNVs in cases compared with controls. Detected CNVs, which were significantly enriched in the ADHD group, were validated by quantitative (q)PCR. Findings were replicated in an independent sample of 386 young patients with ADHD and 781 young population-based healthy controls. We identified rare CNVs within the parkinson protein 2 gene (PARK2) with a significantly higher prevalence in ADHD patients than in controls (P=2.8 × 10(-4) after empirical correction for genome-wide testing). In total, the PARK2 locus (chr 6: 162 659 756-162 767 019) harboured three deletions and nine duplications in the ADHD patients and two deletions and two duplications in the controls. By qPCR analysis, we validated 11 of the 12 CNVs in ADHD patients (P=1.2 × 10(-3) after empirical correction for genome-wide testing). In the replication sample, CNVs at the PARK2 locus were found in four additional ADHD patients and one additional control (P=4.3 × 10(-2)). Our results suggest that copy number variants at the PARK2 locus contribute to the genetic susceptibility of ADHD. Mutations and CNVs in PARK2 are known to be associated with Parkinson disease.Molecular Psychiatry advance online publication, 20 November 2012; doi:10.1038/mp.2012.161.

New method for quantification and statistical analysis of nociceptive reflex receptive fields in humans

A method for quantifying nociceptive withdrawal reflex receptive fields in human volunteers and patients is described. The reflex receptive field (RRF) for a specific muscle denotes the cutaneous area from which a muscle contraction can be evoked by a nociceptive stimulus. The method is based on random stimulations presented in a blinded sequence to 10 stimulation sites. The sensitivity map is derived by interpolating the reflex responses evoked from the 10 sites. A set of features describing the size and location of the RRF is presented based on statistical analysis of the sensitivity map within every subject. The features include RRF area, volume, peak location and center of gravity. The method was applied to 30 healthy volunteers. Electrical stimuli were applied to the sole of the foot evoking reflexes in the ankle flexor tibialis anterior. The RRF area covered a fraction of 0.57+/-0.06 (S.E.M.) of the foot and was located on the medial, distal part of the sole of the foot. An intramuscular injection into flexor digitorum brevis of capsaicin was performed in one spinal cord injured subject to attempt modulation of the reflex receptive field. The RRF area, RRF volume and location of the peak reflex response appear to be the most sensitive measures for detecting modulation of spinal nociceptive processing. This new method has important potential applications for exploring aspects of central plasticity in volunteers and patients. It may be utilized as a new diagnostic tool for central hypersensitivity and quantification of therapeutic interventions.

Phylogenetic analysis of Prevotella nigrescens, Prevotella intermedia and Porphyromonas gingivalis clinical strains reveals a clear species clustering

Prevotella nigrescens, Prevotella intermedia and Porphyromonas gingivalis are oral pathogens from the family Bacteroidaceae, regularly isolated from cases of gingivitis and periodontitis. In this study, the phylogenetic variability of these three bacterial species was investigated by means of 16S rRNA (rrs) gene sequence comparisons of a set of epidemiologically and geographically diverse isolates. For each of the three species, the rrs gene sequences of 11 clinical isolates as well as the corresponding type strains was determined. Comparison of all rrs sequences obtained with those of closely related species revealed a clear clustering of species, with only a little intraspecies variability but a clear difference in the rrs gene with respect to the next related taxon. The results indicate that the three species form stable, homogeneous genetic groups, which favours an rrs-based species identification of these oral pathogens. This is especially useful given the 7% sequence divergence between Prevotella intermedia and Prevotella nigrescens, since phenotypic distinction between the two Prevotella species is inconsistent or involves techniques not applicable in routine identification.

Phylogenetic analysis of Pasteurella multocida subspecies and molecular identification of feline P. multocida subsp. septica by 16S rRNA gene sequencing

Pasteurella multocida is commonly found in the oral cavity of cats and dogs. In humans it is known as an opportunistic pathogen after bites from these animals. Phenotypic identification of P. multocida based on biochemical reactions is often limited and usually only done on a species level, even though 3 subspecies are described. For molecular taxonomy and diagnostic purposes a phylogenetic analysis of the three subspecies of P. multocida based on their 16S rRNA (rrs) gene sequence was therefore carried out. We found P. multocida subsp. septica on a distinguished branch on the phylogenetic tree of Pasteurellaceae, due to a 1.5% divergence of its rrs gene compared to the two other, more closely related subspecies multocida and gallicida. This phylogenetic divergence can be used for the identification of P. multocida subsp. septica by rrs gene determination since they form a phylogenetically well isolated and defined group as shown with a set of feline isolates. Comparison to routine phenotypic identification shows the advantage of the sequence-based identification over conventional methods. It is therefore helpful for future unambiguous identification and molecular taxonomy of P. multocida as well as for epidemiological investigations.

Genome-wide analysis in German shepherd dogs reveals association of a locus on CFA 27 with atopic dermatitis

Humans and dogs are both affected by the allergic skin disease atopic dermatitis (AD), caused by an interaction between genetic and environmental factors. The German shepherd dog (GSD) is a high-risk breed for canine AD (CAD). In this study, we used a Swedish cohort of GSDs as a model for human AD. Serum IgA levels are known to be lower in GSDs compared to other breeds. We detected significantly lower IgA levels in the CAD cases compared to controls (p = 1.1 × 10(-5)) in our study population. We also detected a separation within the GSD cohort, where dogs could be grouped into two different subpopulations. Disease prevalence differed significantly between the subpopulations contributing to population stratification (λ = 1.3), which was successfully corrected for using a mixed model approach. A genome-wide association analysis of CAD was performed (n cases = 91, n controls = 88). IgA levels were included in the model, due to the high correlation between CAD and low IgA levels. In addition, we detected a correlation between IgA levels and the age at the time of sampling (corr = 0.42, p = 3.0 × 10(-9)), thus age was included in the model. A genome-wide significant association was detected on chromosome 27 (praw = 3.1 × 10(-7), pgenome = 0.03). The total associated region was defined as a ~1.5-Mb-long haplotype including eight genes. Through targeted re-sequencing and additional genotyping of a subset of identified SNPs, we defined 11 smaller haplotype blocks within the associated region. Two blocks showed the strongest association to CAD. The ~209-kb region, defined by the two blocks, harbors only the PKP2 gene, encoding Plakophilin 2 expressed in the desmosomes and important for skin structure. Our results may yield further insight into the genetics behind both canine and human AD.

Genome-wide analysis reveals selection for important traits in domestic horse breeds.

Intense selective pressures applied over short evolutionary time have resulted in homogeneity within, but substantial variation among, horse breeds. Utilizing this population structure, 744 individuals from 33 breeds, and a 54,000 SNP genotyping array, breed-specific targets of selection were identified using an F(ST)-based statistic calculated in 500-kb windows across the genome. A 5.5-Mb region of ECA18, in which the myostatin (MSTN) gene was centered, contained the highest signature of selection in both the Paint and Quarter Horse. Gene sequencing and histological analysis of gluteal muscle biopsies showed a promoter variant and intronic SNP of MSTN were each significantly associated with higher Type 2B and lower Type 1 muscle fiber proportions in the Quarter Horse, demonstrating a functional consequence of selection at this locus. Signatures of selection on ECA23 in all gaited breeds in the sample led to the identification of a shared, 186-kb haplotype including two doublesex related mab transcription factor genes (DMRT2 and 3). The recent identification of a DMRT3 mutation within this haplotype, which appears necessary for the ability to perform alternative gaits, provides further evidence for selection at this locus. Finally, putative loci for the determination of size were identified in the draft breeds and the Miniature horse on ECA11, as well as when signatures of selection surrounding candidate genes at other loci were examined. This work provides further evidence of the importance of MSTN in racing breeds, provides strong evidence for selection upon gait and size, and illustrates the potential for population-based techniques to find genomic regions driving important phenotypes in the modern horse.

Genome-wide analysis of genetic and epigenetic control of programmed DNA deletion

During the development of the somatic genome from the Paramecium germline genome the bulk of the copies of ∼45 000 unique, internal eliminated sequences (IESs) are deleted. IES targeting is facilitated by two small RNA (sRNA) classes: scnRNAs, which relay epigenetic information from the parental nucleus to the developing nucleus, and iesRNAs, which are produced and used in the developing nucleus. Why only certain IESs require sRNAs for their removal has been enigmatic. By analyzing the silencing effects of three genes: PGM (responsible for DNA excision), DCL2/3 (scnRNA production) and DCL5 (iesRNA production), we identify key properties required for IES elimination. Based on these results, we propose that, depending on the exact combination of their lengths and end bases, some IESs are less efficiently recognized or excised and have a greater requirement for targeting by scnRNAs and iesRNAs. We suggest that the variation in IES retention following silencing of DCL2/3 is not primarily due to scnRNA density, which is comparatively uniform relative to IES retention, but rather the genetic properties of IESs. Taken together, our analyses demonstrate that in Paramecium the underlying genetic properties of developmentally deleted DNA sequences are essential in determining the sensitivity of these sequences to epigenetic control.

Towards Faster Method Search Through Static Ecosystem Analysis

Software developers are often unsure of the exact name of the method they need to use to invoke the desired behavior in a given context. This results in a process of searching for the correct method name in documentation, which can be lengthy and distracting to the developer. We can decrease the method search time by enhancing the documentation of a class with the most frequently used methods. Usage frequency data for methods is gathered by analyzing other projects from the same ecosystem - written in the same language and sharing dependencies. We implemented a proof of concept of the approach for Pharo Smalltalk and Java. In Pharo Smalltalk, methods are commonly searched for using a code browser tool called "Nautilus", and in Java using a web browser displaying HTML based documentation - Javadoc. We developed plugins for both browsers and gathered method usage data from open source projects, in order to increase developer productivity by reducing method search time. A small initial evaluation has been conducted showing promising results in improving developer productivity.