Anthropogenic contamination sources drive differences in antimicrobial-resistant Escherichia coli in three urban lakes

Antimicrobial resistance (AMR) is an ever-present threat to the treatment of infectious diseases. However, the potential relevance of this phenomenon in environmental reservoirs still raises many questions. Detection of antimicrobial-resistant bacteria in the environment is a critical aspect for understanding the prevalence of resistance outside of clinical settings, as detection in the environment indicates that resistance is likely already widespread. We isolated antimicrobial-resistant Escherichia coli from three urban waterbodies over a 15-month time series, determined their antimicrobial susceptibilities, investigated their population structure, and identified genetic determinants of resistance. We found that E. coli populations at each site were composed of different dominant phylotypes and showed distinct patterns of antimicrobial and multidrug resistance, despite close geographic proximity. Many strains that were genome-sequenced belonged to sequence types of international concern, particularly the ST131 clonal complex. We found widespread resistance to clinically important antimicrobials such as amoxicillin, cefotaxime, and ciprofloxacin, but found that all strains were susceptible to amikacin and the last-line antimicrobials meropenem and fosfomycin. Resistance was most often due to acquirable antimicrobial resistance genes, while chromosomal mutations in gyrA, parC, and parE conferred resistance to quinolones. Whole-genome analysis of a subset of strains further revealed the diversity of the population of E. coli present, with a wide array of AMR and virulence genes identified, many of which were present on the chromosome, including blaCTX-M. Finally, we determined that environmental persistence, transmission between sites, most likely mediated by wild birds, and transfer of mobile genetic elements likely contributed significantly to the patterns observed.IMPORTANCEA One Health perspective is crucial to understand the extent of antimicrobial resistance (AMR) globally, and investigation of AMR in the environment has been increasing in recent years. However, most studies have focused on waterways that are directly polluted by sewage, industrial manufacturing, or agricultural activities. Therefore, there remains a lack of knowledge about more natural, less overtly impacted environments. Through phenotypic and genotypic investigation of AMR in Escherichia coli, this study adds to our understanding of the extent and patterns of resistance in these types of environments, including over a time series, and showed that complex biotic and abiotic factors contribute to the patterns observed. Our study further emphasizes the importance of incorporating the surveillance of microbes in freshwater environments in order to better comprehend potential risks for both human and animal health and how the environment may serve as a sentinel for potential future clinical infections.

Recurring Trans-Atlantic Incursion of Clade 2.3.4.4b H5N1 Viruses by Long Distance Migratory Birds from Northern Europe to Canada in 2022/2023

In December 2022 and January 2023, we isolated clade 2.3.4.4b H5N1 high-pathogenicity avian influenza (HPAI) viruses from six American crows (Corvus brachyrhynchos) from Prince Edward Island and a red fox (Vulpes vulpes) from Newfoundland, Canada. Using full-genome sequencing and phylogenetic analysis, these viruses were found to fall into two distinct phylogenetic clusters: one group containing H5N1 viruses that had been circulating in North and South America since late 2021, and the other one containing European H5N1 viruses reported in late 2022. The transatlantic re-introduction for the second time by pelagic/Icelandic bird migration via the same route used during the 2021 incursion of Eurasian origin H5N1 viruses into North America demonstrates that migratory birds continue to be the driving force for transcontinental dissemination of the virus. This new detection further demonstrates the continual long-term threat of H5N1 viruses for poultry and mammals and the subsequent impact on various wild bird populations wherever these viruses emerge. The continual emergence of clade 2.3.4.4b H5Nx viruses requires vigilant surveillance in wild birds, particularly in areas of the Americas, which lie within the migratory corridors for long-distance migratory birds originating from Europe and Asia. Although H5Nx viruses have been detected at higher rates in North America since 2021, a bidirectional flow of H5Nx genes of American origin viruses to Europe has never been reported. In the future, coordinated and systematic surveillance programs for HPAI viruses need to be launched between European and North American agencies.

Extracellular Polysaccharide Receptor and Receptor-Binding Proteins of the Rhodobacter capsulatus Bacteriophage-like Gene Transfer Agent RcGTA

A variety of prokaryotes produce a bacteriophage-like gene transfer agent (GTA), and the alphaproteobacterial Rhodobacter capsulatus RcGTA is a model GTA. Some environmental isolates of R. capsulatus lack the ability to acquire genes transferred by the RcGTA (recipient capability). In this work, we investigated the reason why R. capsulatus strain 37b4 lacks recipient capability. The RcGTA head spike fiber and tail fiber proteins have been proposed to bind extracellular oligosaccharide receptors, and strain 37b4 lacks a capsular polysaccharide (CPS). The reason why strain 37b4 lacks a CPS was unknown, as was whether the provision of a CPS to 37b4 would result in recipient capability. To address these questions, we sequenced and annotated the strain 37b4 genome and used BLAST interrogations of this genome sequence to search for homologs of genes known to be needed for R. capsulatus recipient capability. We also created a cosmid-borne genome library from a wild-type strain, mobilized the library into 37b4, and used the cosmid-complemented strain 37b4 to identify genes needed for a gain of function, allowing for the acquisition of RcGTA-borne genes. The relative presence of CPS around a wild-type strain, 37b4, and cosmid-complemented 37b4 cells was visualized using light microscopy of stained cells. Fluorescently tagged head spike fiber and tail fiber proteins of the RcGTA particle were created and used to measure the relative binding to wild-type and 37b4 cells. We found that strain 37b4 lacks recipient capability because of an inability to bind RcGTA; the reason it is incapable of binding is that it lacks CPS, and the absence of CPS is due to the absence of genes previously shown to be needed for CPS production in another strain. In addition to the head spike fiber, we found that the tail fiber protein also binds to the CPS.

Patterns of abundance, chromosomal localization, and domain organization among c-di-GMP-metabolizing genes revealed by comparative genomics of five alphaproteobacterial orders

BACKGROUND

Bis-(3'-5')-cyclic dimeric guanosine monophosphate (c-di-GMP) is a bacterial second messenger that affects diverse processes in different bacteria, including the cell cycle, motility, and biofilm formation. Its cellular levels are controlled by the opposing activities of two types of enzymes, with synthesis by diguanylate cyclases containing a GGDEF domain and degradation by phosphodiesterases containing either an HD-GYP or an EAL domain. These enzymes are ubiquitous in bacteria with up to 50 encoded in some genomes, the specific functions of which are mostly unknown.

RESULTS

We used comparative analyses to identify genomic patterns among genes encoding proteins with GGDEF, EAL, and HD-GYP domains in five orders of the class Alphaproteobacteria. GGDEF-containing sequences and GGDEF-EAL hybrids were the most abundant and had the highest diversity of co-occurring auxiliary domains while EAL and HD-GYP containing sequences were less abundant and less diverse with respect to auxiliary domains. There were striking patterns in the chromosomal localizations of the genes found in two of the orders. The Rhodobacterales' EAL-encoding genes and Rhizobiales' GGDEF-EAL-encoding genes showed opposing patterns of distribution compared to the GGDEF-encoding genes. In the Rhodobacterales, the GGDEF-encoding genes showed a tri-modal distribution with peaks mid-way between the origin (ori) and terminus (ter) of replication and at ter while the EAL-encoding genes peaked near ori. The patterns were more complex in the Rhizobiales, but the GGDEF-encoding genes were biased for localization near ter.

CONCLUSIONS

The observed patterns in the chromosomal localizations of these genes suggest a coupling of synthesis and hydrolysis of c-di-GMP with the cell cycle. Moreover, the higher proportions and diversities of auxiliary domains associated with GGDEF domains and GGDEF-EAL hybrids compared to EAL or HD-GYP domains could indicate that more stimuli affect synthesis compared to hydrolysis of c-di-GMP.

Purification of Functional Gene Transfer Agents Using Two-Step Preparative Monolithic Chromatography

Background

Gene transfer agents (GTAs) are phage-like particles that transfer cellular genomic DNA between cells. A hurdle faced in studying GTA function and interactions with cells is the difficulty in obtaining pure and functional GTAs from cultures.

