Detection and characterization of novel luchacoviruses, genus Alphacoronavirus, in saliva and feces of meso-carnivores in the northeastern United States.

IMPORTANCE: Several coronaviruses (CoVs) have been detected in domesticated, farmed, and wild meso-carnivores, causing a wide range of diseases and infecting diverse species, highlighting their important but understudied role in the epidemiology of these viruses. Assessing the viral diversity hosted in wildlife species is essential to understand their significance in the cross-species transmission of CoVs. Our focus here was on CoV discovery in meso-carnivores in the Northeast United States as a potential "hotspot" area with high density of humans and urban wildlife. This study identifies novel alphacoronaviruses circulating in multiple free-ranging wild and domestic species in this area and explores their potential epidemiological importance based on regions of the Spike gene, which are relevant for virus-host interactions.

Evolutionary genomic analyses of canine E. coli infections identify a relic capsular locus associated with resistance to multiple classes of antimicrobials.

Infections caused by antimicrobial-resistant Escherichia coli are the leading cause of death attributed to antimicrobial resistance (AMR) worldwide, and the known AMR mechanisms involve a range of functional proteins. Here, we employed a pan-genome wide association study (GWAS) approach on over 1,000 E. coli isolates from sick dogs collected across the US and Canada and identified a strong statistical association (empirical P < 0.01) of AMR, involving a range of antibiotics to a group 1 capsular (CPS) gene cluster. This cluster included genes under relaxed selection pressure, had several loci missing, and had pseudogenes for other key loci. Furthermore, this cluster is widespread in E. coli and Klebsiella clinical isolates across multiple host species. Earlier studies demonstrated that the octameric CPS polysaccharide export protein Wza can transmit macrolide antibiotics into the E. coli periplasm. We suggest that the CPS in question, and its highly divergent Wza, functions as an antibiotic trap, preventing antimicrobial penetration. We also highlight the high diversity of lineages circulating in dogs across all regions studied, the overlap with human lineages, and regional prevalence of resistance to multiple antimicrobial classes. IMPORTANCE: Much of the human genomic epidemiology data available for E. coli mechanism discovery studies has been heavily biased toward shiga-toxin producing strains from humans and livestock. E. coli occupies many niches and produces a wide variety of other significant pathotypes, including some implicated in chronic disease. We hypothesized that since dogs tend to share similar strains with their owners and are treated with similar antibiotics, their pathogenic isolates will harbor unexplored AMR mechanisms of importance to humans as well as animals. By comparing over 1,000 genomes with in vitro antimicrobial susceptibility data from sick dogs across the US and Canada, we identified a strong multidrug resistance association with an operon that appears to have once conferred a type 1 capsule production system.

White shark genome reveals ancient elasmobranch adaptations associated with wound healing and the maintenance of genome stability.

The white shark (Carcharodon carcharias; Chondrichthyes, Elasmobranchii) is one of the most publicly recognized marine animals. Here we report the genome sequence of the white shark and comparative evolutionary genomic analyses to the chondrichthyans, whale shark (Elasmobranchii) and elephant shark (Holocephali), as well as various vertebrates. The 4.63-Gbp white shark genome contains 24,520 predicted genes, and has a repeat content of 58.5%. We provide evidence for a history of positive selection and gene-content enrichments regarding important genome stability-related genes and functional categories, particularly so for the two elasmobranchs. We hypothesize that the molecular adaptive emphasis on genome stability in white and whale sharks may reflect the combined selective pressure of large genome sizes, high repeat content, high long-interspersed element retrotransposon representation, large body size, and long lifespans, represented across these two species. Molecular adaptation for wound healing was also evident, with positive selection in key genes involved in the wound-healing process, as well as Gene Ontology enrichments in fundamental wound-healing pathways. Sharks, particularly apex predators such as the white shark, are believed to have an acute sense of smell. However, we found very few olfactory receptor genes, very few trace amine-associated receptors, and extremely low numbers of G protein-coupled receptors. We did however, identify 13 copies of vomeronasal type 2 (V2R) genes in white shark and 10 in whale shark; this, combined with the over 30 V2Rs reported previously for elephant shark, suggests this gene family may underlie the keen odorant reception of chondrichthyans.

