Generation of an equine biobank to be used for Functional Annotation of Animal Genomes project.

The Functional Annotation of Animal Genomes (FAANG) project aims to identify genomic regulatory elements in both sexes across multiple stages of development in domesticated animals. This study represents the first stage of the FAANG project for the horse, Equus caballus. A biobank of 80 tissue samples, two cell lines and six body fluids was created from two adult Thoroughbred mares. Ante-mortem assessments included full physical examinations, lameness, ophthalmologic and neurologic evaluations. Complete blood counts and serum biochemistries were also performed. At necropsy, in addition to tissue samples, aliquots of serum, ethylenediaminetetraacetic acid (EDTA) plasma, heparinized plasma, cerebrospinal fluid, synovial fluid, urine and microbiome samples from all regions of the gastrointestinal and urogenital tracts were collected. Epidermal keratinocytes and dermal fibroblasts were cultured from skin samples. All tissues were grossly and histologically evaluated by a board-certified veterinary pathologist. The results of the clinical and pathological evaluations identified subclinical eosinophilic and lymphocytic infiltration throughout the length of the gastrointestinal tract as well as a mild clinical lameness in both animals. Each sample was cryo-preserved in multiple ways, and nuclei were extracted from selected tissues. These samples represent the first published systemically healthy equine-specific biobank with extensive clinical phenotyping ante- and post-mortem. The tissues in the biobank are intended for community-wide use in the functional annotation of the equine genome. The use of the biobank will improve the quality of the reference annotation and allow all equine researchers to elucidate unknown genomic and epigenomic causes of disease.

Identification of long non-coding RNA in the horse transcriptome.

BACKGROUND: Efforts to resolve the transcribed sequences in the equine genome have focused on protein-coding RNA. The transcription of the intergenic regions, although detected via total RNA sequencing (RNA-seq), has yet to be characterized in the horse. The most recent equine transcriptome based on RNA-seq from several tissues was a prime opportunity to obtain a concurrent long non-coding RNA (lncRNA) database. RESULTS: This lncRNA database has a breadth of eight tissues and a depth of over 20 million reads for select tissues, providing the deepest and most expansive equine lncRNA database. Utilizing the intergenic reads and three categories of novel genes from a previously published equine transcriptome pipeline, we better describe these groups by annotating the lncRNA candidates. These lncRNA candidates were filtered using an approach adapted from human lncRNA annotation, which removes transcripts based on size, expression, protein-coding capability and distance to the start or stop of annotated protein-coding transcripts. CONCLUSION: Our equine lncRNA database has 20,800 transcripts that demonstrate characteristics unique to lncRNA including low expression, low exon diversity and low levels of sequence conservation. These candidate lncRNA will serve as a baseline lncRNA annotation and begin to describe the RNA-seq reads assigned to the intergenic space in the horse.

Tissue resolved, gene structure refined equine transcriptome.

BACKGROUND: Transcriptome interpretation relies on a good-quality reference transcriptome for accurate quantification of gene expression as well as functional analysis of genetic variants. The current annotation of the horse genome lacks the specificity and sensitivity necessary to assess gene expression especially at the isoform level, and suffers from insufficient annotation of untranslated regions (UTR) usage. We built an annotation pipeline for horse and used it to integrate 1.9 billion reads from multiple RNA-seq data sets into a new refined transcriptome. RESULTS: This equine transcriptome integrates eight different tissues from 59 individuals and improves gene structure and isoform resolution, while providing considerable tissue-specific information. We utilized four levels of transcript filtration in our pipeline, aimed at producing several transcriptome versions that are suitable for different downstream analyses. Our most refined transcriptome includes 36,876 genes and 76,125 isoforms, with 6474 candidate transcriptional loci novel to the equine transcriptome. CONCLUSIONS: We have employed a variety of descriptive statistics and figures that demonstrate the quality and content of the transcriptome. The equine transcriptomes that are provided by this pipeline show the best tissue-specific resolution of any equine transcriptome to date and are flexible for several downstream analyses. We encourage the integration of further equine transcriptomes with our annotation pipeline to continue and improve the equine transcriptome.

