Coding complete genome of a strain of ovine herpesvirus-2 associated with a clinical case of malignant catarrhal fever in a domestic lamb

Ovine herpesvirus-2 causes sheep-associated malignant catarrhal fever, a fatal disease of ruminants and pigs. The virus is carried by sheep, and infection is typically subclinical. Here, we report the coding complete genome sequence of a strain of OvHV-2 obtained from a clinically affected domestic lamb.

A U.S. Isolate of Theileria orientalis Ikeda Is Not Transstadially Transmitted to Cattle by Rhipicephalus microplus

Theileria orientalis Ikeda has caused an epidemic of bovine anemia and abortion across several U.S. states. This apicomplexan hemoparasite is transmitted by Haemaphysalis longicornis ticks; however, it is unknown if other North American ticks are competent vectors. Since the disease movement is largely determined by the host tick range(s), the prediction of the T. orientalis spread among U.S. cattle populations requires determination of additional competent tick vectors. Although Rhipicephalus microplus has mostly been eradicated from the U.S., outbreaks in populations occur frequently, and the U.S. remains at risk for reintroduction. Since R. microplus is a vector of Theileria equi and T. orientalis DNA has been detected in R. microplus, the goal of this study was to determine whether R. microplus is a competent vector of T. orientalis. Larval R. microplus were applied to a splenectomized, T. orientalis Ikeda-infected calf for parasite acquisition, removed as molted adults, and applied to two T. orientalis naïve, splenectomized calves for transmission. After 60 days, the naïve calves remained negative for T. orientalis by PCR and cytology. Additionally, T. orientalis was not detected in the salivary glands or larval progeny of acquisition-fed adults. These data suggest that R. microplus is not a competent vector of the U.S. T. orientalis Ikeda isolate.

Development of a real-time PCR assay for detection and differentiation of Mycoplasma ovipneumoniae and a novel respiratory-associated Mycoplasma species in domestic sheep and goats

A novel respiratory-associated Mycoplasma species (M. sp. nov.) of unknown clinical significance was recently identified that causes false positive results with multiple published PCR methods reported to specifically detect Mycoplasma ovipneumonaie, a well-known respiratory pathogen in small ruminants. This necessitates our objective to develop a real-time PCR (qPCR) assay for improved specificity and sensitivity, and more rapid detection and differentiation of M. ovipneumoniae and the M. sp. nov. in domestic sheep (DS) and domestic goat (DG) samples, as compared to a conventional PCR and sequencing (cPCR-seq) assay. Primers and probes were designed based on available M. ovipneumoniae 16S rRNA gene sequences in the GenBank database, and partial 16S rRNA gene sequences provided by the United States Department of Agriculture, Agricultural Research Service (USDA-ARS) for M. ovipneumoniae and M. sp. nov. USDA-ARS provided DS (n = 153) and DG (n = 194) nasal swab nucleic acid that previously tested positive for either M. ovipneumoniae (n = 117) or M. sp. nov. (n = 138), or negative for both targets (n = 92) by cPCR-seq. A host 18S rRNA gene was included as an internal control to monitor for the failure of nucleic acid extraction and possible PCR inhibition. For samples positive by cPCR-seq, qPCR agreement was 88.0% (103/117; κ = 0.81) and 89.9% (124/138; κ = 0.84) for M. ovipneumoniae and M. sp. nov., respectively; 12 of 255 (4.7%) cPCR-seq positive samples were qPCR positive for both targets. Of samples negative by cPCR for both mycoplasmas, qPCR detected M. ovipneumoniae and M. sp. nov. in 6.5% (6/92) and 4.3% (4/92), respectively. Samples with discordant results between the cPCR and sequencing assay and the new qPCR were analyzed by target sequencing; successfully sequenced samples had identity matches that confirmed the qPCR result. The increased target specificity of this qPCR is predicted to increase testing accuracy as compared to other published assays.

A U.S. isolate of Theileria orientalis, Ikeda genotype, is transmitted to cattle by the invasive Asian longhorned tick, Haemaphysalis longicornis

Background

