Sequence analysis of a 212 kb defensin gene cluster on ECA 27q17

Defensins are a family of evolutionary ancient antimicrobial peptides consisting of three sub-families: alpha-, beta- and theta-defensins. This investigation was focused on the genomic characterization of equine beta-defensins and the investigation of the potential clustering of beta-defensin genes in the equine genome. Six genomic BAC clones were isolated from the CHORI-241 library and one of these was mapped by FISH to ECA 27q17. This location was confirmed by RH-mapping. The contiguous 212 kb sequence of this clone was determined. Sequence analysis revealed the identification of ten pseudogenes and nine genes, six of which were highly homologous to human beta-defensin DEFB4. Clustering of the beta-defensin genes was confirmed and the order of the genes on the analyzed BAC was related to the corresponding defensin cluster on HSA 8. The knowledge about the sequence and the genomic structure of the equine beta-defensin genes will improve the classification of different paralogous defensin genes and is a prerequisite for subsequent functional studies. Additionally, the first alpha-defensin-like sequence outside the groups of primates, lagomorphs and rodents (glires) was identified.

Vulnerability to risk among small farmers in Tajikistan: results of a 2011 survey

Tajikistan is judged to be highly vulnerable to risk, including food insecurity risks and climate change risks. By some vulnerability measures it is the most vulnerable among all 28 countries in the World Bank’s Europe and Central Asia Region – ECA (World Bank 2009). The rural population, with its relatively high incidence of poverty, is particularly vulnerable. The Pilot Program for Climate Resilience (PPCR) in Tajikistan (2011) provided an opportunity to conduct a farm-level survey with the objective of assessing various dimensions of rural population’s vulnerability to risk and their perception of constraints to farming operations and livelihoods. The survey should be accordingly referred to as the 2011 PPCR survey. The rural population in Tajikistan is highly agrarian, with about 50% of family income deriving from agriculture (see Figure 4.1; also LSMS 2007 – own calculations). Tajikistan’s agriculture basically consists of two groups of producers: small household plots – the successors of Soviet “private agriculture” – and dehkan (or “peasant”) farms – new family farming structures that began to be created under relevant legislation passed after 1992 (Lerman and Sedik, 2008). The household plots manage 20% of arable land and produce 65% of gross agricultural output (GAO). Dehkan farms manage 65% of arable land and produce close to 30% of GAO. The remaining 15% of arable land is held in agricultural enterprises – the rapidly shrinking sector of corporate farms that succeeded the Soviet kolkhozes and sovkhozes and today produces less than 10% of GAO (TajStat 2011) The survey conducted in May 2011 focused on dehkan farms, as budgetary constraints precluded the inclusion of household plots. A total of 142 dehkan farms were surveyed in face-to-face interviews. They were sampled from 17 districts across all four regions – Sughd, Khatlon, RRP, and GBAO. The districts were selected so as to represent different agro-climatic zones, different vulnerability zones (based on the World Bank (2011) vulnerability assessment), and different food-insecurity zones (based on WFP/IPC assessments). Within each district, 3-4 jamoats were chosen at random and 2-3 farms were selected in each jamoat from lists provided by jamoat administration so as to maximize the variability by farm characteristics. The sample design by region/district is presented in Table A, which also shows the agro-climatic zone and the food security phase for each district. The sample districts are superimposed on a map of food security phases based on IPC April 2011.

IgE, IgGa, IgGb and IgG(T) serum antibody levels in offspring of two sires affected with equine recurrent airway obstruction

Equine recurrent airway obstruction (RAO) is a chronic lower airway disease of the horse caused by hypersensitivity reactions to inhaled stable dust, including mould spores such as Aspergillus fumigatus. The goals of this study were to investigate whether total serum IgE levels and allergen-specific IgE and IgG subclasses are influenced by genetic factors and/or RAO and whether quantitative trait loci (QTL) could be identified for these parameters. The offspring of two RAO-affected sires (S1: n=56 and S2: n=65) were grouped by stallion and disease status, and total serum IgE levels and specific IgE, IgGa, IgGb and IgG(T) levels against recombinant Aspergillus fumigatus 7 (rAspf7) were measured by ELISA. A panel of 315 microsatellite markers covering the 31 equine autosomes were used to genotype the stallions and their offspring. A whole-genome scan using half-sib regression interval mapping was performed for each of the IgG and IgE subclasses. There was no significant effect of disease status or sire on total IgE levels, but there was a significant effect of gender and age. rAspf7-specific IgGa levels were significantly higher in RAO-affected than in healthy horses. The offspring of S1 had significantly higher rAspf7-specific IgGa and IgE levels than those of S2. Five QTLs were significant chromosome-wide (P<0.01). QTLs for rAspf7-specific IgGa and IgE were identified on ECA 1, for rAspf7-specific IgGa and IgGb on ECA 24 and for rAspf7 IgGa on ECA 26. These results provide evidence for effects of disease status and genetics on allergen-specific IgGa and IgE.

