Juxtaposed regions of extensive and minimal linkage disequilibrium in human Xq25 and Xq28

The association of gene polymorphisms with milk production and mastitis resistance phenotypic traits in dairy cattle

The aim of this study was to evaluate the association between gene polymorphisms (SNPs) and mastitis indicators and their relationship with milk production profitability in dairy herd.A functional analysis was also performed of five genes containing the studied SNPs and those located close by. DNA was isolated from the hair bulb of 320 dairy cows kept in three herds and SNP-microarray analysis was performed. The data on 299 cows was subjected to final statistical analysis using AI-REML method with one-trait repeatability test-day animal model and pedigree information using the DMU4 package. Five from 35 SNPs significantly associated with mastitis indicators or production traits and located within a gene or no more than 500,000 nucleotides from the gene were selected for the functional and economic analysis. A questionnaire was also developed to collect associated economic data of 219 cows from three herds, such as the value of milk production and direct costs incurred over three years; this allowed the gross margin, direct profitability index and direct costs incurred to produce one liter of milk to be determined, among others. None of the five studied SNPs were related to protein content. The rs110785912(T/A), found near CXCR4, and rs136813430(T/C), located in the TLR4 gene exon, were associated with lnSCC, while rs110455063(C/G), located near IGFI, was associated with milk yield, fat and total solid contents. rs109421300(T/C), associated with fat/protein content ratio, as well as fat and total solid content, is located in the DGAT1 gene intron. rs41587003(A/C), located in the DLG2 gene intron, was associated with lactose content. The economic analysis revealed differences between the variants of the three tested SNPs. The T/C variant of the rs136813430(T/C) SNP was characterized by the highest gross margin, the highest direct profitability index and the lowest costs incurred to produce 1 liter of milk. The T/A variant of rs110785912(T/A) was related to low lnSCC and was characterized by the highest direct profitability index. In turn, the C/C variant of the rs41587003(T/C) was related to the lowest level of lactose and the highest costs of milk production. It appears that rs136813430(T/C) may be the most promising of the tested SNPs for increasing the profitability of milk production. To our knowledge, it is the first effort to assess directly a correlation between the DNA polymorphism and economic output of a dairy enterprise.

Joint Association Analysis Identified 18 New Loci for Bone Mineral Density

Bone mineral density (BMD) at various skeletal sites have shared genetic determinants. In the present study, aiming to identify shared loci associated with BMD, we conducted a joint association study of a genomewide association study (GWAS) and a meta-analysis of BMD at different skeletal sites: (i) a single GWAS of heel BMD in 142,487 individuals from the UK Biobank, and (ii) a meta-analysis of 30 GWASs of total body (TB) BMD in 66,628 individuals from the Genetic Factors for Osteoporosis (GEFOS) Consortium. The genetic correlation coefficient of the two traits was estimated to be 0.57. We performed joint association analysis with a recently developed statistical method multi-trait analysis of GWAS (MTAG) to account for trait heterogeneity and sample overlap. The joint association analysis combining samples of up to 209,115 individuals identified 18 novel loci associated with BMD at the genomewide significance level (α = 5.0 × 10^-8 ), explaining an additional 0.43% and 0.60% of heel-BMD and TB-BMD heritability, respectively. The vast majority of the identified lead SNPs or their proxies exerted local expression quantitative trait loci (cis-eQTL) activity. Credible risk variants, defined as those SNPs located within 500 kilobases (kb) of the lead SNP and with p values within two orders of magnitude of the lead SNP, were enriched in transcription factor binding sites (p = 3.58 × 10^-4 ) and coding regions (p = 5.71 × 10^-4 ). Fifty-six candidate genes were prioritized at these novel loci using multiple sources of information, including several genes being previously reported to play a role in bone biology but not reported in previous GWASs (PPARG, FBN2, DEF6, TNFRSF19, and NFE2L1). One newly identified gene, SCMH1, was shown to upregulate the expression of several bone biomarkers, including alkaline phosphatase (ALP), collagen type 1 (COL-I), osteocalcin (OCN), osteopontin (OPN), and runt-related transcription factor 2 (RUNX2), in mouse osteoblastic MC3T3-E1 cells, highlighting its regulatory role in bone formation. Our results may provide useful candidate genes for future functional investigations. © 2019 American Society for Bone and Mineral Research.

The Genomic Footprints of the Fall and Recovery of the Crested Ibis

Human-induced environmental change and habitat fragmentation pose major threats to biodiversity and require active conservation efforts to mitigate their consequences. Genetic rescue through translocation and the introduction of variation into imperiled populations has been argued as a powerful means to preserve, or even increase, the genetic diversity and evolutionary potential of endangered species [1-4]. However, factors such as outbreeding depression [5, 6] and a reduction in available genetic diversity render the success of such approaches uncertain. An improved evaluation of the consequence of genetic restoration requires knowledge of temporal changes to genetic diversity before and after the advent of management programs. To provide such information, a growing number of studies have included small numbers of genomic loci extracted from historic and even ancient specimens [7, 8]. We extend this approach to its natural conclusion, by characterizing the complete genomic sequences of modern and historic population samples of the crested ibis (Nipponia nippon), an endangered bird that is perhaps the most successful example of how conservation effort has brought a species back from the brink of extinction. Though its once tiny population has today recovered to >2,000 individuals [9], this process was accompanied by almost half of ancestral loss of genetic variation and high deleterious mutation load. We furthermore show how genetic drift coupled to inbreeding following the population bottleneck has largely purged the ancient polymorphisms from the current population. In conclusion, we demonstrate the unique promise of exploiting genomic information held within museum samples for conservation and ecological research.

19p13.3-GADD45B common variants and 19q13.3-PPP1R13L and 19q13.3-CD3EAP in lung cancer risk among Chinese

Lung cancer is the most common cause of cancer-related mortality worldwide. The GADD45 gene family plays important roles in a variety of the responses to cell injury including cell cycle checkpoints, apoptosis, DNA repair and anti-tumor immunity. The 19p13.3-GADD45B encoded protein product is involved in apoptosis and inhibiting tumor growth. To evaluate the association of 19p13.3-GADD45B common variants and lung cancer risk, the present study containing 544 Chinese lung cancer cases and 550 cancer-free controls was conducted. Three htSNPs (haplotype-tagging single nucleotide polymorphism) (rs7354, rs14384, and rs3783501) covering 95% of the common haplotype diversity in 19p13.3-GADD45B and interaction of 19p13.3-GADD45B and 19q13.3-PPP1R13L and 19q13.3-CD3EAP variants and smoking-duration were explored. Genotype and allele frequencies and haplotype distributions of the 19p13.3-GADD45B 3 htSNPs were not associated with lung cancer risk after adjustment for smoking status. 19p13.3-GADD45B rs7354 was associated with lung cancer risk among ≤20 (years) smokers [C/A-A/A versus CC, OR (95% CI) = 3.20 (1.11-9.20), P = 0.025] in a dominant model stratified by smoking duration. MDR (multifactor dimensionality reduction) analyses showed that smoking history as main effect and three-way models (smoking duration, 19p13.3-GADD45B rs3783501, 19q13.3-CD3EAP rs967591) (P = 0.001-0.002) indicated statistically significant association with lung cancer risk. The study identified evidence implicating DNA damage response genes on chromosome 19 in etiology of smoke-exposed lung cancer. In conclusion, our findings demonstrate that 19p13.3-GADD45B rs7354 variant and interaction between 19p13.3-GADD45B rs3783501 and 19q13.3-CD3EAP rs967591 may play a role in association with smoke-exposed lung cancer among Chinese. 19p13.3-GADD45B variants should be further evaluated in large prospective studies with molecular pathological annotations of lung cancer.

Combined effects of multiple linked loci on pairwise sibling tests

