MIDAS: software for analysis and visualisation of interallelic disequilibrium between multiallelic markers

IRGM gene polymorphisms and risk of gastric cancer

OBJECTIVE

The study aimed to assess the possible association of polymorphisms in the autophagy gene IRGM (rs13361189 and rs4958847) with the risk of gastric cancer.

METHODS

A total of 102 patients with gastric adenocarcinoma, 52 with chronic gastritis and 351 healthy controls were included in this study. IRGM allelic variants were genotyped by quantitative real-time polymerase chain reaction. The association between polymorphisms and gastric cancer risk was estimated by odds ratios (OR) and 95% confidence interval (CI).

RESULTS

A significant difference was found for rs4958847 A allele. Carriers of the A allele were protected against gastric cancer (OR = 0.58, 95% CI 0.35-0.97, P = 0.038). Moreover, the presence of this allele seems to play an important role in decreasing the risk for the intestinal type of gastric cancer (OR = 0.47, 95% CI 0.23-0.94, P = 0.03). In contrast, the rs13361189 IRGM polymorphism was not associated with susceptibility to gastric cancer. None of the targeted polymorphisms were associated with chronic gastritis.

CONCLUSION

IRGM rs4958847 polymorphism influences susceptibility to gastric cancer, mainly for the intestinal type.

Adaptive introgression across species boundaries in Heliconius butterflies

It is widely documented that hybridisation occurs between many closely related species, but the importance of introgression in adaptive evolution remains unclear, especially in animals. Here, we have examined the role of introgressive hybridisation in transferring adaptations between mimetic Heliconius butterflies, taking advantage of the recent identification of a gene regulating red wing patterns in this genus. By sequencing regions both linked and unlinked to the red colour locus, we found a region that displays an almost perfect genotype by phenotype association across four species, H. melpomene, H. cydno, H. timareta, and H. heurippa. This particular segment is located 70 kb downstream of the red colour specification gene optix, and coalescent analysis indicates repeated introgression of adaptive alleles from H. melpomene into the H. cydno species clade. Our analytical methods complement recent genome scale data for the same region and suggest adaptive introgression has a crucial role in generating adaptive wing colour diversity in this group of butterflies.

Hybridisation occurs between many animal species, however its evolutionary relevance is still a matter of great debate. While some argue that hybridisation leads to maladaptive gene combinations, and therefore to an evolutionary dead end, others consider interspecific hybridisation as a process with great potential to fuel evolution. We examine this question by exploring the origins of red wing colouration, a trait under natural selection, in the adaptive radiation of closely related species of Heliconius butterflies. By sequencing genetic regions both linked and unlinked to the red wing pattern locus, we found experimental evidence supporting multiple hybridisation events that have mediated the acquisition of colour adaptations from H. melpomene to H. timareta. This introgression has allowed H. timareta to colonise new fitness peaks in the Müllerian mimicry landscape. In this way, our results support the idea that interspecific hybridisation in animals constitutes a source of genetic variation that promotes diversification.

Genetic variation and hitchhiking between structurally polymorphic Alu insertions and HLA-A, -B, and -C alleles and other retroelements within the MHC class I region

We investigated structurally polymorphic Alu insertions (POALINs) at five loci in the major histocompatibility complex (MHC) class I genomic region to determine their allele and haplotype frequencies and associations with the human leukocyte antigen (HLA)-A, -B, and -C genes in three populations, the Australian Caucasians, Japanese, and African Americans. The POALINs varied in allelic frequency between 0% and 42.3% with significant differences between populations at three of the five loci. The linkage disequilibrium (LD) between Alu insertions and the HLA-A, -B, or -C alleles and previously published polymorphic retroelements (four SVA and human endogenous retrovirus type 9 (HERVK9) loci) within the class I region of the MHC were calculated in pairwise analyses of haplotypes to show strong allelic associations and possible crossing-over events between some loci. Each POALIN was in significant LD with a variety of HLA-A, -B, or -C two-digit alleles probably as a result of hitchhiking. The POALINs helped to further stratify the HLA-A:B:C haplotypes into different POALIN:HLA-A:B:C haplotype frequencies. Of the multilocus haplotype analyses, the seven- and eight-locus haplotypes showed the largest number of differences between the populations, and fewer matched haplotypes between populations that ranged, for example, from 49% for HLA-B:HLA-A haplotypes to 7% for AluMICB:HLA-B:HLA-C:AluTF:AluHJ:HLA-A:AluHG:AluTF haplotypes in the Japanese. This comparative study of multilocus POALINs in the HLA class I region of three ethnic populations shows that POALINs alone or together with the HLA class I alleles and other retroelements are informative ancestral markers for assessing the interrelationship of HLA class I haplotype lineages, LD, and genetic diversity within the same and/or different populations.

Molecular and population analysis of natural selection on the human haptoglobin duplication

Haptoglobin binds free haemoglobin that prevents oxidative damage produced by haemolysis. There is a copy number variant (CNV) in the haptoglobin gene (HP) consisting of two alleles, Hp1 (no duplication), and Hp2 (1.7kb duplication involving two exons). The spread of the Hp2 allele is believed to have taken place under selective pressures conferred by malaria resistance. However, molecular evidence is lacking and Hp did not emerge in genomewide SNPs surveys for evidence of selection. In Europe, there is geographical constancy of Hp2 frequency, indicative of absence of clinal pressures and that modern day European alleles represent a “snapshot” of their out-of-Africa migrations. In this work we test for signatures of natural selection acting on the Hp CNV in a sample from the UK population (Avon Longitudinal Study of Parents and Children, ALSPAC). We present here heterozygosity decay, pairwise F_ST values observed between ALSPAC and 301 populations from all five populated continents, extended haplotype homozygosity analyses involving the CNV and 80 SNPs surrounding the CNV ∼500kb in each direction, and linkage disequilibrium and pairwise haplotypic analyses involving 160 SNPs on chromosome 16q22.1. Taken together, our results represent the first molecular analysis of natural selection in the Hp CNV genetic region.

