Characterization of Escherichia coli harboring colibactin genes (clb) isolated from beef production and processing systems

Certain strains of Escherichia coli possess and express the toxin colibactin (Clb) which induces host mutations identical to the signature mutations of colorectal cancer (CRC) that lead to tumorigenic lesions. Since cattle are a known reservoir of several Enterobacteriaceae including E. coli, this study screened for clb amongst E. coli isolated from colons of cattle-at-harvest (entering beef processing facility; n = 1430), across the beef processing continuum (feedlot to finished subprimal beef; n = 232), and in ground beef (n = 1074). Results demonstrated that clb⁺E. coli were present in cattle and beef. Prevalence of clb⁺E. coli from colonic contents of cattle and ground beef was 18.3% and 5.5%, respectively. clb⁺E. coli were found susceptible to commonly used meat processing interventions. Whole genome sequencing of 54 bovine and beef clb⁺ isolates showed clb occurred in diverse genetic backgrounds, most frequently in phylogroup B1 (70.4%), MLST 1079 (42.6%), and serogroup O49 (40.7%).

Twenty-Four-Month Longitudinal Study Suggests Little to No Horizontal Gene Transfer In Situ between Third-Generation Cephalosporin-Resistant Salmonella and Third-Generation Cephalosporin-Resistant Escherichia coli in a Beef Cattle Feedyard

ABSTRACT

Third-generation cephalosporins (3GCs) are preferred treatments for serious human Salmonella enterica infections. Beef cattle are suspected to contribute to human 3GC-resistant Salmonella infections. Commensal 3GC-resistant Escherichia coli are thought to act as reservoirs of 3GC resistance because these strains are isolated more frequently than are 3GC-resistant Salmonella strains at beef cattle feedyards. During each of 24 consecutive months, four samples of pen surface material were obtained from five pens (N = 480) at a Nebraska feedyard to determine to the contribution of 3GC-resistant E. coli to the occurrence of 3GC-resistant Salmonella. Illumina whole genome sequencing was performed, and susceptibility to 14 antimicrobial agents was determined for 121 3GC-susceptible Salmonella, 121 3GC-resistant Salmonella, and 203 3GC-resistant E. coli isolates. 3GC-susceptible Salmonella isolates were predominantly from serotypes Muenchen (70.2%) and Montevideo clade 1 (23.1%). 3GC-resistant Salmonella isolates were predominantly from serotypes Montevideo clade 2 (84.3%). One bla gene type (blaCMY-2) and the IncC plasmid replicon were present in 100 and 97.5% of the 3GC-resistant Salmonella, respectively. Eleven bla gene types were detected in the 3GC-resistant E. coli, which were distributed across 42 multilocus sequence types. The blaCMY-2 gene and IncC plasmid replicon were present in 37.9 and 9.9% of the 3GC-resistant E. coli, respectively. These results suggest that 3GC resistance in Salmonella was primarily due the persistence of Salmonella Montevideo clade 2 with very minimal or no contribution from 3GC-resistant E. coli via horizontal gene transfer and that 3GC-resistant E. coli may not be a useful indicator for 3GC-resistant Salmonella in beef cattle production environments.

HIGHLIGHTS

Classification of small ruminant lentivirus subtype A2, subgroups 1 and 2 based on whole genome comparisons and complex recombination patterns

Background: Small ruminant lentiviruses (SRLVs) cause a multisystemic chronic wasting disease in sheep across much of the world. SRLV subtype A2 is prevalent in North America and further classified into multiple subgroups based on variation in the group antigens gene (gag) and envelope (env) genes. In sheep, the ovine transmembrane protein 154 (TMEM154) gene is associated with SRLV susceptibility. Ewes with at least one copy of TMEM154 encoding a full-length protein with glutamate at position 35 (E35; haplotypes 2 and 3), are highly susceptible to SRLV infection while ewes with any combination of TMEM154 haplotypes which encodes lysine (K35; haplotype 1), or truncated proteins (haplotypes 4 and 6) are several times less so. A2 subgroups 1 and 2 are associated with host TMEM154 genotypes; subgroup 1 with the K35/K35 genotype and subgroup 2 with the E35/E35 genotype.

