Whole genome sequencing of Moraxella bovis strains from North America reveals two genotypes with different genetic determinants.

BACKGROUND: Moraxella bovis and Moraxella bovoculi both associate with infectious bovine keratoconjunctivitis (IBK), an economically significant and painful ocular disease that affects cattle worldwide. There are two genotypes of M. bovoculi (genotypes 1 and 2) that differ in their gene content and potential virulence factors, although neither have been experimentally shown to cause IBK. M. bovis is a causative IBK agent, however, not all strains carry a complete assortment of known virulence factors. The goals of this study were to determine the population structure and depth of M. bovis genomic diversity, and to compare core and accessory genes and predicted outer membrane protein profiles both within and between M. bovis and M. bovoculi. RESULTS: Phylogenetic trees and bioinformatic analyses of 36 M. bovis chromosomes sequenced in this study and additional available chromosomes of M. bovis and both genotype 1 and 2 M. bovoculi, showed there are two genotypes (1 and 2) of M. bovis. The two M. bovis genotypes share a core of 2015 genes, with 121 and 186 genes specific to genotype 1 and 2, respectively. The two genotypes differ by their chromosome size and prophage content, encoded protein variants of the virulence factor hemolysin, and by their affiliation with different plasmids. Eight plasmid types were identified in this study, with types 1 and 6 observed in 88 and 56% of genotype 2 strains, respectively, and absent from genotype 1 strains. Only type 1 plasmids contained one or two gene copies encoding filamentous haemagglutinin-like proteins potentially involved with adhesion. A core of 1403 genes was shared between the genotype 1 and 2 strains of both M. bovis and M. bovoculi, which encoded a total of nine predicted outer membrane proteins. CONCLUSIONS: There are two genotypes of M. bovis that differ in both chromosome content and plasmid profiles and thus may not equally associate with IBK. Immunological reagents specifically targeting select genotypes of M. bovis, or all genotypes of M. bovis and M. bovoculi together could be designed from the outer membrane proteins identified in this study.

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.

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.

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.

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.

The GEA pipeline for characterizing Escherichia coli and Salmonella genomes.

Salmonella enterica and Escherichia coli are major food-borne human pathogens, and their genomes are routinely sequenced for clinical surveillance. Computational pipelines designed for analyzing pathogen genomes should both utilize the most current information from annotation databases and increase the coverage of these databases over time. We report the development of the GEA pipeline to analyze large batches of E. coli and S. enterica genomes. The GEA pipeline takes as input paired Illumina raw reads files which are then assembled followed by annotation. Alternatively, assemblies can be provided as input and directly annotated. The pipeline provides predictive genome annotations for E. coli and S. enterica with a focus on the Center for Genomic Epidemiology tools. Annotation results are provided as a tab delimited text file. The GEA pipeline is designed for large-scale E. coli and S. enterica genome assembly and characterization using the Center for Genomic Epidemiology command-line tools and high-performance computing. Large scale annotation is demonstrated by an analysis of more than 14,000 Salmonella genome assemblies. Testing the GEA pipeline on E. coli raw reads demonstrates reproducibility across multiple compute environments and computational usage is optimized on high performance computers.

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 Extremely Heat-Resistant Escherichia coli in the Beef Production and Processing Continuum.

Understanding the dynamics of stress-resistant Escherichia coli (E. coli) across the meat production and processing continuum is important for tracking sources of such microbes and devising effective modes of control. The Locus of Heat Resistance (LHR) is a ∼14-19 Kb genetic element imparting extreme heat resistance (XHR) in Enterobacteriaceae. It has been hypothesized that thermal and antimicrobial interventions applied during meat processing may select for LHR+E. coli. Thus, our goal was to study the prevalence and molecular biology of LHR+E. coli among lots of beef cattle (n = 3) from production through processing. Two hundred thirty-two generic E. coli isolated from the same animals through seven stages of the beef processing continuum (cattle in feedyards to packaged strip loins) were examined. LHR+E. coli were rare (0.6%; 1 of 180) among the early stages of the beef continuum (feces and hides at feedlot, feces and hides at harvest, and preevisceration carcasses), whereas the prevalence of LHR+E. coli on final carcasses and strip loins was remarkably higher. Half (14 of 28) of the final carcass E. coli possessed the LHR, while 79.2% (19 of 24) of the strip loin E. coli did. Eighty-five percent (29 of 34) of the LHR+E. coli presented with the XHR phenotype. The selection or enrichment of LHR+E. coli from harvest steps to the final products appeared unlikely as the LHR+E. coli isolates were effectively controlled by antimicrobial interventions typically used during beef processing. Further, whole-genome sequencing of the isolates suggested LHR+E. coli are persisting in the chilled processing environment and that horizontal LHR transfer among E. coli isolates may take place.

Probiotics in milk replacer affect the microbiome of the lung in neonatal dairy calves.