Reduced interferon (IFN)-α conditioned by IFNA2 (-173) and IFNA8 (-884) haplotypes is associated with enhanced susceptibility to severe malarial anemia and longitudinal all-cause mortality

Severe malarial anemia (SMA) is a leading cause of pediatric morbidity and mortality in holoendemic Plasmodium falciparum transmission areas. Although dysregulation in cytokine production is an important etiology of SMA, the role of IFN-α in SMA has not been reported. As such, we investigated the relationship between IFN-α promoter polymorphisms [i.e., IFNA2 (A-173T) and IFNA8 (T-884A)], SMA, and functional changes in IFN-α production in children (n = 663; <36 months) residing in a holoendemic P. falciparum transmission region of Kenya. Children with SMA had lower circulating IFN-α than malaria-infected children without severe anemia (P = 0.025). Multivariate logistic regression analyses revealed that heterozygosity at -884 (TA) was associated with an increased risk of SMA [OR 2.80 (95 % CI 1.22-6.43); P = 0.015] and reduced IFN-α relative to wild type (TT; P = 0.038). Additional analyses demonstrated that carriage of the -173T/-884A (TA) haplotype was associated with increased susceptibility to SMA [OR 3.98 (95 % CI 1.17-13.52); P = 0.026] and lower IFN-α (P = 0.031). Follow-up of these children for 36 months revealed that carriers of TA haplotype had greater all-cause mortality than non-carriers (P < 0.001). Generation of reporter constructs showed that the IFNA8 wild-type -884TT exhibited higher levels of luciferase expression than the variant alleles (P < 0.001). Analyses of malaria-associated inflammatory mediators demonstrated that carriers of TA haplotype had altered production of IL-1β, MIG, and IL-13 compared to non-carriers (P < 0.050). Thus, variation at IFNA2 -173 and IFNA8 -884 conditions reduced IFN-α production, and increased susceptibility to SMA and mortality.

Dopamine D4 and D5 receptor gene variant effects on clozapine response in schizophrenia: replication and exploration

OBJECTIVES

This study aimed to: 1) replicate previously reported associations between dopamine D4 receptor gene (DRD4) polymorphisms and antipsychotic (AP) response in a clozapine (CLZ) response sample; and 2) explore possible associations of polymorphisms across dopamine D5 receptor gene (DRD5) as well as other DRD4 regions.

METHODS

DRD4 exon III 48-bp, intron I (G)(n), and 120-bp repeat polymorphisms, and three DRD4 single nucleotide polymorphisms (SNPs); and DRD5 (CA/CT/GT)(n) microsatellite and four DRD5 SNPs were assessed using standard genotyping and statistical procedures.

RESULTS

We report evidence, which does not survive correction for multiple testing, supporting previous DRD4 findings. Findings of interest include the 120-bp 1-copy allele, intron I (G)(n) 142-bp/140-bp genotype, and exon III 4R allele with CLZ response. All DRD5 tests were negative.

CONCLUSIONS

Overall, these results suggest a possible minor contribution of DRD4 variants, but not DRD5 variants, towards the AP/CLZ response phenotype.

Common genetic variation in the IGF1 associates with maximal force output

PURPOSE

We clarified the effect of insulin-like growth factor-1 (IGF1), IGF-binding protein-3 (IGFBP3), interleukin-6 (IL6), and its receptor (IL6R) gene variants on muscular and aerobic performance, body composition, and on circulating levels of IGF-1 and IL-6. Single nucleotide polymorphisms (SNPs) may, in general, influence gene regulation or its expression, or the structure and function of the corresponding protein, and modify its biological effects. IGF-1 is involved in the anabolic pathways of skeletal muscle. IL-6 plays an important role in muscle energy homeostasis during strenuous physical exercise.

METHODS

Eight hundred forty-one healthy Finnish male subjects of Caucasian origin were genotyped for IGF1 (rs6220 and rs7136446), IGFBP3 (rs2854744), IL6 (rs1800795), and IL6R (rs4537545) SNPs, and studied for associations with maximal force of leg extensor muscles, maximal oxygen consumption, body fat percent, and IGF-1 and IL-6 levels. Analytic methods included dynamometer, bicycle ergometer, bioimpedance, ELISA, and polymerase chain reaction assays.

RESULTS

All investigated SNPs conformed to Hardy-Weinberg equilibrium with allele frequencies validated against CEU population. Genotype CC of rs7136446 associated with higher body fat and increased maximal force production. Genotype CC of the IGFBP3 SNP rs2854744 and TT genotype of the IL6R SNP rs4537545 associated with higher IL-6 levels. In logistic regression analysis, allele C of the rs2854744 decreased odds for lower body fat. None of the studied SNPs associated with aerobic performance.

CONCLUSIONS

Our data suggest that common variation in the IGF1 gene may affect maximal force production, which can be explained by the role of IGF-1 in the anabolic pathways of muscle and neurotrophy. Variations in the IGF1 and IGFBP3 gene may result in higher body fat and be related to alterations of IGF-1-mediated tissue growth.

Analytical methods for immunogenetic population data

In this chapter, we describe analyses commonly applied to immunogenetic population data, along with software tools that are currently available to perform those analyses. Where possible, we focus on tools that have been developed specifically for the analysis of highly polymorphic immunogenetic data. These analytical methods serve both as a means to examine the appropriateness of a dataset for testing a specific hypothesis, as well as a means of testing hypotheses. Rather than treat this chapter as a protocol for analyzing any population dataset, each researcher and analyst should first consider their data, the possible analyses, and any available tools in light of the hypothesis being tested. The extent to which the data and analyses are appropriate to each other should be determined before any analyses are performed.

