GARP and EARP are required for efficient BoHV-1 replication as identified by a genome wide CRISPR knockout screen.

The advances in gene editing bring unprecedented opportunities in high throughput functional genomics to animal research. Here we describe a genome wide CRISPR knockout library, btCRISPRko.v1, targeting all protein coding genes in the cattle genome. Using it, we conducted genome wide screens during Bovine Herpes Virus type 1 (BoHV-1) replication and compiled a list of pro-viral and anti-viral candidates. These candidates might influence multiple aspects of BoHV-1 biology such as viral entry, genome replication and transcription, viral protein trafficking and virion maturation in the cytoplasm. Some of the most intriguing examples are VPS51, VPS52 and VPS53 that code for subunits of two membrane tethering complexes, the endosome-associated recycling protein (EARP) complex and the Golgi-associated retrograde protein (GARP) complex. These complexes mediate endosomal recycling and retrograde trafficking to the trans Golgi Network (TGN). Simultaneous loss of both complexes in MDBKs resulted in greatly reduced production of infectious BoHV-1 virions. We also found that viruses released by these deficient cells severely lack VP8, the most abundant tegument protein of BoHV-1 that are crucial for its virulence. In combination with previous reports, our data suggest vital roles GARP and EARP play during viral protein packaging and capsid re-envelopment in the cytoplasm. It also contributes to evidence that both the TGN and the recycling endosomes are recruited in this process, mediated by these complexes. The btCRISPRko.v1 library generated here has been controlled for quality and shown to be effective in host gene discovery. We hope it will facilitate efforts in the study of other pathogens and various aspects of cell biology in cattle.

Changes in hospital mortality in patients with cancer during the COVID-19 pandemic (ISARIC-CCP-UK): a prospective, multicentre cohort study.

BACKGROUND: Patients with cancer are at greater risk of dying from COVID-19 than many other patient groups. However, how this risk evolved during the pandemic remains unclear. We aimed to determine, on the basis of the UK national pandemic protocol, how factors influencing hospital mortality from COVID-19 could differentially affect patients undergoing cancer treatment. We also examined changes in hospital mortality and escalation of care in patients on cancer treatment during the first 2 years of the COVID-19 pandemic in the UK. METHODS: We conducted a prospective cohort study of patients aged older than 19 years and admitted to 306 health-care facilities in the UK with confirmed SARS-CoV-2 infection, who were enrolled in the International Severe Acute Respiratory and emerging Infections Consortium (ISARIC) WHO Clinical Characterisation Protocol (CCP) across the UK from April 23, 2020, to Feb 28, 2022; this analysis included all patients in the complete dataset when the study closed. The primary outcome was 30-day in-hospital mortality, comparing patients on cancer treatment and those without cancer. The study was approved by the South Central-Oxford C Research Ethics Committee in England (Ref: 13/SC/0149) and the Scotland A Research Ethics Committee (Ref 20/SS/0028), and is registered on the ISRCTN Registry (ISRCTN66726260). FINDINGS: 177 871 eligible adult patients either with no history of cancer (n=171 303) or on cancer treatment (n=6568) were enrolled; 93 205 (52·4%) were male, 84 418 (47·5%) were female, and in 248 (13·9%) sex or gender details were not specified or data were missing. Patients were followed up for a median of 13 (IQR 6-21) days. Of the 6568 patients receiving cancer treatment, 2080 (31·7%) died at 30 days, compared with 30 901 (18·0%) of 171 303 patients without cancer. Patients aged younger than 50 years on cancer treatment had the highest age-adjusted relative risk (hazard ratio [HR] 5·2 [95% CI 4·0-6·6], p<0·0001; vs 50-69 years 2·4 [2·2-2·6], p<0·0001; 70-79 years 1·8 [1·6-2·0], p<0·0001; and >80 years 1·5 [1·3-1·6], p<0·0001) but a lower absolute risk (51 [6·7%] of 763 patients <50 years died compared with 459 [30·2%] of 1522 patients aged >80 years). In-hospital mortality decreased for all patients during the pandemic but was higher for patients on cancer treatment than for those without cancer throughout the study period. INTERPRETATION: People with cancer have a higher risk of mortality from COVID-19 than those without cancer. Patients younger than 50 years with cancer treatment have the highest relative risk of death. Continued action is needed to mitigate the poor outcomes in patients with cancer, such as through optimising vaccination, long-acting passive immunisation, and early access to therapeutics. These findings underscore the importance of the ISARIC-WHO pandemic preparedness initiative. FUNDING: National Institute for Health Research and the Medical Research Council.

