Administration and detection of a multi-target rAAV gene doping vector in horses using multiple matrices and molecular techniques.

Gene doping, which includes the non-therapeutic use of genes or genetic elements that have the capacity to enhance athletic performance, is prohibited in horseracing and equestrian sports. To provide a comprehensive assessment of matrix and detection techniques, a custom adeno-associated virus serotype 8 vector was designed to include PCR binding sites for multiple target genes and assay types. The vector was injected via an intramuscular route into two Thoroughbred horses and matrices collected at defined timepoints. DNA was analysed using 3 detection methods: qPCR, digital PCR, and NGS. Overall, there was a strong correlation across the different detection methods employed, although digital PCR was less sensitive at lower concentrations. High concentrations of vector were detected at early timepoints in plasma and whole blood, which rapidly dropped after 0.5 d to trace levels by 4 d and 9 d post-administration respectively, following a similar pattern to previous studies. Vector was detected in dried blood spots at lower levels than whole blood, but with a similar detection time. Detection in hair root bulbs in one horse was observed at over a month post-administration, which opens new avenues for future gene doping testing in humans and animals.

Detection of transgenes in equine dried blood spots using digital PCR and qPCR for gene doping control.

Due to the ease of collection, transport and storage, the use of dried blood spots (DBS) offers an attractive alternative matrix for detection of the abuse of gene therapy, otherwise known as gene doping. This study evaluated the recovery, extraction efficiency and resulting detection capability of DNA from DBS by evaluating different target types, DNA extraction kits, the number of punches and blood tube preservatives. The long-term storage stability of low-copy-number transgene targets in DBS was not assessed in this study but would be noteworthy to investigate further. DNA was quantified using two detection methods: qPCR and digital PCR (dPCR). Using six punches with the Qiagen Investigator kit gave the best overall DNA yield compared with other extraction methods. Including three punches, however, gave better DNA extraction efficiency. Reference material could be detected using qPCR and dPCR in DBS spiked with 5000 copies/mL of blood (approximately 15 copies per 3 mm of punch). The optimal DNA extraction protocol was used on DBS samples from a custom recombinant adeno-associated virus administration study and showed successful detection of vector targets in DBS samples.

Analytical advances in horseracing medication and doping control from 2018 to 2023.

The analytical approaches taken by laboratories to implement robust and efficient regulation of horseracing medication and doping control are complex and constantly evolving. Each laboratory's approach will be dictated by differences in regulatory, economic and scientific drivers specific to their local environment. However, in general, laboratories will all be undertaking developments and improvements to their screening strategies in order to meet new and emerging threats as well as provide improved service to their customers. In this paper, the published analytical advances in horseracing medication and doping control since the 22nd International Conference of Racing Analysts and Veterinarians will be reviewed. Due to the unprecedented impact of COVID-19 on the worldwide economy, the normal 2-year period of this review was extended to over 5 years. As such, there was considerable ground to cover, resulting in an increase in the number of relevant publications included from 107 to 307. Major trends in publications will be summarised and possible future directions highlighted. This will cover developments in the detection of 'small' and 'large' molecule drugs, sample preparation procedures and the use of alternative matrices, instrumental advances/applications, drug metabolism and pharmacokinetics, the detection and prevalence of 'endogenous' compounds and biomarker and OMICs approaches. Particular emphasis will be given to research into the potential threat of gene doping, which is a significant area of new and continued research for many laboratories. Furthermore, developments in analytical instrumentation relevant to equine medication and doping control will be discussed.

NUDCD3 deficiency disrupts V(D)J recombination to cause SCID and Omenn syndrome.

Inborn errors of T cell development present a pediatric emergency in which timely curative therapy is informed by molecular diagnosis. In 11 affected patients across four consanguineous kindreds, we detected homozygosity for a single deleterious missense variant in the gene NudC domain-containing 3 (NUDCD3). Two infants had severe combined immunodeficiency with the complete absence of T and B cells (T -B- SCID), whereas nine showed classical features of Omenn syndrome (OS). Restricted antigen receptor gene usage by residual T lymphocytes suggested impaired V(D)J recombination. Patient cells showed reduced expression of NUDCD3 protein and diminished ability to support RAG-mediated recombination in vitro, which was associated with pathologic sequestration of RAG1 in the nucleoli. Although impaired V(D)J recombination in a mouse model bearing the homologous variant led to milder immunologic abnormalities, NUDCD3 is absolutely required for healthy T and B cell development in humans.