Methods:  Sequence variation within and among full-length assemblies of SRLV subtype A2 subgroups 1 and 2 was analyzed to identify genome-scale recombination patterns and subgroup-specific variants.

Results:  Consensus viral genomes were assembled from 23 infected sheep, including animals of assorted TMEM154 genotypes comprised of haplotypes 1, 2, or 3. Viral genome analysis identified viral subgroups 1 and 2 among the samples, and revealed additional sub-structure within subgroup 2 based on models predicting complex patterns of recombination between the two subgroups in several genomes. Animals with evidence of dual subgroup infection also possessed the most diverse quasi-species and the most highly recombined consensus genomes. After accounting for recombination, 413 subgroup diagnostic single nucleotide polymorphisms (SNPs) were identified.

Conclusions:  The viral subgroup framework developed to classify SRLV consensus genomes along a continuum of recombination suggests that animals with the TMEM154 E35/K35 genotype may represent a reservoir for producing viral genomes representing recombination between A2 subgroups 1 and 2.

MALDI-TOF MS Biomarker Detection Models to Distinguish RTX Toxin Phenotypes of Moraxella bovoculi Strains Are Enhanced Using Calcium Chloride Supplemented Agar

Moraxella bovoculi is the bacterium most often cultured from ocular lesions of cattle with infectious bovine keratoconjunctivitis, also known as bovine pinkeye. Some strains of M. bovoculi contain operons encoding for a repeats-in-toxin (RTX) toxin, which is a known virulence factor of multiple veterinary pathogens. We explored the utility of MALDI-TOF MS and biomarker detection models to classify the presence or absence of an RTX phenotype in M. bovoculi. Ninety strains that had undergone whole genome sequencing were classified by the presence or absence of complete RTX operons and confirmed with a visual assessment of hemolysis on blood agar. Strains were grown on Tryptic Soy Agar (TSA) with 5% sheep blood, TSA with 5% bovine blood that was supplemented with 10% fetal bovine serum, 10 mmol/LCaCl₂, or both. The formulations were designed to determine the influence of growth media on toxin production or activity, as calcium ions are required for toxin secretion and activity. Mass spectra were obtained for strains grown on each agar formulation and biomarker models were developed using ClinProTools 3.0 software. The most accurate model was developed using spectra from strains grown on TSA with 5% bovine blood and supplemented with CaCl₂, which had a sensitivity and specificity of 93.3% and 73.3%, respectively, regarding RTX phenotype classification. The same biomarker model algorithm developed from strains grown on TSA with 5% sheep blood had a substantially lower sensitivity and specificity of 68.0% and 52.0%, respectively. Our results indicate that MALDI-TOF MS biomarker models can accurately classify strains of M. bovoculi regarding the presence or absence of RTX toxin operons and that agar media modifications improve the accuracy of these models.

Classification of small ruminant lentivirus subtype A2, subgroups 1 and 2 based on whole genome comparisons and complex recombination patterns

Background: Small ruminant lentiviruses (SRLVs) cause a multisystemic chronic wasting disease in sheep across much of the world. SRLV subtype A2 is prevalent in North America and further classified into multiple subgroups based on variation in the group antigens gene (gag) and envelope (env) genes. In sheep, the ovine transmembrane protein 154 (TMEM154) gene is associated with SRLV susceptibility. Ewes with at least one copy of TMEM154 encoding a full-length protein with glutamate at position 35 (E35; haplotypes 2 and 3), are highly susceptible to SRLV infection while ewes with any combination of TMEM154 haplotypes which encodes lysine (K35; haplotype 1), or truncated proteins (haplotypes 4 and 6) are several times less so. A2 subgroups 1 and 2 are associated with host TMEM154 genotypes; subgroup 1 with the K35/K35 genotype and subgroup 2 with the E35/E35 genotype. Methods:  Sequence variation within and among full-length assemblies of SRLV subtype A2 subgroups 1 and 2 was analyzed to identify genome-scale recombination patterns and subgroup-specific variants. Results:  Consensus viral genomes were assembled from 23 infected sheep, including animals of assorted TMEM154 genotypes comprised of haplotypes 1, 2, or 3. Viral genome analysis identified viral subgroups 1 and 2 among the samples, and revealed additional sub-structure within subgroup 2 based on models predicting complex patterns of recombination between the two subgroups in several genomes. Animals with evidence of dual subgroup infection also possessed the most diverse quasi-species and the most highly recombined consensus genomes. After accounting for recombination, 413 subgroup diagnostic single nucleotide polymorphisms (SNPs) were identified. Conclusions:  The viral subgroup framework developed to classify SRLV consensus genomes along a continuum of recombination suggests that animals with the TMEM154 E35/K35 genotype may represent a reservoir for producing viral genomes representing recombination between A2 subgroups 1 and 2.

