The pelvic flexure separates distinct microbial communities in the equine hindgut.

As hindgut fermenters, horses are especially dependent on the microbiota residing in their cecum and large intestines. Interactions between these microbial populations and the horse are critical for maintaining gut homeostasis, which supports proper digestion. The current project was motivated to determine if any features of the fecal microbiota are informative of the microbial communities from the cecum, ventral colon, or dorsal colon. Digesta from the cecum, ventral colon, dorsal colon and feces were collected from 6 yearling miniature horses. Microbial DNA was isolated and the microbiota from each sample was characterized by profiling the V4 region of the 16S rRNA. Principal coordinate analysis of the beta diversity results revealed significant (p = 0.0001; F = 5.2393) similarities between the microbial populations from cecal and ventral colon and the dorsal colon and fecal samples, however, there was little overlap between the proximal and distal ends of the hindgut. These distinct population structures observed in our results coincide with the pelvic flexure, which itself separates intestinal compartments with distinct roles in digestive physiology. An indicator species analysis confirmed the population differences, supported by the identification of several microbial families characteristic of the compartments upstream of the pelvic flexure that were not represented following it. Our data suggest that the fecal microbiota is not informative of the proximal hindgut but can provide insight into communities of the distal compartments. Further, our results suggest that the pelvic flexure might be an important anatomical landmark relative to the microbial communities in the equine large intestine.

Next-Generation Sequencing in Equine Genomics.

The sequencing and assembly of a reference genome for the horse has been revolutionary for investigation of horse health and performance. Next-generation sequencing (NGS) methods represent a second revolution in equine genomics. Researchers can align and compare DNA and RNA sequencing data to the reference genome to explore variation that may contribute or be attributed to disease. NGS has also facilitated the translation of research discovery to clinically relevant applications. This article discusses the history and development of NGS, details some of the available sequencing platforms, and describes currently available applications in the context of both discovery and clinical settings.

Transfer of MicroRNAs From Epididymal Epithelium to Equine Spermatozoa.

All epididymal regions are lined with multiple epithelial cell types, each with different functions to provide the luminal environment for spermatozoal maturation. Epithelial cells also create apical blebs, which are released from the apical surface via apocrine secretion and disintegrate in the lumen, thereby releasing epididymosomes. Epididymosomes transport proteins to spermatozoa and contain microRNAs. We hypothesized that epididymosomes also transfer miRNA from epididymal epithelium to spermatozoa. Quantitative real-time polymerase chain reaction was used to determine miRNA profiles of epididymal tissue from caput and cauda, epididymal spermatozoa from caput and cauda, and epididymosomes and from caput, proximal corpus, distal corpus, and cauda. Pathway analysis was performed using DIANA tools on the miRNA unique to caudal spermatozoa. We found 66 newly acquired miRNAs in spermatozoa located in the caudal epididymis. Predicted pathways targeted by these miRNAs suggest a role in cell motility and viability and factors in oocyte and embryo maturation and development. These findings suggest that miRNAs are transported to spermatozoa from epididymal epithelium via epididymosomes.

An isothermal shift assay for proteome scale drug-target identification.

Most small molecule drugs act on living systems by physically interacting with specific proteins and modulating target function. Identification of drug binding targets, within the complex milieu of the human proteome, remains a challenging task of paramount importance in drug discovery. Existing approaches for target identification employ complex workflows with limited throughput. Here, we present the isothermal shift assay (iTSA), a mass spectrometry method for proteome-wide identification of drug targets within lysates or living cells. Compared with prevailing methods, iTSA uses a simplified experimental design with increased statistical power to detect thermal stability shifts that are induced by small molecule binding. Using a pan-kinase inhibitor, staurosporine, we demonstrate improved performance over commonly used thermal proteome profiling methods, identifying known targets in cell lysates and living cells. We also demonstrate the identification of both known targets and additional candidate targets for the kinase inhibitor harmine in cell and tissue lysates.

Non-Coding RNA Sequencing of Equine Endometrium During Maternal Recognition of Pregnancy.