Detection of adeno-associated viral DNA in equine post-administration frozen blood and plasma samples after long-term storage.

Gene doping in horses is a threat to the fairness in sport and has serious implications for animal welfare. To investigate the effect of long-term storage on the detection of AAV in plasma and whole blood, samples from an administration study using an adeno-associated virus serotype 6 expressing green fluorescence protein (AAV6-GFP) were stored at -20°C for 8 months before analysis. The AAV vector was detected in stored plasma samples, following the same detection profile as the fresh plasma samples. The stored blood showed lower overall DNA detection but followed the same detection profile as the plasma samples. This study provides confidence that re-analysing plasma samples and/or analysing a frozen 'B' sample with different matrix such as whole blood after prolonged storage will still result in the detection of gene doping material.

Long-read sequencing assays designed to detect potential gene editing events in the myostatin gene revealed distinct haplotype signatures in the Thoroughbred horse population.

We present here the use of targeted, long-read sequencing of the myostatin (MSTN) gene as a model to detect potential gene editing events in Thoroughbred horses. MSTN is a negative regulator of muscle development, making the gene a prime candidate target for gene doping. By sequencing the complete gene in one PCR product, we can catalogue all mutations without the need to produce short-fragment libraries. A panel of reference material fragments with defined mutations was constructed and successfully sequenced by both Oxford Nanopore and Illumina-based methods, showing that gene doping editing events can be detected using this technology. To ascertain the normal variation within the population, we sequenced the MSTN gene in 119 UK Thoroughbred horses. Variants from the reference genome were assigned to haplotypes and eight distinct patterns, designated Hap1 (reference genome) to Hap8, were determined with haplotypes Hap2 and Hap3 (which includes the 'speed gene' variant) being far the most prevalent. Hap3 was most abundant in flat-racing horses, whereas Hap2 was most abundant in jump-racing. Within this data set, results for 105 racehorses from out-of-competition sampling were compared between matrices of extracted DNA and direct PCR of whole blood from lithium heparin gel tubes, and strong agreement was found between the two methods. The direct-blood PCR was achieved without compromising the sample prior to plasma separation for analytical chemistry, and could thus be used as part of a routine screening workflow for gene editing detection.

Use of mitochondrial sequencing to detect gene doping in horses via gene editing and somatic cell nuclear transfer.

Gene editing and subsequent cloning techniques offer great potential not only in genetic disease correction in domestic animals but also in livestock production by enhancement of desirable traits. The existence of the technology, however, leaves it open to potential misuse in performance-led sports such as horseracing and other equestrian events. Recent advances in equine gene editing, regarding the generation of gene-edited embryos using CRISPR/Cas9 technology and somatic cell nuclear transfer, have highlighted the need to develop tools to detect potential prohibited use of the technology. One possible method involves the characterisation of the mitochondrial genome (which is not routinely preserved during cloning) and comparing it with the sequence of the registered dam. We present here our approach to whole-mitochondrial sequencing using tiled long-range PCR and next-generation sequencing. To determine whether the background mutation rate in the mitochondrial genome could potentially confound results, we sequenced 10 sets of dam and foal duos. We found variation between duos but none within duos, indicating that this method is feasible for future screening systems. Analysis of WGS data from over 100 Thoroughbred horses revealed wide variation in the mitochondria sequence within the breed, further displaying the utility of this approach.

Direct sequence confirmation of qPCR products for gene doping assay validation in horses.