Differences between predicted outer membrane proteins of genotype 1 and 2 Mannheimia haemolytica

BACKGROUND

Mannheimia haemolytica strains isolated from North American cattle have been classified into two genotypes (1 and 2). Although members of both genotypes have been isolated from the upper and lower respiratory tracts of cattle with or without bovine respiratory disease (BRD), genotype 2 strains are much more frequently isolated from diseased lungs than genotype 1 strains. The mechanisms behind the increased association of genotype 2 M. haemolytica with BRD are not fully understood. To address that, and to search for interventions against genotype 2 M. haemolytica, complete, closed chromosome assemblies for 35 genotype 1 and 34 genotype 2 strains were generated and compared. Searches were conducted for the pan genome, core genes shared between the genotypes, and for genes specific to either genotype. Additionally, genes encoding outer membrane proteins (OMPs) specific to genotype 2 M. haemolytica were identified, and the diversity of their protein isoforms was characterized with predominantly unassembled, short-read genomic sequences for up to 1075 additional strains.

RESULTS

The pan genome of the 69 sequenced M. haemolytica strains consisted of 3111 genes, of which 1880 comprised a shared core between the genotypes. A core of 112 and 179 genes or gene variants were specific to genotype 1 and 2, respectively. Seven genes encoding predicted OMPs; a peptidase S6, a ligand-gated channel, an autotransporter outer membrane beta-barrel domain-containing protein (AOMB-BD-CP), a porin, and three different trimeric autotransporter adhesins were specific to genotype 2 as their genotype 1 homologs were either pseudogenes, or not detected. The AOMB-BD-CP gene, however, appeared to be truncated across all examined genotype 2 strains and to likely encode dysfunctional protein. Homologous gene sequences from additional M. haemolytica strains confirmed the specificity of the remaining six genotype 2 OMP genes and revealed they encoded low isoform diversity at the population level.

CONCLUSION

Genotype 2 M. haemolytica possess genes encoding conserved OMPs not found intact in more commensally prone genotype 1 strains. Some of the genotype 2 specific genes identified in this study are likely to have important biological roles in the pathogenicity of genotype 2 M. haemolytica, which is the primary bacterial cause of BRD.

213 Transcriptome profiling of olfactory epithelium in normal cycling and acyclic gilts

A significant proportion of gilts that enter the herd never farrow a litter and are culled because of anestrus or failure to conceive. Genome-wide association studies for pubertal traits have identified loci involved in neuronal and olfactory pathways, and olfaction is critical for expression of reproductive behavior in mammalian females. We evaluated major olfactory epithelium (MOE) for differential gene expression in nonpubertal, behavioral anestrus and normal cycling early follicular and luteal phase gilts (n = 8/group; average age of 259 days). An average of 50 million paired-end RNA-seq reads were collected from each of the 32 RNA libraries and mapped to Sscrofa 11.1. Differential gene expression was determined using DESeq2. A total of 18,484 genes were expressed with a mean normalized expression value greater than 5. Only four genes were differentially expressed between nonpubertal or behavioral anestrus MOE and their cycling controls (early follicular and luteal, respectively). Comparing cycling follicular and luteal phase gilts showed that 1146 genes were more highly expressed in MOE from follicular phase gilts, whereas 1351 genes were more highly expressed in MOE from luteal phase gilts. Pathways for transmembrane receptor protein tyrosine kinase/growth factor signaling, cell junction and epithelium development were overrepresented in follicular phase MOE and cell cycle, chromatin remodeling, inflammatory response and sensory perception were overrepresented in luteal phase MOE. While 1348 locus IDs were identified for olfactory receptors in MOE, only 160 were expressed at an appreciable level (base mean > 5) in MOE and 16 were more highly expressed in MOE from luteal phase than follicular phase gilts. While few genes were differentially expressed in MOE between prepubertal and anestrus gilts and cycling gilts at the same ovarian stage, the comparison between ovarian stages indicates that MOE gene expression is under hormonal control. USDA is an equal opportunity provider and employer.