Evolutionary genomic and bacteria GWAS analysis of Mycobacterium avium subsp. paratuberculosis and dairy cattle Johne's disease phenotypes.

Mycobacterium avium subsp. paratuberculosis (MAP) is the causative agent of Johne's disease in ruminants, which has important health consequences for dairy cattle. The Regional Dairy Quality Management Alliance (RDQMA) project is a multistate research program involving MAP isolates taken from three intensively studied commercial dairy farms in the northeastern United States, which emphasized longitudinal data collection of both MAP isolates and animal health in three regional dairy herds for a period of about 7 years. This paper reports the results of a pan-GWAS analysis involving 318 MAP isolates and dairy cow Johne's disease phenotypes, taken from these three farms. Based on our highly curated accessory gene count the pan-GWAS analysis identified several MAP genes associated with bovine Johne's disease phenotypes scored from these three farms, with some of the genes having functions suggestive of possible cause/effect relationships to these phenotypes. This paper reports a pan-genomic comparative analysis between MAP and Mycobacterium tuberculosis, assessing functional Gene Ontology category enrichments between these taxa. Finally, we also provide a population genomic perspective on the effectiveness of herd isolation, involving closed dairy farms, in preventing MAP inter-farm cross infection on a micro-geographic scale.IMPORTANCE Mycobacterium avium subsp. paratuberculosis (MAP) is the causative agent of Johne's disease in ruminants, which has important health consequences for dairy cattle, and enormous economic consequences for the dairy industry. Understanding which genes in this bacterium are correlated with key disease phenotypes can lead to functional experiments targeting these genes and ultimately lead to improved control strategies. This study represents a rare example of a prolonged longitudinal study of dairy cattle where the disease was measured and the bacteria were isolated from the same cows. The genome sequences of over 300 MAP isolates were analyzed for genes that were correlated with a wide range of Johne's disease phenotypes. A number of genes were identified that were significantly associated with several aspects of the disease and suggestive of further experimental follow-up.

Genomic Assessment of Global Population Structure in a Highly Migratory and Habitat Versatile Apex Predator, the Tiger Shark (Galeocerdo cuvier).

Understanding the population dynamics of highly mobile, widely distributed, oceanic sharks, many of which are overexploited, is necessary to aid their conservation management. We investigated the global population genomics of tiger sharks (Galeocerdo cuvier), a circumglobally distributed, apex predator displaying remarkable behavioral versatility in its diet, habitat use (near coastal, coral reef, pelagic), and individual movement patterns (spatially resident to long-distance migrations). We genotyped 242 tiger sharks from 10 globally distributed locations at more than 2000 single nucleotide polymorphisms. Although this species often conducts massive distance migrations, the data show strong genetic differentiation at both neutral (FST = 0.125-0.144) and candidate outlier loci (FST = 0.570-0.761) between western Atlantic and Indo-Pacific sharks, suggesting the potential for adaptation to the environments specific to these oceanic regions. Within these regions, there was mixed support for population differentiation between northern and southern hemispheres in the western Atlantic, and none for structure within the Indian Ocean. Notably, the results demonstrate a low level of population differentiation of tiger sharks from the remote Hawaiian archipelago compared with sharks from the Indian Ocean (FST = 0.003-0.005, P < 0.01). Given concerns about biodiversity loss and marine ecosystem impacts caused by overfishing of oceanic sharks in the midst of rapid environmental change, our results suggest it imperative that international fishery management prioritize conservation of the evolutionary potential of the highly genetically differentiated Atlantic and Indo-Pacific populations of this unique apex predator. Furthermore, we suggest targeted management attention to tiger sharks in the Hawaiian archipelago based on a precautionary biodiversity conservation perspective.

Population Gene Introgression and High Genome Plasticity for the Zoonotic Pathogen Streptococcus agalactiae.

