Search

We report whole-genome sequencing of influenza A virus (IAV) with 100% diagnostic sensitivity and results available in <24-48 h using amplicon-based nanopore sequencing technology (MinION) on clinical material from wild waterfowl ( n = 19), commercial poultry ( n = 4), and swine ( n = 3). All 8 gene segments of IAV including those from 14 of the 18 recognized hemagglutinin subtypes and 9 of the 11 neuraminidase subtypes were amplified in their entirety at >500× coverage from each of 16 reference virus isolates evaluated. Subgenomic viral sequences obtained in 3 cases using Sanger sequencing as the reference standard were identical to those obtained when sequenced using the MinION approach. An inter-laboratory comparison demonstrated reproducibility when comparing 2 independent laboratories at ≥99.8% across the entirety of the IAV genomes sequenced.

The Longlin goat is one of the most valuable livestock species in Guangxi Autonomous Region of China, but its genomic diversity and selective signals are not clearly elucidated. Here we compared 20 genomes of Longlin goat to 66 genomes of other seven goat breeds worldwide to analyze patterns of Longlin goat genetic variation. We found the lowest linkage disequilibrium at the large distances between SNPs associated with the highest effective population size in the recent generations ago in Longlin goat. The eight goat breeds could be divided into Euro-African and East Asian goat population. Interestingly, like East Asian taurine, the same two migration phases might have occurred in the history of East Asian goat. More importantly, we identified selective signals implicated in immune resistance to disease, especially for skin disease, in Longlin goat. Our findings will not only help understand the evolutionary history and breed characteristic but can provide valuable resources for conservation of germplasm resources and implementation of crossbreeding programs., (Copyright © 2020 Elsevier B.V. All rights reserved.)

The laboratory identification of antibacterial resistance is a cornerstone of infectious disease medicine. In vitro antimicrobial susceptibility testing has long been based on the growth response of organisms in pure culture to a defined concentration of antimicrobial agents. By comparing individual isolates to wild-type susceptibility patterns, strains with acquired resistance can be identified. Acquired resistance can also be detected genetically. After many decades of research, the inventory of genes underlying antimicrobial resistance is well known for several pathogenic genera including zoonotic enteric organisms such as Salmonella and Campylobacter and continues to grow substantially for others. With the decline in costs for large scale DNA sequencing, it is now practicable to characterize bacteria using whole genome sequencing, including the carriage of resistance genes in individual microorganisms and those present in complex biological samples. With genomics, we can generate comprehensive, detailed information on the bacterium, the mechanisms of antibiotic resistance, clues to its source, and the nature of mobile DNA elements by which resistance spreads. These developments point to a new paradigm for antimicrobial resistance detection and tracking for both clinical and public health purposes., (© Published 2020. This article is a U.S. Government work and is in the public domain in the USA.)

The genomic surveillance of porcine reproductive and respiratory syndrome virus (PRRSV) is based on sequencing of the ORF5 gene of the virus, which covers only 4% of the entire viral genome. It is expected that PRRSV whole-genome sequencing (WGS) will improve PRRSV genomic data and allow better understanding of clinical discrepancies observed in the field when using ORF5 sequencing. Our main objective was to implement an efficient method for WGS of PRRSV from clinical samples. The viral genome was purified using a poly(A)-tail viral genome purification method and sequenced using Illumina technology. We tested 149 PRRSV-positive samples: 80 sera, 33 lungs, 33 pools of tissues, 2 oral fluids, and 1 processing fluid (i.e., castration liquid). Overall, WGS of 67.1% of PRRSV-positive cases was successful. The viral load, in particular for tissues, had a major impact on the PRRSV WGS success rate. Serum was the most efficient type of sample to conduct PRRSV WGS poly(A)-tail assays, with a success rate of 76.3%, and this result can be explained by improved sequencing reads dispersion matching throughout the entire viral genome. WGS was unsuccessful for all pools of tissue and lung samples with Cq values > 26.5, whereas it could still be successful with sera at Cq ≤ 34.1. Evaluation of results of highly qualified personnel confirmed that laboratory skills could affect PRRSV WGS efficiency. Oral fluid samples seem very promising and merit further investigation because, with only 2 samples of low viral load (Cq = 28.8, 32.8), PRRSV WGS was successful.

Epidemics of H3N8 and H3N2 influenza A viruses (IAVs) in dogs, along with recognition of spillover infections from IAV strains typically found in humans or other animals, have emphasized the importance of efficient laboratory testing. Given the lack of active IAV surveillance or immunization requirements for dogs, cats, or horses imported into the United States, serotype prediction and whole-genome sequencing of positive specimens detected at veterinary diagnostic laboratories are also needed. The conserved sequences at the ends of the viral genome segments facilitate universal amplification of all segments of viral genomes directly from respiratory specimens. Although several methods for genomic analysis have been reported, no optimization focusing on companion animal strains has been described, to our knowledge. We compared 2 sets of published universal amplification primers using 26 IAV-positive specimens from dogs, horses, and a cat. Libraries prepared from the resulting amplicons were sequenced using Illumina chemistry, and reference-based assemblies were generated from the data produced by both methods. Although both methods produced high-quality data, coverage profiles and base calling differed between the 2 methods. The sequence data were also used to identify the subtype of the IAV strains sequenced and then compared to standard PCR assays for neuraminidase types N2 and N8.

Protozoan parasite isolation and purification are laborious and time-consuming processes required for high quality genomic DNA used in whole genome sequencing. The objective of this study was to capture whole Theileria parva genomes directly from cell cultures and blood samples using RNA baits. Cell culture material was bait captured or sequenced directly, while blood samples were all captured. Baits had variable success in capturing T. parva genomes from blood samples but were successful in cell cultures. Genome mapping uncovered extensive host contamination in blood samples compared to cell cultures. Captured cell cultures had over 81 fold coverage for the reference genome compared to 0-33 fold for blood samples. Results indicate that baits are specific to T. parva, are a good alternative to conventional methods and thus ideal for genomic studies. This study also reports the first whole genome sequencing of South African T. parva., (Copyright © 2020 Elsevier Inc. All rights reserved.)

Salmon pancreas disease virus, more commonly known as salmonid alphavirus (SAV), is a single-stranded positive sense RNA virus and the causative agent of pancreas disease and sleeping disease in salmonids. In this study, a unique strain of SAV previously isolated from ballan wrasse was subjected to whole genome sequencing using nanopore sequencing. In order to accurately examine the evolutionary history of this strain in comparison to other SAV strains, a partitioned phylogenetic analysis was performed to account for variation in the rate of evolution for both individual genes and codon positions. Partitioning the genome alignments almost doubled the observed branch lengths in the phylogenetic tree when compared to the more common approach of applying one model of substitution across the genome and significantly increased the statistical fit of the best-fitting models of nucleotide substitution. Based on the genomic data, a valid case can be made for the viral strain examined in this study to be considered a new SAV genotype. In addition, this study adds to a growing number of studies in which SAV has been found to infect non-salmonid fish, and as such we have suggested that the viral species name be amended to the more inclusive 'piscine alphavirus'.

Shewanella algae is not only the most commonly reported species in Shewanella human infections but also capable to inhabit a wide variety of habitats. Although there is evidence that quorum sensing is associated with bacterial adaptation to changing environmental conditions, little is known of the quorum sensing system in S. algae. In this study, we conducted the whole genome sequencing of S. algae strains and applied comparative genomics to reveal the core genome. Genes related to the quorum sensing system were identified by integrated bioinformatics analysis. S. algae harbor genes involved in all three main types of autoinducer systems. This study provides insights into the quorum sensing systems in S. algae, which might be valuable in the future study of cell behavior in S. algae., (Copyright © 2020 Elsevier Ltd. All rights reserved.)

Large phenotypic differences have been observed between Tongcheng and Large White pigs. However, little is known about their genetic basis. This study performed a genome-wide comparison of CNVs between Tongcheng and Large White pigs using genome sequencing data. By combining the advantages of three different strategies (read depth, paired-end mapping and split read), we detected in total 18 687 CNVs that covered approximately 3.5% of the pig genome length for Tongcheng and Large White pigs. We identified 1864 breed-stratified CNVs (top 10%) by performing V ST statistics. Functional enrichment analyses for genes located in breed-stratified CNVs were found to be involved in pigmentation, behavior, immune system and reproductive processes, which coincide with phenotypic differences between the two breeds. Using a systematic analysis of the genome and transcriptome data, we further identified four novel breed-differential CNVs on the functional genes (disease-resistant, DCUN1D2 and SPARCL1; lipid metabolism, PLEKHA2 and SLCO1A2). Subsequent PCR validation confirmed their accurate breakpoint positions in 33 Tongcheng pigs and 33 Large White pigs. This study provides essential information on differential CNVs for further research on the genetic basis of phenotypic differences between Tongcheng and Large White pigs., (© 2020 Stichting International Foundation for Animal Genetics.)