Comparative human-horse sequence analysis of the CYP3A subfamily gene cluster

Cytochrome P450 enzymes (CYP450s) represent a superfamily of haem-thiolate proteins. CYP450s are most abundant in the liver, a major site of drug metabolism, and play key roles in the metabolism of a variety of substrates, including drugs and environmental contaminants. Interaction of two or more different drugs with the same enzyme can account for adverse effects and failure of therapy. Human CYP3A4 metabolizes about 50% of all known drugs, but little is known about the orthologous CYP450s in horses. We report here the genomic organization of the equine CYP3A gene cluster as well as a comparative analysis with the human CYP3A gene cluster. The equine CYP450 genes of the 3A family are located on ECA 13 between 6.97-7.53 Mb, in a region syntenic to HSA 7 99.05-99.35 Mb. Seven potential, closely linked equine CYP3A genes were found, in contrast to only four genes in the human genome. RNA was isolated from an equine liver sample, and the approximately 1.5-kb coding sequence of six CYP3A genes could be amplified by RT-PCR. Sequencing of the RT-PCR products revealed numerous hitherto unknown single nucleotide polymorphisms (SNPs) in these six CYP3A genes, and one 6-bp deletion compared to the reference sequence (EquCab2.0). The presence of the variants was confirmed in a sample of genomic DNA from the same horse. In conclusion, orthologous genes for the CYP3A family exist in horses, but their number differs from those of the human CYP3A gene family. CYP450 genes of the same family show high homology within and between mammalian species, but can be highly polymorphic.

Allelic Heterogeneity at the Equine KIT Locus in Dominant White (W) Horses

White coat color has been a highly valued trait in horses for at least 2,000 years. Dominant white (W) is one of several known depigmentation phenotypes in horses. It shows considerable phenotypic variation, ranging from approximately 50% depigmented areas up to a completely white coat. In the horse, the four depigmentation phenotypes roan, sabino, tobiano, and dominant white were independently mapped to a chromosomal region on ECA 3 harboring the KIT gene. KIT plays an important role in melanoblast survival during embryonic development. We determined the sequence and genomic organization of the approximately 82 kb equine KIT gene. A mutation analysis of all 21 KIT exons in white Franches-Montagnes Horses revealed a nonsense mutation in exon 15 (c.2151C>G, p.Y717X). We analyzed the KIT exons in horses characterized as dominant white from other populations and found three additional candidate causative mutations. Three almost completely white Arabians carried a different nonsense mutation in exon 4 (c.706A>T, p.K236X). Six Camarillo White Horses had a missense mutation in exon 12 (c.1805C>T, p.A602V), and five white Thoroughbreds had yet another missense mutation in exon 13 (c.1960G>A, p.G654R). Our results indicate that the dominant white color in Franches-Montagnes Horses is caused by a nonsense mutation in the KIT gene and that multiple independent mutations within this gene appear to be responsible for dominant white in several other modern horse populations.

A chromosomal region on ECA13 is associated with maxillary prognathism in horses.

Hereditary variations in head morphology and head malformations are known in many species. The most common variation encountered in horses is maxillary prognathism. Prognathism and brachygnathism are syndromes of the upper and lower jaw, respectively. The resulting malocclusion can negatively affect teeth wear, and is considered a non-desirable trait in breeding programs. We performed a case-control analysis for maxillary prognathism in horses using 96 cases and 763 controls. All horses had been previously genotyped with a commercially available 50 k SNP array. We analyzed the data with a mixed-model considering the genomic relationships in order to account for population stratification. Two SNPs within a region on the distal end of chromosome ECA 13 reached the Bonferroni corrected genome-wide significance level. There is no known prognathism candidate gene located within this region. Therefore, our findings in the horse offer the possibility of identifying a novel gene involved in the complex genetics of prognathism that might also be relevant for humans and other livestock species.

The equine DNAH3 gene: SNP discovery and exclusion of an involvement in recurrent airway obstruction (RAO) in European Warmblood horses

Recurrent airway obstruction is one of the most common airway diseases affecting mature horses. Increased bronchoalveolar mucus, neutrophil accumulation in airways, and airway obstruction are the main features of this disease. Mucociliary clearance is a key component of pulmonary defense mechanisms. Cilia are the motile part of this system and a complex linear array of dynein motors is responsible for their motility by moving along the microtubules in the axonemes of cilia and flagella. We previously detected a QTL for RAO on ECA 13 in a half-sib family of European Warmblood horses. The gene encoding DNAH3 is located in the peak of the detected QTL and encodes a dynein subunit. Therefore, we analysed this gene as a positional and functional candidate gene for RAO. In a mutation analysis of all 62 exons we detected 53 new polymorphisms including 7 non-synonymous variants. We performed an association study using 38 polymorphisms in a cohort of 422 animals. However, after correction for multiple testing we did not detect a significant association of any of these polymorphisms with RAO (P>0.05). Therefore, it seems unlikely that variants at the DNAH3 gene are responsible for the RAO QTL in European Warmblood horses.