The influence that bacterial adaptation (or niche partitioning) within species has on gene spillover and transmission among bacterial populations occupying different niches is not well understood. Streptococcus agalactiae is an important bacterial pathogen that has a taxonomically diverse host range making it an excellent model system to study these processes. Here, we analyze a global set of 901 genome sequences from nine diverse host species to advance our understanding of these processes. Bayesian clustering analysis delineated 12 major populations that closely aligned with niches. Comparative genomics revealed extensive gene gain/loss among populations and a large pan genome of 9,527 genes, which remained open and was strongly partitioned among niches. As a result, the biochemical characteristics of 11 populations were highly distinctive (significantly enriched). Positive selection was detected and biochemical characteristics of the dispensable genes under selection were enriched in ten populations. Despite the strong gene partitioning, phylogenomics detected gene spillover. In particular, tetracycline resistance (which likely evolved in the human-associated population) from humans to bovine, canines, seals, and fish, demonstrating how a gene selected in one host can ultimately be transmitted into another, and biased transmission from humans to bovines was confirmed with a Bayesian migration analysis. Our findings show high bacterial genome plasticity acting in balance with selection pressure from distinct functional requirements of niches that is associated with an extensive and highly partitioned dispensable genome, likely facilitating continued and expansive adaptation.

Transcriptome-Derived Microsatellites Demonstrate Strong Genetic Differentiation in Pacific White Sharks.

Recent advances in genome-scale sequencing technology have allowed the development of high resolution genetic markers for the study of nonmodel taxa. In particular, transcriptome sequencing has proven to be highly useful in generating genomic markers for use in population genetic studies, allowing for insight into species connectivity, as well as local adaptive processes as many transcriptome-derived markers are found within or associated with functional genes. Herein, we developed a set of 30 microsatellite markers from a heart transcriptome for the white shark (Carcharodon carcharias), a widely distributed and globally vulnerable marine predator. Using these markers as well as 10 published anonymous genomic microsatellite loci, we provide 1) the first nuclear genetic assessment of the cross-Pacific connectivity of white sharks, and 2) a comparison of the levels of inferred differentiation across microsatellite marker sets (i.e., transcriptome vs. anonymous) to assess their respective utility to elucidate the population genetic dynamics of white sharks. Significant (FST = 0.083, P = 0.05; G″ST = 0.200; P = 0.001) genetic differentiation was found between Southwestern Pacific (n = 19) and Northeastern Pacific (n = 20) white sharks, indicating restricted, cross Pacific gene flow in this species. Transcriptome-derived microsatellite marker sets identified much higher (up to 2×) levels of genetic differentiation than anonymous genomic markers, underscoring potential utility of transcriptome markers in identifying subtle population genetic differences within highly vagile, globally distributed marine species.Subject areas: Population structure and phylogeography; Conservation genetics and biodiversity.

Comparative transcriptomics of elasmobranchs and teleosts highlight important processes in adaptive immunity and regional endothermy.

BACKGROUND: Comparative genomic and/or transcriptomic analyses involving elasmobranchs remain limited, with genome level comparisons of the elasmobranch immune system to that of higher vertebrates, non-existent. This paper reports a comparative RNA-seq analysis of heart tissue from seven species, including four elasmobranchs and three teleosts, focusing on immunity, but concomitantly seeking to identify genetic similarities shared by the two lamnid sharks and the single billfish in our study, which could be linked to convergent evolution of regional endothermy. RESULTS: Across seven species, we identified an average of 10,877 Swiss-Prot annotated genes from an average of 32,474 open reading frames within each species' heart transcriptome. About half of these genes were shared between all species while the remainder included functional differences between our groups of interest (elasmobranch vs. teleost and endotherms vs. ectotherms) as revealed by Gene Ontology (GO) and selection analyses. A repeatedly represented functional category, in both the uniquely expressed elasmobranch genes (total of 259) and the elasmobranch GO enrichment results, involved antibody-mediated immunity, either in the recruitment of immune cells (Fc receptors) or in antigen presentation, including such terms as "antigen processing and presentation of exogenous peptide antigen via MHC class II", and such genes as MHC class II, HLA-DPB1. Molecular adaptation analyses identified three genes in elasmobranchs with a history of positive selection, including legumain (LGMN), a gene with roles in both innate and adaptive immunity including producing antigens for presentation by MHC class II. Comparisons between the endothermic and ectothermic species revealed an enrichment of GO terms associated with cardiac muscle contraction in endotherms, with 19 genes expressed solely in endotherms, several of which have significant roles in lipid and fat metabolism. CONCLUSIONS: This collective comparative evidence provides the first multi-taxa transcriptomic-based perspective on differences between elasmobranchs and teleosts, and suggests various unique features associated with the adaptive immune system of elasmobranchs, pointing in particular to the potential importance of MHC Class II. This in turn suggests that expanded comparative work involving additional tissues, as well as genome sequencing of multiple elasmobranch species would be productive in elucidating the regulatory and genome architectural hallmarks of elasmobranchs.

