Multilocus sequence analysis and rpoB sequencing of Mycobacterium abscessus (sensu lato) strains

Mycobacterium abscessus, Mycobacterium bolletii, and Mycobacterium massiliense (Mycobacterium abscessus sensu lato) are closely related species that currently are identified by the sequencing of the rpoB gene. However, recent studies show that rpoB sequencing alone is insufficient to discriminate between these species, and some authors have questioned their current taxonomic classification. We studied here a large collection of M. abscessus (sensu lato) strains by partial rpoB sequencing (752 bp) and multilocus sequence analysis (MLSA). The final MLSA scheme developed was based on the partial sequences of eight housekeeping genes: argH, cya, glpK, gnd, murC, pgm, pta, and purH. The strains studied included the three type strains (M. abscessus CIP 104536(T), M. massiliense CIP 108297(T), and M. bolletii CIP 108541(T)) and 120 isolates recovered between 1997 and 2007 in France, Germany, Switzerland, and Brazil. The rpoB phylogenetic tree confirmed the existence of three main clusters, each comprising the type strain of one species. However, divergence values between the M. massiliense and M. bolletii clusters all were below 3% and between the M. abscessus and M. massiliense clusters were from 2.66 to 3.59%. The tree produced using the concatenated MLSA gene sequences (4,071 bp) also showed three main clusters, each comprising the type strain of one species. The M. abscessus cluster had a bootstrap value of 100% and was mostly compact. Bootstrap values for the M. massiliense and M. bolletii branches were much lower (71 and 61%, respectively), with the M. massiliense cluster having a fuzzy aspect. Mean (range) divergence values were 2.17% (1.13 to 2.58%) between the M. abscessus and M. massiliense clusters, 2.37% (1.5 to 2.85%) between the M. abscessus and M. bolletii clusters, and 2.28% (0.86 to 2.68%) between the M. massiliense and M. bolletii clusters. Adding the rpoB sequence to the MLSA-concatenated sequence (total sequence, 4,823 bp) had little effect on the clustering of strains. We found 10/120 (8.3%) isolates for which the concatenated MLSA gene sequence and rpoB sequence were discordant (e.g., M. massiliense MLSA sequence and M. abscessus rpoB sequence), suggesting the intergroup lateral transfers of rpoB. In conclusion, our study strongly supports the recent proposal that M. abscessus, M. massiliense, and M. bolletii should constitute a single species. Our findings also indicate that there has been a horizontal transfer of rpoB sequences between these subgroups, precluding the use of rpoB sequencing alone for the accurate identification of the two proposed M. abscessus subspecies.

Comparative human-horse sequence analysis of the CYP3A subfamily gene cluster

Cytochrome P450 enzymes (CYP450s) represent a superfamily of haem-thiolate proteins. CYP450s are most abundant in the liver, a major site of drug metabolism, and play key roles in the metabolism of a variety of substrates, including drugs and environmental contaminants. Interaction of two or more different drugs with the same enzyme can account for adverse effects and failure of therapy. Human CYP3A4 metabolizes about 50% of all known drugs, but little is known about the orthologous CYP450s in horses. We report here the genomic organization of the equine CYP3A gene cluster as well as a comparative analysis with the human CYP3A gene cluster. The equine CYP450 genes of the 3A family are located on ECA 13 between 6.97-7.53 Mb, in a region syntenic to HSA 7 99.05-99.35 Mb. Seven potential, closely linked equine CYP3A genes were found, in contrast to only four genes in the human genome. RNA was isolated from an equine liver sample, and the approximately 1.5-kb coding sequence of six CYP3A genes could be amplified by RT-PCR. Sequencing of the RT-PCR products revealed numerous hitherto unknown single nucleotide polymorphisms (SNPs) in these six CYP3A genes, and one 6-bp deletion compared to the reference sequence (EquCab2.0). The presence of the variants was confirmed in a sample of genomic DNA from the same horse. In conclusion, orthologous genes for the CYP3A family exist in horses, but their number differs from those of the human CYP3A gene family. CYP450 genes of the same family show high homology within and between mammalian species, but can be highly polymorphic.

Sequence analysis of the porcine IFNAR1 and IFNGR2 genes