The advanced multiplex STR system, PowerPlex Fusion, includes four linked locus pairs. The conventional Identifiler system has one pair of linked loci. Therefore, sibling tests conducted using the advanced system might be more affected by linkage than those conducted using the conventional system. This study simulated single and combined effects of the four linked locus pairs on pairwise sibling tests. Simulated genotypes of 100,000 pairs of full siblings and nonrelatives were constructed according to allele frequencies of the Japanese population. The single linkage effect was evaluated for simulated genotype data by calculating both the likelihood ratio accounting for the linkage between two loci and the likelihood ratio ignoring the linkage. The combined effect was obtained by multiplication of the respective single effects. Furthermore, we investigated the possibility that ignoring the linkage affects subject classification by introducing a scale of the likelihood ratio into sibling tests. The single effect in the Identifiler analysis was 0.645-1.746 times if the linkage was ignored. Overestimations and underestimations were predictable from the identical-by-state status at two linked loci. The combined effect in the PowerPlex Fusion analysis was 0.217-7.390 times. Ignoring the linkage rarely caused a false conclusive or inconclusive result, even from PowerPlex Fusion analysis. Application of the advanced system improved sibling tests considerably. The additional examined loci were more beneficial than the adverse effect of the linkage derived from the four linked locus pairs.

Genome-Wide Association Study Identifies Candidate Loci Underlying Agronomic Traits in a Middle American Diversity Panel of Common Bean

Common bean ( L.) breeding programs aim to improve both agronomic and seed characteristics traits. However, the genetic architecture of the many traits that affect common bean production are not completely understood. Genome-wide association studies (GWAS) provide an experimental approach to identify genomic regions where important candidate genes are located. A panel of 280 modern bean genotypes from race Mesoamerica, referred to as the Middle American Diversity Panel (MDP), were grown in four US locations, and a GWAS using >150,000 single-nucleotide polymorphisms (SNPs) (minor allele frequency [MAF] ≥ 5%) was conducted for six agronomic traits. The degree of inter- and intrachromosomal linkage disequilibrium (LD) was estimated after accounting for population structure and relatedness. The LD varied between chromosomes for the entire MDP and among race Mesoamerica and Durango-Jalisco genotypes within the panel. The LD patterns reflected the breeding history of common bean. Genome-wide association studies led to the discovery of new and known genomic regions affecting the agronomic traits at the entire population, race, and location levels. We observed strong colocalized signals in a narrow genomic interval for three interrelated traits: growth habit, lodging, and canopy height. Overall, this study detected ∼30 candidate genes based on a priori and candidate gene search strategies centered on the 100-kb region surrounding a significant SNP. These results provide a framework from which further research can begin to understand the actual genes controlling important agronomic production traits in common bean.

Genetic structure of the Newfoundland and Labrador population: founder effects modulate variability

The population of the province of Newfoundland and Labrador (NL) has been a resource for genetic studies because of its historical isolation and increased prevalence of several monogenic disorders. Controversy remains regarding the genetic substructure and the extent of genetic homogeneity, which have implications for disease gene mapping. Population substructure has been reported from other isolated populations such as Iceland, Finland and Sardinia. We undertook this study to further our understanding of the genetic architecture of the NL population. We enrolled 494 individuals randomly selected from NL. Genome-wide SNP data were analyzed together with that from 14 other populations including HapMap3, Ireland, Britain and Native American samples from the Human Genome Diversity Project. Using multidimensional scaling and admixture analysis, we observed that the genetic structure of the NL population resembles that of the British population but can be divided into three clusters that correspond to religious/ethnic origins: Protestant English, Roman Catholic Irish and North American aboriginals. We observed reduced heterozygosity and an increased inbreeding coefficient (mean=0.005), which corresponds to that expected in the offspring of third-cousin marriages. We also found that the NL population has a significantly higher number of runs of homozygosity (ROH) and longer lengths of ROH segments. These results are consistent with our understanding of the population history and indicate that the NL population may be ideal for identifying recessive variants for complex diseases that affect populations of European origin.

Vascular endothelial growth factor (VEGF) gene polymorphisms and breast cancer risk in Punjabi population from North West India

The purpose of this study was to evaluate the association of seven VEGF promoter polymorphisms with breast cancer risk in Punjabi population from North West India. We screened DNA samples of 102 sporadic breast cancer patients and 102 unrelated healthy, gender, and age-matched individuals for seven VEGF promoter polymorphisms [-417 C/T (rs833062), -172 C/A (rs59260042), -165 C/T (rs79469752), -160 C/T, -152 G/A (rs13207351), -141 A/C (rs28357093) and -116 G/A (rs1570360)] by direct sequencing. The frequency of GG, GA, and AA genotype of -152 G/A polymorphism was 26.47 vs 38.34%, 46.08 vs 51.96%, and 27.45 vs 9.80%, in patients and controls, respectively. VEGF -152 AA genotype was significantly associated with increased risk for breast cancer (OR = 4.04, 95%CI, 1.69-9.68, p = 0.001; recessive model OR = 3.48, 95%CI, 1.59-7.63, p = 0.001). For VEGF -116 G/A polymorphism, G and A allele frequencies were 65.2 vs 76.47% and 34.8 vs 23.53% in patients and controls, respectively. Individuals having -116 AA genotype (OR = 3.40; 95%CI, 1.24-9.37; p = 0.014) and A allele (OR = 1.73; 95%CI, 1.12-2.67; p = 0.012) were associated with increased risk for breast cancer. VEGF -165 C/T and -141 A/C polymorphisms were associated with reduced risk for breast cancer. There was significantly decreased frequency of CT genotype (4.90 vs 18.63%; p = 0.002) and T allele (2.45 vs 9.31%; p = 0.003) of -165 C/T polymorphism among breast cancer patients as compared to controls. VEGF -141 A and C allele frequency were 96.57 vs 91.18% and 3.43 vs 8.82% in patients and controls, respectively. Significant reduced risk for breast cancer was observed with AC genotype (OR = 0.34, 95%CI, 0.14-0.86; p = 0.019) and C allele (OR = 0.37; 95%CI, 0.15-0.89; p = 0.023) of -141 A/C polymorphism. We did not observe association of VEGF -417 T/C, -172 C/A, -160 C/T polymorphisms with breast cancer risk in the studied subjects (p > 0.05). The VEGF -152 G/A and -116 G/A polymorphisms were found to be significantly associated with increased risk for breast cancer while -165 C/T and -141 A/C polymorphisms were found to be associated with decreased risk for breast cancer in Punjabi population from North West India.

Association of single-nucleotide polymorphisms in a metabotropic glutamate receptor GRM3 gene subunit to alcohol-dependent male subjects

AIMS

The purpose of this study was to investigate the association between the metabotropic glutamate receptor 3 (GRM3) subunit gene and alcohol dependence by the single-nucleotide polymorphisms (SNPs).

METHODS

Two hundred and forty-eight male alcohol-dependent patients and 235 male control subjects were recruited. Ten SNPs in the GRM3 region were studied, and genotyping of SNPs was performed by ligase detection reactions.

RESULTS

We found highly significant differences in allele and genotype frequencies of rs6465084 between the alcohol-dependent and control group, with the greater frequency of A allele of SNP rs6465084 in alcohol-dependent group. We also found significant differences of haplotype frequencies in five combinations (including TAATATT, CAGTATT, TCGTATT, CAATAGC, TAATATC) in the linkage disequilibrium constructed by seven SNPs between the groups.

CONCLUSION

Our results supplied the first evidence that the polymorphism of GRM3 gene associates with the morbidity of alcohol dependence in human being, which may support a new potential target for alcoholism treatment.

Genetic mapping of high caries experience on human chromosome 13

Background

Our previous genome-wide linkage scan mapped five loci for caries experience. The purpose of this study was to fine map one of these loci, the locus 13q31.1, in order to identify genetic contributors to caries.

Methods

Seventy-two pedigrees from the Philippines were studied. Caries experience was recorded and DNA was extracted from blood samples obtained from all subjects. Sixty-one single nucleotide polymorphisms (SNPs) in 13q31.1 were genotyped. Association between caries experience and alleles was tested. We also studied 1,481 DNA samples obtained from saliva of subjects from the USA, 918 children from Brazil, and 275 children from Turkey, in order to follow up the results found in the Filipino families. We used the AliBaba2.1 software to determine if the nucleotide changes of the associated SNPs changed the prediction of the presence of transcription-binding site sequences and we also analyzed the gene expression of the genes selected based on binding predictions. Mutation analysis was also performed in 33 Filipino individuals of a segment of 13q31.1 that is highly conserved in mammals.

Results

