Viral proteogenomic and expression profiling during productive replication of a skin-tropic herpesvirus in the natural host.

Efficient transmission of herpesviruses is essential for dissemination in host populations; however, little is known about the viral genes that mediate transmission, mostly due to a lack of natural virus-host model systems. Marek's disease is a devastating herpesviral disease of chickens caused by Marek's disease virus (MDV) and an excellent natural model to study skin-tropic herpesviruses and transmission. Like varicella zoster virus that causes chicken pox in humans, the only site where infectious cell-free MD virions are efficiently produced is in epithelial skin cells, a requirement for host-to-host transmission. Here, we enriched for heavily infected feather follicle epithelial skin cells of live chickens to measure both viral transcription and protein expression using combined short- and long-read RNA sequencing and LC/MS-MS bottom-up proteomics. Enrichment produced a previously unseen breadth and depth of viral peptide sequencing. We confirmed protein translation for 84 viral genes at high confidence (1% FDR) and correlated relative protein abundance with RNA expression levels. Using a proteogenomic approach, we confirmed translation of most well-characterized spliced viral transcripts and identified a novel, abundant isoform of the 14 kDa transcript family via IsoSeq transcripts, short-read intron-spanning sequencing reads, and a high-quality junction-spanning peptide identification. We identified peptides representing alternative start codon usage in several genes and putative novel microORFs at the 5' ends of two core herpesviral genes, pUL47 and ICP4, along with strong evidence of independent transcription and translation of the capsid scaffold protein pUL26.5. Using a natural animal host model system to examine viral gene expression provides a robust, efficient, and meaningful way of validating results gathered from cell culture systems.

Proteome-wide Analysis of Protein Thermal Stability in the Model Higher Plant Arabidopsis thaliana.

Modern tandem MS-based sequencing technologies allow for the parallel measurement of concentration and covalent modifications for proteins within a complex sample. Recently, this capability has been extended to probe a proteome's three-dimensional structure and conformational state by determining the thermal denaturation profile of thousands of proteins simultaneously. Although many animals and their resident microbes exist under a relatively narrow, regulated physiological temperature range, plants take on the often widely ranging temperature of their surroundings, possibly influencing the evolution of protein thermal stability. In this report we present the first in-depth look at the thermal proteome of a plant species, the model organism Arabidopsis thaliana By profiling the melting curves of over 1700 Arabidopsis proteins using six biological replicates, we have observed significant correlation between protein thermostability and several known protein characteristics, including molecular weight and the composition ratio of charged to polar amino acids. We also report on a divergence of the thermostability of the core and regulatory domains of the plant 26S proteasome that may reflect a unique property of the way protein turnover is regulated during temperature stress. Lastly, the highly replicated database of Arabidopsis melting temperatures reported herein provides baseline data on the variability of protein behavior in the assay. Unfolding behavior and experiment-to-experiment variability were observed to be protein-specific traits, and thus this data can serve to inform the design and interpretation of future targeted assays to probe the conformational status of proteins from plants exposed to different chemical, environmental and genetic challenges.

Runting and Stunting Syndrome in Broiler Chickens: Histopathology and Association With a Novel Picornavirus.

Runting stunting syndrome (RSS) in commercial chickens has been reported worldwide, and although several studies have attempted to clarify the cause and describe the lesions, there are gaps in knowledge of the epidemiology, pathogenesis, and etiology. The study objective was to use commercial chicks naturally affected by RSS to describe the histologic changes of RSS in all segments of the small intestine in chicks of different ages and to identify viral gene sequences in affected chicks and their association with histologic lesions. Chicks lacking clinical signs but from the same houses and from unaffected houses were used as controls. The average weight of affected chicks was significantly lower than expected for their flocks. Macroscopically, the small intestines had paler serosa, with watery, mucoid, or foamy contents and poorly digested food. Histologic lesions were characterized by necrotic crypts, crypt dilation, and flattening of the crypt epithelium. Histomorphometry of the intestines revealed villous atrophy especially in the jejunum and ileum. Histologic changes in other organs were not observed. Random next-generation sequencing of total RNA extracted from formalin-fixed paraffin-embedded tissues detected avian nephritis virus, avian rotavirus, and picornavirus in jejunal segments from 7-day-old chicks. No viruses were detected in the jejunum of 1-day-old chicks. Detection of picornaviral reads was significantly associated (P < .05) with histologic lesions of RSS. Sequence analysis of the picornavirus revealed genetic similarity with the genus Gallivirus. Using in situ hybridization for galliviral nucleic acid sequences, the signal was associated with crypt lesion severity, although signal was detected both in chicks with and without RSS.

