Full-length Genome Characterization Research Articles

Full-length genome characterization and phylogenetic analysis of SARS-CoV-2 virus strains from Yogyakarta and Central Java, Indonesia.

BackgroundRecently, SARS-CoV-2 virus with the D614G mutation has become a public concern due to rapid dissemination of this variant across many countries. Our study aims were (1) to report full-length genome sequences of SARS-CoV-2 collected from four COVID-19 patients in the Special Region of Yogyakarta and Central Java provinces, Indonesia; (2) to compare the clade distribution of full-length genome sequences from Indonesia (n = 60) from March to September 2020 and (3) to perform phylogenetic analysis of SARS-CoV-2 complete genomes from different countries, including Indonesia.MethodsWhole genome sequencing (WGS) was performed using next-generation sequencing (NGS) applied in the Illumina MiSeq instrument. Full-length virus genomes were annotated using the reference genome of hCoV-19/Wuhan/Hu-1/2019 (NC_045512.2) and then visualized in UGENE v. 1.30. For phylogenetic analysis, a dataset of 88 available SARS-CoV-2 complete genomes from different countries, including Indonesia, was retrieved from GISAID.ResultsAll patients were hospitalized with various severities of COVID-19. Phylogenetic analysis revealed that one and three virus samples belong to clade L and GH. These three clade GH virus samples (EPI_ISL_525492, EPI_ISL_516800 and EPI_ISL_516829) were not only located in a cluster with SARS-CoV-2 genomes from Asia but also those from Europe, whereas the clade L virus sample (EPI_ISL_516806) was located amongst SARS-CoV-2 genomes from Asia. Using full-length sequences available in the GISAID EpiCoV Database, 39 of 60 SARS-CoV-2 (65%) from Indonesia harbor the D614G mutation.ConclusionThese findings indicate that SARS-CoV-2 with the D614G mutation appears to become the major circulating virus in Indonesia, concurrent with the COVID-19 situation worldwide.

Insights into the genetic diversity, recombination, and systemic infections with evidence of intracellular maturation of hepadnavirus in cats.

Hepatitis B virus (HBV) is a human pathogen of global concern, while a high diversity of viruses related to HBV have been discovered in other animals during the last decade. Recently, the novel mammalian hepadnavirus, tentatively named domestic cat hepadnavirus (DCH), was detected in an immunocompromised cat. Herein, a collection of 209 cat sera and 15 hepato-diseased cats were screened for DCH using PCR, resulting in 12.4% and 20% positivity in the tested sera and necropsied cats, respectively. Among the DCH-positive sera, a significantly high level of co-detection with retroviral infection was found, with the highest proportion being co-detection with feline immunodeficiency virus (FIV). Full-length genome characterization of DCH revealed the genetic diversity between the nine Thai DCH sequences obtained, and that they phylogenetically formed three distinct monophyletic clades. A putative DCH recombinant strain was found, suggesting a possible role of recombination in DCH evolution. Additionally, quantitative PCR was used to determine the viral copy number in various organs of the DCH-moribund cats, while the pathological findings were compared to the viral localization in hepatocytes, adjacent to areas of hepatic fibrosis, by immunohistochemical (IHC) and western blot analysis. In addition to the liver, positive-DCH immunoreactivity was found in various other organs, including kidneys, lung, heart, intestine, brain, and lymph nodes, providing evidence of systemic infection. Ultrastructure of infected cells revealed electron-dense particles in the nucleus and cytoplasm of hepatocytes, bronchial epithelial cells, and fibroblasts. We propose the intracellular development mechanism of this virus. Although the definitive roles of pathogenicity of DCH remains undetermined, a contributory role of the virus associated with systemic diseases is possible.

Human pegivirus (HPgV, GBV-C) RNA in volunteer blood donors from a public hemotherapy service in Northern Brazil