Statistically significant association with high caries experience was found for 11 markers in 13q31.1 in the Filipino families. Haplotype analysis also confirmed these results. In the populations used for follow-up purposes, associations were found between high caries experience and a subset of these markers. Regarding the prediction of the transcription-binding site, the base change of the SNP rs17074565 was found to change the predicted-binding of genes that could be involved in the pathogenesis of caries. When the sequence has the allele C of rs17074565, the potential transcription factors binding the sequence are GR and GATA1. When the subject carries the G allele of rs17074565, the potential transcription factor predicted to bind to the sequence is GATA3. The expression of GR in whole saliva was higher in individuals with low caries experience when compared to individuals with high caries experience (p = 0.046). No mutations were found in the highly conserved sequence.

Conclusions

Genetic factors contributing to caries experience may exist in 13q31.1. The rs17074565 is located in an intergenic region and is predicted to disrupt the binding sites of two different transcription factors that might be involved with caries experience. GR expression in saliva may be a biomarker for caries risk and should be further explored.

Two- and three-locus haplotypes of the paraoxonase (PON1) gene are associated with coronary artery disease in Asian Indians

INTRODUCTION

In view of the reported association of SNPs in the paraoxonase (PON1) gene with coronary artery disease (CAD), and the absence of conclusive data from India, we investigated the relationship of three SNPs at different loci (-108C/T, L55M and Q192R) of the PON1 gene and their haplotypes with CAD among people residing in the northern plains of India.

MATERIALS AND METHODS

One hundred and seventy-eight healthy controls and two hundred and four angiographically-proven CAD patients were genotyped using PCR-RFLP.

RESULTS

Of the three SNPs, only the R allele of Q192R polymorphism was associated with CAD (p<0.05). Two locus haplotypes QT (OR 0.55, p=0.0004, 95% CI 0.39-0.77, significant) and LQ (odds ratio 0.73, p=0.03, 95% CI 0.55-0.97, trend) showed protective effects, while haplotypes MR (OR=5.36, p=0.0001, 95% CI 2.045-14.049) and MC (OR=2.71, p=0.011, 95% CI 1.221-6.046) were associated with increased risk of CAD. MRT, a minor three-locus haplotype also displayed significant association (OR 4.93, 95% CI 1.7-13.5) with the disease. Significance was assessed after applying Bonferroni's correction.

CONCLUSIONS

Our study revealed that only one SNP at a single locus but several haplotype combinations of PON1 coding and promoter-region polymorphisms were associated with the risk of or protection against CAD. Thus, haplotype analysis brought better insights into the association of PON1 gene polymorphisms with CAD in Asian Indians.

Genome-Wide Linkage Disequilibrium from 100,000 SNPs in the East Finland Founder Population

Information about linkage disequilibrium (LD) is important in understanding the genome structure and has its applications in association studies. Here we present the first genome-wide LD study based on a founder population (East Finland). The LD data consist of 118 unrelated individuals and around 480,000 SNP pairs genotyped with the Affymetrix 100K genotyping assay. Using the minor allele frequency (MAF) limit of .05, the squared correlation coefficient between two loci (r2) was .48, .37, .28, and .20 for distances of 5, 10, 20, and 40 kb respectively. MAF had a significant effect on the mean r2 so that the extent of useful LD (r2 > .3) varied from 17 kb to 80 kb depending on the limit set for the MAF. For D' the effect of MAF was smaller but reflected the possible age of the mutation: SNPs with high MAF had lower D' than those with low MAF. The X chromosome showed higher D' values than autosomes and the extent of useful LD (r2 > .3) was twice as long on the X chromosome than on the autosomes. Based on the results, LD varies across the genome and is correlated to local recombination rate between and within chromosomes. However, the recombination rate does not explain all the variation found in LD. We also report a number of long chromosomal regions where exceptionally high or low LD were detected.

APOE and AβPP gene variation in cortical and cerebrovascular amyloid-β pathology and Alzheimer's disease: a population-based analysis

Cortical and cerebrovascular amyloid-β (Aβ) deposition is a hallmark of Alzheimer's disease (AD), but also occurs in elderly people not affected by dementia. The apolipoprotein E (APOE) ε4 is a major genetic modulator of Aβ deposition and AD risk. Variants of the amyloid-β protein precursor (AβPP) gene have been reported to contribute to AD and cerebral amyloid angiopathy (CAA). We analyzed the role of APOE and AβPP variants in cortical and cerebrovascular Aβ deposition, and neuropathologically verified AD (based on modified NIA-RI criteria) in a population-based autopsy sample of Finns aged ≥ 85 years (Vantaa85 + Study; n = 282). Our updated analysis of APOE showed strong associations of the ε4 allele with cortical (p = 4.91 × 10-17) and cerebrovascular (p = 9.87 × 10-11) Aβ deposition as well as with NIA-RI AD (p = 1.62 × 10-8). We also analyzed 60 single nucleotide polymorphisms (SNPs) at the AβPP locus. In single SNP or haplotype analyses there were no statistically significant AβPP locus associations with cortical or cerebrovascular Aβ deposition or with NIA-RI AD. We sequenced the promoter of the AβPP gene in 40 subjects with very high Aβ deposition, but none of these subjects had any of the previously reported or novel AD-associated mutations. These results suggest that cortical and cerebrovascular Aβ depositions are useful quantitative traits for genetic studies, as highlighted by the strong associations with the APOE ε4 variant. Promoter mutations or common allelic variation in the AβPP gene do not have a major contribution to cortical or cerebrovascular Aβ deposition, or very late-onset AD in this Finnish population based study.

Population-genetic comparison of the Sorbian isolate population in Germany with the German KORA population using genome-wide SNP arrays

Background

The Sorbs are an ethnic minority in Germany with putative genetic isolation, making the population interesting for disease mapping. A sample of N = 977 Sorbs is currently analysed in several genome-wide meta-analyses. Since genetic differences between populations are a major confounding factor in genetic meta-analyses, we compare the Sorbs with the German outbred population of the KORA F3 study (N = 1644) and other publically available European HapMap populations by population genetic means. We also aim to separate effects of over-sampling of families in the Sorbs sample from effects of genetic isolation and compare the power of genetic association studies between the samples.

Results

The degree of relatedness was significantly higher in the Sorbs. Principal components analysis revealed a west to east clustering of KORA individuals born in Germany, KORA individuals born in Poland or Czech Republic, Half-Sorbs (less than four Sorbian grandparents) and Full-Sorbs. The Sorbs cluster is nearest to the cluster of KORA individuals born in Poland. The number of rare SNPs is significantly higher in the Sorbs sample. FST between KORA and Sorbs is an order of magnitude higher than between different regions in Germany. Compared to the other populations, Sorbs show a higher proportion of individuals with runs of homozygosity between 2.5 Mb and 5 Mb. Linkage disequilibrium (LD) at longer range is also slightly increased but this has no effect on the power of association studies.

Oversampling of families in the Sorbs sample causes detectable bias regarding higher FST values and higher LD but the effect is an order of magnitude smaller than the observed differences between KORA and Sorbs. Relatedness in the Sorbs also influenced the power of uncorrected association analyses.

Conclusions

Sorbs show signs of genetic isolation which cannot be explained by over-sampling of relatives, but the effects are moderate in size. The Slavonic origin of the Sorbs is still genetically detectable.

Regarding LD structure, a clear advantage for genome-wide association studies cannot be deduced. The significant amount of cryptic relatedness in the Sorbs sample results in inflated variances of Beta-estimators which should be considered in genetic association analyses.

Mendelian breeding units versus standard sampling strategies: Mitochondrial DNA variation in southwest Sardinia

We report a sampling strategy based on Mendelian Breeding Units (MBUs), representing an interbreeding group of individuals sharing a common gene pool. The identification of MBUs is crucial for case-control experimental design in association studies. The aim of this work was to evaluate the possible existence of bias in terms of genetic variability and haplogroup frequencies in the MBU sample, due to severe sample selection. In order to reach this goal, the MBU sampling strategy was compared to a standard selection of individuals according to their surname and place of birth. We analysed mitochondrial DNA variation (first hypervariable segment and coding region) in unrelated healthy subjects from two different areas of Sardinia: the area around the town of Cabras and the western Campidano area. No statistically significant differences were observed when the two sampling methods were compared, indicating that the stringent sample selection needed to establish a MBU does not alter original genetic variability and haplogroup distribution. Therefore, the MBU sampling strategy can be considered a useful tool in association studies of complex traits.

Genomic Diversity in Pig (Sus scrofa) and its Comparison with Human and other Livestock

