Ian Mackay

From mutations to MAGIC: resources for gene discovery, validation and delivery in crop plants.

The dissection of gene-trait associations and its translation into practice through plant breeding is a central aspect of modern plant biology. The identification of genes underlying simply inherited traits has been very successful. However, the identification of gene-trait associations for complex (multi-genic) traits in crop plants with large, often polyploid genomes has been limited by the absence of appropriate genetic resources that allow quantitative trait loci (QTL) and causal genes to be identified and localised. There has also been a tendency for genetic resources to be developed in germplasm not directly relevant to the breeding community limiting effective implementation. In this review, we discuss approaches to mapping genes and the development of Multi-parent Advanced Generation Inter-cross (MAGIC) populations derived from breeder-relevant germplasm as a platform for a new generation of gene-trait analysis in crop species.

Genome-wide association mapping to candidate polymorphism resolution in the unsequenced barley genome

Although commonplace in human disease genetics, genome-wide association (GWA) studies have only relatively recently been applied to plants. Using 32 phenotypes in the inbreeding crop barley, we report GWA mapping of 15 morphological traits across ∼500 cultivars genotyped with 1,536 SNPs. In contrast to the majority of human GWA studies, we observe high levels of linkage disequilibrium within and between chromosomes. Despite this, GWA analysis readily detected common alleles of high penetrance. To investigate the potential of combining GWA mapping with comparative analysis to resolve traits to candidate polymorphism level in unsequenced genomes, we fine-mapped a selected phenotype (anthocyanin pigmentation) within a 140-kb interval containing three genes. Of these, resequencing the putative anthocyanin pathway gene HvbHLH1 identified a deletion resulting in a premature stop codon upstream of the basic helix-loop-helix domain, which was diagnostic for lack of anthocyanin in our association and biparental mapping populations. The methodology described here is transferable to species with limited genomic resources, providing a paradigm for reducing the threshold of map-based cloning in unsequenced crops.

Alpha-B crystallin gene (CRYAB) mutation causes dominant congenital posterior polar cataract in humans.

Congenital cataracts are an important cause of bilateral visual impairment in infants. In a four-generation family of English descent, we mapped dominant congenital posterior polar cataract to chromosome 11q22-q22.3. The maximum LOD score, 3.92 at recombination fraction 0, was obtained for marker D11S898, near the gene that encodes crystallin alpha-B protein (CRYAB). By sequencing the coding regions of CRYAB, we found in exon 3 a deletion mutation, 450delA, that is associated with cataract in this family. The mutation resulted in a frameshift in codon 150 and produced an aberrant protein consisting of 184 residues. This is the first report of a mutation, in this gene, resulting in isolated congenital cataract.

Genomewide Linkage Study in 1,176 Affected Sister Pair Families Identifies a Significant Susceptibility Locus for Endometriosis on Chromosome 10q26

Endometriosis is a common gynecological disease that affects up to 10% of women in their reproductive years. It causes pelvic pain, severe dysmenorrhea, and subfertility. The disease is defined as the presence of tissue resembling endometrium in sites outside the uterus. Its cause remains uncertain despite >50 years of hypothesis-driven research, and thus the therapeutic options are limited. Disease predisposition is inherited as a complex genetic trait, which provides an alternative route to understanding the disease. We seek to identify susceptibility loci, using a positional-cloning approach that starts with linkage analysis to identify genomic regions likely to harbor these genes. We conducted a linkage study of 1,176 families (931 from an Australian group and 245 from a U.K. group), each with at least two members--mainly affected sister pairs--with surgically diagnosed disease. We have identified a region of significant linkage on chromosome 10q26 (maximum LOD score [MLS] of 3.09; genomewide P = .047) and another region of suggestive linkage on chromosome 20p13 (MLS = 2.09). Minor peaks (with MLS > 1.0) were found on chromosomes 2, 6, 7, 8, 12, 14, 15, and 17. This is the first report of linkage to a major locus for endometriosis. The findings will facilitate discovery of novel positional genetic variants that influence the risk of developing this debilitating disease. Greater understanding of the aberrant cellular and molecular mechanisms involved in the etiology and pathophysiology of endometriosis should lead to better diagnostic methods and targeted treatments.

