RNA-seq transcriptome profiling of equine inner cell mass and trophectoderm

Unveiling the population genetic structure of Iranian horses breeds by whole-genome resequencing analysis

Preserving genetic diversity is pivotal for enhancing genetic improvement and facilitating adaptive responses to selection. This study focuses on identifying key genetic variants, including single nucleotide polymorphisms (SNPs), insertion/deletion polymorphisms (INDELs), and copy number variants (CNVs), while exploring the genomic evolutionary connectedness among seven Iranian horses representing five indigenous breeds: Caspian, Turkemen, DareShuri, Kurdish, and Asil. Using whole-genome resequencing, we generated 2.7 Gb of sequence data, with raw reads ranging from 1.2 Gb for Caspian horses to 0.38 Gb for Turkoman horses. Post-filtering, approximately 1.9 Gb of reads remained, with ~ 1.5 Gb successfully mapped to the horse reference genome (EquCab3.0), achieving mapping rates between 76.4% (Caspian) and 98.35% (Turkoman). We identified 2,909,816 SNPs in Caspian horses, constituting around 0.1% of the genome. Notably, 71% of these SNPs were situated in intergenic regions, while 8.5 and 6.8% were located upstream and downstream, respectively. A comparative analysis of SNPs between Iranian and non-Iranian horse breeds showed that Caspian horses had the lowest number of shared SNPs with Turkoman horses. Instead, they showed a closer genetic relationship with DareShuri, Quarter, Arabian, Standardbred, and Asil breeds. Hierarchical clustering highlighted Caspian horses as a distinct cluster, underscoring their distinctive genomic signature. Caspian horses exhibit a unique genetic profile marked by an enrichment of private mutations in neurological genes, influencing sensory perception and awareness. This distinct genetic makeup shapes mating preferences and signifies a separate evolutionary trajectory. Additionally, significant non-synonymous single nucleotide polymorphisms (nsSNPs) in reproductive genes offer intervention opportunities for managing Caspian horses. These findings reveal the population genetic structure of Iranian horse breeds, contributing to the advancement of knowledge in areas such as conservation, performance traits, climate adaptation, reproduction, and resistance to diseases in equine science.

MUC1 regulation in the left and right uterine horns and conceptus trophectoderm during the peri-implantation period of dromedary camel

Pregnancy maintenance in dromedary camels poses significant challenges, including early embryonic loss in the left uterine horn (LH) and unsuccessful pregnancy in the right uterine horn (RH), suggesting a potential asynchrony between conceptus signaling and uterine receptivity. The transition of the uterine epithelium from a pre-receptive to a receptive state requires a delicate balance of adhesion-promoting and anti-adhesion molecules. Mucin-1 (MUC1) acts as an anti-adhesive molecule on the uterine luminal (LE) and glandular (GE) epithelium. Downregulation of MUC1 is believed to be crucial for successful embryo attachment in various mammals. This study aimed to investigate the temporospatial expression of MUC1 in the LH and RH on Days 8, 10, and 12 pregnant dromedaries and their conceptuses. Quantitative real-time polymerase chain reaction (qrt-PCR), Western blot analysis, immunohistochemistry, and immunofluorescence techniques were employed to assess MUC1 expression at the mRNA and protein levels. The results demonstrated a reduction in MUC1 mRNA expression on Day 8, then increased on Day 10, followed by a decrease on Day 12 in LH. While the RH exhibited progressive increases, peaking on Day 12. However, MUC1 expression constantly exhibited higher levels in RH than in LH in all days. Two bands were detected at 150-kDa and 180-kDa, with the highest intensity observed on Day 10. Spatially, MUC1 was localized in the apical, cytoplasmic, and lumen of uterine glands only. MUC1 was barely detectable on Day 8 but gradually increased on Days 10 and 12 in both horns. Likewise, the RH exhibited higher MUC1 signals than the LH on Days 10 and 12. In the conceptuses, MUC1 mRNA increased on Day 8, peaked on Day 10, and declined on Day 12. Notably, MUC1 protein was detected in both the trophectoderm and endoderm, with high expression observed on Day 10 and reduced by Day 12. In conclusion, the decrease in MUC1 expression on Day 8 in the LH may be associated with maternal recognition of pregnancy (MRP), and the increase on Day 10 may related to embryo protection and movement, while the subsequent decrease on Day 12 could be linked to the embryo attachment and preparation for the implantation. Conversely, the increase of MUC1 in the RH implies a role in the anti-adhesion mechanism. These findings contribute to understanding MUC1's involvement in reproductive processes and provide insights into the complex mechanisms underlying successful pregnancy establishment and maintenance in dromedary camels.

Can Nicotinamide Adenine Dinucleotide (NAD+) and Sirtuins Be Harnessed to Improve Mare Fertility?

Years of sire and dam selection based on their pedigree and athletic performance has resulted in a reduction in the reproductive capability of horses. Mare age is considered a major barrier to equine reproduction largely due to an increase in the age at which mares are typically bred following the end of their racing career. Nicotinamide adenine dinucleotide (NAD+) and its involvement in the activation of Sirtuins in fertility are an emerging field of study, with the role of NAD+ in oocyte maturation and embryo development becoming increasingly apparent. While assisted reproductive technologies in equine breeding programs are in their infancy compared to other livestock species such as cattle, there is much more to be learnt, from oocyte maturation to early embryo development and beyond in the mare, which are difficult to study given the complexities associated with mare fertility research. This review examines what is already known about the role of NAD+ and Sirtuins in fertility and discusses how NAD+-elevating agents may be used to activate Sirtuin proteins to improve equine breeding and embryo production programs both in vivo and in vitro.

In Vitro-Produced Equine Blastocysts Exhibit Greater Dispersal and Intermingling of Inner Cell Mass Cells than In Vivo Embryos

