Evaluation of the effectiveness of single nucleotide polymorphisms compared to microsatellite markers for parentage verification in Moroccan horses.

The International Society for Animal Genetics (ISAG) currently advocates for a transition towards single nucleotide polymorphism (SNP) markers as a potential alternative for equine parentage verification. To ascertain the efficacy of this transition, it is imperative to evaluate the performance of parentage testing using SNPs in juxtaposition with short tandem repeats (STRs). As per ISAG's recommendation, we used an equine genotyping-by-sequencing panel with 144 SNPs for this purpose. Equine parentage is currently realized using 16 microsatellites (STRs) excluding the LEX3 marker. In this study, 1074 horses were genotyped using the 144 SNPs panel, including 432 foals, 414 mares, and 228 stallions, from five different breeds: 293 Arabians, 167 Barbs, 189 Thoroughbreds, 73 Anglo-Arabians, and 352 Arabian-Barbs. As a result, two SNPs markers were eliminated from the panel system due to inconsistent amplification across all examined individuals leaving 142 SNPs markers for analysis. A comparative analysis between SNPs and STRs markers revealed that the mean expected heterozygosity was 0.457 for SNPs and 0.76 for STRs, while the mean observed heterozygosity stood at 0.472 for SNPs and 0.72 for STRs. Furthermore, the probability of identity was calculated to be 5.722 × 10-57 for SNPs and 1.25 × 10-15 for STRs markers. In alignment with the Hardy-Weinberg equilibrium in polyploids test, 110 out of the total SNPs were consistent with the Hardy-Weinberg equilibrium in polyploids test (p > 0.05). Employing both SNPs and STRs markers, the mean polymorphic information content was discerned to be 0.351 for SNPs and 0.72 for STRs. The cumulative exclusion probabilities for SNP markers exceeded 99.99%, indicating that the 142 SNPs panel might be adequate for parentage testing. In contrast, when utilizing STRs markers, the combined average exclusion probabilities for one and both parents were determined to be 99.8% and 99.9%, respectively. Our comprehensive study underscores the potential of SNPs in equine parentage verification, especially when compared to STRs in terms of exclusion probabilities. As a corollary, the application of SNPs for parentage verification and identification can significantly contribute to the conservation initiative for the five Moroccan horse breeds. Nonetheless, further research is required to address and replace the deficient SNPs within the panel.

Multiplex fluorescent amplification-refractory mutation system PCR method for the detection of 10 genetic defects in Holstein cattle and its comparison with the KASP genotyping assay.

The common deleterious genetic defects in Holstein cattle include haplotypes 1-6 (HH1-HH6), haplotypes for cholesterol deficiency (HCD), bovine leukocyte adhesion deficiency (BLAD), complex vertebral malformation (CVM) and brachyspina syndrome (BS). Recessive inheritance patterns of these genetic defects permit the carriers to function normally, but homozygous recessive genotypes cause embryo loss or neonatal death. Therefore, rapid detection of the carriers is essential to manage these genetic defects. This study was conducted to develop a single-tube multiplex fluorescent amplification-refractory mutation system (mf-ARMS) PCR method for efficient genotyping of these 10 genetic defects and to compare its efficiency with the kompetitive allele specific PCR (KASP) genotyping assay. The mf-ARMS PCR method introduced 10 sets of tri-primers optimized with additional mismatches in the 3' end of wild and mutant-specific primers, size differentiation between wild and mutant-specific primers, fluorescent labeling of universal primers, adjustment of annealing temperatures and optimization of primer concentrations. The genotyping of 484 Holstein cows resulted in 16.12% carriers with at least one genetic defect, while no homozygous recessive genotype was detected. This study found carrier frequencies ranging from 0.0% (HH6) to 3.72% (HH3) for individual defects. The mf-ARMS PCR method demonstrated improved detection, time and cost efficiency compared with the KASP method for these defects. Therefore, the application of mf-ARMS PCR for genotyping Holstein cattle is anticipated to decrease the frequency of lethal alleles and limit the transmission of these genetic defects.

Genetic characterization of a captive marmoset (Callithrix jacchus) colony using genotype-by-sequencing.