We have reviewed the current pig (Sus scrofa) genomic diversity within and between sites and compared them with human and other livestock. The current Porcine 60K single nucleotide polymorphism (SNP) panel has an average SNP distance in a range of 30 - 40 kb. Most of genetic variation was distributed within populations, and only a small proportion of them existed between populations. The average heterozygosity was lower in pig than in human and other livestock. Genetic inbreeding coefficient (F(IS)), population differentiation (F(ST)), and Nei's genetic distance between populations were much larger in pig than in human and other livestock. Higher average genetic distance existed between European and Asian populations than between European or between Asian populations. Asian breeds harboured much larger variability and higher average heterozygosity than European breeds. The samples of wild boar that have been analyzed displayed more extensive genetic variation than domestic breeds. The average linkage disequilibrium (LD) in improved pig breeds extended to 1 - 3 cM, much larger than that in human (~ 30 kb) and cattle (~ 100 kb), but smaller than that in sheep (~ 10 cM). European breeds showed greater LD that decayed more slowly than Asian breeds. We briefly discuss some processes for maintaining genomic diversity in pig, including migration, introgression, selection, and drift. We conclude that, due to the long time of domestication, the pig possesses lower heterozygosity, higher F(IS), and larger LD compared with human and cattle. This implies that a smaller effective population size and less informative markers are needed in pig for genome wide association studies.

Evaluating differences in linkage disequilibrium between populations

We propose two methods to evaluate the statistical significance of differences in linkage disequilibrium (LD) between populations, where LD is measured by the standardised parameter D'. The first method is based on bootstrapping individuals within populations in order to test LD differences for each pair of loci. Using this approach we propose a solution to the problem of testing multiple locus-pairs by means of a single test for the number of pairs that exhibit significant LD differences among populations. The second method provides the Bayesian posterior probability that one population has greater LD than the other for each locus pair. Both methods can handle genotypes with unknown phase, and are demonstrated using two data sets. For the purpose of demonstration, we apply the methods to two different sets of data from humans. First, we explore the issue of LD differences between reproductively isolated populations using a new data set of twelve Xq25 microsatellites, typed in four European populations. Second, we examine evidence for LD differences between Alzheimer cases and controls from the Icelandic population using 19 single nucleotide polymorphisms (SNPs) from a 97 kb region flanking the Apolipoprotein E (APOE) gene on chromosome 19.

Shades of gray: a comparison of linkage disequilibrium between Hutterites and Europeans

Founder or isolated populations have advantages for genetic studies due to decreased genetic and environmental heterogeneity. However, whereas longer-range linkage disequilibrium (LD) in these populations is expected to facilitate gene localization, extensive LD may actually limit the ability for gene discovery. The North American Hutterite population is one of the best characterized young founder populations and members of this isolate have been the subjects of our studies of complex traits, including fertility, asthma and cardiovascular disease, for >20 years. Here, we directly assess the patterns and extent of global LD using single nucleotide polymorphism genotypes with minor allele frequencies (MAFs) > or =5% from the Affymetrix GeneChip Mapping 500 K array in 60 relatively unrelated Hutterites and 60 unrelated Europeans (HapMap CEU). Although LD among some marker pairs extends further in the Hutterites than in Europeans, the pattern of LD and MAF are surprisingly similar. These results indicate that (1) identifying disease genes should be no more difficult in the Hutterites than in outbred European populations, (2) the same common susceptibility alleles for complex diseases should be present in the Hutterites and outbred European populations, and (3) imputation algorithms based on HapMap CEU should be applicable to the Hutterites.

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.

Linkage disequilibrium patterns and genetic structure of Amerindian and non-Amerindian Brazilian populations revealed by long-range X-STR markers

The extent of X-chromosome linkage disequilibrium (LD) was studied in a southern Brazilian population, and in a pool of samples from Amerindian populations. For this purpose, 11 microsatellites, located mostly in a Xq region comprising approximately 86 Mb was investigated. The lower Amerindian gene diversity associated with significant differences between the populations studied indicated population structure as the main cause for the higher LD values in the Amerindian pool. On the other hand, the LD levels of the non-Amerindian Brazilian sample, although less extensive than that of the Amerindians, were probably determined by admixture events. Our results indicated that different demographic histories have significant effects on LD levels of human populations, and provide a first approach to the X-chromosome ancestry of Amerindian and non-Amerindian Brazilian populations, being valuable for future studies involving mapping and population genetic studies.

Common polymorphic variation in the genetically diverse African insulin gene and its association with size at birth

The insulin variable number of tandem repeats (INS VNTR) has been variably associated with size at birth in non-African populations. Small size at birth is a major determinant of neonatal mortality, so the INS VNTR may influence survival. We tested the hypothesis, therefore, that genetic variation around the INS VNTR in a rural Gambian population, who experience seasonal variation in nutrition and subsequently birth weight, may be associated with foetal and early growth. Six polymorphisms flanking the INS VNTR were genotyped in over 2,500 people. Significant associations were detected between the maternally inherited SNP 27 (rs689) allele and birth length [effect size 17.5 (5.2-29.8) mm; P = 0.004; n = 361]. Significant associations were also found between the maternally inherited African-specific SNP 28 (rs5506) allele and post-natal weight gain [effect size 0.19 (0.05-0.32) z score points/year; P = 0.005; n = 728). These results suggest that in the Gambian population studied there are associations between polymorphic variation in the genetically diverse INS gene and foetal and early growth characteristics, which contribute to overall polygenic associations with these traits.

Association of the growth hormone receptor gene polymorphisms with mandibular height in a Korean population

Growth hormone receptor gene (GHR) is one of the likely candidates for determining morphological traits, because GH is a key regulator of bone growth. The genetic association of GHR in exon 10 with mandibular ramus height has been found in different populations, Japanese and Chinese. On the other hand, two common isoforms of GHR, one full-length (fl-GHR) and the other lacking the extracellular domain encoded by exon 3 (d3-GHR), are associated with differences in responsiveness to GH. The purpose of this study involving 159 Korean subjects was to study the associations between a GHR polymorphism (d3/fl-GHR) that results in genomic deletion of exon 3 and craniofacial morphology, and to study the associations between GHR genotypes in exon 10 and craniofacial morphology. Moreover, the allelic frequencies in a multi-ethnic population (24 Han Chinese, 24 African-Americans, 24 European-Americans, and 24 Hispanics) in a GHR polymorphism (d3/fl-GHR) were compared in this study. The five craniofacial linear measurements (cranial base length, maxillary length, overall mandibular length, mandibular corpus length, and mandibular ramus height) obtained from lateral cephalograms were examined as craniofacial morphology. We found that the d3/fl-GHR polymorphism had no association for any measurements, and a statistically significant association (P=0.024) between the GHR polymorphisms P561T and C422F in exon 10 and mandibular ramus height. Neither SNPs besides P561T and C422F polymorphisms in exon 10 nor the measurements besides mandibular ramus height have statistically significances. Both derived alleles at P561T and C422F SNPs were highly associated with only one haplotype, haplotype-4 in Korean population. As quantitative haplotype association, the results showed a significant difference in mandibular ramus height between individuals having one haplotype-4 and others without haplotype-4 (P=0.028). Moreover, we found that the d3/fl-GHR polymorphism showed diverse frequency in different population. Regarding GHR genotypes in exon 10, the present study mostly reflected the results obtained for a Japanese population, although our current study does not replicate the correlation between the I526L polymorphism of GHR and mandibular ramus height as was reported in a previous study of Han Chinese. The results of the present study suggest that the GHR exon 10 SNPs, not d3/fl-GHR, contribute to changes in the mandibular ramus height of Koreans.

A follow-up study of chromosome 19q13 in multiple sclerosis susceptibility

A possible role of allelic variation on chromosome 19q13 in multiple sclerosis (MS) susceptibility has been suggested. We tested association of sixteen 19q13 markers with MS in 459 families. Nominally significant associations were tested in an independent set of 323 families as well as in the pooled set of 782 families. We were not able to confirm previously suggested associations with APOE, GIPR, ZNF45, ILT6 and D19S585. In the screening dataset nominally significant associations were found with D19S867 and with APOE haplotype (p=0.007 in both), but these were not replicated in the independent dataset nor in the pooled analysis of 757 families. Thus, we were not able to detect any statistically significant allelic associations. Re-sequencing based approaches may be required for elucidating the role chromosome 19q13 with MS.

