16S rRNA terminal restriction fragment length polymorphism for the characterization of the nasopharyngeal microbiota

A novel non-culture based 16S rRNA Terminal Restriction Fragment Length Polymorphism (T-RFLP) method using the restriction enzymes Tsp509I and Hpy166II was developed for the characterization of the nasopharyngeal microbiota and validated using recently published 454 pyrosequencing data. 16S rRNA gene T-RFLP for 153 clinical nasopharyngeal samples from infants with acute otitis media (AOM) revealed 5 Tsp509I and 6 Hpy166II terminal fragments (TFs) with a prevalence of >10%. Cloning and sequencing identified all TFs with a prevalence >6% allowing a sufficient description of bacterial community changes for the most important bacterial taxa. The conjugated 7-valent pneumococcal polysaccharide vaccine (PCV-7) and prior antibiotic exposure had significant effects on the bacterial composition in an additive main effects and multiplicative interaction model (AMMI) in concordance with the 16S rRNA 454 pyrosequencing data. In addition, the presented T-RFLP method is able to discriminate S. pneumoniae from other members of the Mitis group of streptococci, which therefore allows the identification of one of the most important human respiratory tract pathogens. This is usually not achieved by current high throughput sequencing protocols. In conclusion, the presented 16S rRNA gene T-RFLP method is a highly robust, easy to handle and a cheap alternative to the computationally demanding next-generation sequencing analysis. In case a lot of nasopharyngeal samples have to be characterized, it is suggested to first perform 16S rRNA T-RFLP and only use next generation sequencing if the T-RFLP nasopharyngeal patterns differ or show unknown TFs.

Genotyping of human and porcine Yersinia enterocolitica, Yersinia intermedia, and Yersinia bercovieri strains from Switzerland by amplified fragment length polymorphism analysis

In this study, 231 strains of Yersinia enterocolitica, 25 strains of Y. intermedia, and 10 strains of Y. bercovieri from human and porcine sources (including reference strains) were analyzed using amplified fragment length polymorphism (AFLP), a whole-genome fingerprinting method for subtyping bacterial isolates. AFLP typing distinguished the different Yersinia species examined. Representatives of Y. enterocolitica biotypes 1A, 1B, 2, 3, and 4 belonged to biotype-related AFLP clusters and were clearly distinguished from each other. Y. enterocolitica biotypes 2, 3, and 4 appeared to be more closely related to each other (83% similarity) than to biotypes 1A (11%) and 1B (47%). Biotype 1A strains exhibited the greatest genetic heterogeneity of the biotypes studied. The biotype 1A genotypes were distributed among four major clusters, each containing strains from both human and porcine sources, confirming the zoonotic potential of this organism. The AFLP technique is a valuable genotypic method for identification and typing of Y. enterocolitica and other Yersinia spp.

Genotyping of Francisella tularensis strains by pulsed-field gel electrophoresis, amplified fragment length polymorphism fingerprinting, and 16S rRNA gene sequencing

We evaluated three molecular methods for identification of Francisella strains: pulsed-field gel electrophoresis (PFGE), amplified fragment length polymorphism (AFLP) analysis, and 16S rRNA gene sequencing. The analysis was performed with 54 Francisella tularensis subsp. holarctica, 5 F. tularensis subsp. tularensis, 2 F. tularensis subsp. novicida, and 1 F. philomiragia strains. On the basis of the combination of results obtained by PFGE with the restriction enzymes XhoI and BamHI, PFGE revealed seven pulsotypes, which allowed us to discriminate the strains to the subspecies level and which even allowed us to discriminate among some isolates of F. tularensis subsp. holarctica. The AFLP analysis technique produced some degree of discrimination among F. tularensis subsp. holarctica strains (one primary cluster with three major subclusters and minor variations within subclusters) when EcoRI-C and MseI-A, EcoRI-T and MseI-T, EcoRI-A and MseI-C, and EcoRI-0 and MseI-CA were used as primers. The degree of similarity among the strains was about 94%. The percent similarities of the AFLP profiles of this subspecies compared to those of F. tularensis subsp. tularensis, F. tularensis subsp. novicida, and F. philomiragia were less than 90%, about 72%, and less than 24%, respectively, thus permitting easy differentiation of this subspecies. 16S rRNA gene sequencing revealed 100% similarity for all F. tularensis subsp. holarctica isolates compared in this study. These results suggest that although limited genetic heterogeneity among F. tularensis subsp. holarctica isolates was observed, PFGE and AFLP analysis appear to be promising tools for the diagnosis of infections caused by different subspecies of F. tularensis and suitable techniques for the differentiation of individual strains.