Species discrimination, population structure and linkage disequilibrium in Eucalyptus camaldulensis and Eucalyptus tereticornis using SSR markers

Eucalyptus camaldulensis and E. tereticornis are closely related species commonly cultivated for pulp wood in many tropical countries including India. Understanding the genetic structure and linkage disequilibrium (LD) existing in these species is essential for the improvement of industrially important traits. Our goal was to evaluate the use of simple sequence repeat (SSR) loci for species discrimination, population structure and LD analysis in these species. Investigations were carried out with the most common alleles in 93 accessions belonging to these two species using 62 SSR markers through cross amplification. The polymorphic information content (PIC) ranged from 0.44 to 0.93 and 0.36 to 0.93 in E. camaldulensis and E. tereticornis respectively. A clear delineation between the two species was evident based on the analysis of population structure and species-specific alleles. Significant genotypic LD was found in E. camaldulensis, wherein out of 135 significant pairs, 17 pairs showed r²≥0.1. Similarly, in E. tereticornis, out of 136 significant pairs, 18 pairs showed r²≥0.1. The extent of LD decayed rapidly showing the significance of association analyses in eucalypts with higher resolution markers. The availability of whole genome sequence for E. grandis and the synteny and co-linearity in the genome of eucalypts, will allow genome-wide genotyping using microsatellites or single nucleotide polymorphims.

Functional promoter haplotypes of interleukin-18 condition susceptibility to severe malarial anemia and childhood mortality

Severe malarial anemia (SMA) is a leading cause of morbidity and mortality in children residing in regions where Plasmodium falciparum transmission is holoendemic. Although largely unexplored in children with SMA, interleukin-18 (IL-18) is important for regulating innate and acquired immunity in inflammatory and infectious diseases. As such, we selected two functional single-nucleotide polymorphisms (SNPs) in the IL-18 promoter (-137G→C [rs187238] and -607C→A [rs1946518]) whose haplotypes encompass significant genetic variation due to the presence of strong linkage disequilibrium among these variants. The relationship between the genotypes/haplotypes, SMA (hemoglobin [Hb], <5.0 g/dl], and longitudinal clinical outcomes were then investigated in Kenyan children (n = 719). Multivariate logistic regression analyses controlling for age, gender, sickle cell trait, glucose-6-phosphate dehydrogenase (G6PD) deficiency, HIV-1, and bacteremia revealed that carriage of the -607AA genotype was associated with protection against SMA (odds ratio [OR] = 0.440 [95% confidence interval {CI} = 0.21 to 0.90], P = 0.031) in children with acute infection. In contrast, carriers of the -137G/-607C (GC) haplotype had increased susceptibility to SMA (OR = 2.050 [95% CI = 1.04 to 4.05], P = 0.039). Measurement of IL-18 gene expression in peripheral blood leukocytes demonstrated that elevated IL-18 transcripts were associated with reduced hemoglobin concentrations (ρ = -0.293, P = 0.010) and that carriers of the "susceptible" GC haplotype had elevated IL-18 transcripts (P = 0.026). Longitudinal investigation of clinical outcomes over a 3-year follow-up period revealed that carriers of the rare CC haplotype (∼1% frequency) had 5.76 times higher mortality than noncarriers (P = 0.001). Results presented here demonstrate that IL-18 promoter haplotypes that condition elevated IL-18 gene products during acute infection are associated with increased risk of SMA. Furthermore, carriage of the rare CC haplotype significantly increases the risk of childhood mortality.

Association of heme oxygenase-1 gene polymorphisms with essential hypertension and blood pressure in the Chinese Han population

Many studies have suggested that heme oxygenase-1 (HMOX1) might be implicated in blood pressure (BP) regulation. We hypothesized that this gene might be responsible for the variation of susceptibility to essential hypertension (EH) and BP and investigated three polymorphisms in HMOX1 (i.e., the (GT)n repeat in the HMOX1 promoter and single-nucleotide polymorphisms (SNPs) rs2071746 and rs2071749) in population-based samples of 789 Han Chinese from Xinjiang, China. The GT repeat numbers ≥33, 27-32, and <27 in the HMOX1 promoter were defined as long (L), middle (M), and short (S) alleles, respectively. The participants carrying SS or SM genotype were sorted into one group, and the participants carrying SL, MM, ML, or LL genotype were sorted into another group. The (GT)n repeat in the HMOX1 promoter showed significant association with EH and BP, whereas SNPs rs2071746 and rs2071749 did not. Compared with the SS+SM (GT)n group, the MM+SL+ML+LL (GT)n group had a lower risk of EH (adjusted odds ratio, 0.58; 95% confidence interval, 0.41-0.82; p = 0.002) and lower systolic (128.3 ± 1.3 vs. 132.2 ± 1.0 mm Hg, adjusted p = 0.014) and diastolic BPs (81.7 ± 0.8 vs. 84.5 ± 0.7 mm Hg, adjusted p = 0.009). These findings provide the first genetic evidence for the role of the (GT)n repeat in the HMOX1 promoter in the susceptibility to human EH and the variation of BP.

ABCB1 genotype and CSF beta-amyloid in Alzheimer disease

The ABCB1 gene, coding for the efflux transporter P-glycoprotein (PGP), is a candidate gene for Alzheimer disease (AD). P-glycoprotein is heavily expressed at the blood-brain barrier, where it mediates the efflux of β-amyloid (Aβ) from the brain. In this study, we investigated a possible association between 2 common ABCB1 polymorphisms, G2677T/A (Ala893Ser/Thr) and C3435T, AD, and cerebrospinal fluid (CSF) levels of Aβ. No strong evidence for association was found.

Gametic phase disequilibrium between the syntenic multiallelic HTG4 and HMS3 markers widely used for parentage testing in Thoroughbred horses