The misuse of gene therapy by the introduction of transgenes via plasmid or viral vectors as a doping agent is an increasing concern in human and animal sports, not only in consideration to fair competition but also in potential detrimental effects to welfare. Doping events can be detected by polymerase chain reaction (PCR) amplification of a transgene-specific region of DNA. Quantitative real-time PCR (qPCR) is particularly suited to confirmatory investigations where precise limits of detection can be calculated. To fully validate a qPCR experiment, it is highly desirable to confirm the identity of the amplicon. Although post-PCR techniques such as melt curve and fragment size analysis can provide strong evidence that the amplicon is as expected, sequence identity confirmation may be beneficial as part of regulatory proceedings. We present here our investigation into two alternative processes for the direct assessment of qPCR products for five genes using next-generation sequencing: ligation of sequence-ready adapters to qPCR products and qPCR assays performed with primers tailed with Illumina flow cell binding sites. To fully test the robustness of the techniques at concentrations required for gene doping detection, we also calculated a putative limit of detection for the assays. Both ligated adapters and tailed primers were successful in producing sequence data for the qPCR products without further amplification. Ligated adapters are preferred, however, as they do not require re-optimisation of existing qPCR assays.

Screening for gene doping transgenes in horses via the use of massively parallel sequencing.

Throughout the history of horse racing, doping techniques to suppress or enhance performance have expanded to match the technology available. The next frontier in doping, both in the equine and human sports areas, is predicted to be genetic manipulation; either by prohibited use of genome editing, or gene therapy via transgenes. By using massively-parallel sequencing via a two-step PCR method we can screen for multiple doping targets at once in pooled primer sets. This method has the advantages of high scalability through combinational indexing, and the use of reference standards with altered sequences as controls. Custom software produces transgene-specific amplicons from any Ensembl-annotated genome to facilitate rapid assay design. Additional scripts batch-process FASTQ data from experiments, automatically quality-filtering sequences and assigning hits based on discriminatory motifs. We report here our experiences in establishing the workflow with an initial 31 transgene and vector feature targets. To evaluate the sensitivity of parallel sequencing in a real-world setting, we performed an intramuscular (IM) administration of a control rAAV vector into two horses and compared the detection sensitivity between parallel sequencing and real-time qPCR. Vector was detected by all assays on both methods up to 79 h post-administration, becoming sporadic after 96 h.

High-throughput genotyping of high-homology mutant mouse strains by next-generation sequencing.

Genotyping of knockout alleles in mice is commonly performed by end-point PCR or gene-specific/universal cassette qPCR. Both have advantages and limitations in terms of assay design and interpretation of results. As an alternative method for high-throughput genotyping, we investigated next generation sequencing (NGS) of PCR amplicons, with a focus on CRISPR-mediated exon deletions where antibiotic selection markers are not present. By multiplexing the wild type and mutant-specific PCR reactions, the genotype can be called by the relative sequence counts of each product. The system is highly scalable and can be applied to a variety of different allele types, including those produced by the International Mouse Phenotyping Consortium and associated projects. One potential challenge with any assay design is locating unique areas of the genome, especially when working with gene families or regions of high homology. These can result in misleading or ambiguous genotypes for either qPCR or end-point assays. Here, we show that genotyping by NGS can negate these issues by simple, automated filtering of undesired sequences. Analysis and genotype calls can also be fully automated, using FASTQ or FASTA input files and an in-house Perl script and SQL database.

Tumors induce de novo steroid biosynthesis in T cells to evade immunity.

Tumors subvert immune cell function to evade immune responses, yet the complex mechanisms driving immune evasion remain poorly understood. Here we show that tumors induce de novo steroidogenesis in T lymphocytes to evade anti-tumor immunity. Using a transgenic steroidogenesis-reporter mouse line we identify and characterize de novo steroidogenic immune cells, defining the global gene expression identity of these steroid-producing immune cells and gene regulatory networks by using single-cell transcriptomics. Genetic ablation of T cell steroidogenesis restricts primary tumor growth and metastatic dissemination in mouse models. Steroidogenic T cells dysregulate anti-tumor immunity, and inhibition of the steroidogenesis pathway is sufficient to restore anti-tumor immunity. This study demonstrates T cell de novo steroidogenesis as a mechanism of anti-tumor immunosuppression and a potential druggable target.

Common and distinct transcriptional signatures of mammalian embryonic lethality.