Outcome of COVID-19 in hospitalised immunocompromised patients: An analysis of the WHO ISARIC CCP-UK prospective cohort study.

BACKGROUND: Immunocompromised patients may be at higher risk of mortality if hospitalised with Coronavirus Disease 2019 (COVID-19) compared with immunocompetent patients. However, previous studies have been contradictory. We aimed to determine whether immunocompromised patients were at greater risk of in-hospital death and how this risk changed over the pandemic. METHODS AND FINDINGS: We included patients > = 19 years with symptomatic community-acquired COVID-19 recruited to the ISARIC WHO Clinical Characterisation Protocol UK prospective cohort study. We defined immunocompromise as immunosuppressant medication preadmission, cancer treatment, organ transplant, HIV, or congenital immunodeficiency. We used logistic regression to compare the risk of death in both groups, adjusting for age, sex, deprivation, ethnicity, vaccination, and comorbidities. We used Bayesian logistic regression to explore mortality over time. Between 17 January 2020 and 28 February 2022, we recruited 156,552 eligible patients, of whom 21,954 (14%) were immunocompromised. In total, 29% (n = 6,499) of immunocompromised and 21% (n = 28,608) of immunocompetent patients died in hospital. The odds of in-hospital mortality were elevated for immunocompromised patients (adjusted OR 1.44, 95% CI [1.39, 1.50], p < 0.001). Not all immunocompromising conditions had the same risk, for example, patients on active cancer treatment were less likely to have their care escalated to intensive care (adjusted OR 0.77, 95% CI [0.7, 0.85], p < 0.001) or ventilation (adjusted OR 0.65, 95% CI [0.56, 0.76], p < 0.001). However, cancer patients were more likely to die (adjusted OR 2.0, 95% CI [1.87, 2.15], p < 0.001). Analyses were adjusted for age, sex, socioeconomic deprivation, comorbidities, and vaccination status. As the pandemic progressed, in-hospital mortality reduced more slowly for immunocompromised patients than for immunocompetent patients. This was particularly evident with increasing age: the probability of the reduction in hospital mortality being less for immunocompromised patients aged 50 to 69 years was 88% for men and 83% for women, and for those >80 years was 99% for men and 98% for women. The study is limited by a lack of detailed drug data prior to admission, including steroid doses, meaning that we may have incorrectly categorised some immunocompromised patients as immunocompetent. CONCLUSIONS: Immunocompromised patients remain at elevated risk of death from COVID-19. Targeted measures such as additional vaccine doses, monoclonal antibodies, and nonpharmaceutical preventive interventions should be continually encouraged for this patient group. TRIAL REGISTRATION: ISRCTN 66726260.

Systematic comparison of ranking aggregation methods for gene lists in experimental results.

MOTIVATION: A common experimental output in biomedical science is a list of genes implicated in a given biological process or disease. The gene lists resulting from a group of studies answering the same, or similar, questions can be combined by ranking aggregation methods to find a consensus or a more reliable answer. Evaluating a ranking aggregation method on a specific type of data before using it is required to support the reliability since the property of a dataset can influence the performance of an algorithm. Such evaluation on gene lists is usually based on a simulated database because of the lack of a known truth for real data. However, simulated datasets tend to be too small compared to experimental data and neglect key features, including heterogeneity of quality, relevance and the inclusion of unranked lists. RESULTS: In this study, a group of existing methods and their variations that are suitable for meta-analysis of gene lists are compared using simulated and real data. Simulated data were used to explore the performance of the aggregation methods as a function of emulating the common scenarios of real genomic data, with various heterogeneity of quality, noise level and a mix of unranked and ranked data using 20 000 possible entities. In addition to the evaluation with simulated data, a comparison using real genomic data on the SARS-CoV-2 virus, cancer (non-small cell lung cancer) and bacteria (macrophage apoptosis) was performed. We summarize the results of our evaluation in a simple flowchart to select a ranking aggregation method, and in an automated implementation using the meta-analysis by information content algorithm to infer heterogeneity of data quality across input datasets. AVAILABILITY AND IMPLEMENTATION: The code for simulated data generation and running edited version of algorithms: https://github.com/baillielab/comparison_of_RA_methods. Code to perform an optimal selection of methods based on the results of this review, using the MAIC algorithm to infer the characteristics of an input dataset, can be downloaded here: https://github.com/baillielab/maic. An online service for running MAIC: https://baillielab.net/maic. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Biomechanical Testing and Histologic Examination of Intradermal Skin Closure in Dogs Using Barbed Suture Device and Non-Barbed Monofilament Suture.