Next-Generation Sequencing Identifies Polyunsaturated Fatty Acid Responsive Genes in the Juvenile Rat Cerebellum

Dietary n-3 polyunsaturated fatty acids (PUFA) influence postnatal brain growth and development. However, little data exist regarding the impacts of dietary n-3 PUFA in juvenile animals post weaning, which is a time of rapid growth. We tested the hypothesis that depleting dietary n-3 PUFA would result in modifications to the cerebellar transcriptome of juvenile rats. To test this hypothesis, three week old male rats (an age that roughly corresponds to an 11 month old child in brain development) were fed diets containing either soybean oil (SO) providing 1.1% energy from α-linolenic acid (ALA; 18:3n-3; ALA-sufficient) or corn oil (CO) providing 0.13% energy from ALA (ALA-deficient) for four weeks. Fatty acids (FAs) in the cerebellum were analyzed and revealed a 4-fold increase in n-6 docosapentaenoic acid (DPA; 22:5n-6), increases in arachidonic acid (AA; 20:4n-6) and docosatetraenoic acid (DTA; 22:4n-6), but no decrease in docosahexaenoic acid (DHA; 22:6n-3), in animals fed CO versus SO. Transcript abundance was then characterized to identify differentially expressed genes (DEGs) between the two diets. Upper quartile (UQ) scaling and transcripts per million (TPM) data normalization identified 100 and 107 DEGs, respectively. Comparison of DEGs from the two normalization methods identified 70 genes that overlapped, with 90% having abundance differences less than 2-fold. Nr4a3, a transcriptional activator that plays roles in neuroprotection and learning, was elevated over 2-fold from the CO diet. These data indicate that expression of Nr4a3 in the juvenile rat cerebellum is responsive to dietary n-3 PUFA, but additional studies are needed clarify the neurodevelopmental relationships between n-3 PUFA and Nr4a3 and the resulting impacts.

Whole genome sequencing of Moraxella bovoculi reveals high genetic diversity and evidence for interspecies recombination at multiple loci

Moraxella bovoculi is frequently cultured from the ocular secretions and conjunctiva of cattle with Infectious Bovine Keratoconjunctivitis (IBK). Previous work has shown that single nucleotide polymorphism (SNP) diversity in this species is quite high with 81,284 SNPs identified in eight genomes representing two distinct genotypes isolated from IBK affected eyes (genotype 1) and the nasopharynx of cattle without clinical IBK signs (genotype 2), respectively. The goals of this study were to identify SNPs from a collection of geographically diverse and epidemiologically unlinked M. bovoculi strains from the eyes of IBK positive cattle (n = 183) and another from the eyes of cattle (most from a single population at a single time-point) without signs of IBK (n = 63) and to characterize the genetic diversity. Strains of both genotypes were identified from the eyes of cattle without IBK signs. Only genotype 1 strains were identified from IBK affected eyes, however, these strains were isolated before the discovery of genotype 2, and the protocol for their isolation would have preferentially selected genotype 1 M. bovoculi. The core genome comprised ~74% of the whole and contained >127,000 filtered SNPs. More than 80% of these characterize diversity within genotype 1 while 23,611 SNPs (~18%) delimit the two major genotypes. Genotype 2 strains lacked a repeats-in-toxin (RTX) putative pathogenesis factor and any of ten putative antibiotic resistance genes carried within a genomic island. Within genotype 1, prevalence of these elements was 0.85 and 0.12 respectively in strains from eyes that were IBK positive. Recombination appears to be an important source of genetic diversity for genotype 1 and undermines the utility of ribosomal-locus-based species identification. The extremely high genetic diversity in genotype 1 presents a challenge to the development of an efficacious vaccine directed against them, however, several low-diversity pilin-like genes were identified. Finally, the genotype-defining SNPs described in this study are a resource that can facilitate the development of more accurate M. bovoculi diagnostic tests.

