Search

While gammacoronaviruses mainly comprise infectious bronchitis virus (IBV) and its closely related bird coronaviruses (CoVs), the only mammalian gammacoronavirus was discovered from a white beluga whale (beluga whale CoV [BWCoV] SW1) in 2008. In this study, we discovered a novel gammacoronavirus from fecal samples from three Indo-Pacific bottlenose dolphins (Tursiops aduncus), which we named bottlenose dolphin CoV (BdCoV) HKU22. All the three BdCoV HKU22-positive samples were collected on the same date, suggesting a cluster of infection, with viral loads of 1 × 10(3) to 1 × 10(5) copies per ml. Clearance of virus was associated with a specific antibody response against the nucleocapsid of BdCoV HKU22. Complete genome sequencing and comparative genome analysis showed that BdCoV HKU22 and BWCoV SW1 have similar genome characteristics and structures. Their genome size is about 32,000 nucleotides, the largest among all CoVs, as a result of multiple unique open reading frames (NS5a, NS5b, NS5c, NS6, NS7, NS8, NS9, and NS10) between their membrane (M) and nucleocapsid (N) protein genes. Although comparative genome analysis showed that BdCoV HKU22 and BWCoV SW1 should belong to the same species, a major difference was observed in the proteins encoded by their spike (S) genes, which showed only 74.3 to 74.7% amino acid identities. The high ratios of the number of synonymous substitutions per synonymous site (Ks) to the number of nonsynonymous substitutions per nonsynonymous site (Ka) in multiple regions of the genome, especially the S gene (Ka/Ks ratio, 2.5), indicated that BdCoV HKU22 may be evolving rapidly, supporting a recent transmission event to the bottlenose dolphins. We propose a distinct species, Cetacean coronavirus, in Gammacoronavirus, to include BdCoV HKU22 and BWCoV SW1, whereas IBV and its closely related bird CoVs represent another species, Avian coronavirus, in Gammacoronavirus.

Streptococcus suis is an important zoonotic bacterial pathogen posing a threat to the pig industry as well as public health, for which the mechanisms of growth and cell division remain largely unknown. Developing convenient genetic tools that can achieve strictly controlled gene expression is of great value for investigating these fundamental physiological processes of S. suis. In this study, we first identified three strong constitutive promoters, P

Mutations in the viral genome of SARS-CoV-2 can impact the performance of molecular diagnostic assays. In some cases, such as S gene target failure, the impact can serve as a unique indicator of a particular SARS-CoV-2 variant and provide a method for rapid detection. Here we describe partial ORF1ab gene target failure (pOGTF) on the cobas® SARS-CoV-2 assays, defined by a ≥2 thermocycles delay in detection of the ORF1ab gene compared to the E gene. We demonstrate that pOGTF is 97% sensitive and 99% specific for SARS-CoV-2 lineage BA.2.12.1, an emerging variant in the United States with spike L452Q and S704L mutations that may impact transmission, infectivity, and/or immune evasion. Increasing rates of pOGTF closely mirrored rates of BA.2.12.1 sequences uploaded to public databases, and, importantly increasing local rates of pOGTF also mirrored increasing overall test positivity. Use of pOGTF as a proxy for BA.2.12.1 provides faster tracking of the variant than whole-genome sequencing and can benefit laboratories without sequencing capabilities.

The whole genomic sequencing (WGS) of SARS-CoV-2 has been performed extensively and is playing a crucial role in fighting against COVID-19 pandemic. Obtaining sufficient WGS data from clinical samples is often challenging especially from the samples with low viral load. We evaluated two SARS-CoV-2 sequencing protocols for their efficiency/accuracy and limitations. Sequence coverage of >95% was obtained by Swift normalase amplicon SARS-CoV-2 panels (SNAP) protocol for all the samples with Ct ≤ 35 and by COVIDSeq protocol for 97% of samples with Ct ≤ 30. Sample RNA quantitation obtained using digital PCR provided more precise cutoff values. The quantitative digital PCR cutoff values for obtaining 95% coverage are 10.5 copies/μL for SNAP protocol and 147 copies/μL for COVIDSeq protocol. Combining FASTQ files obtained from 2 protocols improved the outcome of sequence analysis by compensating for missing amplicon regions. This process resulted in an increase of sequencing coverage and lineage call precision.

Long noncoding RNAs (lncRNAs) function as microregulatory factors that influence gene expression after a variety of pathogenic infections, and they have been extensively studied in the past few years. Although less attention has been paid to lncRNAs in lower vertebrates than in mammals, current studies reveals that lncRNAs play a vital role in fish stimulated by pathogens. Here, we discovered a new lncRNA, termed MIR2187HG, which can function as a precursor of a small RNA, miR-2187-3p, with regulatory functions in the miiuy croaker (Miichthys miiuy). Upon Siniperca chuatsi rhabdovirus (SCRV) virus infection, the expression levels of MIR2187HG were remarkably enhanced. Elevated MIR2187HG expression can act as a pivotally negative regulator that participates in the innate immune response of teleost fish to inhibit the intracellular TANK-binding kinase 1 (TBK1)-mediated antiviral signaling pathways, which can effectively avoid excessive immunity. In addition, we found that SCRV could also utilize MIR2187HG to enhance its own numbers. Our results not only provide evidence regarding the involvement of the lncRNAs in response to viruses in fish but also broaden our understanding of the function of lncRNAs as precursor microRNAs (miRNAs) in teleost fish for the first time. IMPORTANCE SCRV infection upregulates MIR2187HG levels, which in turn suppresses SCRV-triggered type I interferon production, thus promoting viral replication in miiuy croaker. Notably, MIR2187HG regulates the release of miR-2187-3p, and TBK1 is a target of miR-2187-3p. MIR2187HG could acquire from miR-2187-3p the function of inhibiting TBK1 expression and subsequently modulate TBK1-mediated NF-κB and IRF3 signaling. The collective results suggest that the novel regulation mechanism of TBK1-mediated antiviral response during RNA viral infection was regulated by MIR2187HG. Therefore, a new regulation mechanism for lncRNAs to regulate antiviral immune responses in fish is proposed.

We provide a novel single restriction enzyme (RE; BsaHI) digestion approach for detecting distinct pathotypes of Newcastle disease virus (NDV). After scanning 4,000 F gene nucleotide sequences in the NCBI database, we discovered a single RE (BsaHI) digestion site in the cleavage site. APMV-I “F gene” class II-specific primer-based reverse transcriptase PCR was utilized to amplify a 535-bp fragment, which was then digested with the RE (BsaHI) for pathotyping avian NDV field isolates and pigeon paramyxovirus-1 isolates. The avirulent (lentogenic and mesogenic strains) produced 189- and 346-bp fragments, respectively, but the result in velogenic strains remained undigested with 535-bp fragments. In addition, 45 field NDV isolates and 8 vaccine strains were used to confirm the approach. The sequence-based analysis also agrees with the data obtained utilizing the single RE (BsaHI) digestion approach. The proposed technique has the potential to distinguish between avirulent and virulent strains in a short time span, making it valuable in NDV surveillance and monitoring research. IMPORTANCE The extensive use of the NDV vaccine strain and the existence of avirulent NDV strains in wild birds makes it difficult to diagnose Newcastle Disease virus (NDV). The intracerebral pathogenicity index (ICPI) and/or sequencing-based identification, which are required to determine virulent NDV, are time-consuming, costly, difficult, and cruel techniques. We evaluated 4,000 F gene nucleotide sequences and discovered a restriction enzyme (RE; BsaHI) digestion technique for detecting NDV and vaccine pathotypes in a short time span, which is cost-effective and useful for field cases as well as for large-scale NDV monitoring and surveillance. The data acquired using the single RE BsaHI digestion technique agree with the sequence-based analysis.