Whole genome scan identifies several chromosomal regions linked to equine sarcoids.

Despite the evidence for a genetic predisposition to develop equine sarcoids (ES), no whole genome scan for ES has been performed to date. The objective of this explorative study was to identify chromosome regions associated with ES. The studied population was comprised of two half-sibling sire families, involving a total of 222 horses. Twenty-six of these horses were affected with ES. All horses had been previously genotyped with 315 microsatellite markers. Quantitative trait locus (QTL) signals were suggested where the F statistic exceeded chromosome-wide significance at P < 0.05. The QTL analyses revealed significant signals reaching P < 0.05 on equine chromosome (ECA) 20, 23 and 25, suggesting a polygenic character for this trait. The candidate regions identified on ECA 20, 23 and 25 include genes regulating virus replication and host immune response. Further investigation of the chromosome regions associated with ES and of genes potentially responsible for the development of ES could form the basis for early identification of susceptible animals, breeding selection or the development of new therapeutic targets.

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.

A Genome-Wide Association Study Reveals Loci Influencing Height and Other Conformation Traits in Horses

The molecular analysis of genes influencing human height has been notoriously difficult. Genome-wide association studies (GWAS) for height in humans based on tens of thousands to hundreds of thousands of samples so far revealed ∼200 loci for human height explaining only 20% of the heritability. In domestic animals isolated populations with a greatly reduced genetic heterogeneity facilitate a more efficient analysis of complex traits. We performed a genome-wide association study on 1,077 Franches-Montagnes (FM) horses using ∼40,000 SNPs. Our study revealed two QTL for height at withers on chromosomes 3 and 9. The association signal on chromosome 3 is close to the LCORL/NCAPG genes. The association signal on chromosome 9 is close to the ZFAT gene. Both loci have already been shown to influence height in humans. Interestingly, there are very large intergenic regions at the association signals. The two detected QTL together explain ∼18.2% of the heritable variation of height in horses. However, another large fraction of the variance for height in horses results from ECA 1 (11.0%), although the association analysis did not reveal significantly associated SNPs on this chromosome. The QTL region on ECA 3 associated with height at withers was also significantly associated with wither height, conformation of legs, ventral border of mandible, correctness of gaits, and expression of the head. The region on ECA 9 associated with height at withers was also associated with wither height, length of croup and length of back. In addition to these two QTL regions on ECA 3 and ECA 9 we detected another QTL on ECA 6 for correctness of gaits. Our study highlights the value of domestic animal populations for the genetic analysis of complex traits.

A Non-Synonymous HMGA2 Variant Decreases Height in Shetland Ponies and Other Small Horses

The identification of quantitative trait loci (QTL) such as height and their underlying causative variants is still challenging and often requires large sample sizes. In humans hundreds of loci with small effects control the heritable portion of height variability. In domestic animals, typically only a few loci with comparatively large effects explain a major fraction of the heritability. We investigated height at withers in Shetland ponies and mapped a QTL to ECA 6 by genome-wide association (GWAS) using a small cohort of only 48 animals and the Illumina equine SNP70 BeadChip. Fine-mapping revealed a shared haplotype block of 793 kb in small Shetland ponies. The HMGA2 gene, known to be associated with height in horses and many other species, was located in the associated haplotype. After closing a gap in the equine reference genome we identified a non-synonymous variant in the first exon of HMGA2 in small Shetland ponies. The variant was predicted to affect the functionally important first AT-hook DNA binding domain of the HMGA2 protein (c.83G>A; p.G28E). We assessed the functional impact and found impaired DNA binding of a peptide with the mutant sequence in an electrophoretic mobility shift assay. This suggests that the HMGA2 variant also affects DNA binding in vivo and thus leads to reduced growth and a smaller stature in Shetland ponies. The identified HMGA2 variant also segregates in several other pony breeds but was not found in regular-sized horse breeds. We therefore conclude that we identified a quantitative trait nucleotide for height in horses.

Selection signatures in Shetland ponies.

Shetland ponies were selected for numerous traits including small stature, strength, hardiness and longevity. Despite the different selection criteria, Shetland ponies are well known for their small stature. We performed a selection signature analysis including genome-wide SNPs of 75 Shetland ponies and 76 large-sized horses. Based upon this dataset, we identified a selection signature on equine chromosome (ECA) 1 between 103.8 Mb and 108.5 Mb. A total of 33 annotated genes are located within this interval including the IGF1R gene at 104.2 Mb and the ADAMTS17 gene at 105.4 Mb. These two genes are well known to have a major impact on body height in numerous species including humans. Homozygosity mapping in the Shetland ponies identified a region with increased homozygosity between 107.4 Mb and 108.5 Mb. None of the annotated genes in this region have so far been associated with height. Thus, we cannot exclude the possibility that the identified selection signature on ECA1 is associated with some trait other than height, for which Shetland ponies were selected.