Matrix-assisted laser desorption/ionization time-of-flight mass spectrometry identification of Moraxella bovoculi and Moraxella bovis isolates from cattle

Infectious bovine keratoconjunctivitis (IBK) is an economically significant disease caused by Moraxella bovis. Moraxella bovoculi, although not reported to cause IBK, has been isolated from the eyes of cattle diagnosed with IBK. Identification of M. bovis and M. bovoculi can be performed using biochemical or DNA-based approaches, both of which may be time consuming and inconsistent between laboratories. We conducted a comparative evaluation of M. bovoculi and M. bovis identification using matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOF MS) with a database provided by Bruker Daltonics (termed the BDAL database), the BDAL database supplemented with spectra generated in our study (termed the UNLVDC database), and with PCR-restriction-fragment length polymorphism (PCR-RFLP) typing. M. bovoculi ( n = 250) and M. bovis ( n = 18) isolates from cattle with or without IBK were used. MALDI-TOF MS using the UNLVDC database correctly identified 250 of 250 (100%) of M. bovoculi and 17 of 18 (94%) of M. bovis isolates. With the BDAL database, MALDI-TOF MS correctly identified 249 of 250 (99%) of M. bovoculi and 7 of 18 (39%) of M. bovis isolates. In comparison, the PCR-RFLP test correctly identified 210 of 250 (84%) of M. bovoculi and 12 of 18 (66%) of M. bovis isolates. Thus, MALDI-TOF MS with the UNLVDC database was the most effective identification methodology for M. bovis and M. bovoculi isolates from cattle.

Unexpected High Intragenomic Variation in Two of Three Major Pest Thrips Species Does Not Affect Ribosomal Internal Transcribed Spacer 2 (ITS2) Utility for Thrips Identification

The mitochondrial cytochrome oxidase I gene (mtCO1) and the ribosomal internal transcribed spacer 2 region (ITS2) are among the most widely used molecular markers for insect taxonomic characterization. Three economically important species of thrips, Scirtothripsdorsalis, Thripspalmi, and Frankliniellaoccidentalis were selected to examine the extent of intragenomic variation within these two marker regions in the family Thripidae, and determine if this variation would affect the utility of markers in thrips molecular diagnostics. For each species, intragenomic (within individual) variation and intergenomic (among individuals) variation was assessed by cloning and sequencing PCR-amplified copies. Intergenomic variation was generally higher than intragenomic variation except in cases where intergenomic variation was very low, as in mtCO1 from S.dorsalis and F.occidentalis. Intragenomic variation was detected in both markers in all three of the thrips species, however, 2-3 times more intragenomic variation was observed for ITS2 than mtCO1 in both S.dorsalis and T.palmi. Furthermore, levels of intragenomic variation were low for both of the genes in F.occidentalis. In all of the three thrips species, no sex-based clustering of haplotypes was observed in either marker. Unexpected high intragenomic variation in ITS2 for two of three thrips species did not interfere with thrips diagnostics. However, caution should be taken in applying ITS2 to certain studies of S.dorsalis and T.palmi when high levels of intragenomic variation could be problematic or confounding. In such studies, mtCO1 may be a preferable marker. Possible reasons for discrepancies in intragenomic variation among genomic regions are discussed.

Rediscovering digitules in Aphidomorpha and the question of homology among Sternorrhyncha (Insecta, Hemiptera)

We explore and expand on the morphological term digitule. The term was originally proposed for toe-like setae on a species of Phylloxera Boyer de Fonscolombe, 1834 (Hemiptera, Sternorrhyncha, Aphidomorpha) by Henry Shimer, an American naturalist. While it is standard terminology in scale systematics (Hemiptera, Sternorrhyncha, Coccidomorpha), the term digitule was ignored by aphid specialists despite being the original taxon for which the term was described. Similar setae occur on many arthropod groups, so the homology is poorly understood even within any superfamily of Hemiptera. We provide the etymology of the term, a proposed explanation for why it was used among scale taxonomists and not aphid taxonomists, and discuss briefly options to progress beyond the confusion between terminology for morphology and homology in Sternorrhyncha.