Materials and Methods

We used a novel two-step method for purification of GTAs from R. capsulatus by monolithic chromatography.

Results

Our efficient and simple process had advantages compared to previous approaches. The purified GTAs retained gene transfer activity and the packaged DNA could be used for further studies.

Conclusions

This method is applicable to GTAs produced by other species and small phages, and could be useful for therapeutic applications.

Shared properties of gene transfer agent and core genes revealed by comparative genomics of Alphaproteobacteria

Gene transfer agents (GTAs) are phage-like particles that transfer pieces of cellular genomic DNA to other cells. Homologues of the Rhodobacter capsulatus GTA (RcGTA) structural genes are widely distributed in the Alphaproteobacteria and particularly well conserved in the order Rhodobacterales. Possible reasons for their widespread conservation are still being discussed. It has been suggested that these alphaproteobacterial elements originate from a prophage that was present in an ancestral bacterium and subsequently evolved into a GTA that is now widely maintained in extant descendant lineages. Here, we analysed genomic properties that might relate to the conservation of these alphaproteobacterial GTAs. This revealed that the chromosomal locations of the GTA gene clusters are biased. They primarily occur on the leading strand of DNA replication, at large distances from long repetitive elements, and thus are in regions of lower plasticity, and in areas of extreme GC skew, which also accumulate core genes. These extreme GC skew regions arise from the preferential use of codons with an excess of G over C, a distinct phenomenon from the elevated GC content that has previously been found to be associated with GTA genes. The observed properties, along with their high level of conservation, show that GTA genes share multiple features with core genes in the examined lineages of the Alphaproteobacteria.

Virus-derived gene transfer agents benefit host cells by providing templates for DNA repair

The potential benefits of producing gene transfer agents (GTAs) have long been speculated. A new study in PLOS Biology shows that DNA transfer by these phage-like elements allows cells with DNA damage to perform DNA repair and survive.

Formal recognition and classification of gene transfer agents as viriforms

Morphological and genetic features strongly suggest that gene transfer agents (GTAs) are caudoviricete-derived entities that have evolved in concert with cellular genomes to such a degree that they should not be considered viruses. Indeed, GTA particles resemble caudoviricete virions but, in contrast to caudoviricetes (or any viruses), GTAs can encapsidate at best only part of their own genomes, are induced solely in small subpopulations of prokaryotic host cells and are transmitted vertically as part of cellular genomes during replication and division. Therefore, the lifecycles of GTAs are analogous to virus-derived entities found in parasitoid wasps, which have recently been recognized as non-virus entities and therefore reclassified as viriforms We evaluated three distinct, independently exapted GTA groups for which the genetic basis for GTA particle production has been established. Based on the evidence, we outline a classification scheme for these viriforms

A new perspective on the evolution and diversity of the genus Amdoparvovirus (family Parvoviridae) through genetic characterization, structural homology modeling, and phylogenetics

Amdoparvoviruses (genus Amdoparvovirus, family Parvoviridae) are primarily viruses of carnivorans, but recent studies have indicated that their host range might also extend to rodents and chiropterans. While their classification is based on the full sequence of the major nonstructural protein (NS1), several studies investigating amdoparvoviral diversity have been focused on partial sequences, leading to difficulties in accurately determining species demarcations and leaving several viruses unclassified. In this study, while reporting the complete genomic sequence of a novel amdoparvovirus identified in an American mink (British Columbia amdoparvovirus, BCAV), we studied the phylogenetic relationships of all amdoparvovirus-related sequences and provide a comprehensive reevaluation of their diversity and evolution. After excluding recombinant sequences, phylogenetic and pairwise sequence identity analyses allowed us to define fourteen different viruses, including the five currently classified species, BCAV, and four additional viruses that fulfill the International Committee on Taxonomy of Viruses criteria to be classified as species. We show that the group of viruses historically known as Aleutian mink disease virus (species Carnivore amdoparvovirus 1) should be considered as a cluster of at least four separate viral species that have been co-circulating in mink farms, facilitating the occurrence of inter-species recombination. Genome organization, splicing donor and acceptor sites, and protein sequence motifs were surprisingly conserved within the genus. The sequence of the major capsid protein virus protein 2 (VP2) was significantly more conserved between and within species compared to NS1, a phenomenon possibly linked to antibody-dependent enhancement (ADE). Homology models suggest a remarkably high degree of conservation of the spikes located near the icosahedral threefold axis of the capsid, comprising the surface region associated with ADE. A surprisingly high number of divergent amino acid positions were found in the luminal threefold and twofold axes of the capsid, regions of hitherto unknown function. We emphasize the importance of complete genome analyses and, given the marked phylogenetic inconsistencies across the genome, advise to obtain the complete coding sequences of divergent strains. Further studies on amdoparvovirus biology and structure as well as epidemiological and virus discovery investigations are required to better characterize the ecology and evolution of this important group of viruses.

Distribution and diversity of dog parvoviruses in wild, free-roaming and domestic canids of Newfoundland and Labrador, Canada

Some parvoviruses of carnivorans can infect multiple host species. Since many canine parvoviruses were only discovered recently, their host-range is still unexplored. We examined the host distribution and diversity of five dog parvoviruses in four canine populations from Newfoundland and Labrador, Canada, and investigated the potential for these viruses to cross the species barriers. Canine bocavirus 2 (CBoV-2) and the minute virus of canines were detected in stool from free-roaming dogs from Labrador (5/48 [10.4%] and 3/48 [6.3%], respectively) and two different CBoV-2 variants were identified. Canine bufavirus was identified in stool from free-roaming dogs (1/48, 2.1%) and foxes (3/80, 3.8%) from Labrador, but two different variants were observed in the two host species. The variant found in foxes was highly divergent from previously identified strains. Two cachavirus 1 variants, genetically similar to those circulating in other Canadian wildlife, were found in spleens from Newfoundland coyotes (3/87, 3.5%). Canine parvovirus type 2 (CPV-2) was found in stool from free-roaming dogs from Labrador (2/48, 4.2%) and in spleens from Newfoundland coyotes (3/87, 3.5%). Comparing CPV-2 sequences from these hosts to those retrieved from local symptomatic domestic dogs revealed the presence of a highly heterogeneous viral population as detected strains belonged to five different clades. The close relationship between CPV-2a strains from a dog and a coyote suggests the occurrence of viral transfer between wild and domestic canids. The identification of highly related strains with a similar molecular signature characteristic of older CPV-2 strains in free-roaming and domestic dogs suggests a probable common ancestry and that older CPV-2 strains, which have not been identified in dogs since the 1990s, persist in this part of Canada. Follow-up studies should evaluate samples from a larger number of animals and host species to extensively investigate the possible occurrence of cross-species transmission for recently discovered parvoviruses.

Comparative transmissibility of SARS-CoV-2 variants Delta and Alpha in New England, USA

The SARS-CoV-2 Delta variant rose to dominance in mid-2021, likely propelled by an estimated 40%-80% increased transmissibility over Alpha. To investigate if this ostensible difference in transmissibility is uniform across populations, we partner with public health programs from all six states in New England in the United States. We compare logistic growth rates during each variant's respective emergence period, finding that Delta emerged 1.37-2.63 times faster than Alpha (range across states). We compute variant-specific effective reproductive numbers, estimating that Delta is 63%-167% more transmissible than Alpha (range across states). Finally, we estimate that Delta infections generate on average 6.2 (95% CI 3.1-10.9) times more viral RNA copies per milliliter than Alpha infections during their respective emergence. Overall, our evidence suggests that Delta's enhanced transmissibility can be attributed to its innate ability to increase infectiousness, but its epidemiological dynamics may vary depending on underlying population attributes and sequencing data availability.

Transmission from vaccinated individuals in a large SARS-CoV-2 Delta variant outbreak