Porcine parvovirus genotype 1 (PPV1) causes reproductive disorder in swine and is prevalent in China. Recently, six new genotypes of PPVs (PPV2 through PPV7) have also been detected in Chinese swine herds. However, the coinfection status of all these seven genotypes of PPVs (PPV1-7) in China was not clarified yet. In this study, we developed a panel of PPV1–7 PCR assays with satisfied specificity, sensitivity and reproducibility and then applied to the detection of PPV1–7 in 435 clinical samples collected from eight provinces of China in 2016–2020. A total of 55.40% samples (241 out of 435) were PPV positive, while PPV2 and PPV3 (both 22.53%) belonging to the genus of Tetraparvovirus were the most prevalent genotypes. Noticeably, PPV1–7 strains were more prevalent in nursery and finishing pigs than in suckling pigs. In addition, coinfection could be detected in all eight provinces and 27.36% (119/435) samples were coinfected with two to five genotypes of PPVs. Meanwhile, the coinfection of PPVs with PCV2 was 22.30% (97/435). Twenty complete genomes of representative PPV1–7 were determined, and phylogenetic analysis confirmed the genotyping results by sequence comparisons and PCR assays. Remarkably, the PPV7 HBTZ20180519-152 strain from domestic pig was recombined from parental JX15-like and JX38-like isolates from wild boars. Selective pressure analysis based on VP2 sequences of PPV1–7 showed that they were predominantly under negative selection, while few positive selection sites could be detected in VP2 of PPV7. Overall, this systematic investigation unveils high prevalence and coinfection of PPV1–7 in China from 2016 to 2020. IMPORTANCE Porcine parvoviruses (PPVs) are prevalent in China associating with reproductive failure in swine. The coinfection of seven genotypes of PPVs (PPV1-7) might have synergistic effects on PPV1 associated SMEDI syndrome. However, the coinfection status of PPV1–7 in China is not clear yet. This study showed that PPV1–7 strains are highly prevalent (55.40%) in China and mainly in nursery and finishing pigs in recent years. In addition, the coinfections of different genotypes of PPVs (27.36%) and PPVs with PCV2 (22.30%) are common. Geographic analysis indicated that different genotypes of PPVs are widely cocirculating in China. Intriguingly, a PPV7 strain from the domestic pig was detected as a recombinant from two wild boar isolates. Selective pressure analyses showed that PPV1–7 are mainly under purifying selection. Our findings provide the first systematic investigation on the prevalence, coinfection, and evolution of PPV1 through PPV7 in Chinese swineherds from 2016 to 2020.

Histophilus somni is a Gram-negative bacterial organism that acts as an opportunistic pathogen and is a fastidious member of the Pasteurellaceae family associated with diseases of respiratory, reproductive, cardiac, and other tissues of ruminants. We identified an intervening sequence (IVS) embedded in all five copies of the 23S rRNA gene in the closed genome sequence of the H. somni isolate USDA-ARS-USMARC-63250 that may play an important role in affecting the biology of the organism. Sequencing the RNA from this isolate shows that most of the IVS is cleaved from the transcript, resulting in independent fragments of this structural rRNA that remain functional within the bacterial ribosome. The IVS lies between positions 1170 and 1278 bp of the 3,017-bp gene and exhibits self-complementarity between its 5′ and 3′ ends that predicts a stem-loop structure interrupting helix-45 in the transcribed 23S rRNA. Excision removes a 94-nucleotide (nt) stem-loop structure that displays an unusual 1-nt 3′ end overhang instead of the more typical 2-nt overhang commonly observed at the ends of other excised IVS stem-loops. A comparison with genomes of other H. somni isolates indicates that this IVS is highly conserved, with 31 of 32 complete genomes having similar interruptions of canonical 23S rRNA genes. The potential biological effects of either the released IVS or the fragmentation of the functional 23S rRNA are unknown, but fragmentation may enhance rRNA degradation in ways that contribute to the regulation of gene expression. IMPORTANCE The genome biology underlying H. somni virulence, pathogenicity, environmental adaptability, and broad tissue tropism is understood poorly. We identified a novel H. somni 109-nt IVS stem-loop structure, of which the central portion is excised from the 23S rRNA transcript, resulting in the fragmentation of this rRNA in the H. somni isolate USDA-ARS-USMARC-63250 and the release of a 94-nt structured RNA of unknown function. We determined that this peculiar rRNA biology is widespread among sequenced H. somni isolates, suggesting it has importance to organism biology. The fragmented 23S rRNA molecules remain functional in the ribosome, given that the isolate grows in culture. The structured excised portion of the IVS, presumably due to the action of the endoribonuclease III, has an unusual 1-nt 3′ end overhang. This newly discovered H. somni 23S rRNA fragmentation may enhance rRNA degradation providing a previously unrecognized avenue for regulating H. somni biological processes.

Azotobacter vinelandii is a nitrogen-fixing free-living soil microbe that has been studied for decades in relation to biological nitrogen fixation (BNF). It is highly amenable to genetic manipulation, helping to unravel the intricate importance of different proteins involved in the process of BNF, including the biosynthesis of cofactors that are essential to assembling the complex metal cofactors that catalyze the difficult reaction of nitrogen fixation. Additionally, A. vinelandii accomplishes this feat while growing as an obligate aerobe, differentiating it from many of the nitrogen-fixing bacteria that are associated with plant roots. The ability to function in the presence of oxygen makes A. vinelandii suitable for application in various potential biotechnological schemes. In this study, we employed transposon sequencing (Tn-seq) to measure the fitness defects associated with disruptions of various genes under nitrogen-fixing dependent growth, versus growth with extraneously provided urea as a nitrogen source. The results allowed us to probe the importance of more than 3,800 genes, revealing that many genes previously believed to be important, can be successfully disrupted without impacting cellular fitness. IMPORTANCE These results provide insights into the functional redundancy in A. vinelandii, while also providing a direct measure of fitness for specific genes associated with the process of BNF. These results will serve as a valuable reference tool in future studies to uncover the mechanisms that govern this process.

Sin Nombre orthohantavirus (SNV), a negative-sense, single-stranded RNA virus that is carried and transmitted by the North American deer mouse Peromyscus maniculatus, can cause infection in humans through inhalation of aerosolized excreta from infected rodents. This infection can lead to hantavirus cardiopulmonary syndrome (HCPS), which has a ∼36% case-fatality rate. We used reverse transcriptase-quantitative polymerase chain reaction (RT-qPCR) to confirm SNV infection in a patient and identified SNV in lung tissue in wild-caught rodents from potential sites of exposure. Using viral whole-genome sequencing (WGS), we identified the likely site of transmission and discovered SNV in multiple rodent species not previously known to carry the virus. Here we report, for the first time, the use of SNV WGS to pinpoint a likely site of human infection and identify SNV simultaneously in multiple rodent species in an area of known host-to-human transmission. These results will impact epidemiology and infection control for hantaviruses by tracing zoonotic transmission and investigating possible novel host reservoirs. Importance Orthohantaviruses cause severe disease in humans and can be lethal in up to 40% of cases. Sin Nombre orthohantavirus (SNV) is the main cause of hantavirus disease in North America. In this study, we sequenced SNV from an infected patient and wild-caught rodents to trace location of infection. We also discovered SNV in rodent species not previously known to carry SNV. These studies demonstrate for the first time the use of virus sequencing to trace transmission of SNV, and describe infection in novel rodent species.