Theileria orientalis is a tick-borne hemoparasite that causes anemia, ill thrift, and death in cattle globally. The Ikeda strain of T.orientalis is more virulent than other strains, leading to severe clinical signs and death of up to 5% of affected animals. Within the Asia–Pacific region, where it affects 25% of Australian cattle, T.orientalis Ikeda has a significant economic impact on the cattle industry. In 2017, T.orientalis Ikeda was detected in a cattle herd in Albermarle County, Virginia, United States. Months earlier, the U.S. was alerted to the invasion of the Asian longhorned tick, Haemaphysalis longicornis, throughout the eastern U.S. Abundant H.longicornis ticks were identified on cattle in the T.orientalis-affected herd in VA, and a subset of ticks from the environment were PCR-positive for T.orientalis Ikeda. A strain of T.orientalis from a previous U.S. outbreak was not transmissible by H.longicornis; however, H.longicornis is the primary tick vector of T.orientalis Ikeda in other regions of the world. Thus, the objective of this study was to determine whether invasive H.longicornis ticks in the U.S. are competent vectors of T.orientalis Ikeda.

Methods

Nymphal H.longicornis ticks were fed on a splenectomized calf infected with the VA-U.S.-T.orientalis Ikeda strain. After molting, a subset of adult ticks from this cohort were dissected, and salivary glands assayed for T.orientalis Ikeda via qPCR. The remaining adult ticks from the group were allowed to feed on three calves. Calves were subsequently monitored for T.orientalis Ikeda infection via blood smear cytology and PCR.

Results

After acquisition feeding on a VA-U.S.-T.orientalis Ikeda-infected calf as nymphs, a subset of molted adult tick salivary glands tested positive by qPCR for T.orientalis Ikeda. Adult ticks from the same cohort successfully transmitted T.orientalis Ikeda to 3/3 naïve calves, each of which developed parasitemia reaching 0.4–0.9%.

Conclusions

Our findings demonstrate that U.S. H.longicornis ticks are competent vectors of the VA-U.S.-T.orientalis Ikeda strain. This data provides important information for the U.S. cattle industry regarding the potential spread of this parasite and the necessity of enhanced surveillance and control measures.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13071-021-04659-9.

Differential pulmonary immunopathology of domestic sheep (Ovis aries) and bighorn sheep (Ovis canadensis) with Mycoplasma ovipneumoniae infection: A retrospective study

Mycoplasma ovipneumoniae is a respiratory pathogen that impacts domestic sheep (Ovis aries; DS) and bighorn sheep (Ovis canadensis; BHS). BHS are reported to be more susceptible than DS to developing polymicrobial pneumonia associated with M. ovipneumoniae infection. Using formalin-fixed paraffin-embedded tissues, we performed a retrospective study investigating the pulmonary immune response of DS and BHS to M. ovipneumoniae infection. M. ovipneumoniae infected DS exhibited a more robust and well-organized BALT formation as compared to BHS. Digital analysis of immunohistochemical chromogen deposition in lung tissue was used to quantitate T cell marker CD3, B cell markers CD20 and CD79a, macrophage markers CD163 and Iba1, and cytokine IL-17. A significant interaction of species and infection status was identified for CD3, CD163, and IL-17. BHS had a greater increase in bronchiolar CD3 and bronchiolar and alveolar CD163 with infection, as compared to DS. BHS had an increase in bronchiolar associated lymph tissue (BALT) and alveolar IL-17 with infection, while these remained similar in DS regardless of infection status. IL-17 in respiratory epithelium of bronchi and bronchioles comparatively decreased in DS and increased in BHS with infection. These data begin to define the interspecies differential immune response to pulmonary M. ovipneumoniae infection in DS and BHS and provide the first investigations of respiratory epithelium-associated IL-17 in ovine.

Genome-Wide Histone Modifications and CTCF Enrichment Predict Gene Expression in Sheep Macrophages

Alveolar macrophages function in innate and adaptive immunity, wound healing, and homeostasis in the lungs dependent on tissue-specific gene expression under epigenetic regulation. The functional diversity of tissue resident macrophages, despite their common myeloid lineage, highlights the need to study tissue-specific regulatory elements that control gene expression. Increasing evidence supports the hypothesis that subtle genetic changes alter sheep macrophage response to important production pathogens and zoonoses, for example, viruses like small ruminant lentiviruses and bacteria like Coxiella burnetii. Annotation of transcriptional regulatory elements will aid researchers in identifying genetic mutations of immunological consequence. Here we report the first genome-wide survey of regulatory elements in any sheep immune cell, utilizing alveolar macrophages. We assayed histone modifications and CTCF enrichment by chromatin immunoprecipitation with deep sequencing (ChIP-seq) in two sheep to determine cis-regulatory DNA elements and chromatin domain boundaries that control immunity-related gene expression. Histone modifications included H3K4me3 (denoting active promoters), H3K27ac (active enhancers), H3K4me1 (primed and distal enhancers), and H3K27me3 (broad silencers). In total, we identified 248,674 reproducible regulatory elements, which allowed assignment of putative biological function in macrophages to 12% of the sheep genome. Data exceeded the FAANG and ENCODE standards of 20 million and 45 million useable fragments for narrow and broad marks, respectively. Active elements showed consensus with RNA-seq data and were predictive of gene expression in alveolar macrophages from the publicly available Sheep Gene Expression Atlas. Silencer elements were not enriched for expressed genes, but rather for repressed developmental genes. CTCF enrichment enabled identification of 11,000 chromatin domains with mean size of 258 kb. To our knowledge, this is the first report to use immunoprecipitated CTCF to determine putative topological domains in sheep immune cells. Furthermore, these data will empower phenotype-associated mutation discovery since most causal variants are within regulatory elements.