An outbreak of over 1,000 COVID-19 cases in Provincetown, Massachusetts (MA), in July 2021-the first large outbreak mostly in vaccinated individuals in the US-prompted a comprehensive public health response, motivating changes to national masking recommendations and raising questions about infection and transmission among vaccinated individuals. To address these questions, we combined viral genomic and epidemiological data from 467 individuals, including 40% of outbreak-associated cases. The Delta variant accounted for 99% of cases in this dataset; it was introduced from at least 40 sources, but 83% of cases derived from a single source, likely through transmission across multiple settings over a short time rather than a single event. Genomic and epidemiological data supported multiple transmissions of Delta from and between fully vaccinated individuals. However, despite its magnitude, the outbreak had limited onward impact in MA and the US overall, likely due to high vaccination rates and a robust public health response.

Prolactin and prolactin receptor expression in the HPG axis and crop during parental care in both sexes of a biparental bird (Columba livia)

During breeding, multiple circulating hormones, including prolactin, facilitate reproductive transitions in species that exhibit parental care. Prolactin underlies parental behaviors and related physiological changes across many vertebrates, including birds and mammals. While circulating prolactin levels often fluctuate across breeding, less is known about how relevant target tissues vary in their prolactin responsiveness via prolactin receptor (PRLR) expression. Recent studies have also investigated prolactin (PRL) gene expression outside of the pituitary (i.e., extra-pituitary PRL), but how PRL gene expression varies during parental care in non-pituitary tissue (e.g., hypothalamus, gonads) remains largely unknown. Further, it is unclear if and how tissue-specific PRL and PRLR vary between the sexes during biparental care. To address this, we measured PRL and PRLR gene expression in tissues relevant to parental care, the endocrine reproductive hypothalamic-pituitary- gonadal (HPG) axis and the crop (a tissue with a similar function as the mammalian mammary gland), across various reproductive stages in both sexes of a biparental bird, the rock dove (Columba livia). We also assessed how these genes responded to changes in offspring presence by adding chicks mid-incubation, simulating an early hatch when prolactin levels were still moderately low. We found that pituitary PRL expression showed similar increases as plasma prolactin levels, and detected extra-pituitary PRL in the hypothalamus, gonads and crop. Hypothalamic and gonadal PRLR expression also changed as birds began incubation. Crop PRLR expression correlated with plasma prolactin, peaking when chicks hatched. In response to replacing eggs with a novel chick mid-incubation, hypothalamic and gonadal PRL and PRLR gene expression differed significantly compared to mid-incubation controls, even when plasma prolactin levels did not differ. We also found sex differences in PRL and PRLR that suggest gene expression may allow males to compensate for lower levels in prolactin by upregulating PRLR in all tissues. Overall, this study advances our understanding of how tissue-specific changes in responsiveness to parental hormones may differ across key reproductive transitions, in response to offspring cues, and between the sexes.

Evidence of transmission from fully vaccinated individuals in a large outbreak of the SARS-CoV-2 Delta variant in Provincetown, Massachusetts

Multiple summer events, including large indoor gatherings, in Provincetown, Massachusetts (MA), in July 2021 contributed to an outbreak of over one thousand COVID-19 cases among residents and visitors. Most cases were fully vaccinated, many of whom were also symptomatic, prompting a comprehensive public health response, motivating changes to national masking recommendations, and raising questions about infection and transmission among vaccinated individuals. To characterize the outbreak and the viral population underlying it, we combined genomic and epidemiological data from 467 individuals, including 40% of known outbreak-associated cases. The Delta variant accounted for 99% of sequenced outbreak-associated cases. Phylogenetic analysis suggests over 40 sources of Delta in the dataset, with one responsible for a single cluster containing 83% of outbreak-associated genomes. This cluster was likely not the result of extensive spread at a single site, but rather transmission from a common source across multiple settings over a short time. Genomic and epidemiological data combined provide strong support for 25 transmission events from, including many between, fully vaccinated individuals; genomic data alone provides evidence for an additional 64. Together, genomic epidemiology provides a high-resolution picture of the Provincetown outbreak, revealing multiple cases of transmission of Delta from fully vaccinated individuals. However, despite its magnitude, the outbreak was restricted in its onward impact in MA and the US, likely due to high vaccination rates and a robust public health response.

Comparative transmissibility of SARS-CoV-2 variants Delta and Alpha in New England, USA

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) Delta variant quickly rose to dominance in mid-2021, displacing other variants, including Alpha. Studies using data from the United Kingdom and India estimated that Delta was 40-80% more transmissible than Alpha, allowing Delta to become the globally dominant variant. However, it was unclear if the ostensible difference in relative transmissibility was due mostly to innate properties of Delta's infectiousness or differences in the study populations. To investigate, we formed a partnership with SARS-CoV-2 genomic surveillance programs from all six New England US states. By comparing logistic growth rates, we found that Delta emerged 37-163% faster than Alpha in early 2021 (37% Massachusetts, 75% New Hampshire, 95% Maine, 98% Rhode Island, 151% Connecticut, and 163% Vermont). We next computed variant-specific effective reproductive numbers and estimated that Delta was 58-120% more transmissible than Alpha across New England (58% New Hampshire, 68% Massachusetts, 76% Connecticut, 85% Rhode Island, 98% Maine, and 120% Vermont). Finally, using RT-PCR data, we estimated that Delta infections generate on average ∼6 times more viral RNA copies per mL than Alpha infections. Overall, our evidence indicates that Delta's enhanced transmissibility could be attributed to its innate ability to increase infectiousness, but its epidemiological dynamics may vary depending on the underlying immunity and behavior of distinct populations.

Newlavirus, a Novel, Highly Prevalent, and Highly Diverse Protoparvovirus of Foxes (Vulpes spp.)

The genus Protoparvovirus (family Parvoviridae) includes several viruses of carnivores. We describe a novel fox protoparvovirus, which we named Newlavirus as it was discovered in samples from Newfoundland and Labrador, Canada. Analysis of the full non-structural protein (NS1) sequence indicates that this virus is a previously uncharacterized species. Newlavirus showed high prevalence in foxes from both the mainland (Labrador, 54/137, 39.4%) and the island of Newfoundland (22/50, 44%) but was not detected in samples from other carnivores, including coyotes (n = 92), lynx (n = 58), martens (n = 146), mink (n = 47), ermines (n = 17), dogs (n = 48), and ringed (n = 4), harp (n = 6), bearded (n = 6), and harbor (n = 2) seals. Newlavirus was found at similar rates in stool and spleen (24/80, 30% vs. 59/152, 38.8%, p = 0.2) but at lower rates in lymph nodes (2/37, 5.4%, p < 0.01). Sequencing a fragment of approximately 750 nt of the capsid protein gene from 53 samples showed a high frequency of co-infection by more than one strain (33.9%), high genetic diversity with 13 genotypes with low sequence identities (70.5-87.8%), and no geographic segregation of strains. Given the high prevalence, high diversity, and the lack of identification in other species, foxes are likely the natural reservoir of Newlavirus, and further studies should investigate its distribution.

ICTV Virus Taxonomy Profile: Marnaviridae 2021

The family Marnaviridae comprises small non-enveloped viruses with positive-sense RNA genomes of 8.6-9.6 kb. Isolates infect marine single-celled eukaryotes (protists) that come from diverse lineages. Some members are known from metagenomic studies of ocean virioplankton, with additional unclassified viruses described from metagenomic datasets derived from marine and freshwater environments. This is a summary of the International Committee on Taxonomy of Viruses (ICTV) Report on the family Marnaviridae, which is available at ictv.global/report/marnaviridae.

A SARS-CoV-2 Cluster in an Acute Care Hospital

BACKGROUND

Little is known about clusters of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection in acute care hospitals.

OBJECTIVE

To describe the detection, mitigation, and analysis of a large cluster of SARS-CoV-2 infections in an acute care hospital with mature infection control policies.

DESIGN

Descriptive study.

SETTING

Brigham and Women's Hospital, Boston, Massachusetts.

PARTICIPANTS

Patients and staff with cluster-related SARS-CoV-2 infections.

INTERVENTION

