Explore the vast range of titles available.

MicroRNAs (miRNAs) are essential regulators of gene expression, yet few comprehensive databases exist for miRNA expression in non-model species, limiting our ability to characterize their roles in gene regulation, development, and disease. Similarly, isomiRs - length and sequence isoforms of canonical miRNAs with potentially altered regulatory targets and functions - have received even less attention in non-model species, including the horse, leaving a critical gap in our understanding of their biological significance. To address these challenges, we developed an open-source, containerized pipeline for identifying and quantifying miRNAs and isomiRs (FARmiR: Framework for Analysis and Refinement of miRNAs), and an associated interactive browser (AIMEE: Animal IsomiR and MiRNA Expression Explorer). AIMEE was developed to make miRNA expression data more accessible and user-friendly, a feature often lacking from other expression atlases. These tools were developed using equine data but can be readily extended to other species. Using these tools, we aggregated 461 small RNA-seq datasets, spanning 61 distinct tissues, integrating data from public repositories, an American Quarter Horse cohort, and the Functional Annotation of ANimal Genome (FAANG) consortium Thoroughbred samples, predicting 5,781 miRNAs and isomiRs. This work represents the largest systematically curated atlas of equine miRNA expression to date, providing a valuable resource that will enhance our understanding of miRNA and isomiR functions in tissue-specific regulation and ultimately improve biomarker discovery, functional genomics, and precision veterinary medicine.

Aberrant microRNA expression contributes to breast cancer progression and endocrine resistance. We reported that although tamoxifen stimulated miR-29b-1/a transcription in tamoxifen (TAM)-resistant breast cancer cells, ectopic expression of miR-29b-1/a did not drive TAM-resistance in MCF-7 breast cancer cells. However, miR-29b-1/a overexpression significantly repressed TAM-resistant LCC9 cell proliferation, suggesting that miR-29b-1/a is not mediating TAM resistance but acts as a tumor suppressor in TAM-resistant cells. The target genes mediating this tumor suppressor activity were unknown. Here, we identify miR-29b-1 and miR-29a target transcripts in both MCF-7 and LCC9 cells. We find that miR-29b-1 and miR-29a regulate common and unique transcripts in each cell line. The cell-specific and common downregulated genes were characterized using the MetaCore Gene Ontology (GO) enrichment analysis algorithm. LCC9-sepecific miR-29b-1/a-regulated GO processes include oxidative phosphorylation, ATP metabolism, and apoptosis. Extracellular flux analysis of cells transfected with anti- or pre- miR-29a confirmed that miR-29a inhibits mitochondrial bioenergetics in LCC9 cells. qPCR,luciferase reporter assays, and western blot also verified the ATP synthase subunit genes ATP5G1 and ATPIF1 as bone fide miR29b-1/a targets. Our results suggest that miR-29 repression of TAM-resistant breast cancer cell proliferation is mediated in part through repression of genes important in mitochondrial bioenergetics.

Centromeres are epigenetically specified by the histone H3 variant CENP-A. Although mammalian centromeres are typically associated with satellite DNA, we previously demonstrated that the centromere of horse chromosome 11 (ECA11) is completely devoid of satellite DNA. We also showed that the localization of its CENP-A binding domain is not fixed but slides within an about 500 kb region in different individuals, giving rise to positional alleles. These epialleles are inherited as Mendelian traits but their position can move in one generation. It is still unknown whether centromere sliding occurs during meiosis or during development. Here, we first improve the sequence of the ECA11 centromeric region in the EquCab3.0 assembly. Then, to test whether centromere sliding may occur during development, we map the CENP-A binding domains of ECA11 using ChIP-seq in five tissues of different embryonic origin from the four horses of the equine FAANG (Functional Annotation of ANimal Genomes) consortium. Our results demonstrate that the centromere is localized in the same region in all tissues, suggesting that the position of the centromeric domain is maintained during development. Sequence improvement of the centromeric region of horse chromosome 11 and mapping of the histone H3 variant CENP-A binding domains show that the centromere is localized in the same region in horse tissues from different embryonic origin.

The genomic sequence of the horse has been available since 2009, providing critical resources for discovering important genomic variants regarding both animal health and population structures. However, to fully understand the functional implications of these variants, detailed annotation of the horse genome is required. Due to the limited availability of functional data for the equine genome, as well as the technical limitations of short-read RNA-seq, existing annotation of the equine genome contains limited information about important aspects of gene regulation, such as alternate isoforms and regulatory elements, which are either not transcribed or transcribed at a very low level. To solve above problems, the Functional Annotation of the Animal Genomes (FAANG) project proposed a systemic approach to tissue collection, phenotyping, and data generation, adopting the blueprint laid out by the Encyclopedia of DNA Elements (ENCODE) project. Here we detail the first comprehensive overview of gene expression and regulation in the horse, presenting 39,625 novel transcripts, 84,613 candidate cis-regulatory elements (CRE) and their target genes, 332,115 open chromatin regions genome wide across a diverse set of tissues. We showed substantial concordance between chromatin accessibility, chromatin states in different genic features and gene expression. This comprehensive and expanded set of genomics resources will provide the equine research community ample opportunities for studies of complex traits in the horse.