Transcriptome dataset of Babesia bovis life stages within vertebrate and invertebrate hosts

Babesia bovis is a hemoprotozoan parasite of cattle that has a complex life cycle within vertebrate and invertebrate hosts. In the mammalian host, B. bovis undergoes asexual reproduction while in the tick midgut, gametes are induced, fuse, and form zygotes. The zygote infects tick gut epithelial cells and transform into kinetes that are released into the hemolymph and invade other tick tissues such as the ovaries, resulting in transovarial transmission to tick offspring. To compare gene regulation between different B. bovis life stages, we collected parasites infecting bovine erythrocytes and tick hemolymph. Total RNA samples were isolated, and multiplexed libraries sequenced using paired-end 100 cycle reads of a HiSeq 2500. The data was normalized using the TMM method and analysed for significant differential expression using the generalized linear model likelihood ratio test (GLM LRT) in edgeR. To validate our datasets, ten genes were selected using NormFinder. Genes that had no significant fold change between the blood and tick stages in the RNA-Seq datasets were tested by quantitative PCR to determine their suitability as “housekeeping” genes. The normalized RNA-Seq data revealed genes upregulated during infection of the mammalian host or tick vector and six upregulated genes were validated by quantitative PCR. These datasets can help identify useful targets for controlling bovine babesiosis.

Comparative analysis of gene expression between Babesia bovis blood stages and kinetes allowed by improved genome annotation

Throughout their life cycle, Babesia parasites alternate between a mammalian host, where they cause babesiosis, and the tick vector. Transition between hosts results in distinct environmental signals that influence patterns of gene expression, consistent with the morphological and functional changes operating in the parasites during their life stages. In addition, comparing differential patterns of gene expression among mammalian and tick parasite stages can provide clues for developing improved methods of control. Hereby, we upgraded the genome assembly of Babesia bovis, a bovine hemoparasite, closing a 139 kbp gap, and used RNA-Seq datasets derived from mammalian blood and tick kinete stages to update the genome annotation. Of the originally annotated genes, 1,254 required structural changes, and 326 new genes were identified, leading to a different predicted proteome compared to the original annotation. Next, the RNA-Seq data was used to identify B. bovis genes that were differentially expressed in the vertebrate and arthropod hosts. In blood stages, 28% of the genes were upregulated up to 300 fold, whereas 26% of the genes in kinetes, a tick stage, were upregulated up to >19,000 fold. We thus discovered differentially expressed genes that may play key biological roles and serve as suitable targets for the development of vaccines to control bovine babesiosis.

Mycoplasma ovis infection in domestic sheep (Ovis aries) in the United States: Prevalence, distribution, associated risk factors, and associated outcomes

Mycoplasma ovis is a hemotropic bacterium reported to infect sheep, goats, and deer species. Infection in these species can result in anemia, jaundice, and ill-thrift. Although of worldwide distribution, only rare reports of this bacterium in the United States exist. The objectives of this retrospective study were to identify the prevalence and distribution of M. ovis, and identify associated demographic and management risk factors, and reproductive and production outcomes associated with infection on domestic sheep (Ovis aries) operations in the United States. As part of the United States Department of Agriculture (USDA), Animal Plant Health Inspection Service, Veterinary Services' National Animal Health Monitoring System (NAHMS) Sheep 2001 and 2011 studies, blood was collected and sera banked from 21,369 ewes in 2001 and 13,128 ewes in 2011. Participating premises were located in 22 states across the United States for each sample year. In 2015 the USDA, Agricultural Research Service, Animal Disease Research Unit received aliquots of these sera, and DNA was extracted and analyzed by PCR for the presence of M. ovis genomic DNA. Flock presence and mean within-flock prevalence of M. ovis were 73.3% and 23.2%, respectively. Model selection using Mallow's C_p Criterion was used to determine which variables significantly affected flock presence and within-flock prevalence. The final flock presence model included flock size, year of blood collection, region, and vaccine administration. The final within-flock prevalence model included year of blood collection, interaction between flock size and region, and interaction between reported abortions and grazing with sheep from other operations. Medium and large operations had a higher flock presence and within-flock prevalence. Flock presence was higher in operations that administered any vaccines. Operations that reported any abortions and grazed with sheep from other operations had a higher within-flock prevalence.

