Long-read sequencing reveals atypical mitochondrial genome structure in a New Zealand marine isopod

Most animal mitochondrial genomes are small, circular and structurally conserved. However, recent work indicates that diverse taxa possess unusual mitochondrial genomes. In Isopoda, species in multiple lineages have atypical and rearranged mitochondrial genomes. However, more species of this speciose taxon need to be evaluated to understand the evolutionary origins of atypical mitochondrial genomes in this group. In this study, we report the presence of an atypical mitochondrial structure in the New Zealand endemic marine isopod, Isocladus armatus. Data from long- and short-read DNA sequencing suggest that I. armatus has two mitochondrial chromosomes. The first chromosome consists of two mitochondrial genomes that have been inverted and fused together in a circular form, and the second chromosome consists of a single mitochondrial genome in a linearized form. This atypical mitochondrial structure has been detected in other isopod lineages, and our data from an additional divergent isopod lineage (Sphaeromatidae) lends support to the hypothesis that atypical structure evolved early in the evolution of Isopoda. Additionally, we find that an asymmetrical site previously observed across many species within Isopoda is absent in I. armatus, but confirm the presence of two asymmetrical sites recently reported in two other isopod species.

Tracking the international spread of SARS-CoV-2 lineages B.1.1.7 and B.1.351/501Y-V2 with grinch

Late in 2020, two genetically-distinct clusters of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) with mutations of biological concern were reported, one in the United Kingdom and one in South Africa. Using a combination of data from routine surveillance, genomic sequencing and international travel we track the international dispersal of lineages B.1.1.7 and B.1.351 (variant 501Y-V2). We account for potential biases in genomic surveillance efforts by including passenger volumes from location of where the lineage was first reported, London and South Africa respectively. Using the software tool grinch (global report investigating novel coronavirus haplotypes), we track the international spread of lineages of concern with automated daily reports, Further, we have built a custom tracking website (cov-lineages.org/global_report.html) which hosts this daily report and will continue to include novel SARS-CoV-2 lineages of concern as they are detected.

Use of Genomics to Track Coronavirus Disease Outbreaks, New Zealand

Real-time genomic sequencing has played a major role in tracking the global spread of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), contributing greatly to disease mitigation strategies. In August 2020, after having eliminated the virus, New Zealand experienced a second outbreak. During that outbreak, New Zealand used genomic sequencing in a primary role, leading to a second elimination of the virus. We generated genomes from 78% of the laboratory-confirmed samples of SARS-CoV-2 from the second outbreak and compared them with the available global genomic data. Genomic sequencing rapidly identified that virus causing the second outbreak in New Zealand belonged to a single cluster, thus resulting from a single introduction. However, successful identification of the origin of this outbreak was impeded by substantial biases and gaps in global sequencing data. Access to a broader and more heterogenous sample of global genomic data would strengthen efforts to locate the source of any new outbreaks.

Tracking the international spread of SARS-CoV-2 lineages B.1.1.7 and B.1.351/501Y-V2

Late in 2020, two genetically-distinct clusters of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) with mutations of biological concern were reported, one in the United Kingdom and one in South Africa. Using a combination of data from routine surveillance, genomic sequencing and international travel we track the international dispersal of lineages B.1.1.7 and B.1.351 (variant 501Y-V2). We account for potential biases in genomic surveillance efforts by including passenger volumes from location of where the lineage was first reported, London and South Africa respectively. Using the software tool grinch (global report investigating novel coronavirus haplotypes), we track the international spread of lineages of concern with automated daily reports, Further, we have built a custom tracking website (cov-lineages.org/global_report.html) which hosts this daily report and will continue to include novel SARS-CoV-2 lineages of concern as they are detected.

Concordant geographic and genetic structure revealed by genotyping-by-sequencing in a New Zealand marine isopod

