Genome‐wide association study of osteochondrosis in the tarsocrural joint of Dutch Warmblood horses identifies susceptibility loci on chromosomes 3 and 10

Equine osteochondrosis is a developmental joint disease that is a significant source of morbidity affecting multiple breeds of horse. The genetic variants underlying osteochondrosis susceptibility have not been established. Here, we describe the results of a genome‐wide association study of osteochondrosis using 90 cases and 111 controls from a population of Dutch Warmblood horses. We report putative associations between osteochondrosis and loci on chromosome 3 (BIEC2‐808543; P = 5.03 × 10^?7) and chromosome 10 (BIEC2‐121323; P = 2.62 × 10^?7).     

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.     

Novel variants in the KIT and PAX3 genes in horses with white‐spotted coat colour phenotypes

Variants in the EDNRB, KIT, MITF, PAX3 and TRPM1 genes are known to cause white spotting phenotypes in horses, which can range from the common white markings up to completely white horses. In this study, we investigated these candidate genes in 169 horses with white spotting phenotypes not explained by the previously described variants. We identified a novel missense variant, PAX3:p.Pro32Arg, in Appaloosa horses with a splashed white phenotype in addition to their leopard complex spotting patterns. We also found three novel variants in the KIT gene. The splice site variant c.1346+1G>A occurred in a Swiss Warmblood horse with a pronounced depigmentation phenotype. The missense variant p.Tyr441Cys was present in several part‐bred Arabians with sabino‐like depigmentation phenotypes. Finally, we provide evidence suggesting that the common and widely distributed KIT:p.Arg682His variant has a very subtle white‐increasing effect, which is much less pronounced than the effect of the other described KIT variants. We termed the new KIT variants W18–W20 to provide a simple and unambiguous nomenclature for future genetic testing applications.     

A candidate gene approach identifies variants in SLC45A2 that explain dilute phenotypes,pearl and sunshine,in compound heterozygote horses

Variations in the SLC45A2 gene are responsible for the dilution phenotypes cream and pearl in domestic horses. Cream dilution is inherited in an incomplete dominant manner, diluting only red in the heterozygous state but both red and black pigments when two alleles are present. The pearl dilution is recessive and dilutes only the red and black pigment in the homozygous state or when paired with a cream allele. Horses that inherit one copy of pearl (C^prl) and one copy of the dominant cream allele (C^Cr) display a dilution phenotype similar to that of homozygous cream, suggesting that pearl is the result of a different variation in the same gene responsible for cream. We sequenced SLC45A2 in two ‘false double dilute’ horses that appeared phenotypically homozygous cream but tested as possessing only a single C^Cr allele. We also sequenced one known pearl carrier to screen for putative causal variants. The missense variant ECA21:SLC45A2:c.985G>A; p.Ala329Thr (C^prl) was present in one false double dilute and the pearl carrier and was also genotyped in an additional 126 horses for statistical evaluation. The genotype matched the expected phenotype in all horses (P‐value = 6.5 × 10^?41) and is identical to a pearl variant found previously. The second false double dilute horse and one non‐dilute offspring genotyped as heterozygous for a novel missense variant ECA21:SLC45A2:c.568G>A (p.Gly190Arg), the proposed C^sun variant (for the name of the horse). This variant produces a recessive dilution similar to pearl and indicates that multiple alleles of SLC45A2 result in dilution phenotypes in the domestic horse.     

Novel insights into Sabino1 and splashed white coat color patterns in horses

Within the framework of genome‐wide analyses using the novel Axiom^® genotyping array, we investigated the distribution of two previously described coat color patterns, namely sabino1 (SBI), associated with the KIT gene (KI16+1037A), and splashed white, associated with the PAX3 gene (ECA6:g.11429753C>T; PAX3^C70Y), including a total of 899 horses originating from eight different breeds (Achal Theke, Purebred Arabian, Partbred Arabian, Anglo‐Arabian, Shagya Arabian, Haflinger, Lipizzan and Noriker). Based on the data we collected we were able to demonstrate that, besides Quarter horses, the PAX3^C70Y allele is also present in Noriker (seven out of 189) and Lipizzan (three out of 329) horses. The SB1 allele was present in three breeds (Haflinger, 14 out of 98; Noriker, four out of 189; Lipizzan one out of 329). Furthermore, we examined the phenotypes of SB1‐ and PAX3^C70Y‐carrier horses for their characteristic white spotting patterns. None of the SB1/sb1‐carrier horses met the criteria defining the Sabino1 pattern according to current applied protocols. From 10 heterozygous PAX3^C70Y‐carrier horses, two had nearly a splashed white phenotype. The results of this large‐scale experiment on the genetic association of white spotting patterns in horses underline the influence of gene interactions and population differences on complex traits such as Sabino1 and splashed white.     

