Pacific lion-paw scallops were collected from natural aggregations in Laguna Ojo de Liebre (Pacific Ocean), the Gulf of California, and from aquaculture facilities for genetic diversity analyses. Mitochondrial DNA sequencing uncovered two highly supported clades separated by 2.5% divergence. Data from ten microsatellite markers suggest individuals from these mitogroups are introgressed, raising questions about the mitotype origin. Some evidence suggests gene flow between La Paz and Ojo de Liebre; otherwise the Gulf of California and Ojo de Liebre are acting as two distinct populations. It is unclear whether translocations between sites have influenced the observed genetic structure or whether gene flow has been facilitated by past geologic events. Finally, scallops spawned for aquaculture are unique from the wild and have significantly less diversity. These results warrant the attention of managers and producers who should work to monitor and conserve genetic diversity in both wild and aquaculture populations.

We characterized 37 new microsatellite markers in the Pacific lion-paw scallop (Nodipecten subnodosus) and tested for cross-amplification in four other species. Genetic diversity was estimated using 24 individuals from the Lagoon Ojo de Liebre, B.C.S., Mexico. Allelic richness varied from 5 to 27 alleles per locus and the average expected heterozygosity was 0.76. Ten loci exhibited significant departure from Hardy–Weinberg equilibrium likely due to the presence of null alleles. Sixteen of these markers cross-amplified in closely related N. nodosus, while little or no amplification was observed in three Argopecten species.

An assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq) has become an increasingly popular method to assess genome-wide chromatin accessibility in isolated nuclei from fresh tissues. However, many biobanks contain only snap-frozen tissue samples. While ATAC-seq has been applied to frozen brain tissues in human, its applicability in a wide variety of tissues in horse remains unclear. The Functional Annotation of Animal Genome (FAANG) project is an international collaboration aimed to provide high quality functional annotation of animal genomes. The equine FAANG initiative has generated a biobank of over 80 tissues from two reference female animals and experiments to begin to characterize tissue specificity of genome function for prioritized tissues have been performed. Due to the logistics of tissue collection and storage, extracting nuclei from a large number of tissues for ATAC-seq at the time of collection is not always practical. To assess the feasibility of using stored frozen tissues for ATAC-seq and to provide a guideline for the equine FAANG project, we compared ATAC-seq results from nuclei isolated from frozen tissue to cryopreserved nuclei (CN) isolated at the time of tissue harvest in liver, a highly cellular homogenous tissue, and lamina, a relatively acellular tissue unique to the horse. We identified 20,000-33,000 accessible chromatin regions in lamina and 22-61,000 in liver, with consistently more peaks identified using CN isolated at time of tissue collection. Our results suggest that frozen tissues are an acceptable substitute when CN are not available. For more challenging tissues such as lamina, nuclei extraction at the time of tissue collection is still preferred for optimal results. Therefore, tissue type and accessibility to intact nuclei should be considered when designing ATAC-seq experiments.

Background

Genetic tests for variants in MYOT (P2; rs1138656462), FLNC (P3a; rs1139799323 or P3b; rs1142918816) and MYOZ3 (P4; rs1142544043) genes are offered commercially to diagnose myofibrillar myopathy (MFM) and type 2 polysaccharide storage myopathy (PSSM2) in Quarter Horses (QH).

Objectives

To determine if PSSM2-QH has histopathological features of MFM. To compare genotype and allele frequencies of variants P2, P3, P4 between control-QH and PSSM2-QH diagnosed by histopathology.

Study design

Retrospective cross-sectional.

Methods

The study includes a total of 229 healthy control-QH, 163 PSSM2-QH GYS1 mutation negative. Desmin stains of gluteal/semimembranosus muscle were evaluated. Purported disease alleles P2, P3a, P3b, P4 were genotyped by pyrosequencing. Genotype, allele frequency and total number of variant alleles or loci were compared between phenotypes using additive/genotypic and dominant models and quantitative effects evaluated by multivariable logistic regression.

