Detecting genes for variation in parasite burden and immunological traits in a wild population: testing the candidate gene approach.

Identifying the genes underlying phenotypic variation in natural populations can provide novel insight into the evolutionary process. The candidate gene approach has been applied to studies of a number of traits in various species, in an attempt to elucidate their genetic basis. Here, we test the application of the candidate gene approach to identify the loci involved in variation in gastrointestinal parasite burden, a complex trait likely to be controlled by many loci, in a wild population of Soay sheep. A comprehensive literature review, Gene Ontology databases, and comparative genomics resources between cattle and sheep were used to generate a list of candidate genes. In a pilot study, these candidates, along with 50 random genes, were then sequenced in two pools of Soay sheep; one with low gastrointestinal nematode burden and the other high, using a NimbleGen sequence capture experiment. Further candidates were identified from single nucleotide polymorphisms (SNPs) that were highly differentiated between high- and low-resistance sheep breeds. A panel of 192 candidate and control SNPs were then typed in 960 individual Soay sheep to examine whether they individually explained variation in parasite burden, as measured as faecal egg count, as well as two immune measures (Teladorsagia circumcincta-specific antibodies and antinuclear antibodies). The cumulative effect of the candidate and control SNPs were estimated by fitting genetic relationship matrices (GRMs) as random effects in animal models of the three traits. No more significant SNPs were identified in the pilot sequencing experiment and association study than expected by chance. Furthermore, no significant difference was found between the proportions of candidate or control SNPs that were found to be significantly associated with parasite burden/immune measures. No significant effect of the candidate or control gene GRMs was found. There is thus little support for the candidate gene approach to the identification of loci explaining variation in parasitological and immunological traits in this population. However, a number of SNPs explained significant variation in multiple traits and significant correlations were found between the proportions of variance explained by individual SNPs across multiple traits. The significant SNPs identified in this study may still, therefore, merit further investigation.

Digital gene expression analysis of gastrointestinal helminth resistance in Scottish blackface lambs.

Digital gene expression (DGE) analysis offers a route to gene discovery which by-passes the need to develop bespoke arrays for nonmodel species, and is therefore a potentially valuable tool for molecular ecologists. Scottish blackface sheep, which vary in resistance to the common abomasal parasitic nematode Teladorsagia circumcincta, were trickle-infected with L3 larvae over 3 months to mimic the natural progression of infection. DGE was performed on abomasal lymph node tissue after the resolution of infection in resistant animals. Susceptible (low resistance) animals showed a large number of differentially expressed genes associated with inflammation and cell activation, but generally few differentially regulated genes in either the susceptible or the resistant group were directly involved in the adaptive immune function. Our results are consistent with the hypothesis that both resistance and susceptibility are active responses to infection and that susceptibility is associated with dysfunction in T cell differentiation and regulation.

Mapping quantitative trait loci in a wild population using linkage and linkage disequilibrium analyses.

Historical information can be used, in addition to pedigree, traits and genotypes, to map quantitative trait locus (QTL) in general populations via maximum likelihood estimation of variance components. This analysis is known as linkage disequilibrium (LD) and linkage mapping, because it exploits both linkage in families and LD at the population level. The search for QTL in the wild population of Soay sheep on St. Kilda is a proof of principle. We analysed the data from a previous study and confirmed some of the QTLs reported. The most striking result was the confirmation of a QTL affecting birth weight that had been reported using association tests but not when using linkage-based analyses.

Horn type and horn length genes map to the same chromosomal region in Soay sheep.

The evolution of male weaponry in animals is driven by sexual selection, which is predicted to reduce the genetic variability underlying such traits. Soay sheep have an inherited polymorphism for horn type in both sexes, with males presenting with either large, normal horns or small, deformed horns (scurs). In addition, there is additive genetic variation in horn length among males with normal horns. Given that scurred males cannot win conflicts with normal-horned males, it is unusual that genes conferring scurs should persist in the population. Identifying the genetic basis of these traits should help us in understanding their evolution. We developed microsatellite markers in a targeted region of the Soay sheep genome and refined the location of the Horns locus (Ho) to a approximately 7.4 cM interval on chromosome 10 (LOD=8.78). We then located quantitative trait loci spanning a 34 cM interval with a peak centred close to Ho, which explained the majority of the genetic variation for horn length and base circumference in normal-horned males (LOD=2.51 and LOD=1.04, respectively). Therefore, the genetic variation in both horn type and horn length is attributable to the same chromosomal region. Understanding the maintenance of horn type and length variation will require an investigation of selection on genotypes that (co)determine both traits.

The genetic basis of recessive self-colour pattern in a wild sheep population.