MinION sequencing to genotype US strains of infectious laryngotracheitis virus.

Over the last decade the US broiler industry has fought long-lasting outbreaks of infectious laryngotracheitis (ILTV). Previously, nine genotypes (I-IX) of ILTVs have been recognized using the polymerase chain reaction-restriction fragment length polymorphisms (PCR-RFLP) method with three viral alleles (gB, gM and UL47/gG). In this study, the genotyping system was simplified to six genotypes by amplicon sequencing and examining discriminating single nucleotide polymorphisms (SNPs) within these open reading frames. Using phylogenomic analysis of 27 full genomes of ILTV, a single allele (ORF A/ORF B) was identified containing SNPs that could differentiate ILTVs into genotypes congruent with the phylogenetic partitioning. The allelic variations allowed for the cataloging of the 27 strains into 5 genotypes: vaccinal TCO, vaccinal CEO, virulent CEO-like, virulent US and virulent US backyard flocks from 1980 to 1990, correlating with the PCR-RFLP genotypes I/ II/ III (TCO), IV (CEO), V (virulent CEO-like), VI (virulent US) and VII/VIII/IX (virulent US backyard flock isolates). With the unique capabilities of third generation sequencing, we investigated the application of Oxford Nanopore MinION technology for rapid sequencing of the amplicons generated in the single-allele assay. This technology was an improvement over Sanger-based sequencing of the single allele amplicons due to a booster amplification step in the MinION sequencing protocol. Overall, there was a 90% correlation between the genotyping results of the single-allele assay and the multi-allele assay. Surveillance of emerging ILTV strains could greatly benefit from real-time amplicon sequencing using the single-allele assay and MinION sequencing. RESEARCH HIGHLIGHTS A multi-allelic assay identified nine ILTV genotypes circulating in the US Single-allele genotyping is congruent with whole genome phylogenetic partitioning US ILTV strains can be grouped into five genotypes using the single-allele assay The single-allele assay can be done using MinION sequencing of barcoded amplicons.

Rapid virulence prediction and identification of Newcastle disease virus genotypes using third-generation sequencing.

BACKGROUND: Newcastle disease (ND) outbreaks are global challenges to the poultry industry. Effective management requires rapid identification and virulence prediction of the circulating Newcastle disease viruses (NDV), the causative agent of ND. However, these diagnostics are hindered by the genetic diversity and rapid evolution of NDVs. METHODS: An amplicon sequencing (AmpSeq) workflow for virulence and genotype prediction of NDV samples using a third-generation, real-time DNA sequencing platform is described here. 1D MinION sequencing of barcoded NDV amplicons was performed using 33 egg-grown isolates, (15 NDV genotypes), and 15 clinical swab samples collected from field outbreaks. Assembly-based data analysis was performed in a customized, Galaxy-based AmpSeq workflow. MinION-based results were compared to previously published sequences and to sequences obtained using a previously published Illumina MiSeq workflow. RESULTS: For all egg-grown isolates, NDV was detected and virulence and genotype were accurately predicted. For clinical samples, NDV was detected in ten of eleven NDV samples. Six of the clinical samples contained two mixed genotypes as determined by MiSeq, of which the MinION method detected both genotypes in four samples. Additionally, testing a dilution series of one NDV isolate resulted in NDV detection in a dilution as low as 101 50% egg infectious dose per milliliter. This was accomplished in as little as 7 min of sequencing time, with a 98.37% sequence identity compared to the expected consensus obtained by MiSeq. CONCLUSION: The depth of sequencing, fast sequencing capabilities, accuracy of the consensus sequences, and the low cost of multiplexing allowed for effective virulence prediction and genotype identification of NDVs currently circulating worldwide. The sensitivity of this protocol was preliminary tested using only one genotype. After more extensive evaluation of the sensitivity and specificity, this protocol will likely be applicable to the detection and characterization of NDV.

Algal ancestor of land plants was preadapted for symbiosis.

