Genomic Epidemiology of Salmonella enterica Circulating in Surface Waters Used in Agriculture and Aquaculture in Central Mexico.

Salmonella enterica can survive in surface waters (SuWa), and the role of nonhost environments in its transmission has acquired increasing relevance. In this study, we conducted comparative genomic analyses of 172 S. enterica isolates collected from SuWa across 3 months in six states of central Mexico during 2019. S. enterica transmission dynamics were assessed using 87 experimental and 112 public isolates from Mexico collected during 2002 through 2019. We also studied genetic relatedness between SuWa isolates and human clinical strains collected in North America during 2005 through 2020. Among experimental isolates, we identified 41 S. enterica serovars and 56 multilocus sequence types (STs). Predominant serovars were Senftenberg (n = 13), Meleagridis, Agona, and Newport (n = 12 each), Give (n = 10), Anatum (n = 8), Adelaide (n = 7), and Infantis, Mbandaka, Ohio, and Typhimurium (n = 6 each). We observed a high genetic diversity in the sample under study, as well as clonal dissemination of strains across distant regions. Some of these strains are epidemiologically important (ST14, ST45, ST118, ST132, ST198, and ST213) and were genotypically close to those involved in clinical cases in North America. Transmission network analysis suggests that SuWa are a relevant source of S. enterica (0.7 source/hub ratio) and contribute to its dissemination as isolates from varied sources and clinical cases have SuWa isolates as common ancestors. Overall, the study shows that SuWa act as reservoirs of various S. enterica serovars of public health significance. Further research is needed to better understand the mechanisms involved in SuWa contamination by S. enterica, as well as to develop interventions to contain its dissemination in food production settings. IMPORTANCE Surface waters are heavily used in food production worldwide. Several human pathogens can survive in these waters for long periods and disseminate to food production environments, contaminating our food supply. One of these pathogens is Salmonella enterica, a leading cause of foodborne infections, hospitalizations, and deaths in many countries. This research demonstrates the role of surface waters as a vehicle for the transmission of Salmonella along food production chains. It also shows that some strains circulating in surface waters are very similar to those implicated in human infections and harbor genes that confer resistance to multiple antibiotics, posing a risk to public health. This study contributes to expand our current knowledge on the ecology and epidemiology of Salmonella in surface waters.

Complete DNA sequence analysis of enterohemorrhagic Escherichia coli plasmid pO157_2 in ß-glucuronidase-positive E. coli O157:H7 reveals a novel evolutionary path.

Strains of enterohemorragic Escherichia coli (EHEC) O157:H7 that are non-sorbitol fermenting (NSF) and ß-glucuronidase negative (GUD(-)) carry a large virulence plasmid, pO157 (>90,000 bp), whereas closely related sorbitol-fermenting (SF) E. coli O157:H(-) strains carry plasmid pSFO157 (>120,000 bp). GUD(+) NSF O157:H7 strains are presumed to be precursors of GUD(-) NSF O157:H7 strains that also carry pO157. In this study, we report the complete sequence of a novel virulence plasmid, pO157-2 (89,762 bp), isolated from GUD(+) NSF O157:H7 strain G5101. PCR analysis confirmed the presence of pO157-2 in six other strains of GUD(+) NSF O157:H7. pO157-2 carries genes associated with virulence (e.g., hemolysin genes) and conjugation (tra and trb genes) but lacks katP and espP present in pO157. Comparative analysis of the three EHEC plasmids shows that pO157-2 is highly related to pO157 and pSFO157 but not ancestral to pO157. These results indicated that GUD(+) NSF O157:H7 strains might not be direct precursors to GUD(-) NSF O157:H7 as previously proposed but rather have evolved independently from a common ancestor.

Draft genome sequences of six Escherichia coli isolates from the stepwise model of emergence of Escherichia coli O157:H7.

Enterohemorrhagic Escherichia coli (EHEC) of serotype O157:H7 has been implicated in food-borne illnesses worldwide. An evolutionary model was proposed in which the highly pathogenic EHEC O157:H7 serotype arose from its ancestor, enteropathogenic E. coli (EPEC) O55:H7 (sorbitol fermenting [SOR(+)] and ß-glucuronidase positive [GUD(+)]), through sequential gain of virulence, phenotypic traits, and serotype change. Here we report six draft genomes of strains belonging to this evolutionary model: two EPEC O55:H7 (SOR(+) GUD(+)) strains, two nonmotile EHEC O157:H(-) strains (SOR(+) GUD(+)) containing plasmid pSFO157, one EHEC O157:H7 (SOR(-) GUD(+)) strain, and one O157:H7 strain containing plasmid pSFO157 (SOR(+) GUD(+)).

Phylogenetics and Mitochondrial DNA.