OBJECTIVE: To compare the biomechanical strength and histologic features of 3-0 Glycomer™ 631 barbed suture (V-LOC™ 90 Absorbable Wound Closure Device, Covidien, Mansfield, MA) to non-barbed 3-0 Glycomer™ 631 suture (Biosyn™, Covidien) for intradermal skin wound closure in the dog. STUDY DESIGN: Randomized, factorial, in vivo. ANIMALS: Eighteen purpose-bred, mature male, and female hound dogs. METHODS: Eighteen adult hound dogs were randomly assigned to 1 of 3 groups designated by postoperative day of assessment. Six skin incisions were made along the dorsum in the thoracolumbar region of each dog with an equal number (n=3) randomly assigned to closure with barbed or non-barbed suture. Six dogs were euthanatized on postoperative days 3, 10, and 14, respectively. Two additional incisions were made on each dog after euthanasia for baseline data (Day 0). The skin incision specimens were harvested for biomechanical testing and histologic evaluation. RESULTS: Non-barbed closure had significantly higher maximum load at failure (P<.001) and stiffness (P<.001) than barbed closure regardless of day. The average tissue reaction score was significantly higher for barbed closure (P=.008), regardless of day. Suturing time for barbed closures was significantly shorter. There was no significant difference in frequency of complications between closures. CONCLUSION: Barbed Glycomer™ 631 closures had a significantly lower maximum load at failure and stiffness, and higher average tissue reaction scores, but showed no difference in short term outcome for intradermal closure of dorsally located skin incisions in dogs.

Rapid identification of bovine MHCI haplotypes in genetically divergent cattle populations using next-generation sequencing.

The major histocompatibility complex (MHC) region contains many genes that are key regulators of both innate and adaptive immunity including the polymorphic MHCI and MHCII genes. Consequently, the characterisation of the repertoire of MHC genes is critical to understanding the variation that determines the nature of immune responses. Our current knowledge of the bovine MHCI repertoire is limited with only the Holstein-Friesian breed having been studied in any depth. Traditional methods of MHCI genotyping are of low resolution and laborious and this has been a major impediment to a more comprehensive analysis of the MHCI repertoire of other cattle breeds. Next-generation sequencing (NGS) technologies have been used to enable high throughput and much higher resolution MHCI typing in a number of species. In this study we have developed a MiSeq platform approach and requisite bioinformatics pipeline to facilitate typing of bovine MHCI repertoires. The method was validated initially on a cohort of Holstein-Friesian animals and then demonstrated to enable characterisation of MHCI repertoires in African cattle breeds, for which there was limited or no available data. During the course of these studies we identified >140 novel classical MHCI genes and defined 62 novel MHCI haplotypes, dramatically expanding the known bovine MHCI repertoire.

Genetic loci inherited from hens lacking maternal behaviour both inhibit and paradoxically promote this behaviour.