Cutinases are enzymes produced by phytopathogenic fungi like Moniliophthora roreri. The three genome-located cutinase genes of M. roreri were amplified from cDNA of fungi growing in different induction culture media for cutinase production. The mrcut1 gene was expressed in the presence of a cacao cuticle, while the mrcut2 and mrcut3 genes were expressed when an apple cuticle was used as the inducer. The sequences of all genes were obtained and analyzed by bioinformatics tools to determine the presence of signal peptides, introns, glycosylation, and regulatory sequences. Also, the theoretical molecular weight and pI were obtained and experimentally confirmed. Finally, cutinase 1 from M. roreri (MRCUT1) was selected for heterologous expression in Escherichia coli. Successful overexpression of MRCUT1 was observed with the highest enzyme activity of 34,036 U/mg under the assay conditions at 40°C and pH 8. Furthermore, the degradation of different synthetic polyesters was evaluated; after 21 days, 59% of polyethylene succinate (PES), 43% of polycaprolactone (PCL), and 31% of polyethylene terephthalate (PET) from plastic residues were degraded. IMPORTANCE Plastic pollution is exponentially increasing; even the G20 has recognized an urgent need to implement actions to reduce it. In recent years, searching for enzymes that can degrade plastics, especially those based on polyesters such as PET, has been increasing as they can be a green alternative to the actual plastic degradation process. A promising option in recent years refers to biological tools such as enzymes involved in stages of partial and even total degradation of some plastics. In this context, the MRCUT1 enzyme can degrade polyesters contained in plastic residues in a short time. Besides, there is limited knowledge about the biochemical properties of cutinases from M. roreri. Commonly, fungal enzymes are expressed as inclusion bodies in E. coli with reduced activity. Interestingly, the successful expression of one cutinase of M. roreri in E. coli with enhanced activity is described.

Shigella infection remains a public health problem in much of the world. Classic models of Shigella pathogenesis suggest that microfold epithelial cells in the small intestine are the preferred initial site of invasion. However, recent evidence supports an alternative model in which Shigella primarily infects a much wider range of epithelial cells that reside primarily in the colon. Here, we investigated whether the luminal pH difference between the small intestine and the colon could provide evidence in support of either model of Shigella flexneri pathogenesis. Because virulence factors culminating in cellular invasion are linked to biofilms in S. flexneri, we examined the effect of pH on the ability of S. flexneri to form and maintain adherent biofilms induced by deoxycholate. We showed that a basic pH (as expected in the small intestine) inhibited formation of biofilms and dispersed preassembled mature biofilms, while an acidic pH (similar to the colonic environment) did not permit either of these effects. To further elucidate this phenomenon at the molecular level, we probed the transcriptomes of biofilms and S. flexneri grown under different pH conditions. We identified specific amino acid (cysteine and arginine) metabolic pathways that were enriched in the bacteria that formed the biofilms but decreased when the pH increased. We then utilized a type III secretion system reporter strain to show that increasing pH reduced deoxycholate-induced virulence of S. flexneri in a dose-dependent manner. Taken together, these experiments support a model in which Shigella infection is favored in the colon because of the local pH differences in these organs.

Viruses of the phylum Nucleocytoviricota, or nucleo-cytoplasmic large DNA viruses (NCLDVs), undergo a cytoplasmic or nucleo-cytoplasmic cycle, the latter of which involves both nuclear and cytoplasmic compartments to proceed viral replication. Medusavirus, a recently isolated NCLDV, has a nucleo-cytoplasmic replication cycle in amoebas during which the host nuclear membrane apparently remains intact, a unique feature among amoeba-infecting NCLDVs. The medusavirus genome lacks most transcription genes but encodes a full set of histone genes. To investigate its infection strategy, we performed a time course RNA sequencing (RNA-seq) experiment. All viral genes were transcribed and classified into five temporal expression clusters. The immediate early genes (cluster 1, 42 genes) were mostly (83%) of unknown functions, frequently (95%) associated with a palindromic promoter-like motif, and often (45%) encoded putative nucleus-localized proteins. These results suggest massive reshaping of the host nuclear environment by viral proteins at an early stage of infection. Genes in other expression clusters (clusters 2 to 5) were assigned to various functional categories. The virally encoded core histone genes were in cluster 3, whereas the viral linker histone H1 gene was in cluster 1, suggesting they have distinct roles during the course of the virus infection. The transcriptional profile of the host Acanthamoeba castellanii genes was greatly altered postinfection. Several encystment-related host genes showed increased representation levels at 48 h postinfection, which is consistent with the previously reported amoeba encystment upon medusavirus infection. IMPORTANCE Medusavirus is an amoeba-infecting giant virus that was isolated from a hot spring in Japan. It belongs to the proposed family “Medusaviridae” in the phylum Nucleocytoviricota. Unlike other amoeba-infecting giant viruses, medusavirus initiates its DNA replication in the host nucleus without disrupting the nuclear membrane. Our RNA sequencing (RNA-seq) analysis of its infection course uncovered ordered viral gene expression profiles. We identified temporal expression clusters of viral genes and associated putative promoter motifs. The subcellular localization prediction showed a clear spatiotemporal correlation between gene expression timing and localization of the encoded proteins. Notably, the immediate early expression cluster was enriched in genes targeting the nucleus, suggesting the priority of remodeling the host intranuclear environment during infection. The transcriptional profile of amoeba genes was greatly altered postinfection.

One of the many challenges faced by RNA viruses is the maintenance of their genomes during infections of host cells. Members of the family Tombusviridae are plus-strand RNA viruses with unmodified triphosphorylated genomic 5′ termini. The tombusvirus Carnation Italian ringspot virus was used to investigate how it protects its RNA genome from attack by 5′-end-targeting degradation enzymes. In vivo and in vitro assays were employed to determine the role of genomic RNA structure in conferring protection from the 5′-to-3′ exoribonuclease Xrn. The results revealed that (i) the CIRV RNA genome is more resistant to Xrn than its sg mRNAs, (ii) the genomic 5′-untranslated region (UTR) folds into a compact RNA structure that effectively and independently prevents Xrn access, (iii) the RNA structure limiting 5′ access is formed by secondary and tertiary interactions that function cooperatively, (iv) the structure is also able to block access of RNA pyrophosphohydrolase to the genomic 5′ terminus, and (v) the RNA structure does not stall an actively digesting Xrn. Based on its proficiency at impeding Xrn 5′ access, we have termed this 5′-terminal structure an Xrn-evading RNA, or xeRNA. These and other findings demonstrate that the 5′UTR of the CIRV RNA genome folds into a complex structural conformation that helps to protect its unmodified 5′ terminus from enzymatic decay during infections. IMPORTANCE The plus-strand RNA genomes of plant viruses in the large family Tombusviridae are not 5′ capped. Here, we explored how a species in the type genus Tombusvirus protects its genomic 5′ end from cellular nuclease attack. Our results revealed that the 5′-terminal sequence of the CIRV genome folds into a complex RNA structure that limits access of the 5′-to-3′ exoribonuclease Xrn, thereby protecting it from processive degradation. The RNA conformation also impeded access of RNA pyrophosphohydrolase, which converts 5′-triphosphorylated RNA termini into 5′-monophosphorylated forms, the preferred substrate for Xrn. This study represents the first report of a higher-order RNA structure in an RNA plant virus genome independently conferring resistance to 5′-end-attacking cellular enzymes.