Close contacts of infected patients and staff were identified and tested every 3 days, patients on affected units were preemptively isolated and repeatedly tested, affected units were cleaned, room ventilation was measured, and specimens were sent for whole-genome sequencing. A case-control study was done to compare clinical interactions, personal protective equipment use, and breakroom and workroom practices in SARS-CoV-2-positive versus negative staff.

MEASUREMENTS

Description of the cluster, mitigation activities, and risk factor analysis.

RESULTS

Fourteen patients and 38 staff members were included in the cluster per whole-genome sequencing and epidemiologic associations. The index case was a symptomatic patient in whom isolation was discontinued after 2 negative results on nasopharyngeal polymerase chain reaction testing. The patient subsequently infected multiple roommates and staff, who then infected others. Seven of 52 (13%) secondary infections were detected only on second or subsequent tests. Eight of 9 (89%) patients who shared rooms with potentially contagious patients became infected. Potential contributing factors included high viral loads, nebulization, and positive pressure in the index patient's room. Risk factors for transmission to staff included presence during nebulization, caring for patients with dyspnea or cough, lack of eye protection, at least 15 minutes of exposure to case patients, and interactions with SARS-CoV-2-positive staff in clinical areas. Whole-genome sequencing confirmed that 2 staff members were infected despite wearing surgical masks and eye protection.

LIMITATION

Findings may not be generalizable.

CONCLUSION

SARS-CoV-2 clusters can occur in hospitals despite robust infection control policies. Insights from this cluster may inform additional measures to protect patients and staff.

PRIMARY FUNDING SOURCE

None.

Connection Between Chromosomal Location and Function of CtrA Phosphorelay Genes in Alphaproteobacteria

Most bacterial chromosomes are circular, with replication starting at one origin (ori) and proceeding on both replichores toward the terminus (ter). Several studies have shown that the location of genes relative to ori and ter can have profound effects on regulatory networks and physiological processes. The CtrA phosphorelay is a gene regulatory system conserved in most alphaproteobacteria. It was first discovered in Caulobacter crescentus where it controls replication and division into a stalked and a motile cell in coordination with other factors. The locations of the ctrA gene and targets of this response regulator on the chromosome affect their expression through replication-induced DNA hemi-methylation and specific positioning along a CtrA activity gradient in the dividing cell, respectively. Here we asked to what extent the location of CtrA regulatory network genes might be conserved in the alphaproteobacteria. We determined the locations of the CtrA phosphorelay and associated genes in closed genomes with unambiguously identifiable ori from members of five alphaproteobacterial orders. The location of the phosphorelay genes was the least conserved in the Rhodospirillales followed by the Sphingomonadales. In the Rhizobiales a trend toward certain chromosomal positions could be observed. Compared to the other orders, the CtrA phosphorelay genes were conserved closer to ori in the Caulobacterales. In contrast, the genes were highly conserved closer to ter in the Rhodobacterales. Our data suggest selection pressure results in differential positioning of CtrA phosphorelay and associated genes in alphaproteobacteria, particularly in the orders Rhodobacterales, Caulobacterales and Rhizobiales that is worth deeper investigation.

Transmission of SARS-CoV-2 from asymptomatic and presymptomatic individuals in healthcare settings despite medical masks and eye protection

We describe 3 instances of SARS-CoV-2 transmission despite medical masks and eye protection, including transmission despite the source person being masked, transmission despite the exposed person being masked, and transmission despite both parties being masked. Whole genome sequencing confirmed perfect homology between source and exposed persons’ viruses in all cases.

Isolating the Role of Corticosterone in the Hypothalamic-Pituitary-Gonadal Transcriptomic Stress Response

Investigation of the negative impacts of stress on reproduction has largely centered around the effects of the adrenal steroid hormone, corticosterone (CORT), and its influence on a system of tissues vital for reproduction-the hypothalamus of the brain, the pituitary gland, and the gonads (the HPG axis). Research on the action of CORT on the HPG axis has predominated the stress and reproductive biology literature, potentially overshadowing other influential mediators. To gain a more complete understanding of how elevated CORT affects transcriptomic activity of the HPG axis, we experimentally examined its role in male and female rock doves (Columba livia). We exogenously administrated CORT to mimic circulating levels during the stress response, specifically 30 min of restraint stress, an experimental paradigm known to increase circulating CORT in vertebrates. We examined all changes in transcription within each level of the HPG axis as compared to both restraint-stressed birds and vehicle-injected controls. We also investigated the differential transcriptomic response to CORT and restraint-stress in each sex. We report causal and sex-specific effects of CORT on the HPG transcriptomic stress response. Restraint stress caused 1567 genes to uniquely differentially express while elevated circulating CORT was responsible for the differential expression of 304 genes. Only 108 genes in females and 8 in males differentially expressed in subjects that underwent restraint stress and those who were given exogenous CORT. In response to elevated CORT and restraint-stress, both sexes shared the differential expression of 5 genes, KCNJ5, CISH, PTGER3, CEBPD, and ZBTB16, all located in the pituitary. The known functions of these genes suggest potential influence of elevated CORT on immune function and prolactin synthesis. Gene expression unique to each sex indicated that elevated CORT affected more gene transcription in females than males (78 genes versus 3 genes, respectively). To our knowledge, this is the first study to isolate the role of CORT in HPG genomic transcription during a stress response. We present an extensive and openly accessible view of the role corticosterone in the HPG transcriptomic stress response. Because the HPG system is well conserved across vertebrates, these data have the potential to inspire new therapeutic strategies for reproductive dysregulation in multiple vertebrate systems, including our own.

Microbiology in the Field: Construction and Validation of a Portable Incubator for Real-Time Quantification of Coliforms and Other Bacteria

Performing microbiological assays on environmental samples in field settings poses logistical challenges with respect to the availability of suitable equipment or the ability to get samples to the laboratory in a timely fashion. For example, the viability of some bacteria can decrease greatly between sampling and arrival to the laboratory for processing. We developed and constructed rugged, reliable, and cost-effective portable incubators that were used by 10 independent field teams to perform microbiological assays on surface water samples from lakes across Canada. Rigorous testing and validation of our incubators ensured that incubation conditions were consistent within and across all 10 field teams and 2 sampling years. Samples from all sites were processed in duplicate and bacterial counts were highly repeatable within and across sampling teams. Bacterial counts were also found to be statistically equivalent to counts obtained with standard laboratory techniques using a conventional incubator. Using this method, thermotolerant coliforms (TTCs) and Escherichia coli were quantified from 432 lakes, allowing comparison to both historical datasets that relied on TTCs and those following current guidelines that use E. coli counts. We found higher loads at the shoreline than the middle of lakes and different patterns between ecozones. E. coli was not frequently detected, but many lakes exceeded Canadian guideline values for activities such as swimming and some even exceeded the guideline value for secondary recreational activities such as boating. To the best of our knowledge, this is the largest bacteriological water quality assessment of freshwater lakes to date in terms of both spatial scale and the number of lakes sampled. Our incubator design can be easily adapted for a wide variety of researcher goals and represents a robust platform for field studies and other applications, including those in remote or low-resources settings.

Multi-host dispersal of known and novel carnivore amdoparvoviruses