Campylobacter is among the most common causes of gastroenteritis worldwide. Campylobacter jejuni and Campylobacter coli are the most common species causing human disease. DNA sequence-based methods for strain characterization have focused largely on C. jejuni , responsible for 80 to 90% of infections, meaning that C. coli epidemiology has lagged behind. Here, we have analyzed the genome of 450 C. coli isolates to determine genetic markers that can discriminate isolates sampled from 3 major reservoir hosts (chickens, cattle, and pigs). These markers then were applied to identify the source of infection of 147 C. coli strains from French clinical cases. Using STRUCTURE software, 259 potential host-segregating markers were revealed by probabilistic characterization of single-nucleotide polymorphism (SNP) frequency variation in strain collections from three different hosts. These SNPs were found in 41 genes or intergenic regions, mostly coding for proteins involved in motility and membrane functions. Source attribution of clinical isolates based on the differential presence of these markers confirmed chickens as the most common source of C. coli infection in France. IMPORTANCE Genome-wide and source attribution studies based on Campylobacter species have shown their importance for the understanding of foodborne infections. Although the use of multilocus sequence typing based on 7 genes from C. jejuni is a powerful method to structure populations, when applied to C. coli , results have not clearly demonstrated its robustness. Therefore, we aim to provide more accurate data based on the identification of single-nucleotide polymorphisms. Results from this study reveal an important number of host-segregating SNPs, found in proteins involved in motility, membrane functions, or DNA repair systems. These findings offer new, interesting opportunities for further study of C. coli adaptation to its environment. Additionally, the results demonstrate that poultry is potentially the main reservoir of C. coli in France., (Copyright © 2020 American Society for Microbiology.)

A growing number of researchers studying horse domestication come to a conclusion that this process happened in multiple locations and involved multiple wild maternal lines. The most promising approach to address this problem involves mitochondrial haplotype comparison of wild and domestic horses from various locations coupled with studies of possible migration routes of the ancient shepherds. Here, we sequenced complete mitochondrial genomes of six horses from burials of the Ukok plateau (Russia, Altai Mountains) dated from 2.7 to 1.4 thousand years before present and a single late Pleistocene wild horse from the neighboring region (Denisova cave). Sequencing data indicates that the wild horse belongs to an extinct pre-domestication lineage. Integration of the domestic horse data with known Eurasian haplotypes of a similar age revealed two distinct groups: the first one widely distributed in Europe and presumably imported to Altai, and the second one specific for Altai Mountains and surrounding area., Competing Interests: The authors have declared that no competing interests exist.

A survey of influenza A viruses (IAVs) in the Mingalar Taung Nyunt live bird market (MTN-LBM), Yangon, Myanmar, was conducted from December 2017 to December 2018. During the survey, 455 swab samples were collected from broilers, layers, backyard chickens and ducks from the MTN-LBM. Ninety-one pooled samples were screened for IAVs by real-time RT-PCR specific to the M gene. Positive pooled samples were individually retested for IAVs. In total, 2.63% of individual samples (12/455) were positive for IAVs. Out of 12 samples, seven samples from layer chickens and the environment were identified as IAV subtype H5N1. In this study, four IAVs were successfully isolated and further characterized by whole genome sequencing. Whole genome sequence analysis revealed that the viruses were characterized as highly pathogenic avian influenza virus subtype H5N1 (HPAIV-H5N1) of clade 2.3.2.1c. Phylogenetic and genetic analyses showed that Myanmar HPAIV-H5N1 was closely related to HPAIV-H5N1 clade 2.3.2.1c isolated from China and Vietnam in 2014. Our results suggested that the live bird market in Myanmar represents a significant risk of HPAIV-H5N1 transmission in poultry and humans. Moreover, HPAIV-H5N1 clade 2.3.2.1c is widely distributed in South-East Asia including Myanmar., (© 2020 Blackwell Verlag GmbH.)

Milking speed (MS) and temperament (MT) are 2 workability traits of great importance in dairy cattle production and breeding. This is mainly due to an increased intensification of the worldwide production systems and greater adoption of precision technologies with less human-cattle interaction. Both MS and MT are heritable traits and thus, genomic selection is a promising tool to expedite their genetic progress. However, the genetic architecture and biological mechanisms underlying the phenotypic expression of these traits remain underexplored. In this study, we investigated the association of >5.7 million imputed whole-genome sequence variants with MT and MS in 4,381 and 4,219 North American Holstein cattle, respectively. The statistical analyses were performed using a mixed linear model fitting a polygenic effect. We detected 40 and 35 significant SNPs independently associated with MT and MS, respectively, which were distributed across 26 chromosomes. Eight candidate genes (GRIN3A, KCNJ3, BOSTAUV1R417, BOSTAUV1R419, MAP2K5, KCTD3, GAP43, and LSAMP) were suggested to play an important role in MT as they are involved in biologically relevant pathways, such as glutamatergic synapse, vomeronasal receptor and oxytocin signaling. Within their coding and upstream sequences, we used an independent data set to further detect or validate significantly differentiated SNP between cattle breeds with known differences in MT. There were fewer candidate genes potentially implicated in MS, but immunity-related genes (e.g., BOLA-NC1 and LOC512672), also identified in other populations, were validated in this study. The significant SNP and novel candidate genes identified contribute to a better understanding of the biological mechanisms underlying both traits in dairy cattle. This information will also be useful for the optimization of prediction of genomic breeding values by giving greater weights to SNP located in the genomic regions identified., (Copyright © 2020 American Dairy Science Association. Published by Elsevier Inc. All rights reserved.)

Coronaviruses (CoVs) is showing obvious interspecies transmission, such as the SARS-CoV, MERS-CoV and SARS-CoV-2. Here, the emerging porcine deltacoronavirus (PDCoV) strain, isolated from Shanghai, China, broadly infects porcine, human and chicken cells in vitro. Previously studies by our group and others have confirmed that PDCoV nucleocapsid (N) protein performs an important role in antagonizing retinoic acid-induced gene I-like receptor (RLR) activation. However, the mechanism of PDCoV N protein suppressing porcine type I IFN production remains unclear, especially the downstream of porcine RLR signaling pathway. In the present study, porcine IRF7 (poIRF7) was identified as the interaction protein of PDCoV N protein through LC-MS/MS. The poIRF7 (268-487aa) was the key region of binding PDCoV N protein. Although IRF7 is a conserved functional protein in species, the PDCoV N protein has been confirmed to interact with only poIRF7 and significantly decrease poIRF7-induced type I IFN production, but not human or chicken IRF7. Furthermore, PDCoV N protein can promote poIRF7 degradation via the ubiquitin-proteasome pathway, which directly increased the K6, K11, and K29-linked polyubiquitination of poIRF7. Lysine 359 of poIRF7 was a key site in PDCoV N protein inducing poIRF7 degradation. Taken together, our results reveal a novel mechanism that PDCoV N protein could species-specifically interact with poIRF7 and then promote its degradation to suppress porcine type I IFN production. The novel findings provide a new insight into PDCoV and other zoonotic coronavirus evading the innate immune response of different species., (Copyright © 2020 Elsevier B.V. All rights reserved.)

As genotypic data are moving from SNP chip toward whole-genome sequence, the accuracy of genomic prediction (GP) exhibits a marginal gain, although all genetic variation, including causal genes, are contained in whole-genome sequence data. Meanwhile, genetic analyses on complex traits, such as genome-wide association studies, have identified an increasing number of genomic regions, including potential causal genes, which would be reliable prior knowledge for GP. Many studies have tried to improve the performance of GP by modifying the prediction model to incorporate prior knowledge. Although several plausible results have been obtained from model modification or strategy optimization, most of them were validated in a specific empirical population with a limited variety of genetic architecture for complex traits. An alternative approach is to use simulated genetic architecture with known causal genes (e.g., simulated causative SNP) to evaluate different GP models with given causal genes. Our objectives were to (1) evaluate the performance of GP under a variety of genetic architectures with a subset of known causal genes and (2) compare different GP models modified by highlighting causal genes and different strategies to weight causal genes. In this study, we simulated pseudo-phenotypes under a variety of genetic architectures based on the real genotypes and phenotypes of a dairy cattle population. Besides classical genomic best linear unbiased prediction, we evaluated 3 modified GP models that highlight causal genes as follows: (1) by treating them as fixed effects, (2) by treating them as a separate random component, and (3) by combining them into the genomic relationship matrix as random effects. Our results showed that highlighting the known causal genes, which explained a considerable proportion of genetic variance in the GP models, increased the predictive accuracy. Combining all given causal genes into the genomic relationship matrix was the optimal strategy under all the scenarios validated, and treating causal genes as a separate random component is also recommended, when more than 20% of genetic variance was explained by known causal genes. Moreover, assigning differential weights to each causal gene further improved the predictive accuracy., (Copyright © 2020 American Dairy Science Association. Published by Elsevier Inc. All rights reserved.)