In vitro production (IVP) of equine embryos is increasingly popular in clinical practice but suffers from higher incidences of early embryonic loss and monozygotic twin development than transfer of in vivo derived (IVD) embryos. Early embryo development is classically characterized by two cell fate decisions: (1) first, trophectoderm (TE) cells differentiate from inner cell mass (ICM); (2) second, the ICM segregates into epiblast (EPI) and primitive endoderm (PE). This study examined the influence of embryo type (IVD versus IVP), developmental stage or speed, and culture environment (in vitro versus in vivo) on the expression of the cell lineage markers, CDX-2 (TE), SOX-2 (EPI) and GATA-6 (PE). The numbers and distribution of cells expressing the three lineage markers were evaluated in day 7 IVD early blastocysts (n = 3) and blastocysts (n = 3), and in IVP embryos first identified as blastocysts after 7 (fast development, n = 5) or 9 (slow development, n = 9) days. Furthermore, day 7 IVP blastocysts were examined after additional culture for 2 days either in vitro (n = 5) or in vivo (after transfer into recipient mares, n = 3). In IVD early blastocysts, SOX-2 positive cells were encircled by GATA-6 positive cells in the ICM, with SOX-2 co-expression in some presumed PE cells. In IVD blastocysts, SOX-2 expression was exclusive to the compacted presumptive EPI, while GATA-6 and CDX-2 expression were consistent with PE and TE specification, respectively. In IVP blastocysts, SOX-2 and GATA-6 positive cells were intermingled and relatively dispersed, and co-expression of SOX-2 or GATA-6 was evident in some CDX-2 positive TE cells. IVP blastocysts had lower TE and total cell numbers than IVD blastocysts and displayed larger mean inter-EPI cell distances; these features were more pronounced in slower-developing IVP blastocysts. Transferring IVP blastocysts into recipient mares led to the compaction of SOX-2 positive cells into a presumptive EPI, whereas extended in vitro culture did not. In conclusion, IVP equine embryos have a poorly compacted ICM with intermingled EPI and PE cells; features accentuated in slowly developing embryos but remedied by transfer to a recipient mare.

Runs of homozygosity and signatures of selection for number of oocytes and embryos in the Gir Indicine cattle

Runs of homozygosity (ROH) and signatures of selection are the results of selection processes in livestock species that have been shown to affect several traits in cattle. The aim of the current work was to verify the profile of ROH and inbreeding depression in the number of total (TO) and viable oocytes (VO) and the number of embryos (EMBR) in Gir Indicine cattle. In addition, we aim to identify signatures of selection, genes, and enriched regions between Gir subpopulations sorted by breeding value for these traits. The genotype file contained 2093 animals and 420,718 SNP markers. Breeding values used to sort Gir animals were previously obtained. ROH and signature of selection analyses were performed using PLINK software, followed by ROH-based (F_ROH) and pedigree-based inbreeding (F_ped) and a search for genes and their functions. An average of 50 ± 8.59 ROHs were found per animal. ROHs were separated into classes according to size, ranging from 1 to 2 Mb (ROH_1-2Mb: 58.17%), representing ancient inbreeding, ROH_2-4Mb (22.74%), ROH_4-8Mb (11.34%), ROH_8-16Mb (5.51%), and ROH_>16Mb (2.24%). Combining our results, we conclude that the increase in general F_ROH and F_ped significantly decreases TO and VO; however, in different chromosomes traits can increase or decrease with F_ROH. In the analysis for signatures of selection, we identified 15 genes from 47 significant genomic regions, indicating differences in populations with high and low breeding value for the three traits.

Dynamic regulation of the transcriptome and proteome of the equine embryo during maternal recognition of pregnancy

During initial maternal recognition of pregnancy (MRP), the equine embryo displays a series of unique events characterized by rapid blastocyst expansion, secretion of a diverse array of molecules, and transuterine migration to interact with the uterine surface. Up to date, the intricate transcriptome and proteome changes of the embryo underlying these events have not been critically studied in horses. Thus, the objective of this study was to perform an integrative transcriptomic (including mRNA, miRNAs, and other small non-coding RNAs) and proteomic analysis of embryos collected from days 10 to 13 of gestation. The results revealed dynamic transcriptome profiles with a total of 1311 differentially expressed genes, including 18 microRNAs (miRNAs). Two main profiles for mRNAs and miRNAs were identified, one with higher expression in embryos ≤5 mm and the second with higher expression in embryos ≥7 mm. At the protein level, similar results were obtained, with 259 differentially abundant proteins between small and large embryos. Overall, the findings demonstrated fine-tuned transcriptomic and proteomic regulations in the developing embryo associated with embryo growth. The identification of specific regulation of mRNAs, proteins, and miRNAs on days 12 and 13 of gestation suggested these molecules as pivotal for embryo development and as involved in MRP, and in establishment of pregnancy in general. In addition, the results revealed new insights into prostaglandin synthesis by the equine embryo, miRNAs and genes potentially involved in modulation of the maternal immune response, regulation of endometrial receptivity and of late implantation in the mare.

The transcriptome landscapes of allantochorion and vitelline-chorion in equine day 30 conceptus

During equine early gestation, trophectoderm forms chorion tissue, which is composed of two parts that one is covering allantoin, called allantochorion (AC) and another is covering yolk sac, which here we call vitelline-chorion (VC). Given that little is known about the equine trophoblast-derived chorion differentiation at an early stage, we first compared the transcriptome of AC and VC of day 30 equine conceptus based on RNA-sequencing. As a result, we found that compared to VC, there are 484 DEGs, including 305 up- and 179 down-regulated genes in AC. GO and KEGG analysis indicated that up-regulated genes in AC are mainly cell proliferation and cell adhesion-related genes, participating in allantois expansion and allantochorionic-placenta formation; dominant genes in VC are extracellular exosome and other cell adhesion-related genes implicated in direct and indirect conceptus-maternal communication. Additionally, as for the progenitor chorion tissue of equine chorionic gonadotropin secreting endometrium cup-the chorionic girdle (CG), which locates at the junction of the dilating AC and regressing VC, we revealed its unique gene expression pattern and the gene regulation during its further differentiation in vitro. Collectively, this study sheds light on the molecular events regarding the trophoblast differentiation and function at an early stage of the equine preimplantation conceptus.

Maternal age affects equine day 8 embryo gene expression both in trophoblast and inner cell mass

BACKGROUND

Breeding a mare until she is not fertile or even until her death is common in equine industry but the fertility decreases as the mare age increases. Embryo loss due to reduced embryo quality is partly accountable for this observation. Here, the effect of mare's age on blastocysts' gene expression was explored. Day 8 post-ovulation embryos were collected from multiparous young (YM, 6-year-old, N = 5) and older (OM, > 10-year-old, N = 6) non-nursing Saddlebred mares, inseminated with the semen of one stallion. Pure or inner cell mass (ICM) enriched trophoblast, obtained by embryo bisection, were RNA sequenced. Deconvolution algorithm was used to discriminate gene expression in the ICM from that in the trophoblast. Differential expression was analyzed with embryo sex and diameter as cofactors. Functional annotation and classification of differentially expressed genes and gene set enrichment analysis were also performed.