Validation of parentage and horse breed registries through DNA typing relies on estimates of random match probabilities with DNA profiles generated from multiple polymorphic loci. Of the twenty-seven microsatellite loci recommended by the International Society for Animal Genetics for parentage testing in Thoroughbred horses, eleven are located on five chromosomes. An important aspect in determining combined exclusion probabilities is the ascertainment of the genetic linkage status of syntenic markers, which may affect reliable use of the product rule in estimating random match probabilities. In principle, linked markers can be in gametic phase disequilibrium (GD). We aimed at determining the extent, by frequency and strength, of GD between the HTG4 and HMS3 multiallelic loci, syntenic on chromosome 9. We typed the qualified offspring (n (1) = 27; n (2) = 14) of two Quarter Bred stallions (registered by the Brazilian Association of Quarter Horse Breeders) and 121 unrelated horses from the same breed. In the 41 informative meioses analyzed, the frequency of recombination between the HTG4 and HMS3 loci was 0.27. Consistent with genetic map distances, this recombination rate does not fit to the theoretical distribution for independently segregated markers. We estimated sign-based D' coefficients as a measure of GD, and showed that the HTG4 and HMS3 loci are in significant, yet partial and weak, disequilibrium, with two allele pairs involved (HTG4 M/HMS3 P, D'(+) = 0.6274; and HTG4 K/HMS3 P, D'(-) = -0.6096). These results warn against the inadequate inclusion of genetically linked markers in the calculation of combined power of discrimination for Thoroughbred parentage validation.

Thymidylate synthase expression and genotype have no major impact on the clinical outcome of colorectal cancer patients treated with 5-fluorouracil

BACKGROUND AND OBJECTIVES

Thymidylate synthase (TS) expression levels appear to be related to response to 5-fluorouracil-(5-FU)-based chemotherapy in colorectal cancer (CRC) patients. Three polymorphisms have been proposed as modulators of TS expression: a tandemly repeated sequence (2R/3R) in the 5' UTR, a SNP (G>C) within the 3R allele and a 6bp deletion in the 3' UTR. To evaluate the influence of TS expression and polymorphisms on clinical outcome of 5-FU-treated patients we performed a comprehensive genetic analysis on 63 CRC patients.

METHODS

TS expression levels were analyzed in normal and tumor tissues. TS coding sequence and UTR polymorphisms were investigated on DNA from normal tissue. LOH analysis was performed to determine tumor genotype.

RESULTS

A difference in disease-free survival (DFS), although not statistically significant, was observed between high and low mRNA expression levels: patients with low levels showed longer DFS. The 2R2R genotype showed significantly lower expression than the 3R3R and 2R3R genotypes in normal tissue. No other TS polymorphism was associated with mRNA expression or clinical outcome.

CONCLUSIONS

The results obtained in this pilot study indicate that the number of 5' UTR repeats is the major genetic determinant of TS expression. The lack of association with other polymorphisms might be partially explained by the existence of linkage disequilibrium in the TS gene. Our data support the growing evidence that TS control may require multiple mechanisms acting in close coordination with one another and suggest that TS genotyping alone in tumor samples is not sufficient to accurately predict response to 5-FU.

Support for a bipolar affective disorder susceptibility locus on chromosome 12q24.3

OBJECTIVE

Linkage and association studies of bipolar affective disorder (BAD) point out chromosome 12q24 as a region of interest.

METHODS

To investigate this region further, we conducted an association study of 22 DNA markers within a 1.14 Mb region in a Danish sample of 166 patients with BAD and 311 control individuals. Two-hundred and four Danish patients with schizophrenia were also included in the study.

RESULTS

We observed highly significant allelic and genotypic association between BAD and two highly correlated markers. The risk allele of both markers considered separately conferred an odds ratio of 2 to an individual carrying one risk allele and an odds ratio of 4 for individuals carrying both risk alleles assuming an additive genetic model. These findings were supported by the haplotype analysis. In addition, we obtained a replication of four markers associated with BAD in an earlier UK study. The most significantly associated marker was also analyzed in a Scottish case-control sample and was earlier associated with BAD in the UK cohort. The association of that particular marker was strongly associated with BAD in a meta-analysis of the Danish, Scottish and UK sample (P=0.0003). The chromosome region confined by our most distant markers is gene-poor and harbours only a few predicted genes. This study implicates the Slynar locus. We confirmed one annotated Slynar transcript and identified a novel transcript in human brain cDNA.

CONCLUSION

This study confirms 12q24.3 as a region of functional importance in the pathogenesis of BAD and highlights the importance of focused genotyping.

Informativeness of a novel multiallelic marker-set comprising an F8 intron 21 and three tightly linked loci for haemophilia A carriership analysis

The extraordinary heterogeneous nature of the deleterious mutations in the F8 gene that lead to functional deficiency of clotting factor VIII in haemophilia A makes routine direct mutation profiling difficult. When direct mutation analysis cannot be performed or a causative/candidate mutation is not found, a second-line approach to track the defective F8 gene within at-risk families is linkage genetic analysis with, tried-and-tested, F8-intragenic and/or extragenic non-recombining multiallelic short tandem repeats (STR). Although several typing STR loci within and around F8 have been described, there is need for improving assessment, because the combined informativeness of available assays rarely reaches 100%. Here, we characterized a newly identified 0.28 cM-resolution marker-set, consisting of a dinucleotide STR located on F8 intron 21 (F8Int21; [AC](n)) and three extragenic tetranucleotide STR located on GAB3 intron 1 (GAB3Int1; [TAAA](n)) and TMLHE intron 1 (TMLHEInt1.1; [GAAA](n) and TMLHEInt1.3; [ATTC](n)). Heterozygosity rates determined in 100 unrelated females ranged from 0.25 (GAB3Int1) to 0.63 (F8Int21). The set rendered a combined informativeness of 0.91 for at least one marker and 0.60 for a minimum of two loci, with at least one F8-intragenic. Multiallelic interlocus non-random association analysis revealed that GAB3Int1 is not in significant gametic disequilibrium (GD) with F8Int21, F8Int9.2, TMLHEInt1.3 or TMLHEInt1.1. Gametic disequilibrium breakdown attests historical recombination between GAB3Int1 and the F8 gene. Through computational analysis of reference assembly sequence data, we note in the GD breakdown region and in the F8 gene a higher than average density of the 13-mer CCNCCNTNNCCNC consensus motif, commonly associated with recombination hotspots.