Mycoplasma ovipneumoniae in Wildlife Species beyond Subfamily Caprinae

Elucidating the emergence of Mycoplasma ovipneumoniae-associated respiratory disease in ruminants requires identification of the pathogen host range. This bacterium was thought to be host restricted to subfamily Caprinae, but we describe its identification in healthy moose, caribou, and mule deer and diseased mule and white-tailed deer, all species in subfamily Capreolinae.

The Babesia bovis hap2 gene is not required for blood stage replication, but expressed upon in vitro sexual stage induction

Babesia bovis, is a tick borne apicomplexan parasite responsible for important cattle losses globally. Babesia parasites have a complex life cycle including asexual replication in the mammalian host and sexual reproduction in the tick vector. Novel control strategies aimed at limiting transmission of the parasite are needed, but transmission blocking vaccine candidates remain undefined. Expression of HAP2 has been recognized as critical for the fertilization of parasites in the Babesia-related Plasmodium, and is a leading candidate for a transmission blocking vaccine against malaria. Hereby we identified the B. bovis hap2 gene and demonstrated that it is widely conserved and differentially transcribed during development within the tick midgut, but not by blood stage parasites. The hap2 gene was disrupted by transfecting B. bovis with a plasmid containing the flanking regions of the hap2 gene and the GPF-BSD gene under the control of the ef-1α-B promoter. Comparison of in vitro growth between a hap2-KO B. bovis clonal line and its parental wild type strain showed that HAP2 is not required for the development of B. bovis in erythrocytes. However, xanthurenic acid-in vitro induction experiments of sexual stages of parasites recovered after tick transmission resulted in surface expression of HAP2 exclusively in sexual stage induced parasites. In addition, hap2-KO parasites were not able to develop such sexual stages as defined both by morphology and by expression of the B. bovis sexual marker genes 6-Cys A and B. Together, the data strongly suggests that tick midgut stage differential expression of hap2 is associated with the development of B. bovis sexual forms. Overall these studies are consistent with a role of HAP2 in tick stages of the parasite and suggest that HAP2 is a potential candidate for a transmission blocking vaccine against bovine babesiosis.

Babesia bovis, is a tick borne apicomplexan parasite responsible for important cattle losses globally. Babesia parasites have a complex life cycle including asexual replication in the mammalian host and sexual reproduction in the tick vector. Novel control strategies aimed at limiting transmission of the parasite are needed, but transmission blocking vaccine candidates remain undefined. In this study we analyze the conservation and role of the hap2 gene in the erythrocyte stage of the life cycle of the parasite and found that expression of the gene is not required for the development of the parasite in erythrocytic stages, using a hap2 mutated parasite line. In addition, we developed an in vitro system for the induction of sexual forms of B. bovis and found expression of the hap2 gene and surface localization of the protein. However, hap2-KO parasites are unable to develop sexual stages. We concluded that HAP2 is a leading candidate for a transmission blocking vaccine against bovine babesiosis due of the high level of conservation, surface exposure, and specific expression in tick stage and in in vitro induced sexual stages parasites.

Geno- and phenotypic characteristics of a transfected Babesia bovis 6-Cys-E knockout clonal line

BACKGROUND

Babesia bovis is an intra-erythrocytic tick-transmitted apicomplexan protozoan parasite. It has a complex lifestyle including asexual replication in the mammalian host and sexual replication occurring in the midgut of host tick vector, typically, Rhipicephalus microplus. Previous evidence showed that certain B. bovis genes, including members of 6-Cys gene family, are differentially expressed during tick and mammalian stages of the parasite's life cycle. Moreover, the 6-Cys E gene is differentially expressed in the T3Bo strain of B. bovis tick stages, and anti 6-Cys E antibodies were shown to be able to inhibit in vitro growth of the phenotypically distinct B. bovis Mo7clonal line.

METHODS

In this study, the 6-Cys E gene of B. bovis T3Bo strain was disrupted by transfection using a plasmid containing 6-Cys gene E 5' and 3' regions to guide homologous recombination, and the egfp-bsd fusion gene under control of a ef-1α promoter, yielding a B. bovis clonal line designated 6-Cys EKO-cln. Full genome sequencing of 6-Cys EKO-cln parasites was performed and in vitro inhibition assays using anti 6-Cys E antibodies.