The Deciphering the Mechanisms of Developmental Disorders programme has analysed the morphological and molecular phenotypes of embryonic and perinatal lethal mouse mutant lines in order to investigate the causes of embryonic lethality. Here we show that individual whole-embryo RNA-seq of 73 mouse mutant lines (>1000 transcriptomes) identifies transcriptional events underlying embryonic lethality and associates previously uncharacterised genes with specific pathways and tissues. For example, our data suggest that Hmgxb3 is involved in DNA-damage repair and cell-cycle regulation. Further, we separate embryonic delay signatures from mutant line-specific transcriptional changes by developing a baseline mRNA expression catalogue of wild-type mice during early embryogenesis (4-36 somites). Analysis of transcription outside coding sequence identifies deregulation of repetitive elements in Morc2a mutants and a gene involved in gene-specific splicing. Collectively, this work provides a large scale resource to further our understanding of early embryonic developmental disorders.

FBXO7 sensitivity of phenotypic traits elucidated by a hypomorphic allele.

FBXO7 encodes an F box containing protein that interacts with multiple partners to facilitate numerous cellular processes and has a canonical role as part of an SCF E3 ubiquitin ligase complex. Mutation of FBXO7 is responsible for an early onset Parkinsonian pyramidal syndrome and genome-wide association studies have linked variants in FBXO7 to erythroid traits. A putative orthologue in Drosophila, nutcracker, has been shown to regulate the proteasome, and deficiency of nutcracker results in male infertility. Therefore, we reasoned that modulating Fbxo7 levels in a murine model could provide insights into the role of this protein in mammals. We used a targeted gene trap model which retained 4-16% residual gene expression and assessed the sensitivity of phenotypic traits to gene dosage. Fbxo7 hypomorphs showed regenerative anaemia associated with a shorter erythrocyte half-life, and male mice were infertile. Alterations to T cell phenotypes were also observed, which intriguingly were both T cell intrinsic and extrinsic. Hypomorphic mice were also sensitive to infection with Salmonella, succumbing to a normally sublethal challenge. Despite these phenotypes, Fbxo7 hypomorphs were produced at a normal Mendelian ratio with a normal lifespan and no evidence of neurological symptoms. These data suggest that erythrocyte survival, T cell development and spermatogenesis are particularly sensitive to Fbxo7 gene dosage.

Mammalian Hbs1L deficiency causes congenital anomalies and developmental delay associated with Pelota depletion and 80S monosome accumulation.

Hbs1 has been established as a central component of the cell's translational quality control pathways in both yeast and prokaryotic models; however, the functional characteristics of its human ortholog (Hbs1L) have not been well-defined. We recently reported a novel human phenotype resulting from a mutation in the critical coding region of the HBS1L gene characterized by facial dysmorphism, severe growth restriction, axial hypotonia, global developmental delay and retinal pigmentary deposits. Here we further characterize downstream effects of the human HBS1L mutation. HBS1L has three transcripts in humans, and RT-PCR demonstrated reduced mRNA levels corresponding with transcripts V1 and V2 whereas V3 expression was unchanged. Western blot analyses revealed Hbs1L protein was absent in the patient cells. Additionally, polysome profiling revealed an abnormal aggregation of 80S monosomes in patient cells under baseline conditions. RNA and ribosomal sequencing demonstrated an increased translation efficiency of ribosomal RNA in Hbs1L-deficient fibroblasts, suggesting that there may be a compensatory increase in ribosome translation to accommodate the increased 80S monosome levels. This enhanced translation was accompanied by upregulation of mTOR and 4-EBP protein expression, suggesting an mTOR-dependent phenomenon. Furthermore, lack of Hbs1L caused depletion of Pelota protein in both patient cells and mouse tissues, while PELO mRNA levels were unaffected. Inhibition of proteasomal function partially restored Pelota expression in human Hbs1L-deficient cells. We also describe a mouse model harboring a knockdown mutation in the murine Hbs1l gene that shared several of the phenotypic elements observed in the Hbs1L-deficient human including facial dysmorphism, growth restriction and retinal deposits. The Hbs1lKO mice similarly demonstrate diminished Pelota levels that were rescued by proteasome inhibition.

No unexpected CRISPR-Cas9 off-target activity revealed by trio sequencing of gene-edited mice.