A microsatellite-based analysis for the detection of selection on BTA1 and BTA20 in northern Eurasian cattle (Bos taurus) populations

BACKGROUND

Microsatellites surrounding functionally important candidate genes or quantitative trait loci have received attention as proxy measures of polymorphism level at the candidate loci themselves. In cattle, selection for economically important traits is a long-term strategy and it has been reported that microsatellites are linked to these important loci.

METHODS

We have investigated the variation of seven microsatellites on BTA1 (Bos taurus autosome 1) and 16 on BTA20, using bovine populations of typical production types and horn status in northern Eurasia. Genetic variability of these loci and linkage disequilibrium among these loci were compared with those of 28 microsatellites on other bovine chromosomes. Four different tests were applied to detect molecular signatures of selection.

RESULTS

No marked difference in locus variability was found between microsatellites on BTA1, BTA20 and the other chromosomes in terms of different diversity indices. Average D' values of pairwise syntenic markers (0.32 and 0.28 across BTA 1 and BTA20 respectively) were significantly (P < 0.05) higher than for non-syntenic markers (0.15). The Ewens-Watterson test, the Beaumont and Nichol's modified frequentist test and the Bayesian FST-test indicated elevated or decreased genetic differentiation, at SOD1 and AGLA17 markers respectively, deviating significantly (P < 0.05) from neutral expectations. Furthermore, lnRV, lnRH and lnRtheta' statistics were used for the pairwise population comparison tests and were significantly less variable in one population relative to the other, providing additional evidence of selection signatures for two of the 51 loci. Moreover, the three Finnish native populations showed evidence of subpopulation divergence at SOD1 and AGLA17. Our data also indicate significant intergenic linkage disequilibrium around the candidate loci and suggest that hitchhiking selection has played a role in shaping the pattern of observed linkage disequilibrium.

CONCLUSION

Hitchhiking due to tight linkage with alleles at candidate genes, e.g. the POLL gene, is a possible explanation for this pattern. The potential impact of selective breeding by man on cattle populations is discussed in the context of selection effects. Our results also suggest that a practical approach to detect loci under selection is to simultaneously apply multiple neutrality tests based on different assumptions and estimations.

Shifting dispersal modes at an expanding species' range margin

While it is generally recognized that noncontiguous (long-distance) dispersal of small numbers of individuals is important for range expansion over large geographic areas, it is often assumed that colonization on more local scales proceeds by population expansion and diffusion dispersal (larger numbers of individuals colonizing adjacent sites). There are few empirical studies of dispersal modes at the front of expanding ranges, and very little information is available on dispersal dynamics at smaller geographic scales where we expect contiguous (diffusion) dispersal to be prevalent. We used highly polymorphic genetic markers to characterize dispersal modes at a local geographic scale for populations at the edge of the range of a newly invasive grass species (Brachypodium sylvaticum) that is undergoing rapid range expansion in the Pacific Northwest of North America. Comparisons of Bayesian clustering of populations, patterns of genetic diversity, and gametic disequilibrium indicate that new populations are colonized ahead of the invasion front by noncontiguous dispersal from source populations, with admixture occurring as populations age. This pattern of noncontiguous colonization was maintained even at a local scale. Absence of evidence for dispersal among adjacent pioneer sites at the edge of the expanding range of this species suggests that pioneer populations undergo an establishment phase during which they do not contribute emigrants for colonization of neighbouring sites. Our data indicate that dispersal modes change as the invasion matures: initial colonization processes appear to be dominated by noncontiguous dispersal from only a few sources, while contiguous dispersal may play a greater role once populations become established.

Genetic evidence for hybrid trait speciation in heliconius butterflies

Homoploid hybrid speciation is the formation of a new hybrid species without change in chromosome number. So far, there has been a lack of direct molecular evidence for hybridization generating novel traits directly involved in animal speciation. Heliconius butterflies exhibit bright aposematic color patterns that also act as cues in assortative mating. Heliconius heurippa has been proposed as a hybrid species, and its color pattern can be recreated by introgression of the H. m. melpomene red band into the genetic background of the yellow banded H. cydno cordula. This hybrid color pattern is also involved in mate choice and leads to reproductive isolation between H. heurippa and its close relatives. Here, we provide molecular evidence for adaptive introgression by sequencing genes across the Heliconius red band locus and comparing them to unlinked wing patterning genes in H. melpomene, H. cydno, and H. heurippa. 670 SNPs distributed among 29 unlinked coding genes (25,847bp) showed H. heurippa was related to H. c. cordula or the three species were intermixed. In contrast, among 344 SNPs distributed among 13 genes in the red band region (18,629bp), most showed H. heurippa related with H. c. cordula, but a block of around 6,5kb located in the 3' of a putative kinesin gene grouped H. heurippa with H. m. melpomene, supporting the hybrid introgression hypothesis. Genealogical reconstruction showed that this introgression occurred after divergence of the parental species, perhaps around 0.43Mya. Expression of the kinesin gene is spatially restricted to the distal region of the forewing, suggesting a mechanism for pattern regulation. This gene therefore constitutes the first molecular evidence for adaptive introgression during hybrid speciation and is the first clear candidate for a Heliconius wing patterning locus.