RESULTS

Full genome sequencing of 6-Cys EKO-cln B. bovis demonstrated single insertion of egfp-bsd gene that disrupts the integrity of 6-Cys gene E. Undistinguishable growth rate of 6-Cys EKO-cln line compared to wild-type 6-Cys E intact T3Bo B. bovis strain in in vitro cultures indicates that expression of gene 6-Cys E is not essential for blood stage replication in this strain. In vitro inhibition assays confirmed the ability of anti-6 Cys E antibodies to inhibit the growth of the wild-type Mo7 and T3Bo B. bovis parasites, but no significant inhibition was found for 6-Cys EKO-cln line parasites.

CONCLUSIONS

Overall, the data suggest that the anti-6 Cys E antibody neutralising effect on the wild type strains is likely due to mechanical hindrance, or cross-reactivity, rather than due to functional requirements of 6-Cys gene E product for survival and development of the erythrocyte stages. Further investigation is underway to determine if the 6-Cys E protein is required for replication and sexual stage development of B. bovis during tick stages.

Expression of 6-Cys Gene Superfamily Defines Babesia bovis Sexual Stage Development within Rhipicephalus microplus

Babesia bovis, an intra-erythrocytic tick-borne apicomplexan protozoan, is one of the causative agents of bovine babesiosis. Its life cycle includes sexual reproduction within cattle fever ticks, Rhipicephalus spp. Six B. bovis 6-Cys gene superfamily members were previously identified (A, B, C, D, E, F) where their orthologues in Plasmodium parasite have been shown to encode for proteins required for the development of sexual stages. The current study identified four additional 6-Cys genes (G, H, I, J) in the B. bovis genome. These four genes are described in the context of the complete ten 6-Cys gene superfamily. The proteins expressed by this gene family are predicted to be secreted or surface membrane directed. Genetic analysis comparing the 6-Cys superfamily among five distinct B. bovis strains shows limited sequence variation. Additionally, A, B, E, H, I and J genes were transcribed in B. bovis infected tick midgut while genes A, B and E were also transcribed in the subsequent B. bovis kinete stage. Transcription of gene C was found exclusively in the kinete. In contrast, transcription of genes D, F and G in either B. bovis infected midguts or kinetes was not detected. None of the 6-Cys transcripts were detected in B. bovis blood stages. Subsequent protein analysis of 6-Cys A and B is concordant with their transcript profile. The collective data indicate as in Plasmodium parasite, certain B. bovis 6-Cys family members are uniquely expressed during sexual stages and therefore, they are likely required for parasite reproduction. Within B. bovis specifically, proteins encoded by 6-Cys genes A and B are markers for sexual stages and candidate antigens for developing novel vaccines able to interfere with the development of B. bovis within the tick vector.

Integration of a transfected gene into the genome of Babesia bovis occurs by legitimate homologous recombination mechanisms

This study examines the patterns of gene integration of gfp-bsd upon stable transfection into the T3Bo strain of Babesia bovis using a plasmid designed to integrate homologous sequences of the parasite's two identical ef-1α A and B genes. While the transfected BboTf-149-6 cell line displayed two distinct patterns of gene integration, clonal lines derived from this strain by cell sorting contained only single gfp-bsd insertions. Whole genome sequencing of two selected clonal lines, E9 and C6, indicated two distinct patterns of gfp-bsd insertion occurring by legitimate homologous recombination mechanisms: one into the expected ef-1α orf B, and another into the ef-1α B promoter. The data suggest that expression of the ef-1α orf B is not required for development of B. bovis in cultured erythrocyte stages. Use of legitimate homologous recombination mechanisms in transfected B. bovis supports the future use of transfection methods for developing efficient gene function assignment experiments using gene knockout techniques.

Malignant catarrhal fever in American bison (Bison bison) experimentally infected with alcelaphine herpesvirus 2

