A cohort study of racing performance in Japanese Thoroughbred racehorses using genome information on ECA18

Using 1710 Thoroughbred racehorses in Japan, a cohort study was performed to evaluate the influence of genotypes at four single nucleotide polymorphisms (SNPs) on equine chromosome 18 (ECA18), which were associated in a previous genome‐wide association study for racing performance with lifetime earnings and performance rank. In males, both g.65809482T>C and g.65868604G>T were related to performance rank (P = 0.005). In females, g.65809482T>C (P = 1.76E‐6), g.65868604G>T (P = 6.81E‐6) and g.66493737C>T (P = 4.42E‐5) were strongly related to performance rank and also to lifetime earnings (P < 0.05). When win‐race distance (WRD) among all winning racehorses and best race distance (BRD) among elite racehorses were considered as the phenotypes, significant associations (P < 0.001) were observed for all four SNPs. The favourable race distance of both elite (BRD) and novice racehorses (WRD) was also associated with genotypes in the ECA18 region, indicating the presence of a gene in this region influencing optimum race distance in Thoroughbred racehorses. Therefore, the association with performance rank is likely due to the bias in the race distances. The location of the SNPs within and proximal to the gene encoding myostatin (MSTN) strongly suggests that regulation of the MSTN gene affects racing performance. In particular, the g.65809482T>C, g.65868604G>T and g.66493737C>T SNPs, or their combinations, may be genetic diagnostic markers for racing performance indicators such as WRD and BRD.     

Analysis of genetic variation contributing to measured speed in Thoroughbreds identifies genomic regions involved in the transcriptional response to exercise

Despite strong selection for athletic traits in Thoroughbred horses, there is marked variation in speed and aptitude for racing performance within the breed. Using global positioning system monitoring during exercise training, we measured speed variables and temporal changes in speed with age to derive phenotypes for GWAS. The aim of the study was to test the hypothesis that genetic variation contributes to variation in end‐point physiological traits, in this case galloping speed measured during field exercise tests. Standardisation of field‐measured phenotypes was attempted by assessing horses exercised on the same gallop track and managed under similar conditions by a single trainer. PCA of six key speed indices captured 73.9% of the variation with principal component 1 (PC1). Verifying the utility of the phenotype, we observed that PC1 (median) in 2‐year‐old horses was significantly different among elite, non‐elite and unraced horses (P < 0.001) and the temporal change with age in PC1 varied among horses with different myostatin (MSTN) g.66493737C>T SNP genotypes. A GWAS for PC1 in 2‐year‐old horses (n = 122) identified four SNPs reaching the suggestive threshold for association (P < 4.80 × 10^?5), defining a 1.09 Mb candidate region on ECA8 containing the myosin XVIIIB (MYO18B) gene. In a GWAS for temporal change in PC1 with age (n = 168), five SNPs reached the suggestive threshold for association and defined candidate regions on ECA2 and ECA11. Both regions contained genes that are significantly differentially expressed in equine skeletal muscle in response to acute exercise and training stimuli, including MYO18A. As MYO18A plays a regulatory role in the skeletal muscle response to exercise, the identified genomic variation proximal to the myosin family genes may be important for the regulation of the response to exercise and training.     

A new single nucleotide polymorphism in the rabbit (Oryctolagus cuniculus) myostatin (MSTN) gene is associated with carcass composition traits

This study aimed at the identification of genetic variations in the myostatin (MSTN) gene and testing their effects on carcass quality traits. We comparatively sequenced Giant Grey (GG) and New Zealand White (NZW) rabbits that were founders of a cross‐bred population. Alignment of our sequence data with the GenBank sequence of the rabbit MSTN gene (Ensembl Gene ID ENSOCUG00000012663) identified three single nucleotide polymorphisms (SNPs). The two novel SNPs (c.?125T>C, c.373+234G>A) and one known SNP (c.747+34C>T) were subsequently analysed for linkage with carcass composition traits in 363 F₂ animals of the cross GG × NZW. Significant linkage was found between c.373+234G>A and nine carcass composition traits (P < 0.05). No significant effects were found for c.?125T>C and c.747+34C>T. Because the linked SNP is located in intron 1 and no genetic variation was found in the coding region, further investigations are necessary to understand the functional effect of the c.373+234G>A variant on the variability of the traits.     

