Establishment and partial characterisation of a new cell line derived from adult tissues of the tsetse fly Glossina morsitans morsitans.

BACKGROUND: Insect cell lines play a vital role in many aspects of research on disease vectors and agricultural pests. The tsetse fly Glossina morsitans morsitans is an important vector of salivarian trypanosomes in sub-Saharan Africa and, as such, is a major constraint on human health and agricultural development in the region. METHODS: Here, we report establishment and partial characterisation of a cell line, GMA/LULS61, derived from tissues of adult female G. m. morsitans. GMA/LULS61 cells, grown at 28 °C in L-15 (Leibovitz) medium supplemented with foetal bovine serum and tryptose phosphate broth, have been taken through 23 passages to date and can be split 1:1 at 2-week intervals. Karyotyping at passage 17 revealed a predominantly haploid chromosome complement. Species origin and absence of contaminating bacteria were confirmed by PCR amplification and sequencing of fragments of the COI gene and pan-bacterial 16S rRNA gene respectively. However, PCR screening of RNA extracted from GMA/LULS61 cells confirmed presence of the recently described Glossina morsitans morsitans iflavirus and Glossina morsitans morsitans negevirus, but absence of Glossina pallipides salivary gland hypertrophy virus. GMA/LULS61 cells supported infection and growth of 6/7 different insect-derived strains of the intracellular bacterial symbiont Wolbachia. CONCLUSIONS: The GMA/LULS61 cell line has potential for application in a variety of studies investigating the biology of G. m. morsitans and its associated pathogenic and symbiotic microorganisms.

Emerging Variants of Canine Enteric Coronavirus Associated with Outbreaks of Gastroenteric Disease.

A 2022 canine gastroenteritis outbreak in the United Kingdom was associated with circulation of a new canine enteric coronavirus closely related to a 2020 variant with an additional spike gene recombination. The variants are unrelated to canine enteric coronavirus-like viruses associated with human disease but represent a model for coronavirus population adaptation.

Microbiome and mitogenomics of the chigger mite Pentidionis agamae: potential role as an Orientia vector and associations with divergent clades of Wolbachia and Borrelia.

BACKGROUND: Trombiculid mites are globally distributed, highly diverse arachnids that largely lack molecular resources such as whole mitogenomes for the elucidation of taxonomic relationships. Trombiculid larvae (chiggers) parasitise vertebrates and can transmit bacteria (Orientia spp.) responsible for scrub typhus, a zoonotic febrile illness. Orientia tsutsugamushi causes most cases of scrub typhus and is endemic to the Asia-Pacific Region, where it is transmitted by Leptotrombidium spp. chiggers. However, in Dubai, Candidatus Orientia chuto was isolated from a case of scrub typhus and is also known to circulate among rodents in Saudi Arabia and Kenya, although its vectors remain poorly defined. In addition to Orientia, chiggers are often infected with other potential pathogens or arthropod-specific endosymbionts, but their significance for trombiculid biology and public health is unclear. RESULTS: Ten chigger species were collected from rodents in southwestern Saudi Arabia. Chiggers were pooled according to species and screened for Orientia DNA by PCR. Two species (Microtrombicula muhaylensis and Pentidionis agamae) produced positive results for the htrA gene, although Ca. Orientia chuto DNA was confirmed by Sanger sequencing only in P. agamae. Metagenomic sequencing of three pools of P. agamae provided evidence for two other bacterial associates: a spirochaete and a Wolbachia symbiont. Phylogenetic analysis of 16S rRNA and multi-locus sequence typing genes placed the spirochaete in a clade of micromammal-associated Borrelia spp. that are widely-distributed globally with no known vector. For the Wolbachia symbiont, a genome assembly was obtained that allowed phylogenetic localisation in a novel, divergent clade. Cytochrome c oxidase I (COI) barcodes for Saudi Arabian chiggers enabled comparisons with global chigger diversity, revealing several cases of discordance with classical taxonomy. Complete mitogenome assemblies were obtained for the three P. agamae pools and almost 50 SNPs were identified, despite a common geographic origin. CONCLUSIONS: P. agamae was identified as a potential vector of Ca. Orientia chuto on the Arabian Peninsula. The detection of an unusual Borrelia sp. and a divergent Wolbachia symbiont in P. agamae indicated links with chigger microbiomes in other parts of the world, while COI barcoding and mitogenomic analyses greatly extended our understanding of inter- and intraspecific relationships in trombiculid mites.