Bridging the genotype-phenotype gap for traits of ecological and evolutionary importance in natural populations can provide a novel insight into the origin and maintenance of variation. Here, we identify the gene and putative causal mutations underlying a recessive colour pattern phenotype ('self' or uniform colour) in a wild population of primitive Soay sheep. We targeted the agouti signalling protein (ASIP) gene, a positional candidate based on previous study that mapped the Coat pattern locus to a presumptive region on chromosome 13. We found evidence for three recessive mutations, including two functional changes in the coding sequence and a putative third cis-regulatory mutation that inactivates the promoter. These mutations define up to five haplotypes in Soays, which collectively explained the coat pattern in all but one member of a complex multi-generational pedigree containing 621 genotyped individuals. The functional mutations are in strong linkage disequilibrium in the study population, and are identical to those known to underlie the self phenotype in domestic sheep. This is indicative of a recent (and simultaneous) origin in Soay sheep, possibly as a consequence of past interbreeding with modern domestic breeds. This is only the second study in which ASIP has been linked to variation in pigmentation in a natural population. Knowledge of the genetic basis of self-colour pattern in Soay sheep, and the recognition that several mutations are segregating in the population, will aid future studies investigating the role of selection in the maintenance of the polymorphism.

Compelling evidence that a single nucleotide substitution in TYRP1 is responsible for coat-colour polymorphism in a free-living population of Soay sheep.

Identifying the genes that underlie phenotypic variation in natural populations is a central objective of evolutionary genetics. Here, we report the identification of the gene and causal mutation underlying coat colour variation in a free-living population of Soay sheep (Ovis aries). We targeted tyrosinase-related protein 1 (TYRP1), a positional candidate gene based on previous work that mapped the Coat colour locus to an approximately 15cM window on chromosome 2. We identified a non-synonymous substitution in exon IV that was perfectly associated with coat colour. This polymorphism is predicted to cause the loss of a cysteine residue that is highly evolutionarily conserved and likely to be of functional significance. We eliminated the possibility that this association is due to the presence of strong linkage disequilibrium with an unknown regulatory mutation by demonstrating that there is no difference in relative TYRP1 expression between colour morphs. Analysis of this putative causal mutation in a complex pedigree of more than 500 sheep revealed almost perfect co-segregation with coat colour (chi2-test, p<0.0001, LOD=110.20), and very tight linkage between Coat colour and TYRP1 (LOD=29.50).

A localized negative genetic correlation constrains microevolution of coat color in wild sheep.

The evolutionary changes that occur over a small number of generations in natural populations often run counter to what is expected on the basis of the heritability of traits and the selective forces acting upon them. In Soay sheep, dark coat color is associated with large size, which is heritable and positively correlated with fitness, yet the frequency of dark sheep has decreased. This unexpected microevolutionary trend is explained by genetic linkage between the causal mutation underlying the color polymorphism and quantitative trait loci with antagonistic effects on size and fitness. As a consequence, homozygous dark sheep are large, but have reduced fitness relative to phenotypically indistinguishable dark heterozygotes and light sheep. This result demonstrates the importance of understanding the genetic basis of fitness variation when making predictions about the microevolutionary consequences of selection.

Phenotypic analysis of host-parasite interactions in lambs infected with Teladorsagia circumcincta.

Female Blackface lambs expected to exhibit genetic variability for resistance to gastrointestinal nematodes, were either exposed to continuous experimental infections of Teladorsagia circumcincta or were sham-dosed to monitor phenotypic responses to infection. As a measure of parasitism and host response, worm-eggs in faeces (faecal egg count, FEC) were counted over a 3-month period and worm burdens were ascertained at post-mortem. The host response to the infection was also measured by differential counts of white blood cells, anti-T. circumcincta IgA antibody levels and body weight. Results suggest that nematode abundance (mean number of parasites per host) and prevalence (proportion of infected animals) were maximal shortly after the beginning of infection (21 days p.i.) when virtually all the infected animals were shedding worm eggs. Increasing anti-T. circumcincta IgA antibody and eosinophil concentrations were associated with a reduction in total numbers of adult worms and an increase in the frequency of early L4s. The data also suggest that genetic selection for an enhanced anti-T. circumcincta IgA response might complement selection based on a reduced FEC as a strategy to select for resistance to gastrointestinal nematodes.

Fine mapping of the parthenocarpic fruit ( pat) mutation in tomato.

The parthenocarpic fruit ( pat) gene of tomato is a recessive mutation conferring parthenocarpy, which is the capability of a plant to set seedless fruits in the absence of pollination and fertilization. Parthenocarpic mutants offer a useful method to regulate fruit production and a suitable experimental system to study ovary and fruit development. In order to map the Pat locus, two populations segregating from the interspecific cross Lycopersicon esculentum x Lycopersicon pennellii were grown, and progeny plants were classified as parthenocarpic or wild-type by taking into account some characteristic aberrations affecting mutant anthers and ovules. Through bulk segregant analysis, we searched for both random and mapped AFLPs linked to the target gene. In this way, the Pat locus was assigned to the long arm of chromosome 3, as also confirmed by the analysis of a set of L. pennellii substitution and introgression lines. Afterwards, the Pat position was refined by using simple sequence repeats (SSRs) and conserved ortholog set (COS) markers mapping in the target region. The tightest COSs were converted into CAPS or SCAR markers. At present, two co-dominant SCAR markers encompassing a genetic window of 1.2 cM flank the Pat locus. Considering that these markers are orthologous to Arabidopsis genes, a positional cloning exploiting the tomato- Arabidopsis microsynteny seems to be a short-term objective.