Haplotype block partitioning as a tool for dimensionality reduction in SNP association studies

BACKGROUND

Identification of disease-related genes in association studies is challenged by the large number of SNPs typed. To address the dilution of power caused by high dimensionality, and to generate results that are biologically interpretable, it is critical to take into consideration spatial correlation of SNPs along the genome. With the goal of identifying true genetic associations, partitioning the genome according to spatial correlation can be a powerful and meaningful way to address this dimensionality problem.

RESULTS

We developed and validated an MCMC Algorithm To Identify blocks of Linkage DisEquilibrium (MATILDE) for clustering contiguous SNPs, and a statistical testing framework to detect association using partitions as units of analysis. We compared its ability to detect true SNP associations to that of the most commonly used algorithm for block partitioning, as implemented in the Haploview and HapBlock software. Simulations were based on artificially assigning phenotypes to individuals with SNPs corresponding to region 14q11 of the HapMap database. When block partitioning is performed using MATILDE, the ability to correctly identify a disease SNP is higher, especially for small effects, than it is with the alternatives considered. Advantages can be both in terms of true positive findings and limiting the number of false discoveries. Finer partitions provided by LD-based methods or by marker-by-marker analysis are efficient only for detecting big effects, or in presence of large sample sizes. The probabilistic approach we propose offers several additional advantages, including: a) adapting the estimation of blocks to the population, technology, and sample size of the study; b) probabilistic assessment of uncertainty about block boundaries and about whether any two SNPs are in the same block; c) user selection of the probability threshold for assigning SNPs to the same block.

CONCLUSION

We demonstrate that, in realistic scenarios, our adaptive, study-specific block partitioning approach is as or more efficient than currently available LD-based approaches in guiding the search for disease loci.

Nonrecurrent MECP2 duplications mediated by genomic architecture-driven DNA breaks and break-induced replication repair

Recurrent submicroscopic genomic copy number changes are the result of nonallelic homologous recombination (NAHR). Nonrecurrent aberrations, however, can result from different nonexclusive recombination-repair mechanisms. We previously described small microduplications at Xq28 containing MECP2 in four male patients with a severe neurological phenotype. Here, we report on the fine-mapping and breakpoint analysis of 16 unique microduplications. The size of the overlapping copy number changes varies between 0.3 and 2.3 Mb, and FISH analysis on three patients demonstrated a tandem orientation. Although eight of the 32 breakpoint regions coincide with low-copy repeats, none of the duplications are the result of NAHR. Bioinformatics analysis of the breakpoint regions demonstrated a 2.5-fold higher frequency of Alu interspersed repeats as compared with control regions, as well as a very high GC content (53%). Unexpectedly, we obtained the junction in only one patient by long-range PCR, which revealed nonhomologous end joining as the mechanism. Breakpoint analysis in two other patients by inverse PCR and subsequent array comparative genomic hybridization analysis demonstrated the presence of a second duplicated region more telomeric at Xq28, of which one copy was inserted in between the duplicated MECP2 regions. These data suggest a two-step mechanism in which part of Xq28 is first inserted near the MECP2 locus, followed by breakage-induced replication with strand invasion of the normal sister chromatid. Our results indicate that the mechanism by which copy number changes occur in regions with a complex genomic architecture can yield complex rearrangements.

Genetic association of low-density lipoprotein receptor-related protein 2 (LRP2) with plasma lipid levels

AIM

Not all genetic factors predisposing phenotypic features of dyslipidemia have been identified. We studied the association between the low density lipoprotein-related protein 2 gene (LRP2) and levels of plasma total cholesterol (T-Cho) and LDL-cholesterol (LDL-C) among 352 adults in Japan.

METHODS

Subjects were obtained from among participants in a cohort study that was carried out with health-check screening in an area of east-central Japan. We selected 352 individuals whose LDL-C levels were higher than 140 mg/dL from the initially screened 22,228 people. We assessed the relation between plasma cholesterol levels and single-nucleotide polymorphisms (SNPs) in the LRP2 gene.

RESULTS

We identified significant correlations between plasma cholesterol levels and two of 19 examined SNPs in LRP2, c.+193826T/C and IVS55 - 147A/G. In particular, the association of c.+193826T/C with the T-Cho level was prominent (p=0.003), showing a co-dominant effect of the minor C-allele on lowering T-Cho and LDL-C levels: for 24 homozygous C-allele carriers, T-Cho=240.7 +/- 24.2 mg/dL and LDL-C=166.1 +/- 21.0 mg/dL); for 130 heterozygous carriers, 248.5 +/- 23.5 mg/dL and 166.6 +/- 19.3 mg/dL; and for 196 homozygous T-allele carriers, 253.9 +/- 23.5 mg/dL and 172.0 +/- 21.0 mg/dL. Linkage disequilibrium (LD) analyses based on 19 selected SNPs showed that c.+193826T/C and IVS55 - 147A/G were in tight LD and that both were located in an LD block covering the genomic sequence from exon 55 to exon 61.

CONCLUSION

We confirm the association between LRP2 and levels of T-Cho and LDL-C in human plasma. The results suggest that genetic variations in LRP2 are important factors affecting lipoprotein phenotypes of patients with hypercholesterolemia.

Population genetics of anthelmintic resistance in parasitic nematodes

A key aim of anthelmintic resistance research is to identify molecular markers that could form the basis of sensitive and accurate diagnostic tests. These would provide powerful tools to study the origin and spread of anthelmintic resistance in the field and to monitor strategies aimed at preventing and managing resistance. Molecular markers could also form the basis of routine diagnostic tests for use in surveillance and clinical veterinary practice. Much of the research conducted to date has focused on the investigation of possible associations of particular candidate genes with the resistance phenotype. In the future, as full parasite genome sequences become available, there will be an opportunity to apply genome-wide approaches to identify the genetic loci that underlie anthelmintic resistance. Both the interpretation of candidate gene studies and the application of genome-wide approaches require a good understanding of the genetics and population biology of the relevant parasites as well as knowledge of how resistance mutations arise and are selected in populations. Unfortunately, much of this information is lacking for parasitic nematodes. This review deals with a number of aspects of genetics and population biology that are pertinent to these issues. We discuss the possible origins of resistance mutations and the likely effects of subsequent selection on the genetic variation at the resistance-conferring locus. We also review some of the experimental approaches that have been used to test associations between candidate genes and anthelmintic resistance phenotypes and highlight implications for future genome-wide studies.

Association of a single-nucleotide variation (A1330V) in the low-density lipoprotein receptor-related protein 5 gene (LRP5) with bone mineral density in adult Japanese women

Low-density lipoprotein receptor-related protein 5 (LRP5), a co-receptor of Wnt signaling, is an important regulator of bone development and maintenance. Recently we identified correlation between an intronic single-nucleotide polymorphism (SNP) in the LRP5 gene and vertebral bone mineral density (BMD), indicating that a genetic ground exists at this locus for determination of BMD. In the study reported here, we searched for nucleotide variation(s) that might confer susceptibility to osteoporosis among an extended panel of 387 healthy subjects recruited from the same hospital (Group-A), as well as among 384 subjects from the general population in eastern Japan (Group-B). We basically focused on two potentially functional variations, Q89R (c.266A > G) and A1330V (c.3989C > T), whose functional effects by the amino-acid changes were estimated by the SIFT software program; it predicted the 1330 V allele as deleterious ("intolerant") although the minor allele of Q89R was questionable. By analyzing associations between the variant alleles and the BMD, reproducible association of the minor variant of A1330V to lower adjusted BMD levels was detected; i.e., In Group-A subjects 1330-V significantly associated with the spinal BMD Z-score (P = 0.034), and in Group-B it associated with low radial BMD (P = 0.019). From haplotype and linkage disequilibrium (LD) analysis for 29 SNPs, we detected two separate LD blocks within the entire 137-kb LRP5 locus, basically consistent with a previous report on Caucasians. One of the second block haplotype significantly associated with adjusted BMD (r = 0.15, P = 0.004). Possible combined effect of Q89R and A1330V belonging to different LD blocks was denied by multiple regression analyses. Our results indicate that genetic variations in LRP5 are important factors affecting BMD in adult women and that 1330 V may contribute to osteoporosis susceptibility, at least in Japanese.