Genome‐wide association analysis identifies quantitative trait loci for growth in a Landrace purebred population

Growth‐related traits are complex and economically important in the livestock industry. The aim of this study was to identify quantitative trait loci (QTL) and the associated positional candidate genes affecting growth in pigs. A genome‐wide association study (GWAS) was performed using the porcine single‐nucleotide polymorphism (SNP) 60K bead chip. A mixed‐effects model and linear regression approach were used for the GWAS. The data used in the study included 490 purebred Landrace pigs. All experimental animals were genotyped with 39 438 SNPs located throughout the pig autosomes. We identified a strong association between a SNP marker on chromosome 16 and body weight at 71 days of age (ALGA0092396, P = 5.35 × 10^?9, Bonferroni adjusted P < 0.05). The SNP marker was located near the genomic region containing IRX4, which encodes iroquois homeobox 4. This SNP marker could be useful in the selective breeding program after validating its effect on other populations.     

Genetic mapping of recurrent exertional rhabdomyolysis in a population of North American Thoroughbreds

Recurrent exertional rhabdomyolysis is a heritable disorder that results in painful skeletal muscle cramping with exercise in up to 10% of all Thoroughbred racehorses. Here, we report a genome‐wide association study with 48 282 SNPs analyzed among 48 case and 37 control Thoroughbreds. The most significant SNPs spanned approximately 13 Mb on ECA16, and the P‐value of the most significant SNP after correcting for population structure was 8.0 × 10^?6. This region on ECA16 was further evaluated by genotyping 247 SNPs in both the initial population and a second population of 34 case and 98 control Thoroughbreds. Several SNPs across the 13‐Mb region on ECA16 showed significance when each population was analyzed separately; however, the exact positions of the most significant SNPs within this region on ECA16 varied between populations. This variability in location may be attributed to lack of power owing to insufficient sample sizes within each population individually, or to the relative distribution of long, conserved haplotypes, characteristic of the Thoroughbred breed. Future genome‐wide association studies with additional horses would likely improve the power to resolve casual loci located on ECA16 and increase the likelihood of detecting any additional loci on other chromosomes contributing to disease susceptibility.     

Leukocyte telomere length in the Thoroughbred racehorse

Thoroughbred racehorses possess superior cardiorespiratory fitness levels and are at the pinnacle of athletic performance compared to other breeds of horses. Although equine athletes have undergone years of artificial selection for racing performance, musculoskeletal injuries and illnesses are common and concerns relating to animal welfare have been proposed. Leukocyte telomere length is indicative of biological age, and accelerated telomere shortening occurs with excess physical and psychological stress. This study was designed to explore the association between leukocyte telomere length, biological factors (age, sex and coat colour), training status, winnings and race history parameters. Blood was collected from 146 Thoroughbred racehorses from around Geelong, Victoria, Australia. DNA was extracted from leukocytes; telomere length was measured using qPCR and analysed in context with traits obtained from the Racing Australia website. Age was inversely correlated with telomere length (r = ?0.194, P = 0.019). The oldest horses (≥11 years) in the highest age quartile possessed shorter telomeres compared to younger horses in the first, second and third quartiles (≤2, 3–5 and 6–10 years respectively; P < 0.05). No statistically significant associations were observed between telomere length and biological factors, training status, winnings or race history parameters in age‐adjusted analyses. The study findings suggest that Thoroughbred horses may undergo age‐related telomere shortening similar to other mixed breeds and humans. Despite concerns from some quarters regarding the welfare of racehorses, there was a lack of accelerated biological ageing observed in the present study, as indicated by leukocyte telomere length.     

A candidate-SNP retrospective cohort study for fracture risk in Japanese Thoroughbred racehorses