A porcine BAC clone harboring the tightly linked IFNAR1 and IFNGR2 genes was identified by comparative analysis of the publicly available porcine BAC end sequences. The complete 168,835 bp insert sequence of this clone was determined. Sequence comparisons of the genomic sequence with EST sequences from public databases were performed and allowed a detailed annotation of the IFNAR1 and IFNGR2 genes. The analyzed genes showed a conserved genomic organization with their known mammalian orthologs, however the sequence conservation of these genes across species was relatively low. In addition to the IFNAR1 and IFNGR2 genes, which were completely sequenced, the analyzed BAC clone also contained parts of an orphan gene encoding a putative transmembrane protein (TMEM50B). In contrast to the IFNAR1 and IFNGR2 genes the sequence conservation of the TMEM50B gene across different mammalian species was extremely high.

Sequence analysis of a 212 kb defensin gene cluster on ECA 27q17

Defensins are a family of evolutionary ancient antimicrobial peptides consisting of three sub-families: alpha-, beta- and theta-defensins. This investigation was focused on the genomic characterization of equine beta-defensins and the investigation of the potential clustering of beta-defensin genes in the equine genome. Six genomic BAC clones were isolated from the CHORI-241 library and one of these was mapped by FISH to ECA 27q17. This location was confirmed by RH-mapping. The contiguous 212 kb sequence of this clone was determined. Sequence analysis revealed the identification of ten pseudogenes and nine genes, six of which were highly homologous to human beta-defensin DEFB4. Clustering of the beta-defensin genes was confirmed and the order of the genes on the analyzed BAC was related to the corresponding defensin cluster on HSA 8. The knowledge about the sequence and the genomic structure of the equine beta-defensin genes will improve the classification of different paralogous defensin genes and is a prerequisite for subsequent functional studies. Additionally, the first alpha-defensin-like sequence outside the groups of primates, lagomorphs and rodents (glires) was identified.

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.

Quantitative sequence analysis of FBN1 premature termination codons provides evidence for incomplete NMD in leukocytes

We improved, evaluated, and used Sanger sequencing for quantification of single nucleotide polymorphism (SNP) variants in transcripts and gDNA samples. This improved assay resulted in highly reproducible relative allele frequencies (e.g., for a heterozygous gDNA 50.0+/-1.4%, and for a missense mutation-bearing transcript 46.9+/-3.7%) with a lower detection limit of 3-9%. It provided excellent accuracy and linear correlation between expected and observed relative allele frequencies. This sequencing assay, which can also be used for the quantification of copy number variations (CNVs), methylations, mosaicisms, and DNA pools, enabled us to analyze transcripts of the FBN1 gene in fibroblasts and blood samples of patients with suspected Marfan syndrome not only qualitatively but also quantitatively. We report a total of 18 novel and 19 known FBN1 sequence variants leading to a premature termination codon (PTC), 26 of which we analyzed by quantitative sequencing both at gDNA and cDNA levels. The relative amounts of PTC-containing FBN1 transcripts in fresh and PAXgene-stabilized blood samples were significantly higher (33.0+/-3.9% to 80.0+/-7.2%) than those detected in affected fibroblasts with inhibition of nonsense-mediated mRNA decay (NMD) (11.0+/-2.1% to 25.0+/-1.8%), whereas in fibroblasts without NMD inhibition no mutant alleles could be detected. These results provide evidence for incomplete NMD in leukocytes and have particular importance for RNA-based analyses not only in FBN1 but also in other genes.

Phylogeny and prediction of genetic similarity of Cronobacter and related taxa by multilocus sequence analysis (MLSA)

Multilocus sequence analysis (MLSA) based on recN, rpoA and thdF genes was done on more than 30 species of the family Enterobacteriaceae with a focus on Cronobacter and the related genus Enterobacter. The sequences provide valuable data for phylogenetic, taxonomic and diagnostic purposes. Phylogenetic analysis showed that the genus Cronobacter forms a homogenous cluster related to recently described species of Enterobacter, but distant to other species of this genus. Combining sequence information on all three genes is highly representative for the species' %GC-content used as taxonomic marker. Sequence similarity of the three genes and even of recN alone can be used to extrapolate genetic similarities between species of Enterobacteriaceae. Finally, the rpoA gene sequence, which is the easiest one to determine, provides a powerful diagnostic tool to identify and differentiate species of this family. The comparative analysis gives important insights into the phylogeny and genetic relatedness of the family Enterobacteriaceae and will serve as a basis for further studies and clarifications on the taxonomy of this large and heterogeneous family.

Streaming Support for Data Intensive Cloud-Based Sequence Analysis