A second generation radiation hybrid map to aid the assembly of the bovine genome sequence

BACKGROUND: Several approaches can be used to determine the order of loci on chromosomes and hence develop maps of the genome. However, all mapping approaches are prone to errors either arising from technical deficiencies or lack of statistical support to distinguish between alternative orders of loci. The accuracy of the genome maps could be improved, in principle, if information from different sources was combined to produce integrated maps. The publicly available bovine genomic sequence assembly with 6x coverage (Btau_2.0) is based on whole genome shotgun sequence data and limited mapping data however, it is recognised that this assembly is a draft that contains errors. Correcting the sequence assembly requires extensive additional mapping information to improve the reliability of the ordering of sequence scaffolds on chromosomes. The radiation hybrid (RH) map described here has been contributed to the international sequencing project to aid this process. RESULTS: An RH map for the 30 bovine chromosomes is presented. The map was built using the Roslin 3000-rad RH panel (BovGen RH map) and contains 3966 markers including 2473 new loci in addition to 262 amplified fragment-length polymorphisms (AFLP) and 1231 markers previously published with the first generation RH map. Sequences of the mapped loci were aligned with published bovine genome maps to identify inconsistencies. In addition to differences in the order of loci, several cases were observed where the chromosomal assignment of loci differed between maps. All the chromosome maps were aligned with the current 6x bovine assembly (Btau_2.0) and 2898 loci were unambiguously located in the bovine sequence. The order of loci on the RH map for BTA 5, 7, 16, 22, 25 and 29 differed substantially from the assembled bovine sequence. From the 2898 loci unambiguously identified in the bovine sequence assembly, 131 mapped to different chromosomes in the BovGen RH map. CONCLUSION: Alignment of the BovGen RH map with other published RH and genetic maps showed higher consistency in marker order and chromosome assignment than with the current 6x sequence assembly. This suggests that the bovine sequence assembly could be significantly improved by incorporating additional independent mapping information.

Differential diagnosis of the honey bee trypanosomatids Crithidia mellificae and Lotmaria passim

Trypanosomatids infecting honey bees have been poorly studied with molecular methods until recently. After the description of Crithidia mellificae (Langridge and McGhee, 1967) it took about forty years until molecular data for honey bee trypanosomatids became available and were used to identify and describe a new trypanosomatid species from honey bees, Lotmaria passim (Evans and Schwarz, 2014). However, an easy method to distinguish them without sequencing is not yet available. Research on the related bumble bee parasites Crithidia bombi and Crithidia expoeki revealed a fragment length polymorphism in the internal transcribed spacer 1 (ITS1), which enabled species discrimination. In search of fragment length polymorphisms for differential diagnostics in honey bee trypanosomatids, we studied honey bee trypanosomatid cell cultures of C. mellificae and L. passim. This research resulted in the identification of fragment length polymorphisms in ITS1 and ITS1-2 markers, which enabled us to develop a diagnostic method to differentiate both honey bee trypanosomatid species without the need for sequencing. However, the amplification success of the ITS1 marker depends probably on the trypanosomatid infection level. Further investigation confirmed that L. passim is the dominant species in Belgium, Japan and Switzerland. We found C. mellificae only rarely in Belgian honey bee samples, but not in honey bee samples from other countries. C. mellificae was also detected in mason bees (Osmia bicornis and Osmia cornuta) besides in honey bees. Further, the characterization and comparison of additional markers from L. passim strain SF (published as C. mellificae strain SF) and a Belgian honey bee sample revealed very low divergence in the 18S rRNA, ITS1-2, 28S rRNA and cytochrome b sequences. Nevertheless, a variable stretch was observed in the gp63 virulence factor.