Fractures are medical conditions that compromise the athletic potential of horses and/or the safety of jockeys. Therefore, the reduction of fracture risk is an important horse and human welfare issue. The present study used molecular genetic approaches to determine the effect of genetic risk for fracture at four candidate SNPs spanning the myostatin (MSTN) gene on horse chromosome 18. Among the 3706 Japanese Thoroughbred racehorses, 1089 (29.4%) had experienced fractures in their athletic life, indicating the common occurrence of this injury in Thoroughbreds. In the case/control association study, fractures of the carpus (carpal bones and distal radius) were statistically associated with g.65809482T/C (P = 1.17 x 10^-8), g.65868604G/T (P = 2.66 x 10^-9), and g.66493737C/T (P = 6.41 x 10^-8). In the retrospective cohort study using 1710 racehorses born in 2000, the relative risk (RR) was highest for male horses at g.65868604G/T, based on the dominant allele risk model (RR = 2.251, 95% confidence interval 1.407–3.604, P = 0.00041), and for female horses at g.65868604G/T, based on the recessive allele risk model (RR = 2.313, 95% confidence interval 1.380–3.877, P = 0.00163). Considering the association of these SNPs with racing performance traits such as speed, these genotypes may affect the occurrence of carpus fractures in Japanese Thoroughbred racehorses as a consequence of the non-genetic influence of the genotype on the distance and/or intensity of racing and training. The genetic information presented here may contribute to the development of strategic training programs and racing plans for racehorses that improve their health and welfare.     

Genetic risk for squamous cell carcinoma of the nictitating membrane parallels that of the limbus in Haflinger horses

Squamous cell carcinoma (SCC) is the most common cancer affecting the equine eye, with a higher incidence documented in Haflinger horses. Recently, a missense variant in the gene damage specific DNA binding protein 2 (DDB2, p.Thr338Met) on ECA12 was identified as a risk factor for the development of limbal SCC in Haflinger horses. SCC also occurs on the nictitating membrane; therefore, investigating the role of this missense variant in nictitating membrane SCC is warranted. In this study, a common ancestor was identified among Haflinger horses affected with limbal SCC or with nictitating membrane SCC, thus supporting a recessive risk factor for the development of cancer at both ocular locations. Analysis of genotype data from Haflinger horses with and without nictitating membrane SCC revealed that the same region on ECA12 associated with limbal SCC was also associated with nictitating membrane SCC (P < 2.04 × 10^?5). Fine mapping of this locus using 25 cases and 49 controls supported the hypothesis that DDB2:c.1013C>T, p.Thr338Met, is a risk factor for nictitating membrane SCC, as 88% of our cases were homozygous for this variant and no other polymorphism was more strongly associated (P = 4.13 × 10^?14). These data indicate that the genetic risk is the same for the development of both limbal and nictitating membrane SCC in Haflinger horses and validates utilization of genetic testing of the DDB2 variant for both clinical management and the guidance of mating decisions.     

A genome-wide association study for racing performances in Thoroughbreds clarifies a candidate region near the MSTN gene

Using 1400 microsatellites, a genome-wide association study (GWAS) was performed to identify genomic regions associated with lifetime earnings and performance ranks, as determined by the Japan Racing Association (JRA). The minimum heritability (h(2) ) was estimated at 7-8% based on the quantitative trait model, suggesting that the racing performance is heritable. Following GWAS with microsatellites, fine mapping led to identification of three SNPs on ECA18, namely, g.65809482T>C (P=1.05E-18), g.65868604G>T (P=6.47E-17), and g.66539967A>G (P=3.35E-14) associated with these performance measures. The haplotype of these SNPs, together with a recently published nearby SNP, g.66493737C>T (P=9.06E-16) in strong linkage disequilibrium, also showed a very clear association with the performance (P<1E-05). The candidate genomic region contained eight genes annotated by ENSEMBL, including the myostatin gene (MSTN). These findings suggest the presence of a gene affecting the racing performance in Thoroughbred racehorses in this region on ECA18.     

A Sequence Polymorphism in MSTN Predicts Sprinting Ability and Racing Stamina in Thoroughbred Horses