Polymorphic major histocompatibility complex class II Alu insertions at five loci and their association with HLA-DRB1 and -DQB1 in Japanese and Caucasians

We investigated polymorphic Alu insertion (POALIN) frequencies at five loci in the major histocompatibility complex (MHC) class II genomic region to determine their allele and haplotype frequencies and associations with the human leukocyte antigen (HLA)-DRB1 and -DQB1 genes for 100 Japanese, 174 Australian Caucasians and 67 HLA reference cell lines obtained from different ethnic groups. The POALINs varied in frequency between 11% and 57% with significant differences between the Japanese and Caucasians at three loci. One POALIN locus deviated significantly from Hardy-Weinberg equilibrium (HWE) and four POALIN loci were in significant linkage disequilibrium and had a high percentage association with a variety of HLA-DRB1 or -DQB1 two-digit alleles. Inferred haplotype analysis among two-locus, five-locus and seven-locus haplotype structures showed maximum differences between the Japanese and Caucasians with the seven-locus haplotypes. The most common multilocus haplotype in Caucasians was DRB1*1501/DQB1*0602/AluDQ1/AluDRB1/AluORF10/AluDPB2 (6.7%), whereas the second most common allele HLA-DRB1*15 (17.5%) in Japanese was associated with three or four Alu insertions. The HLA class II POALINs also differentiated within and between HLA-DRB1 super-haplotypes DR1, DR8, DR51, DR52 and DR53. This is the first comparative population study of multilocus POALINs in the HLA class II region, which shows that POALINs whether investigated alone or together with the HLA class II alleles are informative genetic markers for the identification of allele and haplotype lineages and variations within the same and/or different populations.

Sex-specific population structure, natural selection, and linkage disequilibrium in a wild bird population as revealed by genome-wide microsatellite analyses

Background

Sexual dimorphism in ecologically important traits is widespread, yet the differences in the genomic architecture between the two sexes are largely unexplored. We employed a genome-wide multilocus approach to examine the sexual differences in population subdivision, natural selection and linkage disequilibrium (LD) in a wild Siberian jay (Perisoreus infaustus) population, using genotypes at a total of 107 autosomal and Z-chromosomal microsatellites.

Results

Mean observed heterozygosity was significantly higher in females (H_O = 0.567) than in males (H_O = 0.532), and autosomal markers (H_O = 0.561) were more variable than Z-chromosomal markers (H_O = 0.512). Genetic differentiation (F_ST = 0.002, P < 0.05) between the two sexes was low but significant and males were on average significantly more genetically related to each other than females. Genomescan analyses revealed that 3 out of 101 (3%) autosomal loci were under directional selection, while 4 out of 6 (67%) Z-chromosomal markers were indicated to be under balancing selection. This suggests a significantly greater but contrasting selection force on the Z-chromosome in comparison to autosomes, which is consistent with an overall significantly (P < 0.05) lower F_STvalue for Z-chromosomal (-0.014, 95% CI: -0.025 - -0.011) than for the autosomal loci (0.003, 95% CI: 0.001 - 0.004). Analysis of syntenic marker pairs revealed high levels of LD in both sexes but significantly (P < 0.05) lower levels of LD in the females both on autosomes and Z-chromosome, probably due to the higher rate of dispersal and the higher recombination rates on autosomes, as well as the pseudoautosomal markers. In both sexes LD decayed rapidly with genetic distance in a similar fashion on autosomes, while a more rapid decay of LD in Z-chromosome was detected in females than in males.

Conclusion

We conclude that there are many clear differences in genomic architecture between the sexes studied here which can be at least partly understood in the light of higher dispersal rate of females as compared to males and the unusual structure of the Z-chromosome of the species.

Signatures of natural selection and coevolution between killer cell immunoglobulin-like receptors (KIR) and HLA class I genes

Natural killer (NK) cells are lymphocytes of the innate immune system. In humans, NK cell activities are partly controlled by the diverse killer immunoglobulin-like receptor (KIR) gene family. The importance of NK cells in both immunity to infection and reproduction makes KIR strong candidates for genes undergoing dynamic evolution in the human genome. Using high-resolution allelic typing, we investigated the potential role of natural selection in the diversification of KIR in the Irish population. Higher diversity than expected is observed at several loci, consistent with a history of balancing selection acting to maintain several allelic variants at high frequency in the population. KIR diversity is enhanced further at the haplotype level with functional polymorphisms at KIR2DL4, KIR3DL1 and KIR2DS4 defining nine 'core' haplotypes. Analysis of these core haplotypes in combination with human leukocyte antigen (HLA) class I ligands revealed several nonrandom associations. In particular, the KIR:HLA association for the core haplotype defined by KIR3DL1(*)01502 was female specific and a likely consequence of negative selection acting against KIR3DL1(*)01502 on an HLA-C1/C1 background. Many of the associations between KIR and HLA in the Irish differ from those previously reported, which argues against universal selective pressures for specific KIR:HLA combinations in diverse human populations.

Polymorphic SVA retrotransposons at four loci and their association with classical HLA class I alleles in Japanese, Caucasians and African Americans

Polymorphic insertion frequencies of the retrotransposons known as the "SVA" elements were investigated at four loci in the MHC class I genomic region to determine their allele and haplotype frequencies and associations with the HLA-A, -B or -C genes for 100 Japanese, 100 African Americans, 174 Australian Caucasians and 66 reference cell lines obtained from different ethnic groups. The SVA insertions representing different subfamily members varied in frequency between none for SVA-HF in Japanese and 65% for SVA-HB in Caucasians or African Americans with significant differences in frequencies between the three populations at least at three loci. The SVA loci were in Hardy-Weinberg equilibrium except for the SVA-HA locus which deviated significantly in African Americans and Caucasians possibly because of a genomic deletion of this locus in individuals with the HLA-A*24 allele. Strong linkage disequilibria and high percentage associations between the human leucocyte antigen (HLA) class I gene alleles and some of the SVA insertions were detected in all three populations in spite of significant frequency differences for the SVA and HLA class I alleles between the three populations. The highest percentage associations (>86%) were between SVA-HB and HLA-B*08, -B*27, -B*37 to -B*41, -B*52 and -B*53; SVA-HC and HLA-B*07; SVA-HA and HLA-A*03, -A*11 and -A*30; and SVA-HF and HLA-A*03 and HLA-B*47. From pairwise associations in the three populations and the homozygous cell line results, it was possible to deduce the SVA and HLA class I allelic combinations (haplotypes), population differences and the identity by descent of several common HLA-A allelic lineages.