BackgroundHuman pegivirus (HPgV)—formerly known as GBV-C—is a member of the Flaviviridae family and belongs to the species Pegivirus C. It is a non-pathogenic virus and is transmitted among humans mainly through the exposure to contaminated blood and is often associated with human immunodeficiency virus (HIV) infection, among other viruses. This study aimed to determine the prevalence of HPgV viremia, its association with HIV and clinical epidemiological factors, as well as the full-length sequencing and genome characterization of HPgV recovered from blood donors of the HEMOPA Foundation in Belém-PA-Brazil.MethodsPlasma samples were obtained from 459 donors, tested for the presence of HPgV RNA by the RT-qPCR. From these, a total of 26 RT-qPCR positive samples were submitted to the NGS sequencing approach in order to obtain the full genome. Genome characterization and phylogenetic analysis were conducted.ResultsThe prevalence of HPgV was 12.42%. We observed the highest prevalences among donors aged between 18 and 30 years old (16.5%), with brown skin color (13.2%) and men (15.8%). The newly diagnosed HIV-1 prevalence was 26.67%. The HPgV genotype 2 (2a and 2b) was identified. No data on viral load value was found to corroborate the protective effect of HPgV on HIV evolution.ConclusionsThis study provided information regarding the HPgV infection among blood donors from HEMOPA Foundation. Furthermore, we genetically characterized the HPgV circulating strains and described by the first time nearly complete genomes of genotype 2 in Brazilian Amazon.

Near Full-Length Genomic Characterization of a Novel HIV-1 Circulating Recombinant Form (CRF107_01B) Identified Among Men Who Have Sex with Men in Heilongjiang Province of China.

We reported a novel human immunodeficiency virus type 1 (HIV-1) circulating recombinant form (CRF) among four epidemiologically unlinked patients through men having sex with men in Heilongjiang Province, China. It was named CRF107_01B (this is temporary as we have not received the CRF number from HIV databases). A near full-length genome phylogenetic tree showed that CRF107_01B was generated by two CRF01_AE segments that was described as cluster 5 lineage of CRF01_AE (CRF01-5) inserted into the B (Western origin) backbone. The emergence of CRF107_01B increased the complexity of HIV-1.

Near Full-Length Genomic Characterization of a Novel HIV-1 Circulating Recombinant Form (CRF106_cpx) Identified Among Heterosexuals in China.

In our study, we reported a novel CRF01_AE/B/C human immunodeficiency virus type 1 (HIV-1) recombinant form among six epidemiologically unlinked heterosexual patients in Yunnan, China. It was named CRF106_cpx (this is temporary as we have not received the circulating recombinant form [CRF] number from HIV databases). After sequencing and analyzing of the near full-length genome, we found that CRF106_cpx was generated by three B segments and two CRF01_AE segments inserted into the C backbone. There was quite a bit of diversity between sequences of CRF106_cpx, so this seems to be a CRF that has been spreading for many years.

Full-length genome characterization of a novel recombinant vaccine-like lumpy skin disease virus strain detected during the climatic winter in Russia, 2019.

An uncharacteristic outbreak of lumpy skin disease was reported in the Republic of Udmurtiya, Russia, during the climatic winter of March 2019. The causative lumpy skin disease virus (LSDV_Udmurtiya_Russia_2019) was shown to be a recombinant composed of a live attenuated Neethling-type vaccine strain as the dominant parental strain and a Kenyan KSGP/NI-2490-like virus as its minor parental strain, with 24 statistically significant recombination events that are not identical to those in LSDV Saratov/2017, in which 27 events were identified.

The full-length genome characterization and diversity of faba bean necrotic stunt virus in Iran