Variants of the MSTN gene encoding myostatin are associated with muscle hypertrophy phenotypes in a range of mammalian species, most notably cattle, dogs, mice, and humans. Using a sample of registered Thoroughbred horses (n = 148), we have identified a novel MSTN sequence polymorphism that is strongly associated (g.66493737C>T, P = 4.85×10⁻⁸) with best race distance among elite racehorses (n = 79). This observation was independently validated (P = 1.91×10⁻⁶) in a resampled group of Thoroughbreds (n = 62) and in a cohort of Thoroughbreds (n = 37, P = 0.0047) produced by the same trainer. We observed that C/C horses are suited to fast, short-distance races; C/T horses compete favorably in middle-distance races; and T/T horses have greater stamina. Evaluation of retrospective racecourse performance (n = 142) and stallion progeny performance predict that C/C and C/T horses are more likely to be successful two-year-old racehorses than T/T animals. Here we describe for the first time the identification of a gene variant in Thoroughbred racehorses that is predictive of genetic potential for an athletic phenotype.     

Heritability of racing performance in the Australian Thoroughbred racing population

Performance data for 164 046 Thoroughbreds entered in a race or official barrier trial in Australia were provided by Racing Information Services Australia. Analyses estimating the heritability for a range of racing performance traits using a single‐trait animal model were performed using asreml ‐r . Log of cumulative earnings (LCE; 0.19 ± 0.01), log of earnings per race start (0.23 ± 0.02) and best race distance (0.61 ± 0.03) were all significantly heritable. Fixed effects for sex were significant (P < 0.001) for all performance traits aside from LCE (P = 0.382). With the exception of annual earnings, trainer was also significant for all performance traits. As the application of modern genetic selection methodologies continues to gain popularity in the racing industry, contemporary heritability estimates from the current population of Thoroughbreds will play a vital role in identifying which traits are better suited to selection and in the development of more accurate genomic evaluations for racing performance.     

Evaluation of the causality of the low‐density lipoprotein receptor gene (LDLR) for serum lipids in pigs

A significant quantitative trait locus (QTL) for low‐density lipoprotein cholesterol (LDL‐C) and total cholesterol (TC) was identified around the LDLR gene on chromosome 2 (SSC2) in a White Duroc × Erhualian F₂ resource population and Sutai pigs in our previous study. However, in previous reports, the causality of LDLR with serum lipids is controversial in pigs. To systematically assess the causality of LDLR with serum lipids, association analyses were successively performed in three populations: Sutai pigs, a White Duroc × Erhualian F₂ resource population and a Duroc × (Landrace × Large White) population. We first performed a haplotype‐based association study with 60K SNP genotyping data and evidenced the significant association with LDL‐C and TC around the LDLR gene region. We also found that there is more than one QTL for LDL‐C and TC on SSC2. Then, we evaluated the causalities of two missense mutations, c.1812C>T and c.1520A>G, with LDL‐C and TC. We revealed that the c.1812C>T SNP showed the strongest association with LDL‐C (P = 5.40 × 10^?11) and TC (P = 3.64 × 10^?8) and explained all the QTL effect in Sutai pigs. Haplotype analysis found that two missense SNPs locate within a 1.93‐Mb haplotype block. One major haplotype showed the strongest significant association with LDL‐C (P = 4.62 × 10^?18) and TC (P = 1.06 × 10^?9). However, the c.1812C>T SNP was not identified in the White Duroc × Erhualian intercross, and the association of c.1520A>G with both LDL‐C and TC did not achieve significance in this F₂ population, suggesting population heterogeneity. Both missense mutations were identified in the Duroc × (Landrace × Large White) population and showed significant associations with LDL‐C and TC. Our data give evidence that the LDLR gene should be a candidate causative gene for LDL‐C and TC in pigs, but heterogeneity exists in different populations.     

A Genome-Wide Association Study Identifies Risk Loci to Equine Recurrent Uveitis in German Warmblood Horses