It has been hypothesized that idiopathic hypertriglyceridemia in Miniature Schnauzers is hereditary, but the gene responsible has yet to be identified. The objective of this study was to determine if there were coding variants in the apolipoprotein C2 (APOC2) gene in Miniature Schnauzers with idiopathic hypertriglyceridemia. Blood samples from 12 Miniature Schnauzers with idiopathic hypertriglyceridemia were analyzed. Genomic DNA was extracted from whole blood, and the three coding exons of APOC2 were amplified by PCR. The PCR amplicons were sequenced and analyzed for variants relative to the canine reference genome (CanFam3.1 assembly). A second objective was to determine the extent of variation in coding exons of APOC2 in a large and diverse canine population using the Dog Biomedical Variant Database Consortium variant catalog, comprised of whole genome sequencing variant calls from 582 dogs of 126 breeds and eight wolves. There were no variants detected in the coding exons of APOC2 for any of the 12 Miniature Schnauzers with idiopathic hypertriglyceridemia. Variants in the coding exons of APOC2 were also rare in the Dog Biomedical Variant Database Consortium variant catalog; a single synonymous variant was identified in a heterozygous state in a Tibetan Mastiff. Thus, we concluded that coding variants in APOC2 are unlikely to be a major cause of idiopathic hypertriglyceridemia in North American Miniature Schnauzers and furthermore, that such coding variants are rare in the canine population., (Copyright © 2020 Elsevier Ltd. All rights reserved.)

Dogs are a potential source of zoonotic Salmonella transmission. We had previously estimated the prevalence of Salmonella shedding among shelter dogs throughout Texas using a repeated cross-sectional study design. Our current objectives were to fully characterize the Salmonella isolates and to assess their relatedness, using whole-genome sequencing. Antimicrobial resistance (AMR) genes were detected in 4/27 (15%) of the isolates. The fosfomycin resistance gene fosA7 was identified in two isolates; to our knowledge, there are no published reports of this gene in canine Salmonella isolates. The biocide resistance gene qacEdelta1, conferring resistance to quaternary ammonium compounds, was detected in an isolate that had four additional AMR genes. The most frequently identified serotypes were Newport (6/27, 22%) and Javiana (4/27, 15%), both of which were widespread among animal shelters. For these serotypes, there was evidence of both transmission of Salmonella within the shelter environment and separate introductions of Salmonella into a shelter. Several canine Salmonella isolates were closely related to human clinical isolates (four canine isolates within 10 SNPs and six more within 20 SNPs), suggesting a shared pathogen population. Educational outreach programmes targeting animal shelter workers would be useful for optimizing knowledge of Salmonella and other canine-associated zoonotic pathogens., (© 2020 The Authors. Veterinary Medicine and Science Published by John Wiley & Sons Ltd.)

Swine influenza virus (SIV) causes respiratory diseases in pigs and has impacts on both animal and human health. In this study, we conducted swine influenza surveillance in pig farms in the Yangon and Bago regions, Myanmar, during 2017-2019. Nasal swabs (n = 500) were collected from pigs in 10 swine farms. Our results showed that 11 out of 100 pooled samples (11%) were positive for influenza A virus (IAV) by real-time RT-PCR. Five SIVs could be isolated and could be subtyped as SIV-H1N1 (n = 4) or SIV-H3N2 (n = 1). The viruses were further characterized by whole-genome sequencing and classified as pdmH1N1-2009 (n = 3), reassortant H1N1 (n = 1) or reassortant H3N2 (n = 1). Phylogenetic analysis of Myanmar SIVs showed that all genes of the three SIV-H1N1 (pdmH1N1-2009) were clustered with viruses of the pdm/09 lineage. For one SIV-H1N1 (rH1N1), the HA1 gene was clustered with those of endemic SIVs of the classical swine lineage, and seven genes were clustered with those of viruses of the pdm/09 lineage. For SIV-H3N2 (rH3N2), the HA3 and NA2 genes were clustered with those of endemic SIVs of the human-like swine lineage, while six internal genes were clustered with those of viruses of the pdm/09 lineage. Genetic analysis indicated that all the Myanmar SIVs possessed amino acids that favour binding to the human receptor. All the Myanmar SIVs contained amino acids related to amantadine resistance but not oseltamivir resistance. Notably, the pdmH1N1-2009 virus might have been circulating in the Myanmar pig population for a period of time after pdmH1N1-2009 outbreaks in humans. Then, reassortment between endemic SIV-H1N1 or SIV-H3N2 and pdmH1N1-2009 in pig farms in Myanmar could have occurred. Our findings ascertained the genetic diversity of SIVs, especially pdmH1N1-2009, in the pig population in Myanmar, with zoonotic and reverse zoonotic potentials., (© 2020 Blackwell Verlag GmbH.)

Sockeye salmon (Oncorhynchus nerka) is a commercially and culturally important species to the people that live along the northern Pacific Ocean coast. There are two main sockeye salmon ecotypes-the ocean-going (anadromous) ecotype and the fresh-water ecotype known as kokanee. The goal of this study was to better understand the population structure of sockeye salmon and identify possible genomic differences among populations and between the two ecotypes. In pursuit of this goal, we generated the first reference sockeye salmon genome assembly and an RNA-seq transcriptome data set to better annotate features of the assembly. Resequenced whole-genomes of 140 sockeye salmon and kokanee were analyzed to understand population structure and identify genomic differences between ecotypes. Three distinct geographic and genetic groups were identified from analyses of the resequencing data. Nucleotide variants in an immunoglobulin heavy chain variable gene cluster on chromosome 26 were found to differentiate the northwestern group from the southern and upper Columbia River groups. Several candidate genes were found to be associated with the kokanee ecotype. Many of these genes were related to ammonia tolerance or vision. Finally, the sex chromosomes of this species were better characterized, and an alternative sex-determination mechanism was identified in a subset of upper Columbia River kokanee., Competing Interests: The authors have declared that no competing interests exist.

Background: Disease prevention and control is a significant part in the ex-situ conservation of the endangered red panda (Ailurus fulgens), being bacterial infection is one of the most important health threats to the captive population. To date, studies about the infection caused by Escherichia coli in the red panda are scarce. This study was conducted to determine the cause of death of a captive red panda through clinical symptoms, complete blood count, biochemical analysis, pathological diagnosis and bacterial whole genome sequencing., Case Presentation: The following report describes a case of a 1.5 year old captive red panda (Ailurus fulgens) that was found lethargic and anorectic. She was moved to the quarantine area for daily treatment with 50 mg of Cefpodoxime Proxetil. During the three-day treatment, she did not eat or defecate, and then died. Clinical hematology revealed the values of neutrophils, alanine aminotransferase (ALT), aspartate aminotransferase (AST) and blood urea nitrogen (BUN) were significantly higher. Histological analysis demonstrated major pathological damage in the kidneys, liver and lungs, characterized by hyperemia, parenchymal cell degeneration and necrosis and inflammatory cell infiltration which were predominantly neutrophilic. A bacterial strain confirmed as Escherichia coli was isolated post mortem. Whole genome sequencing of the E. coli showed the complete genome size was 4.99 Mbp. PapA, PapC, OmpA, OmpU and other virulence factors which specific to Uropathogenic Escherichia coli (UPEC) were found in the isolate. Among the virulence factors, P pili, type I pili and related factors of the iron uptake system were associated with nephrotoxicity., Conclusion: The red panda died of bacterial infection caused by an uropathogenic strain of Escherichia coli. The pathogenic mechanisms of the strain are closely related to the expression of specific virulence genes.