Population genetic structure in the marine environment can be influenced by life-history traits such as developmental mode (biphasic, with distinct adult and larval morphology, and direct development, in which larvae resemble adults) or habitat specificity, as well as geography and selection. Developmental mode is thought to significantly influence dispersal, with direct developers expected to have much lower dispersal potential. However, this prediction can be complicated by the presence of geophysical barriers to dispersal. In this study, we use a panel of 8,020 SNPs to investigate population structure and biogeography over multiple spatial scales for a direct-developing species, the New Zealand endemic marine isopod Isocladus armatus. Because our sampling range is intersected by two well-known biogeographic barriers (the East Cape and the Cook Strait), our study provides an opportunity to understand how such barriers influence dispersal in direct developers. On a small spatial scale (20 km), gene flow between locations is extremely high, suggestive of an island model of migration. However, over larger spatial scales (600 km), populations exhibit a clear pattern of isolation-by-distance. Our results indicate that I. armatus exhibits significant migration across the hypothesized barriers and suggest that large-scale ocean currents associated with these locations do not present a barrier to dispersal. Interestingly, we find evidence of a north-south population genetic break occurring between Māhia and Wellington. While no known geophysical barrier is apparent in this area, it coincides with the location of a proposed border between bioregions. Analysis of loci under selection revealed that both isolation-by-distance and adaption may be contributing to the degree of population structure we have observed here. We conclude that developmental life history largely predicts dispersal in the intertidal isopod I. armatus. However, localized biogeographic processes can disrupt this expectation, and this may explain the potential meta-population detected in the Auckland region.

Rapid and inexpensive whole-genome sequencing of SARS-CoV-2 using 1200 bp tiled amplicons and Oxford Nanopore Rapid Barcoding

Rapid and cost-efficient whole-genome sequencing of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the virus that causes coronavirus disease 2019, is critical for understanding viral transmission dynamics. Here we show that using a new multiplexed set of primers in conjunction with the Oxford Nanopore Rapid Barcode library kit allows for faster, simpler, and less expensive SARS-CoV-2 genome sequencing. This primer set results in amplicons that exhibit lower levels of variation in coverage compared to other commonly used primer sets. Using five SARS-CoV-2 patient samples with Cq values between 20 and 31, we show that high-quality genomes can be generated with as few as 10 000 reads (∼5 Mbp of sequence data). We also show that mis-classification of barcodes, which may be more likely when using the Oxford Nanopore Rapid Barcode library prep, is unlikely to cause problems in variant calling. This method reduces the time from RNA to genome sequence by more than half compared to the more standard ligation-based Oxford Nanopore library preparation method at considerably lower costs.

Testing the advantages and disadvantages of short- and long- read eukaryotic metagenomics using simulated reads

BACKGROUND

The first step in understanding ecological community diversity and dynamics is quantifying community membership. An increasingly common method for doing so is through metagenomics. Because of the rapidly increasing popularity of this approach, a large number of computational tools and pipelines are available for analysing metagenomic data. However, the majority of these tools have been designed and benchmarked using highly accurate short read data (i.e. Illumina), with few studies benchmarking classification accuracy for long error-prone reads (PacBio or Oxford Nanopore). In addition, few tools have been benchmarked for non-microbial communities.

RESULTS

Here we compare simulated long reads from Oxford Nanopore and Pacific Biosciences (PacBio) with high accuracy Illumina read sets to systematically investigate the effects of sequence length and taxon type on classification accuracy for metagenomic data from both microbial and non-microbial communities. We show that very generally, classification accuracy is far lower for non-microbial communities, even at low taxonomic resolution (e.g. family rather than genus). We then show that for two popular taxonomic classifiers, long reads can significantly increase classification accuracy, and this is most pronounced for non-microbial communities.

CONCLUSIONS

This work provides insight on the expected accuracy for metagenomic analyses for different taxonomic groups, and establishes the point at which read length becomes more important than error rate for assigning the correct taxon.

Creation of a Dense Transposon Insertion Library Using Bacterial Conjugation in Enterobacterial Strains Such As Escherichia Coli or Shigella flexneri

