Multiplex Detection of Transgenes Using πCode Technology for Gene Doping Control.

To ensure fair competition and sports integrity, gene doping is prohibited in horseracing and equine sports. One gene doping method is by administering exogenous genes, called transgenes, to postnatal animals. Although several transgene detection methods have been developed for horses, many are unsuitable for multiplex detection. In this proof-of-concept study, we developed a highly sensitive and multiplex transgene detection method using multiple πCode with identification patterns printed on the surface. The following steps were employed: (1) multiplex polymerase chain reaction amplification of 12 targeted transgenes in a single tube, (2) detection using a mixture of 12 probes labeled with different πCodes, and (3) median fluorescence intensity measurement of fluorescent πCodes. Twelve transgenes cloned into plasmid vectors were targeted, and 1500 copies of each plasmid were spiked into 1.5 mL of horse plasma. Subsequently, a novel method using πCode succeeded in detecting all the transgenes using their DNA extracts. Additionally, we detected the erythropoietin (EPO) transgene in blood samples from a horse administered solely with the EPO transgene using this method. Therefore, the πCode detection method is considered suitable for multitarget gene detection in gene doping tests.

Fpr1, a primary target of rapamycin, functions as a transcription factor for ribosomal protein genes cooperatively with Hmo1 in Saccharomyces cerevisiae.

Fpr1 (FK506-sensitive proline rotamase 1), a protein of the FKBP12 (FK506-binding protein 12 kDa) family in Saccharomyces cerevisiae, is a primary target for the immunosuppressive agents FK506 and rapamycin. Fpr1 inhibits calcineurin and TORC1 (target of rapamycin complex 1) when bound to FK506 and rapamycin, respectively. Although Fpr1 is recognised to play a crucial role in the efficacy of these drugs, its physiological functions remain unclear. In a hmo1Δ (high mobility group family 1-deleted) yeast strain, deletion of FPR1 induced severe growth defects, which could be alleviated by increasing the copy number of RPL25 (ribosome protein of the large subunit 25), suggesting that RPL25 expression was affected in hmo1Δfpr1Δ cells. In the current study, extensive chromatin immunoprecipitation (ChIP) and ChIP-sequencing analyses revealed that Fpr1 associates specifically with the upstream activating sequences of nearly all RPG (ribosomal protein gene) promoters, presumably in a manner dependent on Rap1 (repressor/activator site binding protein 1). Intriguingly, Fpr1 promotes the binding of Fhl1/Ifh1 (forkhead-like 1/interacts with forkhead 1), two key regulators of RPG transcription, to certain RPG promoters independently of and/or cooperatively with Hmo1. Furthermore, mutation analyses of Fpr1 indicated that for transcriptional function on RPG promoters, Fpr1 requires its N-terminal domain and the binding surface for rapamycin, but not peptidyl-prolyl isomerase activity. Notably, Fpr1 orthologues from other species also inhibit TORC1 when bound to rapamycin, but do not regulate transcription in yeast, which suggests that these two functions of Fpr1 are independent of each other.

Identification of Potential miRNA Biomarkers to Detect Hydrocortisone Administration in Horses.

Circulating microRNAs (miRNAs) are stable in bodily fluids and are potential biomarkers of various diseases and physiological states. Although several studies have been conducted on humans to detect drug doping by miRNAs, research on drugs and miRNAs in horses is limited. In this study, circulating miRNAs in horses after hydrocortisone administration were profiled and variations in miRNAs affected by hydrocortisone administration during endogenous hydrocortisone elevation were examined. The miRNAs were extracted from thoroughbred horse plasma before and after hydrocortisone administration and subjected to small RNA sequencing and reverse transcription quantitative PCR (RT-qPCR). RT-qPCR validation was performed for the 20 miRNAs that were most affected by hydrocortisone administration. The effects of elevated endogenous hydrocortisone levels due to exercise and adrenocorticotropic hormone administration were also confirmed. The validation results showed that approximately half of the miRNAs showed the same significant differences as those obtained using small RNA sequencing. Among the twenty miRNAs, two novel miRNAs and miR-133a were found to vary differently between exogenous hydrocortisone administration and endogenous hydrocortisone elevation. This study provides basic knowledge regarding the circulating miRNA profile of horses after hydrocortisone administration and identifies three miRNAs that could potentially be used as biomarkers to detect hydrocortisone administration.