Cloud computing provides a promising solution to the genomics data deluge problem resulting from the advent of next-generation sequencing (NGS) technology. Based on the concepts of “resources-on-demand” and “pay-as-you-go”, scientists with no or limited infrastructure can have access to scalable and cost-effective computational resources. However, the large size of NGS data causes a significant data transfer latency from the client’s site to the cloud, which presents a bottleneck for using cloud computing services. In this paper, we provide a streaming-based scheme to overcome this problem, where the NGS data is processed while being transferred to the cloud. Our scheme targets the wide class of NGS data analysis tasks, where the NGS sequences can be processed independently from one another. We also provide the elastream package that supports the use of this scheme with individual analysis programs or with workflow systems. Experiments presented in this paper show that our solution mitigates the effect of data transfer latency and saves both time and cost of computation.

Assessing the Potential of Multiple Imputation in Sequence Analysis

Sequence analysis and optimal matching are useful heuristic tools for the descriptive analysis of heterogeneous individual pathways such as educational careers, job sequences or patterns of family formation. However, to date it remains unclear how to handle the inevitable problems caused by missing values with regard to such analysis. Multiple Imputation (MI) offers a possible solution for this problem but it has not been tested in the context of sequence analysis. Against this background, we contribute to the literature by assessing the potential of MI in the context of sequence analyses using an empirical example. Methodologically, we draw upon the work of Brendan Halpin and extend it to additional types of missing value patterns. Our empirical case is a sequence analysis of panel data with substantial attrition that examines the typical patterns and the persistence of sex segregation in school-to-work transitions in Switzerland. The preliminary results indicate that MI is a valuable methodology for handling missing values due to panel mortality in the context of sequence analysis. MI is especially useful in facilitating a sound interpretation of the resulting sequence types.

Amino acid sequence of the alpha subunit and computer modelling of the alpha and beta subunits of echicetin from the venom of Echis carinatus (saw-scaled viper)

Echicetin, a heterodimeric protein from the venom of Echis carinatus, binds to platelet glycoprotein Ib (GPIb) and so inhibits platelet aggregation or agglutination induced by various platelet agonists acting via GPIb. The amino acid sequence of the beta subunit of echicetin has been reported and found to belong to the recently identified snake venom subclass of the C-type lectin protein family. Echicetin alpha and beta subunits were purified. N-terminal sequence analysis provided direct evidence that the protein purified was echicetin. The paper presents the complete amino acid sequence of the alpha subunit and computer models of the alpha and beta subunits. The sequence of alpha echicetin is highly similar to the alpha and beta chains of various heterodimeric and homodimeric C-type lectins. Neither of the fully reduced and alkylated alpha or beta subunits of echicetin inhibited the platelet agglutination induced by von Willebrand factor-ristocetin or alpha-thrombin. Earlier reports about the inhibitory activity of reduced and alkylated echicetin beta subunit might have been due to partial reduction of the protein.

Molecular genetic analysis of Dichelobacter nodosus proteases AprV2/B2, AprV5/B5 and BprV/B in clinical material from European sheep flocks

Dichelobacter nodosus, the etiological agent of ovine footrot, exists both as virulent and as benign strains, which differ in virulence mainly due to subtle differences in the three subtilisin-like proteases AprV2, AprV5 and BprV found in virulent, and AprB2, AprB5 and BprB in benign strains of D. nodosus. Our objective was a molecular genetic epidemiological analysis of the genes of these proteases by direct sequence analysis from clinical material of sheep from herds with and without history of footrot from 4 different European countries. The data reveal the two proteases known as virulent AprV2 and benign AprB2 to correlate fully to the clinical status of the individuals or the footrot history of the herd. In samples taken from affected herds, the aprV2 gene was found as a single allele whereas in samples from unaffected herds several alleles with minor modifications of the aprB2 gene were detected. The different alleles of aprB2 were related to the herds. The aprV5 and aprB5 genes were found in the form of several alleles scattered without distinction between affected and non-affected herds. However, all different alleles of aprV5 and aprB5 encode the same amino acid sequences, indicating the existence of a single protease isoenzyme 5 in both benign and virulent strains. The genes of the basic proteases BprV and BprB also exist as various alleles. However, differences found in samples from affected versus non-affected herds do not reflect the currently known epitopes that are attributed to differences in biochemical activity. The data of the study confirm the prominent role of AprV2 in the virulence of D. nodosus and shed a new light on the presence of the other protease genes and their allelic variants in clinical samples.