Explicit evidence for a missense mutation in exon 4 of SLC45A2 gene causing the pearl coat dilution in horses

Four loci seem responsible for the dilution of the basic coat colours in horse: Dun (D), Silver Dapple (Z), Champagne (CH) and Cream (C). Apart from the current phenotypes ascribed to these loci, pearl has been described as yet another diluted coat colour in this species. To date, this coat colour seems to segregate only in the Iberian breeds Purebred Spanish horse and Lusitano and has also been described in breeds of Iberian origin, such as Quarter Horses and Paint Horse, where it is referred to as the ‘Barlink Factor’. This phenotype segregates in an autosomal recessive manner and resembles some of the coat colours produced by the champagne CH^CH and cream C^Cr alleles, sometimes being difficult to distinguish among them. The interaction between compound heterozygous for the pearl C^prl and cream C^Cr alleles makes SLC45A2 the most plausible candidate gene for the pearl phenotype in horses. Our results provide documented evidence for the missense variation in exon 4 [SLC45A2:c.985G>A; SLC45A2:p.(Ala329Thr)] as the causative mutation for the pearl coat colour. In addition, it is most likely involved as well in the cremello, perlino and smoky cream like phenotypes associated with the compound C^Cr and C^prl heterozygous genotypes (known as cream pearl in the Purebred Spanish horse breed). The characterization of the pearl mutation allows breeders to identify carriers of the C^prl allele and to select this specific coat colour according to personal preferences, market demands or studbook requirements as well as to verify segregation within particular pedigrees.     

A genome‐wide association study in Caucasian women suggests the involvement of HLA genes in the severity of facial solar lentigines

Solar lentigines are a common feature of sun‐induced skin ageing. Little is known, however, about the genetic factors contributing to their development. In this genome‐wide association study, we aimed to identify genetic loci associated with solar lentigines on the face in 502 middle‐aged French women. Nine SNPs, gathered in two independent blocks on chromosome 6, exhibited a false discovery rate below 25% when looking for associations with the facial lentigine score. The first block, in the 6p22 region, corresponded to intergenic SNPs and also exhibited a significant association with forehead lentigines (P = 1.37 × 10^?8). The second block, within the 6p21 HLA region, was associated with decreased HLA‐C expression according to several eQTL databases. Interestingly, these SNPs were also in high linkage disequilibrium with the HLA‐C*0701 allele (r² = 0.95). We replicated an association recently found by GWAS in the IRF4 gene. Finally, a complementary study on 44 selected candidate SNPs revealed novel associations in the MITF gene. Overall, our results point to several mechanisms involved in the severity of facial lentigines, including HLA/immunity and the melanogenesis pathway.     

Understanding the stoichiometric limitation of herbivore growth: the importance of feeding and assimilation flexibilities

Ecological stoichiometry suggests that herbivore growth is limited by phosphorus when this element in the diet is < 8.6 μg P mg C^?1 (C : P atomic ratio > 300). However, in nature, it is not necessarily related to the relative phosphorus content in diets. This may be the result of complex feeding and assimilation responses to diets. We examined these possibilities using herbivorous plankton fed mono‐specific and mixed algae varying in phosphorus content of 1.6 to 8.1 μg P mg C^?1. The herbivores showed a 10‐fold growth rate difference among the diets. Growth rates related poorly with phosphorus content in the diets (r² = 0.07), better with P ingestion rate (r² = 0.41) and best with phosphorus assimilation rate (r² = 0.69). Inclusion of assimilation rates for carbon and fatty acids increased 7% of the explained growth variance. These results indicate that the feeding and assimilation flexibilities play pivotal roles in acquiring a deficient element and in regulating growth rate.     

Carbon isotope discrimination in leaves of the broad‐leaved paperbark tree,Melaleuca quinquenervia,as a tool for quantifying past tropical and subtropical rainfall