Equine recurrent uveitis (ERU) is a common eye disease affecting up to 3–15% of the horse population. A genome-wide association study (GWAS) using the Illumina equine SNP50 bead chip was performed to identify loci conferring risk to ERU. The sample included a total of 144 German warmblood horses. A GWAS showed a significant single nucleotide polymorphism (SNP) on horse chromosome (ECA) 20 at 49.3 Mb, with IL-17A and IL-17F being the closest genes. This locus explained a fraction of 23% of the phenotypic variance for ERU. A GWAS taking into account the severity of ERU, revealed a SNP on ECA18 nearby to the crystalline gene cluster CRYGA-CRYGF. For both genomic regions on ECA18 and 20, significantly associated haplotypes containing the genome-wide significant SNPs could be demonstrated. In conclusion, our results are indicative for a genetic component regulating the possible critical role of IL-17A and IL-17F in the pathogenesis of ERU. The associated SNP on ECA18 may be indicative for cataract formation in the course of ERU.     

Haplotype diversity in the equine myostatin gene with focus on variants associated with race distance propensity and muscle fiber type proportions

Two variants in the equine myostatin gene (MSTN), including a T/C SNP in the first intron and a 227‐bp SINE insertion in the promoter, are associated with muscle fiber type proportions in the Quarter Horse (QH) and with the prediction of race distance propensity in the Thoroughbred (TB). Genotypes from these loci, along with 18 additional variants surrounding MSTN, were examined in 301 horses of 14 breeds to evaluate haplotype relationships and diversity. The C allele of intron 1 was found in 12 of 14 breeds at a frequency of 0.27; the SINE was observed in five breeds, but common in only the TB and QH (0.73 and 0.48 respectively). Haplotype data suggest the SINE insertion is contemporary to and arose upon a haplotype containing the intron 1 C allele. Gluteal muscle biopsies of TBs showed a significant association of the intron 1 C allele and SINE with a higher proportion of Type 2B and lower proportion of Type 1 fibers. However, in the Belgian horse, in which the SINE is not present, the intron 1 SNP was not associated with fiber type proportions, and evaluation of fiber type proportions across the Belgian, TB and QH breeds shows the significant effect of breed on fiber type proportions is negated when evaluating horses without the SINE variant. These data suggest the SINE, rather than the intron 1 SNP, is driving the observed muscle fiber type characteristics and is the variant targeted by selection for short‐distance racing.     

A genome‐wide association study for the incidence of persistent bovine viral diarrhea virus infection in cattle

Bovine viral diarrhea viruses (BVDV) comprise a diverse group of viruses that cause disease in cattle. BVDV may establish both transient and persistent infections depending on the developmental stage of the animal at exposure. The objective was to determine whether genomic regions harboring single nucleotide polymorphisms (SNPs) could be associated with the presence or absence of persistent BVDV infection. A genome‐wide association approach based on 777 000 SNP markers was used. Samples of animals identified as positive (n = 1200) or negative (n = 1200) for the presence of BVDV in skin samples (n = 1200) were used. DNA samples were combined in 24 pools (100 animals per pool). One SNP, significant at the 5 percent genome‐wide level (P = 9.41 × 10^?8), was detected on chromosome 14, located at position 80 675 176 bp. Fifteen SNPs, residing on chromosomes 1, 2, 6, 8, 10, 15 and 18, were moderately associated (P < 1 × 10^?5) with persistent BVDV infection. Results show that genes harboring or neighboring significant SNPs are involved in leucopenia, signal transduction, RNA splicing and DNA methylation processes.     

Length‐weight relationship and biological data of a threatened fish,Ptychobarbus chungtienensis (Tsao, 1964) in Bita Lake of Shangri‐La,Yunnan, China