Amdoparvoviruses (family Parvoviridae) are ssDNA viruses that cause an immune complex-mediated wasting syndrome in carnivores. They are multi-host pathogens and cross-species infection is facilitated by the fact that viral entry is mediated by cellular Fc receptors recognizing antibody-coated viruses. We developed a pan-amdoparvovirus PCR and screened tissue samples from 666 wild carnivores (families Felidae, Canidae, and Mustelidae) from Newfoundland or Labrador (Canada) and molecularly characterized the identified strains. Fifty-four out of 666 (8.1%) animals were amdoparvovirus-positive. Infection rate was the highest in American mink (34/47, 72.3%), followed by foxes (Arctic and red foxes, 13/311, 4.2%), lynx (2/58, 3.5%), and American martens (5/156, 3.4%). No virus was detected in samples from 87 coyotes and 17 ermines. Viruses from Newfoundland were classified as Aleutian mink disease virus (AMDV). Mink harvested near AMDV-affected fur farms had higher prevalence (24/24, 100%) than other mink (10/23, 43.5%; P < 0.001) and their viruses were phylogenetically closely related to those from farms, while most viruses from other mink were in other clades. Strains from three foxes and two lynx were highly related to mink strains. This proves that farms disperse AMDV that subsequently spreads among wild mink (maintenance host) and transmits to other spillover carnivore hosts. In Labrador two novel viruses were identified, Labrador amdoparvovirus 1 (LaAV-1) found in foxes (9/261, 3.5%) and martens (5/156, 3.4%), and LaAV-2 found in one fox (0.4%). LaAV-1 fulfills all requirements to be classified as a novel species. LaAV-1 was most similar to viruses of mink and skunks (AMDV and skunk amdoparvovirus (SKAV)) while LaAV-2 was more closely related to other viruses infecting canids. LaAV-1 capsid proteins were almost indistinguishable from those of AMDV in some regions, suggesting that LaAV-1 could be a virus of mustelids that can infect foxes. While intensive farming practices provide occasions for inter-species transmission in farms, niche overlap or predation could explain cross-species transmission in the wild, but competition among sympatric species reduces the chances of direct contacts, making this an infrequent event. Pan-amdoparvovirus detection methods in wide epidemiological investigations can play a crucial role in defining amdoparvoviral ecology and evolution and discovering novel viruses.

Influenza A viruses remain infectious for more than seven months in northern wetlands of North America

In this investigation, we used a combination of field- and laboratory-based approaches to assess if influenza A viruses (IAVs) shed by ducks could remain viable for extended periods in surface water within three wetland complexes of North America. In a field experiment, replicate filtered surface water samples inoculated with duck swabs were tested for IAVs upon collection and again after an overwintering period of approximately 6-7 months. Numerous IAVs were molecularly detected and isolated from these samples, including replicates maintained at wetland field sites in Alaska and Minnesota for 181-229 days. In a parallel laboratory experiment, we attempted to culture IAVs from filtered surface water samples inoculated with duck swabs from Minnesota each month during September 2018-April 2019 and found monthly declines in viral viability. In an experimental challenge study, we found that IAVs maintained in filtered surface water within wetlands of Alaska and Minnesota for 214 and 226 days, respectively, were infectious in a mallard model. Collectively, our results support surface waters of northern wetlands as a biologically important medium in which IAVs may be both transmitted and maintained, potentially serving as an environmental reservoir for infectious IAVs during the overwintering period of migratory birds.

A Two-Component System Acquired by Horizontal Gene Transfer Modulates Gene Transfer and Motility via Cyclic Dimeric GMP

Bis-(3'-5')-cyclic dimeric guanosine monophosphate (c-di-GMP) is an important intracellular signaling molecule that affects diverse physiological processes in bacteria. The intracellular levels of c-di-GMP are controlled by proteins acting as diguanylate cyclase (DGC) and phosphodiesterase (PDE) enzymes that synthesize and degrade c-di-GMP, respectively. In the alphaproteobacterium Rhodobacter capsulatus, flagellar motility and gene exchange via production of the gene transfer agent RcGTA are regulated by c-di-GMP. One of the R. capsulatus proteins involved in this regulation is Rcc00620, which contains an N-terminal two-component system response regulator receiver (REC) domain and C-terminal DGC and PDE domains. We demonstrate that the enzymatic activity of Rcc00620 is regulated through the phosphorylation status of its REC domain, which is controlled by a cognate histidine kinase protein, Rcc00621. In this system, the phosphorylated form of Rcc00620 is active as a PDE enzyme and stimulates gene transfer and motility. In addition, we discovered that the rcc00620 and rcc00621 genes are present in only one lineage within the genus Rhodobacter and were acquired via horizontal gene transfer from a distantly related alphaproteobacterium in the order Sphingomonadales. Therefore, a horizontally acquired regulatory system regulates gene transfer in the recipient organism.

Gene transfer across species boundaries in bryophytes: evidence from major life cycle stages in Homalothecium lutescens and H. sericeum

BACKGROUND AND AIMS

The mosses Homalothecium lutescens and H. sericeum are genetically, morphologically and ecologically differentiated; mixed populations sometimes occur. In sympatric populations, intermediate character states among gametophytes and sporophytes have been observed, suggesting hybridization and introgression in such populations.

METHODS

We determined genotypes using bi-allelic co-dominant single nucleotide polymorphism (SNP) markers, specific to either H. lutescens or H. sericeum, to estimate the degree of genetic mixing in 449 moss samples collected from seven sympatric and five allopatric populations on the island of Öland, south Sweden. The samples represented three generations: haploid maternal gametophytes; diploid sporophytes; and haploid sporelings.

KEY RESULTS

Admixture analyses of SNP genotypes identified a majority as pure H. lutescens or H. sericeum, but 76 samples were identified as mildly admixed (17 %) and 17 samples (3.8 %) as strongly admixed. Admixed samples were represented in all three generations in several populations. Hybridization and introgression were bidirectional.

CONCLUSIONS

Our results demonstrate that admixed genomes are transferred between the generations, so that the populations behave as true hybrid zones. Earlier studies of sympatric bryophyte populations with admixed individuals have not been able to show that admixed alleles are transferred beyond the first generation. The presence of true hybrid zones has strong evolutionary implications because genetic material transferred across species boundaries can be directly exposed to selection in the long-lived haploid generation of the bryophyte life cycle, and contribute to local adaptation, long-term survival and speciation.

Stress-mediated convergence of splicing landscapes in male and female rock doves

BACKGROUND

The process of alternative splicing provides a unique mechanism by which eukaryotes are able to produce numerous protein products from the same gene. Heightened variability in the proteome has been thought to potentiate increased behavioral complexity and response flexibility to environmental stimuli, thus contributing to more refined traits on which natural and sexual selection can act. While it has been long known that various forms of environmental stress can negatively affect sexual behavior and reproduction, we know little of how stress can affect the alternative splicing associated with these events, and less still about how splicing may differ between sexes. Using the model of the rock dove (Columba livia), our team previously uncovered sexual dimorphism in the basal and stress-responsive gene transcription of a biological system necessary for facilitating sexual behavior and reproduction, the hypothalamic-pituitary-gonadal (HPG) axis. In this study, we delve further into understanding the mechanistic underpinnings of how changes in the environment can affect reproduction by testing the alternative splicing response of the HPG axis to an external stressor in both sexes.

RESULTS

This study reveals dramatic baseline differences in HPG alternative splicing between males and females. However, after subjecting subjects to a restraint stress paradigm, we found a significant reduction in these differences between the sexes. In both stress and control treatments, we identified a higher incidence of splicing activity in the pituitary in both sexes as compared to other tissues. Of these splicing events, the core exon event is the most abundant form of splicing and more frequently occurs in the coding regions of the gene. Overall, we observed less splicing activity in the 3'UTR (untranslated region) end of transcripts than the 5'UTR or coding regions.

CONCLUSIONS

Our results provide vital new insight into sex-specific aspects of the stress response on the HPG axis at an unprecedented proximate level. Males and females uniquely respond to stress, yet exhibit splicing patterns suggesting a convergent, optimal splicing landscape for stress response. This information has the potential to inform evolutionary theory as well as the development of highly-specific drug targets for stress-induced reproductive dysfunction.

Ecology and Infection Dynamics of Multi-Host Amdoparvoviral and Protoparvoviral Carnivore Pathogens