SARS-CoV-2 variants of concern (VOCs) have emerged as a global threat to the COVID-19 pandemic response. We implemented a combined approach to quickly detect known VOCs while continuously monitoring for evolving mutations of the virus. To rapidly detect VOCs, two real-time reverse transcriptase PCR assays were designed and implemented, targeting the spike gene H69/V70 deletion and the N501Y mutation. The H69/V70 deletion and N501Y mutation assays demonstrated accuracies of 98.3% (95% CI 93.8 to 99.8) and 100% (95% CI 96.8 to 100), limits of detection of 1,089 and 294 copies/ml, and percent coefficients of variation of 0.08 to 1.16% and 0 to 2.72% for the two gene targets, respectively. No cross-reactivity with common respiratory pathogens was observed with either assay. Implementation of these tests allowed the swift escalation in testing for VOCs from 2.2% to ∼100% of all SARS-CoV-2-positive samples over 12 January to 9 February 2021, and resulted in the detection of a rapid rise of B.1.1.7 cases within the province of Alberta, Canada. A prospective comparison of the VOC assays to genome sequencing for the detection of B.1.1.7, combined detection of P.1 and B.1.351, and wild-type (i.e., non-VOC) lineages showed sensitivities of 98.2 to 100%, specificities of 98.9 to 100%, positive predictive values of 76.9% to 100%, and negative predictive values of 96 to 100%. Variant screening results inform sampling strategies for regular surveillance by genome sequencing, thus allowing rapid identification of known VOCs while continuously monitoring the evolution of SARS-CoV-2 in the province. IMPORTANCE Different strains, or variants, of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2, the virus that causes COVID-19) have emerged that have higher levels of transmission, less susceptibility to our immune response, and possibly cause more severe disease than previous strains of the virus. Rapid detection of these variants of concern is important to help contain them and prevent them from spreading widely within the population. This study describes two newly developed tests that are able to identify and differentiate the variants of concern from regular strains of SARS-CoV-2. These tests are faster and simpler than the main, gold standard method of identifying variants of concern (genome sequencing). These tests also demonstrated a high correlation with genome sequencing and allowed for the rapid and accurate detection of the rise of B.1.1.7 (one of the variants of concern) in the province of Alberta, Canada.

The dynamics of quasispecies afford RNA viruses a great fitness on cell tropism and host range. To study the quasispecies features and the intra-host evolution of SARS-CoV-2, we collected nine confirmed patients and sequenced the haplotypes of spike gene using a single-molecule real-time platform. Fourteen samples were extracted from sputum, nasopharyngeal swabs, or stool, which in total produced 283,655 high-quality circular consensus sequences. We observed a stable quasispecies structure that one master mutant (mean abundance ∼0.70), followed by numerous minor mutants (mean abundance ∼1.21 × 10−3). Under high selective pressure, minor mutants may obtain a fitness advantage and become the master ones. The later predominant substitution D614G existed in the minor mutants of more than one early patient. An epidemic variant had a possibility to be independently originated from multiple hosts. The mutant spectrums covered ∼85% amino acid variations of public genomes (GISAID; frequency ≥ 0.1) and likely provided an advantage mutation pool for the current/future epidemic variants. Notably, 32 of 35 collected antibody escape substitutions were preexistent in the early quasispecies. Virus populations in different tissues/organs revealed potentially independent replications. The quasispecies complexity of sputum samples was significantly lower than that of nasopharyngeal swabs (P = 0.02). Evolution analysis revealed that three continuous S2 domains (HR1, CH, and CD) had undergone a positive selection. Cell fusion-related domains may play a crucial role in adapting to the intrahost immune system. Our findings suggested that future epidemiologic investigations and clinical interventions should consider the quasispecies information that has missed by routine single consensus genome. IMPORTANCE RNA virus population in a host does not consist of a consensus single haplotype but rather an ensemble of related sequences termed quasispecies. The dynamics of quasispecies afford SARS-CoV-2 a great ability on genetic fitness during intrahost evolution. The process is likely achieved by changing the genetic characteristics of key functional genes, such as the spike glycoprotein. Previous studies have applied the next-generation sequencing (NGS) technology to evaluate the quasispecies of SARS-CoV-2, and results indicated a low genetic diversity of the spike gene. However, the NGS platform cannot directly obtain the full haplotypes without assembling, and it is also difficult to predict the extremely low-frequency variations. Therefore, we introduced a single-molecule real-time technology to directly obtain the haplotypes of the RNA population and further study the quasispecies features and intrahost evolution of the spike gene.

Alphaviruses have positive-strand RNA genomes containing two open reading frames (ORFs). The first ORF encodes the nonstructural (ns) polyproteins P123 and P1234 that act as precursors for the subunits of the viral RNA replicase (nsP1 to nsP4). Processing of P1234 leads to the formation of a negative-strand replicase consisting of nsP4 (RNA polymerase) and P123 components. Subsequent processing of P123 results in a positive-strand replicase. The second ORF encoding the structural proteins is expressed via the synthesis of a subgenomic RNA. Alphavirus replicase is capable of using template RNAs that contain essential cis-active sequences. Here, we demonstrate that the replicases of nine alphaviruses, expressed in the form of separate P123 and nsP4 components, are active. Their activity depends on the abundance of nsP4. The match of nsP4 to its template strongly influences efficient subgenomic RNA synthesis. nsP4 of Barmah Forest virus (BFV) formed a functional replicase only with matching P123, while nsP4s of other alphaviruses were compatible also with several heterologous P123s. The P123 components of Venezuelan equine encephalitis virus and Sindbis virus (SINV) required matching nsP4s, while P123 of other viruses could form active replicases with different nsP4s. Chimeras of Semliki Forest virus, harboring the nsP4 of chikungunya virus, Ross River virus, BFV, or SINV were viable. In contrast, chimeras of SINV, harboring an nsP4 from different alphaviruses, exhibited a temperature-sensitive phenotype. These findings highlight the possibility for formation of new alphaviruses via recombination events and provide a novel approach for the development of attenuated chimeric viruses for vaccination strategies. IMPORTANCE A key element of every virus with an RNA genome is the RNA replicase. Understanding the principles of RNA replicase formation and functioning is therefore crucial for understanding and responding to the emergence of new viruses. Reconstruction of the replicases of nine alphaviruses from nsP4 and P123 polyproteins revealed that the nsP4 of the majority of alphaviruses, including the mosquito-specific Eilat virus, could form a functional replicase with P123 originating from a different virus, and the corresponding chimeric viruses were replication-competent. nsP4 also had an evident role in determining the template RNA preference and the efficiency of RNA synthesis. The revealed broad picture of the compatibility of the replicase components of alphaviruses is important for understanding the formation and functioning of the alphavirus RNA replicase and highlights the possibilities for recombination between different alphavirus species.

Bioinformatic analysis showed previously that a majority of promoters in the photoheterotrophic alphaproteobacterium Rhodobacter sphaeroides lack the thymine at the last position of the −10 element (−7T), a base that is very highly conserved in promoters in bacteria other than alphaproteobacteria. The absence of −7T was correlated with low promoter activity using purified R. sphaeroides RNA polymerase (RNAP), but the transcription factor CarD compensated by activating almost all promoters lacking −7T tested in vitro, including rRNA promoters. Here, we show that a previously uncharacterized R. sphaeroides promoter, the promoter for carD itself, has high basal activity relative to other tested R. sphaeroides promoters despite lacking −7T, and its activity is inhibited rather than activated by CarD. This high basal activity is dependent on a consensus-extended −10 element (TGn) and specific features in the spacer immediately upstream of the extended −10 element. CarD negatively autoregulates its own promoter by producing abortive transcripts, limiting promoter escape, and reducing full-length mRNA synthesis. This mechanism of negative regulation differs from that employed by classical repressors, in which the transcription factor competes with RNA polymerase for binding to the promoter, and with the mechanism of negative regulation used by transcription factors like DksA/ppGpp and TraR that allosterically inhibit the rate of open complex formation. IMPORTANCE R. sphaeroides CarD activates many promoters by binding directly to RNAP and DNA just upstream of the −10 element. In contrast, we show here that CarD inhibits its own promoter using the same interactions with RNAP and DNA used for activation. Inhibition results from increasing abortive transcript formation, thereby decreasing promoter escape and full-length RNA synthesis. We propose that the combined interactions of RNAP with CarD, with the extended −10 element and with features in the adjacent −10/−35 spacer DNA, stabilize the promoter complex, reducing promoter clearance. These findings support previous predictions that the effects of CarD on transcription can be either positive or negative, depending on the kinetic properties of the specific promoter.

