SARS-CoV-2 Delta variant isolates from vaccinated individuals.

BACKGROUND: The SARS-CoV-2 Delta variant was first identified in the U.S. in March 2021 and has rapidly become the predominant lineage across the U.S. due to increased transmissibility, immune evasion and vaccine breakthrough. The aim of this study was to better understand the genetic diversity and the potential impact of mutations observed in SARS-CoV-2 viruses circulating in the U.S. in vaccinated individuals. RESULTS: Whole genome sequencing was performed on thirty-four SARS-CoV-2 positive samples using the Oxford Nanopore MinION. Evolutionary genomic analysis revealed two novel mutations, ORF1b:V2354F and a premature stop codon, ORF7a:Q94*, identified in a cluster of SARS-CoV-2 Delta isolates collected from vaccinated individuals in Colorado. The ORF1b:V2354F mutation, corresponding to NSP15:V303F, may induce a conformational change and result in a disruption to a flanking beta-sheet structure. The premature stop codon, ORF7a:Q94*, truncates the transmembrane protein and cytosolic tail used to mediate protein transport. This may affect protein localization to the ER-Golgi. In addition to these novel mutations, the cluster of vaccinated isolates contain an additional mutation in the spike protein, at position 112, compared to the Delta variant defining mutations. This mutation, S112L, exists in isolates previously obtained in the U.S. The S112L mutation substitutes a bulky hydrophobic side chain for a polar side chain, which results in a non-conservative substitution within the protein that may affect antibody-binding affinity. Additionally, the vaccinated cluster of isolates contains non-synonymous mutations within ORF8 and NSPs which further distinguish this cluster from the respective ancestral Delta variant. CONCLUSIONS: These results show there is an emerging sub-lineage of the ancestral Delta variant circulating in the U.S. As mutations emerge in constellations, those with a potentially beneficial advantage to the virus may continue to circulate while others will cease.

Clinical evaluation of a multiplex real-time RT-PCR assay for detection of SARS-CoV-2 in individual and pooled upper respiratory tract samples.

The aim of this study was to identify and validate a sensitive, high-throughput, and cost-effective SARS-CoV-2 real-time RT-PCR assay to be used as a surveillance and diagnostic tool for SARS-CoV-2 in a university surveillance program. We conducted a side-by-side clinical evaluation of a newly developed SARS-CoV-2 multiplex assay (EZ-SARS-CoV-2 Real-Time RT-PCR) with the commercial TaqPath COVID-19 Combo Kit, which has an Emergency Use Authorization from the FDA. The EZ-SARS-CoV-2 RT-PCR incorporates two assays targeting the SARS-CoV-2 N gene, an internal control targeting the human RNase P gene, and a PCR inhibition control in a single reaction. Nasopharyngeal (NP) and anterior nares (AN) swabs were tested as individuals and pools with both assays and in the ABI 7500 Fast and the QuantStudio 5 detection platforms. The analytical sensitivity of the EZ-SARS-CoV-2 RT-PCR assay was 250 copies/ml or approximately 1.75 genome copy equivalents per reaction. The clinical performance of the EZ-SARS-CoV-2 assay was evaluated using NP and AN samples tested in other laboratories. The diagnostic sensitivity of the assay ranged between 94 and 96% across the detection platforms, and the diagnostic specificity was 94.06%. The positive predictive value was 94%, and the negative predictive value ranged from 94 to 96%. Pooling five NP or AN specimens yielded 93% diagnostic sensitivity. The overall agreement between these SARS-CoV-2 RT-PCR assays was high, supported by a Cohen's kappa value of 0.93. The EZ-SARS-CoV-2 RT-PCR assay performance attributes of high sensitivity and specificity with AN sample matrix and pooled upper respiratory samples support its use in a high-throughput surveillance testing program.

Analyzing Tumor Gene Expression Factors with the CorExplorer Web Portal.

Differential gene expression analysis is an important technique for understanding disease states. The machine learning algorithm CorEx has shown utility in analyzing differential expression of groups of genes in tumor RNA-seq in a way that may be helpful for advancing precision oncology. However, CorEx produces many factors that can be challenging to analyze and connect to existing understanding. To facilitate such connections, we have built a website, CorExplorer, that allows users to interactively explore the data and answer common questions related to its analysis. We trained CorEx on RNA-seq gene expression data for four tumor types: ovarian, lung, melanoma, and colorectal. We then incorporated corresponding survival, protein-protein interactions, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichments, and heatmaps into the website for association with the factor graph visualization. Here we employ example protocols to illustrate use of the database for comprehending the significance of the learned tumor factors in the context of this external data.