A unique open reading frame within the comX gene of Streptococcus mutans regulates genetic competence and oxidative stress tolerance.

Streptococcus mutans displays complex regulation of genetic competence, with ComX controlling late competence gene transcription. The rcrRPQ operon has been shown to link oxidative stress tolerance, (p)ppGpp metabolism and competence in S. mutans. Importantly, an rcrR polar (ΔrcrR-P) mutant is hyper-transformable, but an rcrR non-polar (ΔrcrR-NP) mutant cannot be transformed. Transcriptome comparisons of the rcrR mutants using RNA-Seq and quantitative real-time polymerase chain reaction revealed little expression in the 5' region of comX in ΔrcrR-NP, but high level expression in the 3' region. Northern blotting with comX probes revealed two distinct transcripts in the ΔrcrR-P and ΔrcrR-NP strains, and 5' Rapid Amplification of cDNA Ends mapped the 5' terminus of the shorter transcript to nt +140 of the comX structural gene, where a unique 69-aa open reading frame, termed XrpA, was encoded in a different reading frame than ComX. Two single-nucleotide substitution mutants (comX::T162C; comX::T210A) were introduced to disrupt XrpA without affecting the sequence of ComX. When the mutations were in the ΔrcrR-NP genetic background, ComX production and transformation were restored. Overexpression of xrpA led to impaired growth in aerobic conditions and decreased transformability. These results reveal an unprecedented mechanism for competence regulation and stress tolerance by a gene product encoded within the comX gene that appears unique to S. mutans.

Genome Content and Phylogenomics Reveal both Ancestral and Lateral Evolutionary Pathways in Plant-Pathogenic Streptomyces Species.

Streptomyces spp. are highly differentiated actinomycetes with large, linear chromosomes that encode an arsenal of biologically active molecules and catabolic enzymes. Members of this genus are well equipped for life in nutrient-limited environments and are common soil saprophytes. Out of the hundreds of species in the genus Streptomyces, a small group has evolved the ability to infect plants. The recent availability of Streptomyces genome sequences, including four genomes of pathogenic species, provided an opportunity to characterize the gene content specific to these pathogens and to study phylogenetic relationships among them. Genome sequencing, comparative genomics, and phylogenetic analysis enabled us to discriminate pathogenic from saprophytic Streptomyces strains; moreover, we calculated that the pathogen-specific genome contains 4,662 orthologs. Phylogenetic reconstruction suggested that Streptomyces scabies and S. ipomoeae share an ancestor but that their biosynthetic clusters encoding the required virulence factor thaxtomin have diverged. In contrast, S. turgidiscabies and S. acidiscabies, two relatively unrelated pathogens, possess highly similar thaxtomin biosynthesis clusters, which suggests that the acquisition of these genes was through lateral gene transfer.

Diverse sampling of East African haemosporidians reveals chiropteran origin of malaria parasites in primates and rodents.