The marmoset is a fundamental nonhuman primate model for the study of aging, neurobiology, and many other topics. Genetic management of captive marmoset colonies is complicated by frequent chimerism in the blood and other tissues, a lack of tools to enable cost-effective, genome-wide interrogation of variation, and historic mergers and migrations of animals between colonies. We implemented genotype-by-sequencing (GBS) of hair follicle derived DNA (a minimally chimeric DNA source) of 82 marmosets housed at the Southwest National Primate Research Center (SNPRC). Our primary goals were the genetic characterization of our marmoset population for pedigree verification and colony management and to inform the scientific community of the functional genetic makeup of this valuable resource. We used the GBS data to reconstruct the genetic legacy of recent mergers between colonies, to identify genetically related animals whose relationships were previously unknown due to incomplete pedigree information, and to show that animals in the SNPRC colony appear to exhibit low levels of inbreeding. Of the >99,000 single-nucleotide variants (SNVs) that we characterized, >9800 are located within gene regions known to harbor pathogenic variants of clinical significance in humans. Overall, we show the combination of low-resolution (sparse) genotyping using hair follicle DNA is a powerful strategy for the genetic management of captive marmoset colonies and for identifying potential SNVs for the development of biomedical research models.

Development of a high-resolution typing method for SLA-3, swine MHC class I antigen 3.

We developed a high-resolution and comprehensive typing method for swine leukocyte antigen 3 (SLA-3), an MHC class I gene, employing locus-specific genomic PCR followed by subsequent direct sequencing. A total of 292 individuals from nine pure, one cross-breed and six cell lines were successfully typed. A total of 21 SLA-3 alleles were identified, of which four were found to be novel alleles. However, the allelic diversity of SLA-3 was lower than that of previously reported class I genes, SLA-1 and -2. More SLA-3 alleles were observed in the Landrace and Yorkshire breeds than the other breeds. SLA-3*04:01 was identified in seven out of nine breeds and was the most widely distributed allele across all breeds. Therefore, the typing method reported in this study completes our efforts to develop high-resolution typing methods for major SLA molecules, facilitating the combined analysis of major SLA genes from field samples, which is important to understand the relationship between the adaptive immune responses against pathogens and the immunogenetic makeup of an individual.

Sero- and apx-typing of German Actinobacillus pleuropneumoniae field isolates from 2010 to 2019 reveals a predominance of serovar 2 with regular apx-profile.

Serotyping is the most common method to characterize field isolates of Actinobacillus (A.) pleuropneumoniae, the etiological agent of porcine pleuropneumonia. Based on serology, many farms seem to be infected and antibodies against a wide variety of serovars are detectable, but, so far it is unknown to what degree respective serovars contribute to outbreaks of clinical manifest disease. In this study, 213 German A. pleuropneumoniae field isolates retrieved for diagnostic purposes from outbreaks of porcine pleuropneumonia between 2010 and 2019 were genetically serotyped and analyzed regarding their apx-toxin gene profile using molecular methods. Serotyping revealed a prominent role of serovar 2 in clinical cases (64% of all isolates) and an increase in the detection of this serovar since 2010 in German isolates. Serovar 9/11 followed as the second most frequent serovar with about 15% of the isolates. Furthermore, very recently described serovars 16 (n = 2) and 18 (n = 8) were detected. Most isolates (93.4%) showed apx-profiles typical for the respective serovar. However, this does not hold true for isolates of serovar 18, as 75% (n = 6) of all isolates of this serovar deviated uniformly from the "typical" apx-gene profile of the reference strain 7311555. Notably, isolates from systemic lesions such as joints or meninges did not harbor the complete apxICABD operon which is considered typical for highly virulent strains. Furthermore, the extremely low occurrence (n = 1) of NAD independent (biovar II) isolates in German A. pleuropneumoniae was evident in our collection of clinical isolates.

Investigation of giardiasis in captive animals in zoological gardens with strain typing of assemblages in China.

Giardia duodenalis is a common zoonotic intestinal pathogen. It has been increasingly reported in humans and animals; however, genotyping information for G. duodenalis in captive animals is still limited. This study was conducted to assess the prevalence and multilocus genotyping of G. duodenalis in captive animals in zoological gardens in Shanghai, China. A total of 678 fresh fecal samples were randomly collected from captive animals including non-human primates (NHPs) (n = 190), herbivores (n = 190), carnivores (n = 151), birds (n = 138) and reptiles (n = 9) in a zoo and were examined for the presence of G. duodenalis using nested polymerase chain reaction (nested PCR). All G. duodenalis positive samples were assayed with PCR followed by sequencing at ß-giardin (bg), glutamate dehydrogenase (gdh) and triose phosphate isomerase (tpi) genes. In this study, 42 specimens (6.2%) were tested G. duodenalis-positive of the 678 fecal samples examined based on a single locus. A total of 30 (4.4%), 30 (4.4%) and 22 (3.2%) specimens were successfully amplified and sequenced at gdh, tpi and bg loci, respectively. Assemblages A and B were identified with assemblage B dominating in NHPs. Sequence analysis demonstrated that one, two and five new isolates were identified at bg, gdh and tpi loci. DNA sequences and new assemblage-subtypes of zoonotic G. duodenalis assemblages A and B were identified in the current study. Our data indicate the occurrence and molecular diversity of G. duodenalis and the potential zoonotic transmission in captive animals in China.