BACKGROUND: A major step towards the success of chickens as a domesticated species was the separation between maternal care and reproduction. Artificial incubation replaced the natural maternal behaviour of incubation and, thus, in certain breeds, it became possible to breed chickens with persistent egg production and no incubation behaviour; a typical example is the White Leghorn strain. Conversely, some strains, such as the Silkie breed, are prized for their maternal behaviour and their willingness to incubate eggs. This is often colloquially known as broodiness. RESULTS: Using an F2 linkage mapping approach and a cross between White Leghorn and Silkie chicken breeds, we have mapped, for the first time, genetic loci that affect maternal behaviour on chromosomes 1, 5, 8, 13, 18 and 19 and linkage group E22C19W28. Paradoxically, heterozygous and White Leghorn homozygous genotypes were associated with an increased incidence of incubation behaviour, which exceeded that of the Silkie homozygotes for most loci. In such cases, it is likely that the loci involved are associated with increased egg production. Increased egg production increases the probability of incubation behaviour occurring because egg laying must precede incubation. For the loci on chromosomes 8 and 1, alleles from the Silkie breed promote incubation behaviour and influence maternal behaviour (these explain 12 and 26% of the phenotypic difference between the two founder breeds, respectively). CONCLUSIONS: The over-dominant locus on chromosome 5 coincides with the strongest selective sweep reported in chickens and together with the loci on chromosomes 1 and 8, they include genes of the thyrotrophic axis. This suggests that thyroid hormones may play a critical role in the loss of incubation behaviour and the improved egg laying behaviour of the White Leghorn breed. Our findings support the view that loss of maternal incubation behaviour in the White Leghorn breed is the result of selection for fertility and egg laying persistency and against maternal incubation behaviour.

Cementless total hip replacement in 20 juveniles using BFX™ arthroplasty.

OBJECTIVE: To describe outcome after a minimum of 1 year for total hip replacement (THR) using BioMedtrix BFX™ biologic fixation implants in skeletally immature dogs (6-10 months of age). STUDY DESIGN: Case series. ANIMALS: Dogs (n = 20). METHODS: Medical records (November 2007-June 2010) of 20 dogs, 6-10 months old that had cementless THR were reviewed. Preoperative, immediate, 6-week and >1-year postoperative radiographs were compared. Clinical examination was performed at 6 weeks and >1 year postoperatively. Owner questionnaire was obtained at final follow-up. RESULTS: Cementless THR using BFX™ implants was performed because of debilitating coxofemoral pain, resulting from canine hip dysplasia, after unsatisfactory outcome with medical management. Radiographs taken immediately, 6 weeks and at a mean of 29.8 months (range, 12-48 months) postoperatively revealed satisfactory implant positioning and stability. Significant change in measured cranial-caudal femoral stem fill and level (P < .001 and P = .006, respectively) were recorded at 6 weeks postoperatively without clinical significance. No further change in stem positioning occurred. Acetabular cup orientation remained unchanged throughout follow-up. Complications requiring further surgical intervention were not encountered. One minor superficial infection was recorded and treated. Lameness and pain on manipulation of the affected hip had greatly improved in all dogs by 6 weeks postoperatively. CONCLUSION: Cementless THR can safely be performed in skeletally immature dogs, providing satisfactory return to normal joint function and implant survivability for at least a mean of 29.8 months.

Validation of Candidate Host Cell Entry Factors for Bovine Herpes Virus Type-1 Based on a Genome-Wide CRISPR Knockout Screen.

To identify host factors that affect Bovine Herpes Virus Type 1 (BoHV-1) infection we previously applied a genome wide CRISPR knockout screen targeting all bovine protein coding genes. By doing so we compiled a list of both pro-viral and anti-viral proteins involved in BoHV-1 replication. Here we provide further analysis of those that are potentially involved in viral entry into the host cell. We first generated single cell knockout clones deficient in some of the candidate genes for validation. We provide evidence that Polio Virus Receptor-related protein (PVRL2) serves as a receptor for BoHV-1, mediating more efficient entry than the previously identified Polio Virus Receptor (PVR). By knocking out two enzymes that catalyze HSPG chain elongation, HST2ST1 and GLCE, we further demonstrate the significance of HSPG in BoHV-1 entry. Another intriguing cluster of candidate genes, COG1, COG2 and COG4-7 encode six subunits of the Conserved Oligomeric Golgi (COG) complex. MDBK cells lacking COG6 produced fewer but bigger plaques compared to control cells, suggesting more efficient release of newly produced virions from these COG6 knockout cells, due to impaired HSPG biosynthesis. We further observed that viruses produced by the COG6 knockout cells consist of protein(s) with reduced N-glycosylation, potentially explaining their lower infectivity. To facilitate candidate validation, we also detailed a one-step multiplex CRISPR interference (CRISPRi) system, an orthogonal method to KO that enables quick and simultaneous deployment of three CRISPRs for efficient gene inactivation. Using CRISPR3i, we verified eight candidates that have been implicated in the synthesis of surface heparan sulfate proteoglycans (HSPGs). In summary, our experiments confirmed the two receptors PVR and PVRL2 for BoHV-1 entry into the host cell and other factors that affect this process, likely through the direct or indirect roles they play during HSPG synthesis and glycosylation of viral proteins.