Antimicrobial resistance in Neisseria gonorrhoeae has reached an alarming level, severely impacting the effective treatment of gonorrhea. Belonging to the resistance-nodulation-cell division (RND) superfamily of efflux transporters, the MtrD membrane protein of N. gonorrhoeae provides resistance to a broad range of antimicrobial compounds. A unique feature of MtrD is an 11-residue sequence (from N917 to P927 [N917-P927]) that connects transmembrane helices (TMS) 9 and 10; this sequence is not present in homologous RND proteins. This study explores the structural and functional roles of the N917-P927 region by means of mutant analysis and molecular dynamics simulations. We show that N917-P927 plays a key role in modulating substrate access to the binding cleft and influences the overall orientation of the protein within the inner membrane necessary for optimal functioning. Removal of N917-P927 significantly reduced MtrD-mediated resistance to a range of antimicrobials and mutations of three single amino acids impacted MtrD-mediated multidrug resistance. Furthermore, molecular dynamics simulations showed deletion of N917-P927 in MtrD may dysregulate access of the substrate to the binding cleft and closure of the substrate-binding pocket during the transport cycle. These findings indicate that N917-P927 is a key region for interacting with the inner membrane, conceivably influencing substrate capture from the membrane-periplasm interface and thus is essential for full multidrug resistance capacity of MtrD. IMPORTANCE The historical sexually transmitted infection gonorrhea continues to be a major public health concern with an estimated global annual incidence of 86.9 million cases. N. gonorrhoeae has been identified by the World Health Organization as one of the 12 antimicrobial-resistant bacterial species that poses the greatest risk to human health. As the major efflux pump in gonococci, the MtrD transporter contributes to the cell envelope barrier in this organism and pumps antimicrobials from the periplasm and inner membrane, resulting in resistance. This study demonstrates that a unique region of the MtrD protein that connects TMS 9 and TMS 10 forms a structure that may interact with the inner membrane positioning TMS 9 and stabilizing the protein facilitating substrate capture from the inner membrane-periplasm interface. Analysis of mutants of this region identified that it was essential for MtrD-mediated multidrug resistance. Characterization of the structure and function of this unique local region of MtrD has implications for drug efflux mechanisms used by related proteins and is important knowledge for development of antibiotics that bypass efflux.

Gene drives are genetic systems designed to efficiently spread a modification through a population. They have been designed almost exclusively in eukaryotic species, especially in insects. We recently developed a CRISPR-based gene drive system in herpesviruses that relies on similar mechanisms and could efficiently spread into a population of wild-type viruses. A common consequence of gene drives in insects is the appearance and selection of drive-resistant sequences that are no longer recognized by CRISPR-Cas9. In this study, we analyzed in cell culture experiments the evolution of resistance in a viral gene drive against human cytomegalovirus. We report that after an initial invasion of the wild-type population, a drive-resistant population is positively selected over time and outcompetes gene drive viruses. However, we show that targeting evolutionarily conserved sequences ensures that drive-resistant viruses acquire long-lasting mutations and are durably attenuated. As a consequence, and even though engineered viruses do not stably persist in the viral population, remaining viruses have a replication defect, leading to a long-term reduction of viral levels. This marks an important step toward developing effective gene drives in herpesviruses, especially for therapeutic applications. IMPORTANCE The use of defective viruses that interfere with the replication of their infectious parent after coinfecting the same cells—a therapeutic strategy known as viral interference—has recently generated a lot of interest. The CRISPR-based system that we recently reported for herpesviruses represents a novel interfering strategy that causes the conversion of wild-type viruses into new recombinant viruses and drives the native viral population to extinction. In this study, we analyzed how targeted viruses evolved resistance against the technology. Through numerical simulations and cell culture experiments with human cytomegalovirus, we showed that after the initial propagation, a resistant viral population is positively selected and outcompetes engineered viruses over time. We show, however, that targeting evolutionarily conserved sequences ensures that resistant viruses are mutated and attenuated, which leads to a long-term reduction of viral levels. This marks an important step toward the development of novel therapeutic strategies against herpesviruses.