Genomic signatures of Mannheimia haemolytica that associate with the lungs of cattle with respiratory disease, an integrative conjugative element, and antibiotic resistance genes

BACKGROUND

Mannheimia haemolytica typically resides in cattle as a commensal member of the upper respiratory tract microbiome. However, some strains can invade their lungs and cause respiratory disease and death, including those with multi-drug resistance. A nucleotide polymorphism typing system was developed for M. haemolytica from the genome sequences of 1133 North American isolates, and used to identify genetic differences between isolates from the lungs and upper respiratory tract of cattle with and without clinical signs of respiratory disease.

RESULTS

A total of 26,081 nucleotide polymorphisms were characterized after quality control filtering of 48,403 putative polymorphisms. Phylogenetic analyses of nucleotide polymorphism genotypes split M. haemolytica into two major genotypes (1 and 2) that each were further divided into multiple subtypes. Multiple polymorphisms were identified with alleles that tagged genotypes 1 or 2, and their respective subtypes. Only genotype 2 M. haemolytica associated with the lungs of diseased cattle and the sequence of a particular integrative and conjugative element (ICE). Additionally, isolates belonging to one subtype of genotype 2 (2b), had the majority of antibiotic resistance genes detected in this study, which were assorted into seven combinations that ranged from 1 to 12 resistance genes.

CONCLUSIONS

Typing of diverse M. haemolytica by nucleotide polymorphism genotypes successfully identified associations with diseased cattle lungs, ICE sequence, and antibiotic resistance genes. Management of cattle by their carriage of M. haemolytica could be an effective intervention strategy to reduce the prevalence of respiratory disease and supplemental needs for antibiotic treatments in North American herds.

Large genomic differences between Moraxella bovoculi isolates acquired from the eyes of cattle with infectious bovine keratoconjunctivitis versus the deep nasopharynx of asymptomatic cattle

Moraxella bovoculi is a recently described bacterium that is associated with infectious bovine keratoconjunctivitis (IBK) or “pinkeye” in cattle. In this study, closed circularized genomes were generated for seven M. bovoculi isolates: three that originated from the eyes of clinical IBK bovine cases and four from the deep nasopharynx of asymptomatic cattle. Isolates that originated from the eyes of IBK cases profoundly differed from those that originated from the nasopharynx of asymptomatic cattle in genome structure, gene content and polymorphism diversity and consequently placed into two distinct phylogenetic groups. These results suggest that there are genetically distinct strains of M. bovoculi that may not associate with IBK.

A novel mitochondrial genome architecture in thrips (Insecta: Thysanoptera): extreme size asymmetry among chromosomes and possible recent control region duplication

Background

Multipartite mitochondrial genomes are very rare in animals but have been found previously in two insect orders with highly rearranged genomes, the Phthiraptera (parasitic lice), and the Psocoptera (booklice/barklice).

Results

We provide the first report of a multipartite mitochondrial genome architecture in a third order with highly rearranged genomes: Thysanoptera (thrips). We sequenced the complete mitochondrial genomes of two divergent members of the Scirtothrips dorsalis cryptic species complex. The East Asia 1 species has the single circular chromosome common to animals while the South Asia 1 species has a genome consisting of two circular chromosomes. The fragmented South Asia 1 genome exhibits extreme chromosome size asymmetry with the majority of genes on the large, 14.28 kb, chromosome and only nad6 and trnC on the 0.92 kb mini-circle chromosome. This genome also features paralogous control regions with high similarity suggesting a very recent origin of the nad6 mini-circle chromosome in the South Asia 1 cryptic species.

Conclusions

Thysanoptera, along with the other minor paraenopteran insect orders should be considered models for rapid mitochondrial genome evolution, including fragmentation. Continued use of these models will facilitate a greater understanding of recombination and other mitochondrial genome evolutionary processes across eukaryotes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1672-4) contains supplementary material, which is available to authorized users.

The Scirtothrips dorsalis Species Complex: Endemism and Invasion in a Global Pest