Sequencing and Analysis of Lumpy Skin Disease Virus Whole Genomes Reveals a New Viral Subgroup in West and Central Africa.

Lumpy skin disease virus (LSDV) is a member of the capripoxvirus (CPPV) genus of the Poxviridae family. LSDV is a rapidly emerging, high-consequence pathogen of cattle, recently spreading from Africa and the Middle East into Europe and Asia. We have sequenced the whole genome of historical LSDV isolates from the Pirbright Institute virus archive, and field isolates from recent disease outbreaks in Sri Lanka, Mongolia, Nigeria and Ethiopia. These genome sequences were compared to published genomes and classified into different subgroups. Two subgroups contained vaccine or vaccine-like samples ("Neethling-like" clade 1.1 and "Kenya-like" subgroup, clade 1.2.2). One subgroup was associated with outbreaks of LSD in the Middle East/Europe (clade 1.2.1) and a previously unreported subgroup originated from cases of LSD in west and central Africa (clade 1.2.3). Isolates were also identified that contained a mix of genes from both wildtype and vaccine samples (vaccine-like recombinants, grouped in clade 2). Whole genome sequencing and analysis of LSDV strains isolated from different regions of Africa, Europe and Asia have provided new knowledge of the drivers of LSDV emergence, and will inform future disease control strategies.

ArkMAP: integrating genomic maps across species and data sources.

BACKGROUND: The visualisation of genetic and genomic maps aligned within and between species and across data sources can be used to inform studies of genome evolution, assist genome assembly projects and aid gene discovery and identification. Whilst annotation, integration and exploration of assembled genome sequences is well supported, there are fewer tools available which can display genetic maps for less well-characterized species, and integrate these maps with annotated reference genomes to support cross species comparisons. RESULTS: We have developed a desktop application to draw and align genetic and genomic maps, retrieved from remote data sources or loaded as local files. Maps can be retrieved from our public map database ArkDB or from any Ensembl data source (i.e. Ensembl and Ensembl Genomes). By using the JEnsembl API, maps can be drawn for any release version of any of the thousands of species present in Ensembl data sources, allowing not only inter-specific comparisons, but also comparisons between different versions/revisions of assembled genomes. Maps can be aligned by relating identical or synonymous markers across maps, or through the gene homology/orthology relationship data stored in the Ensembl Compara databases, allowing ready visualization of regions of conserved synteny between species. The map drawing canvas is highly configurable, supports interactive exploration of maps, markers and relationships and allows export of publication quality graphics. CONCLUSIONS: ArkMAP allows users to draw and interactively explore gene and variation maps for any version of any annotated genome curated in the Ensembl data sources, and to integrate local mapping data. The maps and inter-map relationships drawn are highly configurable and ArkMAP may be used to produce publication quality graphics. ArkMAP is freely available as an auto-updating Java 'Web Start' application, or as a standalone archived application.

Decreased expression of the satiety signal receptor CCKAR is responsible for increased growth and body weight during the domestication of chickens.