To convert external light into internal neural signal, vertebrates rely on a special group of proteins, the visual opsins. Four of the five types of visual opsins-short-wavelength sensitive 1 (Sws1), short-wavelength sensitive 2 (Sws2), medium-wavelength sensitive (Rh2), and long-wavelength sensitive (Lws)-are expressed in cone cells for scotopic vision, with the fifth, rhodopsin (Rh1), being expressed in rod cells for photopic vision. Fish often display differing ontogenetic cone opsin expression profiles, which may be related to dietary and/or habitat ontogenetic shift. The western mosquitofish (Gambusia affinis) is an aggressive invader that has successfully colonized every continent except Antarctica. The strong invasiveness of this species may be linked to its visual acuity since it can inhabit turbid waters better than other fishes. By genome screening and transcriptome analysis, we identify seven cone opsin genes in the western mosquitofish, including one sws1, two sws2, one rh2, and three lws. The predicted maximal absorbance wavelength (λmax) values of the respective proteins are 353 nm for Sws1, 449 nm for Sws2a, 408 nm for Sws2b, 516 nm for Rh2-1, 571 nm for Lws-1, and 519 nm for Lws-3. Retention of an intron in the lws-r transcript likely renders this visual opsin gene non-functional. Our real-time quantitative PCR demonstrates that adult male and female western mosquitofish do not differ in their cone opsin expression profiles, but we do reveal an ontogenetic shift in cone opsin expression. Compared to adults, larvae express proportionally more sws1 and less lws-1, suggesting that the western mosquitofish is more sensitive to shorter wavelengths in the larval stage, but becomes more sensitive to longer wavelengths in adulthood., Competing Interests: The authors have declared that no competing interests exist.

Cattle pastoralism plays a central role in human livelihood in Africa. However, the genetic history of its success remains unknown. Here, through whole-genome sequence analysis of 172 indigenous African cattle from 16 breeds representative of the main cattle groups, we identify a major taurine × indicine cattle admixture event dated to circa 750-1,050 yr ago, which has shaped the genome of today's cattle in the Horn of Africa. We identify 16 loci linked to African environmental adaptations across crossbred animals showing an excess of taurine or indicine ancestry. These include immune-, heat-tolerance- and reproduction-related genes. Moreover, we identify one highly divergent locus in African taurine cattle, which is putatively linked to trypanotolerance and present in crossbred cattle living in trypanosomosis-infested areas. Our findings indicate that a combination of past taurine and recent indicine admixture-derived genetic resources is at the root of the present success of African pastoralism.

Limousin, a renowned beef breed originating from central France, has been selectively bred over the last 100 years to improve economically important traits. We used whole-genome sequencing data from 10 unrelated Limousin bull calves to detect polymorphisms and identify regions under selection. A total of 13 943 766 variants were identified. Moreover, 311 852 bi-allelic SNPs and 92 229 indels located on autosomes were fixed for the alternative allele in all sequenced animals, including the previously reported missense deleterious F94L mutation in MSTN. We performed a whole-genome screen to discover genomic regions with excess homozygosity, using the pooled heterozygosity score and identified 171 different candidate selective sweeps. In total, 68 candidate genes were found in only 57 of these regions, indicating that a large fraction of the genome under selection might lie in non-coding regions and suggesting that a majority of adaptive mutations might be regulatory in nature. Many QTL were found within candidate selective sweep regions, including QTL associated with shear force or carcass weight. Among the putative selective sweeps, we located genes (MSTN, NCKAP5, RUNX2) that potentially contribute to important phenotypes in Limousin. Several candidate regions and genes under selection were also found in previous genome-wide selection scans performed in Limousin. In addition, we were able to pinpoint candidate causative regulatory polymorphisms in GRIK3 and RUNX2 that might have been under selection. Our results will contribute to improved understanding of the mechanisms and targets of artificial selection and will facilitate the interpretation of GWASs performed in Limousin., (© 2020 Stichting International Foundation for Animal Genetics.)

The complete genomic sequence along with phylogenetic analyses of an adenovirus (AdV), isolated from a dead captive pygmy marmoset (Callithrix pygmaea) from a Hungarian zoo is reported. Earlier, based on the phylogenetic analysis of the sequence of a PCR-amplified fragment from the DNA polymerase gene, the pygmy marmoset AdV (PMAdV) has been reported to cluster closest to certain chiropteran AdVs. In the following years similar AdVs were discovered in additional mammalian hosts, including a skunk (Mephitis mephitis), African pygmy hedgehogs (Atelerix albiventris), North American porcupines (Erethizon dorsatum) and grey fox (Urocyon cinereoargenteus). After the full genome analysis of the skunk adenovirus (SkAdV-1), a novel species Skunk mastadenovirus A (SkAdV-A) has been established. The AdVs, originating from the African pygmy hedgehogs, have been found to belong to virus species SkAdV-A. Partial gene sequences from the porcupine AdVs have also implied their very close genetic relatedness to SkAdV-A. The complete genomic sequence of PMAdV, examined in this study, was found to share 99.83% nucleotide identity with SkAdV-1, thus unequivocally represents a genomic variant of SkAdV-1. The observation that viruses classifiable as SkAdV-A are able to infect and cause diseases in several, distantly related mammals seems to deserve further studies to elucidate the infection biology of this intriguing AdV.

The production of male and female offspring is often determined by the presence of specific sex chromosomes which control sex-specific expression, and sex chromosomes evolve through reduced recombination and specialized gene content. Here we present the genomes of Chrysomya rufifacies, a monogenic blow fly (females produce female or male offspring, exclusively) by separately sequencing and assembling each type of female and the male. The genomes (> 25X coverage) do not appear to have any sex-linked Muller F elements (typical for many Diptera) and exhibit little differentiation between groups supporting the morphological assessments of C. rufifacies homomorphic chromosomes. Males in this species are associated with a unimodal coverage distribution while females exhibit bimodal coverage distributions, suggesting a potential difference in genomic architecture. The presence of the individual-sex draft genomes herein provides new clues regarding the origination and evolution of the diverse sex-determining mechanisms observed within Diptera. Additional genomic analysis of sex chromosomes and sex-determining genes of other blow flies will allow a refined evolutionary understanding of how flies with a typical X/Y heterogametic amphogeny (male and female offspring in similar ratios) sex determination systems evolved into one with a dominant factor that results in single sex progeny in a chromosomally monomorphic system.

Background: Brachygnathia, cardiomegaly and renal hypoplasia syndrome (BCRHS, OMIA 001595-9940 ) is a previously reported recessively inherited disorder in Australian Poll Merino/Merino sheep. Affected lambs are stillborn with various congenital defects as reflected in the name of the disease, as well as short stature, a short and broad cranium, a small thoracic cavity, thin ribs and brachysternum. The BCRHS phenotype shows similarity to certain human short stature syndromes, in particular the human 3M syndrome-2. Here we report the identification of a likely disease-causing variant and propose an ovine model for human 3M syndrome-2., Results: Eight positional candidate genes were identified among the 39 genes in the approximately 1 Mb interval to which the disease was mapped previously. Obscurin like cytoskeletal adaptor 1 (OBSL1) was selected as a strong positional candidate gene based on gene function and the resulting phenotypes observed in humans with mutations in this gene. Whole genome sequencing of an affected lamb (BCRHS3) identified a likely causal variant ENSOARG00000020239:g.220472248delC within OBSL1. Sanger sequencing of seven affected, six obligate carrier, two phenotypically unaffected animals from the original flock and one unrelated control animal validated the variant. A genotyping assay was developed to genotype 583 animals from the original flock, giving an estimated allele frequency of 5%., Conclusions: The identification of a likely disease-causing variant resulting in a frameshift (p.(Val573Trpfs*119)) in the OBSL1 protein has enabled improved breeding management of the implicated flock. The opportunity for an ovine model for human 3M syndrome and ensuing therapeutic research is promising given the availability of carrier ram semen for BCRHS.