Role of the NFKB1 -94ins/delATTG promoter polymorphism in IBD and potential interactions with polymorphisms in the CARD15/NOD2, IKBL, and IL-1RN genes

BACKGROUND: Recently, an association of the NFKB1 polymorphism -94ins/delATTG with ulcerative colitis (UC) has been reported. This 4-bp insertion/deletion polymorphism is localized in the promoter region of the NFKB1 gene and appears to be functionally relevant. The aim of the present study was to confirm the association of the -94ins/delATTG (W/D) NFKB1 promoter polymorphism with UC in a population of German origin and to test for a potential association with Crohn's disease (CD). Furthermore, potential interactions of the -94ins/delATTG polymorphism with the IKBL and the IL-1RN genes should be determined. MATERIALS AND METHODS: The study population comprised 630 patients with CD, 365 patients with UC, and 974 healthy controls. Genotyping was performed using polymerase chain reaction and restriction fragment length polymorphism analysis. For statistical evaluation, the chi-square test and the Fisher exact test were used. RESULTS: No significant association of the W/D NFKB1 polymorphism with CD or UC was detected. In addition, no significant interactions between the -94ins/delATTG NFKB1 polymorphism and polymorphisms within the IKBL and the IL-1RN genes, respectively, were found in CD or UC. Also, no significant interactions of the NFKB1 polymorphism with mutations of the CARD15/NOD2 gene and with clinical phenotypes were detected in CD. Moreover, no associations of the NFKB1 polymorphism were found in UC depending on disease localization. CONCLUSIONS: The present study could not confirm the reported association of the -94ins/delATTG NFKB1 polymorphism with UC and also found no evidence for a role of this polymorphism in CD. The results do not give evidence for a role of this NFKB1 polymorphism in the pathogenesis of UC and CD.

Genetic polymorphisms of manganese-superoxide dismutase and glutathione-S-transferase in chronic alcoholic pancreatitis

Chronic alcohol consumption is a major risk factor for the development of chronic pancreatitis. However, chronic pancreatitis occurs only in a minority of heavy drinkers. This variability may be due to yet unidentified genetic factors. Several enzymes involved in the degradation of reactive oxidants and xenobiotics, such as glutathione-S-transferase P1 (GSTP1) and manganese-superoxide dismutase (MnSOD) reveal functional polymorphisms that affect the antioxidative capacity and may therefore modulate the development of chronic pancreatitis and long-term complications like endocrine and exocrine pancreatic insufficiency. Two functional polymorphisms of the MnSOD and the GSTP1 gene were assessed by polymerase chain reaction and restriction fragment length polymorphism in 165 patients with chronic alcoholic pancreatitis, 140 alcoholics without evidence of pancreatic disease and 160 healthy control subjects. The distribution of GSTP1 and MnSOD genotypes were in Hardy-Weinberg equilibrium in the total cohort. Genotype and allele frequencies for both genes were not statistically different between the three groups. Although genotype MnSOD Ala/Val was seemingly associated with the presence of exocrine pancreatic insufficiency, this subgroup was too small and the association statistically underpowered. None of the tested genotypes affected the development of endocrine pancreatic insufficiency. Polymorphisms of MnSOD and GSTP1 are not associated with chronic alcoholic pancreatitis. The present data emphasize the need for stringently designed candidate gene association studies with well-characterized cases and controls and sufficient statistical power to exclude chance observations.

Genetic polymorphisms associated with the inflammatory response in bacterial meningitis.