Monitoring the positive conversion of anti-erythrocyte antibodies in blood transfusion donor horses.

To confirm the positive conversion of antibodies against erythrocyte antigens in horses, possible blood transfusion donor horses selected from draft horse populations were periodically monitored with an indirect antiglobulin (Coombs) test for approximately 3 years. In this study, 19 horses (16 females and 3 males) were investigated, and five mares showed alloantibodies during the monitoring period. Four mares were typically pregnant when positive conversion was detected, whereas no particular cause of conversion could be observed for one mare based on its clinical records. In the analyzed horses, most positive conversions were possibly due to pregnancy, as conversion occurred more often during this period than after parturition. Pregnancy is considered a key event for positive conversion. Additionally, in cases in which unknown causative sensitization is confirmed, continuous monitoring with a test to detect antibodies should be performed, even if the possible donor is selected and maintained.

Construction of an individual identification panel for horses using insertion and deletion markers.

Individual identification and paternity testing are important for avoiding inbreeding in the management of small populations of wild and domestic animals. In horse racing industries, they are extremely important for identifying and registering individuals and doping control to ensure fair competition. In this study, we constructed an individual identification panel for horses by using insertion and deletion (INDEL) markers. The panel included 39 INDEL markers selected from a whole-genome INDEL database. Genotyping of 89 Thoroughbreds showed polymorphisms with minor allele frequencies (MAFs) of 0.180-0.489 in all markers. The total probability of exclusion for paternity testing, power of discrimination, and probability of identity were 0.9994271269, >0.9999999999, and 0.9999999987, respectively. The panel was applied to 13 trios (sires, dams, and foals), and no contradictions were observed in genetic inheritance among the trios. When this panel was applied to the trios (52 trios) containing false fathers, an average of 7.3 markers excluded parentage relationships. In addition, genomic DNA extracted from the urine of six horses was partially genotyped for 39 markers, and 6-28 markers were successfully genotyped. The newly constructed panel has two advantages: a low marker mutation rate compared with short tandem repeats and a genotyping procedure that is as simple as short tandem repeat typing compared with single nucleotide variant typing. This panel can be applied for individual identification, paternity determination, and urine-sample identification in Thoroughbred horses.

Investigation of optimal procedures for storage and use of plasma samples suitable for gene doping tests.

Gene doping, which is prohibited in horseracing and equestrian sports, can be performed by introducing exogenous genes, known as transgenes, into the bodies of postnatal animals. To detect exogenous genes, a method utilizing quantitative polymerase chain reaction (qPCR) with a hydrolysis probe was developed to test whole blood and plasma samples, thereby protecting the fairness of competition and the rights of stakeholders in horseracing and equestrian sports. Therefore, we aimed to develop sample storage methods suitable for A and B samples in gene doping tests using blood. For sample A, sufficient qPCR detection was demonstrated after refrigeration for 1 to 2 weeks post collection. For sample B, the following procedures were confirmed to be suitable for storage: 1) centrifugation after sample receipt, 2) frozen storage, 3) natural thawing at room temperature, and 4) centrifugation without mixing blood cell components. Our results indicated that long-term cryopreservation yielded good plasma components from frozen blood samples even though it destroyed blood cells, indicating its applicability to the gene doping test using sample B, which can be stored for later use. Sample storage procedures are as important as detection methods in doping tests. Therefore, the series of procedures that we evaluated in this study will contribute to the efficient performance of gene doping tests through qPCR using blood samples.

Identification of processed pseudogenes in the genome of Thoroughbred horses: Possibility of gene-doping detection considering the presence of pseudogenes.