Overcoming the challenges of pen-side molecular diagnosis of African swine fever to support outbreak investigations under field conditions.

African swine fever (ASF) is a devastating disease of pigs. Without a vaccine, early detection and rapid diagnosis of ASF is a crucial step towards effective disease control. In many countries where ASF is endemic, laboratory infrastructure including sampling and sample shipment is inadequate, and a rapid laboratory confirmation would require that the diagnosis is performed at regional laboratories close to the pig farms of concern, or even at the farm-side. This study intended to evaluate measures including sample preparation methods, a dried-down assay, and a portable, battery-powered real-time PCR instrument, to improve molecular diagnosis under field conditions. A simple dilution of blood samples, either in Phosphate-buffered saline or a commercial buffer, worked similarly to beads-based nucleic acid extraction using a magnet as the core equipment; the latter method did work as well for those samples with low viral load or high Ct values. The real-time PCR assay using a Universal ProbeLibrary (UPL) probe tolerated suspected inhibitory substances present in the prepared samples better, whereas the dried-down assay had a higher diagnostic sensitivity. Additionally, an inhibition control assay proved to be helpful in avoiding false negative results when interpreting negative results of samples that might be of low quality or with inadequate reduction in inhibitory substances. When tested with synthetic DNA standards, the portable instrument performed at a level approaching stationary thermocyclers. In summary, the developments of suitable sample preparation methods, robust and thermal-stable real-time PCR assays with inhibition control, and battery-powered portable thermocyclers with middle-throughput offer one way forward to provide rapid, reliable molecular diagnosis under challenging field conditions.

Clinical, Serological, and Molecular Observations from a Case Series Study during the Asian Lineage Zika Virus Outbreak in Grenada during 2016.

This paper describes the spatial and temporal distribution of cases, demographic characteristics of patients, and clinical manifestations of Zika virus (ZIKV) during the 2016 outbreak in Grenada. The first reported case was recorded in St. Andrew Parish in April, and the last reported case was seen in November, with peak transmission occurring in the last week of June, based on test results. Data were collected from a total of 514 patients, of whom 207 (40%) tested positive for ZIKV. No evidence was found that testing positive for ZIKV infection was related to age, gender, or pregnancy status. Clinical presentation with rash (OR = 2.4, 95% CI = 1.5 to 3.7) or with lymphadenopathy (OR = 1.7, 95% CI = 1.0 to 2.9) were the only reported symptoms consistent with testing positive for ZIKV infection. During the Zika outbreak, the infection rate was 20 clinical cases per 10,000 in the population compared to 41 cases per 10,000 during the chikungunya outbreak in Grenada in 2014 and 17 cases per 10,000 during the dengue outbreak in 2001-2002. Even though the country has employed vector control programs, with no apparent decrease in infection rates, it appears that new abatement approaches are needed to minimize morbidity in future arbovirus outbreaks.

Allele Workbench: transcriptome pipeline and interactive graphics for allele-specific expression.

Sequencing the transcriptome can answer various questions such as determining the transcripts expressed in a given species for a specific tissue or condition, evaluating differential expression, discovering variants, and evaluating allele-specific expression. Differential expression evaluates the expression differences between different strains, tissues, and conditions. Allele-specific expression evaluates expression differences between parental alleles. Both differential expression and allele-specific expression have been studied for heterosis (hybrid vigor), where the hybrid has improved performance over the parents for one or more traits. The Allele Workbench software was developed for a heterosis study that evaluated allele-specific expression for a mouse F1 hybrid using libraries from multiple tissues with biological replicates. This software has been made into a distributable package, which includes a pipeline, a Java interface to build the database, and a Java interface for query and display of the results. The required input is a reference genome, annotation file, and one or more RNA-Seq libraries with optional replicates. It evaluates allelic imbalance at the SNP and transcript level and flags transcripts with significant opposite directional allele-specific expression. The Java interface allows the user to view data from libraries, replicates, genes, transcripts, exons, and variants, including queries on allele imbalance for selected libraries. To determine the impact of allele-specific SNPs on protein folding, variants are annotated with their effect (e.g., missense), and the parental protein sequences may be exported for protein folding analysis. The Allele Workbench processing results in transcript files and read counts that can be used as input to the previously published Transcriptome Computational Workbench, which has a new algorithm for determining a trimmed set of gene ontology terms. The software with demo files is available from https://code.google.com/p/allele-workbench. Additionally, all software is ready for immediate use from an Atmosphere Virtual Machine Image available from the iPlant Collaborative (www.iplantcollaborative.org).

