NucAmino: a nucleotide to amino acid alignment optimized for virus gene sequences

CNCA aligns small annotated genomes

BACKGROUND

To explore the evolutionary history of sequences, a sequence alignment is a first and necessary step, and its quality is crucial. In the context of the study of the proximal origins of SARS-CoV-2 coronavirus, we wanted to construct an alignment of genomes closely related to SARS-CoV-2 using both coding and non-coding sequences. To our knowledge, there is no tool that can be used to construct this type of alignment, which motivated the creation of CNCA.

RESULTS

CNCA is a web tool that aligns annotated genomes from GenBank files. It generates a nucleotide alignment that is then updated based on the protein sequence alignment. The output final nucleotide alignment matches the protein alignment and guarantees no frameshift. CNCA was designed to align closely related small genome sequences up to 50 kb (typically viruses) for which the gene order is conserved.

CONCLUSIONS

CNCA constructs multiple alignments of small genomes by integrating both coding and non-coding sequences. This preserves regions traditionally ignored in conventional back-translation methods, such as non-coding regions.

Development and validation of a next-generation sequencing assay with open-access analysis software for detecting resistance-associated mutations in CMV

Cytomegalovirus (CMV) resistance testing by targeted next-generation sequencing (NGS) allows for the simultaneous analysis of multiple genes. We developed and validated an amplicon-based Ion Torrent NGS assay to detect CMV resistance mutations in UL27, UL54, UL56, and UL97 and compared the results to standard Sanger sequencing. NGS primers were designed to generate 83 overlapping amplicons of four CMV genes (~10 kb encompassing 138 mutation sites). An open-access software plugin was developed to perform read alignment, call variants, and interpret drug resistance. Plasmids were tested to determine NGS error rate and minor variant limit of detection. NGS limit of detection was determined using the CMV WHO International Standard and quantified clinical specimens. Reproducibility was also assessed. After establishing quality control metrics, 185 patient specimens previously tested using Sanger were reanalyzed by NGS. The NGS assay had a low error rate (<0.05%) and high accuracy (95%) for detecting CMV-associated resistance mutations present at ≥5% in contrived mixed populations. Mutation sites were reproducibly sequenced with 40× coverage when plasma viral loads were ≥2.6 log IU/mL. NGS detected the same resistance-associated mutations identified by Sanger in 68/69 (98.6%) specimens. In 16 specimens, NGS detected 18 resistance mutations that Sanger failed to detect; 14 were low-frequency variants (<20%), and six would have changed the drug resistance interpretation. The NGS assay showed excellent agreement with Sanger and generated high-quality sequence from low viral load specimens. Additionally, the higher resolution and analytic sensitivity of NGS potentially enables earlier detection of antiviral resistance.

Improved DNA-Versus-Protein Homology Search for Protein Fossils

Protein fossils, i.e., noncoding DNA descended from coding DNA, arise frequently from transposable elements (TEs), decayed genes, and viral integrations. They can reveal, and mislead about, evolutionary history and relationships. They have been detected by comparing DNA to protein sequences, but current methods are not optimized for this task. We describe a powerful DNA-protein homology search method. We use a 64×21 substitution matrix, which is fitted to sequence data, automatically learning the genetic code. We detect subtly homologous regions by considering alternative possible alignments between them, and calculate significance (probability of occurring by chance between random sequences). Our method detects TE protein fossils much more sensitively than blastx, and faster. Of the  ∼ 7 major categories of eukaryotic TE, three were long thought absent in mammals: we find two of them in the human genome, polinton and DIRS/Ngaro. This method increases our power to find ancient fossils, and perhaps to detect non-standard genetic codes. The alternative-alignments and significance paradigm is not specific to DNA-protein comparison, and could benefit homology search generally. This is an extended version of a conference paper (Yao & Frith, 2021).

Characteristics of hepatitis C virus resistance in an international cohort after a decade of direct-acting antivirals

Background & Aims