Background: Extensive drug resistance (XDR) is an emerging concern with Mannheimia haemolytica, and a variety of testing methods are available for characterizing in vitro antimicrobial susceptibility., Objectives: To compare the concordance among disk diffusion, broth microdilution, and whole genome sequencing (WGS) for susceptibility testing of M. haemolytica before and after mass treatment using tulathromycin., Animals: Forty-eight M. haemolytica isolates collected from high-risk beef stocker calves before and after mass treatment (metaphylaxis) using tulathromycin (Draxxin, Zoetis, Parsippany, NJ) given at the label dosage of 2.5 mg/kg body weight SC in the neck., Methods: In vitro antimicrobial susceptibility was determined for all 48 isolates using disk diffusion, broth microdilution, and WGS. Concordance was calculated between pairs of susceptibility testing methods as follows: number of isolates classified identically by the 2 testing methods for each timepoint, divided by the number of isolates tested at that timepoint. Discordance was calculated as follows: number of isolates classified differently by the 2 testing methods for each timepoint, divided by the number of isolates tested at that timepoint., Results: Concordance between testing methods ranged from 42.3% to 100%, depending on antimicrobial evaluated, timing of sample collection, and testing method used. Very major errors were identified in up to 7.7% of classifications whereas minor errors were seen in up to 50% of classifications depending on antimicrobial evaluated, timing of sample collection, and testing method used., Conclusions and Clinical Importance: Our results show that discrepancies in the results of different susceptibility testing methods occur and suggest a need for greater harmonization of susceptibility testing methods., (© 2020 The Authors. Journal of Veterinary Internal Medicine published by Wiley Periodicals LLC on behalf of American College of Veterinary Internal Medicine.)

Although geese possess an adaptive physiological capacity for lipid storage, few candidate genes contributing to this ability are characterised. By comparing the genomes of individuals with extremely high and low fatty liver weights (FLW), candidate genes were identified, including ARAP2, GABRE, and IL6. Single-nucleotide polymorphisms in or near these genes were significantly (p < 0.05) associated with carcass traits (FLW) and biochemical indexes (very-low-density lipoprotein and N-terminal procollagen III), suggesting contribution to trait variation. A common variant at the 5'-end of LCORL explained ~ 18% and ~ 26% of the phenotypic variance in body weight with/without overfeeding and had significant effects on FLW (p < 0.01). ZFF36L1, ARHGEF1 and IQCJ, involved in bile acid metabolism, blood pressure, and lipid concentration modulation, were also identified. The presence of highly divergent haplotypes within these genes suggested involvement in protection against negative effects from excessive lipids in the liver or circulatory system. Based on this and transcriptomic data, we concluded that geese hepatosteatosis results from severe imbalance between lipid accumulation and secretion, comparable to human non-alcohol fatty liver disease but involving other genes. Our results provided valuable insights into the genesis of geese fatty liver and detected potential target genes for treatment of lipid-related diseases.

Pigs are one of the earliest domesticated animals and multiple breeds have been developed to meet the various demands of consumers. EigenGWAS is a novel strategy to identify candidate genes that underlying population genetic differences and to infer candidate regions under selection as well. In this study, EigenGWAS and F st analyses were performed using the public re-sequencing data of three typical commercial pig breeds, Duroc, Landrace and Yorkshire. The intersection of genome-wide significant SNPs detected by EigenGWAS and top-ranked 1% SNPs of F st results were treated as signals under selection. Using the data of all three breeds, 3062 signals under selection were detected and the nearby genomic regions within 300 kb upstream and downstream covered 6.54% of whole genome. Pairs of breeds were analysed along with the pathway analysis. The gene function enrichment results indicated that many candidate genes located in the genomic regions of the signals under selection were associated with biological processes related to growth, metabolism, reproduction, sensory perception, etc. Among the candidate genes, the FSHB, AHR, PTHLH, KDR and FST genes were reported to be associated with reproductive performance; the KIT, KITLG, MITF, MC1R and EDNRB genes were previously identified to affect coat colour; the RETREG1, TXNIP, BMP5, PPARD and RBP4 genes were reported to be associated with lipid metabolism and growth traits. The identified genetic differences across the three commercial breeds will advance understanding of the artificial selection history of pigs and the signals under selection will suggest potential uses in pig genomic breeding programmes., (© 2020 Stichting International Foundation for Animal Genetics.)

In this study, we identified copy number variants (CNVs) in 19 European autochthonous pig breeds and in two commercial breeds (Italian Large White and Italian Duroc) that represent important genetic resources for this species. The genome of 725 pigs was sequenced using a breed-specific DNA pooling approach (30-35 animals per pool) obtaining an average depth per pool of 42×. This approach maximised CNV discovery as well as the related copy number states characterising, on average, the analysed breeds. By mining more than 17.5 billion reads, we identified a total of 9592 CNVs (~683 CNVs per breed) and 3710 CNV regions (CNVRs; 1.15% of the reference pig genome), with an average of 77 CNVRs per breed that were considered as private. A few CNVRs were analysed in more detail, together with other information derived from sequencing data. For example, the CNVR encompassing the KIT gene was associated with coat colour phenotypes in the analysed breeds, confirming the role of the multiple copies in determining breed-specific coat colours. The CNVR covering the MSRB3 gene was associated with ear size in most breeds. The CNVRs affecting the ELOVL6 and ZNF622 genes were private features observed in the Lithuanian Indigenous Wattle and in the Turopolje pig breeds respectively. Overall, the genome variability unravelled here can explain part of the genetic diversity among breeds and might contribute to explain their origin, history and adaptation to a variety of production systems., (© 2020 Stichting International Foundation for Animal Genetics.)

Enteric pathogens, such as non-typhoidal Salmonella, Campylobacter and Escherichia coli, can reside in the intestinal tract of many animals, including livestock, companion animals, small mammals and reptiles. Often, these animals can appear healthy; nonetheless, humans can become infected after direct or indirect contact, resulting in a substantial illness burden. An estimated 14% of the 3.2 million illnesses that occur in the United States of America (USA) each year from such enteric pathogens are attributable to animal contact. Surveillance for enteric pathogens in the USA includes the compilation and interpretation of both laboratory and epidemiologic data. However, the authors feel that a collaborative, multisectoral and transdisciplinary - or One Health - approach is needed for data collection and analysis, at every level. In addition, they suggest that the future of enteric illness surveillance lies in the development of improved technologies for pathogen detection and characterisation, such as genomic sequencing and metagenomics. In particular, using whole-genome sequencing to compare genetic sequences of enteric pathogens from humans, food, animals and the environment, can help to predict antimicrobial resistance among these pathogens, determine their genetic relatedness and identify outbreaks linked to a common source. In this paper, the authors describe three recent, multi-state human enteric illness outbreaks linked to animal contact in the USA and discuss how integrated disease surveillance was essential to outbreak detection and response. Additional datasharing between public health and animal health laboratories and epidemiologists at the local, national, regional and international level may help to improve surveillance for emerging animal and human health threats and lead to new opportunities for prevention.

Natural selection and domestication have shaped modern horse populations, resulting in a vast range of phenotypically diverse breeds. Horse breeds are classified into three types (pony, light, and draft) generally based on their body type. Understanding the genetic basis of horse type variation and selective pressures related to the evolutionary trend can be particularly important for current selection strategies. Whole-genome sequences were generated for 14 pony and 32 light horses to investigate the genetic signatures of selection of the horse type in pony and light horses. In the overlapping extremes of the fixation index and nucleotide diversity results, we found novel genomic signatures of selective sweeps near key genes previously implicated in body measurements including C4ORF33 , CRB1 , CPN1 , FAM13A , and FGF12 that may influence variation in pony and light horse types. This study contributes to a better understanding of the genetic background of differences between pony and light horse types.

Polycystic kidney disease (PKD) is one of the most common hereditary diseases in cats, with high prevalence in Persian and Persian-related cats. PKD is caused mainly by an inherited autosomal dominant (AD) mutation, and animals may be asymptomatic for years. We screened 16 cats from various breeds exhibiting a renal abnormality by ultrasound examination and genotyped them for the c.10063C>A transversion on exon 29 of the polycystin-1 ( PKD1 ) gene, by PCR-restriction fragment length polymorphism (PCR-RFLP). Among these cats, a Siamese nuclear family of 4 cats with ancestral hereditary renal failure were screened by whole-genome sequencing (WGS) to determine novel variations in genes associated with both AD and autosomal recessive PKD in humans. During the study period, one cat died as a result of renal failure and was forwarded for autopsy. Additionally, we screened 294 cats asymptomatic for renal disease (Angora, Van, Persian, Siamese, Scottish Fold, Exotic Shorthair, British Shorthair, and mixed breeds) to determine the prevalence of the mutation in cats in Turkey. Ten of the symptomatic and 2 of the asymptomatic cats carried the heterozygous C → A transversion, indicating a prevalence of 62.5% and 0.68%, respectively. In the WGS analysis of 4 cats in the Siamese nuclear family, novel variations were determined in the fibrocystin gene ( PKHD1 ), which was not compatible with dominant inheritance of PKD.