Phylogenies of parasites provide hypotheses on the history of their movements between hosts, leading to important insights regarding the processes of host switching that underlie modern-day epidemics. Haemosporidian (malaria) parasites lack a well resolved phylogeny, which has impeded the study of evolutionary processes associated with host-switching in this group. Here we present a novel phylogenetic hypothesis that suggests bats served as the ancestral hosts of malaria parasites in primates and rodents. Expanding upon current taxon sampling of Afrotropical bat and bird parasites, we find strong support for all major nodes in the haemosporidian tree using both Bayesian and maximum likelihood approaches. Our analyses support a single transition of haemosporidian parasites from saurian to chiropteran hosts, and do not support a monophyletic relationship between Plasmodium parasites of birds and mammals. We find, for the first time, that Hepatocystis and Plasmodium parasites of mammals represent reciprocally monophyletic evolutionary lineages. These results highlight the importance of broad taxonomic sampling when analyzing phylogenetic relationships, and have important implications for our understanding of key host switching events in the history of malaria parasite evolution.

Laser capture microdissection microscopy and genome sequencing of the avian malaria parasite, Plasmodium relictum.

Acquiring genomic material from avian malaria parasites for genome sequencing has proven problematic due to the nucleation of avian erythrocytes, which produces a large ratio of host to parasite DNA (∼1 million to 1 bp). We tested the ability of laser capture microdissection microscopy to isolate parasite cells from individual avian erythrocytes for four avian Plasmodium species, and subsequently applied whole genome amplification and Illumina sequencing methods to Plasmodium relictum (lineage pSGS1) to produce sequence reads of the P. relictum genome. We assembled ∼335 kbp of parasite DNA from this species, but were unable to completely avoid contamination by host DNA and other sources. However, it is clear that laser capture microdissection holds promise for the isolation of genomic material from haemosporidian parasites in intracellular life stages. In particular, laser capture microdissection may prove useful for isolating individual parasite species from co-infected hosts. Although not explicitly tested in this study, laser capture microdissection may also have important applications for isolation of rare parasite lineages and museum specimens for which no fresh material exists.

Evolutionary and population genomics of the cavity causing bacteria Streptococcus mutans.

Streptococcus mutans is widely recognized as one of the key etiological agents of human dental caries. Despite its role in this important disease, our present knowledge of gene content variability across the species and its relationship to adaptation is minimal. Estimates of its demographic history are not available. In this study, we generated genome sequences of 57 S. mutans isolates, as well as representative strains of the most closely related species to S. mutans (S. ratti, S. macaccae, and S. criceti), to identify the overall structure and potential adaptive features of the dispensable and core components of the genome. We also performed population genetic analyses on the core genome of the species aimed at understanding the demographic history, and impact of selection shaping its genetic variation. The maximum gene content divergence among strains was approximately 23%, with the majority of strains diverging by 5-15%. The core genome consisted of 1,490 genes and the pan-genome approximately 3,296. Maximum likelihood analysis of the synonymous site frequency spectrum (SFS) suggested that the S. mutans population started expanding exponentially approximately 10,000 years ago (95% confidence interval [CI]: 3,268-14,344 years ago), coincidental with the onset of human agriculture. Analysis of the replacement SFS indicated that a majority of these substitutions are under strong negative selection, and the remainder evolved neutrally. A set of 14 genes was identified as being under positive selection, most of which were involved in either sugar metabolism or acid tolerance. Analysis of the core genome suggested that among 73 genes present in all isolates of S. mutans but absent in other species of the mutans taxonomic group, the majority can be associated with metabolic processes that could have contributed to the successful adaptation of S. mutans to its new niche, the human mouth, and with the dietary changes that accompanied the origin of agriculture.

Positive selection, genetic recombination, and intra-host evolution in novel equine coronavirus genomes and other members of the Embecovirus subgenus.