Processed pseudogenes, also known as retrocopy genes, are copies of messenger RNAs that have been reverse transcribed into DNA and inserted into the genome. In this study, we identified 62 processed pseudogene candidates as intron-less genes from whole-genome sequencing (WGS) data of Thoroughbred horses using delly structural variation software. The 62 processed pseudogene candidates were confirmed by PCR amplification of intron-less products. A total of 11 processed pseudogenes were confirmed in the genome of all 23 analysed horses, whereas three processed pseudogenes with structures of ATP11B, DPH3 and RPL17 were detected in only one of 115 horses by PCR amplification of intron-less products. Currently, most of the gene doping tests proposed in human and horse sports are adapted PCR-based methods using hydrolysis probes to detect exon/exon junctions in transgenes because the operation is simple and economical. However, when the pseudogene is present in the host genome, the PCR-based methods may have a potential risk of detecting false positives. In this study, because processed pseudogenes that exist less frequently in the horse genome may affect PCR-based transgene detection in gene-doping tests, we propose and demonstrate that PCR amplification and sequencing using primers designed on transgene and promotors and/or polyadenylation signal for gene expression are useful for gene-doping detection as an additional confirmatory test to prevent false positives.

Detection of non-targeted transgenes by whole-genome resequencing for gene-doping control.

Gene doping has raised concerns in human and equestrian sports and the horseracing industry. There are two possible types of gene doping in the sports and racing industry: (1) administration of a gene-doping substance to postnatal animals and (2) generation of genetically engineered animals by modifying eggs. In this study, we aimed to identify genetically engineered animals by whole-genome resequencing (WGR) for gene-doping control. Transgenic cell lines, in which the erythropoietin gene (EPO) cDNA form was inserted into the genome of horse fibroblasts, were constructed as a model of genetically modified horse. Genome-wide screening of non-targeted transgenes was performed to find structural variation using DELLY based on split-read and paired-end algorithms and Control-FREEC based on read-depth algorithm. We detected the EPO transgene as an intron deletion in the WGR data by the split-read algorithm of DELLY. In addition, single-nucleotide polymorphisms and insertions/deletions artificially introduced in the EPO transgene were identified by WGR. Therefore, genome-wide screening using WGR can contribute to gene-doping control even if the targets are unknown. This is the first study to detect transgenes as intron deletions for gene-doping detection.

Robustness of Digital PCR and Real-Time PCR in Transgene Detection for Gene-Doping Control.

Gene doping is banned in human sports, horseracing, and equestrian sports. One possible form of gene doping is to administer exogenous genes, called transgenes. Several transgene detection methods based on quantitative PCR have been developed. In this study, we investigated the robustness of digital PCR and real-time PCR in transgene detection using primers and probes that matched (P-true) or incompletely matched (P-false) the template DNA. Fluorescence intensity was significantly reduced when substituted probes were used compared to that using the matched probe in both digital and real-time PCR assays. Digital PCR yielded a similar copy number regardless of the probe (P-true: 1230.7, P-false: 1229.7), whereas real-time PCR revealed a decrease in sensitivity based on Cq values (P-true: 23.5, P-false: 29.7). When substituted primers were used, the detected copy number decreased in the digital PCR assay, and the Cq value in real-time PCR was much higher. Interestingly, digital PCR copy numbers improved by performing PCR at a low annealing temperature, even if a substituted probe was used. Thus, when primer and probe sequences did not completely match the template transgene, digital PCR was relatively robust, but real-time PCR was less sensitive. Although PCR specificity may be reduced, PCR sensitivity can be improved by lowering the annealing temperature. If the target sequence is substituted to escape doping detection, it may be desirable to set the annealing temperature lower and use a more robust method, such as digital PCR, to increase the detection of positive cases, which will also result in fewer false-negative results.

Distribution of the Water-Soluble Astaxanthin Binding Carotenoprotein (AstaP) in Scenedesmaceae.

Photooxidative stress-inducible water-soluble astaxanthin-binding proteins, designated as AstaP, were identified in two Scenedesmaceae strains, Coelastrella astaxanthina Ki-4 and Scenedesmus obtusus Oki-4N; both strains were isolated under high light conditions. These AstaPs are classified as a novel family of carotenoprotein and are useful for providing valuable astaxanthin in water-soluble form; however, the distribution of AstaP orthologs in other microalgae remains unknown. Here, we examined the distribution of AstaP orthologs in the family Scenedesmaceae with two model microalgae, Chlamydomonas reinhardtii and Chlorella variabilis. The expression of AstaP orthologs under photooxidative stress conditions was detected in cell extracts of Scenedesmaceae strains, but not in model algal strains. Aqueous orange proteins produced by Scenedesmaceae strains were shown to bind astaxanthin. The protein from Scenedesmus costatus SAG 46.88 was purified. It was named ScosAstaP and found to bind astaxanthin. The deduced amino acid sequence from a gene encoding ScosAstaP showed 62% identity to Ki-4 AstaP. The expression of the genes encoding AstaP orthologs was shown to be inducible under photooxidative stress conditions; however, the production amounts of AstaP orthologs were estimated to be approximately 5 to 10 times lower than that of Ki-4 and Oki-4N.