Porcine respiratory and reproductive syndrome virus (PRRSV) causes an economically important disease affecting commercial pork production worldwide. NADC34-like PRRSV has had a strong impact on the U.S. and Peruvian pig industries in recent years and also emerged in northeastern China in 2017. However, the endemic status of NADC34-like PRRSV in China is unclear. In this study, we examined 650 tissue samples collected from 16 Provinces in China from 2018 to 2019. Six NADC34-like PRRSV strains were detected in samples from three Provinces, and the complete genomes of four of these strains were sequenced. Phylogenetic analysis showed that these novel PRRSV strains belong to sublineage 1.5 (or NADC34-like PRRSV), forming two groups in China. Sequence alignment suggested that these novel strains share the same 100-aa deletion in the Nsp2 protein that was identified in IA/2014/NADC34 isolated from the United States in 2014. Recombination analysis revealed that five of eight complete genome sequences are derived from recombination between IA/2014/NADC34 and ISU30 or NADC30. The number and distribution of NADC34-like PRRSVs is increasing in China. Importantly, compared with the currently endemic strain NADC30-like PRRSV, NADC34-like PRRSV has the potential to be an endemic strain in China. This study will help us understand the epidemic status of NADC34-like PRRSV in China and provide data for further monitoring this type of PRRSV in China., (© 2020 Blackwell Verlag GmbH.)