There are several examples of coronaviruses in the Betacoronavirus subgenus Embecovirus that have jumped from an animal to the human host. Studying how evolutionary factors shape coronaviruses in non-human hosts may provide insight into the coronavirus host-switching potential. Equids, such as horses and donkeys, are susceptible to equine coronaviruses (ECoVs). With increased testing prevalence, several ECoV genome sequences have become available for molecular evolutionary analyses, especially those from the United States of America (USA). To date, no analyses have been performed to characterize evolution within coding regions of the ECoV genome. Here, we obtain and describe four new ECoV genome sequences from infected equines from across the USA presenting clinical symptoms of ECoV, and infer ECoV-specific and Embecovirus-wide patterns of molecular evolution. Within two of the four data sets analyzed, we find evidence of intra-host evolution within the nucleocapsid (N) gene, suggestive of quasispecies development. We also identify 12 putative genetic recombination events within the ECoV genome, 11 of which fall in ORF1ab. Finally, we infer and compare sites subject to positive selection on the ancestral branch of each major Embecovirus member clade. Specifically, for the two currently identified human coronavirus (HCoV) embecoviruses that have spilled from animals to humans (HCoV-OC43 and HCoV-HKU1), we find that there are 42 and 2 such sites, respectively, perhaps reflective of the more complex ancestral evolutionary history of HCoV-OC43, which involves several different animal hosts.IMPORTANCEThe Betacoronavirus subgenus Embecovirus contains coronaviruses that not only pose a health threat to animals and humans, but also have jumped from animal to human host. Equids, such as horses and donkeys are susceptible to equine coronavirus (ECoV) infections. No studies have systematically examined evolutionary patterns within ECoV genomes. Our study addresses this gap and provides insight into intra-host ECoV evolution from infected horses. Further, we identify and report natural selection pattern differences between two embecoviruses that have jumped from animals to humans [human coronavirus OC43 and HKU1 (HCoV-OC43 and HCoV-HKU1, respectively)], and hypothesize that the differences observed may be due to the different animal host(s) that each virus circulated in prior to its jump into humans. Finally, we contribute four novel, high-quality ECoV genomes to the scientific community.

Core-gene-encoded peptide regulating virulence-associated traits in Streptococcus mutans.

Recently, high-coverage genome sequence of 57 isolates of Streptococcus mutans, the primary etiological agent of human dental caries, was completed. The SMU.1147 gene, encoding a 61-amino-acid (61-aa) peptide, was present in all sequenced strains of S. mutans but absent in all bacteria in current databases. Reverse transcription-PCR revealed that SMU.1147 is cotranscribed with scnK and scnR, which encode the histidine kinase and response regulator, respectively, of a two-component system (TCS). The C terminus of the SMU.1147 gene product was tagged with a FLAG epitope and shown to be expressed in S. mutans by Western blotting with an anti-FLAG antibody. A nonpolar mutant of SMU.1147 formed less biofilm in glucose-containing medium and grew slower than did the wild-type strain under aerobic and anaerobic conditions, at low pH, or in the presence of H2O2. Mutation of SMU.1147 dramatically reduced genetic competence and expression of comX and comY, compared to S. mutans UA159. The competence defect of the SMU.1147 mutant could not be overcome by addition of sigX-inducing peptide (XIP) in defined medium or by competence-stimulating peptide (CSP) in complex medium. Complementation with SMU.1147 on a plasmid restored all phenotypes. Interestingly, mutants lacking either one of the TCS components and a mutant lacking all three genes behaved like the wild-type strain for all phenotypes mentioned above, but all mutant strains grew slower than UA159 in medium supplemented with 0.3 M NaCl. Thus, the SMU.1147-encoded peptide affects virulence-related traits and dominantly controls quorum-sensing pathways for development of genetic competence in S. mutans.

Characterization of the heart transcriptome of the white shark (Carcharodon carcharias).