Transposon mutagenesis is a method that allows gene disruption via the random genomic insertion of a piece of DNA called a transposon. The protocol below outlines a method for high efficiency transfer between bacterial strains of a plasmid harboring a transposon containing a kanamycin resistance marker. The plasmid-borne transposase is encoded by a variant tnp gene that inserts the transposon into the genome of the recipient strain with very low insertional bias. This method thus allows the creation of large mutant libraries in which transposons have been inserted into unique genomic positions in a recipient strain of either Escherichia coli or Shigella flexneri bacteria. By using bacterial conjugation, as opposed to other methods such as electroporation or chemical transformation, large libraries with hundreds of thousands of unique clones can be created. This yields high-density insertion libraries, with insertions occurring as frequently as every 4-6 base pairs in non-essential genes. This method is superior to other methods as it allows for an inexpensive, easy to use, and high efficiency method for the creation of a dense transposon insertion library. The transposon library can be used in downstream applications such as transposon sequencing (Tn-Seq), to infer genetic interaction networks, or more simply, in mutational (forward genetic) screens.

Combining Shigella Tn-seq data with gold-standard E. coli gene deletion data suggests rare transitions between essential and non-essential gene functionality

BACKGROUND

Gene essentiality - whether or not a gene is necessary for cell growth - is a fundamental component of gene function. It is not well established how quickly gene essentiality can change, as few studies have compared empirical measures of essentiality between closely related organisms.

RESULTS

Here we present the results of a Tn-seq experiment designed to detect essential protein coding genes in the bacterial pathogen Shigella flexneri 2a 2457T on a genome-wide scale. Superficial analysis of this data suggested that 481 protein-coding genes in this Shigella strain are critical for robust cellular growth on rich media. Comparison of this set of genes with a gold-standard data set of essential genes in the closely related Escherichia coli K12 BW25113 revealed that an excessive number of genes appeared essential in Shigella but non-essential in E. coli. Importantly, and in converse to this comparison, we found no genes that were essential in E. coli and non-essential in Shigella, implying that many genes were artefactually inferred as essential in Shigella. Controlling for such artefacts resulted in a much smaller set of discrepant genes. Among these, we identified three sets of functionally related genes, two of which have previously been implicated as critical for Shigella growth, but which are dispensable for E. coli growth.

CONCLUSIONS

The data presented here highlight the small number of protein coding genes for which we have strong evidence that their essentiality status differs between the closely related bacterial taxa E. coli and Shigella. A set of genes involved in acetate utilization provides a canonical example. These results leave open the possibility of developing strain-specific antibiotic treatments targeting such differentially essential genes, but suggest that such opportunities may be rare in closely related bacteria.

A simple screen to identify promoters conferring high levels of phenotypic noise

Genetically identical populations of unicellular organisms often show marked variation in some phenotypic traits. To investigate the molecular causes and possible biological functions of this phenotypic noise, it would be useful to have a method to identify genes whose expression varies stochastically on a certain time scale. Here, we developed such a method and used it for identifying genes with high levels of phenotypic noise in Salmonella enterica ssp. I serovar Typhimurium (S. Typhimurium). We created a genomic plasmid library fused to a green fluorescent protein (GFP) reporter and subjected replicate populations harboring this library to fluctuating selection for GFP expression using fluorescent-activated cell sorting (FACS). After seven rounds of fluctuating selection, the populations were strongly enriched for promoters that showed a high amount of noise in gene expression. Our results indicate that the activity of some promoters of S. Typhimurium varies on such a short time scale that these promoters can absorb rapid fluctuations in the direction of selection, as imposed during our experiment. The genomic fragments that conferred the highest levels of phenotypic variation were promoters controlling the synthesis of flagella, which are associated with virulence and host–pathogen interactions. This confirms earlier reports that phenotypic noise may play a role in pathogenesis and indicates that these promoters have among the highest levels of noise in the S. Typhimurium genome. This approach can be applied to many other bacterial and eukaryotic systems as a simple method for identifying genes with noisy expression.