Maternal recognition of pregnancy (MRP) in the mare is not well defined. In a non-pregnant mare, prostaglandin F2α (PGF) is released on day 14 post-ovulation (PO) to cause luteal regression, resulting in loss of progesterone production. Equine MRP occurs prior to day 14 to halt PGF production. Studies have failed to identify a gene candidate for MRP, so attention has turned to small, non-coding RNAs. The objective of this study was to evaluate small RNA (<200 nucleotides) content in endometrium during MRP. Mares were used in a cross-over design with each having a pregnant and non-mated cycle. Each mare was randomly assigned to collection day 11 or 13 PO (n = 3/day) and endometrial biopsies were obtained. Total RNA was isolated and sequencing libraries were prepared using a small RNA library preparation kit and sequenced on a HiSeq 2000. EquCab3 was used as the reference genome and DESeq2 was used for statistical analysis. On day 11, 419 ncRNAs, representing miRNA, snRNA, snoRNA, scaRNA, and vaultRNA, were different between pregnancy statuses, but none on day 13. Equine endometrial ncRNAs with unknown structure and function were also identified. This study is the first to describe ncRNA transcriptome in equine endometrium. Identifying targets of these ncRNAs could lead to determining MRP.

Coding RNA Sequencing of Equine Endometrium during Maternal Recognition of Pregnancy.

Equine maternal recognition of pregnancy (MRP) is a process whose signal remains unknown. During MRP the conceptus and endometrium communicate to attenuate prostaglandin F2α (PGF) secretion, sparing the corpus luteum and maintaining progesterone production. Recognition of a mobile conceptus by the endometrium is critical by days 14-16 post-ovulation (PO), when endometrium produces PGF, initiating luteolysis. The objective of this study was to evaluate endometrial gene expression changes based upon pregnancy status via RNA sequencing. This experiment utilized a cross-over design with each mare serving as both a pregnant and non-mated control on days nine, 11, and 13 PO (n = 3/status/day). Mares were randomly assigned to collection day and pregnancy confirmed by terminal uterine lavage at the time of endometrial biopsy. Total RNA was isolated and libraries prepared using Illumina TruSeq RNA sample preparation kit. Reads were mapped and annotated using HISAT2 and Stringtie. Expression values were evaluated with DESEQ2 (P ≤ 0.05 indicated significance). On day nine, 11, and 13 there were 1435, 1435 and 916 significant transcripts, respectively. Multiple genes with splice variants had different expression patterns within the same day. These are the first data to evaluate the endometrial transcriptome during MRP on days nine, 11, and 13.

Equine Fecal Microbiota Changes Associated With Anthelmintic Administration.

The gastrointestinal microbiota (GIM) plays an essential role in maintaining intestinal homeostasis with disruptions having profound effects on the wellbeing of the host animal. Parasitic infection is a long-standing issue for the equine industry, and the use of anthelmintic drugs for parasite control has long been standard practice. The impact of anthelmintic treatment on the GIM in healthy horses is not well known. This study evaluated the hypothesis that anthelmintic administration will alter the equine fecal microbiota in horses without an observed helminth infection. Ten horses were treated with a single dose of QUEST PLUS (active ingredients: Moxidectin and Praziquantel) (Zoetis), and fecal samples were collected before and after treatment. Amplicon sequencing data were quality filtered, processed, and analyzed using QIIME2. Anthelmintic treatment corresponded with a small but significant decrease in alpha diversity (P-value < .05). Analysis of taxonomic abundances before and after treatment with DESeq2 identified 21 features that were significantly different after treatment (Padj-value < .05). Differences in beta diversity associated with treatment were not significant and potentially suggest factors unique to the individual may play an essential role in the specific responses observed. Overall, the present study does not indicate a broad, large-scale impact on the GIM after anthelmintic treatment. The results do, however, suggest the potential of individualized responses that are based instead on host factors. Identification of these factors and investigation of their impact on the host/microbiota relationship will contribute significantly to our understanding of the role of the microbiome in horse health.

Tissue Restricted Splice Junctions Originate Not Only from Tissue-Specific Gene Loci, but Gene Loci with a Broad Pattern of Expression.

Cellular mechanisms that achieve protein diversity in eukaryotes are multifaceted, including transcriptional components such as RNA splicing. Through alternative splicing, a single protein-coding gene can generate multiple mRNA transcripts and protein isoforms, some of which are tissue-specific. We have conducted qualitative and quantitative analyses of the Bodymap 2.0 messenger RNA-sequencing data from 16 human tissue samples and identified 209,363 splice junctions. Of these, 22,231 (10.6%) were not previously annotated and 21,650 (10.3%) were expressed in a tissue-restricted pattern. Tissue-restricted alternative splicing was found to be widespread, with approximately 65% of expressed multi-exon genes containing at least one tissue-specific splice junction. Interestingly, we observed many tissue-specific splice junctions not only in genes expressed in one or a few tissues, but also from gene loci with a broad pattern of expression.