BACKGROUND: The white shark (Carcharodon carcharias) is a globally distributed, apex predator possessing physical, physiological, and behavioral traits that have garnered it significant public attention. In addition to interest in the genetic basis of its form and function, as a representative of the oldest extant jawed vertebrate lineage, white sharks are also of conservation concern due to their small population size and threat from overfishing. Despite this, surprisingly little is known about the biology of white sharks, and genomic resources are unavailable. To address this deficit, we combined Roche-454 and Illumina sequencing technologies to characterize the first transciptome of any tissue for this species. RESULTS: From white shark heart cDNA we generated 665,399 Roche 454 reads (median length 387-bp) that were assembled into 141,626 contigs (mean length 503-bp). We also generated 78,566,588 Illumina reads, which we aligned to the 454 contigs producing 105,014 454/Illumina consensus sequences. To these, we added 3,432 non-singleton 454 contigs. By comparing these sequences to the UniProtKB/Swiss-Prot database we were able to annotate 21,019 translated open reading frames (ORFs) of ≥ 20 amino acids. Of these, 19,277 were additionally assigned Gene Ontology (GO) functional annotations. While acknowledging the limitations of our single tissue transcriptome, Fisher tests showed the white shark transcriptome to be significantly enriched for numerous metabolic GO terms compared to the zebra fish and human transcriptomes, with white shark showing more similarity to human than to zebra fish (i.e. fewer terms were significantly different). We also compared the transcriptome to other available elasmobranch sequences, for signatures of positive selection and identified several genes of putative adaptive significance on the white shark lineage. The white shark transcriptome also contained 8,404 microsatellites (dinucleotide, trinucleotide, or tetranucleotide motifs ≥ five perfect repeats). Detailed characterization of these microsatellites showed that ORFs with trinucleotide repeats, were significantly enriched for transcription regulatory roles and that trinucleotide frequency within ORFs was lower than for a wide range of taxonomic groups including other vertebrates. CONCLUSION: The white shark heart transcriptome represents a valuable resource for future elasmobranch functional and comparative genomic studies, as well as for population and other biological studies vital for effective conservation of this globally vulnerable species.

Transcriptomic and genomic evidence for Streptococcus agalactiae adaptation to the bovine environment.

BACKGROUND: Streptococcus agalactiae is a major cause of bovine mastitis, which is the dominant health disorder affecting milk production within the dairy industry and is responsible for substantial financial losses to the industry worldwide. However, there is considerable evidence for host adaptation (ecotypes) within S. agalactiae, with both bovine and human sourced isolates showing a high degree of distinctiveness, suggesting differing ability to cause mastitis. Here, we (i) generate RNAseq data from three S. agalactiae isolates (two putative bovine adapted and one human) and (ii) compare publicly available whole genome shotgun sequence data from an additional 202 isolates, obtained from six host species, to elucidate possible genetic factors/adaptations likely important for S. agalactiae growth and survival in the bovine mammary gland. RESULTS: Tests for differential expression showed distinct expression profiles for the three isolates when grown in bovine milk. A key finding for the two putatively bovine adapted isolates was the up regulation of a lactose metabolism operon (Lac.2) that was strongly correlated with the bovine environment (all 36 bovine sourced isolates on GenBank possessed the operon, in contrast to only 8/151 human sourced isolates). Multi locus sequence typing of all genome sequences and phylogenetic analysis using conserved operon genes from 44 S. agalactiae isolates and 16 additional Streptococcus species provided strong evidence for acquisition of the operon via multiple lateral gene transfer events, with all Streptococcus species known to be major causes of mastitis, identified as possible donors. Furthermore, lactose fermentation tests were only positive for isolates possessing Lac.2. Combined, these findings suggest that lactose metabolism is likely an important adaptation to the bovine environment. Additional up regulation in the bovine adapted isolates included genes involved in copper homeostasis, metabolism of purine, pyrimidine, glycerol and glucose, and possibly aminoglycoside antibiotic resistance. CONCLUSION: We detected several genetic factors likely important in S. agalactiae's adaptation to the bovine environment, in particular lactose metabolism. Of concern is the up regulation of a putative antibiotic resistance gene (GCN5-related N-acetyltransferase) that might reflect an adaptation to the use of aminoglycoside antibiotics within this environment.

Replacing and additive horizontal gene transfer in Streptococcus.