Animal domestication has resulted in changes in growth and size. It has been suggested that this may have involved selection for differences in appetite. Divergent growth between chickens selected for egg laying or meat production is one such example. The neurons expressing AGRP and POMC in the basal hypothalamus are important components of appetite regulation, as are the satiety feedback pathways that carry information from the intestine, including CCK and its receptor CCKAR (CCK1 receptor). Using 16 generations of a cross between a fast and a relatively slow growing strain of chicken has identified a region on chromosome 4 downstream of the CCKAR gene, which is responsible for up to a 19% difference in body weight at 12 wk of age. Animals possessing the high-growth haplotype at the locus have lower expression of mRNA and immunoreactive CCKAR in the brain, intestine, and exocrine organs, which is correlated with increased levels of orexigenic AGRP in the hypothalamus. Animals with the high-growth haplotype are resistant to the anorectic effect of exogenously administered CCK, suggesting that their satiety set point has been altered. Comparison with traditional breeds shows that the high-growth haplotype has been present in the founders of modern meat-type strains and may have been selected early in domestication. This is the first dissection of the physiological consequences of a genetic locus for a quantitative trait that alters appetite and gives us an insight into the domestication of animals. This will allow elucidation of how differences in appetite occur in birds and also mammals.

JEnsembl: a version-aware Java API to Ensembl data systems.