Annotation of the Protein Coding Regions of the Equine Genome.

Current gene annotation of the horse genome is largely derived from in silico predictions and cross-species alignments. Only a small number of genes are annotated based on equine EST and mRNA sequences. To expand the number of equine genes annotated from equine experimental evidence, we sequenced mRNA from a pool of forty-three different tissues. From these, we derived the structures of 68,594 transcripts. In addition, we identified 301,829 positions with SNPs or small indels within these transcripts relative to EquCab2. Interestingly, 780 variants extend the open reading frame of the transcript and appear to be small errors in the equine reference genome, since they are also identified as homozygous variants by genomic DNA resequencing of the reference horse. Taken together, we provide a resource of equine mRNA structures and protein coding variants that will enhance equine and cross-species transcriptional and genomic comparisons.

Comparison of the Equine Reference Sequence with Its Sanger Source Data and New Illumina Reads.

The reference assembly for the domestic horse, EquCab2, published in 2009, was built using approximately 30 million Sanger reads from a Thoroughbred mare named Twilight. Contiguity in the assembly was facilitated using nearly 315 thousand BAC end sequences from Twilight's half brother Bravo. Since then, it has served as the foundation for many genome-wide analyses that include not only the modern horse, but ancient horses and other equid species as well. As data mapped to this reference has accumulated, consistent variation between mapped datasets and the reference, in terms of regions with no read coverage, single nucleotide variants, and small insertions/deletions have become apparent. In many cases, it is not clear whether these differences are the result of true sequence variation between the research subjects' and Twilight's genome or due to errors in the reference. EquCab2 is regarded as "The Twilight Assembly." The objective of this study was to identify inconsistencies between the EquCab2 assembly and the source Twilight Sanger data used to build it. To that end, the original Sanger and BAC end reads have been mapped back to this equine reference and assessed with the addition of approximately 40X coverage of new Illumina Paired-End sequence data. The resulting mapped datasets identify those regions with low Sanger read coverage, as well as variation in genomic content that is not consistent with either the original Twilight Sanger data or the new genomic sequence data generated from Twilight on the Illumina platform. As the haploid EquCab2 reference assembly was created using Sanger reads derived largely from a single individual, the vast majority of variation detected in a mapped dataset comprised of those same Sanger reads should be heterozygous. In contrast, homozygous variations would represent either errors in the reference or contributions from Bravo's BAC end sequences. Our analysis identifies 720,843 homozygous discrepancies between new, high throughput genomic sequence data generated for Twilight and the EquCab2 reference assembly. Most of these represent errors in the assembly, while approximately 10,000 are demonstrated to be contributions from another horse. Other results are presented that include the binary alignment map file of the mapped Sanger reads, a list of variants identified as discrepancies between the source data and resulting reference, and a BED annotation file that lists the regions of the genome whose consensus was likely derived from low coverage alignments.

Determinants and implications of mRNA poly(A) tail size--does this protein make my tail look big?

While the phenomenon of polyadenylation has been well-studied, the dynamics of poly(A) tail size and its impact on transcript function and cell biology are less well-appreciated. The goal of this review is to encourage readers to view the poly(A) tail as a dynamic, changeable aspect of a transcript rather than a simple static entity that marks the 3' end of an mRNA. This could open up new angles of regulation in the post-transcriptional control of gene expression throughout development, differentiation and cancer.

Analysis of unannotated equine transcripts identified by mRNA sequencing.

Sequencing of equine mRNA (RNA-seq) identified 428 putative transcripts which do not map to any previously annotated or predicted horse genes. Most of these encode the equine homologs of known protein-coding genes described in other species, yet the potential exists to identify novel and perhaps equine-specific gene structures. A set of 36 transcripts were prioritized for further study by filtering for levels of expression (depth of RNA-seq read coverage), distance from annotated features in the equine genome, the number of putative exons, and patterns of gene expression between tissues. From these, four were selected for further investigation based on predicted open reading frames of greater than or equal to 50 amino acids and lack of detectable homology to known genes across species. Sanger sequencing of RT-PCR amplicons from additional equine samples confirmed expression and structural annotation of each transcript. Functional predictions were made by conserved domain searches. A single transcript, expressed in the cerebellum, contains a putative kruppel-associated box (KRAB) domain, suggesting a potential function associated with zinc finger proteins and transcriptional regulation. Overall levels of conserved synteny and sequence conservation across a 1MB region surrounding each transcript were approximately 73% compared to the human, canine, and bovine genomes; however, the four loci display some areas of low conservation and sequence inversion in regions that immediately flank these previously unannotated equine transcripts. Taken together, the evidence suggests that these four transcripts are likely to be equine-specific.