Investigation of erythrocyte antigen frequencies in draft horse populations in Japan to assess blood donor suitability.

Erythrocyte alloantigen frequencies of draft horses in Japan were investigated to assess blood donor suitability for transfusion. Here, 148 Japanese draft, 69 Percheron, and 65 Breton horses were blood-typed and subjected to an indirect antiglobulin test. Regarding the major immunogenic factors, the rates of Aa- and Qa-negative horses ranged from 0.35 to 0.49 and from 0.82 to 1.00, respectively. The rate of alloantibody-positive horses ranged from 0.12 to 0.35. Although the prevalence of alloantibodies in these horses was higher than that expected naturally, the rates of Aa- and Qa-negative horses were higher than those of some breeds reported previously. The current draft horse population could provide potential candidates for donors, and the obtained information may contribute to the selection of a safe donor for transfusion.

Design and storage stability of reference materials for microfluidic quantitative PCR-based equine gene doping tests.

One method of gene doping in horseracing is administering of exogenous genetic materials, known as transgenes. Several polymerase chain reaction (PCR)-based methods have been developed for detecting transgenes with high sensitivity and specificity. However, novel designs for reference materials (RMs) and/or positive template controls (PTCs) are necessary for simultaneous analysis of multiple transgene targets. In this study, we designed and developed a novel RM for simultaneously detecting multiple targets via microfluidic quantitative PCR (MFQPCR). Twelve equine genes were selected as targets in this study. A sequence region including primers and probes for quantitative PCR was designed, and a 10 bp sequence was inserted to allow the RM to be distinguished from the original transgene sequences. The sequences of individual detection sites were then connected for 12 genes and cloned into a single plasmid vector. We performed fragment size analysis to distinguish between the PCR products of the original transgene sequence and those of the RM, enabling identification of RM contamination. PTCs diluted to 10,000, 1,000, 100, and 10 copies/µl with horse genomic DNA from RM were stably stored at 4°C for 1 year. As digital PCR enabled absolute quantification, the designed substances can serve as an RM. These findings indicate that the RM design and storage conditions were suitable for gene doping tests using MFQPCR.

Analysis of the Diversity of the AvBD Gene Region in Japanese Quail.

The avian ß-defensin (AvBD) gene region is an important component of the innate immune system, encoding a variety of antimicrobial peptides. The AvBD region forms a multigene cluster in a specific chromosomal region. Comparison of the AvBD region among various birds suggests the presence of defects, duplications, and pseudogenization at many loci. The AvBD region in certain galliform birds, namely chicken, turkey, and bobwhite quail, includes AvBD3, -6, and -7, with the latter exhibiting copy number variants (CNVs) in chickens. DNA for genomic analysis was extracted from the peripheral blood of 99 randomly selected quail (Coturnix japonica) from 6 inbred lines. Nine CjAvBD1 and 8 CjAvBD12 alleles were detected. Ten haplotypes, including three that were strain specific, were found in alleles from the quail AvBD1 (CjAvBD1) and -12 (CjAvBD12) loci. Next-generation sequencing was used to determine the nucleotide sequences of the CjAvBD gene region (56-70 kb) for 7 homozygous diplotypes of these 10 haplotypes. These 7 haplotypes contained between 12 and 16 CjAvBD genes and were composed of 11 common loci: CjAvBD1, -2, -4, -5, -8, -9, -10, -11, -12, -13, and -14, but lacked CjAvBD3 and -7. Furthermore, up to 5 CjAvBD101 (AvBD6 ortholog) CNVs were observed among the 7 haplotypes. In addition, we detected amino acid substitutions causing net charge mutations that could affect antimicrobial activity in CjAvBD4, -13, -14, and -101. These results suggest that the CjAvBD region is unique among the Galliformes and that its diversity results in potential functional variation in innate immunity.