Background: Munchkin cats were founded on a naturally occurring mutation segregating into long-legged and short-legged types. Short-legged cats showed disproportionate dwarfism (chondrodysplasia) in which all four legs are short and are referred as standard Munchkin cats. Long-legged animals are referred as non-standard Munchkin cats. A previous study using genome-wide single nucleotide polymorphisms (SNPs) for genome-wide association analysis identified a significantly associated region at 168-184 Mb on feline chromosome (FCA) B1., Results: In this study, we validated the critical region on FCA B1 using a case-control study with 89 cats and 14 FCA B1-SNPs. A structural variant within UGDH (NC&#95;018726.2:g.173294289&#95;173297592delins108, Felis catus 8.0, equivalent to NC&#95;018726.3:g.174882895&#95;174886198delins108, Felis catus 9.0) on FCA B1 was perfectly associated with the phenotype of short-legged standard Munchkin cats., Conclusion: This UGDH structural variant very likely causes the chondrodysplastic (standard) phenotype in Munchkin cats. The lack of homozygous mutant phenotypes and reduced litter sizes in standard Munchkin cats suggest an autosomal recessive lethal trait in the homozygote state. We propose an autosomal dominant mode of inheritance for the chondrodysplastic condition in Munchkin cats.

Three isolates of Ehrlichia ruminantium (Kümm 2, Omatjenne and Riverside), the causative agent of heartwater in domestic ruminants, were isolated in Ixodes scapularis (IDE8) tick cell cultures using the leukocyte fraction of infected sheep blood. All stocks were successfully propagated in IDE8 cells, whereas initiation attempts using endothelial cell cultures were unsuccessful. Therefore, the new technique should be included in any attempt to isolate field strains of E. ruminantium to enhance the probability of getting E. ruminantium isolates which might not be initiated in endothelial cells. Draft genome sequences of all three isolates were generated and compared with published genomes. The data confirmed previous phylogenetic studies that these three isolates are genetically very close to each other, but distinct from previously characterised E. ruminantium isolates. Genome comparisons indicated that the gene content and genomic synteny were highly conserved, with the exception of the membrane protein families. These findings expand our understanding of the genetic diversity of E. ruminantium and confirm the distinct phenotypic and genetic characteristics shared by these three isolates.

Background: African indigenous taurine cattle display unique adaptive traits shaped by husbandry management, regional climate and exposure to endemic pathogens. They are less productive with respect to milk and meat production which has been associated with amongst others, small size, traditional beliefs, husbandry practices, limited feed resources, disease burden and lack of sustained breeding for trait improvement. This resulted in the severe dwindling of their population size rendering them vulnerable to extinction. The Namchi taurine cattle breed is referred to as [Namchi (Doayo)] and shows resistance traits against trypanosome infection and exposure to tick infestation. Nonetheless, the historically later introduced Zebu cattle are the main cattle breeds in Africa today, even though they suffer more from locally prevailing pathogens. By using a whole genome sequencing approach, we sequenced with high depth for the first time the genomes of five cattle breeds from Cameroon in order to provide a valuable genetic resource for future African cattle breeding: the Namchi, an endangered trypano-tolerant taurine breed, the Kapsiki, an indigenous trypano-susceptible taurine breed, and three Zebu (Bos indicus indicus) breeds: Ngaoundere Gudali, White Fulani and Red Fulani., Results: Approximately 167 Gigabases of raw sequencing data were generated for each breed and mapped to the cattle reference genomes ARS-UCD1.2 and UMD3.1.The coverage was 103 to 140-fold when aligning the reads to ARS-UCD1.2 with an average mapping rate of ~ 99%, and 22 to 30-fold when aligning the reads to UMD3.1 with an average mapping rate of ~ 64%. The single nucleotide polymorphisms (SNPs) obtained from analysis using the genome ARS-UCD1.2 were compared with reference genomes of European Bos taurus Holstein, the Asian Bos indicus Brahman, and the African trypanotolerant N'Dama breeds. A total of ~ 100 million (M) SNPs were identified and 7.7 M of those were breed-specific. An approximately 11.1 M constituted of small insertions and deletions. By using only breed-specific non-synonymous variants we identified genes as genetic signatures and associated Gene Ontology (GO) terms that could explain certain cattle-breed specific phenotypes such as increased tolerance against trypanosome parasites in the Namchi breed and heat tolerance in the Kapsiki breed. Phylogenetic analysis grouped, except for Namchi, the Bos taurus breeds Kapsiki, N'Dama and Holstein together while the B. indicus breeds White and Red Fulani, Gudali and Brahman clustered separately. The deviating result for Namchi indicates a hybrid status of the selected animal with a recent introgression of Zebu genes into its genome., Conclusions: The findings provide the first comprehensive set of genome-wide variant data of the most important Cameroonian cattle breeds. The genomic data shall constitute a foundation for breed amelioration whilst exploiting the heritable traits and support conservation efforts for the endangered local cattle breeds.

White-spotting coat colour phenotypes in cattle are either fixed characteristics of specific cattle breeds or occur sporadically owing to germline genetic variation of solid-coloured parents. A Brown Swiss cow showing a piebald pattern resembling colour-sidedness was referred for genetic evaluation. Both parents were normal solid-brown-coloured cattle. The cow was tested negative for the three known DNA variants in KIT, MITF and TWIST2 associated with different depigmentation phenotypes in Brown Swiss cattle. Whole-genome sequencing of the cow was performed and a heterozygous variant affecting the coding sequence of the bovine KIT gene was identified on chromosome 6. The variant is a 40 bp deletion in exon 9, NM&#95;001166484.1:c.1390&#95;1429del, and leads to a frameshift that is predicted to produce a novel 50 amino acid-long C-terminus replacing almost 50% of the wt KIT protein, including the functionally important intracellular tyrosine kinase domain (NP&#95;001159956.1:p.(Asn464AlafsTer50)). Interestingly, among three available offspring, two solid-coloured daughters were genotyped as homozygous wt whereas a single son showing a slightly milder but still obvious depigmentation phenotype inherited a copy of the novel variant allele. The genetic findings provide strong evidence that the identified loss-of-function KIT variant most likely represents a de novo germline mutation that is causative owing to haploinsufficiency., (© 2020 Stichting International Foundation for Animal Genetics.)

Trioza erytreae is the main vector for 'Candidatus Liberibacter africanus', the causative agent of African Citrus Greening disease. The insect is widespread in Africa, and has recently disseminated to Southwestern Europe. This study aimed at generating reference mitogenome sequences for T. erytreae, as a background for future genetic diversity surveys. Complete mitochondrial sequences of three specimens collected in Ethiopia, Uganda and South Africa were recovered using Ion Torrent technology. The mitogenomes of T. erytreae from Uganda and Ethiopia were highly similar, and distinct from that found in South Africa. The phylogeographic structure of T. erytreae was assessed using genetic clustering and pairwise distances, based on a dataset of public COI sequences recorded as T. erytreae. The dataset revealed ten haplotypes with strong phylogeographic structure in Africa and Europe. Three haplotypes found in Kenya on Clausena anisata belonged to pairs separated by distances as high as 11.2%, and were basal to all other sequences. These results indicate that not all sequences identified as T. erytreae belong to the same species, and that some degree of specificity with different plant hosts is likely to exist. This study provides new baseline information on the diversity of T. erytreae, with potential implications for the epidemiology of African Citrus Greening disease.

Duck Tembusu virus (DTMUV), an emerging infectious disease in ducks, belongs to the Flavivirus genus and Flaviviridae family. The transmission of DTUMV involves mosquito vectors; however, the exact role of mosquitoes in the ecology of DTMUV in Thailand remains unclear. This study was conducted to examine DTMUV detection and characterization from mosquitoes in duck farms in central Thailand. Mosquitoes were collected from two duck farms in Sing Buri Province and two duck farms in Ang Thong Province from September 2015 to July 2016 using four CDC-light traps. A total of 30,841 mosquitoes were collected and identified to seven species (Anopheles (An.) barbirostris, An. stephensi, Culex (Cx.) gelidus, Cx. quinquefasciatus, Cx. tritaeniorhynchus, Mansonia (Ma.) annulifera and Ma. uniformis). The most common collected species from each duck farm and each collection time was Cx. tritaeniorhynchus. Mosquitoes were pooled according to species, location, and collection time and then examined for DTMUV by RT-PCR. A total of 273 mosquito pools were examined, with only one pool of Cx. tritaeniorhynchus collected from Sing Buri Province in November 2015 testing positive for DTMUV. Phylogenetic analysis of the polyprotein genes demonstrated that a mosquito-derived Thai DTMUV was grouped into subcluster 2.1 and most closely related to the 2013 Thai DTMUVs. Thus, this study indicated that Cx. tritaeniorhynchus may play a role as a vector in the transmission of DTMUV in Thailand. However, additional studies concerning the vector competence of this mosquito for DTMUV are needed., (© 2020 Blackwell Verlag GmbH.)

Whole-genome next-generation sequencing was used to investigate the local evolution of viral haemorrhagic septicaemia virus, a serious pathogen affecting economically important fish such as rainbow trout and turbot in Europe and olive flounder in Asia. Sequence analysis showed that all isolates were genotype IVa, but could be classified further into four subgroups (K1-K4). In addition, genomic regions encompassing the nucleoprotein, phosphoprotein, matrix protein and non-virion protein genes, as well as the seven non-coding regions, were relatively conserved, whereas glycoprotein and RNA-dependent RNA polymerase genes were variable in the coding region. Taken together, the data demonstrate that whole-genome next-generation sequencing may be useful for future surveillance, prevention and control strategies against viral haemorrhagic septicaemia., (© 2020 John Wiley & Sons Ltd.)

Background: The coupling of appropriate sequencing strategies and imputation methods is critical for assembling large whole-genome sequence datasets from livestock populations for research and breeding. In this paper, we describe and validate the coupling of a sequencing strategy with the imputation method hybrid peeling in real animal breeding settings., Methods: We used data from four pig populations of different size (18,349 to 107,815 individuals) that were widely genotyped at densities between 15,000 and 75,000 markers genome-wide. Around 2% of the individuals in each population were sequenced (most of them at 1× or 2× and 37-92 individuals per population, totalling 284, at 15-30×). We imputed whole-genome sequence data with hybrid peeling. We evaluated the imputation accuracy by removing the sequence data of the 284 individuals with high coverage, using a leave-one-out design. We simulated data that mimicked the sequencing strategy used in the real populations to quantify the factors that affected the individual-wise and variant-wise imputation accuracies using regression trees., Results: Imputation accuracy was high for the majority of individuals in all four populations (median individual-wise dosage correlation: 0.97). Imputation accuracy was lower for individuals in the earliest generations of each population than for the rest, due to the lack of marker array data for themselves and their ancestors. The main factors that determined the individual-wise imputation accuracy were the genotyping status, the availability of marker array data for immediate ancestors, and the degree of connectedness to the rest of the population, but sequencing coverage of the relatives had no effect. The main factors that determined variant-wise imputation accuracy were the minor allele frequency and the number of individuals with sequencing coverage at each variant site. Results were validated with the empirical observations., Conclusions: We demonstrate that the coupling of an appropriate sequencing strategy and hybrid peeling is a powerful strategy for generating whole-genome sequence data with high accuracy in large pedigreed populations where only a small fraction of individuals (2%) had been sequenced, mostly at low coverage. This is a critical step for the successful implementation of whole-genome sequence data for genomic prediction and fine-mapping of causal variants.

After extinction in the wild in the beginning of the 20th century, the European bison has been successfully recovered in 2 distinct genetic lines from only 12 and 7 captive founders. We here aimed at characterizing the levels of realized inbreeding in these 2 restored lines to provide empirical insights into the genomic footprints left by population recovery from a small number of founders. To that end, we genotyped 183 European bison born over the last 40 years with the Illumina BovineHD beadchip that contained 22 602 informative autosomal single-nucleotide polymorphisms after data filtering. We then identified homozygous-by-descent (HBD) segments and classified them into different age-related classes relying on a model-based approach. As expected, we observed that the strong and recent founder effect experienced by the 2 lines resulted in very high levels of recent inbreeding and in the presence of long HBD tracks (up to 120 Mb). These long HBD tracks were associated with ancestors living approximately from 4 to 32 generations in the past, suggesting that inbreeding accumulated over multiple generations after the bottleneck. The contribution to inbreeding of the most recent groups of ancestors was however found to be decreasing in both lines. In addition, comparison of Lowland individuals born at different time periods showed that the levels of inbreeding tended to stabilize, HBD segments being shorter in animals born more recently which indicates efficient control of inbreeding. Monitoring HBD segment lengths over generations may thus be viewed as a valuable genomic diagnostic tool for populations in conservation or recovery programs., (© The American Genetic Association 2020. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

To explore the mechanisms leading to the poor quality of IVF blastocysts, the single-cell whole-genome methylation sequencing technique was used in this study to analyse the methylation patterns of bovine blastocysts derived from invivo, fresh (IVF) or vitrified (V&#95;IVF) oocytes. Genome methylation levels of blastocysts in the IVF and V&#95;IVF groups were significantly lower than those of the invivo group (P<0.05). In all, 1149 differentially methylated regions (DMRs) were identified between the IVF and invivo groups, 1578 DMRs were identified between the V&#95;IVF and invivo groups and 151 DMRs were identified between the V&#95;IVF and IVF groups. For imprinted genes, methylation levels of insulin-like growth factor 2 receptor (IGF2R) and protein phosphatase 1 regulatory subunit 9A (PPP1R9A) were lower in the IVF and V&#95;IVF groups than in the invivo group, and the methylation level of paternally expressed 3 (PEG3) was lower in the V&#95;IVF group than in the IVF and invivo groups. Genes with DMRs between the IVF and invivo and the V&#95;IVF and IVF groups were primarily enriched in oocyte maturation pathways, whereas DMRs between the V&#95;IVF and invivo groups were enriched in fertilisation and vitrification-vulnerable pathways. The results of this study indicate that differences in the methylation of critical DMRs may contribute to the differences in quality between invitro- and invivo-derived embryos.

Background: The yak is a species of livestock which is crucial for local communities of the Qinghai-Tibet Plateau and adjacent regions and naturally owns one more thoracic vertebra than cattle. Recently, a sub-population of yak termed as the Jinchuan yak has been identified with over half its members own a thoracolumbar vertebral formula of T15L5 instead of the natural T14L5 arrangement. The novel T15L5 positioning is a preferred genetic trait leading to enhanced meat and milk production. Selective breeding of this trait would have great agricultural value and exploration of the molecular mechanisms underlying this trait would both accelerate this process and provide us insight into the development and regulation of somitogenesis., Results: Here we investigated the genetic background of the Jinchuan yak through resequencing fifteen individuals, comprising five T15L5 individuals and ten T14L5 individuals with an average sequencing depth of > 10X, whose thoracolumbar vertebral formulae were confirmed by anatomical observation. Principal component analysis, linkage disequilibrium analysis, phylogenetic analysis, and selective sweep analysis were carried out to explore Jinchuan yak's genetic background. Three hundred and thirty candidate markers were identified as associated with the additional thoracic vertebrae and target sequencing was used to validate seven carefully selected markers in an additional 51 Jinchuan yaks. The accuracies of predicting 15 thoracic vertebrae and 20 thoracolumbar vertebrae with these 7 markers were 100.00 and 33.33% despite they both could only represent 20% of all possible genetic diversity. Two genes, PPP2R2B and TBLR1, were found to harbour the most candidate markers associated with the trait and likely contribute to the unique somitic number and identity according to their reported roles in the mechanism of somitogenesis., Conclusions: Our findings provide a clear depiction of the Jinchuan yak's genetic background and a solid foundation for marker-assistant selection. Further exploitation of this unique population and trait could be promoted with the aid of our genomic resource.

Background: Advances in genome technology have simplified a new comprehension of the genetic and historical processes crucial to rapid phenotypic evolution under domestication. To get new insight into the genetic basis of the dog domestication process, we conducted whole-genome sequence analysis of three wolves and three dogs from Iran which covers the eastern part of the Fertile Crescent located in Southwest Asia where the independent domestication of most of the plants and animals has been documented and also high haplotype sharing between wolves and dog breeds has been reported., Results: Higher diversity was found within the wolf genome compared with the dog genome. A total number of 12.45 million SNPs were detected in all individuals (10.45 and 7.82 million SNPs were identified for all the studied wolves and dogs, respectively) and a total number of 3.49 million small Indels were detected in all individuals (3.11 and 2.24 million small Indels were identified for all the studied wolves and dogs, respectively). A total of 10,571 copy number variation regions (CNVRs) were detected across the 6 individual genomes, covering 154.65 Mb, or 6.41%, of the reference genome (canFam3.1). Further analysis showed that the distribution of deleterious variants in the dog genome is higher than the wolf genome. Also, genomic annotation results from intron and intergenic regions showed that the proportion of variations in the wolf genome is higher than that in the dog genome, while the proportion of the coding sequences and 3'-UTR in the dog genome is higher than that in the wolf genome. The genes related to the olfactory and immune systems were enriched in the set of the structural variants (SVs) identified in this work., Conclusions: Our results showed more deleterious mutations and coding sequence variants in the domestic dog genome than those in wolf genome. By providing the first Iranian dog and wolf variome map, our findings contribute to understanding the genetic architecture of the dog domestication.

Background: The AIG (avrRpt2-induced gene) family of GTPases, characterized by the presence of a distinctive AIG1 domain, is mysterious in having a peculiar phylogenetic distribution, a predilection for undergoing expansion and loss, and an uncertain functional role, especially in invertebrates. AIGs are frequently represented as GIMAPs (GTPase of the immunity associated protein family), characterized by presence of the AIG1 domain along with coiled-coil domains. Here we provide an overview of the remarkably expanded AIG repertoire of the freshwater gastropod Biomphalaria glabrata, compare it with AIGs in other organisms, and detail patterns of expression in B. glabrata susceptible or resistant to infection with Schistosoma mansoni, responsible for the neglected tropical disease of intestinal schistosomiasis., Results: We define the 7 conserved motifs that comprise the AIG1 domain in B. glabrata and detail its association with at least 7 other domains, indicative of functional versatility of B. glabrata AIGs. AIG genes were usually found in tandem arrays in the B. glabrata genome, suggestive of an origin by segmental gene duplication. We found 91 genes with complete AIG1 domains, including 64 GIMAPs and 27 AIG genes without coiled-coils, more than known for any other organism except Danio (with > 100). We defined expression patterns of AIG genes in 12 different B. glabrata organs and characterized whole-body AIG responses to microbial PAMPs, and of schistosome-resistant or -susceptible strains of B. glabrata to S. mansoni exposure. Biomphalaria glabrata AIG genes clustered with expansions of AIG genes from other heterobranch gastropods yet showed unique lineage-specific subclusters. Other gastropods and bivalves had separate but also diverse expansions of AIG genes, whereas cephalopods seem to lack AIG genes., Conclusions: The AIG genes of B. glabrata exhibit expansion in both numbers and potential functions, differ markedly in expression between strains varying in susceptibility to schistosomes, and are responsive to immune challenge. These features provide strong impetus to further explore the functional role of AIG genes in the defense responses of B. glabrata, including to suppress or support the development of medically relevant S. mansoni parasites.

Ascaris lumbricoides and Ascaris suum are parasitic nematodes that mainly parasitize the small intestines of people and pigs, respectively. Ascariasis seriously endangers human health and causes huge economic losses in the pig industry. A. lumbricoides and A. suum have similar morphologies and genetic structures, and occasionally these organisms cross-infect the alternate host. Therefore, their taxonomies are controversial. In this study, the whole genomes of A. lumbricoides (n = 6) and A. suum (n = 6) were resequenced using a HiSeq X Ten sequencing platform. Phylogenetic, principal component, and population structure analyses showed clear genetic differentiation between the two Ascaris populations. Linkage disequilibrium analysis indicated that the A. lumbricoides population was more primitive than the A. suum population. In the selective elimination analysis, 160 and 139 candidate regions were screened in A. lumbricoides and A. suum, respectively, and the selected regions were analyzed by Gene Ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses. The A. lumbricoides population had no significant enrichment in GO terms, but two KEGG pathways, the RNA degradation and tyrosine metabolism pathways, were significantly enriched. Five GO entries and one KEGG pathway, the alanine, aspartate, and glutamate metabolism signaling pathway, were significantly enriched in the A. suum population. An analysis of the demographic histories of Ascaris populations revealed that A. lumbricoides and A. suum had similar trends in effective population size in different historical periods. Ascaris populations peaked about 1 million years ago and then began to decline. In the last glacial period, they dropped to a historical low and continued at this level until the last glacial maximum. This phenomenon may be associated with the cold climate at that time. This study provides new information on the genetic differentiation, evolutionary relationships, gene functional enrichment, and population dynamics of Ascaris populations, with implications for host differences, evolution, and classification of A. lumbricoides and A. suum., (Copyright © 2020 Elsevier B.V. All rights reserved.)

Cervidpoxvirus is one of the more recently designated genera within the subfamily Chordopoxvirinae , with Deerpox virus (DPV) as the only recognized species to date. In this study, the authors describe spontaneous disease and infection in the North American moose ( Alces americanus ) by a novel Cervidpoxvirus , here named Moosepox virus (MPV). Three 4-month-old moose calves developed a multifocal subacute-to-chronic, necrotizing, suppurative-to-granulomatous dermatitis that affected the face and the extremities. Ultrastructurally, all stages of MPV morphogenesis-that is, crescents, spherical immature particles, mature particles, and enveloped mature virus-were observed in skin tissue. In vitro infection with MPV confirmed that its morphogenesis was similar to that of the prototype vaccinia virus. The entire coding region, including 170 putative genes of this MPV, was sequenced and annotated. The sequence length was 164,258 bp with 98.5% nucleotide identity with DPV (strain W-1170-84) based on the whole genome. The genome of the study virus was distinct from that of the reference strain (W-1170-84) in certain genes, including the CD30-like protein (83.9% nucleotide, 81.6% amino acid), the endothelin precursor (73.2% nucleotide including some indels, 51.4% amino acid), and major histocompatibility class (MHC) class I-like protein (81.0% nucleotide, 68.2% amino acid). This study provides biological characterization of a new Cervidpoxvirus attained through in vivo and in vitro ultrastructural analyses. It also demonstrates the importance of whole-genome sequencing in the molecular characterization of poxviruses identified in taxonomically related hosts.

Objectives: A multidrug-resistant Morganella morganii strain (CQ-M7), isolated from the kidney of a diseased Chinese giant salamander in China, was examined with whole genome sequencing to better understand drug tolerance and its pathogenicity., Methods: The draft genome of the investigated strain was assembled using HGA assembler and annotated using Rapid Annotations Subsystems Technology (RAST) server. The contigs were annotated by the appropriate bioinformatics tools available on the National Center for Biotechnology Information (NCBI) website. Antibiotic resistance genes were detected by PCR. Pathogenicity of the isolate was performed on 30 healthy Chinese giant salamanders with different infection dosages., Results: The CQ-M7 strain showed resistance to multiple antimicrobials, especially to aminoglycoside and β-lactam antibiotics. Seventeen drug-resistance genes were detected, which were related to β-lactams, aminoglycosides, fluoroquinolones, tetracyclines, peptide antibiotic, and fosfomycin resistance. Sequence analysis showed the assembled genome size to be 4 966 326bp with 51.16% of GC content, containing 4587 protein-coding genes, 71 pseudogenes, five rRNAs, 80 tRNAs, and five noncoding RNAs. The genome sequence was deposited in GenBank under accession number RQIJ00000000. Artificial infection results indicated that the CQ-M7 strain was a low-virulence strain for the Chinese giant salamander., Conclusion: It is believed that this is the first draft genome of Chinese giant salamander original Morganella morganii strain harbouring multiple antibiotic resistance genes in China. The reported genome sequence could provide insights into antibiotic resistance mechanisms and control strategies of Morganella morganii., (Copyright © 2019 International Society for Antimicrobial Chemotherapy. Published by Elsevier Ltd. All rights reserved.)