key: cord-0834136-u66yn560 authors: Hashish, Amro; Sato, Yuko; Li, Ganwu; Zheng, Ying; Gauger, Phillip C.; El-Gazzar, Mohamed title: Near-Complete Genome Sequence of GI-17 Lineage Infectious Bronchitis Virus, Circulating in Iowa date: 2021-05-20 journal: Microbiol Resour Announc DOI: 10.1128/mra.01406-20 sha: 29a28d6ecf56356aee8b6a891de7d9ca665b8be4 doc_id: 834136 cord_uid: u66yn560 Avian infectious bronchitis virus (AvIBV) is the causative agent of a highly contagious respiratory disease in chickens which results in significant economic losses in the poultry industry. Here, we report a near-complete genome sequence of the strain, designated IA1162/2020, identified in tracheal swabs from chickens in Iowa in 2020. A vian infectious bronchitis virus (AvIBV) belongs to the genus Gammacoronavirus, family Coronaviridae, subfamily Coronavirinae, and order Nidovirales (1). AvIBV is characterized by high genetic diversity and the continuous emergence of new variant strains (2) . In 2011, a nephropathogenic strain of AvIBV, Delmarva/1639 (DMV/1639), from genotype I lineage 17 (GI-17) was reported (3) . Despite the recent wide distribution of this strain in the United States and the availability of multiple S1 gene sequences, whole-genome sequences of this strain are limited. In this study, we report the near-complete genome sequence of a new GI-17 strain, designated IA1162/2020, related to the DMV/1639 strain. Two pooled tracheal swabs were collected from 10-day-old layer chicks in Iowa with mild respiratory signs and submitted to the Iowa State University Veterinary Diagnostic Laboratory (ISU-VDL) in January 2020. The samples tested positive by reverse transcription-quantitative PCR (qRT-PCR) (4) with a threshold cycle (C T ) value of 18. Total nucleic acid was reextracted from the tracheal swabs using the MagMAX pathogen RNA/DNA kit (Thermo Fisher Scientific, Massachusetts) with a KingFisher 96 instrument (Thermo Fisher Scientific) (5). Double-stranded cDNA was synthesized using the NextFlex rapid transcriptome sequencing (RNA-seq) kit (Bioo Scientific Corp., Texas). The sequencing library was prepared using the Nextera XT DNA library preparation kit (Illumina, California) with dual indexing. The pooled libraries were sequenced on an Illumina MiSeq platform using the 300-cycle v2 reagent kit (Illumina). The raw reads were analyzed and preprocessed using Trimmomatic v0.36 (6) . The cleaned reads were classified using Kraken v0.10.5-beta (7) with the standard database. Unclassified reads were classified using Kaiju v1.6.2 (8), and KronaTools v2.6 (9) was used to generate the interactive html charts for the hierarchical classification results. Out of a total of 603,724 reads generated for this sample, 8,214 raw reads were identified as IBV. AvIBV reads were extracted from the classification results for de novo assembly using ABySS v1.3.9 (10). The resulting contigs were manually refined and curated to remove contaminated (nonviral) contigs and to trim chimeric (misassembled) contigs in the SeqMan Pro DNASTAR Lasergene 11 core suite. Finally, a near-complete sequence of AvIBV IA1162/2020 (99.9% genome coverage based on the GA9977/2019 reference genome FIG 1 Phylogenetic tree generated from the complete S1 gene sequences of IA1162/2020 and other avian infectious bronchitis virus GI-17 strains. The analysis was done with the neighbor-joining method, Kimura 2 parameters (15) , and MEGA X (16) . Numbers at the nodes correspond to bootstrap values (1,000 replications). Bootstrap values below 70 were removed. The IA1162/2020 strain is marked with a red triangle. All ambiguous positions were removed for each sequence pair (pairwise deletion option). There were a total of 1,773 positions in the final data set. This analysis involved 35 nucleotide sequences. Sequences from the GII-1 lineage are included as an outgroup. [GenBank accession number MK878536.1]) was identified with a genome length of 27,486 nucleotides, excluding the poly(A) tail, with a GC content of 38%. Twelve open reading frames (ORFs) were predicted (Table 1) using the ORFfinder program (https://www.ncbi.nlm.nih.gov/orffinder/), resulting in the following genome organization: 59-1a-1ab-S-3a-3b-E-M-4b-4c-5a-5b-N-39. This genome (IA1162/2020) is missing ORF 6b, which was supposed to be between the nucleocapsid gene (N) and the 39 untranscribed region (UTR). ORF 6b (74 amino acids [aa]) is thought to be an apoptosis inducer (11) and accelerated the replication of murine coronavirus (CoV) in vitro (12) . Full-length genome sequence BLAST comparison of IA1162/2020 showed the highest nucleotide similarity (94.78%) to Georgia strain GA9977/2019 (GenBank accession number MK878536.1) belonging to lineage GI-17 (13) . The sequence of the full S1 of IA1162/2020 clustered within GI-17 (14) ; however, it formed a separate clade away from other GI-17 strains (Fig. 1) . This near-complete genome will be useful for developing new diagnostic assays and new and alternative vaccination strategies and for better understanding the evolution of AvIBV. Data availability. The near-complete genome sequence of the IA1162/2020 AvIBV DMV/1639 strain has been deposited in GenBank under the accession number MW024789. The raw data were deposited under SRA number SRR12575717, BioSample number SAMN15963238, and BioProject number PRJNA660897. Discovery of seven novel mammalian and avian coronaviruses in the genus deltacoronavirus supports bat coronaviruses as the gene source of alphacoronavirus and betacoronavirus and avian coronaviruses as the gene source of gammacoronavirus and deltacoronavirus Genetic analysis of avian coronavirus infectious bronchitis virus in yellow chickens in southern china over the past decade: revealing the changes of genetic diversity, dominant genotypes, and selection pressure Characterization of nephropathogenic infectious bronchitis virus DMV/1639/11 recovered from Delmarva broiler chickens in 2011 Development and evaluation of a real-time Taqman RT-PCR assay for the detection of infectious bronchitis virus from infected chickens High-throughput whole genome sequencing of porcine reproductive and respiratory syndrome virus from cell culture materials and clinical specimens using next-generation sequencing technology Trimmomatic: a flexible trimmer for Illumina sequence data Kraken: ultrafast metagenomic sequence classification using exact alignments Fast and sensitive taxonomic classification for metagenomics with Kaiju Interactive metagenomic visualization in a Web browser ABySS: a parallel assembler for short read sequence data Genomic and single nucleotide polymorphism analysis of infectious bronchitis coronavirus Severe acute respiratory syndrome coronavirus protein 6 accelerates murine coronavirus infections First complete genome sequence of currently circulating infectious bronchitis virus strain DMV/1639 of the GI-17 lineage S1 gene-based phylogeny of infectious bronchitis virus: an attempt to harmonize virus classification A simple method for estimating evolutionary rates of base substitutions through comparative studies of nucleotide sequences MEGA X: molecular evolutionary genetics analysis across computing platforms We thank the ISU-VDL faculty and staff for assistance with sample processing and diagnostic testing.This study was supported by the ISU-VDL.