Whole-genome resequencing using genomic DNA extracted from horsehair roots for gene-doping control in horse sports.

Gene doping is prohibited in horseracing and equestrian sports. In previous studies, we developed non-targeted transgene and genome editing detection methods based on whole genome resequencing (WGR) using genomic DNA extracted from whole blood. In this study, we aimed to develop a WGR method using DNA extracts from hair roots. Hair roots are a preferred substrate because their collection is less invasive than blood collection. Hair is also easier to store for long periods of time. Although almost all genomic DNA extracted from hair root samples stored for years at room temperature was degraded, the quality of genomic DNA from samples stored for years at refrigerated temperatures (4-8°C) was maintained. High-molecular-weight genomic DNA was isolated from hair roots using a magnetic silica beads method of extraction, enabling WGR from horsehair root extracts. Nucleotide sequencing results and numbers of single-nucleotide polymorphisms and insertions/deletions concurred with those previously reported for WGR of DNA extracted from whole blood. Therefore, we consider that storing hair samples at refrigerated temperatures prevents degradation of DNA, allowing the detection of gene doping in these samples based on WGR. It is likely this finding will also have a deterrent effect, as it is now possible to test horses with archived samples even if they or their parents are deceased. To our knowledge, this is the first report employing WGR on horsehair roots stored for a long term.

Suppression of substrate inhibition in phenanthrene-degrading Mycobacterium by co-cultivation with a non-degrading Burkholderia strain.

In natural environments contaminated by recalcitrant organic pollutants, efficient biodegradation of such pollutants has been suggested to occur through the cooperation of different bacterial species. A phenanthrene-degrading bacterial consortium, MixEPa4, from polluted soil was previously shown to include a phenanthrene-degrading strain, Mycobacterium sp. EPa45, and a non-polycyclic aromatic hydrocarbon (PAH)-degrading strain, Burkholderia sp. Bcrs1W. In this study, we show that addition of phenanthrene to rich liquid medium resulted in the transient growth arrest of EPa45 during its degradation of phenanthrene. RNA-sequencing analysis of the growth-arrested cells showed the phenanthrene-dependent induction of genes that were predicted to be involved in the catabolism of this compound, and many other cell systems, such as a ferric iron-uptake, were up-regulated, implying iron deficiency of the cells. This negative effect of phenanthrene became much more apparent when using phenanthrene-containing minimal agar medium; colony formation of EPa45 on such agar was significantly inhibited in the presence of phenanthrene and its intermediate degradation products. However, growth inhibition was suppressed by the co-residence of viable Bcrs1W cells. Various Gram-negative bacterial strains, including the three other strains from MixEPa4, also exhibited varying degrees of suppression of the growth inhibition effect on EPa45, strongly suggesting that this effect is not strain-specific. Growth inhibition of EPa45 was also observed by other PAHs, biphenyl and naphthalene, and these two compounds and phenanthrene also inhibited the growth of another mycobacterial strain, M. vanbaalenii PYR-1, that can use them as carbon sources. These phenomena of growth inhibition were also suppressed by Bcrs1W. Our findings suggest that, in natural environments, various non-PAH-degrading bacterial strains play potentially important roles in the facilitation of PAH degradation by the co-residing mycobacteria.

Isolation of Salmonella enterica serovar Agona strains and their similarities to strains derived from a clone caused a serovar shift in broilers.

Salmonella enterica serovar Agona strains isolated from human cases were compared to strains that were derived from a clone caused a serovar shift in broilers. Pulsed field gel electrophoresis (PFGE) analysis with XbaI or BlnI digestion showed that three of seven strains from human case strains and most of the 81 strains from broilers were clustered in single complex in a minimum spanning tree (MST) reconstructed from the PFGE data. All the strains from human cases and 22 randomly selected strains from broilers were also analyzed by whole genome sequencing (WGS). Analysis of single nucleotide polymorphism (SNP) in the S. Agona core genes showed that four strains from human cases and all the strains from broilers were clustered in a maximum likelihood phylogenetic tree (ML tree) and an MST. These results indicated that the strains derived from the clone caused the serovar shift had already spread to humans. PFGE analysis with XbaI showed that four strains from broilers did not cluster with the other strains in an MST, though all those strains clustered in an ML tree and an MST reconstructed from SNP data. Moreover, three strains from broilers did not cluster in an MST reconstructed from PFGE with BlnI digestion, though those strains clustered in an ML tree and an MST reconstructed from SNP data. Therefore, it was suggested that S. Agona strains derived from a particular clone could not be traced by PFGE analysis but can be investigated by WGS analysis.