Colonization of land by plants was a major transition on Earth, but the developmental and genetic innovations required for this transition remain unknown. Physiological studies and the fossil record strongly suggest that the ability of the first land plants to form symbiotic associations with beneficial fungi was one of these critical innovations. In angiosperms, genes required for the perception and transduction of diffusible fungal signals for root colonization and for nutrient exchange have been characterized. However, the origin of these genes and their potential correlation with land colonization remain elusive. A comprehensive phylogenetic analysis of 259 transcriptomes and 10 green algal and basal land plant genomes, coupled with the characterization of the evolutionary path leading to the appearance of a key regulator, a calcium- and calmodulin-dependent protein kinase, showed that the symbiotic signaling pathway predated the first land plants. In contrast, downstream genes required for root colonization and their specific expression pattern probably appeared subsequent to the colonization of land. We conclude that the most recent common ancestor of extant land plants and green algae was preadapted for symbiotic associations. Subsequent improvement of this precursor stage in early land plants through rounds of gene duplication led to the acquisition of additional pathways and the ability to form a fully functional arbuscular mycorrhizal symbiosis.

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

BACKGROUND: Next-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. METHODS: In 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. RESULTS: Twenty-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. CONCLUSIONS: This 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.

Unique patterns of transcript and miRNA expression in the South American strong voltage electric eel (Electrophorus electricus).

BACKGROUND: With its unique ability to produce high-voltage electric discharges in excess of 600 volts, the South American strong voltage electric eel (Electrophorus electricus) has played an important role in the history of science. Remarkably little is understood about the molecular nature of its electric organs. RESULTS: We present an in-depth analysis of the genome of E. electricus, including the transcriptomes of eight mature tissues: brain, spinal cord, kidney, heart, skeletal muscle, Sachs' electric organ, main electric organ, and Hunter's electric organ. A gene set enrichment analysis based on gene ontology reveals enriched functions in all three electric organs related to transmembrane transport, androgen binding, and signaling. This study also represents the first analysis of miRNA in electric fish. It identified a number of miRNAs displaying electric organ-specific expression patterns, including one novel miRNA highly over-expressed in all three electric organs of E. electricus. All three electric organ tissues also express three conserved miRNAs that have been reported to inhibit muscle development in mammals, suggesting that miRNA-dependent regulation of gene expression might play an important role in specifying an electric organ identity from its muscle precursor. These miRNA data were supported using another complete miRNA profile from muscle and electric organ tissues of a second gymnotiform species. CONCLUSIONS: Our work on the E. electricus genome and eight tissue-specific gene expression profiles will greatly facilitate future research on determining the coding and regulatory sequences that specify the function, development, and evolution of electric organs. Moreover, these data and future studies will be informed by the first comprehensive analysis of miRNA expression in an electric fish presented here.

Rapid phosphoproteomic and transcriptomic changes in the rhizobia-legume symbiosis.

Symbiotic associations between legumes and rhizobia usually commence with the perception of bacterial lipochitooligosaccharides, known as Nod factors (NF), which triggers rapid cellular and molecular responses in host plants. We report here deep untargeted tandem mass spectrometry-based measurements of rapid NF-induced changes in the phosphorylation status of 13,506 phosphosites in 7739 proteins from the model legume Medicago truncatula. To place these phosphorylation changes within a biological context, quantitative phosphoproteomic and RNA measurements in wild-type plants were compared with those observed in mutants, one defective in NF perception (nfp) and one defective in downstream signal transduction events (dmi3). Our study quantified the early phosphorylation and transcription dynamics that are specifically associated with NF-signaling, confirmed a dmi3-mediated feedback loop in the pathway, and suggested "cryptic" NF-signaling pathways, some of them being also involved in the response to symbiotic arbuscular mycorrhizal fungi.

A proteogenomic survey of the Medicago truncatula genome.