According to the conventional view, the characteristics of an organism are determined by nature and nurture—by its genes and by the environment it lives in. Consequently, one would expect that two organisms that share the same genes and live in the same environment have identical characteristics. Recently it has become clear that this expectation is often not borne out; clonal families of simple organisms living under constant conditions often show variation in biological traits and sometimes even have markedly different properties and do different things. In order to investigate molecular causes and possible biological functions of such phenotypic noise, it would be very valuable to have a simple and fast method for identifying biological traits that are particularly noisy. Here, we developed such a method, and screened for noisy traits in the bacterial pathogen Salmonella typhimurium that vary at a time scale of one day. We found that traits involved in interaction with the host are particularly noisy, suggesting that phenotypic noise might be important in pathogenesis. This method can be readily adopted for other organisms and might contribute to elucidating the role of noise in biology.

Diagnostic discrimination of live attenuated influenza vaccine strains and community-acquired pathogenic strains in clinical samples

Live vaccines can generate false-positive results on common influenza assays including reverse transcriptase-PCR (RT-PCR), culture and antigen tests. This threatens the integrity of epidemiological data and may misdirect treatment and control efforts. We report the development of RT-PCR tests that distinguish live FluMist vaccine (FMV) strains from circulating influenza strains in clinical samples. Primers were validated using influenza-positive samples from unvaccinated patients, packaged FMV, and one PCR-positive asymptomatic vaccine. Furthermore, the assay was used to experimentally test our lab's collection of influenza-positive samples from the 2004-05 and 2005-06 influenza seasons and several 2005 preseason isolates to determine the rate of vaccine-derived false-positive results under differing epidemiological conditions. Analytical and clinical validations show that the assay is both sensitive and specific. Experimental results demonstrate that 51 out of 51 influenza-positive samples collected during influenza season from ill, previously-vaccinated military personnel represent real infections with circulating strains. Finally, the assay shows that four preseason influenza-positive samples were false positives resulting from vaccine shedding. The vaccine-discriminatory RT-PCR methods described here provide the first test designed to distinguish FMV strains from circulating strains. The results show that the test is effective, and demonstrate the importance of such tests in the age of live vaccines.

Transmission dynamics and prospective environmental sampling of adenovirus in a military recruit setting

BackgroundHigh levels of morbidity caused by adenovirus among US military recruits have returned since the loss of adenovirus vaccines in 1999. The transmission dynamics of adenovirus have never been well understood, which complicates prevention efforts

MethodsEnrollment and end-of-study samples were obtained and active surveillance for febrile respiratory illnesses (FRIs) was performed for 341 recruits and support personnel. Environmental samples were collected simultaneously. Classic and advanced diagnostic techniques were used

ResultsSeventy-nine percent (213/271) of new recruits were seronegative for either adenovirus serotype 4 (Ad-4) or adenovirus serotype 7 (Ad-7). FRI caused by Ad-4 was observed in 25% (67/271) of enrolled recruits, with 100% of them occurring in individuals with enrollment titers <1:4. The percentage of recruits seropositive for Ad-4 increased from 34% at enrollment to 97% by the end of the study. Adenovirus was most commonly detected in the environment on pillows, lockers, and rifles

ConclusionsPotential sources of adenovirus transmission among US military recruits included the presence of adenovirus on surfaces in living quarters and extended pharyngeal viral shedding over the course of several days. The introduction of new recruits, who were still shedding adenovirus, into new training groups was documented. Serological screening could identify susceptible recruits for the optimal use of available vaccines. New high-throughput technologies show promise in providing valuable data for clinical and research applications

Surveillance for febrile respiratory infections during Cobra Gold 2003