Study on the concordance between different SNP-genotyping platforms in sheep.

Different SNP genotyping technologies are commonly used in multiple studies to perform QTL detection, genotype imputation, and genomic predictions. Therefore, genotyping errors cannot be ignored, as they can reduce the accuracy of different procedures applied in genomic selection, such as genomic imputation, genomic predictions, and false-positive results in genome-wide association studies. Currently, whole-genome resequencing (WGR) also offers the potential for variant calling analysis and high-throughput genotyping. WGR might overshadow array-based genotyping technologies due to the larger amount and precision of the genomic information provided; however, its comparatively higher price per individual still limits its use in larger populations. Thus, the objective of this work was to evaluate the accuracy of the two most popular SNP-chip technologies, namely, Affymetrix and Illumina, for high-throughput genotyping in sheep considering high-coverage WGR datasets as references. Analyses were performed using two reference sheep genome assemblies, the popular Oar_v3.1 reference genome and the latest available version Oar_rambouillet_v1.0. Our results demonstrate that the genotypes from both platforms are suggested to have high concordance rates with the genotypes determined from reference WGR datasets (96.59% and 99.51% for Affymetrix and Illumina technologies, respectively). The concordance results provided in the current study can pinpoint low reproducible markers across multiple platforms used for sheep genotyping data. Comparing results using two reference genome assemblies also informs how genome assembly quality can influence genotype concordance rates among different genotyping platforms. Moreover, we describe an efficient pipeline to test the reliability of markers included in sheep SNP-chip panels against WGR datasets available on public databases. This pipeline may be helpful for discarding low-reliability markers before exploiting genomic information for gene mapping analyses or genomic prediction.

Screening and detection of chromosomal copy number alterations in the domestic horse using SNP-array genotyping data.

Chromosomal abnormalities are a common cause of infertility in horses. However, they are difficult to detect using automated methods. Here, we propose a simple methodology based on single nucleotide polymorphism (SNP)-array data that allows us to detect the main chromosomal abnormalities in horses in a single procedure. As proof of concept, we were able to detect chromosomal abnormalities in 33 out of 268 individuals, including monosomies, chimerisms, and male and female sex-reversions, by analyzing the raw signal intensity produced by an SNP array-based genotyping platform. We also demonstrated that the procedure is not affected by the SNP density of the array employed or by the inbreeding level of the individuals. Finally, the methodology proposed in this study could be performed in an open bioinformatic environment, thus permitting its integration as a flexible screening tool in diagnostic laboratories and genomic breeding programs.

The equine graying with age mutation of the STX17 gene: A copy number study using droplet digital PCR reveals a new pattern.

The equine graying with age causative mutation in the syntaxin-17 gene (STX17) has been known for over a decade, but proper genotyping of this variant remains challenging due to its molecular character (4.6-kb tandem duplication). Precise information on gray mutation status is important for horse breeders and veterinarians, since gray homozygous horses are more prone to developing aggressive melanoma tumors than heterozygotes. Since recent studies have confirmed that droplet digital PCR is a valuable technique for copy number analysis, we decided to investigate whether this method can be used for accurate genotyping of the horse graying-related variant and established the copy numbers of the 4.6-kb fragment in the available cohort (n = 75) of gray and nongray horses of various breeds. Surprisingly, we found that our STX17 genotype results varied from what has been previously published, suggesting that gray phenotype is associated with the presence of six (GG) or four (Gg) copies of studied region. All the examined nongray horses (gg) have the two copies of these fragments. This new pattern and its inheritance were also confirmed by an analysis conducted for the Polish Warmblood horse family. We noted no further copy number variation in the entire tested samples set. Our study confirmed the usefulness and accuracy of droplet digital PCR for genotyping STX17 gene variant. Further studies on a broader range of materials are needed to fully understand the origin and molecular structure of the graying causative mutation in the horse STX17.

A genome-wide association study for the number of teats in European rabbits (Oryctolagus cuniculus) identifies several candidate genes affecting this trait.