Peptide sequencing by computational assignment of tandem mass spectra to a database of putative protein sequences provides an independent approach to confirming or refuting protein predictions based on large-scale DNA and RNA sequencing efforts. This use of mass spectrometrically-derived sequence data for testing and refining predicted gene models has been termed proteogenomics. We report herein the application of proteogenomic methodology to a database of 10.9 million tandem mass spectra collected over a period of two years from proteolytically generated peptides isolated from the model legume Medicago truncatula. These spectra were searched against a database of predicted M. truncatula protein sequences generated from public databases, in silico gene model predictions, and a whole-genome six-frame translation. This search identified 78,647 distinct peptide sequences, and a comparison with the publicly available proteome from the recently published M. truncatula genome supported translation of 9,843 existing gene models and identified 1,568 novel peptides suggesting corrections or additions to the current annotations. Each supporting and novel peptide was independently validated using mRNA-derived deep sequencing coverage and an overall correlation of 93% between the two data types was observed. We have additionally highlighted examples of several aspects of structural annotation for which tandem MS provides unique evidence not easily obtainable through typical DNA or RNA sequencing. Proteogenomic analysis is a valuable and unique source of information for the structural annotation of genomes and should be included in such efforts to ensure that the genome models used by biologists mirror as accurately as possible what is present in the cell.

Expression of chicken parvovirus VP2 in chicken embryo fibroblasts requires codon optimization for production of naked DNA and vectored meleagrid herpesvirus type 1 vaccines.

Meleagrid herpesvirus type 1 (MeHV-1) is an ideal vector for the expression of antigens from pathogenic avian organisms in order to generate vaccines. Chicken parvovirus (ChPV) is a widespread infectious virus that causes serious disease in chickens. It is one of the etiological agents largely suspected in causing Runting Stunting Syndrome (RSS) in chickens. Initial attempts to express the wild-type gene encoding the capsid protein VP2 of ChPV by insertion into the thymidine kinase gene of MeHV-1 were unsuccessful. However, transient expression of a codon-optimized synthetic VP2 gene cloned into the bicistronic vector pIRES2-Ds-Red2, could be demonstrated by immunocytochemical staining of transfected chicken embryo fibroblasts (CEFs). Red fluorescence could also be detected in these transfected cells since the red fluorescent protein gene is downstream from the internal ribosome entry site (IRES). Strikingly, fluorescence could not be demonstrated in cells transiently transfected with the bicistronic vector containing the wild-type or non-codon-optimized VP2 gene. Immunocytochemical staining of these cells also failed to demonstrate expression of wild-type VP2, indicating that the lack of expression was at the RNA level and the VP2 protein was not toxic to CEFs. Chickens vaccinated with a DNA vaccine consisting of the bicistronic vector containing the codon-optimized VP2 elicited a humoral immune response as measured by a VP2-specific ELISA. This VP2 codon-optimized bicistronic cassette was rescued into the MeHV-1 genome generating a vectored vaccine against ChPV disease.

Identification and characterization of the genomic termini and cleavage/packaging signals of gallid herpesvirus type 2.

Herpesvirus replication within host cells results in concatemeric genomic DNA, which is cleaved into unit-length genomes and packaged into the capsid by a complex of proteins. The sites of cleavage have been identified for many herpesviruses, and conserved signaling sequences involved in cleavage and packaging have been characterized. The cleavage/packaging motifs pac-1, pac-2, and DR1 and two distinct groups of telomeric repeat sequences (static TRS and variable TRS) have been identified. By sequencing the termini of the gallid herpesvirus type 2 (GaHV-2) strain CU-2, two different cleavage sites (classical and aberrant) have been identified. Unlike classical cleavage of human herpesvirus type 1, which occurs within the DR1 site, classical cleavage of the GaHV-2 concatemers occurs 8.5 bp upstream of the DR1 site and results in an S-terminus containing telomeric repeats. Aberrant cleavage occurs the same distance from the DR1 site and generates a telomeric S-terminus but an L-terminus lacking an a sequence. These results are consistent with previous findings in other herpesviruses and should prove useful in the future study and manipulation of the GaHV-2 genome.

Virulent Newcastle disease virus genotypes V.3, VII.2, and XIII.1.1 and their coinfections with infectious bronchitis viruses and other avian pathogens in backyard chickens in Tanzania.