Increased bacterial taxonomic and functional diversity is associated with impaired rotavirus vaccine immunogenicity in infants from India and Malawi.

The immunogenicity and effectiveness of oral rotavirus vaccines (ORVs) against severe rotavirus-associated gastroenteritis are impaired in low- and middle-income countries (LMICs) where the burden of disease is highest. Determining risk factors for impaired ORV response may help identify strategies to enhance vaccine effectiveness. In this study, we use metagenomic sequencing to provide a high-resolution taxonomic analysis of stool samples collected at 6 weeks of age (coinciding with the first ORV dose) during a prospective study of ORV immunogenicity in India and Malawi. We then analyse the functional capacity of the developing microbiome in these cohorts. Microbiome composition differed significantly between countries, although functional capacity was more similar than taxonomic composition. Our results confirm previously reported findings that the developing microbiome is more diverse in taxonomic composition in ORV non-seroconverters compared with seroconverters, and we additionally demonstrate a similar pattern in functional capacity. Although taxonomic or functional feature abundances are poor predictors of ORV response, we show that skews in the direction of associations within these microbiome data can be used to identify consistent markers of ORV response across LMIC infant cohorts. We also highlight the systemic under-representation of reference genes from LMICs that limit functional annotation in our study (7% and 13% annotation at pathway and enzyme commission level, respectively). Overall, higher microbiome diversity in early life may act as marker for impaired ORV response in India and Malawi, whilst a holistic perspective of functional capacity may be hidden in the "dark matter" of the microbiome.

Arsenophonus symbiosis with louse flies: multiple origins, coevolutionary dynamics, and metabolic significance.

IMPORTANCE: Insects that live exclusively on vertebrate blood utilize symbiotic bacteria as a source of essential compounds, e.g., B vitamins. In louse flies, the most frequent symbiont originated in genus Arsenophonus, known from a wide range of insects. Here, we analyze genomic traits, phylogenetic origins, and metabolic capacities of 11 Arsenophonus strains associated with louse flies. We show that in louse flies, Arsenophonus established symbiosis in at least four independent events, reaching different stages of symbiogenesis. This allowed for comparative genomic analysis, including convergence of metabolic capacities. The significance of the results is twofold. First, based on a comparison of independently originated Arsenophonus symbioses, it determines the importance of individual B vitamins for the insect host. This expands our theoretical insight into insect-bacteria symbiosis. The second outcome is of methodological significance. We show that the comparative approach reveals artifacts that would be difficult to identify based on a single-genome analysis.

Association Between Maternal Breastmilk Microbiota Composition and Rotavirus Vaccine Response in African, Asian, and European Infants: A Prospective Cohort Study.

BACKGROUND: Maternal breastmilk is a source of pre- and pro-biotics that impact neonatal gut microbiota colonization. Because oral rotavirus vaccines (ORVs) are administered at a time when infants are often breastfed, breastmilk microbiota composition may have a direct or indirect influence on vaccine take and immunogenicity. METHODS: Using standardized methods across sites, we compared breastmilk microbiota composition in relation to geographic location and ORV response in cohorts prospectively followed from birth to 18 weeks of age in India (n = 307), Malawi (n = 119), and the United Kingdom ([UK] n = 60). RESULTS: Breastmilk microbiota diversity was higher in India and Malawi than the UK across 3 longitudinal samples spanning weeks of life 1 to 13. Dominant taxa such as Streptococcus and Staphylococcus were consistent across cohorts; however, significant geographic differences were observed in the prevalence and abundance of common and rare genera throughout follow up. No consistent associations were identified between breastmilk microbiota composition and ORV outcomes including seroconversion, vaccine shedding after dose 1, and postvaccination rotavirus-specific immunoglobulin A level. CONCLUSIONS: Our findings suggest that breastmilk microbiota composition may not be a key factor in shaping trends in ORV response within or between countries.

The P323L substitution in the SARS-CoV-2 polymerase (NSP12) confers a selective advantage during infection.