Introduction: Probiotics have been investigated for their many health benefits and impact on the microbiota of the gut. Recent data have also supported a gut-lung axis regarding the bacterial populations (microbiomes) of the two locations; however, little research has been performed to determine the effects of oral probiotics on the microbiome of the bovine respiratory tract. We hypothesized that probiotic treatment would result in changes in the lung microbiome as measured in lung lavage fluid. Our overall goal was to characterize bacterial populations in the lungs of calves fed probiotics in milk replacer and dry rations from birth to weaning. Methods: A group of 20 dairy calves was split into two treatment groups: probiotic (TRT; N = 10, milk replacer +5 g/d probiotics; Bovamine Dairy, Chr. Hansen, Inc., Milwaukee, WI) and control (CON; N = 10, milk replacer only). On day 0, birth weight was obtained, and calves were provided colostrum as per the dairy SOP. On day 2, probiotics were added to the milk replacer of the treated group and then included in their dry ration. Lung lavages were performed on day 52 on five random calves selected from each treatment group. DNA was extracted from lavage fluid, and 16S ribosomal RNA (rRNA) gene hypervariable regions 1-3 were amplified by PCR and sequenced using next-generation sequencing (Illumina MiSeq) for the identification of the bacterial taxa present. Taxa were classified into both operational taxonomic units (OTUs) and amplicon sequence variants (ASVs). Results: Overall, the evaluation of these samples revealed that the bacterial genera identified in the lung lavage samples of probiotic-fed calves as compared to the control calves were significantly different based on the OTU dataset (p < 0.05) and approached significance for the ASV dataset (p < 0.06). Additionally, when comparing the diversity of taxa in lung lavage samples to nasal and tonsil samples, taxa diversity of lung samples was significantly lower (p < 0.05). Discussion: In conclusion, analysis of the respiratory microbiome in lung lavage samples after probiotic treatment provides insight into the distribution of bacterial populations in response to oral probiotics and demonstrates that oral probiotics affect more than the gut microbiome.

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.

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.

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/LCaCl2, 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 CaCl2, 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.

Rapid differentiation of Moraxella bovoculi genotypes 1 and 2 using MALDI-TOF mass spectrometry profiles.

Moraxella bovoculi is the most frequently isolated bacteria from the eyes of cattle with Infectious Bovine Keratoconjunctivitis (IBK), also known as bovine pinkeye. Two distinct genotypes of M. bovoculi, genotype 1 and genotype 2, were characterized after whole genome sequencing showed a large degree of single nucleotide polymorphism (SNP) diversity within the species. To date, both genotypes have been isolated from the eyes of cattle without clinical signs of IBK while only genotype 1 strains have been isolated from the eyes of cattle with clinical signs of IBK. We used 38 known genotype 1 strains and 26 known genotype 2 strains to assess the ability of matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOF MS) to accurately genotype M. bovoculi strains using mass spectrum biomarkers. Mass spectrum data was analyzed with ClinProTools 3.0 software and six models were developed that classify strain genotypes with accuracies ranging from 90.6% - 100%. Finally, using four of the most genotype-specific peaks that also exhibited high peak intensities from the six automated models, we developed a customized model (UNL assisted model) that had recognition capability, validation, and classification accuracies of 100% for genotype classification. Our results indicate that MALDI-TOF MS biomarkers can be used to accurately discriminate genotypes of M. bovoculi without the need for additional methods.

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.

Digestive tract microbiota of beef cattle that differed in feed efficiency.

We hypothesized cattle that differed in BW gain had different digestive tract microbiota. Two experiments were conducted. In both experiments, steers received a diet that consisted of 8.0% chopped alfalfa hay, 20% wet distillers grain with solubles, 67.75% dry-rolled corn, and 4.25% vitamin/mineral mix (including monensin) on a dry matter basis. Steers had ad libitum access to feed and water. In experiment 1, 144 steers (age = 310 ± 1.5 d; BW = 503 ± 37.2 kg) were individually fed for 105 d. Ruminal digesta samples were collected from eight steers with the greatest (1.96 ± 0.02 kg/d) and eight steers with the least ADG (1.57 ± 0.02 kg/d) that were within ±0.32 SD of the mean (10.1 ± 0.05 kg/d) dry matter. In experiment 2, 66 steers (age = 396 ± 1 d; BW = 456 ± 5 kg) were individually fed for 84 d. Rumen, duodenum, jejunum, ileum, cecum, and colon digesta samples were collected from eight steers with the greatest (2.39 ± 0.06 kg/d) and eight steers with the least ADG (1.85 ± 0.06 kg/d) that were within ±0.55 SD of the mean dry matter intake (11.9 ± 0.1 kg/d). In both studies, DNA was isolated and the V1 to V3 regions of the 16S rRNA gene were sequenced. Operational taxonomic units were classified using 0.03 dissimilarity and identified using the Greengenes 16S rRNA gene database. In experiment 1, there were no differences in the Chao1, Shannon, Simpson, and InvSimpson diversity indexes or the permutation multivariate analysis of variance (PERMANOVA; P = 0.57). The hierarchical test returned six clades as being differentially abundant between steer classifications (P < 0.05). In experiment 2, Chao1, Shannon, Simpson, and InvSimpson diversity indexes and PERMANOVA between steer classified as less or greater ADG did not differ (P > 0.05) for the rumen, duodenum, ileum, cecum, and colon. In the jejunum, there tended to be a difference in the Chao1 (P = 0.09) and Simpson diversity (P = 0.09) indexes between steer classifications, but there was no difference in the Shannon (P = 0.14) and InvSimpson (P = 0.14) diversity indexes. Classification groups for the jejunum differed (P = 0.006) in the PERMANOVA. The hierarchical dependence false discovery rate procedure returned 11 clades as being differentially abundant between steer classifications in the jejunum (P < 0.05). The majority of the OTU were in the Families Corynebacteriaceae and Coriobacteriaceae. This study suggests that intestinal differences in the microbiota of ruminants may be associated with animal performance.

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.

Complete Genome Sequence of Moraxella bovis Strain Epp-63 (300), an Etiologic Agent of Infectious Bovine Keratoconjunctivitis.

We report here the complete closed genome sequence of Moraxella bovis strain Epp-63 (300) (Epp63). This strain was isolated from an infectious bovine keratoconjunctivitis (IBK) case in 1963. Since then, Epp63 has been used extensively for IBK research. Consequently, the genome sequence of Epp63 should help elucidate IBK host-pathogen interactions.