Malignant catarrhal fever (MCF), due to ovine herpesvirus 2 (OvHV-2), causes appreciable death loss in ranched bison (Bison bison) throughout North America. No vaccine exists to protect animals from disease. Since OvHV-2 has not been propagated in vitro, one strategy to develop a modified live vaccine is to use a closely related, non-pathogenic member of the malignant catarrhal fever virus family as a vector expressing potentially protective OvHV-2 epitopes. To date, no controlled experimental challenge studies with alcelaphine herpesvirus 2 (AlHV-2) derived from topi (Damaliscus lunatus jimela) have been reported The unique or light DNA segment of the AlHV-2 genome was sequenced and annotated and the virus was tested for its ability to infect and induce disease in American bison. Yearling bison were inoculated intranasally (n=4) or intramuscularly (n=3) with 2 × 10(-4.7) TCID50 of AlHV-2, and monitored for infection and the development of disease. Six inoculated bison became infected with AlHV-2. Two of the six animals developed clinical signs and had gross and histological lesions consistent with terminal MCF, which differed in distribution from those in bison with MCF due to OvHV-2. One other animal developed minor clinical signs and had gross and histological pulmonary lesions consistent with early (pre-clinical) stages of MCF. Unmodified low cell culture passage AlHV-2 derived from topi is an unsuitable vaccine vector for the prevention of MCF. However, the annotated genome might be useful in identifying genes which could be deleted to potentially attenuate the virus for bison.

A divergent Artiodactyl MYADM-like repeat is associated with erythrocyte traits and weight of lamb weaned in domestic sheep

A genome-wide association study (GWAS) was performed to investigate seven red blood cell (RBC) phenotypes in over 500 domestic sheep (Ovis aries) from three breeds (Columbia, Polypay, and Rambouillet). A single nucleotide polymorphism (SNP) showed genome-wide significant association with increased mean corpuscular hemoglobin concentration (MCHC, P = 6.2×10⁻¹⁴) and genome-wide suggestive association with decreased mean corpuscular volume (MCV, P = 2.5×10⁻⁶). The ovine HapMap project found the same genomic region and the same peak SNP has been under extreme historical selective pressure, demonstrating the importance of this region for survival, reproduction, and/or artificially selected traits. We observed a large (>50 kb) variant haplotype sequence containing a full-length divergent artiodactyl MYADM-like repeat in strong linkage disequilibrium with the associated SNP. MYADM gene family members play roles in membrane organization and formation in myeloid cells. However, to our knowledge, no member of the MYADM gene family has been identified in development of morphologically variant RBCs. The specific RBC differences may be indicative of alterations in morphology. Additionally, erythrocytes with altered morphological structure often exhibit increased structural fragility, leading to increased RBC turnover and energy expenditure. The divergent artiodactyl MYADM-like repeat was also associated with increased ewe lifetime kilograms of lamb weaned (P = 2×10⁻⁴). This suggests selection for normal RBCs might increase lamb weights, although further validation is required before implementation in marker-assisted selection. These results provide clues to explain the strong selection on the artiodactyl MYADM-like repeat locus in sheep, and suggest MYADM family members may be important for RBC morphology in other mammals.

Comparative genomic analysis and phylogenetic position of Theileria equi

Background

Transmission of arthropod-borne apicomplexan parasites that cause disease and result in death or persistent infection represents a major challenge to global human and animal health. First described in 1901 as Piroplasma equi, this re-emergent apicomplexan parasite was renamed Babesia equi and subsequently Theileria equi, reflecting an uncertain taxonomy. Understanding mechanisms by which apicomplexan parasites evade immune or chemotherapeutic elimination is required for development of effective vaccines or chemotherapeutics. The continued risk of transmission of T. equi from clinically silent, persistently infected equids impedes the goal of returning the U. S. to non-endemic status. Therefore comparative genomic analysis of T. equi was undertaken to: 1) identify genes contributing to immune evasion and persistence in equid hosts, 2) identify genes involved in PBMC infection biology and 3) define the phylogenetic position of T. equi relative to sequenced apicomplexan parasites.

Results

The known immunodominant proteins, EMA1, 2 and 3 were discovered to belong to a ten member gene family with a mean amino acid identity, in pairwise comparisons, of 39%. Importantly, the amino acid diversity of EMAs is distributed throughout the length of the proteins. Eight of the EMA genes were simultaneously transcribed. As the agents that cause bovine theileriosis infect and transform host cell PBMCs, we confirmed that T. equi infects equine PBMCs, however, there is no evidence of host cell transformation. Indeed, a number of genes identified as potential manipulators of the host cell phenotype are absent from the T. equi genome. Comparative genomic analysis of T. equi revealed the phylogenetic positioning relative to seven apicomplexan parasites using deduced amino acid sequences from 150 genes placed it as a sister taxon to Theileria spp.

Conclusions

The EMA family does not fit the paradigm for classical antigenic variation, and we propose a novel model describing the role of the EMA family in persistence. T. equi has lost the putative genes for host cell transformation, or the genes were acquired by T. parva and T. annulata after divergence from T. equi. Our analysis identified 50 genes that will be useful for definitive phylogenetic classification of T. equi and closely related organisms.