Faba bean necrotic stunt virus (FBNSV) (genus Nanovirus; family Nanoviridae) is a multipartite virus comprising eight circular single-stranded DNA components, each about 1 kb in length. Also, two alphasatellites have often been reported to be associated with FBNSV infection. The only full-length genome sequence of an Iranian FBNSV isolate, FBNSV-[IR;Ka:FB], so far reported is one infecting faba bean (Vicia faba L.). To obtain more information about Iranian isolates of FBNSV (FBNSV-IRs) and their diversity, the full-length genome of four FBNSV isolates and associated alphasatellites obtained from faba bean (two isolates), common bean (Phaseolus vulgaris L., one isolate) and chickpea (Cicer arietinum L., one isolate) were characterized. All genomic components except U4 had the expected size in accordance with the FBNSV genomic sequences published in the GenBank. The component U4 had an additional stretch of 16 nucleotides (nts 765–780) in its non-coding region (NCR). Alphasatellites of the aforementioned isolates showed 98–100% sequence identity. Intra-component sequence analyses of FBNSV-IRs indicated a range of 98 to 100% nt identity between each single FBNSV-IRs component or between their associated alphasatellites. DNA-M and DNA-U2 exhibited the highest value of diversity, sharing 98–99% identity. Phylogenetic analyses indicated that FBNSV-IRs were clustered with a Moroccan isolate and some Azerbaijani isolates of FBNSV forming a clade distinct from Ethiopian isolates. Inter-component sequence analyses of NCRs revealed the DNA-S and DNA-U4 subgrouped together in all FBNSV-IRs isolates except the chickpea isolate in which NCRs of DNA-S and DNA–U4 formed separate subgroups.

Near Full-Length Genomic Characterization of a Novel Unique Recombinant (CRF01_AE/CRF07_BC) in Fuyang City of China.

Recombinant forms contribute substantially to the genetic diversity of HIV-1. Recent studies have also revealed that three major viral strains (CRF07_BC, CRF01_AE, and subtype B) have been cocirculating among men who have sex with men (MSM) in Anhui, suggesting a high probability of generating new recombinants. In this study, we reported a novel CRF01_AE and CRF07_BC HIV-1 recombinant form in MSM in Fuyang city of China. Two near full-length genome (NFLG) named FY184 and FY208 were successfully obtained. The genomic composition analysis of the NFLG reveals that it was divided into four segments by three breakpoints, with two regions of CRF07_BC inserted into a CRF01_AE backbone's gag and pol regions. The CRF01_AE regions were originated from a subcluster lineage of CRF01_AE, which is mainly circulating among MSM in China. The emergence of a novel recombinant of CRF01_AE/CRF07_BC is indicative of the increasing genetic diversity of the HIV epidemic in MSM in Anhui.

Characterization of a New Isolate of Chilli ringspot virus in Yunnan, China

Virus diseases are the major limiting factors of chilli production in China. A virus isolate of chilli ringspot virus-Yunnan (ChiRSV-YN) from chilli (Capsicum annuum L.) showing mild mosaic, leaf crinkling, ringspot, and vein banding symptoms was recently identified in Yunnan, China. The electron microscopy confirmed the virus as a typical potyvirus having flexuous filamentous particles ca. 780 nm in length. Several chilli and tobacco varieties were tested through mechanical inoculation and found to be susceptible to the virus isolate ChiRSV-YN. The viral genome was sequenced through RNA sequencing. The complete genomic RNA of the isolate (GenBank Acc. No. KX258620) consisted of 9,652 nucleotides (nt) excluding the poly (A) tail at the 3′ end. The isolate contains a large ORF which encodes a polyprotein of 3,086 amino acids (aa) with Mr. of 349.21 kDa. It had a typical genomic organization of potyviruses. The complete nucleotide and polyprotein sequences of ChiRSV-YN had 87.6 and 93.2% identities with the ChiRSV-HN/14, respectively. The result showed ChiRSV-YN is a new isolate of ChiRSV and it is the first full-length genome characterization of ChiRSV in Yunnan, China. Understanding of the characteristics of this virus isolate might provide the basis for disease diagnosis and control in the future.