Quantitative reconstructions of terrestrial climate are highly sought after but rare, particularly in Australia. Carbon isotope discrimination in plant leaves (Δ_leaf) is an established indicator of past hydroclimate because the fractionation of carbon isotopes during photosynthesis is strongly influenced by water stress. Leaves of the evergreen tree Melaleuca quinquenervia have been recovered from the sediments of some perched lakes on North Stradbroke and Fraser Islands, south‐east Queensland, eastern Australia. Here, we examine the potential for using M. quinquenervia ?_leaf as a tracer of past rainfall by analysing carbon isotope ratios (δ¹³C) of modern leaves. We firstly assess Δ_leaf variation at the leaf and stand scale and find no systematic pattern within leaves or between leaves due to their position on the tree. We then examine the relationships between climate and Δ_leaf for a 11‐year time series of leaves collected in a litter tray. M. quinquenervia retains its leaves for 1–4 years; thus, cumulative average climate data are used. There is a significant relationship between annual mean ?_leaf and mean annual rainfall of the hydrological year for 1–4 years (i.e. 365–1460 days) prior to leaf fall (r² = 0.64, P = 0.003, n = 11). This relationship is marginally improved by accounting for the effect of pCO₂ on discrimination (r² = 0.67, P = 0.002, n = 11). The correlation between rainfall and Δ_leaf, and the natural distribution of Melaleuca quinquenervia around wetlands of eastern Australia, Papua New Guinea and New Caledonia offers significant potential to infer past rainfall on a wide range of spatial and temporal scales.     

Genetic contributions to precocity traits in racing Thoroughbreds

Adaptation to early training and racing (i.e. precocity), which is highly variable in racing Thoroughbreds, has implications for the selection and training of horses. We hypothesised that precocity in Thoroughbred racehorses is heritable. Age at first sprint training session (work day), age at first race and age at best race were used as phenotypes to quantify precocity. Using high‐density SNP array data, additive SNP heritability () was estimated to be 0.17, 0.14 and 0.17 for the three traits respectively. In genome‐wide association studies (GWAS) for age at first race and age at best race, a 1.98‐Mb region on equine chromosome 18 (ECA18) was identified. The most significant association was with the myostatin (MSTN) g.66493737C>T SNP (P = 5.46 × 10^?12 and P = 1.89 × 10^?14 respectively). In addition, two SNPs on ECA1 (g.37770220G>A and g.37770305T>C) within the first intron of the serotonin receptor gene HTR7 were significantly associated with age at first race and age at best race. Although no significant associations were identified for age at first work day, the MSTN:g.66493737C>T SNP was among the top 20 SNPs in the GWAS (P = 3.98 × 10^?5). Here we have identified variants with potential roles in early adaptation to training. Although there was an overlap in genes associated with precocity and distance aptitude (i.e. MSTN), the HTR7 variants were more strongly associated with precocity than with distance. Because HTR7 is closely related to the HTR1A gene, previously implicated in tractability in young Thoroughbreds, this suggests that behavioural traits may influence precocity.     

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.     

Whole genome sequencing reveals a novel deletion variant in the KIT gene in horses with white spotted coat colour phenotypes

White spotting phenotypes in horses can range in severity from the common white markings up to completely white horses. EDNRB, KIT, MITF, PAX3 and TRPM1 represent known candidate genes for such phenotypes in horses. For the present study, we re‐investigated a large horse family segregating a variable white spotting phenotype, for which conventional Sanger sequencing of the candidate genes’ individual exons had failed to reveal the causative variant. We obtained whole genome sequence data from an affected horse and specifically searched for structural variants in the known candidate genes. This analysis revealed a heterozygous ~1.9‐kb deletion spanning exons 10–13 of the KIT gene (chr3:77,740,239_77,742,136del1898insTATAT). In continuity with previously named equine KIT variants we propose to designate the newly identified deletion variant W22. We had access to 21 horses carrying the W22 allele. Four of them were compound heterozygous W20/W22 and had a completely white phenotype. Our data suggest that W22 represents a true null allele of the KIT gene, whereas the previously identified W20 leads to a partial loss of function. These findings will enable more precise genetic testing for depigmentation phenotypes in horses.     

Replication and fine-mapping of a QTL for recurrent airway obstruction in European Warmblood horses