BACKGROUND Bacterial meningitis (BM) is an infectious disease that results in high mortality and morbidity. Despite efficacious antibiotic therapy, neurological sequelae are often observed in patients after disease. Currently, the main challenge in BM treatment is to develop adjuvant therapies that reduce the occurrence of sequelae. In recent papers published by our group, we described the associations between the single nucleotide polymorphisms (SNPs) AADAT +401C > T, APEX1 Asn148Glu, OGG1 Ser326Cys and PARP1 Val762Ala and BM. In this study, we analyzed the associations between the SNPs TNF -308G > A, TNF -857C > T, IL-8 -251A > T and BM and investigated gene-gene interactions, including the SNPs that we published previously. METHODS The study was conducted with 54 BM patients and 110 healthy volunteers (as the control group). The genotypes were investigated via primer-introduced restriction analysis-polymerase chain reaction (PIRA-PCR) or polymerase chain reaction-based restriction fragment length polymorphism (PCR-RFLP) analysis. Allelic and genotypic frequencies were also associated with cytokine and chemokine levels, as measured with the x-MAP method, and cell counts. We analyzed gene-gene interactions among SNPs using the generalized multifactor dimensionality reduction (GMDR) method. RESULTS We did not find significant association between the SNPs TNF -857C > T and IL-8 -251A > T and the disease. However, a higher frequency of the variant allele TNF -308A was observed in the control group, associated with changes in cytokine levels compared to individuals with wild type genotypes, suggesting a possible protective role. In addition, combined inter-gene interaction analysis indicated a significant association between certain genotypes and BM, mainly involving the alleles APEX1 148Glu, IL8 -251 T and AADAT +401 T. These genotypic combinations were shown to affect cyto/chemokine levels and cell counts in CSF samples from BM patients. CONCLUSIONS In conclusion, this study revealed a significant association between genetic variability and altered inflammatory responses, involving important pathways that are activated during BM. This knowledge may be useful for a better understanding of BM pathogenesis and the development of new therapeutic approaches.

Population genomic tests of models of adaptive radiation in Lake Victoria Region cichlid fish

Adaptive radiation is usually thought to be associated with speciation, but the evolution of intraspecific polymorphisms without speciation is also possible. The radiation of cichlid fish in Lake Victoria (LV) is perhaps the most impressive example of a recent rapid adaptive radiation, with 600+ very young species. Key questions about its origin remain poorly characterized, such as the importance of speciation versus polymorphism, whether species persist on evolutionary time scales, and if speciation happens more commonly in small isolated or in large connected populations. We used 320 individuals from 105 putative species from Lakes Victoria, Edward, Kivu, Albert, Nabugabo and Saka, in a radiation-wide amplified fragment length polymorphism (AFLP) genome scan to address some of these questions. We demonstrate pervasive signatures of speciation supporting the classical model of adaptive radiation associated with speciation. A positive relationship between the age of lakes and the average genomic differentiation of their species, and a significant fraction of molecular variance explained by above-species level taxonomy suggest the persistence of species on evolutionary time scales, with radiation through sequential speciation rather than a single starburst. Finally the large gene diversity retained from colonization to individual species in every radiation suggests large effective population sizes and makes speciation in small geographical isolates unlikely.

CYP2E1 genotype and isoniazid-induced hepatotoxicity in patients treated for latent tuberculosis