Invasive arthropods pose unique management challenges in various environments, the first of which is correct identification. This apparently mundane task is particularly difficult if multiple species are morphologically indistinguishable but accurate identification can be determined with DNA barcoding provided an adequate reference set is available. Scirtothrips dorsalis is a highly polyphagous plant pest with a rapidly expanding global distribution and this species, as currently recognized, may be comprised of cryptic species. Here we report the development of a comprehensive DNA barcode library for S. dorsalis and seven nuclear markers via next-generation sequencing for identification use within the complex. We also report the delimitation of nine cryptic species and two morphologically distinguishable species comprising the S. dorsalis species complex using histogram analysis of DNA barcodes, Bayesian phylogenetics, and the multi-species coalescent. One member of the complex, here designated the South Asia 1 cryptic species, is highly invasive, polyphagous, and likely the species implicated in tospovirus transmission. Two other species, South Asia 2, and East Asia 1 are also highly polyphagous and appear to be at an earlier stage of global invasion. The remaining members of the complex are regionally endemic, varying in their pest status and degree of polyphagy. In addition to patterns of invasion and endemism, our results provide a framework both for identifying members of the complex based on their DNA barcode, and for future species delimiting efforts.

ESTIMATING BACTERIAL DIVERSITY IN SCIRTOTHRIPS DORSALIS (THYSANOPTERA: THRIPIDAE) VIA NEXT GENERATION SEQUENCING

The last 2 decades have produced a better understanding of insect-microbial associations and yielded some important opportunities for insect control. However, most of our knowledge comes from model systems. Thrips (Thysanoptera: Thripidae) have been understudied despite their global importance as invasive species, plant pests and disease vectors. Using a culture and primer independent next-generation sequencing and metagenomics pipeline, we surveyed the bacteria of the globally important pest, Scirtothrips dorsalis Hood. The most abundant bacterial phyla identified were Actinobacteria and Proteobacteria and the most abundant genera were Propionibacterium, Stenotrophomonas, and Pseudomonas. A total of 189 genera of bacteria were identified. The absence of any vertically transferred symbiont taxa commonly found in insects is consistent with other studies suggesting that thrips primarilly acquire resident microbes from their environment. This does not preclude a possible beneficial/intimate association between S. dorsalis and the dominant taxa identified and future work should determine the nature of these associations.

Evolution and homoplasy at the Bem6 microsatellite locus in three sweetpotato whitefly (Bemisia tabaci) cryptic species

BACKGROUND

The evolution of individual microsatellite loci is often complex and homoplasy is common but often goes undetected. Sequencing alleles at a microsatellite locus can provide a more complete picture of the common evolutionary mechanisms occurring at that locus and can reveal cases of homoplasy. Within species homoplasy can lead to an underestimate of differentiation among populations and among species homoplasy can produce a misleading interpretation regarding shared alleles and hybridization. This is especially problematic with cryptic species.

RESULTS

By sequencing alleles from three cryptic species of the sweetpotato whitefly (Bemisia tabaci), designated MEAM1, MED, and NW, the evolution of the putatively dinucleotide Bem6 (CA₈)imp microsatellite locus is inferred as one of primarily stepwise mutation occurring at four distinct heptaucleotide tandem repeats. In two of the species this pattern yields a compound tandem repeat. Homoplasy was detected both among species and within species.

CONCLUSIONS

In the absence of sequencing, size homoplasious alleles at the Bem6 locus lead to an overestimate of alleles shared and hybridization among cryptic species of Bemisia tabaci. Furthermore, the compound heptanucleotide motif structure of a putative dinucleotide microsatellite has implications for the nomenclature of heptanucleotide tandem repeats with step-wise evolution.

Population genetics of invasive Bemisia tabaci (Hemiptera: Aleyrodidae) cryptic species in the United States based on microsatellite markers