Oropharyngeal (OP) and cloacal (CL) swabs from 2049 adult backyard chickens collected at 12 live bird markets, two each in Arusha, Dar es Salaam, Iringa, Mbeya, Morogoro and Tanga regions of Tanzania were screened for Newcastle disease virus (NDV) using reverse transcription real-time PCR (rRT-PCR). The virus was confirmed in 25.23% of the birds (n = 517; rRT-PCR CT ≤ 30), with the highest positivity rates observed in birds from Dar es Salaam region with higher prevalence during the dry season (September-November 2018) compared to the rainy season (January and April-May 2019). Next-generation sequencing of OP/CL samples of 20 out of 32 birds that had high amounts of viral RNAs (CT ≤ 25) resulted in the assembly of 18 complete and two partial genome sequences (15,192 bp and 15,045-15,190 bp in length, respectively) of NDV sub-genotypes V.3, VII.2 and XIII.1.1 (n = 1, 13 and 4 strains, respectively). Two birds had mixed NDV infections (V.3/VII.2 and VII.2/XIII.1.1), and nine were coinfected with viruses of families Astroviridae, Coronaviridae, Orthomyxoviridae, Picornaviridae, Pneumoviridae, and Reoviridae. Of the coinfecting viruses, complete genome sequences of two avastroviruses (a recombinant chicken astrovirus antigenic group-Aii and avian nephritis virus genogroup-5) and two infectious bronchitis viruses (a turkey coronavirus-like recombinant and a GI-19 virus) were determined. The fusion (F) protein F1/F2 cleavage sites of the Tanzanian NDVs have the consensus motifs 112 RRRKR↓F 117 (VII.2 strains) and 112 RRQKR↓F 117 (V.3 and XIII.1.1 strains) consistent with virulent virus; virulence was confirmed by intracerebral pathogenicity index scores of 1.66-1.88 in 1-day-old chicks using nine of the 20 isolates. Phylogenetically, the complete F-gene and full genome sequences regionally cluster the Tanzanian NDVs with, but distinctly from, other strains previously reported in eastern and southern African countries. These data contribute to the understanding of NDV epidemiology in Tanzania and the region.

Genome Sequences and Characterization of Chicken Astrovirus and Avian Nephritis Virus from Tanzanian Live Bird Markets.