CRISPR-Cas9 technologies have transformed genome-editing of experimental organisms and have immense therapeutic potential. Despite significant advances in our understanding of the CRISPR-Cas9 system, concerns remain over the potential for off-target effects. Recent studies have addressed these concerns using whole-genome sequencing (WGS) of gene-edited embryos or animals to search for de novo mutations (DNMs), which may represent candidate changes introduced by poor editing fidelity. Critically, these studies used strain-matched, but not pedigree-matched controls and thus were unable to reliably distinguish generational or colony-related differences from true DNMs. Here we used a trio design and whole genome sequenced 8 parents and 19 embryos, where 10 of the embryos were mutagenised with well-characterised gRNAs targeting the coat colour Tyrosinase (Tyr) locus. Detailed analyses of these whole genome data allowed us to conclude that if CRISPR mutagenesis were causing SNV or indel off-target mutations in treated embryos, then the number of these mutations is not statistically distinguishable from the background rate of DNMs occurring due to other processes.

Placentation defects are highly prevalent in embryonic lethal mouse mutants.

Large-scale phenotyping efforts have demonstrated that approximately 25-30% of mouse gene knockouts cause intrauterine lethality. Analysis of these mutants has largely focused on the embryo and not the placenta, despite the crucial role of this extraembryonic organ for developmental progression. Here we screened 103 embryonic lethal and sub-viable mouse knockout lines from the Deciphering the Mechanisms of Developmental Disorders program for placental phenotypes. We found that 68% of knockout lines that are lethal at or after mid-gestation exhibited placental dysmorphologies. Early lethality (embryonic days 9.5-14.5) is almost always associated with severe placental malformations. Placental defects correlate strongly with abnormal brain, heart and vascular development. Analysis of mutant trophoblast stem cells and conditional knockouts suggests that a considerable number of factors that cause embryonic lethality when ablated have primary gene function in trophoblast cells. Our data highlight the hugely under-appreciated importance of placental defects in contributing to abnormal embryo development and suggest key molecular nodes that govern placenta formation.

Corrigendum: High-throughput discovery of novel developmental phenotypes.

This corrects the article DOI: 10.1038/nature19356.

Genome-wide in vivo screen identifies novel host regulators of metastatic colonization.

Metastasis is the leading cause of death for cancer patients. This multi-stage process requires tumour cells to survive in the circulation, extravasate at distant sites, then proliferate; it involves contributions from both the tumour cell and tumour microenvironment ('host', which includes stromal cells and the immune system). Studies suggest the early steps of the metastatic process are relatively efficient, with the post-extravasation regulation of tumour growth ('colonization') being critical in determining metastatic outcome. Here we show the results of screening 810 mutant mouse lines using an in vivo assay to identify microenvironmental regulators of metastatic colonization. We identify 23 genes that, when disrupted in mouse, modify the ability of tumour cells to establish metastatic foci, with 19 of these genes not previously demonstrated to play a role in host control of metastasis. The largest reduction in pulmonary metastasis was observed in sphingosine-1-phosphate (S1P) transporter spinster homologue 2 (Spns2)-deficient mice. We demonstrate a novel outcome of S1P-mediated regulation of lymphocyte trafficking, whereby deletion of Spns2, either globally or in a lymphatic endothelial-specific manner, creates a circulating lymphopenia and a higher percentage of effector T cells and natural killer (NK) cells present in the lung. This allows for potent tumour cell killing, and an overall decreased metastatic burden.

High-throughput discovery of novel developmental phenotypes.

Approximately one-third of all mammalian genes are essential for life. Phenotypes resulting from knockouts of these genes in mice have provided tremendous insight into gene function and congenital disorders. As part of the International Mouse Phenotyping Consortium effort to generate and phenotypically characterize 5,000 knockout mouse lines, here we identify 410 lethal genes during the production of the first 1,751 unique gene knockouts. Using a standardized phenotyping platform that incorporates high-resolution 3D imaging, we identify phenotypes at multiple time points for previously uncharacterized genes and additional phenotypes for genes with previously reported mutant phenotypes. Unexpectedly, our analysis reveals that incomplete penetrance and variable expressivity are common even on a defined genetic background. In addition, we show that human disease genes are enriched for essential genes, thus providing a dataset that facilitates the prioritization and validation of mutations identified in clinical sequencing efforts.