Phylogenetic analysis can be conducted using a variety of methods, generally classified as distancebased or character-based approaches. Patterns found through phylogenetic analysis of human mitochondrial DNA (mtDNA) sequences have revealed a wealth of information in such disparate fields as the human evolution; the movement of human lineages throughout history (phylogeography); and the susceptibility of certain groups to devastating diseases. Forensic mtDNA analysis has also benefited from the use of these methods. Phylogenetic assessment of forensic mtDNA databases has revealed a consistency with published data at a depth of analysis that is not attainable with basic population genetic methods. The detailed characteristics of specific sites within a mtDNA sequence are best assessed using phylogenetic methods. These studies have identified the most informative sites for individual differentiation, while also providing quality assurance metrics to apply to individual mtDNA profiles or entire databases. Such a level of evaluation and understanding enhances the interpretation of forensic casework.

Sources of incongruence among mammalian mitochondrial sequences: COII, COIII, and ND6 genes are main contributors.

To investigate the origins of incongruence among mammalian mitochondrial protein-coding genes, we compiled a matrix that included 13 protein-coding-genes for 41 mammals from 14 different orders. This matrix was examined for congruence using different partitioning strategies. The incongruence length difference test showed significant incongruence among the 13 gene partitions used simultaneously, and the result was not affected by third codon or transversion weighting. In the pair-wise comparisons, significant incongruence was detected between NADH:ubiquinone oxidoreductase subunit 6 gene (ND6), cytochrome oxidase subunit II (COII), or cytochrome oxidase subunit III (COIII) gene partitioned individually against the rest of the genes. Omission of any of the 14 mammalian orders alone or in combinations from the matrix did not result in a statistically significant improvement of congruence, suggesting that taxonomic sampling will not improve congruence among the data sets. However, omission of the ND6, COII, and COIII significantly improved congruence in our data matrix. Possible origins of unusual phylogenetic properties of the three genes are discussed.

A test of Archonta monophyly and the phylogenetic utility of the mitochondrial gene 12S rRNA.

The relationships within the superorder Archonta, which contains the orders Dermoptera (flying lemurs), Scandentia (tree shrews), Chiroptera (bats), and Primates, were examined through the analysis of five newly derived and complete mitochondrial 12S rRNA sequences. The new data is combined with 83 additional known mammalian sequences to provide a full phylogenetic sampling. Phylogenetic hypotheses are generated using PAUP 3.1.1 (Swofford [1993] Illinois Natural History Survey, Champaign, IL) through analyses of all characters equally weighted, transversions only, and the effect of alignment gaps on phylogeny. The Parsimony Jackknifer (Farris et al. [1996] Cladistics 12:99-124) was used to assess the level of ambiguity present in the sequence data, and therefore the strength of the tree topologies. The conclusions of Springer and Douzery (1996, J. Mol. Evol. 43:357-373) which states that 12S rRNA is reliable to a time depth of 100 mya is unsupported by these analyses. The usefulness of 12S rRNA to aid in solving Archonta relationships and others of similar time depth is found to be suspect.

Support for interordinal eutherian relationships with an emphasis on primates and their archontan relatives.

Current knowledge about mammalian interordinal relationships is growing rapidly; thus this contribution is an attempt to summarize the past 5 years of this literature. We have focused on the recent controversies in mammalian phylogenetics including hypotheses concerning the monophyly of Archonta, the diphyly of Chiroptera, and the polyphyly of Rodentia. All of these issues have been proposed recently, challenging these phylogenetic hypotheses. We have attempted to include all of the comprehensive analyses of eutherian mammal systematics with an emphasis on morphological and molecular data sets where discrete characters are listed so they could be compiled and used in support of interordinal relationships. Particular attention is given to determining which of the living eutherian orders is the closest relative to primates. In reviewing relationships among the mammals, we have focused on collating all of the available evidence so that one could know where each of the specific data sets is in support of the various relationships.

Molecular systematics of hystricognath rodents: evidence from the mitochondrial 12S rRNA gene.

Nucleotide sequence variation among 22 representatives of 14 families of hystricognathid rodents was examined using an 814-bp region of the mitochondrial 12S ribosomal RNA (rRNA) gene composing domains I-III. The purpose of this study was twofold. First, the phylogenetic relationships among Old World phiomorph (primarily African) and New World caviomorph (primarily South American) families were investigated, with a special emphasis on testing hypotheses pertaining to the origin of New World families and the identification of major monophyletic groups. Second, divergence times derived from molecular data were compared to those suggested by the fossil record. The resultant 12S rRNA gene phylogeny, analyzed separately and in combination with other morphological and molecular data, supported a monophyletic Caviomorpha. This finding is counter to the idea of a multiple origin for the South American families. The most strongly supported relationships within the Caviomorpha were a monophyletic Octodontoidea (containing five families) and the placement of New World porcupines (family Erethizontidae) as the most divergent family. Although comparisons to other data were more equivocal, the most parsimonious 12S rRNA trees also supported a monophyletic Phiomorpha that could be subdivided into two major groups, a clade containing the Thryonomyoidea (Thryonomyidae and Petromuridae) plus Bathyergidae and the more divergent Hystricidae (Old World porcupines). No significant differences in rates of 12S rRNA gene divergence were observed for hystricognathids in comparison to other rodent groups. Although time since divergence estimates were influenced by the fossil dates chosen to calibrate absolute rates, the overall divergence times derived from both transversions only and Kimura corrected distances and calibrations using two independent dates revealed a divergence time between Old and New World groups dating in the Eocene.