Whole genome sequences (WGS) of 185 North American Thoroughbred horses were compared to quantify the number and frequency of variants, diversity of mitotypes, and autosomal runs of homozygosity (ROH). Of the samples, 82 horses were born between 1965 and 1986 (Group 1); the remaining 103, selected to maximize pedigree diversity, were born between 2000 and 2020 (Group 2). Over 14.3 million autosomal variants were identified with 4.5–5.0 million found per horse. Mitochondrial sequences associated the North American Thoroughbreds with 9 of 17 clades previously identified among diverse breeds. Individual coefficients of inbreeding, estimated from ROH, averaged 0.266 (Group 1) and 0.283 (Group 2). When SNP arrays were simulated using subsets of WGS markers, the arrays over-estimated lengths of ROH. WGS-based estimates of inbreeding were highly correlated (r > 0.98) with SNP array-based estimates, but only moderately correlated (r = 0.40) with inbreeding based on 5-generation pedigrees. On average, Group 1 horses had more heterozygous variants ( P  < 0.001), more total variants ( P  < 0.001), and lower individual inbreeding (F ROH ; P  < 0.001) than horses in Group 2. However, the distribution of numbers of variants, allele frequency, and extent of ROH overlapped among all horses such that it was not possible to identify the group of origin of any single horse using these measures. Consequently, the Thoroughbred population would be better monitored by investigating changes in specific variants, rather than relying on broad measures of diversity. The WGS for these 185 horses is publicly available for comparison to other populations and as a foundation for modeling changes in population structure, breeding practices, or the appearance of deleterious variants.

Significance The domestication of the horse revolutionized warfare, trade, and the exchange of people and ideas. This at least 5,500-y-long process, which ultimately transformed wild horses into the hundreds of breeds living today, is difficult to reconstruct from archeological data and modern genetics alone. We therefore sequenced two complete horse genomes, predating domestication by thousands of years, to characterize the genetic footprint of domestication. These ancient genomes reveal predomestic population structure and a significant fraction of genetic variation shared with the domestic breeds but absent from Przewalski’s horses. We find positive selection on genes involved in various aspects of locomotion, physiology, and cognition. Finally, we show that modern horse genomes contain an excess of deleterious mutations, likely representing the genetic cost of domestication. The domestication of the horse ∼5.5 kya and the emergence of mounted riding, chariotry, and cavalry dramatically transformed human civilization. However, the genetics underlying horse domestication are difficult to reconstruct, given the near extinction of wild horses. We therefore sequenced two ancient horse genomes from Taymyr, Russia (at 7.4- and 24.3-fold coverage), both predating the earliest archeological evidence of domestication. We compared these genomes with genomes of domesticated horses and the wild Przewalski’s horse and found genetic structure within Eurasia in the Late Pleistocene, with the ancient population contributing significantly to the genetic variation of domesticated breeds. We furthermore identified a conservative set of 125 potential domestication targets using four complementary scans for genes that have undergone positive selection. One group of genes is involved in muscular and limb development, articular junctions, and the cardiac system, and may represent physiological adaptations to human utilization. A second group consists of genes with cognitive functions, including social behavior, learning capabilities, fear response, and agreeableness, which may have been key for taming horses. We also found that domestication is associated with inbreeding and an excess of deleterious mutations. This genetic load is in line with the “cost of domestication” hypothesis also reported for rice, tomatoes, and dogs, and it is generally attributed to the relaxation of purifying selection resulting from the strong demographic bottlenecks accompanying domestication. Our work demonstrates the power of ancient genomes to reconstruct the complex genetic changes that transformed wild animals into their domesticated forms, and the population context in which this process took place.