RESULTS

Maternal aging did not affect embryo recovery rate, embryo diameter nor total RNA quantity. In both compartments, the expression of genes involved in mitochondria and protein metabolism were disturbed by maternal age, although more genes were affected in the ICM. Mitosis, signaling and adhesion pathways and embryo development were decreased in the ICM of embryos from old mares. In trophoblast, ion movement pathways were affected.

CONCLUSIONS

This is the first study showing that maternal age affects gene expression in the equine blastocyst, demonstrating significant effects as early as 10 years of age. These perturbations may affect further embryo development and contribute to decreased fertility due to aging.

Transcript profiling of bovine embryos implicates specific transcription factors in the maternal-to-embryo transition

Full-grown oocytes are transcriptionally quiescent. Following maturation and fertilization, the early stages of embryonic development occur in the absence (or low levels) of transcription that results in a period of development relying on maternally derived products (e.g., mRNAs and proteins). Two critical steps occur during the transition from maternal to embryo control of development: maternal mRNA clearance and embryonic genome activation with an associated dramatic reprogramming of gene expression required for further development. By combining an RNA polymerase II inhibitor with RNA sequencing, we were able not only to distinguish maternally derived from embryonic transcripts in bovine preimplantation embryos but also to establish that embryonic gene activation is required for clearance of maternal mRNAs as well as to identify putative transcription factors that are likely critical for early bovine development.

Paternally expressed retrotransposon Gag-like 1 gene, RTL1, is one of the crucial elements for placental angiogenesis in horses†

RTL1 (retrotransposon Gag-like 1) is an essential gene in the development of the human and murine placenta. Several fetal and placental abnormalities such as intrauterine growth restriction (IUGR) and hydrops conditions have been associated with altered expression of this gene. However, the function of RTL1 has not been identified. RTL1 is located on a highly conserved region in eutherian mammals. Therefore, the genetic and molecular analysis in horses could hold important implications for other species, including humans. Here, we demonstrated that RTL1 is paternally expressed and is localized within the endothelial cells of the equine (Equus caballus) chorioallantois. We developed an equine placental microvasculature primary cell culture and demonstrated that RTL1 knockdown leads to loss of the sprouting ability of these endothelial cells. We further demonstrated an association between abnormal expression of RTL1 and development of hydrallantois. Our data suggest that RTL1 may be essential for placental angiogenesis, and its abnormal expression can lead to placental insufficiency. This placental insufficiency could be the reason for IUGR and hydrops conditions reported in other species, including humans.

Changes in testicular gene expression following reduced estradiol synthesis: A complex pathway to increased porcine Sertoli cell proliferation

Porcine Sertoli cell number including number present at puberty is increased if testicular estradiol synthesis is reduced during the neonatal interval. Evaluating the changes in gene expression during the crucial interval of suppressed estradiol that leads to the increased Sertoli cell population will increase our understanding of Sertoli cell biology but this evaluation first required a more precise determination of the critical interval for treatment and timing of a detectable response. Previously, reduced testicular estrogens from 1 week of age were accompanied by increased Sertoli cell number at 6.5 weeks of age but the age at which Sertoli cell numbers were initially increased was unknown, one of the current objectives. Additional experiments were designed to further delineate the essential timing of treatment for the Sertoli cell response. Finally, changes in gene expression induced by the reduced estradiol synthesis were evaluated to elucidate molecular mechanisms. Experimental design typically consisted of one member of littermate pairs of boars treated with the aromatase inhibitor, letrozole, beginning at 1 week of age and the remaining member treated with canola oil vehicle. Weekly treatments continued through 5 weeks of age or tissue collection, whichever came first. Increases in Sertoli cell numbers were not detectable prior to 6.5 weeks of age and persistent treatment through 5 weeks of age was required to induce the increase in Sertoli cell numbers. This increase resulted from prolonging the first interval of Sertoli cell proliferation in the treated animals. Few genes exhibited dramatically altered transcription and similarities in pathway analysis or principal modified genes were quite limited in 2, 3, and 5-week-old boars. The critical timing and prolonged treatment required and the sequential changes in gene expression suggest a complex mechanism is involved in this model of increased proliferation of Sertoli cells.

A Review of OCT4 Functions and Applications to Equine Embryos

OCT4 is a core transcription factor involved in pluripotency maintenance in the early mammalian embryo. The POU5F1 gene that encodes the OCT4 protein is highly conserved across species, suggesting conserved function. However, studies in several species including mice, cattle, and pigs, suggest that there are differences in where and when OCT4 is expressed. Specifically, in the horse, several studies have shown that exposure to the uterine environment may be necessary to induce OCT4 expression restriction to the inner cell mass (ICM) of the developing embryo, suggesting that there may be equine-specific extrinsic regulators of OCT4 expression that have not yet been investigated. However, an alternative hypothesis is that this restriction may not be evident in equine embryos because of our inability to culture them to the epiblast stage, preventing the observation of this restriction. In vitro studies have identified that OCT4 is expressed in the immature equine oocyte and in the early equine embryo, but OCT4 expression has not been studied after the formation of the ICM in the equine embryo. Despite the gaps in knowledge about equine-specific functions of OCT4, this factor has been used in studies assessing equine embryonic stem cells and to induce pluripotency in equine somatic cells. This review describes the role of OCT4 in the equine embryo and its applications in equine stem cell research.

Placentation in Equids

This chapter focuses on the early stages of placental development in horses and their relatives in the genus Equus and highlights unique features of equid reproductive biology. The equine placenta is classified as a noninvasive, epitheliochorial type. However, equids have evolved a minor component of invasive trophoblast, the chorionic girdle and endometrial cups, which links the equine placenta with the highly invasive hemochorial placentae of rodents and, particularly, with the primate placenta. Two types of fetus-to-mother signaling in equine pregnancy are mediated by the invasive equine trophoblast cells. First, endocrinological signaling mediated by equine chorionic gonadotrophin (eCG) drives maternal progesterone production to support the equine conceptus between days 40 and 100 of gestation. Only in primates and equids does the placenta produce a gonadotrophin, but the evolutionary paths taken by these two groups of mammals to produce this placental signal were very different. Second, florid expression of paternal major histocompatibility complex (MHC) class I molecules by invading chorionic girdle cells stimulates strong maternal anti-fetal antibody responses that may play a role in the development of immunological tolerance that protects the conceptus from destruction by the maternal immune system. In humans, invasive extravillous trophoblasts also express MHC class I molecules, but the loci involved, and their likely function, are different from those of the horse. Comparison of the cellular and molecular events in these disparate species provides outstanding examples of convergent evolution and co-option in mammalian pregnancy and highlights how studies of the equine placenta have produced new insights into reproductive strategies.