Haplotype association of IL-8 gene with Behcet's disease

Interleukin-8 (IL-8), a CXC chemokine that recruits and activates inflammatory cells, plays a critical role in the pathogenesis of Behcet's disease (BD). To investigate the association of the genetic polymorphism of IL-8 and BD, we genotyped IL-8 -845 T/C, -738 T/A, -353 A/T, -251 A/T, +293 G/T, +678 T/C and receptors CXCR-1 +2607 G/C and CXCR-2 +785 C/T polymorphisms in 119 Korean patients with BD and 119 age- and sex-matched healthy blood donors. Then, single nucleotide polymorphisms (SNPs) and haplotypes were analyzed between patients and controls. There were no SNPs associated with BD. However, the frequency of haplotype TAT inferred from SNPs, IL-8 -353 A/T, -251 A/T and +678 T/C, was significantly higher in patients with BD than controls (5.9 vs 0.0%, P = 0.0001), as was haplotype ATC (6.7 vs 0.0%, P < 0.0001). The haplotype difference was still valid in human leukocyte antigen-B51-negative subjects. In conclusion, we found a significant difference in the distribution of IL-8 gene haplotypes between patients with BD and healthy controls. These results suggest that the genetic polymorphisms of proinflammatory chemokine IL-8 can contribute to the pathogenesis of BD.

Power-based, phase-informed selection of single nucleotide polymorphisms for disease association screens

Single nucleotide polymorphisms (SNPs) are becoming widely used as genotypic markers in genetic association studies of common, complex human diseases. For such association screens, a crucial part of study design is determining what SNPs to prioritize for genotyping. We present a novel power-based algorithm to select a subset of tag SNPs for genotyping from a map of available SNPs. Blocks of markers in strong linkage disequilibrium (LD) are identified, and SNPs are selected to represent each block such that power to detect disease association with an underlying disease allele in LD with block members is preserved; all markers outside of blocks are also included in the tagging subset. A key, novel element of this method is that it incorporates information about the phase of LD observed among marker pairs to retain markers likely to be in coupling phase with an underlying disease locus, thus increasing power compared to a phase-blind approach. Power calculations illustrate important issues regarding LD phase and make clear the advantages of our approach to SNP selection. We apply our algorithm to genotype data from the International HapMap Consortium and demonstrate that considerable reduction in SNP genotyping may be attained while retaining much of the available power for a disease association screen. We also demonstrate that these tag SNPs effectively represent underlying variants not included in the LD analysis and SNP selection, by using leave-one-out tests to show that most (approximately 90%) of the "untyped" variants lying in blocks are in coupling-phase LD with a tag SNP. Additional performance tests using the HapMap ENCyclopedia of DNA Elements (ENCODE) regions show that the method compares well with the popular r2 bin tagging method. This work is a concrete example of how empirical LD phase may be used to benefit study design.

No association between the polymorphisms in matrix metalloproteinase-1 and matrix metalloproteinase-3 promoter regions and colorectal cancer in Chinese

PURPOSE

Matrix metalloproteinase-1 and matrix metalloproteinase-3 are implicated in all steps of cancer initiation, invasion, and metastasis. Recently, several genetic studies have demonstrated that matrix metalloproteinase-1-1607 ins/delG (1G/2G) polymorphism and matrix metalloproteinase-3-1612 ins/delA (5A/6A) polymorphism modify each transcriptional activity in allele-specific manners. In this study, we investigated whether these functional polymorphisms are associated with colorectal cancer in a Chinese population.

METHODS

Matrix metalloproteinase-1 and matrix metalloproteinase-3 genotypings were performed on 126 pathologically diagnosed colorectal cancer patients and 126 age-matched and gender-matched controls by polymerase chain reaction-based denaturing high-performance liquid chromatography analysis or restriction fragment length polymorphism, respectively.

RESULTS

The distributions of the matrix metalloproteinase-1 and matrix metalloproteinase-3 genotypes in healthy control subjects were inconsistent with Hardy-Weinberg equilibrium. Neither the genotype frequencies nor allele frequencies of matrix metalloproteinase-1 and matrix metalloproteinase-3 polymorphisms showed significant difference from those in healthy control subjects. There also were no significant associations between matrix metalloproteinase-1 and matrix metalloproteinase-3 genotypes and clinicopathologic features. When we examined the linkage disequilibrium between these two single nucleotide polymorphisms using expectation-maximization algorithm, we found that the two single nucleotide polymorphisms were in a strong linkage disequilibrium, but no significant difference was found in haplotype distribution.

CONCLUSIONS

Our present data suggest that the matrix metalloproteinase-1 and matrix metalloproteinase-3 promoter polymorphisms may not be useful markers to predicate susceptibility of colorectal cancer in Chinese.

Analysis of single nucleotide polymorphisms in genes in the chromosome 12Q24.31 region points to P2RX7 as a susceptibility gene to bipolar affective disorder

Previous results from our genetic analyses using pedigrees from a French Canadian population suggested that the interval delimited by markers on chromosome 12, D12S86 and D12S378, was the most probable genomic region to contain a susceptibility gene for affective disorders. Association studies with microsatellite markers using a case/control sample from the same population (n = 427) revealed significant allelic associations between the bipolar phenotype and marker NBG6. Since this marker is located in intron 9 of the P2RX7 gene, we analyzed the surrounding genomic region for the presence of polymorphisms in regulatory, coding and intron/exon junction sequences. Twenty four (24) SNPs were genotyped in a case/control sample and 12 SNPs in all pedigrees used for linkage analysis. Allelic, genotypic or family-based association studies suggest the presence of two susceptibility loci, the P2RX7 and CaMKK2 genes. The strongest association was observed in bipolar families at the non-synonymous SNP P2RX7-E13A (rs2230912, P-value = 0.000708), which results from an over-transmission of the mutant G-allele to affected offspring. This Gln460Arg polymorphism occurs at an amino acid that is conserved between humans and rodents and is located in the C-terminal domain of the P2X7 receptor, known to be essential for normal P2RX7 function.

Comparative study of the linkage disequilibrium of an ENCODE region, chromosome 7p15, in Korean, Japanese, and Han Chinese samples

The extent and pattern of linkage disequilibrium (LD) in the human genome provide important information for disease gene mapping. Previous studies have shown that LDs vary depending on chromosomal regions and populations. As the Asian samples of the International HapMap Project consisted of Japanese and Chinese populations, it was of interest whether we could use the HapMap data as a reference to carry out association studies of common complex diseases in a closely related population, such as Koreans. We have compared the LD and recombination patterns defined by single-nucleotide polymorphisms (SNPs) in ENCODE region ENm010, chromosome 7p15.2, in Korean, Japanese, and Chinese samples and further tested the robustness of tagSNPs among the Asian samples. We genotyped 792 SNPs in 500 kb (chromosome 7: 26699793-27199792, NCBI build 34) from 90 unrelated Koreans by fluorescence polarization detection and compared the data with Asian data from the HapMap project. Despite some differences in the position of high LD region boundaries, the overall patterns of LD were remarkably similar across the three samples, reflecting strong genetic affinities among them. Furthermore, the haplotype tag SNP transferability across the three samples was greater than 90%. Our results support the initial suggestion that the populations genotyped in the HapMap project might serve as reference populations for the selection of tagSNPs in association studies.

Localized breakdown in linkage disequilibrium does not always predict sperm crossover hot spots in the human MHC class II region

To investigate the relationship between meiotic crossover hot spots and block-like linkage disequilibrium (LD), we have extended our high-resolution studies of the human MHC class II region to a 90-kb segment upstream of the HLA-DOA gene. LD blocks in this region are not as well defined as in the neighboring 210-kb DNA segment but do show two regions of LD breakdown in which coalescent analysis indicates substantial historical recombination. Sperm crossover analysis of one region revealed a novel localized hot spot similar in intensity and morphology to most other MHC hot spots. Crossovers at this hot spot are not obviously affected by a large insertion/deletion polymorphism near the hot spot. The second region of LD breakdown, within the DPB1 gene, shows an extremely low level of sperm crossover activity and does not contain a sperm crossover hot spot. These results highlight the complexity of LD patterns and the importance of experimentally verifying crossover hot spots.

A genome-wide search for alleles and haplotypes associated with autism and related pervasive developmental disorders on the Faroe Islands