BACKGROUND: The mutational landscape of SARS-CoV-2 varies at the dominant viral genome sequence and minor genomic variant population. During the COVID-19 pandemic, an early substitution in the genome was the D614G change in the spike protein, associated with an increase in transmissibility. Genomes with D614G are accompanied by a P323L substitution in the viral polymerase (NSP12). However, P323L is not thought to be under strong selective pressure. RESULTS: Investigation of P323L/D614G substitutions in the population shows rapid emergence during the containment phase and early surge phase during the first wave. These substitutions emerge from minor genomic variants which become dominant viral genome sequence. This is investigated in vivo and in vitro using SARS-CoV-2 with P323 and D614 in the dominant genome sequence and L323 and G614 in the minor variant population. During infection, there is rapid selection of L323 into the dominant viral genome sequence but not G614. Reverse genetics is used to create two viruses (either P323 or L323) with the same genetic background. L323 shows greater abundance of viral RNA and proteins and a smaller plaque morphology than P323. CONCLUSIONS: These data suggest that P323L is an important contribution in the emergence of variants with transmission advantages. Sequence analysis of viral populations suggests it may be possible to predict the emergence of a new variant based on tracking the frequency of minor variant genomes. The ability to predict an emerging variant of SARS-CoV-2 in the global landscape may aid in the evaluation of medical countermeasures and non-pharmaceutical interventions.

Genomic and cytogenetic analysis of the Ceratitis capitata temperature-sensitive lethal region.

Genetic sexing strains (GSS) are an important tool in support of sterile insect technique (SIT) applications against insect pests and disease vectors. The yet unknown temperature-sensitive lethal (tsl) gene and the recently identified white pupae (wp) gene have been used as selectable markers in the most successful GSS developed so far, the Ceratitis capitata (medfly) VIENNA 8 GSS. The molecular identification of the tsl gene may open the way for its use as a marker for the development of GSS in other insect pests and disease vectors of SIT importance. Prior studies have already shown that the tsl gene is located on the right arm of chromosome 5, between the wp and Zw loci (tsl genomic region). In the present study, we used genomic, transcriptomic, bioinformatic, and cytogenetic approaches to characterize and analyze this genomic region in wild-type and tsl mutant medfly strains. Our results suggested the presence of 561 genes, with 322 of them carrying SNPs and/or insertion-deletion (indel) mutations in the tsl genomic region. Furthermore, comparative transcriptomic analysis indicated the presence of 32 differentially expressed genes, and bioinformatic analysis revealed the presence of 33 orthologs with a described heat-sensitive phenotype of Drosophila melanogaster in this region. These data can be used in functional genetic studies to identify the tsl gene(s) and the causal mutation(s) responsible for the temperature-sensitive lethal phenotype in medfly, and potentially additional genes causing a similar phenotype.

Timing and delivery route effects of cecal microbiome transplants on Salmonella Typhimurium infections in chickens: potential for in-hatchery delivery of microbial interventions.

BACKGROUND: Exposure to microbes early in life has long-lasting effects on microbial community structure and function of the microbiome. However, in commercial poultry settings chicks are reared as a single-age cohort with no exposure to adult birds which can have profound effects on microbiota development and subsequent pathogen challenge. Microbiota manipulation is a proven and promising strategy to help reduce pathogen load and transmission within broiler flocks. However, administration of microbiota transplant products in a hatchery setting may prove challenging. Effective administration strategies are dependent on key factors, such as; the age of chicks receiving interventions and mode of delivery. This study aimed to assess these two aspects to provide supporting evidence towards microbiome manipulation strategies for use in commercial hatcheries. RESULTS: Manipulation of the microbiota between 4 and 72 h of hatch markedly reduced faecal shedding and colonisation with the foodborne pathogen Salmonella enterica serovar Typhimurium (ST4/74). Administration of transplant material via spray or gel drop delivery systems had minimal effect on the protection conferred with fewer birds in transplant groups shown to shed ST4/74 in the faeces compared to PBS-gavaged control birds. Analysis of the microbiome following transplantation demonstrated that all transplant groups had higher diversity and species richness than non-transplant groups during the first week of life and the early stages of infection with ST47/4.The relative abundance of the bacterium Faecalibacterium prausnitzii was significantly higher in CMT groups compared to PBS controls. The presence of F. prausnitzii was also shown to increase in PBS-challenged birds compared to unchallenged birds potentially indicating a role of this bacterium in limiting Salmonella infections. CONCLUSIONS: This study demonstrated that administration of microbiome transplants, using methods that would align with hatchery practices, effectively reduced colonisation and shedding of Salmonella in chickens. Age of chicks at microbiome administration had limited effect on the diversity and composition of the microbiome and conferred protection against Salmonella infections. Traditional hatchery delivery systems, such as spray or gel-drop, are sufficient to transfer donor material, alter the microbiome and confer protection against Salmonella. This study helps highlight the opportunity for use of microbiome modification methods within the hatchery.