Genetic Relatedness, Antimicrobial and Biocide Susceptibility Comparative Analysis of Methicillin-Resistant and -Susceptible Staphylococcus pseudintermedius from Portugal

Forty methicillin-resistant and -susceptible Staphylococcus pseudintermedius (MRSP and MSSP, respectively) from colonization and infection in dogs and cats were characterized for clonality, antimicrobial, and biocide susceptibility. MSSP were genetically more diverse than MRSP by multi-locus sequence typing and pulsed-field gel electrophoresis. Three different spa types (t06, t02, t05) and two SCCmec types (II-III and V) were detected in the MRSP isolates. All MRSP and two MSSP strains were multidrug-resistant. Several antibiotic resistance genes (mecA, blaZ, tet(M), tet(K), aac(6')-Ie-aph(2')-Ia, aph(3')-III, ant(6)-Ia, sat4, erm(B), lnu(A), dfr(G), and catpC221) were identified by microarray and double mutations in the gyrA and grlA genes and a single mutation in the rpoB gene were detected by sequence analysis. No differences were detected between MSSP and MRSP in the chlorhexidine acetate (CHA) minimum inhibitory concentrations (MICs). However, two MSSP had elevated MIC to triclosan (TCL) and one to benzalkonium chloride and ethidium bromide. One MSSP isolate harboured a qacA gene, while in another a qacB gene was detected. None of the isolates harboured the sh-fabI gene. Three of the biocide products studied had high bactericidal activity (Otodine(®), Clorexyderm Spot Gel(®), Dermocanis Piocure-M(®)), while Skingel(®) failed to achieve a five log reduction in the bacterial counting. S. pseudintermedius have become a serious therapeutic challenge in particular if methicillin- resistance and/or multidrug-resistance are involved. Biocides, like CHA and TCL, seem to be clinically effective and safe topical therapeutic options.

Play in Rats: Association across Contexts and Types , and Analysis of Structure.

Play has been proposed as a promising indicator of positive animal welfare. We aimed to study play in rats across contexts (conspecific/heterospecific) and types (social: pinning, being pinned; solitary: scampering), and we investigated its structure using behavioral sequence analysis. Group-housed (three per cage) adolescent male Lister Hooded rats (n = 21) were subjected to a Play-In-Pairs test: after a 3 hour isolation period, a pair of cage-mates was returned to the home cage and both social and solitary play were scored for 20 min. This procedure was repeated for each pair combination across three consecutive days, and individual play scores were calculated. Heterospecific play was measured using a Tickling test: rats were individually tickled by the experimenter through bouts of gentle, rapid finger movements on their underside, and the number of positive 50 kHz frequency modulated vocalizations and experimenter-directed approach behaviors were recorded. Both of the above tests were compared with social play in the home cage. While conspecific play in both the Play-In-Pairs test and home cage were correlated, both seemed to be unrelated to heterospecific play in the Tickling test. During the Play-In-Pairs test, although both solitary and social play types occurred, they were unrelated, and solitary locomotor play of one rat did not predict the subsequent play behavior of its cage mate. Analysis of play structure revealed that social play occurred more often in bouts of repeated behaviors while solitary play sequences did not follow a specific pattern. If play is to be used as an indicator of positive welfare in rats, context, type and structure differences should be taken into account.

Detection of one VH antibody sequence in both healthy donors and urticaria patients.

We have previously isolated anti-FcepsilonRIalpha autoantibodies from phage libraries of healthy donors and urticaria patients. Strikingly, the same antibody, LTMalpha15, was isolated from both libraries. Sequence analysis revealed a germline configuration of the LTMalpha15 variable heavy (V(H)) chain with a slightly mutated variable light (V(L)) chain supporting its classification as a natural autoantibody. Distribution analysis of anti-FcepsilonRIalpha autoantibodies by functional or serological tests delivered conflicting data. For this reason we have developed a new real-time PCR to analyse the distribution of LTMalpha15V(H) in healthy donors and urticaria patients. Our new bioinformatic program permitted the design of a minor groove binder (MGB) TaqMan probe that specifically detected the LTMalpha15V(H). We were able to demonstrate a broad range of rearranged V(H) gene copy number without any correlation to the state of health. Monitoring LTMalpha15V(H) gene copy number in a single donor over a period of 70 days revealed a time-related fluctuation of circulating B cells carrying LTMalpha15V(H). We propose that our real-time PCR may serve as a model for the quantification of natural antibody sequences at a monoclonal level.