Results

Histopathological features of MFM were absent in all QH. A P variant allele at any locus was not associated (P > .05) with a histopathological diagnosis of PSSM2 and one or more P variants were common in control-QH (57%) and PSSM2-QH (61%). Allele frequencies (control/PSSM2) were: 0.24/0.21 (P2), 0.07/0.12 (P3a), 0.07/0.11 (P3b) and 0.06/0.08 (P4). P3a and P3b loci were not independent (r²  = 0.894); and not associated with PSSM2 histopathology comparing the haplotype of both P3a and P3b variants to other haplotypes. A receiver operator curve did not accurately predict the PSSM2 phenotype (AUC = 0.67, 95% CI 0.62-0.72), and there was no difference in the total number of variant loci or total variant allele count between control-QH and PSSM2-QH.

Main limitations

P3a and P3b were not in complete linkage disequilibrium.

Conclusions

The P2, P3 and P4 variants in genes associated with human MFM were not associated with PSSM2 in 392 QH. Their use would improperly diagnose PSSM2/MFM in 57% of healthy QH and fail to diagnose PSSM2 in 40% of QH with histopathological evidence of PSSM2.

Background

Myofibrillar myopathy (MFM) of unknown aetiology has recently been identified in Warmblood (WB) horses. In humans, 16 genes have been implicated in various MFM-like disorders.

Objectives

To identify variants in 16 MFM candidate genes and compare allele frequencies of all variants between MFM WB and non-MFM WB and coding variants with moderate or severe predicted effects in MFM WB with publicly available data of other breeds. To compare differential gene expression and muscle fibre contractile force between MFM and non-MFM WB.

Study design

Case-control.

Animals

8 MFM WB, 8 non-MFM WB, 33 other WB, 32 Thoroughbreds, 80 Quarter Horses and 77 horses of other breeds in public databases.

Methods

Variants were called within transcripts of 16 candidate genes using gluteal muscle mRNA sequences aligned to EquCab3.0 and allele frequencies compared by Fisher's exact test among MFM WB, non-MFM WB and public sequences across breeds. Candidate gene differential expression was determined between MFM and non-MFM WB by fitting a negative binomial generalised log-linear model per gene (false discovery rate <0.05). The maximal isometric force/cross-sectional area generated by isolated membrane-permeabilised muscle fibres was determined.

Results

None of the 426 variants identified in 16 candidate genes were associated with MFM including 26 missense variants. Breed-specific differences existed in allele frequencies. Candidate gene differential expression and muscle fibre-specific force did not differ between MFM WB (143.1 ± 34.7 kPa) and non-MFM WB (140.2 ± 43.7 kPa) (P = .8).

Main limitations

RNA-seq-only assays transcripts expressed in skeletal muscle. Other possible candidate genes were not evaluated.

Conclusions

Evidence for association of variants with a disease is essential because coding sequence variants are common in the equine genome. Variants identified in MFM candidate genes, including two coding variants offered as commercial MFM equine genetic tests, did not associate with the WB MFM phenotype.

Background

Commercial genetic tests for type 2 polysaccharide storage myopathy (PSSM2) and myofibrillar myopathy (MFM) have not been validated by peer-review, and formal regulation of veterinary genetic testing is lacking.

Objectives

To compare genotype and allele frequencies of commercial test variants (P variants) in MYOT (P2; rs1138656462), FLNC (P3a; rs1139799323), FLNC (P3b; rs1142918816) and MYOZ3 (P4; rs1142544043) between Warmblood (WB) and Arabian (AR) horses diagnosed with PSSM2/MFM by muscle histopathology, and phenotyped breed-matched controls. To quantify variant frequency in public repositories of ancient and modern horse breeds.

Study design

Cross sectional using archived clinical material and publicly available data.

Methods