ZC3H4-a novel Cys-Cys-Cys-His-type zinc finger protein-is essential for early embryogenesis in mice†

Zinc finger domains of the Cys-Cys-Cys-His (CCCH) class are evolutionarily conserved proteins that bind nucleic acids and are involved in various biological processes. Nearly 60 CCCH-type zinc finger proteins have been identified in humans and mice, most have not been functionally characterized. Here, we provide the first in vivo functional characterization of ZC3H4-a novel CCCH-type zinc finger protein. Our results show that although Zc3h4 mutant embryos exhibit normal morphology at E3.5 blastocyst stage, they cannot be recovered at E7.5 early post-gastrulation stage, suggesting implantation failure. Outgrowth assays reveal that mutant blastocysts either fail to hatch from the zona pellucida, or can hatch but do not form a typical inner cell mass colony, the source of embryonic stem cells (ESCs). Although there is no change in levels of reactive oxygen species, Zc3h4 mutants display severe DNA breaks and reduced cell proliferation. Analysis of lineage specification reveals that both epiblast and primitive endoderm lineages are compromised with severe reductions in cell number and/or specification in the mutant blastocysts. In summary, these findings demonstrate the essential role of ZC3H4 during early mammalian embryogenesis.

Uterine infusion of conceptus fragments changes the protein profile from cyclic mares

This experiment aimed to compare at day seven after ovulation, the protein profile of uterine fluid in cyclic mares with mares infused two days before with Day 13 conceptus fragments. Experimental animals were ten healthy cyclic mares, examined daily to detect ovulation (Day 0) as soon as estrus was confirmed. On day seven, after ovulation, uterine fluid was collected, constituting the Cyclic group (n = 10). The same mares were examined in the second cycle until ovulation was detected. On day five, after ovulation, fragments from a previously collected concepti were infused into each mare's uterus. Two days after infusion, uterine fluid was collected, constituting the Fragment group (n = 10). Two-dimensional electrophoresis technique processed uterine fluid samples. A total of 373 spots were detected. MALDI-TOF/TOF and NanoUHPLC-QTOF mass spectrometry identified twenty spots with differences in abundance between the Cyclic and Fragment group. Thirteen proteins were identified, with different abundance between groups. Identified proteins may be related to embryo-maternal communication, which involves adhesion, nutrition, endothelial cell proliferation, transport, and immunological tolerance. In conclusion, conceptus fragments signalized changes in the protein profile of uterine fluid seven days after ovulation in comparison to the observed at Day 7 in the same cyclic mares.

Kinetics of placenta-specific 8 (PLAC8) in equine placenta during pregnancy and placentitis

Placenta-specific 8 (PLAC8) is one of the placenta-regulatory genes which is highly conserved among eutherian mammals. However, little is known about its expression in equine placenta (chorioallantois; CA and endometrium; EN) during normal and abnormal pregnancy. Therefore, the current study was designed to 1) elucidate the expression of PLAC8 in equine embryonic membranes during the preimplantation period, 2) characterize the expression profile of PLAC8 in equine CA (45d, 4mo, 6mo, 10 mo, 11 mo and postpartum) and EN (14d, 4mo, 6mo, 10 mo, and 11 mo) obtained from pregnant mares (n = 4/timepoint), as well as, d14 non-pregnant EN (n = 4), and 3) investigate the expression profile of PLAC8 in ascending placentitis (n = 5) and in nocardioform placentitis (n = 6) in comparison to normal CA. In the preimplantation period, PLAC8 mRNA was not abundant in the trophectoderm of d8 equine embryo and d14 conceptus, while it was abundant later in d 30, 31, 34, and 45 chorion. In normal pregnancy, PLAC8 mRNA expression in CA at 45 d gradually decline to reach nadir at 6mo before gradually increasing to its peak at 11mo and postpartum CA. The mRNA expression of PLAC8 was significantly upregulated in CA from mares with ascending and nocardioform placentitis compared to control mares. Immunohistochemistry revealed that PLAC8 is localized in equine chorionic epithelium and immune cells. Our results revealed that PLAC8 expression in equine chorion is dynamic during pregnancy and is regulated in an implantation-dependent manner. Moreover, PLAC8 is implicated in the immune response in CA during equine ascending placentitis and nocardioform placentitis.

Loss of POLR1D results in embryonic lethality prior to blastocyst formation in mice

In eukaryotic cells, RNA polymerase (Pol) I and Pol III are dedicated to the synthesis of ribosomal RNA precursors and a variety of small RNAs, respectively. Although RNA Pol I and Pol III complexes are crucial for the regulation of cell growth and cell cycle in all cell types, many of the components of the Pol I and Pol III complexes have not been functionally characterized in mammals. Here, we provide the first in vivo functional characterization of POLR1D, a subunit shared by RNA Pol I and Pol III, during early mammalian embryo development. Our results show that Polr1d mutant embryos cannot be recovered at E7.5 early post-gastrulation stage, suggesting failed implantation. Although Polr1d mutants can be recovered at E3.5, they exhibit delayed/stalled development with morula morphology rather than differentiation into blastocysts. Even with extended time in culture, mutant embryos fail to form blastocysts and eventually die. Analysis of E3.0 embryos revealed severe DNA damage in Polr1d mutants. Additionally, lineage assessment reveals that trophectoderm specification is compromised in the absence of Polr1d. In summary, these findings demonstrate the essential role of POLR1D during early mammalian embryogenesis and highlight cell-lethal phenotype without Polr1d function.

Equine non-invasive time-lapse imaging and blastocyst development