OBJECTIVE: To determine whether pharmacogenetic tests such as N-acetyltransferase 2 (NAT2) and cytochrome P450 2E1 (CYP2E1) genotyping are useful in identifying patients prone to antituberculosis drug-induced hepatotoxicity in a cosmopolite population. METHODS: In a prospective study we genotyped 89 patients treated with isoniazid (INH) for latent tuberculosis. INH-induced hepatitis (INH-H) or elevated liver enzymes including hepatitis (INH-ELE) was diagnosed based on the clinical diagnostic scale (CDS) designed for routine clinical practice. NAT2 genotypes were assessed by fluorescence resonance energy transfer probe after PCR analysis, and CYP2E1 genotypes were determined by PCR with restriction fragment length polymorphism analysis. RESULTS: Twenty-six patients (29%) had INH-ELE, while eight (9%) presented with INH-H leading to INH treatment interruption. We report no significant influence of NAT2 polymorphism, but we did find a significant association between the CYP2E1 *1A/*1A genotype and INH-ELE (OR: 3.4; 95% CI:1.1-12; p = 0.02) and a non significant trend for INH-H (OR: 5.9; 95% CI: 0.69-270; p = 0.13) compared with other CYP2E1 genotypes. This test for predicting INH-ELE had a positive predictive value (PPV) of 39% (95% CI: 26-54%) and a negative predictive value (NPV) of 84% (95% CI: 69-94%). CONCLUSION: The genotyping of CYP2E1 polymorphisms may be a useful predictive tool in the common setting of a highly heterogeneous population for predicting isoniazid-induced hepatic toxicity. Larger prospective randomized trials are needed to confirm these results.

Nuclear markers reveal unexpected genetic variation and a Congolese-Nilotic origin of the Lake Victoria cichlid species flock

Phylogenetic analyses based on mitochondrial (mt) DNA have indicated that the cichlid species flock of the Lake Victoria region is derived from a single ancestral species found in East African rivers, closely related to the ancestor of the Lake Malawi cichlid species flock. The Lake Victoria flock contains ten times less mtDNA variation than the Lake Malawi radiation, consistent with current estimates of the ages of the lakes. We present results of a phylogenetic investigation using nuclear (amplified fragment length polymorphism) markers and a wider coverage of riverine haplochromines. We demonstrate that the Lake Victoria–Edward flock is derived from the morphologically and ecologically diverse cichlid genus Thoracochromis from the Congo and Nile, rather than from the phenotypically conservative East African Astatotilapia. This implies that the ability to express much of the morphological diversity found in the species flock may by far pre–date the origin of the flock. Our data indicate that the nuclear diversity of the Lake Victoria–Edward species flock is similar to that of the Lake Malawi flock, indicating that the genetic diversity is considerably older than the 15 000 years that have passed since the lake began to refill. Most of this variation is manifested in trans–species polymorphisms, indicating very recent cladogenesis from a genetically very diverse founder stock. Our data do not confirm strict monophyly of either of the species flocks, but raise the possibility that these flocks have arisen from hybrid swarms.

Continental-Scale Footprint of Balancing and Positive Selection in a Small Rodent (Microtus arvalis)

Genetic adaptation to different environmental conditions is expected to lead to large differences between populations at selected loci, thus providing a signature of positive selection. Whereas balancing selection can maintain polymorphisms over long evolutionary periods and even geographic scale, thus leads to low levels of divergence between populations at selected loci. However, little is known about the relative importance of these two selective forces in shaping genomic diversity, partly due to difficulties in recognizing balancing selection in species showing low levels of differentiation. Here we address this problem by studying genomic diversity in the European common vole (Microtus arvalis) presenting high levels of differentiation between populations (average FST = 0.31). We studied 3,839 Amplified Fragment Length Polymorphism (AFLP) markers genotyped in 444 individuals from 21 populations distributed across the European continent and hence over different environmental conditions. Our statistical approach to detect markers under selection is based on a Bayesian method specifically developed for AFLP markers, which treats AFLPs as a nearly codominant marker system, and therefore has increased power to detect selection. The high number of screened populations allowed us to detect the signature of balancing selection across a large geographic area. We detected 33 markers potentially under balancing selection, hence strong evidence of stabilizing selection in 21 populations across Europe. However, our analyses identified four-times more markers (138) being under positive selection, and geographical patterns suggest that some of these markers are probably associated with alpine regions, which seem to have environmental conditions that favour adaptation. We conclude that despite favourable conditions in this study for the detection of balancing selection, this evolutionary force seems to play a relatively minor role in shaping the genomic diversity of the common vole, which is more influenced by positive selection and neutral processes like drift and demographic history.

No evidence for large differences in genomic methylation between wild and hatchery steelhead (Oncorhynchus mykiss)

When salmonid fish that have been raised in hatcheries spawn in the wild, they often produce fewer surviving adult offspring than wild fish. Recent data from steelhead (Oncorhynchus mykiss) in the Hood River (Oregon, USA) show that even one or two generations of hatchery culture can result in dramatic declines in fitness. Although intense domestication selection could cause such declines, it is worth considering alternative explanations. One possibility is heritable epigenetic changes induced by the hatchery environment. Here, we show, using methylation-sensitive amplified fragment length polymorphism, that hatchery and wild adult steelhead from the Hood River do not appear to differ substantially in overall levels of genomic methylation. Thus, although altered methylation of specific DNA sites or other epigenetic processes could still be important, the hatchery environment does not appear to cause a global hypo- or hypermethylation of the genome or create a large number of sites that are differentially methylated.

Hybridization between distant lineages increases adaptive variation during a biological invasion: stickleback in Switzerland

The three-spined stickleback is a widespread Holarctic species complex that radiated from the sea into freshwaters after the retreat of the Pleistocene ice sheets. In Switzerland, sticklebacks were absent with the exception of the far northwest, but different introduced populations have expanded to occupy a wide range of habitats since the late 19th century. A well-studied adaptive phenotypic trait in sticklebacks is the number of lateral plates. With few exceptions, freshwater and marine populations in Europe are fixed for either the low plated phenotype or the fully plated phenotype, respectively. Switzerland, in contrast, harbours in close proximity the full range of phenotypic variation known from across the continent. We addressed the phylogeographic origins of Swiss sticklebacks using mitochondrial partial cytochrome b and control region sequences. We found only five different haplotypes but these originated from three distinct European regions, fixed for different plate phenotypes. These lineages occur largely in isolation at opposite ends of Switzerland, but co-occur in a large central part. Across the country, we found a strong correlation between a microsatellite linked to the high plate ectodysplasin allele and the mitochondrial haplotype from a region where the fully plated phenotype is fixed. Phylogenomic and population genomic analysis of 481 polymorphic amplified fragment length polymorphism loci indicate genetic admixture in the central part of the country. The same part of the country also carries elevated within-population phenotypic variation. We conclude that during the recent invasive range expansion of sticklebacks in Switzerland, adaptive and neutral between-population genetic variation was converted into within-population variation, raising the possibility that hybridization between colonizing lineages contributed to the ecological success of sticklebacks in Switzerland.

Prevalence and genotypes of Toxoplasma gondii in feline faeces (oocysts) and meat from sheep, cattle and pigs in Switzerland

The protozoan parasite Toxoplasma gondii infects almost all warm blooded animal species including humans, and is one of the most prevalent zoonotic parasites worldwide. Post-natal infection in humans is acquired through oral uptake of sporulated T. gondii oocysts or by ingestion of parasite tissue cysts upon consumption of raw or undercooked meat. This study was undertaken to determine the prevalence of oocyst-shedding by cats and to assess the level of infection with T. gondii in meat-producing animals in Switzerland via detection of genomic DNA (gDNA) in muscle samples. In total, 252 cats (44 stray cats, 171 pet cats, 37 cats with gastrointestinal disorders) were analysed coproscopically, and subsequently species-specific identification of T. gondii oocysts was achieved by Polymerase Chain Reaction (PCR). Furthermore, diaphragm samples of 270 domestic pigs (120 adults, 50 finishing, and 100 free-range animals), 150 wild boar, 250 sheep (150 adults and 100 lambs) and 406 cattle (47 calves, 129 heifers, 100 bulls, and 130 adult cows) were investigated by T. gondii-specific real-time PCR. For the first time in Switzerland, PCR-positive samples were subsequently genotyped using nine PCR-restriction fragment length polymorphism (PCR-RFLP) loci (SAG2, SAG3, BTUB, GRA6, c22-8, c29-2, L358, PK1 and Apico) for analysis. Only one of the cats shed T. gondii oocysts, corresponding to a T. gondii prevalence of 0.4% (95% CI: 0.0-2.2%). In meat-producing animals, gDNA prevalence was lowest in wild boar (0.7%; 95% CI: 0.0-3.7%), followed by sheep (2.0%; 95% CI: 0.1-4.6%) and pigs (2.2%; 95% CI: 0.8-4.8%). The highest prevalence was found in cattle (4.7%; 95% CI: 2.8-7.2%), mainly due to the high prevalence of 29.8% in young calves. With regard to housing conditions, conventional fattening pigs and free-range pigs surprisingly exhibited the same prevalence (2.0%; 95% CI: 0.2-7.0%). Genotyping of oocysts shed by the cat showed T. gondii with clonal Type II alleles and the Apico I allele. T. gondii with clonal Type II alleles were also predominantly observed in sheep, while T. gondii with mixed or atypical allele combinations were very rare in sheep. In pigs and cattle however, genotyping of T. gondii was often incomplete. These findings suggested that cattle in Switzerland might be infected with Toxoplasma of the clonal Types I or III, atypical T. gondii or more than one clonal Type.