Interlocus non-random association of multiallelic polymorphisms spanning the coagulation factor VIII gene on human chromosome distalmost Xq28

The most common severe hereditary bleeding disorder phenotype in humans, the coagulation factor VIII (F8) deficiency haemophilia A (HEMA), maps on Xq28 band, a region that comprises 11.7% of genes and 14.2% of phenotypes on X chromosome. Information about the distribution and extent of gametic disequilibrium (GD) covering the F8 gene is scarce, despite its relevance for linkage and association studies. The aim of this study was to determine the patterns, by frequency and strength, of non-random multiallelic interallelic associations between two-locus combinations of seven microsatellite loci (REN90833, F8Int25.2, F8Int22, F8Int13.2, HEMA154311.3, TMLHEInt5 and HEMA154507.3, in that physical order) spanning 0.813 Mb on distalmost Xq28. We measured sign-based interallelic D' coefficients in 106 men and in 100 women drawn from a single unrelated Brazilian population. Significance and patterns of GD using haploid and phased diploid sample probabilities were close to conformity. Only 9.18% of the variance of D' could be accounted for by changes in length, indicating that GD is not a monotonically decreasing function of length. We defined two regions of overlapping long-range GD extending 698 735 base pairs (bp) (REN90833/TMLHEInt5 block) and 689 900 bp (F8Int13.2/HEMA154507.3 block) The extent of GD overlap is 575 637 bp (F8Int13.2/TMLHEInt5 interstice). Extended haplotype homozygosity analysis centred at the F8 intronic loci revealed that the most frequent core haplotypes decay the least in the flanking GD. The F8 intronic loci attend distinct non-random association forces; F8Int13.2 serves at maintenance of the long-range overlapping pattern of GD, whereas F8Int25.2 and F8Int22 serve at lessening it in force or effect.

Do common genetic variants in endotoxin signaling pathway contribute to predisposition to alcoholic liver cirrhosis?

BACKGROUND

Tumor necrosis factor-alpha (TNF-alpha) and interleukin-1beta (IL-1beta), produced by endotoxin-activated Kupffer cells, play a key role in the pathogenesis of alcoholic liver cirrhosis (ALC). Alleles TNFA -238A, IL1B -31T and variant IL1RN*2 of repeat polymorphism in the gene encoding the IL-1 receptor antagonist increase production of TNF-alpha and IL-1beta, respectively. Alleles CD14 -159T, TLR4 c.896G and TLR4 c.1196T modify activation of Kupffer cells by endotoxin. We confirmed the published associations between these common variants and genetic predisposition to ALC by means of a large case-control association study conducted on two Central European populations.

METHODS

The study population comprised a Czech sample of 198 ALC patients and 370 controls (MONICA project), and a German sample of 173 ALC patients and 331 controls (KORA-Augsburg), and 109 heavy drinkers without liver disease.

RESULTS

Single locus analysis revealed no significant difference between patients and controls in all tested loci. Diplotype [IL1RN 2/ 2; IL1B -31T+] was associated with increased risk of ALC in the pilot study, but not in the validation samples.

CONCLUSIONS

Although cytokine mediated immune reactions play a role in the pathogenesis of ALC, hereditary susceptibility caused by variants in the corresponding genes is low in Central European populations.

Association of MICA with rheumatoid arthritis independent of known HLA-DRB1 risk alleles in a family-based and a case control study

INTRODUCTION

The gene MICA encodes the protein major histocompatibility complex class I polypeptide-related sequence A. It is expressed in synovium of patients with rheumatoid arthritis (RA) and its implication in autoimmunity is discussed. We analyzed the association of genetic variants of MICA with susceptibility to RA.

METHODS

Initially, 300 French Caucasian individuals belonging to 100 RA trio families were studied. An additional 100 independent RA trio families and a German Caucasian case-control cohort (90/182 individuals) were available for replication. As MICA is situated in proximity to known risk alleles of the HLA-DRB1 locus, our analysis accounted for linkage disequilibrium either by analyzing the subgroup consisting of parents not carrying HLA-DRB1 risk alleles with transmission disequilibrium test (TDT) or by implementing a regression model including all available data. Analysis included a microsatellite polymorphism (GCT)n and single-nucleotide polymorphisms (SNPs) rs3763288 and rs1051794.

RESULTS

In contrast to the other investigated polymorphisms, the non-synonymously coding SNP MICA-250 (rs1051794, Lys196Glu) was strongly associated in the first family cohort (TDT: P = 0.014; regression model: odds ratio [OR] 0.46, 95% confidence interval [CI] 0.25 to 0.82, P = 0.007). Although the replication family sample showed only a trend, combined family data remained consistent with the hypothesis of MICA-250 association independent from shared epitope (SE) alleles (TDT: P = 0.027; regression model: OR 0.56, 95% CI 0.38 to 0.83, P = 0.003). We also replicated the protective association of MICA-250A within a German Caucasian cohort (OR 0.31, 95% CI 0.1 to 0.7, P = 0.005; regression model: OR 0.6, 95% CI 0.37 to 0.96, P = 0.032). We showed complete linkage disequilibrium of MICA-250 (D' = 1, r2= 1) with the functional MICA variant rs1051792 (D' = 1, r2= 1). As rs1051792 confers differential allelic affinity of MICA to the receptor NKG2D, this provides a possible functional explanation for the observed association.