The involvement of genetic factors in the etiology of autism has been clearly established. We undertook a genome-wide search for regions containing susceptibility genes for autism in 12 subjects with childhood autism and related pervasive developmental disorders (PDDs) and 44 controls from the relatively isolated population of the Faroe Islands. In total, 601 microsatellite markers distributed throughout the human genome with an average distance of 5.80 cM were genotyped, including 502 markers in the initial scan. The Faroese population structure and genetic relatedness of cases and controls were also evaluated. Based on a combined approach, including an assumption-free test as implemented in CLUMP, Fisher's exact test for specific alleles and haplotypes, and IBD(0) probability calculations, we found association between autism and microsatellite markers in regions on 2q, 3p, 6q, 15q, 16p, and 18q. The most significant finding was on 3p25.3 (P(T1)=0.00003 and P(T4)=0.00007), which was also supported by other genetic studies. Furthermore, no evidence of population substructure was found, and a higher degree of relatedness among cases could not be detected, decreasing the risk of inflated P-values. Our data suggest that markers in these regions are in linkage disequilibrium with genes involved in the etiology of autism, and we hypothesize susceptibility genes for autism and related PDDs to be localized within these regions.

High-density single-nucleotide polymorphism maps of the human genome

Here we report a large, extensively characterized set of single-nucleotide polymorphisms (SNPs) covering the human genome. We determined the allele frequencies of 55,018 SNPs in African Americans, Asians (Japanese-Chinese), and European Americans as part of The SNP Consortium's Allele Frequency Project. A subset of 8333 SNPs was also characterized in Koreans. Because these SNPs were ascertained in the same way, the data set is particularly useful for modeling. Our results document that much genetic variation is shared among populations. For autosomes, some 44% of these SNPs have a minor allele frequency > or =10% in each population, and the average allele frequency differences between populations with different continental origins are less than 19%. However, the several percentage point allele frequency differences among the closely related Korean, Japanese, and Chinese populations suggest caution in using mixtures of well-established populations for case-control genetic studies of complex traits. We estimate that approximately 7% of these SNPs are private SNPs with minor allele frequencies <1%. A useful set of characterized SNPs with large allele frequency differences between populations (>60%) can be used for admixture studies. High-density maps of high-quality, characterized SNPs produced by this project are freely available.

Study of six X-linked tetranucleotide microsatellites: population data from five Spanish regions

We studied six X-linked microsatellites in a large group of Spanish individuals (n=614) from five different regions located in northern, central and southern Spain. All the markers had tetranucleotide repeat units (DXS9895, DXS9898, DXS7130, DXS7132, GATA172D05 and DXS6789). They were amplified in two triplex PCR reactions. There were no significant sex- or region-related differences in allelic frequencies, suggesting that general national databases can be adequate as a reference for X-linked markers. The analysis of those six short tandem repeats combined in 316 males revealed 300 different "temporary haplotypes", 283 of which were found only once. There was no evidence for statistically significant linkage disequilibrium among the loci studied. Therefore these markers are quite polymorphic and useful for forensic purposes.

Perspectives on human genetic variation from the HapMap Project

The completion of the International HapMap Project marks the start of a new phase in human genetics. The aim of the project was to provide a resource that facilitates the design of efficient genome-wide association studies, through characterising patterns of genetic variation and linkage disequilibrium in a sample of 270 individuals across four geographical populations. In total, over one million SNPs have been typed across these genomes, providing an unprecedented view of human genetic diversity. In this review we focus on what the HapMap Project has taught us about the structure of human genetic variation and the fundamental molecular and evolutionary processes that shape it.

Nucleotide variations in genes encoding plasminogen activator inhibitor-2 and serine proteinase inhibitor B10 associated with prostate cancer

Genes encoding the serine proteinase inhibitor B family (SERPINBs) are mainly clustered on human chromosome 18 (18q21). Several serpins are known to affect malignant phenotypes of tumor cells, so aberrant genetic variants in this molecular family are candidates for conferring susceptibility for risk of cancer. We investigated whether eight selected non-synonymous variations within SERPINB loci at 18q21 might be associated with risk of prostate cancer in Japanese men. A case-control study involving 292 prostate-cancer patients and 384 controls revealed significant differences in regard to distribution of four missense variations in genes encoding plasminogen activator inhibitor 2 (PAI2) and SERPINB10. The most significant association was detected for the N120D polymorphism in the PAI2 gene (P = 5.0 x 10(-5)); men carrying the 120-N allele (120-N/N and 120-N/D genotypes) carried a 2.4-fold increased risk of prostate cancer (95% confidence interval 1.45-4.07). Associations were also detected for three other missense polymorphisms in those two genes. Strong linkage disequilibrium in the region encompassing PAI2 and SERPINB10 extended to about 50 kbp. The results suggested that missense variations in one or both of these genes confer important risks for prostate cancer, and may be themselves tumorigenic. Although confirmative replication studies on larger cohorts are awaited, clinical examination of these variations may become useful for identifying individuals at high risk for prostate cancer.

Association study of polymorphisms in synaptic vesicle-associated genes, SYN2 and CPLX2, with schizophrenia

BACKGROUND

The occurrence of aberrant functional connectivity in the neuronal circuit is one of the integrative theories of the etiology of schizophrenia. Previous studies have reported that the protein and mRNA levels of the synapsin 2 (SYN2) and complexin 2 (CPLX2) genes were decreased in patients with schizophrenia. Synapsin 2 and complexin 2 are involved in synaptogenesis and the modulation of neurotransmitter release. This report presents a study of the association of polymorphisms of SYN2 and CPLX2 with schizophrenia in the Korean population.

METHODS

Six single nucleotide polymorphisms (SNPs) and one 5-bp insertion/deletion in SYN2 and five SNPs in CPLX2 were genotyped in 154 Korean patients with schizophrenia and 133 control patients using direct sequencing or restriction fragment length polymorphism analysis. An intermarker linkage disequilibrium map was constructed for each gene.

RESULTS

Although there was no significant difference in the genotypic distributions and allelic frequencies of either SYN2 or CPLX2 polymorphisms between the schizophrenia and control groups, the two-way haplotype analyses revealed significant associations with the disease (P < 0.05 after Bonferroni correction). The three-way haplotype analyses also revealed a significant association of SYN2 with schizophrenia (P < 0.001 after Bonferroni correction).

CONCLUSION

These results suggest that both SYN2 and CPLX2 may confer susceptibility to schizophrenia in the Korean population.

Linkage disequilibrium patterns vary substantially among populations

A major initiative to create a global human haplotype map has recently been launched as a tool to improve the efficiency of disease gene mapping. The 'HapMap' project will study common variants in depth in four (and to a lesser degree in up to 12) populations to catalogue haplotypes that are expected to be common to all populations. A hope of the 'HapMap' project is that much of the genome occurs in regions of limited diversity such that only a few of the SNPs in each region will capture the diversity and be relevant around the world. In order to explore the implications of studying only a limited number of populations, we have analyzed linkage disequilibrium (LD) patterns of three 175-320 kb genomic regions in 16 diverse populations with an emphasis on African and European populations. Analyses of these three genomic regions provide empiric demonstration of marked differences in frequencies of the same few haplotypes, resulting in differences in the amount of LD and very different sets of haplotype frequencies. These results highlight the distinction between the statistical concept of LD and the biological reality of haplotypes and their frequencies. The significant quantitative and qualitative variation in LD among populations, even for populations within a geographic region, emphasizes the importance of studying diverse populations in the HapMap project to assure broad applicability of the results.

Complex HTR2C linkage disequilibrium and promoter associations with body mass index and serum leptin

The occurrence of obesity, eating disorders, and related diseases has increased in many parts of the world. Given that few strong genetic factors have been found, it is clear that these are complex multi-factorial diseases. The serotonin receptor 2C, a member of the 5-HTergic system, has been implicated in the control of phagia and obesity. We report a detailed investigation of linkage disequilibrium (LD) within and between the HTR2C promoter and the flanking sequences around a commonly utilized marker in the second coding exon of HTR2C. We suggest that inconsistent associations between HTR2C and several phenotypes, including obesity, may be due to the LD pattern across the gene in which recombination and gene conversion have been influential. The nucleotide and haplotype distribution is consistent with that of the neutral mutation model. The number of haplotypes suggests demographic influences or over dominant selection that may have a function in HTR2C expression. Using the fine LD pattern, we describe a possible association with promoter haplotypes and diplotypes, including a GT microsatellite, and body mass index (BMI) > or =30 kgm(-2) (P<0.0001). SNP -995G>A heterozygotes, as well as promoter diplotypes, were found to marginally influence higher serum leptin corrected for percentage body fat (P=0.01), which might suggest that these subjects are leptin resistant. Our results complement previous studies of HTR2C in both mice and humans, and suggest the importance of genetic variation and elucidating the fine LD structure in uncovering the genetic factors of obesity.