Amdoparvovirus and Protoparvovirus are monophyletic viral genera that infect carnivores. We performed surveillance for and sequence analyses of parvoviruses in mustelids in insular British Columbia to investigate parvoviral maintenance and cross-species transmission among wildlife. Overall, 19.1% (49/256) of the tested animals were parvovirus-positive. Aleutian mink disease virus (AMDV) was more prevalent in mink (41.6%, 32/77) than martens (3.1%, 4/130), feline panleukopenia virus (FPV) was more prevalent in otters (27.3%, 6/22) than mink (5.2%, 4/77) or martens (2.3%, 3/130), and canine parvovirus 2 (CPV-2) was found in one mink, one otter, and zero ermines (N = 27). Viruses were endemic and bottleneck events, founder effects, and genetic drift generated regional lineages. We identified two local closely related AMDV lineages, one CPV-2 lineage, and five FPV lineages. Highly similar viruses were identified in different hosts, demonstrating cross-species transmission. The likelihood for cross-species transmission differed among viruses and some species likely represented dead-end spillover hosts. We suggest that there are principal maintenance hosts (otters for FPV, raccoons for CPV-2/FPV, mink for AMDV) that enable viral persistence and serve as sources for other susceptible species. In this multi-host system, viral and host factors affect viral persistence and distribution, shaping parvoviral ecology and evolution, with implications for insular carnivore conservation.

Cyclic di-GMP-Mediated Regulation of Gene Transfer and Motility in Rhodobacter capsulatus

Gene transfer agents (GTAs) are bacteriophage-like particles produced by several bacterial and archaeal lineages that contain small pieces of the producing cells' genomes that can be transferred to other cells in a process similar to transduction. One well-studied GTA is RcGTA, produced by the alphaproteobacterium Rhodobacter capsulatus RcGTA gene expression is regulated by several cellular regulatory systems, including the CckA-ChpT-CtrA phosphorelay. The transcription of multiple other regulator-encoding genes is affected by the response regulator CtrA, including genes encoding putative enzymes involved in the synthesis and hydrolysis of the second messenger bis-(3'-5')-cyclic dimeric GMP (c-di-GMP). To investigate whether c-di-GMP signaling plays a role in RcGTA production, we disrupted the CtrA-affected genes potentially involved in this process. We found that disruption of four of these genes affected RcGTA gene expression and production. We performed site-directed mutagenesis of key catalytic residues in the GGDEF and EAL domains responsible for diguanylate cyclase (DGC) and c-di-GMP phosphodiesterase (PDE) activities and analyzed the functions of the wild-type and mutant proteins. We also measured RcGTA production in R. capsulatus strains where intracellular levels of c-di-GMP were altered by the expression of either a heterologous DGC or a heterologous PDE. This adds c-di-GMP signaling to the collection of cellular regulatory systems controlling gene transfer in this bacterium. Furthermore, the heterologous gene expression and the four gene disruptions had similar effects on R. capsulatus flagellar motility as found for gene transfer, and we conclude that c-di-GMP inhibits both RcGTA production and flagellar motility in R. capsulatusIMPORTANCE Gene transfer agents (GTAs) are virus-like particles that move cellular DNA between cells. In the alphaproteobacterium Rhodobacter capsulatus, GTA production is affected by the activities of multiple cellular regulatory systems, to which we have now added signaling via the second messenger dinucleotide molecule bis-(3'-5')-cyclic dimeric GMP (c-di-GMP). Similar to the CtrA phosphorelay, c-di-GMP also affects R. capsulatus flagellar motility in addition to GTA production, with lower levels of intracellular c-di-GMP favoring increased flagellar motility and gene transfer. These findings further illustrate the interconnection of GTA production with global systems of regulation in R. capsulatus, providing additional support for the notion that the production of GTAs has been maintained in this and related bacteria because it provides a benefit to the producing organisms.

Application of a sequence-based taxonomic classification method to uncultured and unclassified marine single-stranded RNA viruses in the order Picornavirales

Metagenomics has altered our understanding of microbial diversity and ecology. This includes its applications to viruses in marine environments that have demonstrated their enormous diversity. Within these are RNA viruses, many of which share genetic features with members of the order Picornavirales; yet, very few of these have been taxonomically classified. The only recognized family of marine RNA viruses is the Marnaviridae, which was founded based on discovery and characterization of the species Heterosigma akashiwo RNA virus. Two additional genera of marine RNA viruses, Labyrnavirus (one species) and Bacillarnavirus (three species), were subsequently defined within the order Picornavirales but not assigned to a family. We have defined a sequence-based framework for taxonomic classification of twenty marine RNA viruses into the family Marnaviridae. Using RNA-dependent RNA polymerase (RdRp) phylogeny and distance-based analyses, we assigned the genera Labyrnavirus and Bacillarnavirus to the family Marnaviridae and created four additional genera in the family: Locarnavirus (four species), Kusarnavirus (one species), Salisharnavirus (four species) and Sogarnavirus (six species). We used pairwise capsid protein comparisons to delineate species within families, with 75 per cent identity as the species demarcation threshold. The family displays high sequence diversities and Jukes-Cantor distances for both the RdRp and capsid genes, suggesting that the classified viruses are not representative of all of the virus diversity within the family and that there are many more extant taxa. Our proposed taxonomic framework provides a sound classification system for this group of viruses that will have broadly applicable principles for other viral groups. It is based on sequence data alone and provides a robust taxonomic framework to include viruses discovered via metagenomic studies, thereby greatly expanding the realm of viruses subject to taxonomic classification.

The NSERC Canadian Lake Pulse Network: A national assessment of lake health providing science for water management in a changing climate

The distribution and quality of water resources vary dramatically across Canada, and human impacts such as land-use and climate changes are exacerbating uncertainties in water supply and security. At the national level, Canada has no enforceable standards for safe drinking water and no comprehensive water-monitoring program to provide detailed, timely reporting on the state of water resources. To provide Canada's first national assessment of lake health, the NSERC Canadian Lake Pulse Network was launched in 2016 as an academic-government research partnership. LakePulse uses traditional approaches for limnological monitoring as well as state-of-the-art methods in the fields of genomics, emerging contaminants, greenhouse gases, invasive pathogens, paleolimnology, spatial modelling, statistical analysis, and remote sensing. A coordinated sampling program of about 680 lakes together with historical archives and a geomatics analysis of over 80,000 lake watersheds are used to examine the extent to which lakes are being altered now and in the future, and how this impacts aquatic ecosystem services of societal importance. Herein we review the network context, objectives and methods.

Induction of Rhodobacter capsulatus Gene Transfer Agent Gene Expression Is a Bistable Stochastic Process Repressed by an Extracellular Calcium-Binding RTX Protein Homologue

Bacteriophage-like gene transfer agents (GTAs) have been discovered in both of the prokaryotic branches of the three-domain phylogenetic tree of life. The production of a GTA (RcGTA) by the phototrophic alphaproteobacterium Rhodobacter capsulatus is regulated by quorum sensing and a phosphorelay homologous to systems in other species that control essential functions such as the initiation of chromosome replication and cell division. In wild-type strains, RcGTA is produced in <3% of cells in laboratory cultures. Mutants of R. capsulatus that exhibit greatly elevated production of RcGTA were created decades ago by chemical mutagenesis, but the nature and molecular consequences of the mutation were unknown. We show that the number of cells in a population that go on to express RcGTA genes is controlled by a stochastic process, in contrast to a genetic process. We used transposon mutagenesis along with a fluorescent protein reporter system and genome sequence data to identify a gene, rcc00280, that encodes an RTX family calcium-binding protein homologue. The Rc280 protein acts as an extracellular repressor of RcGTA gene expression by decreasing the percentage of cells that induce the production of RcGTA.IMPORTANCE GTAs catalyze horizontal gene transfer (HGT), which is important for genomic evolution because the majority of genes found in bacterial genomes have undergone HGT at some point in their evolution. Therefore, it is important to determine how the production of GTAs is regulated to understand the factors that modulate the frequency of gene transfer and thereby specify the tempo of evolution. This work describes a new type of genetic regulation in which an extracellular calcium-binding protein homologue represses the induction of the Rhodobacter capsulatus GTA, RcGTA.

Small extracellular particles with big potential for horizontal gene transfer: membrane vesicles and gene transfer agents

Bacteria are known to release different types of particles that serve various purposes such as the processing of metabolites, communication, and the transfer of genetic material. One of the most interesting aspects of the production of such particles is the biogenesis and trafficking of complex particles that can carry DNA, RNA, proteins or toxins into the surrounding environment to aid in bacterial survival or lead to gene transfer. Two important bacterial extracellular complexes are membrane vesicles and gene transfer agents. In this review, we will discuss the production, contents and functions of these two types of particles as related to their abilities to facilitate horizontal gene transfer.