This paper provides the length‐weight relationship (LWR) and other biological information of a threatened plateau freshwater fish species Ptychobarbus chungtienensis in the subfamily Schizothoracinae of the family Cyprinidae, from the Bita Lake of Shangri‐la, in Yunnan, China. The LWR was BW = 0.00954 TL^2.95 for a total of 392 individuals, where the coefficient b was close to 3.0, suggesting a near‐isometric growth in this species. The total length: standard length relationship was TL = 0.399 + 1.151 SL. Fitting a von Bertalanffy growth function to the observed length‐at‐age data resulted in L_(t) = 53.1 (1‐e^?0.08(t?0.2)). The natural mortality rate (M) was 0.20 year^?1. Length at first maturity (Lm) was 20.57 cm, and 6.3 years the theoretical estimated age at first maturity (Am). Ptychobarbus chungtienensis mature gradually beginning in May, reaching a reproduction peak around August, with most individuals completing reproduction by October. The gonadosomatic index (GSI) was highest in August (GSI = 0.1351) and lowest in October (GSI = 0.0467). The male index (GSI = 0.0866) was lower than in females (GSI = 0.1466) (P < 0.05). Mean condition factor (Km) of P. chungtienensis was 0.8500, and not significantly different between sexes (P = 0.124), but differing significantly between seasons. In October the condition factor (Km = 0.9211) was significantly higher than in May (Km = 0.8379) or in August (Km = 0.8412) (P < 0.05).     

Single-nucleotide polymorphisms in the promoter of the growth hormone-releasing hormone receptor gene are associated with growth and reproduction traits in chickens

Growth hormone‐releasing hormone receptor (GHRHR) plays a critical role in growth hormone (GH) synthesis, release and regulation in animals. The objective of this study was to investigate variations of the chicken GHRHR gene and their associations with growth and reproduction traits in 768 Beijing You chickens. Results revealed three single nucleotide polymorphisms (SNPs) in the promoter region of the gene (g.‐1654A>G, g.‐1411A>G and g.‐142T>C). Association analysis revealed that the novel SNP g.‐1654A>G had significant effects on chicken body weight at 7, 9, 11, 13, 17 weeks of age and the age of first egg as well as egg number at 32, 36 and 40 weeks. Significant association was also observed between g.‐1411A>G and g.‐142T>C with EN24. Moreover, the age of first egg was distinctly related with g.‐142T>C (P < 0.05). Although significant statistical difference was not detected in GHRHR mRNA levels among genotypes of the SNPs (P > 0.05), strong expression variations of the gene were found between the ages 17 and 20 weeks in the population (P < 0.05). These results suggest that the three SNPs in the GHRHR promoter could be used as potential genetic markers to improve the growth and reproductive traits in chickens.     

Genomic analysis of the differential response to experimental infection with porcine circovirus 2b

Porcine circovirus type 2 (PCV2) is the etiological agent of a group of associated diseases (PCVAD) that affect production efficiency and can lead to mortality. Using different crossbred lines of pigs, we analyzed host genetic variation of viral load, immune response and weight change following experimental infection with a PCV2b strain (n = 386). Pigs expressed variation in the magnitude and initiation of viremia and immune response recorded weekly until 28 days post‐infection. A higher viral load was correlated with weight gain (r = ?0.26, P < 0.0001) and presence of PCV2‐specific antibodies (IgM, r = 0.26–0.34, P < 0.0001; IgG, r = 0.17–0.20, P < 0.01). In genome‐wide association analyses of the responses at different time points, the proportions of phenotypic variation explained by combined effects of 56 433 SNPs were 34.8–59.4% for viremia, 10.1–59.5% for antibody response and 5.6–14.9% for weight change. Relationships between genomic prediction of overall viral load and weight gain during the first weeks of challenge were negative (?0.21 and ?0.24 respectively, P < 0.0001). Individuals that carried more favorable alleles across three SNPs on SSC9 (0.60 Mb) and SSC12 (6.8 and 18.2 Mb) partially explained this relationship, having lower viral load (P < 0.0001); lower viremia at day 14 (P < 0.0001), day 21 (P < 0.01) and day 28 (P < 0.05) and greater overall average daily gain during infection (ADG_i; P < 0.01), ADG_i at week 3 (P < 0.001) and week 4 (P < 0.01). These additive genetic relationships could lead to molecular solutions to improve animal health and reduce production costs.