The enteric chicken astrovirus (CAstV) and avian nephritis virus (ANV) are the type species of the genus Avastrovirus (AAstV; Astroviridae family), capable of causing considerable production losses in poultry. Using next-generation sequencing of a cloacal swab from a backyard chicken in Tanzania, we assembled genome sequences of ANV and CAstV (6918 nt and 7318 nt in length, respectively, excluding poly(A) tails, which have a typical AAstV genome architecture (5'-UTR-ORF1a-ORF1b-ORF2-'3-UTR). They are most similar to strains ck/ANV/BR/RS/6R/15 (82.72%) and ck/CAstV/PL/G059/14 (82.23%), respectively. Phylogenetic and sequence analyses of the genomes and the three open reading frames (ORFs) grouped the Tanzanian ANV and CAstV strains with Eurasian ANV-5 and CAstV-Aii viruses, respectively. Compared to other AAstVs, the Tanzanian strains have numerous amino acid variations (substitutions, insertions and deletions) in the spike region of the capsid protein. Furthermore, CAstV-A has a 4018 nt recombinant fragment in the ORF1a/1b genomic region, predicted to be from Eurasian CAstV-Bi and Bvi parental strains. These data should inform future epidemiological studies and options for AAstV diagnostics and vaccines.

Molecular Characterization of Complete Genome Sequence of an Avian Coronavirus Identified in a Backyard Chicken from Tanzania.

A complete genome sequence of an avian coronavirus (AvCoV; 27,663 bp excluding 3' poly(A) tail) was determined using nontargeted next-generation sequencing (NGS) of an oropharyngeal swab from a backyard chicken in a live bird market in Arusha, Tanzania. The open reading frames (ORFs) of the Tanzanian strain TZ/CA127/19 are organized as typical of gammaCoVs (Coronaviridae family): 5'UTR-[ORFs 1a/1b encoding replicase complex (Rep1ab) non-structural peptides nsp2-16]-[spike (S) protein]-[ORFs 3a/3b]-[small envelop (E) protein]-[membrane (M) protein]-[ORFs 4a/4c]-[ORFs 5a/5b]-[nucleocapsid (N) protein]-[ORF6b]-3'UTR. The structural (S, E, M and N) and Rep1ab proteins of TZ/CA127/19 contain features typically conserved in AvCoVs, including the cleavage sites and functional motifs in Rep1ab and S. Its genome backbone (non-spike region) is closest to Asian GI-7 and GI-19 infectious bronchitis viruses (IBVs) with 87.2-89.7% nucleotide (nt) identities, but it has a S gene closest (98.9% nt identity) to the recombinant strain ck/CN/ahysx-1/16. Its 3a, 3b E and 4c sequences are closest to the duck CoV strain DK/GD/27/14 at 99.43%, 100%, 99.65% and 99.38% nt identities, respectively. Whereas its S gene phylogenetically cluster with North American TCoVs and French guineafowl COVs, all other viral genes group monophyletically with Eurasian GI-7/GI-19 IBVs and Chinese recombinant AvCoVs. Detection of a 4445 nt-long recombinant fragment with breakpoints at positions 19,961 and 24,405 (C- and N-terminus of nsp16 and E, respectively) strongly suggested that TZ/CA127/19 acquired its genome backbone from an LX4-type (GI-19) field strain via recombination with an unknown AvCoV. This is the first report of AvCoV in Tanzania and leaves unanswered the questions of its emergence and the biological significance.

Unique Variants of Avian Coronaviruses from Indigenous Chickens in Kenya.

The avian gamma-coronavirus infectious bronchitis virus (AvCoV, IBV; Coronaviridae family) causes upper respiratory disease associated with severe economic losses in the poultry industry worldwide. Here, we report for the first time in Kenya and the Eastern African region two novel AvCoVs, designated IBV/ck/KE/1920/A374/2017 (A374/17) and AvCoV/ck/KE/1922/A376/2017 (A376/17), inadvertently discovered using random nontargeted next-generation sequencing (NGS) of cloacal swabs collected from indigenous chickens. Despite having genome organization (5'UTR-[Rep1a/1ab-S-3a-3b-E-M-4b-4c-5a-5b-N-6b]-3'UTR), canonical conservation of essential genes and size (~27.6 kb) typical of IBVs, the Kenyan isolates do not phylogenetically cluster with any genotypes of the 37 IBV lineages and 26 unique variants (UVs). Excluding the spike gene, genome sequences of A374/17 and A376/17 are only 93.1% similar to each other and 86.7-91.4% identical to genomes of other AvCoVs. All five non-spike genes of the two isolates phylogenetically cluster together and distinctly from other IBVs and turkey coronaviruses (TCoVs), including the indigenous African GI-26 viruses, suggesting a common origin of the genome backbone of the Kenyan isolates. However, isolate A376/17 contains a TCoV-like spike (S) protein coding sequence and is most similar to Asian TCoVs (84.5-85.1%) compared to other TCoVs (75.6-78.5%), whereas isolate A374/17 contains an S1 gene sequence most similar to the globally distributed lineage GI-16 (78.4-79.5%) and the Middle Eastern lineage GI-23 (79.8-80.2%) viruses. Unanswered questions include the actual origin of the Kenyan AvCoVs, the potential pathobiological significance of their genetic variations, whether they have indeed established themselves as independent variants and subsequently spread within Kenya and to the neighboring east/central African countries that have porous live poultry trade borders, and whether the live-attenuated Mass-type (lineage GI-1)-based vaccines currently used in Kenya and most of the African countries provide protection against these genetically divergent field variants.

Complete Genome Sequences of Avian Metapneumovirus Subtype B Vaccine Strains from Brazil.

Avian metapneumovirus (aMPV) causes a highly contagious upper respiratory and reproductive disease in chickens, turkeys, and ducks. Here, complete genome sequences of aMPV-B vaccine strains BR/1890/E1/19 (PL21, Nemovac; Boehringer Ingelheim Animal Health, Brazil) and BR/1891/E2/19 (1062; Hipraviar, France) were sequenced and compared with the pathogenic field strain VCO3/60616.

Dynamic equilibrium of Marek's disease genomes during in vitro serial passage.

Attenuation of Gallid herpesvirus-2 (GaHV-2), the causative agent of Marek's disease, can occur through serial passage of a virulent field isolate in avian embryo fibroblasts. In order to gain a better understanding of the genes involved in attenuation and associate observed changes in phenotype with specific genetic variations, the genomic DNA sequence of a single GaHV-2 virulent strain (648A) was determined at defined passage intervals. Biological characterization of these "interval-isolates" in chickens previously indicated that the ability to induce transient paralysis was lost by passages 40 and the ability to induce persistent neurological disease was lost after passage 80, coincident with the loss of neoplastic lesion formation. Deep sequencing of the interval-isolates allowed for a detailed cataloguing of the mutations that exist within a single passage population and the frequency with which a given mutation occurs across passages. Gross genetic alterations were identified in both novel and well-characterized genes and cis-acting regions involved in replication and cleavage/packaging. Deletions in genes encoding the virulence factors vLipase, vIL8, and RLORF4, as well as a deletion in the promoter of ICP4, appeared between passages 61 and 101. Three mutations in the virus-encoded telomerase which predominated in late passages were also identified. Overall, the frequency of mutations fluctuated greatly during serial passage and few genetic changes were absolute. This indicates that serial passage of GaHV-2 results in the generation of a collection of genomes with limited sequence heterogeneity.

Genome Sequence Variations of Infectious Bronchitis Virus Serotypes From Commercial Chickens in Mexico.

New variants of infectious bronchitis viruses (IBVs; Coronaviridae) continuously emerge despite routine vaccinations. Here, we report genome sequence variations of IBVs identified by random non-targeted next generation sequencing (NGS) of vaccine and field samples collected on FTA cards from commercial flocks in Mexico in 2019-2021. Paired-ended sequencing libraries prepared from rRNA-depleted RNAs were sequenced using Illumina MiSeq. IBV RNA was detected in 60.07% (n = 167) of the analyzed samples, from which 33 complete genome sequences were de novo assembled. The genomes are organized as 5'UTR-[Rep1a-Rep1b-S-3a-3b-E-M-4b-4c-5a-5b-N-6b]-3'UTR, except in eight sequences lacking non-structural protein genes (accessory genes) 4b, 4c, and 6b. Seventeen sequences have auxiliary S2' cleavage site located 153 residues downstream the canonically conserved primary furin-specific S1/S2 cleavage site. The sequences distinctly cluster into lineages GI-1 (Mass-type; n = 8), GI-3 (Holte/Iowa-97; n = 2), GI-9 (Arkansas-like; n = 8), GI-13 (793B; n = 14), and GI-17 (California variant; CAV; n = 1), with regional distribution in Mexico; this is the first report of the presence of 793B- and CAV-like strains in the country. Various point mutations, substitutions, insertions and deletions are present in the S1 hypervariable regions (HVRs I-III) across all 5 lineages, including in residues 38, 43, 56, 63, 66, and 69 that are critical in viral attachment to respiratory tract tissues. Nine intra-/inter-lineage recombination events are present in the S proteins of three Mass-type sequences, two each of Holte/Iowa-97 and Ark-like sequence, and one each of 793B-like and CAV-like sequences. This study demonstrates the feasibility of FTA cards as an attractive, adoptable low-cost sampling option for untargeted discovery of avian viral agents in field-collected clinical samples. Collectively, our data points to co-circulation of multiple distinct IBVs in Mexican commercial flocks, underscoring the need for active surveillance and a review of IBV vaccines currently used in Mexico and the larger Latin America region.

Detection and Genome Sequence Analysis of Avian Metapneumovirus Subtype A Viruses Circulating in Commercial Chicken Flocks in Mexico.

Avian metapneumoviruses (aMPV subtypes A-D) are respiratory and reproductive pathogens of poultry. Since aMPV-A was initially reported in Mexico in 2014, there have been no additional reports of its detection in the country. Using nontargeted next-generation sequencing (NGS) of FTA card-spotted respiratory samples from commercial chickens in Mexico, seven full genome sequences of aMPV-A (lengths of 13,288-13,381 nucleotides) were de novo assembled. Additionally, complete coding sequences of genes N (n = 2), P and M (n = 7 each), F and L (n = 1 each), M2 (n = 6), SH (n = 5) and G (n = 2) were reference-based assembled from another seven samples. The Mexican isolates phylogenetically group with, but in a distinct clade separate from, other aMPV-A strains. The genome and G-gene nt sequences of the Mexican aMPVs are closest to strain UK/8544/06 (97.22-97.47% and 95.07-95.83%, respectively). Various amino acid variations distinguish the Mexican isolates from each other, and other aMPV-A strains, most of which are in the G (n = 38), F (n = 12), and L (n = 19) proteins. Using our sequence data and publicly available aMPV-A data, we revised a previously published rRT-PCR test, which resulted in different cycling and amplification conditions for aMPV-A to make it more compatible with other commonly used rRT-PCR diagnostic cycling conditions. This is the first comprehensive sequence analysis of aMPVs in Mexico and demonstrates the value of nontargeted NGS to identify pathogens where targeted virus surveillance is likely not routinely performed.