Experimental designs for reliable detection of linkage disequilibrium in unstructured random population association studies

A method is given for design of experiments to detect associations (linkage disequilibrium) in a random population between a marker and a quantitative trait locus (QTL), or gene, with a given strength of evidence, as defined by the Bayes factor. Using a version of the Bayes factor that can be linked to the value of an F-statistic with an existing deterministic power calculation makes it possible to rapidly evaluate a comprehensive range of scenarios, demonstrating the feasibility, or otherwise, of detecting genes of small effect. The Bayes factor is advocated for use in determining optimal strategies for selecting candidate genes for further testing or applications. The prospects for fine-scale mapping of QTL are reevaluated in this framework. We show that large sample sizes are needed to detect small-effect genes with a respectable-sized Bayes factor, and to have good power to detect a QTL allele at low frequency it is necessary to have a marker with similar allele frequency near the gene.

Functional analysis of the single nucleotide polymorphism (787T>C) in the tissue-nonspecific alkaline phosphatase gene associated with BMD

Polymorphisms of the TNSALP gene have not previously been studied as a possible determinant for variations in BMD or as a predisposing genetic factor for osteoporosis. This study showed a significantly higher association between the 787T>C (Tyr246His) TNSALP gene and BMD among 501 postmenopausal women. Furthermore, the effects of amino acid substitution on the catalytic property of the protein translated from the 787T>C gene were examined.

INTRODUCTION

Alkaline phosphatase (ALP) is present mainly on the cell membrane in various tissues and hydrolyzes a variety of monophosphate esters into inorganic phosphoric acid and alcohol. Human ALPs are classified into four types: tissue-nonspecific, intestinal, placental, and germ cell types. Based on studies of hypophosphatasia, which is a systemic skeletal disorder resulting from a tissue-nonspecific ALP (TNSALP) deficiency, TNSALP was suggested to be indispensable for bone mineralization.

MATERIALS AND METHODS

We explored the possibility that the TNSALP gene may contribute to age-related bone loss in humans by examining the association between TNSALP gene polymorphisms and BMD in 501 Japanese postmenopausal women. To analyze the protein translated from the TNSALP gene associated with BMD, we constructed a TNSALP cDNA expression plasmid.

RESULTS

We genotyped two single nucleotide polymorphisms (787T>C[Tyr246His] and 876A>G[Pro275Pro]), which proved to be in complete linkage disequilibrium. There was a significant difference in BMD and the BMD score adjusted for age and body weight (Z score) among haplotypes (p = 0.041), which was lowest among 787T/876A homozygotes, highest among 787T>C/876A>G homozygotes, and intermediate among heterozygotes. In subgroups divided by age, haplotypes were significantly associated with BMD in older postmenopausal women (>74 years; p = 0.001), but not in younger postmenopausal women (<74 years; p = 0.964). Expression of the 787T>C TNSALP gene using COS-1 cells showed that the protein translated from 787T>C had ALP-specific activity similar to that of 787T. Interestingly, the K(m) value for TNSALP in cells transfected with the 787T>C TNSALP gene was decreased significantly compared with that of cells bearing the 787T gene, reflecting the higher affinity.

CONCLUSIONS

These results suggest that variation in TNSALP may be an important determinant of age-related bone loss in humans and that the phosphate metabolism pathway may provide a novel target for the prevention and treatment of osteoporosis.

Characterisation of the genomic architecture of human chromosome 17q and evaluation of different methods for haplotype block definition

BACKGROUND

The selection of markers in association studies can be informed through the use of haplotype blocks. Recent reports have determined the genomic architecture of chromosomal segments through different haplotype block definitions based on linkage disequilibrium (LD) measures or haplotype diversity criteria. The relative applicability of distinct block definitions to association studies, however, remains unclear. We compared different block definitions in 6.1 Mb of chromosome 17q in 189 unrelated healthy individuals. Using 137 single nucleotide polymorphisms (SNPs), at a median spacing of 15.5 kb, we constructed haplotype block maps using published methods and additional methods we have developed. Haplotype tagging SNPs (htSNPs) were identified for each map.

RESULTS

Blocks were found to be shorter and coverage of the region limited with methods based on LD measures, compared to the method based on haplotype diversity. Although the distribution of blocks was highly variable, the number of SNPs that needed to be typed in order to capture the maximum number of haplotypes was consistent.

CONCLUSION

For the marker spacing used in this study, choice of block definition is not important when used as an initial screen of the region to identify htSNPs. However, choice of block definition has consequences for the downstream interpretation of association study results.

Association of nucleotide variations in the apolipoprotein B48 receptor gene (APOB48R) with hypercholesterolemia

Factors predisposing to the phenotypic features of high total cholesterol (T-Cho) in human plasma have not been clearly defined. Here we report an association between two variations in the apolipoprotein B48 receptor gene (APOB48R) and plasma T-Cho levels among 352 adult individuals in Japan. By analyzing phenotypic associations between age- and gender-adjusted levels of plasma T-Cho, low-density lipoprotein (LDL) cholesterol (LDL-C), and high-density lipoprotein (HDL) cholesterol (HDL-C), we detected a significant correlation between genotypes of the A419P variation and adjusted T-Cho levels. Among homozygous G-allele carriers (n=265), heterozygous carriers (n=78), and homozygous minor C-allele carriers (n=9), T-Cho levels were 2.43+/-0.21 mg/cm(3), 2.48+/-0.24 mg/cm(3), and 2.63+/-0.21 mg/cm(3), respectively, indicating a codominant T-Cho-elevating effect of the minor C-allele (r=0.15, P=0.007). A similar effect was detected for c.934-960/del (r=0.13, P=0.015). Linkage disequilibrium (LD) analysis detected significant LD among eight variant sites that included neighboring loci. Our results indicate that variations in APOB48R and nearby genes are among the many factors involved in hypercholesterolemia. The etiological studies should now include consideration of this novel aspect of the mechanism(s) leading to hypercholesterolemic disease.

Linkage disequilibrium maps constructed with common SNPs are useful for first-pass disease association screens

To develop an efficient strategy for mapping genetic factors associated with common diseases, we constructed linkage disequilibrium (LD) maps of human chromosomes 5, 7, 17, and X. These maps consist of common single nucleotide polymorphisms at an average intermarker distance of 100 kb. The genotype data from these markers in a panel of American samples of European descent were analyzed to produce blocks of markers in strong pair-wise LD. Power calculations were used to guide block definitions and predicted that high-level LD maps would be useful in initial genome scans for susceptibility alleles in case-control association studies of complex diseases. As anticipated, LD blocks on the X chromosome were larger and covered more of the chromosome than those found on the autosomes.

Linkage disequilibrium and association studies in higher plants: present status and future prospects

During the last two decades, DNA-based molecular markers have been extensively utilized for a variety of studies in both plant and animal systems. One of the major uses of these markers is the construction of genome-wide molecular maps and the genetic analysis of simple and complex traits. However, these studies are generally based on linkage analysis in mapping populations, thus placing serious limitations in using molecular markers for genetic analysis in a variety of plant systems. Therefore, alternative approaches have been suggested, and one of these approaches makes use of linkage disequilibrium (LD)-based association analysis. Although this approach of association analysis has already been used for studies on genetics of complex traits (including different diseases) in humans, its use in plants has just started. In the present review, we first define and distinguish between LD and association mapping, and then briefly describe various measures of LD and the two methods of its depiction. We then give a list of different factors that affect LD without discussing them, and also discuss the current issues of LD research in plants. Later, we also describe the various uses of LD in plant genomics research and summarize the present status of LD research in different plant genomes. In the end, we discuss briefly the future prospects of LD research in plants, and give a list of softwares that are useful in LD research, which is available as electronic supplementary material (ESM).