In the European rabbit (Oryctolagus cuniculus), a polytocous livestock species, the number of teats indirectly impacts the doe reproduction efficiency and, in turn, the sustainable production of rabbit meat. In this study, we carried out a genome-wide association study (GWAS) for the total number of teats in 247 Italian White does included in the Italian White rabbit breed selection program, by applying a selective genotyping approach. Does had either 8 (n = 121) or 10 teats (n = 126). All rabbits were genotyped with the Affymetrix Axiom OrcunSNP Array. Genomic data from the two extreme groups of rabbits were also analysed with the single-marker fixation index statistic and combined with the GWAS results. The GWAS identified 50 significant SNPs and the fixation index analysis identified a total of 20 SNPs that trespassed the 99.98th percentile threshold, 19 of which confirmed the GWAS results. The most significant SNP (P = 4.31 × 10-11 ) was located on OCU1, close to the NUDT2 gene, a breast carcinoma cells proliferation promoter. Another significant SNP identified as candidate gene NR6A1, which is well known to play an important role in affecting the correlated number of vertebrae in pigs. Other significant markers were close to candidate genes involved in determining body length in mice. Markers associated with increased number of teats could be included in selection programmes to speed up the improvement for this trait in rabbit lines that need to increase maternal performances.

Identification of eQTLs and sQTLs associated with meat quality in beef.

BACKGROUND: Transcription has a substantial genetic control and genetic dissection of gene expression could help us understand the genetic architecture of complex phenotypes such as meat quality in cattle. The objectives of the present research were: 1) to perform eQTL and sQTL mapping analyses for meat quality traits in longissimus dorsi muscle; 2) to uncover genes whose expression is influenced by local or distant genetic variation; 3) to identify expression and splicing hot spots; and 4) to uncover genomic regions affecting the expression of multiple genes. RESULTS: Eighty steers were selected for phenotyping, genotyping and RNA-seq evaluation. A panel of traits related to meat quality was recorded in longissimus dorsi muscle. Information on 112,042 SNPs and expression data on 8588 autosomal genes and 87,770 exons from 8467 genes were included in an expression and splicing quantitative trait loci (QTL) mapping (eQTL and sQTL, respectively). A gene, exon and isoform differential expression analysis previously carried out in this population identified 1352 genes, referred to as DEG, as explaining part of the variability associated with meat quality traits. The eQTL and sQTL mapping was performed using a linear regression model in the R package Matrix eQTL. Genotype and year of birth were included as fixed effects, and population structure was accounted for by including as a covariate the first PC from a PCA analysis on genotypic data. The identified QTLs were classified as cis or trans using 1 Mb as the maximum distance between the associated SNP and the gene being analyzed. A total of 8377 eQTLs were identified, including 75.6% trans, 10.4% cis, 12.5% DEG trans and 1.5% DEG cis; while 11,929 sQTLs were uncovered: 66.1% trans, 16.9% DEG trans, 14% cis and 3% DEG cis. Twenty-seven expression master regulators and 13 splicing master regulators were identified and were classified as membrane-associated or cytoskeletal proteins, transcription factors or DNA methylases. These genes could control the expression of other genes through cell signaling or by a direct transcriptional activation/repression mechanism. CONCLUSION: In the present analysis, we show that eQTL and sQTL mapping makes possible positional identification of gene and isoform expression regulators.

SNP discovery and population structure analysis in Lassi and Marecha camel breeds using a genotyping by sequencing method.

Pakistani camels have been classified socio-geographically into 20 breeds, but they have not yet been subjected to substantial selective pressures and the genetic basis for these breeds is not understood. However, it should be possible to distinguish them by use of molecular data. This study investigated the genetic diversity and population structure within and between two major Pakistani camel breeds, Marecha and Lassi. As no SNP array is currently available, we first identified 63 619 SNPs using a genotyping by sequencing approach. After quality control, a panel of 36 926 SNPs was used in the analysis. Population structure was investigated with a principal coordinate analysis as well as a cluster analysis using NetView, and multilocus heterozygosity analysis to explore between- and within-breed genetic variation. In addition, between-breed variation was explored using the fixation index, FST . We also compared relationship matrices computed using the VanRaden SNP-based method and a method developed specifically for genotyping by sequencing data. Among the two camel breeds, Lassi showed a lower level of genetic diversity whereas Marecha showed a higher level. As a genotyping platform has not yet been developed for the camel, the SNPs discovered in this study will be useful in future genetic studies in camels.