Characterisation of SARS-CoV-2 genomic variation in response to molnupiravir treatment in the AGILE Phase IIa clinical trial.

Molnupiravir is an antiviral, currently approved by the UK Medicines and Healthcare products Regulatory Agency (MHRA) for treating at-risk COVID-19 patients, that induces lethal error catastrophe in SARS-CoV-2. How this drug-induced mechanism of action might impact the emergence of resistance mutations is unclear. To investigate this, we used samples from the AGILE Candidate Specific Trial (CST)-2 (clinical trial number NCT04746183). The primary outcomes of AGILE CST-2 were to measure the drug safety and antiviral efficacy of molnupiravir in humans (180 participants randomised 1:1 with placebo). Here, we describe the pre-specified exploratory virological endpoint of CST-2, which was to determine the possible genomic changes in SARS-CoV-2 induced by molnupiravir treatment. We use high-throughput amplicon sequencing and minor variant analysis to characterise viral genomics in each participant whose longitudinal samples (days 1, 3 and 5 post-randomisation) pass the viral genomic quality criteria (n = 59 for molnupiravir and n = 65 for placebo). Over the course of treatment, no specific mutations were associated with molnupiravir treatment. We find that molnupiravir significantly increased the transition:transversion mutation ratio in SARS-CoV-2, consistent with the model of lethal error catastrophe. This study highlights the utility of examining intra-host virus populations to strengthen the prediction, and surveillance, of potential treatment-emergent adaptations.

Recovering individual haplotypes and a contiguous genome assembly from pooled long-read sequencing of the diamondback moth (Lepidoptera: Plutellidae).

The assembly of divergent haplotypes using noisy long-read data presents a challenge to the reconstruction of haploid genome assemblies, due to overlapping distributions of technical sequencing error, intralocus genetic variation, and interlocus similarity within these data. Here, we present a comparative analysis of assembly algorithms representing overlap-layout-consensus, repeat graph, and de Bruijn graph methods. We examine how postprocessing strategies attempting to reduce redundant heterozygosity interact with the choice of initial assembly algorithm and ultimately produce a series of chromosome-level assemblies for an agricultural pest, the diamondback moth, Plutella xylostella (L.). We compare evaluation methods and show that BUSCO analyses may overestimate haplotig removal processing in long-read draft genomes, in comparison to a k-mer method. We discuss the trade-offs inherent in assembly algorithm and curation choices and suggest that "best practice" is research question dependent. We demonstrate a link between allelic divergence and allele-derived contig redundancy in final genome assemblies and document the patterns of coding and noncoding diversity between redundant sequences. We also document a link between an excess of nonsynonymous polymorphism and haplotigs that are unresolved by assembly or postassembly algorithms. Finally, we discuss how this phenomenon may have relevance for the usage of noisy long-read genome assemblies in comparative genomics.

Analysis of SARS-CoV-2 known and novel subgenomic mRNAs in cell culture, animal model, and clinical samples using LeTRS, a bioinformatic tool to identify unique sequence identifiers.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has a complex strategy for the transcription of viral subgenomic mRNAs (sgmRNAs), which are targets for nucleic acid diagnostics. Each of these sgmRNAs has a unique 5' sequence, the leader-transcriptional regulatory sequence gene junction (leader-TRS junction), that can be identified using sequencing. High-resolution sequencing has been used to investigate the biology of SARS-CoV-2 and the host response in cell culture and animal models and from clinical samples. LeTRS, a bioinformatics tool, was developed to identify leader-TRS junctions and can be used as a proxy to quantify sgmRNAs for understanding virus biology. LeTRS is readily adaptable for other coronaviruses such as Middle East respiratory syndrome coronavirus or a future newly discovered coronavirus. LeTRS was tested on published data sets and novel clinical samples from patients and longitudinal samples from animal models with coronavirus disease 2019. LeTRS identified known leader-TRS junctions and identified putative novel sgmRNAs that were common across different mammalian species. This may be indicative of an evolutionary mechanism where plasticity in transcription generates novel open reading frames, which can then subject to selection pressure. The data indicated multiphasic abundance of sgmRNAs in two different animal models. This recapitulates the relative sgmRNA abundance observed in cells at early points in infection but not at late points. This pattern is reflected in some human nasopharyngeal samples and therefore has implications for transmission models and nucleic acid-based diagnostics. LeTRS provides a quantitative measure of sgmRNA abundance from sequencing data. This can be used to assess the biology of SARS-CoV-2 (or other coronaviruses) in clinical and nonclinical samples, especially to evaluate different variants and medical countermeasures that may influence viral RNA synthesis.