Direct-acting antiviral (DAA) regimens provide a cure in >95% of patients with chronic HCV infection. However, in some patients in whom therapy fails, resistance-associated substitutions (RASs) can develop, limiting retreatment options and risking onward resistant virus transmission. In this study, we evaluated RAS prevalence and distribution, including novel NS5A RASs and clinical factors associated with RAS selection, among patients who experienced DAA treatment failure.

Methods

SHARED is an international consortium of clinicians and scientists studying HCV drug resistance. HCV sequence linked metadata from 3,355 patients were collected from 22 countries. NS3, NS5A, and NS5B RASs in virologic failures, including novel NS5A substitutions, were examined. Associations of clinical and demographic characteristics with RAS selection were investigated.

Results

The frequency of RASs increased from its natural prevalence following DAA exposure: 37% to 60% in NS3, 29% to 80% in NS5A, 15% to 22% in NS5B for sofosbuvir, and 24% to 37% in NS5B for dasabuvir. Among 730 virologic failures, most were treated with first-generation DAAs, 94% had drug resistance in ≥1 DAA class: 31% single-class resistance, 42% dual-class resistance (predominantly against protease and NS5A inhibitors), and 21% triple-class resistance. Distinct patterns containing ≥2 highly resistant RASs were common. New potential NS5A RASs and adaptive changes were identified in genotypes 1a, 3, and 4. Following DAA failure, RAS selection was more frequent in older people with cirrhosis and those infected with genotypes 1b and 4.

Conclusions

Drug resistance in HCV is frequent after DAA treatment failure. Previously unrecognized substitutions continue to emerge and remain uncharacterized.

Lay summary

Although direct-acting antiviral medications effectively cure hepatitis C in most patients, sometimes treatment selects for resistant viruses, causing antiviral drugs to be either ineffective or only partially effective. Multidrug resistance is common in patients for whom DAA treatment fails. Older patients and patients with advanced liver diseases are more likely to select drug-resistant viruses. Collective efforts from international communities and governments are needed to develop an optimal approach to managing drug resistance and preventing the transmission of resistant viruses.

RNA genome conservation and secondary structure in SARS-CoV-2 and SARS-related viruses: a first look

As the COVID-19 outbreak spreads, there is a growing need for a compilation of conserved RNA genome regions in the SARS-CoV-2 virus along with their structural propensities to guide development of antivirals and diagnostics. Here we present a first look at RNA sequence conservation and structural propensities in the SARS-CoV-2 genome. Using sequence alignments spanning a range of betacoronaviruses, we rank genomic regions by RNA sequence conservation, identifying 79 regions of length at least 15 nt as exactly conserved over SARS-related complete genome sequences available near the beginning of the COVID-19 outbreak. We then confirm the conservation of the majority of these genome regions across 739 SARS-CoV-2 sequences subsequently reported from the COVID-19 outbreak, and we present a curated list of 30 "SARS-related-conserved" regions. We find that known RNA structured elements curated as Rfam families and in prior literature are enriched in these conserved genome regions, and we predict additional conserved, stable secondary structures across the viral genome. We provide 106 "SARS-CoV-2-conserved-structured" regions as potential targets for antivirals that bind to structured RNA. We further provide detailed secondary structure models for the extended 5' UTR, frameshifting stimulation element, and 3' UTR. Lastly, we predict regions of the SARS-CoV-2 viral genome that have low propensity for RNA secondary structure and are conserved within SARS-CoV-2 strains. These 59 "SARS-CoV-2-conserved-unstructured" genomic regions may be most easily accessible by hybridization in primer-based diagnostic strategies.

VIRULIGN: fast codon-correct alignment and annotation of viral genomes

Summary

Virus sequence data are an essential resource for reconstructing spatiotemporal dynamics of viral spread as well as to inform treatment and prevention strategies. However, the potential benefit of these applications critically depends on accurate and correctly annotated alignments of genetically heterogeneous data. VIRULIGN was built for fast codon-correct alignments of large datasets, with standardized and formalized genome annotation and various alignment export formats.

Availability and implementation

VIRULIGN is freely available at https://github.com/rega-cev/virulign as an open source software project.

Supplementary information

Supplementary data is available at Bioinformatics online.