We studied 54 control-WB, 68 PSSM2/MFM-WB, 30 control-AR, 30 PSSM2/MFM-AR and 205 public genotypes. Variants were genotyped by pyrosequencing archived DNA. Genotype and allele frequency, and number of variant alleles or loci were compared within breed between controls, PSSM2/MFM combined and MFM or PSSM2 horses considered separately using additive/genotypic and dominant models (Fisher's exact tests). Variant frequencies in modern, early domestic and Przewalski horses were determined from a public data repository.

Results

There was no significant association between any P locus and a histopathological diagnosis of PSSM2/MFM, and no difference between control and myopathic horses in total loci with alternative alleles, or total alternate alleles when PSSM2/MFM was considered combined or separately as PSSM2 or MFM. For all tests, sensitivity was <0.33. Allele frequencies in WB (controls/cases) were: 8%/15% (P2), 5%/6% (P3a/b) and 9%/13% (P4); in AR, frequencies were: 12%/17% (P2), 2%/2% (P3a/b) and 7%/12% (P4). All P variants were present in early domestic (400- to 5500-year-old) horses and P2 present in the Przewalski.

Conclusions

Because of the lack of significant association between a histopathological diagnosis of PSSM2 or MFM and the commercial genetic test variants P2, P3 and P4 in WB and AR, we cannot recommend the use of these variant genotypes for selection and breeding, prepurchase examination or diagnosis of a myopathy.

Background

Sarcolipin (SLN), myoregulin (MRLN), and dwarf open reading frame (DWORF) are transmembrane regulators of the sarcoplasmic reticulum calcium transporting ATPase (SERCA) that we hypothesized played a role in recurrent exertional rhabdomyolysis (RER).

Objectives

Compare coding sequences of SLN, MRLN, DWORF across species and between RER and control horses. Compare expression of muscle Ca²⁺ regulatory genes between RER and control horses.

Animals

Twenty Thoroughbreds (TB), 5 Standardbreds (STD), 6 Quarter Horses (QH) with RER and 39 breed-matched controls.

Methods

Sanger sequencing of SERCA regulatory genes with comparison of amino acid (AA) sequences among control, RER horses, human, mouse, and rabbit reference genomes. In RER and control gluteal muscle, quantitative real-time polymerase chain reaction of SERCA regulatory peptides, the calcium release channel (RYR1), and its accessory proteins calsequestrin (CASQ1), and calstabin (FKBP1A).

Results

The SLN gene was the highest expressed horse SERCA regulatory gene with a uniquely truncated AA sequence (29 versus 31) versus other species. Coding sequences of SLN, MRLN, and DWORF were identical in RER and control horses. A sex-by-phenotype effect occurred with lower CASQ1 expression in RER males versus control males (P < .001) and RER females (P = .05) and higher FKBP1A (P = .01) expression in RER males versus control males.

Conclusions and clinical importance

The SLN gene encodes a uniquely truncated peptide in the horse versus other species. Variants in the coding sequence of SLN, MLRN, or DWORF were not associated with RER. Males with RER have differential gene expression that could reflect adaptations to stabilize RYR1.

The horse reference genome assemblies, EquCab2.0 and EquCab3.0, have enabled great advancements in the equine genomics field, from tools to novel discoveries. However, significant gaps of knowledge regarding genome function remain, hindering the study of complex traits in horses. In an effort to address these gaps and with inspiration from the Encyclopedia of DNA Elements (ENCODE) project, the equine Functional Annotation of Animal Genome (FAANG) initiative was proposed to bridge the gap between genome and gene expression, providing further insights into functional regulation within the horse genome. Three years after launching the initiative, the equine FAANG group has generated data from more than 400 experiments using over 50 tissues, targeting a variety of regulatory features of the equine genome. In this review, we examine how valuable lessons learned from the ENCODE project informed our decisions in the equine FAANG project. We report the current state of the equine FAANG project and discuss how FAANG can serve as a template for future expansion of functional annotation in the equine genome and be used as a reference for studies of complex traits in horse. A well-annotated reference functional atlas will also help advance equine genetics in the pan-genome and precision medicine era.

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.