Genotyping more cows increases genetic gain and reduces rate of true inbreeding in a dairy cattle breeding scheme using female reproductive technologies.

Both small dairy cattle populations and dairy cattle populations with a low level of linkage disequilibrium (LD) suffer from low reliability of genomic prediction. In this study, we investigated whether adding more genotyped cows to the reference population influences the rate of genetic gain and rate of inbreeding by affecting the reliability. A standard breeding program with a large reference population and high LD, which mimicked a breeding program for Danish Holstein population, was simulated as a reference. A Danish Jersey population with a small reference population and high LD and a Red Dairy Cattle population with a large reference population and low LD were also simulated. Two additional breeding programs were simulated for Danish Jersey and Red Dairy Cattle populations, where 2,000 additional genotyped cows were included in the population for genomic selection. All 5 simulated breeding programs were initiated by a founder population to generate LD resembling the real LD pattern, followed by a 20-yr conventional progeny-testing scheme with 1,000 or 10,000 genotyped progeny-tested bulls and a 10-yr genomic selection scheme with or without 2,000 additional genotyped cows. Evaluation criteria were annual monetary genetic gain and rate of true inbreeding. Our results showed that adding more genotyped cows to the reference in dairy cattle populations has the potential to increase genetic gain and reduce the rate of inbreeding, regardless of reference population size and level of LD. However, it is still not possible to reach the same genetic gain as in the simulated Danish Holstein population with either a small reference population or low LD. Our results also showed that in a small reference population with high LD, it is difficult to manage inbreeding because of lower accuracy compared with the simulated Danish Holstein population and a smaller number of relevant families to select from. Therefore, breeding strategies need to be chosen to match population size and structure. The rate of true inbreeding is always underestimated by pedigree inbreeding and even more in genomic breeding programs, indicating that some forms of genome-wide inbreeding, instead of pedigree-based inbreeding, should be used to monitor inbreeding when genomic selection is implemented.

ACAA2 and FASN polymorphisms affect the fatty acid profile of Chios sheep milk.

The objective of the research reported in this research communication was the identification and association of single nucleotide polymorphisms (SNP) in the ovine DGAT1, FASN, SCD1 and ACAA2 genes with milk fat percentage and fatty acid (FA) content. Three consecutive monthly milk samplings were obtained from a total of 429 purebred Chios ewes during mid-lactation. Genotypic data were jointly analyzed with 1184 fat content and 37 718 FA percentage records using mixed models. The 3' untranslated region (UTR) of the DGAT1 gene and the 5' and 3'UTRs of the SCD1 gene appeared to be monomorphic. The FASN g.14777C>T SNP on exon 31 was associated with C13:0 and the ACAA2 g.2982T>C SNP on the 3'UTR was associated with C9:0, C11:0, C12:1 cis-9, C13:0 and the ω6/ω3 index, while fat percentage was not affected by the identified SNPs. The results could be useful for breeding programs aiming to improve the quality and nutritional value of ovine milk.

Effects of FABP4 variation on milk fatty-acid composition for dairy cattle grazed on pasture in late lactation.

This research communication describes associations between variation in the fatty acid binding protein 4 gene (FABP4) and milk fat composition in New Zealand Holstein-Friesian × Jersey cross dairy cows. After correcting for the effect of the amino acid substitution p.K232A in diacylglycerol O-acyltransferase 1 (DGAT1), which is associated with variation in many milk fatty acid (FA) component levels, the effect of FABP4 c.328A/G on milk FA levels was typically small. For the five genotypes analysed, the AB cows produced more medium-chain fatty acids than CC cows (P < 0.05), and more C14:0 FA than AA and AC cows (P < 0.05). The AA and AC cows produced less C22:0 FA (P < 0.01) than the BC cows, and the AC cows produced more C24:0 FA (P < 0.05) than was produced by the BC cows. Cows of genotype CC produce more long-chain fatty acids than cows of genotype BC (P < 0.05).

Equine Genotyping Arrays.

High-quality genomic tools have been integral in understanding genomic architecture and function in the modern-day horse. The equine genetics community has a long tradition of pooling resources to develop genomic tools. Since the equine genome was sequenced in 2006, several iterations of high throughput genotyping arrays have been developed and released, enabling rapid and cost-effective genotyping. This review highlights the design considerations of each iteration, focusing on data available during development and outlining considerations in selecting the genetic variants included on each array. Additionally, we outline recent applications of equine genotyping arrays as well as future prospects and applications.