Genomic diversity across the Rickettsia and 'Candidatus Megaira' genera and proposal of genus status for the Torix group.

Members of the bacterial genus Rickettsia were originally identified as causative agents of vector-borne diseases in mammals. However, many Rickettsia species are arthropod symbionts and close relatives of 'Candidatus Megaira', which are symbiotic associates of microeukaryotes. Here, we clarify the evolutionary relationships between these organisms by assembling 26 genomes of Rickettsia species from understudied groups, including the Torix group, and two genomes of 'Ca. Megaira' from various insects and microeukaryotes. Our analyses of the new genomes, in comparison with previously described ones, indicate that the accessory genome diversity and broad host range of Torix Rickettsia are comparable to those of all other Rickettsia combined. Therefore, the Torix clade may play unrecognized roles in invertebrate biology and physiology. We argue this clade should be given its own genus status, for which we propose the name 'Candidatus Tisiphia'.

Impact of maternal antibodies and microbiota development on the immunogenicity of oral rotavirus vaccine in African, Indian, and European infants.

Identifying risk factors for impaired oral rotavirus vaccine (ORV) efficacy in low-income countries may lead to improvements in vaccine design and delivery. In this prospective cohort study, we measure maternal rotavirus antibodies, environmental enteric dysfunction (EED), and bacterial gut microbiota development among infants receiving two doses of Rotarix in India (n = 307), Malawi (n = 119), and the UK (n = 60), using standardised methods across cohorts. We observe ORV shedding and seroconversion rates to be significantly lower in Malawi and India than the UK. Maternal rotavirus-specific antibodies in serum and breastmilk are negatively correlated with ORV response in India and Malawi, mediated partly by a reduction in ORV shedding. In the UK, ORV shedding is not inhibited despite comparable maternal antibody levels to the other cohorts. In both India and Malawi, increased microbiota diversity is negatively correlated with ORV immunogenicity, suggesting that high early-life microbial exposure may contribute to impaired vaccine efficacy.

Isolation in Natural Host Cell Lines of Wolbachia Strains wPip from the Mosquito Culex pipiens and wPap from the Sand Fly Phlebotomus papatasi.

Endosymbiotic intracellular bacteria of the genus Wolbachia are harboured by many species of invertebrates. They display a wide range of developmental, metabolic and nutritional interactions with their hosts and may impact the transmission of arboviruses and protozoan parasites. Wolbachia have occasionally been isolated during insect cell line generation. Here, we report the isolation of two strains of Wolbachia, wPip and wPap, during cell line generation from their respective hosts, the mosquito Culex pipiens and the sand fly Phlebotomus papatasi. wPip was pathogenic for both new C. pipiens cell lines, CPE/LULS50 and CLP/LULS56, requiring tetracycline treatment to rescue the lines. In contrast, wPap was tolerated by the P. papatasi cell line PPL/LULS49, although tetracycline treatment was applied to generate a Wolbachia-free subline. Both Wolbachia strains were infective for a panel of heterologous insect and tick cell lines, including two novel lines generated from the sand fly Lutzomyia longipalpis, LLE/LULS45 and LLL/LULS52. In all cases, wPip was more pathogenic for the host cells than wPap. These newly isolated Wolbachia strains, and the novel mosquito and sand fly cell lines reported here, will add to the resources available for research on host-endosymbiont relationships, as well as on C. pipiens, P. papatasi, L. longipalpis and the pathogens that they transmit.

A gene expression panel for estimating age in males and females of the sleeping sickness vector Glossina morsitans.