Making a living while starving in the dark: metagenomic insights into the energy dynamics of a carbonate cave.

Carbonate caves represent subterranean ecosystems that are largely devoid of phototrophic primary production. In semiarid and arid regions, allochthonous organic carbon inputs entering caves with vadose-zone drip water are minimal, creating highly oligotrophic conditions; however, past research indicates that carbonate speleothem surfaces in these caves support diverse, predominantly heterotrophic prokaryotic communities. The current study applied a metagenomic approach to elucidate the community structure and potential energy dynamics of microbial communities, colonizing speleothem surfaces in Kartchner Caverns, a carbonate cave in semiarid, southeastern Arizona, USA. Manual inspection of a speleothem metagenome revealed a community genetically adapted to low-nutrient conditions with indications that a nitrogen-based primary production strategy is probable, including contributions from both Archaea and Bacteria. Genes for all six known CO2-fixation pathways were detected in the metagenome and RuBisCo genes representative of the Calvin-Benson-Bassham cycle were over-represented in Kartchner speleothem metagenomes relative to bulk soil, rhizosphere soil and deep-ocean communities. Intriguingly, quantitative PCR found Archaea to be significantly more abundant in the cave communities than in soils above the cave. MEtaGenome ANalyzer (MEGAN) analysis of speleothem metagenome sequence reads found Thaumarchaeota to be the third most abundant phylum in the community, and identified taxonomic associations to this phylum for indicator genes representative of multiple CO2-fixation pathways. The results revealed that this oligotrophic subterranean environment supports a unique chemoautotrophic microbial community with potentially novel nutrient cycling strategies. These strategies may provide key insights into other ecosystems dominated by oligotrophy, including aphotic subsurface soils or aquifers and photic systems such as arid deserts.

Portable nucleic acid thermocyclers.

A nucleic acid thermal cycler is considered to be portable if it is under ten pounds, easily carried by one individual, and battery powered. Nucleic acid amplification includes both polymerase chain reaction (e.g. PCR, RT-PCR) and isothermal amplification (e.g. RPA, HDA, LAMP, NASBA, RCA, ICAN, SMART, SDA). There are valuable applications for portable nucleic acid thermocyclers in fields that include clinical diagnostics, biothreat detection, and veterinary testing. A system that is portable allows for the distributed detection of targets at the point of care and a reduction of the time from sample to answer. The designer of a portable nucleic acid thermocycler must carefully consider both thermal control and the detection of amplification. In addition to thermal control and detection, the designer may consider the integration of a sample preparation subsystem with the nucleic acid thermocycler. There are a variety of technologies that can achieve accurate thermal control and the detection of nucleic acid amplification. Important evaluation criteria for each technology include maturity, power requirements, cost, sensitivity, speed, and manufacturability. Ultimately the needs of a particular market will lead to user requirements that drive the decision between available technologies.

Profiling bacterial diversity and taxonomic composition on speleothem surfaces in Kartchner Caverns, AZ.

Caves are relatively accessible subterranean habitats ideal for the study of subsurface microbial dynamics and metabolisms under oligotrophic, non-photosynthetic conditions. A 454-pyrotag analysis of the V6 region of the 16S rRNA gene was used to systematically evaluate the bacterial diversity of ten cave surfaces within Kartchner Caverns, a limestone cave. Results showed an average of 1,994 operational taxonomic units (97 % cutoff) per speleothem and a broad taxonomic diversity that included 21 phyla and 12 candidate phyla. Comparative analysis of speleothems within a single room of the cave revealed three distinct bacterial taxonomic profiles dominated by either Actinobacteria, Proteobacteria, or Acidobacteria. A gradient in observed species richness along the sampling transect revealed that the communities with lower diversity corresponded to those dominated by Actinobacteria while the more diverse communities were those dominated by Proteobacteria. A 16S rRNA gene clone library from one of the Actinobacteria-dominated speleothems identified clones with 99 % identity to chemoautotrophs and previously characterized oligotrophs, providing insights into potential energy dynamics supporting these communities. The robust analysis conducted for this study demonstrated a rich bacterial diversity on speleothem surfaces. Further, it was shown that seemingly comparable speleothems supported divergent phylogenetic profiles suggesting that these communities are very sensitive to subtle variations in nutritional inputs and environmental factors typifying speleothem surfaces in Kartchner Caverns.

Life at the hyperarid margin: novel bacterial diversity in arid soils of the Atacama Desert, Chile.