Nearly full-length genome characterization of canine parvovirus strains circulating in Nigeria.

Canine parvovirus type 2 (CPV‐2) emerged suddenly in the late 1970s as pathogen of dogs, causing a severe and often fatal gastroenteric disease. The original CPV‐2 was replaced by three antigenic variants, CPV‐2a, CPV‐2b and CPV‐2c, which to date have gained a worldwide distribution with different relative proportions. All previous studies conducted in Africa were based on partial VP2 gene sequences. The aim of this study was to provide a genome analysis to characterize the CPV strains collected in Nigeria, Africa. Rectal swab samples (n = 320) were collected in 2018 and tested by means of an immunochromatographic assay. Among the 144 positive samples, 59 were selected for further analyses using different molecular assays. The results revealed a high prevalence of CPV‐2c (91.5%) compared to the CPV‐2a variant (8.5%). The VP2 gene sequences showed a divergence from the strains analysed in 2010 in Nigeria and a closer connection with CPV strains of Asian origin. The non‐structural gene analysis evidenced amino acid changes never previously reported. The molecular analysis based on genomic sequences evidenced a geographical pattern of distribution of the analysed strains, suggesting a potential common evolutionary origin with CPV of Asian origin. This study represents the first CPV molecular characterization including all the encoding gene sequences conducted in the African continent and contributes to define the current geographical spread of the CPV variants worldwide.

Retrospective study, full-length genome characterization and evaluation of viral infectivity and pathogenicity of chimeric porcine deltacoronavirus detected in Vietnam.

Increased evidence of porcine deltacoronavirus (PDCoV) causing diarrhoea in pigs has been reported in several countries worldwide. The virus has currently evolved into three separated groups including US, China and Southeast Asia (SEA) groups. In Vietnam, PDCoV was first reported in 2015. Based on phylogenetic analyses of spike, membrane and nucleocapsid genes, it is suggested that Vietnam PDCoV is chimeric virus. In the present study, we retrospectively investigated the presence of PDCoV in Vietnam and the full‐length genomes of six PDCoV isolates identified in 2014–2016 were further characterized. The results demonstrated that Vietnam PDCoV was first detected as early as 2014. All six Vietnam PDCoV are in the SEA group and further divided into two separated subgroups including SEA‐1 and SEA‐2. Vietnam PDCoV in SEA‐2 was closely related to Thai and Lao PDCoV. Recombination analysis demonstrated that three isolates in SEA‐1 were a chimeric virus of which P12_14_VN_0814, the first Vietnam isolate, and US PDCoV isolates were major and minor parents, respectively. The recombination was further evaluated by phylogenetic construction based on 3 recombinant fragments. The first and third fragments, closely related to P12_14_VN_0814, were associated with ORF1a/1b and N genes, respectively. The second fragment, associated with S, E, and M genes, was closely related to US PDCoV isolates. High antigenic and hydrophobic variations were detected in S1 protein. Three‐day‐old pigs challenged with the chimeric virus displayed clinical diseases and villus atrophy. In conclusion, Vietnam PDCoV is genetically diverse influenced by an external introduction from neighbouring countries. The chimeric Vietnam PDCoV can induce a disease similar to Thai PDCoV.

Full-Length Genomic Characterization of Chicken Anemia Virus in Turkey