In this study we examined the timeline of mitotic events of invitro-produced equine embryos that progressed to blastocyst stage using non-invasive time-lapse microscopy (TLM). Intracytoplasmic sperm injection (ICSI) embryos were cultured using a self-contained imaging incubator system (Miri®TL; Esco Technologies) that captured brightfield images at 5-min intervals that were then generated into video for retrospective analysis. For all embryos that progressed to the blastocyst stage, the initial event of extrusion of acellular debris preceded all first cleavages and occurred at mean (±s.e.m.) time of 20.0±1.1h after ICSI, whereas 19 of 24 embryos that did not reach the blastocyst stage demonstrated debris extrusion that occurred at 23.8±1.1h, on average 4h longer for this initial premitotic event (P<0.05). Embryos that failed to reach the blastocyst stage demonstrated a 4-h delay compared with those that reached the blastocyst stage to reach the 2-cell stage (P<0.05). All embryos that reached the blastocyst stage expressed pulsation of the blastocyst with visible expansion and contraction at approximate 10-min intervals, or five to six times per hour. Using a logit probability method, we determined that 2- and 8-cell stage embryos could reasonably predict which embryos progressed to the blastocyst stage. Together, the results indicate that TLM for equine embryo development is a dynamic tool with promise for predicting successful embryo development.

A genome-wide scan for candidate lethal variants in Thoroughbred horses

Domestic animal populations are often characterised by high rates of inbreeding and low effective population sizes due to selective breeding practices. These practices can result in otherwise rare recessive deleterious alleles drifting to high frequencies, resulting in reduced fertility rates. This study aimed to identify potential recessive lethal haplotypes in the Thoroughbred horse breed, a closed population that has been selectively bred for racing performance. In this study, we identified a haplotype in the LY49B gene that shows strong evidence of being homozygous lethal, despite having high frequencies of heterozygotes in Thoroughbreds and other domestic horse breeds. Variant analysis of whole-genome sequence data identified two SNPs in the 3'UTR of the LY49B gene that may result in loss of function. Analysis of transcriptomic data from equine embryonic tissue revealed that LY49B is expressed in the trophoblast during placentation stage of development. These findings suggest that LY49B may have an essential, but as yet unknown function in the implantation stage of equine development. Further investigation of this region may allow for the development of a genetic test to improve fertility rates in horse populations. Identification of other lethal variants could assist in improving natural levels of fertility in horse populations.

Common principles of early mammalian embryo self-organisation

Pre-implantation mammalian development unites extreme plasticity with a robust outcome: the formation of a blastocyst, an organised multi-layered structure ready for implantation. The process of blastocyst formation is one of the best-known examples of self-organisation. The first three cell lineages in mammalian development specify and arrange themselves during the morphogenic process based on cell-cell interactions. Despite decades of research, the unifying principles driving early mammalian development are still not fully defined. Here, we discuss the role of physical forces, and molecular and cellular mechanisms, in driving self-organisation and lineage formation that are shared between eutherian mammals.

Identification of Novel lncRNAs Differentially Expressed in Placentas of Chinese Ningqiang Pony and Yili Horse Breeds

As a nutrient sensor, the placenta plays a key role in regulating fetus growth and development. Long non-coding RNAs (lncRNAs) have been shown to regulate growth-related traits. However, the biological function of lncRNAs in horse placentas remains unclear. To compare the expression patterns of lncRNAs in the placentas of the Chinese Ningqiang (NQ) and Yili (YL) breeds, we performed a transcriptome analysis using RNA sequencing (RNA-seq) technology. NQ is a pony breed with an average adult height at the withers of less than 106 cm, whereas that of YL is around 148 cm. Based on 813 million high-quality reads and stringent quality control procedures, 3011 transcripts coding for 1464 placental lncRNAs were identified and mapped to the horse reference genome. We found 107 differentially expressed lncRNAs (DELs) between NQ and YL, including 68 up-regulated and 39 down-regulated DELs in YL. Six (TBX3, CACNA1F, EDN3, KAT5, ZNF281, TMED2, and TGFB1) out of the 233 genes targeted by DELs were identified as being involved in limb development, skeletal myoblast differentiation, and embryo development. Two DELs were predicted to target the TBX3 gene, which was found to be under strong selection and associated with small body size in the Chinese Debao pony breed. This finding suggests the potential functional significance of placental lncRNAs in regulating horse body size.

Characterization of the placental transcriptome through mid to late gestation in the mare

The placenta is a dynamic organ which undergoes extensive remodeling throughout pregnancy to support, protect and nourish the developing fetus. Despite the importance of the placenta, very little is known about its gene expression beyond very early pregnancy and post-partum. Therefore, we utilized RNA-sequencing to characterize the transcriptome from the fetal (chorioallantois) and maternal (endometrium) components of the placenta from mares throughout gestation (4, 6, 10, 11 m). Within the endometrium, 47% of genes changed throughout pregnancy, while in the chorioallantois, 29% of genes underwent significant changes in expression. Further bioinformatic analyses of both differentially expressed genes and highly expressed genes help reveal similarities and differences between tissues. Overall, the tissues were more similar than different, with ~ 95% of genes expressed in both tissues, and high similarities between the most highly expressed genes (9/20 conserved), as well as marked similarities between the PANTHER pathways identified. The most highly expressed genes fell under a few broad categories, including endocrine and immune-related transcripts, iron-binding proteins, extracellular matrix proteins, transport proteins and antioxidants. Serine protease inhibitors were particularly abundant, including SERPINA3, 6 and 14, as well as SPINK7 and 9. This paper also demonstrates the ability to effectively separate maternal and fetal components of the placenta, with only a minimal amount of chorioallantoic contamination in the endometrium (~8%). This aspect of equine placentation is a boon for better understanding gestational physiology and allows the horse to be used in areas where a separation of fetal and maternal tissues is essential. Overall, these data represent the first large-scale characterization of placental gene expression in any species and include time points from multiple mid- to late-gestational stages, helping further our understanding of gestational physiology.

Can trophectoderm RNA analysis predict human blastocyst competency?