Brother of CDO (BOC) expression in equine articular cartilage.

Brother of CDO (BOC) is a cell surface receptor that derives its name from the structurally related protein, cell adhesion molecule-related/down-regulated by oncogenes (CDO, sometimes CDON). High levels of BOC mRNA and protein expression have been described in embryonic tissues with active cell proliferation and ongoing cellular differentiation(1,2). A microarray-based screen of RNA isolated from 11 different adult equine tissues unexpectedly identified BOC as having an expression pattern restricted to articular cartilage. The objective of this study was to further investigate BOC expression in adult articular cartilage relative to other tissues. Both RT-qPCR and mRNA sequencing confirmed the microarray data. Steady state BOC mRNA levels in articular cartilage were substantially higher than in the other adult tissues tested, neonatal tendon, placenta, and whole embryo. The expression of BOC displayed a pattern of tissue specificity comparable to well established cartilage matrix protein biomarkers. BOC mRNA levels in articular cartilage increased with age, but were rapidly down-regulated when chondrocytes were enzymatically isolated from the cartilage matrix and expanded in monolayer culture. Relative expression patterns of CDO were broadly similar, but displayed lower fold change differences. A functional role in articular cartilage that involves Hedgehog signaling is suggested by the known binding affinity of BOC for all three Hedgehog ligands. These data also extend BOC and CDO biology to a post-mitotic and highly differentiated cell type within a mature tissue.

Structural annotation of equine protein-coding genes determined by mRNA sequencing.

The horse, like the majority of animal species, has a limited amount of species-specific expressed sequence data available in public databases. As a result, structural models for the majority of genes defined in the equine genome are predictions based on ab initio sequence analysis or the projection of gene structures from other mammalian species. The current study used Illumina-based sequencing of messenger RNA (RNA-seq) to help refine structural annotation of equine protein-coding genes and for a preliminary assessment of gene expression patterns. Sequencing of mRNA from eight equine tissues generated 293,758105 sequence tags of 35 bases each, equalling 10.28 gbp of total sequence data. The tag alignments represent approximately 207 × coverage of the equine mRNA transcriptome and confirmed transcriptional activity for roughly 90% of the protein-coding gene structures predicted by Ensembl and NCBI. Tag coverage was sufficient to refine the structural annotation for 11,356 of these predicted genes, while also identifying an additional 456 transcripts with exon/intron features that are not listed by either Ensembl or NCBI. Genomic locus data and intervals for the protein-coding genes predicted by the Ensembl and NCBI annotation pipelines were combined with 75,116 RNA-seq-derived transcriptional units to generate a consensus equine protein-coding gene set of 20,302 defined loci. Gene ontology annotation was used to compare the functional and structural categories of genes expressed in either a tissue-restricted pattern or broadly across all tissue samples.

Transcriptional regulation of the mosquito vitellogenin gene via a blood meal-triggered cascade.

In anautogenous mosquitoes, a blood meal is required for activation of genes encoding yolk protein precursors (YPP). Vitellogenin (Vg), the major YPP gene, is transcribed at a very high level following blood meal activation. It is expressed exclusively in the female fat body, the tissue producing most of mosquito hemolymph and immune proteins. In this paper, we analyzed the upstream region of the Aedes aegypti Vg gene in order to identify regulatory elements responsible for its unique expression pattern. To achieve this goal, we analyzed the gene using transgenic Drosophila and Aedes as well as DNA-binding assays. These analyses revealed three regulatory regions in the 2.1 kb upstream portion of the Vg gene. The proximal region containing binding sites to EcR/USP, GATA, C/EBP and HNF3/fkh is required for the correct tissue- and stage-specific expression at a low level. The median region carrying sites for early ecdysone response factors E74 and E75 is responsible for hormonal enhancement of Vg expression. Finally, the distal GATA-rich region is necessary for extremely high expression levels characteristic of the Vg gene. The present work elucidates the molecular basis of blood meal-dependent expression of this mosquito gene, laying the foundation for mosquito-specific expression cassettes with predictable stage and tissue specificity.