ABSTRACT Chicken infectious anemia (CIA) is an immune-suppressive disease caused by chicken anemia virus (CAV). It is characterized by lymphoid atrophy, aplastic anemia, especially in chicks. In this study, full-length genomic characterization of CAV DNA from the broiler flocks in Turkey and phylogenetic analysis were aimed. In the study, CAV DNA were found positive for 37 (53%) flocks with PCR studies from thymus tissues of each 70 broiler flocks. And 17 purified CAV DNA PCR products from these 37 CAV isolates were full length sequenced with the NGS method (Illumina MiSeq). Also with the phylogenetic analyses, full length PCR products of 17 purified CAV isolates have been determined as 2298bp genome size and 99% similarity with each other. The highest similarity (99%) has been detected with the isolates from China and Taiwan. Furthermore, a 97-98% similarity has been detected with vaccine strains (Cux-1, 26P4 and Del Ros) and also 88-90 % similarity has been detected with GyV4 and GyV3 isolates. As a result, in the study full length genomic characterization of CAV DNA from the 7 regions of Turkey were determined. And also all Turkish CAV isolates and vaccine strains were in group 2 according to the phylogenetic tree were obtained. But these isolates and vaccine strains were not found in the same group with GyV3 and GyV4 strains. Besides, these CAV isolates were showed more similarity to the isolates reported from Taiwan and China than the vaccine strains.

Oropouche orthobunyavirus: Genetic characterization of full-length genomes and development of molecular methods to discriminate natural reassortments

Oropouche orthobunyavirus (OROV) has significant impact in public health in Amazon region. This arbovirus is one of the most common causes of febrile illness in Brazil, and is responsible for several epidemics since 1960's. In this study, we sequenced and characterized the complete coding sequences (S-, M- and L-RNA) of 35 OROV isolates from Brazil. Here, we classified 20 strains in genotype I from Pará and Maranhão states, nine as genotype II from Pará and Rondônia states confirmed, four classified into genotype III from Acre, Maranhão, Minas Gerais and Rondônia states and two genotype IV from Amazonas State. Also, we did not observe reassortment events involving the OROV isolates. In addition, we developed novel RT-PCR tools to identify reassortment events among OROV strains. These data will be useful to better understand the molecular epidemiology and diagnostic of OROV infections.

Near Full-Length Genomic Characterization of HIV-1 CRF01_AE/B Recombinant Strains Identified in Hebei, China.

This study reports a novel HIV-1 recombinant form of CRF01_AE/B that was isolated from a 22-year-old male individual (HE150308) residing in Hebei province, China. The near full-length genome (NFLG) phylogenetic tree showed that the strain was clustered with CRF01_AE reference sequences and placed at the peripheral position within the branch of CRF01_AE strains. The bootscaning and similarity plot analysis revealed that the NFLG of this novel recombinant was composed of eight interlaced segments, including four CRF01_AE and four subtype B segments separated by seven breakpoints observed in the gag, pol, and nef regions of HIV-1 genome. Phylogenetic analysis demonstrated that the CRF01_AE subregions were from the previously identified CRF01_AE cluster 5, and the B subregions were correlated with the B strains originated from Europe and America. They were all the lineages widely prevalent in men who have sex with men (MSM) population in China. In recent years, a large number of recombinants between CRF01_AE and B strains are constantly emerging in the MSM population in China. This continual and recurrent recombination between CRF01_AE and B in high-risk group people deserves more attention and further monitoring.

Molecular Characterization of Near Full-Length Genomes of Hepatitis B Virus Isolated from Predominantly HIV Infected Individuals in Botswana.

The World Health Organization plans to eliminate hepatitis B and C Infections by 2030. Therefore, there is a need to study and understand hepatitis B virus (HBV) epidemiology and viral evolution further, including evaluating occult (HBsAg-negative) HBV infection (OBI), given that such infections are frequently undiagnosed and rarely treated. We aimed to molecularly characterize HBV genomes from 108 individuals co-infected with human immunodeficiency virus (HIV) and chronic hepatitis B (CHB) or OBI identified from previous HIV studies conducted in Botswana from 2009 to 2012. Full-length (3.2 kb) and nearly full-length (~3 kb) genomes were amplified by nested polymerase chain reaction (PCR). Sequences from OBI participants were compared to sequences from CHB participants and GenBank references to identify OBI-unique mutations. HBV genomes from 50 (25 CHB and 25 OBI) individuals were successfully genotyped. Among OBI participants, subgenotype A1 was identified in 12 (48%), D3 in 12 (48%), and E in 1 (4%). A similar genotype distribution was observed in CHB participants. Whole HBV genome sequences from Botswana, representing OBI and CHB, were compared for the first time. There were 43 OBI-unique mutations, of which 26 were novel. Future studies using larger sample sizes and functional analysis of OBI-unique mutations are warranted.