Genome-wide association study among four horse breeds identifies a common haplotype associated with in vitro CD3+ T cell susceptibility/resistance to equine arteritis virus infection.

Previously, we have shown that horses could be divided into susceptible and resistant groups based on an in vitro assay using dual-color flow cytometric analysis of CD3+ T cells infected with equine arteritis virus (EAV). Here, we demonstrate that the differences in in vitro susceptibility of equine CD3+ T lymphocytes to EAV infection have a genetic basis. To investigate the possible hereditary basis for this trait, we conducted a genome-wide association study (GWAS) to compare susceptible and resistant phenotypes. Testing of 267 DNA samples from four horse breeds that had a susceptible or a resistant CD3+ T lymphocyte phenotype using both Illumina Equine SNP50 BeadChip and Sequenom's MassARRAY system identified a common, genetically dominant haplotype associated with the susceptible phenotype in a region of equine chromosome 11 (ECA11), positions 49572804 to 49643932. The presence of a common haplotype indicates that the trait occurred in a common ancestor of all four breeds, suggesting that it may be segregated among other modern horse breeds. Biological pathway analysis revealed several cellular genes within this region of ECA11 encoding proteins associated with virus attachment and entry, cytoskeletal organization, and NF-κB pathways that may be associated with the trait responsible for the in vitro susceptibility/resistance of CD3+ T lymphocytes to EAV infection. The data presented in this study demonstrated a strong association of genetic markers with the trait, representing de facto proof that the trait is under genetic control. To our knowledge, this is the first GWAS of an equine infectious disease and the first GWAS of equine viral arteritis.

MapSplice: accurate mapping of RNA-seq reads for splice junction discovery.

The accurate mapping of reads that span splice junctions is a critical component of all analytic techniques that work with RNA-seq data. We introduce a second generation splice detection algorithm, MapSplice, whose focus is high sensitivity and specificity in the detection of splices as well as CPU and memory efficiency. MapSplice can be applied to both short (<75 bp) and long reads (≥ 75 bp). MapSplice is not dependent on splice site features or intron length, consequently it can detect novel canonical as well as non-canonical splices. MapSplice leverages the quality and diversity of read alignments of a given splice to increase accuracy. We demonstrate that MapSplice achieves higher sensitivity and specificity than TopHat and SpliceMap on a set of simulated RNA-seq data. Experimental studies also support the accuracy of the algorithm. Splice junctions derived from eight breast cancer RNA-seq datasets recapitulated the extensiveness of alternative splicing on a global level as well as the differences between molecular subtypes of breast cancer. These combined results indicate that MapSplice is a highly accurate algorithm for the alignment of RNA-seq reads to splice junctions. Software download URL: http://www.netlab.uky.edu/p/bioinfo/MapSplice.

Characterization of equine and other vertebrate TLR3, TLR7, and TLR8 genes.

Toll-like receptors 3, 7, and 8 (TLR3, TLR7, and TLR8) were studied in the genomes of the domestic horse and several other mammals. The messenger RNA sequences and exon/intron structures of these TLR genes were determined. An equine bacterial artificial chromosome clone containing the TLR3 gene was assigned by fluorescent in situ hybridization to the horse chromosomal location ECA27q16-q17 and this map location was confirmed using an equine radiation hybrid panel. Direct sequencing revealed 13 single-nucleotide polymorphisms in the coding regions of the equine TLR 3, 7, and 8 genes. Of these polymorphisms, 12 were not previously reported. The allelic frequency was estimated for each single-nucleotide polymorphism from genotyping data obtained for 154 animals from five horse breeds. Some of these frequencies varied significantly among different horse breeds. Domain architecture predictions for the three equine TLR protein sequences revealed several conserved regions within the variable leucine-rich repeats between the corresponding horse and cattle TLR proteins. A phylogenetic analysis did not indicate that any significant exchanges had occurred between paralogous TLR7 and TLR8 genes in 20 vertebrate species analyzed.