Recurrent airway obstruction (RAO), or ‘heaves’, is a common performance‐limiting allergic respiratory disease of mature horses. It is related to sensitization and exposure to mouldy hay and has a familial basis with a complex mode of inheritance. In a previous study, we detected a QTL for RAO on ECA 13 in a half‐sib family of European Warmblood horses. In this study, we genotyped additional markers in the family and narrowed the QTL down to about 1.5 Mb (23.7–25.2 Mb). We detected the strongest association with SNP BIEC2‐224511 (24 309 405 bp). We also obtained SNP genotypes in an independent cohort of 646 unrelated Warmblood horses. There was no genome‐wide significant association with RAO in these unrelated horses. However, we performed a genotypic association study of the SNPs on ECA 13 in these unrelated horses, and the SNP BIEC2‐224511 also showed the strongest association with RAO in the unrelated horses (p_raw = 0.00037). The T allele at this SNP was associated with RAO both in the family and the unrelated horses. Thus, the association study in the unrelated animals provides independent support for the previously detected QTL. The association study allows further narrowing of the QTL interval to about 0.5 Mb (24.0–24.5 Mb). We sequenced the coding regions of the genes in the critical region but did not find any associated coding variants. Therefore, the causative variant underlying this QTL is likely to be a regulatory mutation.     

Variation in Aggressiveness Components in the Hemileia vastatrix Population in Brazil

The Pucciniomycete fungus Hemileia vastatrix causes leaf rust on coffee trees. The pathogen is responsible for considerable yield losses in susceptible coffee cultivars if appropriate management strategies are not implemented. Rapid spread and epidemics of rust fungi are usually associated with the emergence of new races of the pathogen that overcome resistance or with the emergence of more aggressive populations of the pathogen. In Brazil, coffee production is dominated by susceptible cultivars of Coffea arabica and Coffea canephora. We assessed aggressiveness in 46 populations of H. vastatrix from Minas Gerais and Espírito Santo, two of the most important coffee‐producing states in Brazil. We observed a significant difference in the incubation period between the populations from Minas Gerais and Espírito Santo when 183 single‐pustule isolates were inoculated onto Catuaí Vermelho IAC 44, a susceptible C. arabica cultivar. Variation in aggressiveness components was observed between and within localities. Isolates with longer incubation periods also tended to have longer latent periods, although there was only a low correlation between these two aggressiveness components (r² = 0.34, P = 2.2 × 10^?16). Low‐sporulating isolates also had significantly longer incubation and latent periods. The H. vastatrix population from Minas Gerais and Espírito Santo is structured by the formation of groups of individuals with differential level of aggressiveness. Our results indicate that the variation in aggressiveness of the Brazilian H. vastatrix population may be associated with the geographic coffee‐producing areas.     

Age‐associated microRNA expression in human peripheral blood is associated with all‐cause mortality and age‐related traits

Recent studies provide evidence of correlations of DNA methylation and expression of protein‐coding genes with human aging. The relations of microRNA expression with age and age‐related clinical outcomes have not been characterized thoroughly. We explored associations of age with whole‐blood microRNA expression in 5221 adults and identified 127 microRNAs that were differentially expressed by age at P < 3.3 × 10^?4 (Bonferroni‐corrected). Most microRNAs were underexpressed in older individuals. Integrative analysis of microRNA and mRNA expression revealed changes in age‐associated mRNA expression possibly driven by age‐associated microRNAs in pathways that involve RNA processing, translation, and immune function. We fitted a linear model to predict ‘microRNA age’ that incorporated expression levels of 80 microRNAs. MicroRNA age correlated modestly with predicted age from DNA methylation (r = 0.3) and mRNA expression (r = 0.2), suggesting that microRNA age may complement mRNA and epigenetic age prediction models. We used the difference between microRNA age and chronological age as a biomarker of accelerated aging (Δage) and found that Δage was associated with all‐cause mortality (hazards ratio 1.1 per year difference, P = 4.2 × 10^?5 adjusted for sex and chronological age). Additionally, Δage was associated with coronary heart disease, hypertension, blood pressure, and glucose levels. In conclusion, we constructed a microRNA age prediction model based on whole‐blood microRNA expression profiling. Age‐associated microRNAs and their targets have potential utility to detect accelerated aging and to predict risks for age‐related diseases.     

Late‐onset progressive retinal atrophy in the Gordon and Irish Setter breeds is associated with a frameshift mutation in C2orf71