The Naval Health Research Center conducted laboratory-based surveillance for febrile respiratory infections at the 2003 Cobra Gold Exercise in Thailand. Seventeen individuals met the case definition for febrile respiratory illness, and diagnostic specimens were obtained from 16. Laboratory testing identified influenza A for 44%; sequence analysis demonstrated that these were Fujian-like influenza strains, which represented the predominant strain found globally in 2003/2004. Other pathogens identified included coronavirus OC43, respiratory syncytial virus, and rhinovirus. Logistical challenges were overcome as laboratory-supported febrile respiratory illness surveillance was conducted during a military training exercise. With heightened concern over the potential for another global influenza pandemic, such surveillance could prove critical for the detection of emerging influenza and respiratory pathogen strains with potential for importation to the United States.

Co-infections of adenovirus species in previously vaccinated patients

Adenoviral infections associated with respiratory illness in military trainees involve multiple co-infecting species and serotypes.

Despite the success of the adenovirus vaccine administered to US military trainees, acute respiratory disease (ARD) surveillance still detected breakthrough infections (respiratory illnesses associated with the adenovirus serotypes specifically targeted by the vaccine). To explore the role of adenoviral co-infection (simultaneous infection by multiple pathogenic adenovirus species) in breakthrough disease, we examined specimens from patients with ARD by using 3 methods to detect multiple adenoviral species: a DNA microarray, a polymerase chain reaction­ (PCR)–enzyme-linked immunosorbent assay, and a multiplex PCR assay. Analysis of 52 samples (21 vaccinated, 31 unvaccinated) collected from 1996 to 2000 showed that all vaccinated samples had co-infections. Most of these co-infections were community-acquired serotypes of species B1 and E. Unvaccinated samples primarily contained only 1 species (species E) associated with adult respiratory illness. This study highlights the rarely reported phenomenon of adenoviral co-infections in a clinically relevant environment suitable for the generation of new recombinational variants.

Effectiveness of the 2003–2004 influenza vaccine among U.S. military basic trainees: a year of suboptimal match between vaccine and circulating strain

Effectiveness of the 2003-2004 influenza vaccine was evaluated at five military basic training centers throughout the United States. Data from surveillance conducted in December and January 2003-2004 in this highly vaccinated population were evaluated. During this period, 10.6% (37/350) of specimens were positive for influenza A. A 14-day period after vaccination was considered the period prior to immune protection; vaccine effectiveness (VE) was calculated based on febrile respiratory illness presentation and laboratory confirmation of influenza before or after this 14-day period. Thirty-two cases presented within 14 days of vaccination, and five cases presented beyond 14 days from vaccination. VE in this population was estimated to be 94.4% for laboratory-confirmed influenza. In contrast, VE was only 13.9% for influenza-like illness (ILI) without a laboratory confirmation.

Molecular analysis of adenovirus isolates from vaccinated and unvaccinated young adults

Infections of adenovirus type 4 (Ad4) and Ad7 were discovered among previously vaccinated individuals through febrile respiratory illness surveillance at military recruit camps. Genetic analysis was performed on these isolates and a sample of adenovirus isolates from unvaccinated patients. Antigenic regions of the adenovirus hexon gene from 21 vaccinated and 31 unvaccinated patients were sequenced and compared to homologous regions of Ad4 and Ad7 vaccine strains and of other representative hexon sequences archived in GenBank. The phylogenetic distribution of sequences from vaccinated individuals closely resembled those from unvaccinated individuals. The most common Ad7 strain was the Ad7d2 hexon genotype, and the most common Ad4 strain was a genotype nearly identical to the recently discovered Z-G 95-873 Ad4 variant. Near exclusive isolation of Ad4 since 1999 indicates that the Ad4 variant is currently responsible for the vast majority of adenovirus morbidity in military recruit camps. Different ratios of nonsynonymous to synonymous nucleotide substitution rates in known antigenic regions compared to nonantigenic regions indicated positive selection for diversity in the antigenic regions and purifying selection in the nonantigenic regions.