Discovery and Characterization of Novel RNA Viruses in Aquatic North American Wild Birds

Wild birds are recognized viral reservoirs but our understanding about avian viral diversity is limited. We describe here three novel RNA viruses that we identified in oropharyngeal/cloacal swabs collected from wild birds. The complete genome of a novel gull metapneumovirus (GuMPV B29) was determined. Phylogenetic analyses indicated that this virus could represent a novel avian metapneumovirus (AMPV) sub-group, intermediate between AMPV-C and the subgroup of the other AMPVs. This virus was detected in an American herring (1/24, 4.2%) and great black-backed (4/26, 15.4%) gulls. A novel gull coronavirus (GuCoV B29) was detected in great black-backed (3/26, 11.5%) and American herring (2/24, 8.3%) gulls. Phylogenetic analyses of GuCoV B29 suggested that this virus could represent a novel species within the genus Gammacoronavirus, close to other recently identified potential novel avian coronaviral species. One GuMPV-GuCoV co-infection was detected. A novel duck calicivirus (DuCV-2 B6) was identified in mallards (2/5, 40%) and American black ducks (7/26, 26.9%). This virus, of which we identified two different types, was fully sequenced and was genetically closest to other caliciviruses identified in Anatidae, but more distant to other caliciviruses from birds in the genus Anas. These discoveries increase our knowledge about avian virus diversity and host distributions.

Genetic diversity of Borrelia garinii from Ixodes uriae collected in seabird colonies of the northwestern Atlantic Ocean

The occurrence of Borrelia garinii in seabird ticks, Ixodes uriae, associated with different species of colonial seabirds has been studied since the early 1990s. Research on the population structure of this bacterium in ticks from seabird colonies in the northeastern Atlantic Ocean has revealed admixture between marine and terrestrial tick populations. We studied B. garinii genetic diversity and population structure in I. uriae collected from seabird colonies in the northwestern Atlantic Ocean, in Newfoundland and Labrador, Canada. We applied a multi-locus sequence typing (MLST) scheme to B. garinii found in ticks from four species of seabirds. The B. garinii strains found in this seabird colony ecosystem were diverse. Some were very similar to strains from Asia and Europe, including some obtained from human clinical samples, while others formed a divergent group specific to this region of the Atlantic Ocean. Our findings highlight the genetic complexity of B. garinii circulating in seabird ticks and their avian hosts but also demonstrate surprisingly close connections between B. garinii in this ecosystem and terrestrial sources in Eurasia. Genetic similarities among B. garinii from seabird ticks and humans indicate the possibility that B. garinii circulating within seabird tick-avian host transmission cycles could directly, or indirectly via connectivity with terrestrial transmission cycles, have consequences for human health.

Identification of a Novel Gammaherpesvirus in Canada lynx (Lynx canadensis)

Gammaherpesviruses (GHVs) infect many animal species and are associated with lymphoproliferative disorders in some. Previously, we identified several novel GHVs in North American felids; however, a GHV had never been identified in Canada lynx (Lynx canadensis). We, therefore, hypothesized the existence of an unidentified GHV in lynx. Using degenerate nested and subsequently virus-specific PCR, we amplified and sequenced 3.4 kb of DNA from a novel GHV in lynx, which we named Lynx canadensis gammaherpesvirus 1 (LcaGHV1). Phylogenetic analysis determined that LcaGHV1 is a distinct GHV species belonging to the genus Percavirus. We then estimated the prevalence of LcaGHV1 in lynx by developing a PCR-based assay and detected LcaGHV1 DNA in 36% (95% CI: 22-53%) of lynx spleen DNA samples from Maine, USA and 17% (95% CI: 8-31%) from Newfoundland, Canada. The LcaGHV1 DNA sequences from Maine and Newfoundland lynx were nearly identical to each other (two nucleotide substitutions in 3.4 kb), suggesting that the unique lynx subspecies present on the island of Newfoundland (Lynx canadensis subsolanus) is infected with virus that very closely resembles virus found in mainland lynx. The potential ecologic and pathologic consequences of this novel virus for Canada lynx populations warrant further study.

New Insight Into Avian Papillomavirus Ecology and Evolution From Characterization of Novel Wild Bird Papillomaviruses

Viruses in the family Papillomaviridae have circular dsDNA genomes of approximately 5.7–8.6 kb that are packaged within non-enveloped, icosahedral capsids. The known papillomavirus (PV) representatives infect vertebrates, and there are currently more than 130 recognized PV species in more than 50 genera. We identified 12 novel avian papillomavirus (APV) types in wild birds that could represent five distinct species and two genera. Viruses were detected in paired oropharyngeal/cloacal swabs collected from six bird species, increasing the number of avian species known to harbor PVs by 40%. A new duck PV (DuPV-3) was found in mallard and American black duck (27.6% estimated prevalence) that was monophyletic with other known DuPVs. A single viral type was identified in Atlantic puffin (PuPV-1, 9.8% estimated prevalence), while a higher genetic diversity was found in other Charadriiformes. Specifically, three types [gull PV-1 (GuPV-1), -2, and -3] were identified in two gull species (estimated prevalence of 17% and 2.6% in American herring and great black-backed gull, respectively), and seven types [kittiwake PV-1 (KiPV-1) through -7] were found in black-legged kittiwake (81.3% estimated prevalence). Significantly higher DuPV-3 circulation was observed in spring compared to fall and in adults compared to juveniles. The studied host species’ tendencies to be in crowded environments likely affect infection rates and their migratory behaviors could explain the high viral diversity, illustrating how host behavior can influence viral ecology and distribution. For DuPV-3, GuPV-1, PuPV-1, and KiPV-2, we obtained the complete genomic sequences, which showed the same organization as other known APVs. Phylogenetic analyses showed evidence for virus–host co-divergence at the host taxonomic levels of family, order, and inter-order, but we also observed that host-specificity constraints are relaxed among highly related hosts as we found cross-species transmission within ducks and within gulls. Furthermore, the phylogeny of viruses infecting the Charadriiformes did not match the host phylogeny and gull viruses formed distinct monophyletic clades with kittiwake viruses, possibly reflecting past host-switching events. Considering the vast PV genotype diversity in other hosts and the large number of bird species, many more APVs likely remain to be discovered.

Packaging of Dinoroseobacter shibae DNA into Gene Transfer Agent Particles Is Not Random

Gene transfer agents (GTAs) are phage-like particles which contain a fragment of genomic DNA of the bacterial or archaeal producer and deliver this to a recipient cell. GTA gene clusters are present in the genomes of almost all marine Rhodobacteraceae (Roseobacters) and might be important contributors to horizontal gene transfer in the world’s oceans. For all organisms studied so far, no obvious evidence of sequence specificity or other nonrandom process responsible for packaging genomic DNA into GTAs has been found. Here, we show that knock-out of an autoinducer synthase gene of Dinoroseobacter shibae resulted in overproduction and release of functional GTA particles (DsGTA). Next-generation sequencing of the 4.2-kb DNA fragments isolated from DsGTAs revealed that packaging was not random. DNA from low-GC conjugative plasmids but not from high-GC chromids was excluded from packaging. Seven chromosomal regions were strongly overrepresented in DNA isolated from DsGTA. These packaging peaks lacked identifiable conserved sequence motifs that might represent recognition sites for the GTA terminase complex. Low-GC regions of the chromosome, including the origin and terminus of replication, were underrepresented in DNA isolated from DsGTAs. DNA methylation reduced packaging frequency while the level of gene expression had no influence. Chromosomal regions found to be over- and underrepresented in DsGTA-DNA were regularly spaced. We propose that a “headful” type of packaging is initiated at the sites of coverage peaks and, after linearization of the chromosomal DNA, proceeds in both directions from the initiation site. GC-content, DNA-modifications, and chromatin structure might influence at which sides GTA packaging can be initiated.

ViDiT–CACTUS: an inexpensive and versatile library preparation and sequence analysis method for virus discovery and other microbiology applications