SUMMARY The malaria parasite has a complex life cycle exhibiting phenotypic and morphogenic variations in two different hosts by existing in heterogeneous developmental states. To investigate this cellular heterogeneity of the parasite within the human host, we performed single-cell RNA sequencing of synchronized Plasmodium cells under control and temperature treatment conditions. Using the Malaria Cell Atlas (https://www.sanger.ac.uk/science/tools/mca) as a guide, we identified 9 subtypes of the parasite distributed across known intraerythrocytic stages. Interestingly, temperature treatment results in the upregulation of the AP2-G gene, the master regulator of sexual development in a small subpopulation of the parasites. Moreover, we identified a heterogeneous stress-responsive subpopulation (clusters 5, 6, and 7 [∼10% of the total population]) that exhibits upregulation of stress response pathways under normal growth conditions. We also developed an online exploratory tool that will provide new insights into gene function under normal and temperature stress conditions. Thus, our study reveals important insights into cell-to-cell heterogeneity in the parasite population under temperature treatment that will be instrumental toward a mechanistic understanding of cellular adaptation and population dynamics in Plasmodium falciparum. IMPORTANCE The malaria parasite has a complex life cycle exhibiting phenotypic variations in two different hosts accompanied by cell-to-cell variability that is important for stress tolerance, immune evasion, and drug resistance. To investigate cellular heterogeneity determined by gene expression, we performed single-cell RNA sequencing (scRNA-seq) of about 12,000 synchronized Plasmodium cells under physiologically relevant normal (37°C) and temperature stress (40°C) conditions phenocopying the cyclic bouts of fever experienced during malarial infection. In this study, we found that parasites exhibit transcriptional heterogeneity in an otherwise morphologically synchronized culture. Also, a subset of parasites is continually committed to gametocytogenesis and stress-responsive pathways. These observations have important implications for understanding the mechanisms of drug resistance generation and vaccine development against the malaria parasite.

The influenza A virus genome is comprised of eight single-stranded negative-sense viral RNA (vRNA) segments. Each of the eight vRNA segments contains segment-specific nonconserved noncoding regions (NCRs) of similar sequence and length in different influenza A virus strains. However, in the subtype-determinant segments, encoding hemagglutinin (HA) and neuraminidase (NA), the segment-specific noncoding regions are subtype specific, varying significantly in sequence and length at both the 3' and 5' termini among different subtypes. The significance of these subtype-specific noncoding regions (ssNCR) in the influenza virus replication cycle is not fully understood. In this study, we show that truncations of the 3'-end H1-subtype-specific noncoding region (H1-ssNCR) resulted in recombinant viruses with decreased HA vRNA replication and attenuated growth phenotype, although the vRNA replication was not affected in single-template RNP reconstitution assays. The attenuated viruses were unstable, and point mutations at nucleotide position 76 or 56 in the adjacent coding region of HA vRNA were found after serial passage. The mutations restored the HA vRNA replication and reversed the attenuated virus growth phenotype. We propose that the terminal noncoding and adjacent coding regions act synergistically to ensure optimal levels of HA vRNA replication in a multisegment environment. These results provide novel insights into the role of the 3'-end nonconserved noncoding regions and adjacent coding regions on template preference in multiple-segmented negative-strand RNA viruses. IMPORTANCE While most influenza A virus vRNA segments contain segment-specific nonconserved noncoding regions of similar length and sequence, these regions vary considerably both in length and sequence in the segments encoding HA and NA, the two major antigenic determinants of influenza A viruses. In this study, we investigated the function of the 3'-end H1-ssNCR and observed a synergistic effect between the 3'-end H1-ssNCR nucleotides and adjacent coding nucleotide(s) of the HA segment on template preference in a multisegment environment. The results unravel an additional level of complexity in the regulation of RNA replication in multiple-segmented negative-strand RNA viruses.

Genome-wide variation in SARS-CoV-2 reveals evolution and transmission dynamics which are critical considerations for disease control and prevention decisions. Here, we review estimates of the genome-wide viral mutation rates, summarize current COVID-19 case load in the province of Ontario, Canada (5 January 2021), and analyze published SARS-CoV-2 genomes from Ontario (collected prior to 24 November 2020) to test for more infectious genetic variants or lineages., Genome-wide variation in SARS-CoV-2 reveals evolution and transmission dynamics which are critical considerations for disease control and prevention decisions. Here, we review estimates of the genome-wide viral mutation rates, summarize current COVID-19 case load in the province of Ontario, Canada (5 January 2021), and analyze published SARS-CoV-2 genomes from Ontario (collected prior to 24 November 2020) to test for more infectious genetic variants or lineages. The reported mutation rate (∼10−6 nucleotide [nt]−1 cycle−1) for SARS-CoV-2 is typical for coronaviruses. Analysis of published SARS-CoV-2 genomes revealed that the G614 spike protein mutation has dominated infections in Ontario and that SARS-CoV-2 lineages present in Ontario have not differed significantly in their rate of spread. These results suggest that the SARS-CoV-2 population circulating in Ontario has not changed significantly to date. However, ongoing genome monitoring is essential for identification of new variants and lineages that may contribute to increased viral transmission.

The intracellular bacterial pathogen Brucella causes persistent infections in various mammalian species. To survive and replicate within macrophages, these bacteria must be able to withstand oxidative stresses and express the type IV secretion system (T4SS) to evade host immune responses. The extracytoplasmic function (ECF) sigma factor system is a major signal transduction mechanism in bacteria that senses environmental cues and responds by regulating gene expression. In this study, we defined an ECF σ bcrS and its cognate anti-σ factor abcS in Brucella melitensis M28 by conserved domain analysis and a protein interaction assay. BcrS directly activates an adjacent operon, bcrXQP, that encodes a methionine-rich peptide and a putative methionine sulfoxide reductase system, whereas AbcS is a negative regulator of bcrS and bcrXQP. The bcrS-abcS and bcrXQP operons can be induced by hypochlorous acid and contribute to hypochlorous acid resistance in vitro. Next, RNA sequencing analysis and genome-wide recognition sequence search identified the regulons of BcrS and AbcS. Interestingly, we found that BcrS positively influences T4SS expression in an AbcS-dependent manner and that AbcS also affects T4SS expression independently of BcrS. Last, we demonstrate that abcS is required for the maintenance of persistent infection, while bcrS is dispensable in a mouse infection model. Collectively, we conclude that BcrS and AbcS influence expression of multiple genes responsible for Brucella virulence traits. IMPORTANCEBrucella is a notorious intracellular pathogen that induces chronic infections in animals and humans. To survive and replicate within macrophages, these bacteria require a capacity to withstand oxidative stresses and to express the type IV secretion system (T4SS) to combat host immune responses. In this study, we characterized an extracytoplasmic function sigma/anti-sigma factor system that regulates resistance to reactive chlorine species and T4SS expression, thereby establishing a potential link between two crucial virulence traits of Brucella. Furthermore, the anti-sigma factor AbcS contributes to Brucella persistent infection of mice. Thus, this work provides novel insights into Brucella virulence regulation as well as a potential drug target for fighting Brucella infections.

Cyanophages are viruses that target cyanobacteria and directly control their abundance via viral lysis. Cyanobacteria are known to cause large blooms in water bodies, substantially contributing to oxygen depletion in bottom waters resulting in areas called dead zones., Up to 20% of prokaryotic organisms in the oceans are estimated to die every day due to viral infection and lysis. Viruses can therefore alter microbial diversity, community structure, and biogeochemical processes driven by these organisms. Cyanophages are viruses that infect and lyse cyanobacterial cells, adding bioavailable carbon and nutrients into the environment. Cyanobacteria are photosynthesizing bacteria, with some species capable of N2 fixation, which are known to form large blooms as well as resistant resting cells known as akinetes. Here, we investigated cyanophage diversity and community structure plus cyanobacteria in dead zone sediments. We sampled surface sediments and sequenced DNA and RNA, along an oxygen gradient—representing oxic, hypoxic, and anoxic conditions—in one of the world’s largest dead zones located in the Baltic Sea. Cyanophages were detected at all stations and, based on partial genome contigs, had a higher alpha diversity and different beta diversity in the hypoxic-anoxic sediments, suggesting that cyanobacteria in dead zone sediments and/or environmental conditions select for specific cyanophages. Some of these cyanophages can infect cyanobacteria with potential consequences for gene expression related to their photosystem and phosphate regulation. Top cyanobacterial genera detected in the anoxic sediment included Dolichospermum/Anabaena, Synechococcus, and Cyanobium. RNA transcripts classified to cyanobacteria were associated with numerous pathways, including anaerobic carbon metabolism and N2 fixation. Cyanobacterial blooms are known to fuel oxygen-depleted ecosystems with phosphorus (so-called internal loading), and our cyanophage data indicate the potential for viral lysis of cyanobacteria which might explain the high nutrient turnover in these environments. IMPORTANCE Cyanophages are viruses that target cyanobacteria and directly control their abundance via viral lysis. Cyanobacteria are known to cause large blooms in water bodies, substantially contributing to oxygen depletion in bottom waters resulting in areas called dead zones. Our knowledge of cyanophages in dead zones is very scarce, and so far, no studies have assembled partial cyanophage genomes and investigated their associated cyanobacteria in these dark and anoxic sediments. Here, we present the first study using DNA and RNA sequencing to investigate in situ diversity of cyanophages and cyanobacteria in dead zones. Our study shows that dead zone sediments contain different cyanophages compared to oxic sediments and suggest that these viruses are able to affect cyanobacterial photosystem and phosphate regulation. Furthermore, cyanophage-controlled lysis of cyanobacteria might also increase the turnover of carbon, phosphorus, and nitrogen in these oxygen-free environments at the bottom of the sea.

Diverse microbial communities drive biogeochemical cycles in Earth’s ocean, yet studying these organisms and processes is often limited by technological capabilities, especially in the deep ocean. In this study, we used a novel marine microbial incubator instrument capable of in situ experimentation to investigate microbial primary producers at deep-sea hydrothermal vents., Depressurization and sample processing delays may impact the outcome of shipboard microbial incubations of samples collected from the deep sea. To address this knowledge gap, we developed a remotely operated vehicle (ROV)-powered incubator instrument to carry out and compare results from in situ and shipboard RNA stable isotope probing (RNA-SIP) experiments to identify the key chemolithoautotrophic microbes and metabolisms in diffuse, low-temperature venting fluids from Axial Seamount. All the incubations showed microbial uptake of labeled bicarbonate primarily by thermophilic autotrophic Epsilonbacteraeota that oxidized hydrogen coupled with nitrate reduction. However, the in situ seafloor incubations showed higher abundances of transcripts annotated for aerobic processes, suggesting that oxygen was lost from the hydrothermal fluid samples prior to shipboard analysis. Furthermore, transcripts for thermal stress proteins such as heat shock chaperones and proteases were significantly more abundant in the shipboard incubations, suggesting that depressurization induced thermal stress in the metabolically active microbes in these incubations. Together, the results indicate that while the autotrophic microbial communities in the shipboard and seafloor experiments behaved similarly, there were distinct differences that provide new insight into the activities of natural microbial assemblages under nearly native conditions in the ocean. IMPORTANCE Diverse microbial communities drive biogeochemical cycles in Earth’s ocean, yet studying these organisms and processes is often limited by technological capabilities, especially in the deep ocean. In this study, we used a novel marine microbial incubator instrument capable of in situ experimentation to investigate microbial primary producers at deep-sea hydrothermal vents. We carried out identical stable isotope probing experiments coupled to RNA sequencing both on the seafloor and on the ship to examine thermophilic, microbial autotrophs in venting fluids from an active submarine volcano. Our results indicate that microbial communities were significantly impacted by the effects of depressurization and sample processing delays, with shipboard microbial communities being more stressed than seafloor incubations. Differences in metabolism were also apparent and are likely linked to the chemistry of the fluid at the beginning of the experiment. Microbial experimentation in the natural habitat provides new insights into understanding microbial activities in the ocean.

Compared to other human coronaviruses, the genetic diversity and evolution of human coronavirus 229E (HCoV-229E) are relatively understudied. We report a fatal case of COVID-19 pneumonia coinfected with HCoV-229E in Hong Kong. Genome sequencing of SARS-CoV-2 and HCoV-229E from a nasopharyngeal sample of the patient showed that the SARS-CoV-2 strain HK13 was most closely related to SARS-CoV-2 type strain Wuhan-Hu-1 (99.99% nucleotide identity), compatible with his recent history of travel to Wuhan. The HCoV-229E strain HK20-42 was most closely related to HCoV-229E strain SC0865 from the United States (99.86% nucleotide identity). To investigate if it may represent a newly emerged HCoV-229E genotype in Hong Kong, we retrieved 41 archived respiratory samples that tested positive for HCoV-229E from 2004 to 2019. Pneumonia and exacerbations of chronic airway diseases were common among infected patients. Complete RdRp, S, and N gene sequencing of the 41 HCoV-229E strains revealed that our contemporary HCoV-229E strains have undergone significant genetic drift with clustering of strains in chronological order. Two novel genogroups were identified, in addition to previously described genogroups 1 to 4, with recent circulating strains including strain HK20-42 belonging to novel genogroup 6. Positive selection was detected in the spike protein and receptor-binding domain, which may be important for viral evolution at the receptor-binding interphase. Molecular dating analysis showed that HCoV-229E shared the most recent common ancestor with bat and camel/alpaca 229E-related viruses at ∼1884, while camel/alpaca viruses had a relatively recent common ancestor at ∼1999. Further studies are required to ascertain the evolutionary origin and path of HCoV-229E. IMPORTANCE Since its first appearance in the 1960s, the genetic diversity and evolution of human coronavirus 229E (HCoV-229E) have been relatively understudied. In this study, we report a fatal case of COVID-19 coinfected with HCoV-229E in Hong Kong. Genome sequencing revealed that our SARS-CoV-2 strain is highly identical to the SARS-CoV-2 strain from Wuhan, compatible with the patient’s recent travel history, whereas our HCoV-229E strain in this study is highly identical to a recent strain in the United States. We also retrieved 41 archived HCoV-229E strains from 2004 to 2019 in Hong Kong for sequence analysis. Pneumonia and exacerbations of chronic airway diseases were common diagnoses among the 41 patients. The results showed that HCoV-229E was evolving in chronological order. Two novel genogroups were identified in addition to the four preexisting HCoV-229E genogroups, with recent circulating strains belonging to novel genogroup 6. Molecular clock analysis dated bat-to-human and bat-to-camelid transmission to as early as 1884.

Human tetherin is a host restriction factor that inhibits replication of enveloped viruses by blocking viral release. Tetherin has an unusual topology that includes an N-terminal cytoplasmic tail, a single transmembrane domain, an extracellular domain, and a C-terminal glycosylphosphatidylinositol anchor. Tetherin is not well conserved across species, so it inhibits viral replication in a species-specific manner. Thus, studies of tetherin activities from different species provide an important tool for understanding its antiviral mechanism. Here, we report cloning of equine tetherin and characterization of its antiviral activity. Equine tetherin shares 53%, 40%, 36%, and 34% amino acid sequence identity with feline, human, simian, and murine tetherins, respectively. Like the feline tetherin, equine tetherin has a shorter N-terminal domain than human tetherin. Equine tetherin is localized on the cell surface and strongly blocks human immunodeficiency virus type 1 (HIV-1), simian immunodeficiency virus (SIV), and equine infectious anemia virus (EIAV) release from virus-producing cells. The antiviral activity of equine tetherin is neutralized by EIAV envelope protein, but not by the HIV-1 accessory protein Vpu, which is a human tetherin antagonist, and EIAV envelope protein does not counteract human tetherin. These results shed new light on our understanding of the species-specific tetherin antiviral mechanism.

Typical enteropathogenic Escherichia coli (tEPEC) is a leading cause of diarrhea and associated death in children worldwide. Atypical EPEC (aEPEC) lacks the plasmid encoding bundle-forming pili and is considered less virulent, but the molecular mechanism of virulence is poorly understood. We recently identified kittens as a host for aEPEC where intestinal epithelial colonization was associated with diarrheal disease and death. The purposes of this study were to (i) determine the genomic similarity between kitten aEPEC and human aEPEC isolates and (ii) identify genotypic or phenotypic traits associated with virulence in kitten aEPEC. We observed no differences between kitten and human aEPEC in core genome content or gene cluster sequence identities, and no distinguishing genomic content was observed between aEPEC isolates from kittens with nonclinical colonization (NC) versus those with lethal infection (LI). Variation in adherence patterns and ability to aggregate actin in cultured cells mirrored descriptions of human aEPEC. The aEPEC isolated from kittens with LI were significantly more motile than isolates from kittens with NC. Kittens may serve as a reservoir for aEPEC that is indistinguishable from human aEPEC isolates and may provide a needed comparative animal model for the study of aEPEC pathogenesis. Motility seems to be an important factor in pathogenesis of LI associated with aEPEC in kittens.

The IS26 family includes the ISs that are commonly found associated with antibiotic resistance genes in multiply resistant Gram-negative and Gram-positive bacteria. IS26 is most prevalent in Gram-negative species and can generate the clusters of antibiotic resistance genes interspersed with directly oriented IS26 seen in multiply resistant pathogens., IS26 forms cointegrates using two distinct routes, a copy-in mechanism involving one insertion sequence (IS) and a target and a targeted conservative mechanism involving two ISs in different DNA molecules. In this study, the ability of IS26 and some close relatives, IS1006, IS1008, and a natural hybrid, IS1006/IS1008, which are found predominantly in Acinetobacter spp., to interact was examined. IS1006/1008 consists of 175 bp from IS1006 at the left end, with the remainder from IS1008. These ISs all have the same 14-bp terminal inverted repeats, and the Tnp26, Tnp1006, and Tnp1008 transposases, with pairwise identities of 83.7% to 93.1%, should be able to recognize each other’s ends. In a recA-negative Escherichia coli strain, IS1006, IS1008, and IS1006/1008 each formed cointegrates via the copy-in route and via the targeted conservative route, albeit at frequencies for the targeted reaction at least 10-fold lower than for IS26. However, using mixed pairs, targeted cointegration was detected only when IS1008 was paired with the IS1006/1008 hybrid, which also encodes Tnp1008, and the targeted cointegrates formed all arose from a reaction occurring at the end where the DNA sequences are identical. The reaction also occurred at the end with extended DNA identity using IS26 paired with IS26::catA1, an artificially constructed IS26 derivative that includes the catA1 gene. Thus, both identical transposases and identical DNA sequences at the reacting end were required. These features indicate that the targeted conservative pathway proceeds via a single transposase-catalyzed strand transfer, followed by migration and resolution of the Holliday junction formed. IMPORTANCE The IS26 family includes the ISs that are commonly found associated with antibiotic resistance genes in multiply resistant Gram-negative and Gram-positive bacteria. IS26 is most prevalent in Gram-negative species and can generate the clusters of antibiotic resistance genes interspersed with directly oriented IS26 seen in multiply resistant pathogens. This ability relies on the novel dual mechanistic capabilities of IS26 family members. However, the mechanism underlying the recently discovered targeted conservative mode of cointegrate formation mediated by IS26, IS257/IS431, and IS1216, which is unlike any previously studied IS movement mechanism, is not well understood. An important question is what features of the IS and the transposase are key to allowing IS26 family members to undertake targeted conservative reaction. In this study, this question was addressed using mixed-partner crosses involving IS26 and naturally occurring close relatives of IS26 that are found near resistance genes in Acinetobacter baumannii and are widespread in Acinetobacter species.

Viral genome sequencing has guided our understanding of the spread and extent of genetic diversity of SARS-CoV-2 during the COVID-19 pandemic. SARS-CoV-2 viral genomes are usually sequenced from nasopharyngeal swabs of individual patients to track viral spread., Viral genome sequencing has guided our understanding of the spread and extent of genetic diversity of SARS-CoV-2 during the COVID-19 pandemic. SARS-CoV-2 viral genomes are usually sequenced from nasopharyngeal swabs of individual patients to track viral spread. Recently, RT-qPCR of municipal wastewater has been used to quantify the abundance of SARS-CoV-2 in several regions globally. However, metatranscriptomic sequencing of wastewater can be used to profile the viral genetic diversity across infected communities. Here, we sequenced RNA directly from sewage collected by municipal utility districts in the San Francisco Bay Area to generate complete and nearly complete SARS-CoV-2 genomes. The major consensus SARS-CoV-2 genotypes detected in the sewage were identical to clinical genomes from the region. Using a pipeline for single nucleotide variant calling in a metagenomic context, we characterized minor SARS-CoV-2 alleles in the wastewater and detected viral genotypes which were also found within clinical genomes throughout California. Observed wastewater variants were more similar to local California patient-derived genotypes than they were to those from other regions within the United States or globally. Additional variants detected in wastewater have only been identified in genomes from patients sampled outside California, indicating that wastewater sequencing can provide evidence for recent introductions of viral lineages before they are detected by local clinical sequencing. These results demonstrate that epidemiological surveillance through wastewater sequencing can aid in tracking exact viral strains in an epidemic context.

Novel recombinants, generated from inter- and intraspecies recombination of norovirus lineages, often emerge and pose a threat to public health. However, the factors determining emergence of these particular recombinants from all possible combinations of parental lineages remain largely unknown., Genetic recombination plays a pivotal role in the appearance of human norovirus recombinants that cause global epidemics. However, the factors responsible for the appearance of these recombinants remains largely unknown. In this study, we revealed a selective pressure that restricts parental combinations leading to the emergence of norovirus recombinants. To investigate traces of emerging novel recombinants and their parents in the human population, we isolated mass nucleotide sequence clones of human norovirus genogroups I and II in sewage-affected waters over a 4-year sampling period. Fourteen different phylogenetic combinations of recombinants and their parents were defined from the dozens of phylogenetic lineages circulating in the human population. To evaluate the probability of these combinations, parental lineages of each recombinant were categorized into two groups as HP (relatively higher-competitiveness parents) and LP (relatively lower-competitiveness parents), according to their relative detection frequency. Strong categorization of HP and LP was confirmed by tests with modified data and additional variables. An algorithm that was developed in this study to visualize the chance of mixed infection between parents revealed that HP lineages have a higher chance of mixed infection than LP lineages in the human population. Three parental pairing types in recombinants were defined: HP-HP, HP-LP, and LP-LP. Among these, most recombinants were identified as HP-LP, despite the prediction of dominant emergence of HP-HP-type recombinants. These results suggest that nature favors recombinants of human norovirus that originate from parental pairing of heterogeneous competitiveness. IMPORTANCE Novel recombinants, generated from inter- and intraspecies recombination of norovirus lineages, often emerge and pose a threat to public health. However, the factors determining emergence of these particular recombinants from all possible combinations of parental lineages remain largely unknown. Therefore, current investigations on these recombinants are inevitably limited to postepidemic analyses, which merely identify genetic or phenotypic changes in the newly emerged recombinants compared to their parents. Here, we provide a new theoretical concept that emergence of novel recombinants could be explained by a combination of parental noroviruses thriving in the human population and those circulating at lower levels. This study could provide an additional and important rationale for the proactive environmental monitoring of potential future epidemics due to viral recombinants.

Mycoplasma bovis causes serious infections in ruminants, leading to huge economic losses. Lipoproteins are key components of the mycoplasma membrane and are believed to function in nutrient acquisition, adherence, enzymatic interactions with the host, and induction of the host’s immune response to infection. Many genes of M. bovis have not been assigned functions, in part because of their low sequence similarity with other bacteria, making it difficult to extrapolate gene functions. This study examined functions of a surface-localized leucine-rich repeat (LRR) lipoprotein encoded by mbfN of M. bovis PG45. Homologs of MbfN were detected as 48-kDa peptides by Western blotting in all the strains of M. bovis included in this study, with the predicted 70-kDa full-length polypeptide detected in some strains. Sequence analysis of the gene revealed the absence in some strains of a region encoding the carboxyl-terminal 147 amino acids found in strain PG45, which could account for the variation detected by immunoblotting. In silico analysis of MbfN suggested that it may have an adhesion-related function. In vitro binding assays confirmed MbfN to be a fibronectin and heparin-binding protein. Disruption of mbfN in M. bovis PG45 significantly reduced (P = 0.033) the adherence of M. bovis PG45 to MDBK cells in vitro, demonstrating the role of MbfN as an adhesin. IMPORTANCE Experimental validation of the putative functions of genes in M. bovis will advance our understanding of the basic biology of this economically important pathogen and is crucial in developing prevention strategies. This study demonstrated the extracellular matrix binding ability of a novel immunogenic lipoprotein of M. bovis, and the role of this protein in adhesion by M. bovis suggests that it could play a role in virulence.

The TyrR transcription factor controls the expression of genes for the uptake and biosynthesis of aromatic amino acids in Escherichia coli. In the plant-associated and clinically significant proteobacterium Enterobacter ludwigii UW5, the TyrR orthologue was previously shown to regulate genes that encode enzymes for synthesis of the plant hormone indole-3-acetic acid and for gluconeogenesis, indicating a broader function for the transcription factor. This study aimed to delineate the TyrR regulon of E. ludwigii by comparing the transcriptomes of the wild type and a tyrR deletion strain. In E. ludwigii, TyrR positively or negatively regulates the expression of over 150 genes. TyrR downregulated expression of envelope stress response regulators CpxR and CpxP through interaction with a DNA binding site in the intergenic region between divergently transcribed cpxP and cpxR. Repression of cpxP was alleviated by tyrosine. Methyltransferase gene dmpM, which is possibly involved in antibiotic synthesis, was strongly activated in the presence of tyrosine and phenylalanine by TyrR binding to its promoter region. TyrR also regulated expression of genes for aromatic catabolism and anaerobic respiration. Our findings suggest that the E. ludwigii TyrR regulon has diverged from that of E. coli to include genes for survival in the diverse environments that this bacterium inhabits and illustrate the expansion and plasticity of transcription factor regulons. IMPORTANCE Genome-wide RNA sequencing revealed a broader regulatory role for the TyrR transcription factor in the ecologically versatile bacterium Enterobacter ludwigii beyond that of aromatic amino acid synthesis and transport that constitute the role of the TyrR regulon of E. coli. In E. ludwigii, a plant symbiont and human gut commensal, the TyrR regulon is expanded to include genes that are beneficial for plant interactions and response to stresses. Identification of the genes regulated by TyrR provides insight into the mechanisms by which the bacterium adapts to its environment.