MOTIVATION: The Ensembl Project provides release-specific Perl APIs for efficient high-level programmatic access to data stored in various Ensembl database schema. Although Perl scripts are perfectly suited for processing large volumes of text-based data, Perl is not ideal for developing large-scale software applications nor embedding in graphical interfaces. The provision of a novel Java API would facilitate type-safe, modular, object-orientated development of new Bioinformatics tools with which to access, analyse and visualize Ensembl data. RESULTS: The JEnsembl API implementation provides basic data retrieval and manipulation functionality from the Core, Compara and Variation databases for all species in Ensembl and EnsemblGenomes and is a platform for the development of a richer API to Ensembl datasources. The JEnsembl architecture uses a text-based configuration module to provide evolving, versioned mappings from database schema to code objects. A single installation of the JEnsembl API can therefore simultaneously and transparently connect to current and previous database instances (such as those in the public archive) thus facilitating better analysis repeatability and allowing 'through time' comparative analyses to be performed. AVAILABILITY: Project development, released code libraries, Maven repository and documentation are hosted at SourceForge (http://jensembl.sourceforge.net).

VIPER: a visualisation tool for exploring inheritance inconsistencies in genotyped pedigrees.

BACKGROUND: Pedigree genotype datasets are used for analysing genetic inheritance and to map genetic markers and traits. Such datasets consist of hundreds of related animals genotyped for thousands of genetic markers and invariably contain multiple errors in both the pedigree structure and in the associated individual genotype data. These errors manifest as apparent inheritance inconsistencies in the pedigree, and invalidate analyses of marker inheritance patterns across the dataset. Cleaning raw datasets of bad data points (incorrect pedigree relationships, unreliable marker assays, suspect samples, bad genotype results etc.) requires expert exploration of the patterns of exposed inconsistencies in the context of the inheritance pedigree. In order to assist this process we are developing VIPER (Visual Pedigree Explorer), a software tool that integrates an inheritance-checking algorithm with a novel space-efficient pedigree visualisation, so that reported inheritance inconsistencies are overlaid on an interactive, navigable representation of the pedigree structure. METHODS AND RESULTS: This paper describes an evaluation of how VIPER displays the different scales and types of dataset that occur experimentally, with a description of how VIPER's display interface and functionality meet the challenges presented by such data. We examine a range of possible error types found in real and simulated pedigree genotype datasets, demonstrating how these errors are exposed and explored using the VIPER interface and we evaluate the utility and usability of the interface to the domain expert.Evaluation was performed as a two stage process with the assistance of domain experts (geneticists). The initial evaluation drove the iterative implementation of further features in the software prototype, as required by the users, prior to a final functional evaluation of the pedigree display for exploring the various error types, data scales and structures. CONCLUSIONS: The VIPER display was shown to effectively expose the range of errors found in experimental genotyped pedigrees, allowing users to explore the underlying causes of reported inheritance inconsistencies. This interface will provide the basis for a full data cleaning tool that will allow the user to remove isolated bad data points, and reversibly test the effect of removing suspect genotypes and pedigree relationships.

A novel strain of lumpy skin disease virus causes clinical disease in cattle in Hong Kong.

Lumpy skin disease virus (LSDV) is an emerging poxviral pathogen of cattle that is currently spreading throughout Asia. The disease situation is of high importance for farmers and policy makers in Asia. In October 2020, feral cattle in Hong Kong developed multi-focal cutaneous nodules consistent with lumpy skin disease (LSD). Gross and histological pathology further supported the diagnosis and samples were sent to the OIE Reference Laboratory at The Pirbright Institute for confirmatory testing. LSDV was detected using quantitative polymerase chain reaction (qPCR) and additional molecular analyses. This is the first report of LSD in Hong Kong. Whole genome sequencing (WGS) of the strain LSDV/Hong Kong/2020 and phylogenetic analysis were carried out in order to identify connections to previous outbreaks of LSD, and better understand the drivers of LSDV emergence. Analysis of the 90 core poxvirus genes revealed LSDV/Hong Kong/2020 was a novel strain most closely related to the live-attenuated Neethling vaccine strains of LSDV and more distantly related to wildtype LSDV isolates from Africa, the Middle East and Europe. Analysis of the more variable regions located towards the termini of the poxvirus genome revealed genes in LSDV/Hong Kong/2020 with different patterns of grouping when compared to previously published wildtype and vaccine strains of LSDV. This work reveals that the LSD outbreak in Hong Kong in 2020 was caused by a different strain of LSDV than the LSD epidemic in the Middle East and Europe in 2015-2018. The use of WGS is highly recommended when investigating LSDV disease outbreaks.

Evaluation of approaches for identifying population informative markers from high density SNP chips.

BACKGROUND: Genetic markers can be used to identify and verify the origin of individuals. Motivation for the inference of ancestry ranges from conservation genetics to forensic analysis. High density assays featuring Single Nucleotide Polymorphism (SNP) markers can be exploited to create a reduced panel containing the most informative markers for these purposes. The objectives of this study were to evaluate methods of marker selection and determine the minimum number of markers from the BovineSNP50 BeadChip required to verify the origin of individuals in European cattle breeds. Delta, Wright's FST, Weir & Cockerham's FST and PCA methods for population differentiation were compared. The level of informativeness of each SNP was estimated from the breed specific allele frequencies. Individual assignment analysis was performed using the ranked informative markers. Stringency levels were applied by log-likelihood ratio to assess the confidence of the assignment test. RESULTS: A 95% assignment success rate for the 384 individually genotyped animals was achieved with <80, <100, <140 and <200 SNP markers (with increasing stringency threshold levels) across all the examined methods for marker selection. No further gain in power of assignment was achieved by sampling in excess of 200 SNP markers. The marker selection method that required the lowest number of SNP markers to verify the animal's breed origin was Wright's FST (60 to 140 SNPs depending on the chosen degree of confidence). Certain breeds required fewer markers (<100) to achieve 100% assignment success. In contrast, closely related breeds require more markers (~200) to achieve>95% assignment success. The power of assignment success, and therefore the number of SNP markers required, is dependent on the levels of genetic heterogeneity and pool of samples considered. CONCLUSIONS: While all SNP selection methods produced marker panels capable of breed identification, the power of assignment varied markedly among analysis methods. Thus, with effective exploration of available high density genetic markers, a diagnostic panel of highly informative markers can be produced.

Suppression of HDAC nuclear export and cardiomyocyte hypertrophy by novel irreversible inhibitors of CRM1.

Histone deacetylase 5 (HDAC5) represses expression of nuclear genes that promote cardiac hypertrophy. Agonism of a variety of G protein coupled receptors (GPCRs) triggers phosphorylation-dependent nuclear export of HDAC5 via the CRM1 nuclear export receptor, resulting in derepression of pro-hypertrophic genes. A cell-based high-throughput screen of a commercial compound collection was employed to identify compounds with the ability to preserve the nuclear fraction of GFP-HDAC5 in primary cardiomyocytes exposed to GPCR agonists. A hit compound potently inhibited agonist-induced GFP-HDAC5 nuclear export in cultured neonatal rat ventricular myocytes (NRVMs). A small set of related compounds was designed and synthesized to evaluate structure-activity relationship (SAR). The results demonstrated that inhibition of HDAC5 nuclear export was a result of compounds irreversibly reacting with a key cysteine residue in CRM1 that is required for its function. CRM1 inhibition by the compounds also resulted in potent suppression of cardiomyocyte hypertrophy. These studies define a novel class of anti-hypertrophic compounds that function through irreversible inhibition of CRM1-dependent nuclear export.

Dynamic data-driven meta-analysis for prioritisation of host genes implicated in COVID-19.

The increasing body of literature describing the role of host factors in COVID-19 pathogenesis demonstrates the need to combine diverse, multi-omic data to evaluate and substantiate the most robust evidence and inform development of therapies. Here we present a dynamic ranking of host genes implicated in human betacoronavirus infection (SARS-CoV-2, SARS-CoV, MERS-CoV, seasonal coronaviruses). We conducted an extensive systematic review of experiments identifying potential host factors. Gene lists from diverse sources were integrated using Meta-Analysis by Information Content (MAIC). This previously described algorithm uses data-driven gene list weightings to produce a comprehensive ranked list of implicated host genes. From 32 datasets, the top ranked gene was PPIA, encoding cyclophilin A, a druggable target using cyclosporine. Other highly-ranked genes included proposed prognostic factors (CXCL10, CD4, CD3E) and investigational therapeutic targets (IL1A) for COVID-19. Gene rankings also inform the interpretation of COVID-19 GWAS results, implicating FYCO1 over other nearby genes in a disease-associated locus on chromosome 3. Researchers can search and review the gene rankings and the contribution of different experimental methods to gene rank at https://baillielab.net/maic/covid19 . As new data are published we will regularly update the list of genes as a resource to inform and prioritise future studies.

An XML transfer schema for exchange of genomic and genetic mapping data: implementation as a web service in a Taverna workflow.

BACKGROUND: Genomic analysis, particularly for less well-characterized organisms, is greatly assisted by performing comparative analyses between different types of genome maps and across species boundaries. Various providers publish a plethora of on-line resources collating genome mapping data from a multitude of species. Datasources range in scale and scope from small bespoke resources for particular organisms, through larger web-resources containing data from multiple species, to large-scale bioinformatics resources providing access to data derived from genome projects for model and non-model organisms. The heterogeneity of information held in these resources reflects both the technologies used to generate the data and the target users of each resource. Currently there is no common information exchange standard or protocol to enable access and integration of these disparate resources. Consequently data integration and comparison must be performed in an ad hoc manner. RESULTS: We have developed a simple generic XML schema (GenomicMappingData.xsd - GMD) to allow export and exchange of mapping data in a common lightweight XML document format. This schema represents the various types of data objects commonly described across mapping datasources and provides a mechanism for recording relationships between data objects. The schema is sufficiently generic to allow representation of any map type (for example genetic linkage maps, radiation hybrid maps, sequence maps and physical maps). It also provides mechanisms for recording data provenance and for cross referencing external datasources (including for example ENSEMBL, PubMed and Genbank.). The schema is extensible via the inclusion of additional datatypes, which can be achieved by importing further schemas, e.g. a schema defining relationship types. We have built demonstration web services that export data from our ArkDB database according to the GMD schema, facilitating the integration of data retrieval into Taverna workflows. CONCLUSION: The data exchange standard we present here provides a useful generic format for transfer and integration of genomic and genetic mapping data. The extensibility of our schema allows for inclusion of additional data and provides a mechanism for typing mapping objects via third party standards. Web services retrieving GMD-compliant mapping data demonstrate that use of this exchange standard provides a practical mechanism for achieving data integration, by facilitating syntactically and semantically-controlled access to the data.

An introduction to scripting in Ruby for biologists.

The Ruby programming language has a lot to offer to any scientist with electronic data to process. Not only is the initial learning curve very shallow, but its reflection and meta-programming capabilities allow for the rapid creation of relatively complex applications while still keeping the code short and readable. This paper provides a gentle introduction to this scripting language for researchers without formal informatics training such as many wet-lab scientists. We hope this will provide such researchers an idea of how powerful a tool Ruby can be for their data management tasks and encourage them to learn more about it.