Nearly half the earth's surface is occupied by dryland ecosystems, regions susceptible to reduced states of biological productivity caused by climate fluctuations. Of these regions, arid zones located at the interface between vegetated semiarid regions and biologically unproductive hyperarid zones are considered most vulnerable. The objective of this study was to conduct a deep diversity analysis of bacterial communities in unvegetated arid soils of the Atacama Desert, to characterize community structure and infer the functional potential of these communities based on observed phylogenetic associations. A 454-pyrotag analysis was conducted of three unvegetated arid sites located at the hyperarid-arid margin. The analysis revealed communities with unique bacterial diversity marked by high abundances of novel Actinobacteria and Chloroflexi and low levels of Acidobacteria and Proteobacteria, phyla that are dominant in many biomes. A 16S rRNA gene library of one site revealed the presence of clones with phylogenetic associations to chemoautotrophic taxa able to obtain energy through oxidation of nitrite, carbon monoxide, iron, or sulfur. Thus, soils at the hyperarid margin were found to harbor a wealth of novel bacteria and to support potentially viable communities with phylogenetic associations to non-phototrophic primary producers and bacteria capable of biogeochemical cycling.

SyMAP v3.4: a turnkey synteny system with application to plant genomes.

SyMAP (Synteny Mapping and Analysis Program) was originally developed to compute synteny blocks between a sequenced genome and a FPC map, and has been extended to support pairs of sequenced genomes. SyMAP uses MUMmer to compute the raw hits between the two genomes, which are then clustered and filtered using the optional gene annotation. The filtered hits are input to the synteny algorithm, which was designed to discover duplicated regions and form larger-scale synteny blocks, where intervening micro-rearrangements are allowed. SyMAP provides extensive interactive Java displays at all levels of resolution along with simultaneous displays of multiple aligned pairs. The synteny blocks from multiple chromosomes may be displayed in a high-level dot plot or three-dimensional view, and the user may then drill down to see the details of a region, including the alignments of the hits to the gene annotation. These capabilities are illustrated by showing their application to the study of genome duplication, differential gene loss and transitive homology between sorghum, maize and rice. The software may be used from a website or standalone for the best performance. A project manager is provided to organize and automate the analysis of multi-genome groups. The software is freely distributed at http://www.agcol.arizona.edu/software/symap.

A BAC-based physical map of Brachypodium distachyon and its comparative analysis with rice and wheat.

BACKGROUND: Brachypodium distachyon (Brachypodium) has been recognized as a new model species for comparative and functional genomics of cereal and bioenergy crops because it possesses many biological attributes desirable in a model, such as a small genome size, short stature, self-pollinating habit, and short generation cycle. To maximize the utility of Brachypodium as a model for basic and applied research it is necessary to develop genomic resources for it. A BAC-based physical map is one of them. A physical map will facilitate analysis of genome structure, comparative genomics, and assembly of the entire genome sequence. RESULTS: A total of 67,151 Brachypodium BAC clones were fingerprinted with the SNaPshot HICF fingerprinting method and a genome-wide physical map of the Brachypodium genome was constructed. The map consisted of 671 contigs and 2,161 clones remained as singletons. The contigs and singletons spanned 414 Mb. A total of 13,970 gene-related sequences were detected in the BAC end sequences (BES). These gene tags aligned 345 contigs with 336 Mb of rice genome sequence, showing that Brachypodium and rice genomes are generally highly colinear. Divergent regions were mainly in the rice centromeric regions. A dot-plot of Brachypodium contigs against the rice genome sequences revealed remnants of the whole-genome duplication caused by paleotetraploidy, which were previously found in rice and sorghum. Brachypodium contigs were anchored to the wheat deletion bin maps with the BES gene-tags, opening the door to Brachypodium-Triticeae comparative genomics. CONCLUSION: The construction of the Brachypodium physical map, and its comparison with the rice genome sequence demonstrated the utility of the SNaPshot-HICF method in the construction of BAC-based physical maps. The map represents an important genomic resource for the completion of Brachypodium genome sequence and grass comparative genomics. A draft of the physical map and its comparisons with rice and wheat are available at http://phymap.ucdavis.edu/brachypodium/.

Efficacy of clone fingerprinting methodologies.

With the development of new high-information content fingerprinting techniques for constructing BAC-based physical maps, physical map construction is accelerating and it is important to determine which methodologies work best. In a recent publication (Z. Xu et al., 2004, Genomics 84:941-951), Xu et al. evaluated five different techniques (one agarose-based and four using multiple enzymes) and concluded that a two-enzyme technique was superior. In addition, they found that no benefit was gained from fingerprinting more than 10x coverage. In this paper we report our own extensive simulation results, which lead to contrasting conclusions. Our data indicate that the five-enzyme method known as SNaPshot is the most effective and that the assembly can in fact be significantly improved with greater than 10x coverage.