A systematic review of the literature showed that trophectoderm biopsy could assist in the selection of healthy embryos for uterine transfer without affecting implantation rates. However, previous studies attempting to establish the relationship between trophectoderm gene expression profiles and implantation competency using either microarrays or RNA sequencing strategies, were not sufficiently optimized to handle the exceptionally low RNA inputs available from biopsied material. In this pilot study, we report that differential gene expression in human trophectoderm biopsies assayed by an ultra-sensitive next generation RNA sequencing strategy could predict blastocyst implantation competence. RNA expression profiles from isolated human trophectoderm cells were analysed with established clinical pregnancy being the primary endpoint. Following RNA sequencing, a total of 47 transcripts were found to be significantly differentially expressed between the trophectoderm cells from successfully implanted (competent) versus unsuccessful (incompetent) blastocysts. Of these, 36 transcripts were significantly down-regulated in the incompetent blastocysts, including Hydroxysteroid 17-Beta Dehydrogenase 1 (HSD17B1) and Cytochrome P450 Family 11 Subfamily A Member 1 (CYP11A1), while the remaining 11 transcripts were significantly up-regulated, including BCL2 Antagonist/Killer 1 (BAK1) and KH Domain Containing 1 Pseudogene 1 (KHDC1P1) of which the latter was always detected in the incompetent and absent in all competent blastocysts. Ontological analysis of differentially expressed RNAs revealed pathways involved in steroidogenic processes with high confidence. Novel differentially expressed transcripts were also noted by reference to a de novo sequence assembly. The selection of the blastocyst with the best potential to support full-term pregnancy following single embryo transfer could reduce the need for multiple treatment cycles and embryo transfers. The main limitation was the low sample size (N = 8). Despite this shortcoming, the pilot suggests that trophectoderm biopsy could assist with the selection of healthy embryos for embryo transfer. A larger cohort of samples is needed to confirm these findings. Abbreviations: AMA: advanced maternal age; ART: assisted reproductive technology; CP: clinical pregnancy; DE: differential expression; FDR: false discovery rate; IVF: in vitro fertilization; LD PCR: long distance PCR; qRT-PCR: quantitative real-time PCR; SET: single embryo transfer; TE: trophectoderm.

Transcriptome analysis reveals that fertilization with cryopreserved sperm downregulates genes relevant for early embryo development in the horse

Artificial insemination with cryopreserved spermatozoa is a major assisted reproductive technology in many species. In horses, as in humans, insemination with cryopreserved sperm is associated with lower pregnancy rates than those for fresh sperm, however, direct effects of sperm cryopreservation on the development of resulting embryos are largely unexplored. The aim of this study was to investigate differences in gene expression between embryos resulting from fertilization with fresh or cryopreserved sperm. Embryos were obtained at 8, 10 or 12 days after ovulation from mares inseminated post-ovulation on successive cycles with either fresh sperm or frozen-thawed sperm from the same stallion, providing matched embryo pairs at each day. RNA was isolated from two matched pairs (4 embryos) for each day, and cDNA libraries were built and sequenced. Significant differences in transcripts per kilobase million (TPM) were determined using (i) genes for which the expression difference between treatments was higher than 99% of that in the random case (P < 0.01), and (ii) genes for which the fold change was ≥ 2, to avoid expression bias in selection of the candidate genes. Molecular pathways were explored using the DAVID webserver, followed by network analyses using STRING, with a threshold of 0.700 for positive interactions. The transcriptional profile of embryos obtained with frozen-thawed sperm differed significantly from that for embryos derived from fresh sperm on all days, showing significant down-regulation of genes involved in biological pathways related to oxidative phosphorylation, DNA binding, DNA replication, and immune response. Many genes with reduced expression were orthologs of genes known to be embryonic lethal in mice. This study, for the first time, provides evidence of altered transcription in embryos resulting from fertilization with cryopreserved spermatozoa in any species. As sperm cryopreservation is commonly used in many species, including human, the effect of this intervention on expression of developmentally important genes in resulting embryos warrants attention.

Growth factor modulation of equine trophoblast mitosis and prostaglandin gene expression

 To provide insight into maternal recognition of pregnancy control in equids, the mitogenic and developmental effects of endometrium-expressed growth factors (epidermal growth factor (EGF), fibroblast growth factor 2 (FGF2), hepatocyte growth factor (HGF), and insulin-like growth factor-1 (IGF-1)) were examined in equine iTr cells, an equine trophectoderm cell line. Initial western blots revealed that HGF and IGF-1 stimulate phosphorylation of AKT serine/threonine kinase 1 (AKT1) and EGF, FGF2, or HGF resulted in phosphorylation of both extracellular signal-regulated kinase 1 (ERK1) and ERK2. Mitotic activity was stimulated (P < 0.05) by EGF, FGF2, and HGF. Chemical disruption of mitogen-activated protein kinase kinases 1 and 2 (MEK1/2) phosphorylation suppresses (P < 0.05) proliferation in control and growth factor treated cells demonstrating a dependence on ERK1/2 for mitotic activity. Treatment of iTr cells with EGF or HGF in the presence of an AKT1 inhibitor prohibits (P < 0.05) growth factor stimulated proliferation. The effect of EGF, FGF2, HGF, and IGF-1 on steroid biosynthetic enzyme gene expression, including prostaglandin-endoperoxide synthase 2 (PTGS2), was determined by real-time PCR. Neither EGF, FGF2, nor IGF-1 affected PTGS2 expression while HGF caused a two-fold increase (P < 0.05) in expression. Co-supplementation with HGF and an AKT1 inhibitor did not block PTGS2 expression, whereas providing an MEK1/2 inhibitor prevented (P < 0.05) the HGF-mediated increase in PTGS2. These results provide novel evidence of a role for HGF in equine trophectoderm proliferation and prostaglandin biosynthesis.

Targeting αGal epitopes for multi-species embryo immunosurgery

Immunosurgical isolation of the inner cell mass (ICM) from blastocysts is based on complement-mediated lysis of antibody-coated trophectoderm (TE) cells. Conventionally, anti-species antisera, containing antibodies against multiple undefined TE-cell epitopes, have been used as the antibody source. We previously generated α-1,3-galactosyltransferase deficient (GTKO) pigs to prevent hyperacute rejection of pig-to-primate xenotransplants. Since GTKO pigs lack galactosyl-α-1,3-galactose (αGal) but are exposed to this antigen (e.g. αGal on gut bacteria), they produce anti-αGal antibodies. In this study, we examined whether serum from GTKO pigs could be used as a novel antibody source for multi-species embryo immunosurgery. Mouse, rabbit, pig and cattle blastocysts were used for the experiment. Expression of αGal epitopes on the surface of TE cells was detected in blastocysts of all species tested. GTKO pig serum contained sufficient anti-αGal antibodies to induce complement-mediated lysis of TE cells in blastocysts from all species investigated. Intact ICMs could be successfully recovered and the majority showed the desired level of purity. Our study demonstrates that GTKO pig serum is a reliable and effective source of antibodies targeting the αGal epitopes of TE cells for multi-species embryo immunosurgery.

Dynamic changes in gene expression and signalling during trophoblast development in the horse