Many vector-borne diseases are controlled by methods that kill the insect vectors responsible for disease transmission. Recording the age structure of vector populations provides information on mortality rates and vectorial capacity, and should form part of the detailed monitoring that occurs in the wake of control programmes, yet tools for obtaining estimates of individual age remain limited. We investigate the potential of using markers of gene expression to predict age in tsetse flies, which are the vectors of deadly and economically damaging African trypanosomiases. We use RNAseq to identify candidate expression markers, and test these markers using qPCR in laboratory-reared Glossina morsitans morsitans of known age. Measuring the expression of six genes was sufficient to obtain a prediction of age with root mean squared error of less than 8 days, while just two genes were sufficient to classify flies into age categories of ≤15 and >15 days old. Further testing of these markers in field-caught samples and in other species will determine the accuracy of these markers in the field.

Generation and transmission of interlineage recombinants in the SARS-CoV-2 pandemic.

We present evidence for multiple independent origins of recombinant SARS-CoV-2 viruses sampled from late 2020 and early 2021 in the United Kingdom. Their genomes carry single-nucleotide polymorphisms and deletions that are characteristic of the B.1.1.7 variant of concern but lack the full complement of lineage-defining mutations. Instead, the remainder of their genomes share contiguous genetic variation with non-B.1.1.7 viruses circulating in the same geographic area at the same time as the recombinants. In four instances, there was evidence for onward transmission of a recombinant-origin virus, including one transmission cluster of 45 sequenced cases over the course of 2 months. The inferred genomic locations of recombination breakpoints suggest that every community-transmitted recombinant virus inherited its spike region from a B.1.1.7 parental virus, consistent with a transmission advantage for B.1.1.7's set of mutations.

Inter-kingdom relationships in Crohn's disease explored using a multi-omics approach.

The etiology of Crohn's disease (CD) is multifactorial. Bacterial and fungal microbiota are involved in the onset and/or progression of the disease. A bacterial dysbiosis in CD patients is accepted; however, less is known about the mycobiome and the relationships between the two communities. We investigated the interkingdom relationships, their metabolic consequences, and the changes in the fungal community during relapse and remission in CD.Two cohorts were evaluated: a British cohort (n = 63) comprising CD and ulcerative colitis patients, and controls. The fungal and bacterial communities of biopsy and fecal samples were analyzed, with the fecal volatiles; datasets were also integrated; and a Dutch cohort (n = 41) comprising CD patients and healthy controls was analyzed for stability of the gut mycobiome.A dysbiosis of the bacterial community was observed in biopsies and stool. Results suggest Bacteroides is likely key in CD and may modulate Candida colonization. A dysbiosis of the fungal community was observed only in the Dutch cohort; Malassezia and Candida were increased in patients taking immunosuppressants. Longitudinal analysis showed an increase in Cyberlindnera in relapse. Saccharomyces was dominant in all fecal samples, but not in biopsies, some of which did not yield fungal reads; amino acid degradation was the main metabolic change associated with CD and both bacteria and fungi might be implicated.We have shown that Bacteroides and yeasts may play a role in CD; understanding their role and relationship in the disease would shed new light on the development and treatment of CD.

Complete genome characterization of human noroviruses allows comparison of minor alleles during acute and chronic infections.

Human noroviruses (HuNoVs) circulate globally, affect all age groups and place a substantial burden upon health services. High genetic diversity leading to antigenic variation plays a significant role in HuNoV epidemiology, driving periodic global emergence of epidemic variants. Studies have suggested that immunocompromised individuals may be a reservoir for such epidemic variants, but studies investigating the diversity and emergence of HuNoV variants in immunocompetent individuals are underrepresented. To address this, we sequenced the genomes of HuNoVs present in samples collected longitudinally from one immunocompetent (acute infection) and one immunocompromised (chronic infection) patient. A broadly reactive HuNoV capture-based method was used to concentrate the virus present in these specimens prior to massively parallel sequencing to recover near complete viral genomes. Using a novel bioinformatics pipeline, we demonstrated that persistent minor alleles were present in both acute and chronic infections, and that minor allele frequencies represented a larger proportion of the population during chronic infection. In acute infection, minor alleles were more evenly spread across the genome, although present at much lower frequencies, and therefore difficult to discern from error. By contrast, in the chronic infection, more minor alleles were present in the minor structural protein. No non-synonymous minor alleles were detected in the major structural protein over the short sampling period of the HuNoV chronic infection, suggesting where immune pressure is variable or non-existent, epidemic variants could emerge over longer periods of infection by random chance.