Association of the TSHR gene with Graves' disease: the first disease specific locus

The development of autoimmune thyroid disease (AITD) is associated with autoantibodies directed against the thyroid stimulating hormone receptor (TSHR). Previous studies have failed to demonstrate a consistent association between the TSHR and AITD, or any of its sub-phenotypes. In the present study, we analysed the linkage disequilibrium (LD) structure encompassing the TSHR, to identify LD ‘blocks’ and SNPs, which capture the majority of intra-block haplotype diversity. The haplotype tagging SNPs, plus all common SNPs reported in previous studies were genotyped in 1059 AITD Caucasian cases and 971 Caucasian controls. A haplotype, across two LD blocks, showed association (P<1 × 10−6, OR 1.7) with Graves’ disease (GD) but not autoimmune hypothyroidism (AIH). We replicated these findings by genotyping the most associated GD SNP, rs2268458, in a separate UK Caucasian cohort of 1366 AITD cases and 1061 controls (GD, P=2 × 10−6, OR 1.3; AIH, P=NS). These results in two independent Caucasian data sets suggest that the TSHR is the first replicated GD-specific locus meriting further fine mapping and functional analysis to identify the aetiological variants.

Population-Based Resequencing Reveals That the Flowering Time Adaptation of Cultivated Barley Originated East of the Fertile Crescent

Gene resequencing and association analysis present new opportunities to study the evolution of adaptive traits in crop plants. Here we apply these tools to an extensive set of barley accessions to identify a component of the molecular basis of the flowering time adaptation, a trait critical to plant survival. Using an association-based study to relate variation in flowering time to sequence-based polymorphisms in the Ppd-H1 gene, we identify a causative polymorphism (SNP48) that accounts for the observed variation in barley flowering time. This polymorphism also shows latitude-dependent geographical distribution, consistent with the expected clinal variation in phenotype with the nonresponsive form predominating in the north. Networks, genealogies, and phylogenetic trees drawn for the Ppd-H1 haplotypes reveal population structure both in wild barley and in domesticated barley landraces. The spatial distribution of these population groups indicates that phylogeographical analysis of European landraces can provide information relevant to the Neolithic spread of barley cultivation and also has implications for the origins of domesticated barley, including those with the nonresponsive ppd-H1 phenotype. Haplotypes containing the nonresponsive version of SNP48 are present in wild barley accessions, indicating that the nonresponsive phenotype of European landraces originated in wild barley. The wild accessions whose nonresponsive haplotypes are most closely similar to those of landraces are found in Iran, within a region suggested as an area for domestication of barley east of the Fertile Crescent but which has previously been thought to have contributed relatively little to the diversity of European cultivars.

Linkage disequilibrium in young genetically isolated Dutch population

The design and feasibility of genetic studies of complex diseases are critically dependent on the extent and distribution of linkage disequilibrium (LD) across the genome and between different populations. We have examined genomewide and region-specific LD in a young genetically isolated population identified in the Netherlands by genotyping approximately 800 Short Tandem Repeat markers distributed genomewide across 58 individuals. Several regions were analyzed further using a denser marker map. The permutation-corrected measure of LD was used for analysis. A significant (P<0.0004) relation between LD and genetic distance on a genomewide scale was found. Distance explained 4% of the total LD variation. For fine-mapping data, distance accounted for a larger proportion of LD variation (up to 39%). A notable similarity in the genomewide distribution of LD was revealed between this population and other young genetically isolated populations from Micronesia and Costa Rica. Our study population and experiment was simulated in silico to confirm our knowledge of the history of the population. High agreement was observed between results of analysis of simulated and empirical data. We conclude that our population shows a high level of LD similar to that demonstrated previously in other young genetic isolates. In Europe, there may be a large number of young genetically isolated populations that are similar in history to ours. In these populations, a similar degree of LD is expected and thus they may be effectively used for linkage or LD mapping.