Full-length genome and molecular characterization of dengue virus serotype 2 isolated from an imported patient from Myanmar

BackgroundDengue is the most common mosquito-borne infection worldwide and a serious threat to global public health. Sporadic dengue virus serotype 2 (DENV-2) imported cases from Myanmar have been documented almost every year in Yunnan Province of China since 2005. However, the complete genome sequences of DENV-2 isolates imported from Myanmar are not available.MethodsThe full-length genome of the DENV-2 strain (YNPE2), isolated from an imported case from Myanmar in 2013, was identified by the next-generation sequencing. The extreme ends of the viral genome were validated by 5′/3′ RACE and Sanger sequencing. Furthermore, phylogenetic, recombination and selection pressure analyses were conducted for the molecular characterization of YNPE2 strain.ResultsWhole-genome sequencing revealed that the full-length sequence of YNPE2 strain was 10,724 bases, with an open reading frame encoding for 3391 amino acids. The YNPE2 strain had 99.0% nucleotide identity and 99.8% amino acid identity with two closely related strains, ThD2_0078_01 strain (DQ181797) and DENV-2/TH/BID-V2157/200 strain (FJ639832). The phylogenetic analysis suggested that the YNPE2 strain belonged to Asian I genotype and was likely derived from Thailand strain (DQ181797). Moreover, selection pressure analysis revealed two amino acid sites of the NS4B and NS5 proteins, with important evidence of positive selection.ConclusionThis study revealed the first complete genome sequence and molecular characterization of a DENV-2 strain (YNPE2) isolated from an imported case from Myanmar, thus providing a valuable reference genome source for future surveillance, epidemiology and vaccine development of DENV-2 virus in Yunnan, China.

The full-length genome characterization, genetic diversity and evolutionary analyses of Senecavirus A isolated in Thailand in 2016

Senecavirus A (SVA) is a novel picornavirus that causes porcine idiopathic vesicular disease characterized by lameness, coronary band hyperemia, and vesicles on the snout and coronary bands. An increase in the detection rate of SVA in several countries suggests that the disease has become a widespread problem. Herein, we report the detection of SVA in Thailand and the characterization of full-length genomic sequences of six Thai SVA isolates. Phylogenetic, genetic, recombination, and evolutionary analyses were performed. The full-length genome, excluding the poly (A) tail of the Thai SVA isolates, was 7282 nucleotides long, with the genomic organization resembling other previously reported SVA isolates. Phylogenetic and genetic analyses based on full-length genome demonstrated that the Thai SVA isolates were grouped in a novel cluster, separated from SVA isolates from other countries. Although the Thai SVA isolates were closely related to 11-55910-3, the first SVA isolate from Canada, with 97.9–98.2%, but they are different. Evolutionary and recombinant analyses suggested that the Thai SVA isolates shared a common ancestor with the 11-55910-3 isolate. The positive selection in the VP4 and 3D genes suggests that the virus was not externally introduced, but rather continuously evolved in the population prior to the first detection. Addition, the presence of SVA could have been ignored due to the presence of other pathogens causing similar clinical diseases. This study warrants further investigations into molecular epidemiology and genetic evolution of the SVA in Thailand.

Short Communication: Characterization of a New HIV-1 Group N Isolate Originating from a Cameroonian Patient

HIV-1 group N (HIV-1/N) remains rare and mainly restricted to Cameroon. In this study, we report a new HIV-1/N infected case identified during routine HIV screening activities in Yaounde. The genetic characterization of the near full-length genome of this virus strain revealed that it is genetically distinct to all HIV-1/N described to date. However, the Vpu protein responsible for tetherin antagonism displayed the same amino acid substitutions (E15A, V19A, I25L, and V26L) as other HIV-1/N from Cameroon.