Detection of foot-and-mouth disease virus: comparative diagnostic sensitivity of two independent real-time reverse transcription-polymerase chain reaction assays.

Rapid and accurate diagnosis is central to the effective control of foot-and-mouth disease (FMD). It is now recognized that reverse-transcription polymerase chain reaction (RT-PCR) assays can play an important role in the routine detection of FMD virus (FMDV) in clinical samples. The aim of this study was to compare the ability of 2 independent real-time RT-PCR (rRT-PCR) assays targeting the 5' untranslated region (5'UTR) and RNA polymerase (3D) to detect FMDV in clinical samples. There was concordance between the results generated by the 2 assays for 88.1% (347 of 394) of RNA samples extracted from suspensions of epithelial tissue obtained from suspect FMD cases. The comparison between the 2 tests highlighted 19 FMDV isolates (13 for the 5'UTR and 6 for the 3D assay), which failed to produce a signal in 1 assay but gave a positive signal in the other. The sequence of the genomic targets of selected isolates highlighted nucleotide substitutions in the primer or probe regions, thereby providing an explanation for negative results generated in the rRT-PCR assays. These data illustrate the importance of the continuous monitoring of circulating FMDV field strains to ensure the design of the rRT-PCR assay remains fit for purpose and suggest that the use of multiple diagnostic targets could further enhance the sensitivity of molecular methods for the detection of FMDV.

Whole-genome validation of high-information-content fingerprinting.

Fluorescent-based high-information-content fingerprinting (HICF) techniques have recently been developed for physical mapping. These techniques make use of automated capillary DNA sequencing instruments to enable both high-resolution and high-throughput fingerprinting. In this article, we report the construction of a whole-genome HICF FPC map for maize (Zea mays subsp. mays cv B73), using a variant of HICF in which a type IIS restriction enzyme is used to generate the fluorescently labeled fragments. The HICF maize map was constructed from the same three maize bacterial artificial chromosome libraries as previously used for the whole-genome agarose FPC map, providing a unique opportunity for direct comparison of the agarose and HICF methods; as a result, it was found that HICF has substantially greater sensitivity in forming contigs. An improved assembly procedure is also described that uses automatic end-merging of contigs to reduce the effects of contamination and repetitive bands. Several new features in FPC v7.2 are presented, including shared-memory multiprocessing, which allows dramatically faster assemblies, and automatic end-merging, which permits more accurate assemblies. It is further shown that sequenced clones may be digested in silico and located accurately on the HICF assembly, despite size deviations that prevent the precise prediction of experimental fingerprints. Finally, repetitive bands are isolated, and their effect on the assembly is studied.

Use of a portable real-time reverse transcriptase-polymerase chain reaction assay for rapid detection of foot-and-mouth disease virus.

OBJECTIVE: To evaluate a portable real-time reverse transcriptase-polymerase chain reaction (RT-PCR) assay designed to detect all 7 viral serotypes of foot-and-mouth disease virus (FMDV). DESIGN: Laboratory and animal studies. STUDY POPULATION: Viruses grown in tissue culture and animals experimentally infected with FMDV. PROCEDURE: 1 steer, pig, and sheep were infected with serotype O FMDV. Twenty-four hours later, animals were placed in separate rooms that contained 4 FMDV-free, healthy animals of the same species. Oral and nasal swab specimens, oropharyngeal specimens obtained with a probang, and blood samples were obtained at frequent intervals, and animals were observed for fever and clinical signs of foot-and-mouth disease (FMD). Samples from animals and tissue cultures were assayed for infectious virus and viral RNA. RESULTS: The assay detected viral RNA representing all 7 FMDV serotypes grown in tissue culture but did not amplify a panel of selected viruses that included those that cause vesicular diseases similar to FMD; thus, the assay had a specificity of 100%, depending on the panel selected. The assay also met or exceeded sensitivity of viral culture on samples from experimentally infected animals. In many instances, the assay detected viral RNA in the mouth and nose 24 to 96 hours before the onset of clinical disease. CONCLUSIONS AND CLINICAL RELEVANCE: The assay reagents are produced in a vitrified form, which permits storage and transportation at ambient temperatures. The test can be performed in 2 hours or less on a portable instrument, thus providing a rapid, portable, sensitive, and specific method for detection of FMDV.