A seminested RT-PCR for molecular genotyping of the Brazilian BR-I Infectious Bronchitis Virus Strain (GI-11).

Infectious bronchitis (IB) is one of the avian diseases with the greatest impact on poultry farming worldwide. In Brazil, strain BR-I (GI-11) is the most prevalent in poultry flocks. The present study aimed to develop a seminested RT-PCR assay specific for the diagnosis of BR-I IBV in Brazilian samples, targeting subunit 1 of the S gene. The detection limit of this assay was 10 copies of the IBV genome. In this study, 62.24% of 572 organ pools from the 5 regions of Brazil tested positive in a 3'UTR screening, and 84.83% were typed as BR-I IBV. BR-I was detected in the respiratory, digestive and urogenital tracts in pooled samples from all Brazilian geographical regions and in all the breeding systems analyzed. Specificity and sensitivity tests as well as phylogenetic analysis successfully confirmed the expected clustering of the sequences detected by this assay with the BR-I (GI-11) group. The nested PCR described in this study represents a suitable and valuable tool in the diagnosis, epidemiology, monitoring and vaccination decisions of IBV.

Differential expression and characterization of ATP1A1 exon17 gene by high resolution melting analysis and RT-PCR in Indian goats.

The investigation was carried out to analyse the genetic polymorphism and gene expression of ATP1A1 gene in four different Indian goat breeds by using high resolution melting (HRM) and real time-PCR. ATPase is electro-genic ion pump which is maintains the balance of sodium and potassium ions in animal cells. The transport of Na+& K+ is variable at cellular level during extreme hot period. Therefore, susceptible and tolerant animals were selected based on the physiological responses during hot period. Blood samples were collected from individuals, DNA was isolated. The 300 bp fragment of ATP1A1 gene was amplified by PCR and HRM genotyping was performed. The melting curves were analysed, differential temperature-shift plot showed three different genotypes in all the analysed samples. Out of the 135 samples, the distribution percentages were 55.56% (AA/blue), 33.33% (AC/red) and 11.11% (CC/green). The sequence variation revealed a SNP at 143rd position (A>C). The nucleotide diversity was 0.695 ± 0.403, 0.732 ± 0.424, 0.662 ± 0.433 and 0.687 ± 0.398 in Barbari, Jamunapari, Jakharna and Sirohi, respectively. The respiration rate (RR) was significantly different (P < 0.05) between AA and AC (t = 1.875, df = 38) genotype and heart rate (HR) was significantly different (P < 0.05) between AA and CC genotype. The relative expression pattern of ATP1A1 in SNP variants and non-variants animal tissues showed 19.09 and 6.93 fold higher than control (non-variant), respectively. Jamunapari showed higher fold value of ATP1A1 gene in comparison to Barbari, Jakharna and Sirohi. However, the heat stress-susceptible phenotype had significantly higher gene expression than stress-tolerant in all the breeds. The variation may be used as a marker for selection on the basis of physiological parameters and expression of ATP1A1 gene in goats indicating the specificity of expression in each tissue.

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.

Comparative genomic analyses of Mycoplasma synoviae vaccine strain MS-H and its wild-type parent strain 86079/7NS: implications for the identification of virulence factors and applications in diagnosis of M. synoviae.

Mycoplasma synoviae is an economically important avian pathogen worldwide, causing respiratory disease, infectious synovitis, airsacculitis and eggshell apex abnormalities in commercial chickens. Despite the widespread use of MS-H as a live attenuated vaccine over the past two decades, the precise molecular basis for loss of virulence in this vaccine is not yet fully understood. To address this, the whole genome sequence of the vaccine parent strain, 86079/7NS, was obtained and compared to that of the MS-H vaccine. Except for the vlhA expressed region, both genomes were nearly identical. Thirty-two single nucleotide polymorphisms (SNPs) were identified in MS-H, including 11 non-synonymous mutations that were predicted, by bioinformatics analysis, to have changed the secondary structure of the deduced proteins. One of these mutations caused truncation of the oppF-1 gene, which encodes the ATP-binding protein of an oligopeptide permease transporter. Overall, the attenuation of MS-H strain may be caused by the cumulative and complex effects of several mutations. The SNPs identified in MS-H were further analyzed by comparing the MS-H and 86079/7NS sequences with the strains WVU-1853 and MS53. In the genomic regions conserved between all strains, 30 SNPs were found to be unique to MS-H lineage. These results have provided a foundation for developing novel biomarkers for the detection of virulence in M. synoviae and also for designing new genotyping assays for discrimination of MS-H from field strains.