DNA systematics and evolution of the artiodactyl family Bovidae.

Nine additional sequences from representatives of different tribes of the family Bovidae were combined with six published artiodactyl sequences to provide orthologous mtDNA for investigation of bovid phylogeny and evolution. Each species was represented by a homologous 2.7-kilobase-pair stretch of mtDNA for the complete 12S and 16S rRNA genes and three adjacent tRNA genes. These data, when compared to other results, provided evidence for a monophyletic Bovidae and for two clades within the family: one including the tribes Boselaphini, Bovini, and Tragelaphini and another for an Antilopini/Neotragini grouping. All other intrafamilial relationships were only weakly supported. These sequence comparisons suggest that most bovid tribes originated early in the Miocene with all extant lineages present by approximately 16-17 million years ago. Thus, bovid tribes provide an example of rapid cladogenesis, following the origin of families in the infraorder Pecora.

Nucleotide sequence variation in the mitochondrial 12S rRNA gene and the phylogeny of African mole-rats (Rodentia: Bathyergidae).

Mitochondrial DNA (mtDNA) sequence variation was examined in eight taxa of the African rodent family Bathyergidae, as well as in two taxa representative of the Old-World hystricognathid rodent families Petromyidae and Thryonomyidae. A total of 812 bp, constituting domains I-III of the 12S ribosomal rRNA gene, were compared for each taxon. The high levels of intrafamilial mtDNA sequence divergence observed (average 16.8, range 3.5-23.2) support an ancient origin for the five genera, 20-38 Mya. These data do not support the current subfamilial groupings of the Bathyergidae. The eastern African naked mole-rat, Heterocephalus glaber, is the most basal representative of the family, with the silvery mole-rat, Heliophobius, being the next most basal. South African forms [dune, common, and cape mole-rats (Bathyergus, Cryptomys, and Georychus, respectively)] group together. The independent origin of the common mole-rat, relative to the naked mole-rat, suggests that complex social systems evolved in parallel along different bathyergid lineages. The 12S rRNA gene is not evolving at a higher rate within the rodent lineages, relative to that seen for artiodactyls and primates. Bathyergid rodents appear to fall at an extreme end of the spectrum of mammalian variation, with respect to both transition/transversion ratios and divergence, showing much lower transition/transversion ratios than those previously reported for intrafamilial comparisons.

Resolution of the African hominoid trichotomy by use of a mitochondrial gene sequence.

Mitochondrial DNA sequences encoding the cytochrome oxidase subunit II gene have been determined for five primate species, siamang (Hylobates syndactylus), lowland gorilla (Gorilla gorilla), pygmy chimpanzee (Pan paniscus), crab-eating macaque (Macaca fascicularis), and green monkey (Cercopithecus aethiops), and compared with published sequences of other primate and nonprimate species. Comparisons of cytochrome oxidase subunit II gene sequences provide clear-cut evidence from the mitochondrial genome for the separation of the African ape trichotomy into two evolutionary lineages, one leading to gorillas and the other to humans and chimpanzees. Several different tree-building methods support this same phylogenetic tree topology. The comparisons also yield trees in which a substantial length separates the divergence point of gorillas from that of humans and chimpanzees, suggesting that the lineage most immediately ancestral to humans and chimpanzees may have been in existence for a relatively long time.

Ribosomal DNA variation within and between species of rodents, with emphasis on the genus Onychomys.

Patterns of restriction-endonuclease site and length variation at the nuclear rDNA locus (18S + 28S rRNA gene complex) were examined in rodents. Of the 164 restriction sites mapped for seven species, 22 were conserved (mapping to the 18S, 28S, and 5.8S genes and ITS1) in all three Onychomys species as well as in Mus musculus and in three closely related peromyscine rodents, Peromyscus boylii, P. eremicus, and Reithrodontomys megalotis. The nontranscribed spacer (NTS) region revealed most of the variation among these taxa, with the patterns of variation grouping into the following categories, (1) intraindividual variation revealing as many as four site-specific repeat types within an individual, (2) intraspecific and interspecific site variation confined to the NTS, and (3) length variation in both the transcribed and NTS regions. Length variation in the 28S rRNA gene was also examined in 17 additional rodent species, and most size differences mapped to the divergent domain, D8, found in sequence comparisons between Mus and Rattus. The systematic implications of rDNA variation are discussed using the perspective gained from these rodent comparisons.

The production of single-stranded DNA suitable for sequencing using the polymerase chain reaction.

A simple and reliable procedure for the amplification of single-stranded DNA suitable for sequencing is described. This procedure employs the polymerase chain reaction and implements modifications pertaining to the purification of the double-stranded DNA product prior to single-stranded DNA amplification. The most consistent sequencing reactions are obtained when the double-stranded DNA product is purified by centrifugation with a microconcentrator prior to single-stranded DNA amplification and the overall amount of specific primers and number of cycles used, in both single-stranded and double-stranded DNA polymerase chain reactions, are reduced.