Equine chorionic girdle trophoblast cells play important endocrine and immune functions critical in supporting pregnancy. Very little is known about the genes and pathways that regulate chorionic girdle trophoblast development. Our aim was to identify genes and signalling pathways active in vivo in equine chorionic girdle trophoblast within a critical 7-days window. We exploited the late implantation of the equine conceptus to obtain trophoblast tissue. An Agilent equine 44K microarray was performed using RNA extracted from chorionic girdle and chorion (control) from equine pregnancy days 27, 30, 31 and 34 (n = 5), corresponding to the initiation of chorionic girdle trophoblast proliferation, differentiation and migration. Data were analysed using R packages limma and maSigPro, Ingenuity Pathway Analysis and DAVID and verified using qRT-PCR, promoter analysis, western blotting and migration assays. Microarray analysis showed gene expression (absolute log FC >2, FDR-adjusted P < 0.05) was rapidly and specifically induced in the chorionic girdle between days 27 and 34 (compared to day 27, day 30 = 116, day 31 = 317, day 34 = 781 genes). Pathway analysis identified 35 pathways modulated during chorionic girdle development (e.g. FGF, integrin, Rho GTPases, MAPK) including pathways that have limited description in mammalian trophoblast (e.g. IL-9, CD40 and CD28 signalling). Rho A and ERK/MAPK activity was confirmed as was a role for transcription factor ELF5 in regulation of the CGB promoter. The purity and accessibility of chorionic girdle trophoblast proved to be a powerful resource to identify candidate genes and pathways involved in early equine placental development.

Pre-implantation Development of Domestic Animals

During the first days following fertilization, cells of mammalian embryo gradually lose totipotency, acquiring distinct identity. The first three lineages specified in the mammalian embryo are pluripotent epiblast, which later gives rise to the embryo proper, and two extraembryonic lineages, hypoblast (also known as primitive endoderm) and trophectoderm, which form tissues supporting development of the fetus in utero. Most of our knowledge regarding the mechanisms of early lineage specification in mammals comes from studies in the mouse. However, the growing body of evidence points to both similarities and species-specific differences. Understanding molecular and cellular mechanisms of early embryonic development in nonrodent mammals expands our understanding of basic mechanisms of differentiation and is essential for the development of effective protocols for assisted reproduction in agriculture, veterinary medicine, and for biomedical research. This review summarizes the current state of knowledge on key events in epiblast, hypoblast, and trophoblast differentiation in domestic mammals.

Characterization of companion animal pluripotent stem cells

Pluripotent stem cells have the capacity to grow indefinitely in culture and differentiate into derivatives of the three germ layers. These properties underpin their potential to be used in regenerative medicine. Originally derived from early embryos, pluripotent stem cells can now be derived by reprogramming an adult cell back to a pluripotent state. Companion animals such as horses, dogs, and cats suffer from many injuries and diseases for which regenerative medicine may offer new treatments. As many of the injuries and diseases are similar to conditions in humans the use of companion animals for the experimental and clinical testing of stem cell and regenerative medicine products would provide relevant animal models for the translation of therapies to the human field. In order to fully utilize companion animal pluripotent stem cells robust, standardized methods of characterization must be developed to ensure that safe and effective treatments can be delivered. In this review we discuss the methods that are available for characterizing pluripotent stem cells and the techniques that have been applied in cells from companion animals. We describe characteristics which have been described consistently across reports as well as highlighting discrepant results. Significant steps have been made to define the in vitro culture requirements and drive lineage specific differentiation of pluripotent stem cells in companion animal species. However, additional basic research to compare pluripotent stem cell types and define characteristics of pluripotency in companion animal species is still required. © 2017 International Society for Advancement of Cytometry.

Tissue organization alters gene expression in equine induced trophectoderm cells

Rapid morphological and gene expression changes occur during the early formation of a mammalian blastocyst. Critical to successful retention of the blastocyst and pregnancy is a functional trophectoderm (TE) that supplies the developing embryo with paracrine factors and hormones. The contribution of TE conformational changes to gene expression was examined in equine induced trophoblast (iTr) cells. Equine iTr cells were cultured as monolayers or in suspension to form spheres. The spheres are hollow and structurally reminiscent of native equine blastocysts. Total RNA was isolated from iTr monolayers and spheres and analyzed by RNA sequencing. An average of 32.2 and 31million aligned reads were analyzed for the spheres and monolayers, respectively. Forty-four genes were unique to monolayers and 45 genes were expressed only in spheres. Conformation did not affect expression of CDX2, POU5F1, TEAD4, ETS2, ELF3, GATA2 or TFAP2A, the core gene network of native TE. Bioinformatic analysis was used to identify classes of genes differentially expressed in response to changes in tissue shape. In both iTr spheres and monolayers, the majority of the differentially expressed genes were associated with binding activity in cellular, developmental and metabolic processes. Inherent to protein:protein interactions, several receptor-ligand families were identified in iTr cells with enrichment of genes coding for PI3-kinase and MAPK signaling intermediates. Our results provide evidence for ligand initiated kinase signaling pathways that underlie early trophectoderm structural changes.

Patients with Dilated Cardiomyopathy and Sustained Monomorphic Ventricular Tachycardia Show Up-Regulation of KCNN3 and KCNJ2 Genes and CACNG8-Linked Left Ventricular Dysfunction

AIMS

Disruptions in cardiac ion channels have shown to influence the impaired cardiac contraction in heart failure. We sought to determine the altered gene expression profile of this category in dilated cardiomyopathy (DCM) patients and relate the altered gene expression with the clinical signs present in our patients, such as ventricular dysfunction and sustained monomorphic ventricular tachycardia (SMVT).

METHODS AND RESULTS

Left ventricular (LV) tissue samples were used in RNA-sequencing technique to elucidate the transcriptomic changes of 13 DCM patients compared to controls (n = 10). We analyzed the differential gene expression of cardiac ion channels, and we found a total of 34 altered genes. We found that the calcium channel CACNG8 mRNA and protein levels were down-regulated and highly and inversely related with LV ejection fraction (LVEF) (r = -0.78, P<0.01). Furthermore, the potassium channels KCNN3 and KCNJ2 mRNA and protein levels were up-regulated and showed also a significant and inverse correlation with LVEF (r = -0.61, P<0.05; r = -0.60, P<0.05) in patients with SMVT.

CONCLUSION

A broad set of deregulated genes have been identified by RNA-sequencing technique. The relationship of CACNG8, KCNN3 and KCNJ2 with LVEF, and the up-regulation of KCNN3 and KCNJ2 in all patients with SMVT, irrespective of CACNG8 expression, suggest a significant role for these three ion flux related genes in the LV dysfunction present in this cardiomyopathy and an important relationship between KCNN3 and KCNJ2 up-regulation and the presence of SMVT.