The Bemisia tabaci (Gennadius) (Hemiptera: Aleyrodidae) cryptic species complex of whiteflies contains two species, MEAM1 and MED, that are highly invasive in supportive climates the world over. In the United States, MEAM1 occurs both in the field and in the greenhouse, but MED is only found in the greenhouse. To make inferences about the population structure of both species, and the origin and recent spread of MED within the United States, 987 MEAM1 whiteflies and 340 MED whiteflies were genotyped at six and seven microsatellite loci, respectively, for population genetic analyses. Major results of the study are 1) MED exhibits more population structure and genetic differentiation than MEAM1, 2) nuclear microsatellite markers exhibit a high degree of concordance with mitochondrial markers recovering a major east-west phylogeographic break within MED, 3) both eastern and western MED are found throughout the continental United States and eastern MED is present in Hawaii, and 4) MEAM1 contains two greenhouse U.S. populations significantly differentiated from other U.S. MEAM1. The results suggest that MED was introduced into the United States on at least three occasions and rapidly spread throughout the United States, showing no discernible differentiation across 7,000 km. The results further suggest that there is an enhanced role of the protected agricultural environment in promoting genetic differentiation in both invasive B. tabaci cryptic species.

Distribution of Bemisia tabaci (Hemiptera: Aleyrodidae) biotypes in North America after the Q invasion

After the 2004 discovery of the Bemisia tabaci (Gennadius) (Hemiptera Aleyrodidae) Q biotype in the United States, there was a vital need to determine the geographical and host distribution as well as its interaction with the resident B biotype because of its innate ability to rapidly develop high-level insecticide resistance that persists in the absence of exposure. As part of a coordinated country-wide effort, an extensive survey of B. tabaci biotypes was conducted in North America, with the cooperation of growers, industry, local, state, and federal agencies, to monitor the introduction and distribution of the Q biotype. The biotype status of submitted B. tabaci samples was determined either by polymerase chain reaction amplification and sequencing of a mitochondrial cytochrome oxidase I small subunit gene fragment and characterization of two biotype discriminating nuclear microsatellite markers or esterase zymogram analysis. Two hundred and eighty collections were sampled from the United States, Bermuda, Canada, and Mexico during January 2005 through December 2011. Host plants were split between ornamental plant and culinary herb (67%) and vegetable and field crop (33%) commodities. The New World biotype was detected on field-grown tomatoes (Solanum lycopersicum L.) in Mexico (two) and in commercial greenhouses in Texas (three) and represented 100% of these five collections. To our knowledge, the latter identification represents the first report of the New World biotype in the United States since its rapid displacement in the late 1980s after the introduction of biotype B. Seventy-one percent of all collections contained at least one biotype B individual, and 53% of all collections contained only biotype B whiteflies. Biotype Q was detected in 23 states in the United States, Canada (British Columbia and Ontario territories), Bermuda, and Mexico. Forty-five percent of all collections were found to contain biotype Q in samples from ornamentals, herbs and a single collection from tomato transplants located in protected commercial horticultural greenhouses, but there were no Q detections in outdoor agriculture (vegetable or field crops). Ten of the 15 collections (67%) from Canada and a single collection from Bermuda contained biotype Q, representing the first reports of biotype Q for both countries. Three distinct mitochondrial haplotypes of B. tabaci biotype Q whiteflies were detected in North America Our data are consistent with the inference of independent invasions from at least three different locations. Of the 4,641 individuals analyzed from 517 collections that include data from our previous work, only 16 individuals contained genetic or zymogram evidence of possible hybridization of the Q and B biotypes, and there was no evidence that rare hybrid B-Q marker co-occurrences persisted in any populations.

Population structure of Phthorimaea operculella (Lepidoptera: Gelechiidae) in the United States

The potato tuberworm, Phthorimaea operculella Zeller, is one of the most damaging pests of potatoes worldwide. Although tuberworm was recorded in California as early as 1856, it was first reported in damaging numbers in Oregon and Washington in the early 2000s. The objective of this study was to provide baseline information on the population structure of potato tuberworm in the United States. Potato tuberworm adults were collected from potato fields in the major production regions in the United States. Amplified fragment length polymorphisms (AFLP) were used to determine the genetic population structure of potato tuberworm. We used 219 individuals and 335 polymorphic AFLP bands to infer the geographic population structure of potato tuberworm. Posterior probabilities calculated using the model-based clustering method implemented in STRUCTURE, and F(st) values calculated using AFLP-SURV confirmed the existence of two distinct populations in the United States (one mostly associated to the western United States and the other one mostly associated to the eastern United States). This study provides baseline data for the molecular characterization of potato tuberworm populations, which will aid in tracking the origin of future invasions within the United States.