High-throughput sequencing (HTS) technologies are becoming increasingly important within microbiology research, but aspects of library preparation, such as high cost per sample or strict input requirements, make HTS difficult to implement in some niche applications and for research groups on a budget. To answer these necessities, we developed ViDiT, a customizable, PCR-based, extremely low-cost (less than US$5 per sample), and versatile library preparation method, and CACTUS, an analysis pipeline designed to rely on cloud computing power to generate high-quality data from ViDiT-based experiments without the need of expensive servers. We demonstrate here the versatility and utility of these methods within three fields of microbiology: virus discovery, amplicon-based viral genome sequencing, and microbiome profiling. ViDiT–CACTUS allowed the identification of viral fragments from 25 different viral families from 36 oropharyngeal–cloacal swabs collected from wild birds, the sequencing of three almost complete genomes of avian influenza A viruses (>90% coverage), and the characterization and functional profiling of the complete microbial diversity (bacteria, archaea, viruses) within a deep-sea carnivorous sponge. ViDiT–CACTUS demonstrated its validity in a wide range of microbiology applications, and its simplicity and modularity make it easily implementable in any molecular biology laboratory, towards various research goals.

Sex-biased transcriptomic response of the reproductive axis to stress

Stress is a well-known cause of reproductive dysfunction in many species, including birds, rodents, and humans, though males and females may respond differently. A powerful way to investigate how stress affects reproduction is by examining its effects on a biological system essential for regulating reproduction, the hypothalamic-pituitary-gonadal (HPG) axis. Often this is done by observing how a stressor affects the amount of glucocorticoids, such as cortisol or corticosterone, circulating in the blood and their relationship with a handful of known HPG-producing reproductive hormones, like testosterone and estradiol. Until now, we have lacked a full understanding of how stress affects all genomic activity of the HPG axis and how this might differ between the sexes. We leveraged a highly replicated and sex-balanced experimental approach to test how male and female rock doves (Columba livia) respond to restraint stress at the level of their transcriptome. Females exhibit increased genomic responsiveness to stress at all levels of their HPG axis as compared to males, and these responsive genes are mostly unique to females. Reasons for this may be due to fluctuations in the female endocrine environment over the reproductive cycle and/or their evolutionary history, including parental investment and the potential for maternal effects. Direct links between genome to phenome cause and effect cannot be ascertained at this stage; however, the data we report provide a vital genomic foundation on which sex-specific reproductive dysfunction and adaptation in the face of stress can be further experimentally studied, as well as novel gene targets for genetic intervention and therapy investigations.

Dynamics of Heterotrophic Bacterial Assemblages within Synechococcus Cultures

Interactions between photoautotrophic and heterotrophic microorganisms are central to the marine microbial ecosystem. Lab cultures of one of the dominant marine photoautotrophs, Synechococcus, have historically been difficult to render axenic, presumably because these bacteria depend upon other organisms to grow under these conditions. These tight associations between Synechococcus and heterotrophic bacteria represent a good relevant system to study interspecies interactions. Ten individual Synechococcus strains, isolated from eutrophic and oligotrophic waters, were chosen for investigation. Four to six dominant associated heterotrophic bacteria were detected in the liquid cultures of each Synechococcus isolate, comprising members of the Cytophaga-Flavobacteria-Bacteroides (CFB) group (mainly from Flavobacteriales and Cytophagales), Alphaproteobacteria (mainly from the Roseobacter clade), Gammaproteobacteria (mainly from the Alteromonadales and Pseudomonadales), and Actinobacteria The presence of the CFB group, Gammaproteobacteria, and Actinobacteria showed clear geographic patterns related to the isolation environments of the Synechococcus bacteria. An investigation of the population dynamics within a growing culture (XM-24) of one of the isolates, including an evaluation of the proportions of cells that were free-living versus aggregated/attached, revealed interesting patterns for different bacterial groups. In Synechococcus sp. strain XM-24 culture, flavobacteria, which was the most abundant group throughout the culture period, tended to be aggregated or attached to the Synechococcus cells, whereas the actinobacteria demonstrated a free-living lifestyle, and roseobacters displayed different patterns depending on the culture growth phase. Factors contributing to these succession patterns for the heterotrophs likely include interactions among the culture community members, their relative abilities to utilize different compounds produced by Synechococcus cells and changes in the compounds released as culture growth proceeds, and their responses to other changes in the environmental conditions throughout the culture period.IMPORTANCE Marine microbes exist within an interactive ecological network, and studying their interactions is an important part of understanding their roles in global biogeochemical cycling and the determinants of microbial diversity. In this study, the dynamic relationships between Synechococcus spp. and their associated heterotrophic bacteria were investigated. Synechococcus-associated heterotrophic bacteria had similar geographic distribution patterns as their "host" and displayed different lifestyles (free-living versus attached/aggregated) according to the Synechococcus culture growth phases. Combined organic carbon composition and bacterial lifestyle data indicated a potential for succession in carbon utilization patterns by the dominant associated heterotrophic bacteria. Comprehending the interactions between photoautotrophs and heterotrophs and the patterns of organic carbon excretion and utilization is critical to understanding their roles in oceanic biogeochemical cycling.

CD24 induces changes to the surface receptors of B cell microvesicles with variable effects on their RNA and protein cargo

The CD24 cell surface receptor promotes apoptosis in developing B cells, and we recently found that it induces B cells to release plasma membrane-derived, CD24-bearing microvesicles (MVs). Here we have performed a systematic characterization of B cell MVs released from WEHI-231 B lymphoma cells in response to CD24 stimulation. We found that B cells constitutively release MVs of approximately 120 nm, and that CD24 induces an increase in phosphatidylserine-positive MV release. RNA cargo is predominantly comprised of 5S rRNA, regardless of stimulation; however, CD24 causes a decrease in the incorporation of protein coding transcripts. The MV proteome is enriched with mitochondrial and metabolism-related proteins after CD24 stimulation; however, these changes were variable and could not be fully validated by Western blotting. CD24-bearing MVs carry Siglec-2, CD63, IgM, and, unexpectedly, Ter119, but not Siglec-G or MHC-II despite their presence on the cell surface. CD24 stimulation also induces changes in CD63 and IgM expression on MVs that is not mirrored by the changes in cell surface expression. Overall, the composition of these MVs suggests that they may be involved in releasing mitochondrial components in response to pro-apoptotic stress with changes to the surface receptors potentially altering the cell type(s) that interact with the MVs.

Full genetic characterization and epidemiology of a novel amdoparvovirus in striped skunk (Mephitis mephitis)

Amdoparvovirus is a newly defined parvoviral genus that contains four species (Carnivore amdoparvovirus 1-4), including the well-known Aleutian mink disease virus (AMDV). Amdoparvoviruses cause an immune-associated and often lethal wasting syndrome in Mustelidae and Caninae hosts. In this study, we molecularly investigated amdoparvoviruses detected in 44 striped skunks (Mephitis mephitis) found dead in and around Vancouver, British Columbia, Canada. Some of the animals exhibited pathological changes compatible with amdoparvovirus-associated disease. The nearly complete genomic sequence was obtained for seven different strains and our analyses show how this virus, which we named skunk amdoparvovirus (SKAV), should be classified as a separate species within the genus (proposed Carnivore amdoparvovirus 5). We detected co-infections, recombinant genomes, at least three separate viral lineages, and preliminary evidence for geographic segregation of lineages. Furthermore, we proved that similar viruses, only partially characterized in previous studies and labeled as AMDV, circulate in skunks from other distant areas of North America (Ontario and California) and found evidence for spillover events in mink (Neovison vison). Although SKAVs are capable of causing disease in infected animals, a high proportion of sub-clinical infections has been observed, suggesting these animals might act as asymptomatic carriers and pose a threat to wild and captive carnivores. Finally, we highlight the need for more specific diagnostic tests and further molecular investigations to clarify the epidemiology and host- and geographical distributions of amdoparvoviruses in terrestrial carnivores, especially because the whole spectrum of viral diversity in this group is likely still unknown.