RNA sequencing as a powerful tool in searching for genes influencing health and performance traits of horses

RNA sequencing (RNA-seq) by next-generation technology is a powerful tool which creates new possibilities in whole-transcriptome analysis. In recent years, with the use of the RNA-seq method, several studies expanded transcriptional gene profiles to understand interactions between genotype and phenotype, supremely contributing to the field of equine biology. To date, in horses, massive parallel sequencing of cDNA has been successfully used to identify and quantify mRNA levels in several normal tissues, as well as to annotate genes. Moreover, the RNA-seq method has been applied to identify the genetic basis of several diseases or to investigate organism adaptation processes to the training conditions. The use of the RNA-seq approach has also confirmed that horses can be useful as a large animal model for human disease, especially in the field of immune response. The presented review summarizes the achievements of profiling gene expression in horses (Equus caballus).

Generation of a de novo transcriptome from equine lamellar tissue

BACKGROUND

Laminitis, the structural failure of interdigitated tissue that suspends the distal skeleton within the hoof capsule, is a devastating disease that is the second leading cause of both lameness and euthanasia in the horse. Current transcriptomic research focuses on the expression of known genes. However, as this tissue is quite unique and equine gene annotation is largely derived from computational predictions, there are likely yet uncharacterized transcripts that may be involved in the etiology of laminitis. In order to create a novel annotation resource, we performed whole transcriptome sequencing of sagittal lamellar sections from one control and two laminitis affected horses.

RESULTS

Whole transcriptome sequencing of the three samples resulted in 113 million reads. Overall, 88 % of the reads mapped to the equCab2 reference genome, allowing for the identification of 119,430 SNPs. The de novo assembly generated around 75,000 transcripts, of which 36,000 corresponded to known annotations. Annotated transcript models are hosted in a public data repository and thus can be easily accessed or loaded into genome browsers. RT-PCR of 12 selected assemblies confirmed structure and expression in lamellar tissue.

CONCLUSIONS

Transcriptome sequencing represents a powerful tool to expand on equine annotation and identify novel targets for further laminitis research.

Genetic Susceptibility to Rhodococcus equi

Rhodococcus equi pneumonia is a major cause of morbidity and mortality in neonatal foals. Much effort has been made to identify preventative measures and new treatments for R. equi with limited success. With a growing focus in the medical community on understanding the genetic basis of disease susceptibility, investigators have begun to evaluate the interaction of the genetics of the foal with R. equi. This review describes past efforts to understand the genetic basis underlying R. equi susceptibility and tolerance. It also highlights the genetic technology available to study horses and describes the use of this technology in investigating R. equi. This review provides readers with a foundational understanding of candidate gene approaches, single nucleotide polymorphism‐based, and copy number variant‐based genome‐wide association studies, and next generation sequencing (both DNA and RNA).

The transcriptome of equine peripheral blood mononuclear cells

Complete transcriptomic data at high resolution are available only for a few model organisms with medical importance. The gene structures of non-model organisms are mostly computationally predicted based on comparative genomics with other species. As a result, more than half of the horse gene models are known only by projection. Experimental data supporting these gene models are scarce. Moreover, most of the annotated equine genes are single-transcript genes. Utilizing RNA sequencing (RNA-seq) the experimental validation of predicted transcriptomes has become accessible at reasonable costs. To improve the horse genome annotation we performed RNA-seq on 561 samples of peripheral blood mononuclear cells (PBMCs) derived from 85 Warmblood horses. The mapped sequencing reads were used to build a new transcriptome assembly. The new assembly revealed many alternative isoforms associated to known genes or to those predicted by the Ensembl and/or Gnomon pipelines. We also identified 7,531 transcripts not associated with any horse gene annotated in public databases. Of these, 3,280 transcripts did not have a homologous match to any sequence deposited in the NCBI EST database suggesting horse specificity. The unknown transcripts were categorized as coding and noncoding based on predicted coding potential scores. Among them 230 transcripts had high coding potential score, at least 2 exons, and an open reading frame of at least 300 nt. We experimentally validated 9 new equine coding transcripts using RT-PCR and Sanger sequencing. Our results provide valuable detailed information on many transcripts yet to be annotated in the horse genome.

DPPA3 prevents cytosine hydroxymethylation of the maternal pronucleus and is required for normal development in bovine embryos

Dppa3 has been described in mice as an important maternal factor contributed by the oocyte that participates in protecting the maternal genome from oxidation of methylated cytosines (5mC) to hydroxymethylated cytosines (5hmC). Dppa3 is also required for normal mouse preimplantation development. This gene is poorly conserved across mammalian species, with less than 32% of protein sequence shared between mouse, cow and human. RNA-seq analysis of bovine oocytes and preimplantation embryos revealed that DPPA3 transcripts are some of the most highly abundant mRNAs in the oocyte, and their levels gradually decrease toward the time of embryonic genome activation (EGA). Knockdown of DPPA3 by injection of siRNA in germinal vesicle (GV) stage oocytes was used to assess its role in epigenetic remodeling and embryo development. DPPA3 knockdown resulted in increased intensity of 5hmC staining in the maternal pronucleus (PN), demonstrating a role for this factor in the asymmetric remodeling of the maternal and paternal PN in bovine zygotes. Also, DPPA3 knockdown decreased the developmental competence of parthenogenetic and in vitro fertilized embryos. Finally, DPPA3 knockdown embryos that reached the blastocyst stage had significantly fewer ICM cells as compared with control embryos. We conclude that DPPA3 is a maternal factor important for correct epigenetic remodeling and normal embryonic development in cattle, indicating that the role of DPPA3 during early development is conserved between species.

Characterisation of the horse transcriptome from immunologically active tissues

The immune system of the horse has not been well studied, despite the fact that the horse displays several features such as sensitivity to bacterial lipopolysaccharide that make them in many ways a more suitable model of some human disorders than the current rodent models. The difficulty of working with large animal models has however limited characterisation of gene expression in the horse immune system with current annotations for the equine genome restricted to predictions from other mammals and the few described horse proteins. This paper outlines sequencing of 184 million transcriptome short reads from immunologically active tissues of three horses including the genome reference “Twilight”. In a comparison with the Ensembl horse genome annotation, we found 8,763 potentially novel isoforms.