Progressive retinal atrophy (PRA) in dogs is characterised by the degeneration of the photoreceptor cells of the retina, resulting in vision loss and eventually complete blindness. The condition affects more than 100 dog breeds and is known to be genetically heterogeneous between breeds. Around 14 mutations have now been identified that are associated with PRA in around 49 breeds, but for the majority of breeds the mutation(s) responsible have yet to be identified. Using genome‐wide association with 16 Gordon Setter PRA cases and 22 controls, we identified a novel PRA locus, termed rod–cone degeneration 4 (rcd4), on CFA17 (P_raw = 2.22 × 10^?8, P_genome = 2.00 × 10^?5), where a 3.2‐Mb region was homozygous within cases. A frameshift mutation was identified in C2orf71, a gene located within this region. This variant was homozygous in 19 of 21 PRA cases and was at a frequency of approximately 0.37 in the Gordon Setter population. Approximately 10% of cases in our study (2 of 21) are not associated with this C2orf71 mutation, indicating that PRA in this breed is genetically heterogeneous and caused by at least two mutations. This variant is also present in a number of Irish Setter dogs with PRA and has an estimated allele frequency of 0.26 in the breed. The function of C2orf71 remains unknown, but it is important for retinal development and function and has previously been associated with autosomal recessive retinitis pigmentosa in humans.     

Immunochip analysis identifies association of the RAD50/IL13 region with human longevity

Human longevity is characterized by a remarkable lack of confirmed genetic associations. Here, we report on the identification of a novel locus for longevity in the RAD50/IL13 region on chromosome 5q31.1 using a combined European sample of 3208 long‐lived individuals (LLI) and 8919 younger controls. First, we performed a large‐scale association study on 1458 German LLI (mean age 99.0 years) and 6368 controls (mean age 57.2 years) by targeting known immune‐associated loci covered by the Immunochip. The analysis of 142 136 autosomal single nucleotide polymorphisms (SNPs) revealed an Immunochip‐wide significant signal (P_Immunochip = 7.01 × 10^–9) for the SNP rs2075650 in the TOMM40/APOE region, which has been previously described in the context of human longevity. To identify novel susceptibility loci, we selected 15 markers with P_Immunochip < 5 × 10^–4 for replication in two samples from France (1257 LLI, mean age 102.4 years; 1811 controls, mean age 49.1 years) and Denmark (493 LLI, mean age 96.2 years; 740 controls, mean age 63.1 years). The association at SNP rs2706372 replicated in the French study collection and showed a similar trend in the Danish participants and was also significant in a meta‐analysis of the combined French and Danish data after adjusting for multiple testing. In a meta‐analysis of all three samples, rs2706372 reached a P‐value of P_{Immunochip+Repl} = 5.42 × 10^?7 (OR = 1.20; 95% CI = 1.12–1.28). SNP rs2706372 is located in the extended RAD50/IL13 region. RAD50 seems a plausible longevity candidate due to its involvement in DNA repair and inflammation. Further studies are needed to identify the functional variant(s) that predispose(s) to a long and healthy life.     

Spatial and Temporal Patterns of Frugivorous Hornbill Movements in Central Africa and their Implications for Rain Forest Conservation

Tropical forest conservation and restoration require an understanding of the movements and habitat preferences of important seed dispersers. With forests now being altered at an unprecedented rate, avian frugivores are becoming increasingly vital for forest regeneration. Seed movement, however, is highly dependent on the behavioral characteristics of their dispersers. Here, we examined the movements, habitat preferences, and range sizes of two African frugivores: the Black‐casqued (Ceratogymna atrata) and the White‐thighed (Bycanistes albotibialis) Hornbill, in the lowland rain forests of southern Cameroon. Using satellite transmitters, we tracked eight hornbills for 3 yr to characterize their movements and relate them to environmental landscape features. Hornbill movements differed significantly, with B. albotibialis ranging over larger areas (mean = 20,274 ha) than C. atrata (mean = 5604 ha), and females of both species covering over 15 times the area of males. Evidence suggests that movements are irruptive during particular periods, perhaps driven by low resource availability. In addition, hornbills often returned to the same localities within a year, although movements were not characterized as migratory. Both species displayed significant differences in habitat preference, with B. albotibialis utilizing disturbed habitat more frequently than C. atrata (t = ?22.04, P = 2.2 × 10^?16). Major roads were found to act as barriers for C. atrata, but not for B. albotibialis. The ability of both hornbill species to move large distances suggests hornbills will play a vital role in the maintenance and regeneration of rain forests in Central Africa as forest fragmentation increases and terrestrial vertebrates decline in numbers.