Archetype JC Polyomavirus (JCPyV) Prevails in a Rare Case of JCPyV Nephropathy and in Stable Renal Transplant Recipients With JCPyV Viruria

JC polyomavirus (JCPyV) is reactivated in approximately 20% of renal transplant recipients, and it may rarely cause JCPyV-associated nephropathy (JCPyVAN). Whereas progressive multifocal leukoencephalopathy of the brain is caused by rearranged neurotropic JCPyV, little is known about viral sequence variation in JCPyVAN owing to the rarity of this condition. Using single-molecule real-time sequencing, characterization of full-length JCPyV genomes in urine and plasma samples from 1 patient with JCPyVAN and 20 stable renal transplant recipients with JCPyV viruria was attempted. Sequence analysis of JCPyV strains was performed, with emphasis on the noncoding control region, the major capsid protein gene VP1, and the large T antigen gene. Exclusively archetype strains were identified in urine from the patient with JCPyVAN. Full-length JCPyV sequences were not retrieved from plasma. Archetype strains were found in urine samples from 19 stable renal transplant recipients, with JCPyV quasispecies detected in 5 samples. In a patient with minor graft dysfunction, a strain with an archetype-like noncoding cont rol region was discovered. Individual point mutations were detected in both VP1 and large T antigen genes. Archetype JCPyV was dominant in the patient with JCPyVAN and in stable renal transplant recipients. Archetype rather than rearranged JCPyV seems to drive the pathogenesis of JCPyVAN.

A robust and cost-effective approach to sequence and analyze complete genomes of small RNA viruses

BackgroundNext-generation sequencing (NGS) allows ultra-deep sequencing of nucleic acids. The use of sequence-independent amplification of viral nucleic acids without utilization of target-specific primers provides advantages over traditional sequencing methods and allows detection of unsuspected variants and co-infecting agents. However, NGS is not widely used for small RNA viruses because of incorrectly perceived cost estimates and inefficient utilization of freely available bioinformatics tools.MethodsIn this study, we have utilized NGS-based random sequencing of total RNA combined with barcode multiplexing of libraries to quickly, effectively and simultaneously characterize the genomic sequences of multiple avian paramyxoviruses. Thirty libraries were prepared from diagnostic samples amplified in allantoic fluids and their total RNAs were sequenced in a single flow cell on an Illumina MiSeq instrument. After digital normalization, data were assembled using the MIRA assembler within a customized workflow on the Galaxy platform.ResultsTwenty-eight avian paramyxovirus 1 (APMV-1), one APMV-13, four avian influenza and two infectious bronchitis virus complete or nearly complete genome sequences were obtained from the single run. The 29 avian paramyxovirus genomes displayed 99.6% mean coverage based on bases with Phred quality scores of 30 or more. The lower and upper quartiles of sample median depth per position for those 29 samples were 2984 and 6894, respectively, indicating coverage across samples sufficient for deep variant analysis. Sample processing and library preparation took approximately 25–30 h, the sequencing run took 39 h, and processing through the Galaxy workflow took approximately 2–3 h. The cost of all steps, excluding labor, was estimated to be 106 USD per sample.ConclusionsThis work describes an efficient multiplexing NGS approach, a detailed analysis workflow, and customized tools for the characterization of the genomes of RNA viruses. The combination of multiplexing NGS technology with the Galaxy workflow platform resulted in a fast, user-friendly, and cost-efficient protocol for the simultaneous characterization of multiple full-length viral genomes. Twenty-nine full-length or near-full-length APMV genomes with a high median depth were successfully sequenced out of 30 samples. The applied de novo assembly approach also allowed identification of mixed viral populations in some of the samples.