The prominent role of Horizontal Gene Transfer (HGT) in the evolution of bacteria is now well documented, but few studies have differentiated between evolutionary events that predominantly cause genes in one lineage to be replaced by homologs from another lineage ("replacing HGT") and events that result in the addition of substantial new genomic material ("additive HGT"). Here in, we make use of the distinct phylogenetic signatures of replacing and additive HGTs in a genome-wide study of the important human pathogen Streptococcus pyogenes (SPY) and its close relatives S. dysgalactiae subspecies equisimilis (SDE) and S. dysgalactiae subspecies dysgalactiae (SDD). Using recently developed statistical models and computational methods, we find evidence for abundant gene flow of both kinds within each of the SPY and SDE clades and of reduced levels of exchange between SPY and SDD. In addition, our analysis strongly supports a pronounced asymmetry in SPY-SDE gene flow, favoring the SPY-to-SDE direction. This finding is of particular interest in light of the recent increase in virulence of pathogenic SDE. We find much stronger evidence for SPY-SDE gene flow among replacing than among additive transfers, suggesting a primary influence from homologous recombination between co-occurring SPY and SDE cells in human hosts. Putative virulence genes are correlated with transfer events, but this correlation is found to be driven by additive, not replacing, HGTs. The genes affected by additive HGTs are enriched for functions having to do with transposition, recombination, and DNA integration, consistent with previous findings, whereas replacing HGTs seen to influence a more diverse set of genes. Additive transfers are also found to be associated with evidence of positive selection. These findings shed new light on the manner in which HGT has shaped pathogenic bacterial genomes.

Genomes of endangered great hammerhead and shortfin mako sharks reveal historic population declines and high levels of inbreeding in great hammerhead.

Despite increasing threats of extinction to Elasmobranchii (sharks and rays), whole genome-based conservation insights are lacking. Here, we present chromosome-level genome assemblies for the Critically Endangered great hammerhead (Sphyrna mokarran) and the Endangered shortfin mako (Isurus oxyrinchus) sharks, with genetic diversity and historical demographic comparisons to other shark species. The great hammerhead exhibited low genetic variation, with 8.7% of the 2.77 Gbp genome in runs of homozygosity (ROH) > 1 Mbp and 74.4% in ROH >100 kbp. The 4.98 Gbp shortfin mako genome had considerably greater diversity and <1% in ROH > 1 Mbp. Both these sharks experienced precipitous declines in effective population size (Ne) over the last 250 thousand years. While shortfin mako exhibited a large historical Ne that may have enabled the retention of higher genetic variation, the genomic data suggest a possibly more concerning picture for the great hammerhead, and a need for evaluation with additional individuals.

Natural selection differences detected in key protein domains between non-pathogenic and pathogenic Feline Coronavirus phenotypes.

Feline Coronaviruses (FCoVs) commonly cause mild enteric infections in felines worldwide (termed Feline Enteric Coronavirus [FECV]), with around 12% developing into deadly Feline Infectious Peritonitis (FIP; Feline Infectious Peritonitis Virus [FIPV]). Genomic differences between FECV and FIPV have been reported, yet the putative genotypic basis of the highly pathogenic phenotype remains unclear. Here, we used state-of-the-art molecular evolutionary genetic statistical techniques to identify and compare differences in natural selection pressure between FECV and FIPV sequences, as well as to identify FIPV and FECV specific signals of positive selection. We analyzed full length FCoV protein coding genes thought to contain mutations associated with FIPV (Spike, ORF3abc, and ORF7ab). We identified two sites exhibiting differences in natural selection pressure between FECV and FIPV: one within the S1/S2 furin cleavage site, and the other within the fusion domain of Spike. We also found 15 sites subject to positive selection associated with FIPV within Spike, 11 of which have not previously been suggested as possibly relevant to FIP development. These sites fall within Spike protein subdomains that participate in host cell receptor interaction, immune evasion, tropism shifts, host cellular entry, and viral escape. There were 14 sites (12 novel) within Spike under positive selection associated with the FECV phenotype, almost exclusively within the S1/S2 furin cleavage site and adjacent C domain, along with a signal of relaxed selection in FIPV relative to FECV, suggesting that furin cleavage functionality may not be needed for FIPV. Positive selection inferred in ORF7b was associated with the FECV phenotype, and included 24 positively selected sites, while ORF7b had signals of relaxed selection in FIPV. We found evidence of positive selection in ORF3c in FCoV wide analyses, but no specific association with the FIPV or FECV phenotype. We hypothesize that some combination of mutations in FECV may contribute to FIP development, and that is unlikely to be one singular "switch" mutational event. This work expands our understanding of the complexities of FIP development and provides insights into how evolutionary forces may alter pathogenesis in coronavirus genomes.