Non-maternal transmission is the major mode of ovine lentivirus transmission in a ewe flock: a molecular epidemiology study

Transmission of ovine progressive pneumonia virus (OPPV), a lentivirus of sheep, occurs through both maternal and non-maternal means. Currently, the contribution of each route to the overall flock OPPV prevalence is poorly understood since previous serological epidemiologic studies lacked the ability to accurately track routes of transmission within an infected flock. In this study, the amount of maternal OPP transmission was assessed in a naturally infected ewe flock by applying molecular analyses to proviral sequences derived from peripheral blood leukocytes of OPP positive dam-daughter pairs (N=40). Both proviral envelope (env) and long terminal repeat (LTR) sequences, separately and combined, were utilized in the following 2 sequence analysis methods: phylogenetic analysis and pairwise distance calculations. True maternal transmission events were defined as agreement in 2 out of the 2 sequence analysis methods. Using this criterion, proviral env sequences resulted in a 14.3% maternal transmission frequency, and proviral LTR sequences resulted in a 10% maternal transmission frequency. Both proportions of maternal transmission varied significantly from equality (P<0.0001). This indicates that the remaining 85.7-90% of daughters are infected via non-maternal transmission. This is also the first study to calculate the OPP proviral rate of change for the env gene and LTR promoter. Accurately defining the routes of OPPV transmission provides critical epidemiological data supporting management intended to reduce flock transmission and viral dose.

Conservation in the face of diversity: multistrain analysis of an intracellular bacterium

BACKGROUND

With the recent completion of numerous sequenced bacterial genomes, notable advances have been made in understanding the level of conservation between various species. However, relatively little is known about the genomic diversity among strains. We determined the complete genome sequence of the Florida strain of Anaplasma marginale, and near complete (>96%) sequences for an additional three strains, for comparative analysis with the previously fully sequenced St. Maries strain genome.

RESULTS

These comparisons revealed that A. marginale has a closed-core genome with few highly plastic regions, which include the msp2 and msp3 genes, as well as the aaap locus. Comparison of the Florida and St. Maries genome sequences found that SNPs comprise 0.8% of the longer Florida genome, with 33.5% of the total SNPs between all five strains present in at least two strains and 3.0% of SNPs present in all strains except Florida. Comparison of genomes from three strains of Mycobacterium tuberculosis, Bacillus anthracis, and Nessieria meningiditis, as well as four Chlamydophila pneumoniae strains found that 98.8%-100% of SNPs are unique to each strain, suggesting A. marginale, with 76.0%, has an intermediate level of strain-specific SNPs. Comparison of genomes from other organisms revealed variation in diversity that did not segregate with the environmental niche the bacterium occupies, ranging from 0.00% to 8.00% of the larger pairwise-compared genome.

CONCLUSION

Analysis of multiple A. marginale strains suggests intracellular bacteria have more variable SNP retention rates than previously reported, and may have closed-core genomes in response to the host organism environment and/or reductive evolution.

Genomic analysis of Ovis aries (Ovar) MHC class IIa loci

Determining the genomic organization of the Ovis aries (Ovar) major histocompatibility complex class IIa region is essential for future functional studies related to antigen presentation. In this study, a bacterial artificial chromosome (BAC) library of genomic DNA from peripheral blood leukocytes (PBL) of a Rambouillet ram was constructed, and BAC clone consisting of the major histocompatibility complex (MHC) class II DQB2, DQA2, DQB1, DQA1, and DRB1 loci was identified and completely sequenced. The BAC clone consists of 160,889 bp of finished sequence with the loci arranged in the following order: DQB2, DQA2, DQB1, DQA1, and DRB1 with 14.3, 25, 6.6, and 40.9 Kb spanning between the loci, respectively. All five of these loci were transcribed in the animal used to generate the MHC class II BAC clone. Repeat or retrotransposable elements along with MHC class II cis promoter elements consisting of S, X, and Y boxes were identified in the sequence. In addition, 16 non-coding conserved sequences amongst primates, carnivores, and ruminants were identified (p < 0.001). These conserved sequences include binding sites for transcription factors with known roles in immune cells, and they provide a basis for further functional investigation of the genes in this region. This is the first ruminant finished sequence of the DQB2-DRB1 region, and this sequence information will aid in whole genome and transcriptome analyses of MHC class II.

Comparison of ovine herpesvirus 2 genomes isolated from domestic sheep (Ovis aries) and a clinically affected cow (Bos bovis)