The genetic diversity of UK, US and Australian cultivars of Triticum aestivum measured by DArT markers and considered by genome

The genetic diversity of UK, US and Australian wheat varieties over the period of modern plant breeding is estimated using diversity array technology markers. Diversity is assessed by both genetic distance between varieties, by AMOVA and as the volumes of multi-dimensional convex hulls estimated from principal co-ordinate analysis. At the whole genome level the three populations are genetically distinct; this is also true of the B genome. However, the US and Australian D genomes are found to occupy the same region of diversity space and the A genomes for these countries are partially overlapping. The use of high-density genotyping with a common marker set allows an unprecedented direct comparison between the diversities of the national populations, between individual genomes and the fluctuation of diversity over time. The highest genetic diversity amongst varieties is reported in the Australian population followed by the US, which in turn is more diverse than the UK. However the average diversity of loci is higher in the US set than in the Australian. Non-random fluctuations in genetic diversity over time are observed.

Association mapping of partitioning loci in barley

BackgroundAssociation mapping, initially developed in human disease genetics, is now being applied to plant species. The model species Arabidopsis provided some of the first examples of association mapping in plants, identifying previously cloned flowering time genes, despite high population sub-structure. More recently, association genetics has been applied to barley, where breeding activity has resulted in a high degree of population sub-structure. A major genotypic division within barley is that between winter- and spring-sown varieties, which differ in their requirement for vernalization to promote subsequent flowering. To date, all attempts to validate association genetics in barley by identifying major flowering time loci that control vernalization requirement (VRN-H1 and VRN-H2) have failed. Here, we validate the use of association genetics in barley by identifying VRN-H1 and VRN-H2, despite their prominent role in determining population sub-structure.ResultsBy taking barley as a typical inbreeding crop, and seasonal growth habit as a major partitioning phenotype, we develop an association mapping approach which successfully identifies VRN-H1 and VRN-H2, the underlying loci largely responsible for this agronomic division. We find a combination of Structured Association followed by Genomic Control to correct for population structure and inflation of the test statistic, resolved significant associations only with VRN-H1 and the VRN-H2 candidate genes, as well as two genes closely linked to VRN-H1 (HvCSFs1 and HvPHYC).ConclusionWe show that, after employing appropriate statistical methods to correct for population sub-structure, the genome-wide partitioning effect of allelic status at VRN-H1 and VRN-H2 does not result in the high levels of spurious association expected to occur in highly structured samples. Furthermore, we demonstrate that both VRN-H1 and the candidate VRN-H2 genes can be identified using association mapping. Discrimination between intragenic VRN-H1 markers was achieved, indicating that candidate causative polymorphisms may be discerned and prioritised within a larger set of positive associations. This proof of concept study demonstrates the feasibility of association mapping in barley, even within highly structured populations. A major advantage of this method is that it does not require large numbers of genome-wide markers, and is therefore suitable for fine mapping and candidate gene evaluation, especially in species for which large numbers of genetic markers are either unavailable or too costly.

Replication and extension studies of inflammatory bowel disease susceptibility regions confirm linkage to chromosome 6p ( IBD3 )

Inflammatory bowel disease (IBD) is a chronic inflammatory disease of the intestine, commonly diagnosed as either ulcerative colitis (UC) or Crohns disease (CD). Epidemiological studies have consistently shown that both genetic and environmental factors influence the pathogenesis of IBD. A number of genome scans have been conducted in cohorts of IBD families with affected sibling pairs (ASPs) to identify chromosomal regions that harbour IBD susceptibility genes. Several putative linked loci have been identified, including two loci on chromosomes 16 and 12, IBD1 and IBD2, which have subsequently been replicated by independent region-specific studies. We have conducted both a replication study on another linkage region, chromosome 6p (IBD3), and extension studies on two other regions, chromosomes 3p and 7q. Microsatellite markers across each region were genotyped in 284 IBD ASPs from 234 families. A nonparametric peak multipoint LOD score of 3.0 was observed near D6S291, replicating the previous linkage to chromosome 6p (IBD3). Nominal evidence of linkage was observed at both the 3p and 7q regions.