Background To date, genome-scale analyses in the domestic horse have been limited by suboptimal single nucleotide polymorphism (SNP) density and uneven genomic coverage of the current SNP genotyping arrays. The recent availability of whole genome sequences has created the opportunity to develop a next generation, high-density equine SNP array. Results Using whole genome sequence from 153 individuals representing 24 distinct breeds collated by the equine genomics community, we cataloged over 23 million de novo discovered genetic variants. Leveraging genotype data from individuals with both whole genome sequence, and genotypes from lower-density, legacy SNP arrays, a subset of ~5 million high-quality, high-density array candidate SNPs were selected based on breed representation and uniform spacing across the genome. Considering probe design recommendations from a commercial vendor (Affymetrix, now Thermo Fisher Scientific) a set of ~2 million SNPs were selected for a next-generation high-density SNP chip (MNEc2M). Genotype data were generated using the MNEc2M array from a cohort of 332 horses from 20 breeds and a lower-density array, consisting of ~670 thousand SNPs (MNEc670k), was designed for genotype imputation. Conclusions Here, we document the steps taken to design both the MNEc2M and MNEc670k arrays, report genomic and technical properties of these genotyping platforms, and demonstrate the imputation capabilities of these tools for the domestic horse.

The domestic yak (Bos grunniens) has experienced introgression with domestic cattle in its native Qinghai‐Tibetan Plateau and after introduction to North America (NA), although the extent to which the latter has occurred is not well documented. We used complete mitochondrial (mtDNA) sequences of 12 NA yak and aligned them to the B. taurus reference genome for annotation. Identified variation among the NA haplotypes included a total of 982 variants, of which 99 were nonsynonymous single nucleotide polymorphisms. Among the NA yak, we identified nine unique mitotypes, which a haplotype network separated into two distinct clusters. A maximum likelihood tree including 86 Bovidae taxa revealed six NA yak haplotypes formed a clade with B. indicus; the other three haplotypes grouped with B. grunniens and fell as a sister clade to bison, gaur, and gayal. These data demonstrate two mitochondrial origins of NA yak, likely dating prior to their importation to NA. We isolated satellite cells from seven yak that represented both major mitotypes (B. indicus [N = 4] or yak [N = 3]) to investigate possible differences in ATP production. Oxidative consumption rates and extracellular acidification rates were quantified as measures of mitochondrial respiration and glycolysis using the Seahorse ATP Rate Assay. Cells with the B. grunniens mitotype had less total energy metabolism (p = 0.016). This difference may reflect adaptations to ancestral environments and selective pressures associated with husbandry practices and breeding. Little is known about the origin of yak that were imported to North America. Using complete mtDNA sequence, we identified two distinct lineages in North American yak; one of these appears to be the consequence of past introgression from indicine cattle. Variation in the mtDNA included many putatively functional variants. In a pilot study of mitochondrial function, we observed a difference in ATP production of satellite cells associated with the mtDNA haplotype.

Equine arteritis virus (EAV) is the causative agent of equine viral arteritis (EVA), a respiratory, systemic, and reproductive disease of horses and other equid species. Following natural infection, 10-70% of the infected stallions can become persistently infected and continue to shed EAV in their semen for periods ranging from several months to life. Recently, we reported that some stallions possess a subpopulation(s) of CD3+ T lymphocytes that are susceptible to in vitro EAV infection and that this phenotypic trait is associated with long-term carrier status following exposure to the virus. In contrast, stallions not possessing the CD3+ T lymphocyte susceptible phenotype are at less risk of becoming long-term virus carriers. A genome wide association study (GWAS) using the Illumina Equine SNP50 chip revealed that the ability of EAV to infect CD3+ T lymphocytes and establish long-term carrier status in stallions correlated with a region within equine chromosome 11. Here we identified the gene and mutations responsible for these phenotypes. Specifically, the work implicated three allelic variants of the equine orthologue of CXCL16 (EqCXCL16) that differ by four non-synonymous nucleotide substitutions (XM_00154756; c.715 A → T, c.801 G → C, c.804 T → A/G, c.810 G → A) within exon 1. This resulted in four amino acid changes with EqCXCL16S (XP_001504806.1) having Phe, His, Ile and Lys as compared to EqCXL16R having Tyr, Asp, Phe, and Glu at 40, 49, 50, and 52, respectively. Two alleles (EqCXCL16Sa, EqCXCL16Sb) encoded identical protein products that correlated strongly with long-term EAV persistence in stallions (P<0.000001) and are required for in vitro CD3+ T lymphocyte susceptibility to EAV infection. The third (EqCXCL16R) was associated with in vitro CD3+ T lymphocyte resistance to EAV infection and a significantly lower probability for establishment of the long-term carrier state (viral persistence) in the male reproductive tract. EqCXCL16Sa and EqCXCL16Sb exert a dominant mode of inheritance. Most importantly, the protein isoform EqCXCL16S but not EqCXCL16R can function as an EAV cellular receptor. Although both molecules have equal chemoattractant potential, EqCXCL16S has significantly higher scavenger receptor and adhesion properties compared to EqCXCL16R.

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.