The rhadinovirus Ovine herpesvirus 2 (OvHV-2) is the causative agent of sheep-associated malignant catarrhal fever. OvHV-2 primarily affects ruminants and has a worldwide distribution. In this study, a composite sequence of OvHV-2 genomic DNA isolated from nasal secretions of sheep experiencing virus-shedding episodes was determined and compared with the sequence of OvHV-2 DNA isolated from a lymphoblastoid cell line derived from a clinically affected cow. The study confirmed the OvHV-2 sequence information determined for the cell line-isolated DNA and showed no apparently significant changes in the OvHV-2 genome during passage through a clinically susceptible species with subsequent maintenance in vitro. Amino acid identity between the predicted open reading frames (ORFs) of the two genomes was 94-100%, except for ORF73, which had an identity of 83%. Polymorphism in ORF73 was due primarily to variability in the G/E-rich repetitive central region of the ORF.

A novel 78-kDa fatty acyl-CoA synthetase (ACS1) of Babesia bovis stimulates memory CD4+ T lymphocyte responses in B. bovis-immune cattle

Antigen-specific CD4+ T lymphocyte responses contribute to protective immunity against Babesia bovis, however the antigens that induce these responses remain largely unknown. A proteomic approach was used to identify novel B. bovis antigens recognized by memory CD4+ T cells from immune cattle. Fractions obtained from merozoites separated by continuous-flow electrophoresis (CFE) that contained proteins ranging from 20 to 83 kDa were previously shown to stimulate memory CD4+ lymphocyte responses in B. bovis-immune cattle. Expression library screening with rabbit antiserum raised against an immunostimulatory CFE fraction identified a clone encoding a predicted 78 kDa protein. BLAST analysis revealed sequence identity of this B. bovis protein with Plasmodium falciparum fatty acyl coenzyme A synthetase (ACS) family members (PfACS1-PfACS11), and the protein was designated B. bovis acyl-CoA synthetase 1 (ACS1). Southern blot analysis indicated that B. bovis ACS1 is encoded by a single gene, although BLAST analysis of the preliminary B. bovis genome sequence identified two additional family members, ACS2 and ACS3. Peripheral blood lymphocytes and CD4+ T cell lines from B. bovis-immune cattle proliferated significantly against recombinant ACS1 protein, consistent with its predicted involvement in protective immunity. However, immune sera from cattle recovered from B. bovis infection did not react with ACS1, indicating that epitopes may be conformationally dependent.

Complete genome sequencing of Anaplasma marginale reveals that the surface is skewed to two superfamilies of outer membrane proteins

The rickettsia Anaplasma marginale is the most prevalent tick-borne livestock pathogen worldwide and is a severe constraint to animal health. A. marginale establishes lifelong persistence in infected ruminants and these animals serve as a reservoir for ticks to acquire and transmit the pathogen. Within the mammalian host, A. marginale generates antigenic variants by changing a surface coat composed of numerous proteins. By sequencing and annotating the complete 1,197,687-bp genome of the St. Maries strain of A. marginale, we show that this surface coat is dominated by two families containing immunodominant proteins: the msp2 superfamily and the msp1 superfamily. Of the 949 annotated coding sequences, just 62 are predicted to be outer membrane proteins, and of these, 49 belong to one of these two superfamilies. The genome contains unusual functional pseudogenes that belong to the msp2 superfamily and play an integral role in surface coat antigenic variation, and are thus distinctly different from pseudogenes described as byproducts of reductive evolution in other Rickettsiales.

End-sequence profiling: sequence-based analysis of aberrant genomes

Genome rearrangements are important in evolution, cancer, and other diseases. Precise mapping of the rearrangements is essential for identification of the involved genes, and many techniques have been developed for this purpose. We show here that end-sequence profiling (ESP) is particularly well suited to this purpose. ESP is accomplished by constructing a bacterial artificial chromosome (BAC) library from a test genome, measuring BAC end sequences, and mapping end-sequence pairs onto the normal genome sequence. Plots of BAC end-sequences density identify copy number abnormalities at high resolution. BACs spanning structural aberrations have end pairs that map abnormally far apart on the normal genome sequence. These pairs can then be sequenced to determine the involved genes and breakpoint sequences. ESP analysis of the breast cancer cell line MCF-7 demonstrated its utility for analysis of complex genomes. End sequencing of approximately 8,000 clones (0.37-fold haploid genome clonal coverage) produced a comprehensive genome copy number map of the MCF-7 genome at better than 300-kb resolution and identified 381 genome breakpoints, a subset of which was verified by fluorescence in situ hybridization mapping and sequencing.