CONCLUSIONS

We present evidence for linkage and association of MICA-250 (rs1051794) with RA independent of known HLA-DRB1 risk alleles, suggesting MICA as an RA susceptibility gene. However, more studies within other populations are necessary to prove the general relevance of this polymorphism for RA.

SNP HiTLink: a high-throughput linkage analysis system employing dense SNP data

BACKGROUND

During this recent decade, microarray-based single nucleotide polymorphism (SNP) data are becoming more widely used as markers for linkage analysis in the identification of loci for disease-associated genes. Although microarray-based SNP analyses have markedly reduced genotyping time and cost compared with microsatellite-based analyses, applying these enormous data to linkage analysis programs is a time-consuming step, thus, necessitating a high-throughput platform.

RESULTS

We have developed SNP HiTLink (SNP High Throughput Linkage analysis system). In this system, SNP chip data of the Affymetrix Mapping 100 k/500 k array set and Genome-Wide Human SNP array 5.0/6.0 can be directly imported and passed to parametric or model-free linkage analysis programs; MLINK, Superlink, Merlin and Allegro. Various marker-selecting functions are implemented to avoid the effect of typing-error data, markers in linkage equilibrium or to select informative data.

CONCLUSION

The results using the 100 k SNP dataset were comparable or even superior to those obtained from analyses using microsatellite markers in terms of LOD scores obtained. General personal computers are sufficient to execute the process, as runtime for whole-genome analysis was less than a few hours. This system can be widely applied to linkage analysis using microarray-based SNP data and with which one can expect high-throughput and reliable linkage analysis.

Heritable variation of ERBB2 and breast cancer risk

Amplification of the epithelial growth factor receptor gene ERBB2 (HER2, NEU) in breast cancer is associated with a poor clinical prognosis. In mammary gland development, this receptor plays a role in ductal and lobuloalveolar differentiation. We conducted a systematic investigation of the role of genetic variation of the ERBB2 gene in breast cancer risk in a study of 842 histologically confirmed invasive breast cancer cases and 1,108 controls from the Shanghai Breast Cancer Study. We observed that the ERBB2 gene resides within a locus of high linkage disequilibrium, composed of three major ancestral haplotypes in the study population. These haplotypes are marked by simple tandem repeat and single nucleotide polymorphisms, including the missense variants I655V and P1170A. We observed a risk-modifying effect of a highly polymorphic simple tandem repeat within an evolutionarily conserved region, 4.4 kb upstream from the ERBB2 transcription start site. Under a dominant genetic model, the age-adjusted odds ratio was 1.74 (95% confidence interval, 1.27-2.37). Its association with breast cancer, and with breast cancer stratified by histology, by histologic grade, and by stage, remained significant after correction for multiple comparisons. In contrast, we observed no association of ERBB2 single nucleotide polymorphism haplotypes with breast cancer predisposition.

Cubic exact solutions for the estimation of pairwise haplotype frequencies: implications for linkage disequilibrium analyses and a web tool 'CubeX'

BACKGROUND

The frequency of a haplotype comprising one allele at each of two loci can be expressed as a cubic equation (the 'Hill equation'), the solution of which gives that frequency. Most haplotype and linkage disequilibrium analysis programs use iteration-based algorithms which substitute an estimate of haplotype frequency into the equation, producing a new estimate which is repeatedly fed back into the equation until the values converge to a maximum likelihood estimate (expectation-maximisation).

RESULTS

We present a program, "CubeX", which calculates the biologically possible exact solution(s) and provides estimated haplotype frequencies, D', r2 and chi2 values for each. CubeX provides a "complete" analysis of haplotype frequencies and linkage disequilibrium for a pair of biallelic markers under situations where sampling variation and genotyping errors distort sample Hardy-Weinberg equilibrium, potentially causing more than one biologically possible solution. We also present an analysis of simulations and real data using the algebraically exact solution, which indicates that under perfect sample Hardy-Weinberg equilibrium there is only one biologically possible solution, but that under other conditions there may be more.

CONCLUSION

Our analyses demonstrate that lower allele frequencies, lower sample numbers, population stratification and a possible |D'| value of 1 are particularly susceptible to distortion of sample Hardy-Weinberg equilibrium, which has significant implications for calculation of linkage disequilibrium in small sample sizes (eg HapMap) and rarer alleles (eg paucimorphisms, q < 0.05) that may have particular disease relevance and require improved approaches for meaningful evaluation.

Haplotype analysis of CYP11A1 identifies promoter variants associated with breast cancer risk

The CYP11A1 gene encodes the cholesterol side chain cleavage enzyme that catalyzes the initial and rate-limiting step of steroidogenesis. A large number of epidemiologic studies have implicated the duration and degree of endogenous estrogen exposure in the development of breast cancer in women. Here, we conduct a systematic investigation of the role of genetic variation of the CYP11A1 gene in breast cancer risk in a study of 1193 breast cancer cases and 1310 matched controls from the Shanghai Breast Cancer Study. We characterize the genetic architecture of the CYP11A1 gene in a Chinese study population. We then genotype tagging polymorphisms to capture common variation at the locus for tests of association. Variants designating a haplotype encompassing the gene promoter are significantly associated with both increased expression (P = 1.6e-6) and increased breast cancer risk: heterozygote age-adjusted odds ratio (OR), 1.51 [95% confidence interval (95% CI), 1.19-1.91]; homozygote age-adjusted OR, 2.94 (95% CI, 1.22-7.12), test for trend, P = 5.0e-5. Among genes controlling endogenous estrogen metabolism, CYP11A1 harbors common variants that may influence expression to significantly modify risk of breast cancer.