Isolation of colonization-defective Escherichia coli mutants reveals critical requirement for fatty acids in bacterial colony formation.

Most bacterial cells in nature exhibit extremely low colony-forming activity, despite showing various signs of viability, impeding the isolation and utilization of many bacterial resources. However, the general causes responsible for this state of low colony formation are largely unknown. Because liquid cultivation typically yields more bacterial cell cultures than traditional solid cultivation, we hypothesized that colony formation requires one or more specific gene functions that are dispensable or less important for growth in liquid media. To verify our hypothesis and reveal the genetic background limiting colony formation among bacteria in nature, we isolated Escherichia coli mutants that had decreased frequencies of colony formation but could grow in liquid medium from a temperature-sensitive mutant collection. Mutations were identified in fabB, which is essential for the synthesis of long unsaturated fatty acids. We then constructed a fabB deletion mutant in a wild-type background. Detailed behavioural analysis of the mutant revealed that under fatty acid-limited conditions, colony formation on solid media was more sensitively and seriously impaired than growth in liquid media. Furthermore, growth under partial inhibition of fatty acid synthesis with cerulenin or triclosan brought about similar phenotypes, not only in E. coli but also in Bacillus subtilis and Corynebacterium glutamicum. These results indicate that fatty acids have a critical importance in colony formation and that depletion of fatty acids in the environment partly accounts for the low frequency of bacterial colony formation.

Lactobacillus metriopterae sp. nov., a novel lactic acid bacterium isolated from the gut of grasshopper Metrioptera engelhardti.

A strain (Hime 5-1T) of lactic acid bacterium was isolated from the gut of the grasshopper Metrioptera engelhardti from a mountainous area of Nagano Prefecture, Japan. Strain Hime 5-1T had a low 16S rRNA gene sequence similarity to known lactic acid bacteria, with the closest recognized relatives being Lactobacillus tucceti (96.7 %), Lactobacillus furfuricola (96.5 %), Lactobacillus versmoldensis (96.3 %) and Lactobacillus nodensis (96.1 %). Comparative analyses of the rpoA and pheS gene sequences indicated that Hime 5-1T is not closely related to other Lactobacillus species. Strain Hime 5-1T is a Gram-stain-positive, catalase-negative and homofermentative bacterium with yellowish colonies, which contrasts with the whitish colonies of its closest recognized relatives. Based on phenotypic and genotypic properties, we conclude that the isolated bacterium represents a novel species of the genus Lactobacillus, for which the name Lactobacillus metriopterae sp. nov. is proposed. The type strain is Hime 5-1T (=JCM 31635T=DSM 103730T). 16S rRNA gene based high-throughput sequencing revealed that L. metriopterae is the dominant microbiota in the gut of Metrioptera engelhardti.

Transcriptomic characterization of gall tissue of Japanese elm tree (Ulmus davidiana var. japonica) induced by the aphid Tetraneura nigriabdominalis.

Insect galls are abnormal plant tissues induced by parasitic insect(s) for use as their habitat. In previous work, we suggested that gall tissues induced by the aphid Tetraneura nigriabdominalis on Japanese elm trees are less responsive than leaf tissues to jasmonic acid (JA), which is involved in the production of volatile organic compounds as a typical defensive reaction of plants against attack by insect pests. A comprehensive analysis of gene expression by RNA sequencing indicated that the number of JA responsive genes was markedly lower in gall tissues than in leaf tissues. This suggests that gall tissues are mostly defective in JA signaling, although JA signaling is not entirely compromised in gall tissue. Gene ontology analysis sheds light on some stress-related unigenes with higher expression levels in gall tissues, suggesting that host plants sense aphids as a biotic stress but are defective in the JA-mediated defense response in gall tissues.