Spectrum of nonrandom associations between microsatellite loci on human chromosome 11p15

Genetic Diversity and Population Structure of a Large USDA Sesame Collection

Sesame, Sesamum indicum L., is one of the oldest domesticated crops used for its oil and protein in many parts of the world. To build genomic resources for sesame that could be used to improve sesame productivity and responses to stresses, a USDA sesame germplasm collection of 501 accessions originating from 36 countries was used in this study. The panel was genotyped using genotyping-by-sequencing (GBS) technology to explore its genetic diversity and population structure and the relatedness among its accessions. A total of 24,735 high-quality single-nucleotide polymorphism (SNP) markers were identified over the 13 chromosomes. The marker density was 1900 SNP per chromosome, with an average polymorphism information content (PIC) value of 0.267. The marker polymorphisms and heterozygosity estimators indicated the usefulness of the identified SNPs to be used in future genetic studies and breeding activities. The population structure, principal components analysis (PCA), and unrooted neighbor-joining phylogenetic tree analyses classified two distinct subpopulations, indicating a wide genetic diversity within the USDA sesame collection. Analysis of molecular variance (AMOVA) revealed that 29.5% of the variation in this population was due to subpopulations, while 57.5% of the variation was due to variation among the accessions within the subpopulations. These results showed the degree of differentiation between the two subpopulations as well as within each subpopulation. The high fixation index (FST) between the distinguished subpopulations indicates a wide genetic diversity and high genetic differentiation among and within the identified subpopulations. The linkage disequilibrium (LD) pattern averaged 161 Kbp for the whole sesame genome, while the LD decay ranged from 168 Kbp at chromosome LG09 to 123 Kbp in chromosome LG05. These findings could explain the complications of linkage drag among the traits during selections. The selected accessions and genotyped SNPs provide tools to enhance genetic gain in sesame breeding programs through molecular approaches.

Population genetic structure of Guiana dolphin (Sotalia guianensis) from the southwestern Atlantic coast of Brazil

Sotalia guianensis is a small dolphin that is vulnerable to anthropogenic impacts. Along the Brazilian Atlantic coast, this species is threatened with extinction. A prioritized action plan for conservation strategies relies on increased knowledge of the population. The scarcity of studies about genetic diversity and assessments of population structure for this animal have precluded effective action in the region. Here, we assessed, for the first time, the genetic differentiation at 14 microsatellite loci in 90 S. guianensis specimens stranded on the southeastern Atlantic coast of the State of Espírito Santo, Brazil. We estimated population parameters and structure, measured the significance of global gametic disequilibrium and the intensity of non-random multiallelic interallelic associations and constructed a provisional synteny map using Bos taurus, the closest terrestrial mammal with a reference genome available. All microsatellite loci were polymorphic, with at least three and a maximum of ten alleles each. Allele frequencies ranged from 0.01 to 0.97. Observed heterozygosity ranged from 0.061 to 0.701. The mean inbreeding coefficient was 0.103. Three loci were in Hardy-Weinberg disequilibrium even when missing genotypes were inferred. Although 77 of the 91 possible two-locus associations were in global gametic equilibrium, we unveiled 13 statistically significant, sign-based, non-random multiallelic interallelic associations in 10 two-locus combinations with either coupling (D' values ranging from 0.782 to 0.353) or repulsion (D' values -0.517 to -1.000) forces. Most of the interallelic associations did not involve the major alleles. Thus, for either physically or non-physically linked loci, measuring the intensity of non-random interallelic associations is important for defining the evolutionary forces at equilibrium. We uncovered a small degree of genetic differentiation (F_ST = 0.010; P-value = 0.463) with a hierarchical clustering into one segment containing members from the southern and northern coastal regions. The data thus support the scenario of little genetic structure in the population of S. guianensis in this geographic area.

Genetic epidemiology of osteoporosis across four microsatellite markers near the VDR gene

The large amount of positive genetic association data in a number of bone diseases suggests functional consequences of Vitamin D receptor (VDR) gene polymorphism. In the present study, four microsatellite markers viz., D12S1633, D12S1635, D12S347, and D12S96, that lie in the vicinity of the VDR gene on chromosome 12 were selected to assess the allele distribution pattern and diversity among three groups of individuals - normal, osteopenia and osteoporosis. Genetic association study was performed using allele frequency data. Total genomic DNA was isolated from the whole blood of 226 individuals, after recording their bone mineral density (BMD) using Dual X-ray absorptiometry (DXA). All DNA samples were subjected to multiplex Polymerase Chain Reaction (PCR) - genotyping. Allele frequencies and genetic diversity parameters like - number of alleles, average variance and average heterozygosity across all the four markers among three groups were computed. Effect of population stratification was excluded by investigating population structure. A trend of decreasing genetic diversity across four loci from normal to pre- and post-disease condition has been observed. Lesser recombination rate (θ) indicates linkage between studied microsatellite markers and VDR gene. Statistically significant linkage disequilibrium was detected for the allele - 22 of locus D12S96 with osteoporosis. A positive association of allele - 22 suggests susceptibility to disease whereas predominance of allele - 27 among non - diseased group implicates its association with normal bone health.

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.

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.

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.

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.

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.

Multilocus genomics of outcrossing plant populations

The structure and organization of natural plant populations can be understood by estimating the genetic parameters related to mating behavior, recombination frequency, and gene associations with DNA-based markers typed throughout the genome. We developed a statistical and computational model for estimating and testing these parameters from multilocus data collected in a natural population. This model, constructed by a maximum likelihood approach and implemented within the EM algorithm, is shown to be robust for simultaneously estimating the outcrossing rate, recombination frequencies and linkage disequilibria. The algorithm built with three or more markers allows the characterization of crossover interference in meiosis and high-order disequilibria among different genes, thus providing a powerful tool for illustrating a detailed picture of genetic diversity and organization in natural populations. Computer simulations demonstrate the statistical properties of the proposed model. This multilocus model will be useful for studying the pattern and amount of genetic variation within and among populations to further infer the evolutionary history of a plant species.

Mapping of susceptibility and protective loci for acute GVHD in unrelated HLA-matched bone marrow transplantation donors and recipients using 155 microsatellite markers on chromosome 22

Despite matching donors and recipients for the human leukocyte antigens (HLAs) expressed by the major histocompatibility genomic region of the short arm of chromosome 6, several recipients still develop acute graft-versus-host disease (aGVHD) after bone marrow transplantation (BMT). This is possibly due to non-HLA gene polymorphisms, such as minor histocompatibility antigens (mHas) and genes coding for cytokines. However, a detailed genetic background for aGVHD has not yet been established. To find novel susceptibility and/or protective loci for aGVHD, a whole genome-wide association study of donors and recipients needs to be performed. As the first step to such a study, we retrospectively analyzed polymorphisms of 155 microsatellite markers spread across the long arm of chromosome 22 in 70 pairs of HLA-matched unrelated BMT donors and recipients. We performed individual typing and then compared the markers' allele frequencies (1) between all the aGVHD (grades III and IV GVHD) and GVHD-free (grade 0 GVHD) groups in donors and recipients and (2) between the aGVHD and aGVHD-free groups in donor/recipient pairs that were matched and mismatched for the microsatellite marker's allele. Screening of the microsatellite markers revealed five loci with a significant difference between the aGVHD and GVHD-free groups and revealed eight loci on chromosome 22, where the microsatellite allele mismatched markers were associated with aGVHD. This screening analysis suggests that several aGVHD-associated susceptible and protective loci exist on chromosome 22, which may encompass novel gene regions that need to be elucidated for their role in aGVHD.

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

Background

Various software tools are available for the display of pairwise linkage disequilibrium across multiple single nucleotide polymorphisms. The HapMap project also presents these graphics within their website. However, these approaches are limited in their use of data from multiallelic markers and provide limited information in a graphical form.

Results

We have developed a software package (MIDAS – Multiallelic Interallelic Disequilibrium Analysis Software) for the estimation and graphical display of interallelic linkage disequilibrium. Linkage disequilibrium is analysed for each allelic combination (of one allele from each of two loci), between all pairwise combinations of any type of multiallelic loci in a contig (or any set) of many loci (including single nucleotide polymorphisms, microsatellites, minisatellites and haplotypes). Data are presented graphically in a novel and informative way, and can also be exported in tabular form for other analyses. This approach facilitates visualisation of patterns of linkage disequilibrium across genomic regions, analysis of the relationships between different alleles of multiallelic markers and inferences about patterns of evolution and selection.

Conclusion

MIDAS is a linkage disequilibrium analysis program with a comprehensive graphical user interface providing novel views of patterns of linkage disequilibrium between all types of multiallelic and biallelic markers.

Availability

Available from and

A study of TH01 and IGF2-INS-TH haplotypes in relation to smoking initiation in three independent surveys

OBJECTIVES

Recent studies suggest an association between a microsatellite locus (TH01) located in intron 1 of the tyrosine hydroxylase gene (TH) and nicotine dependence. We aimed here to study whether both TH01 and haplotypes of the wider IGF2-INS-TH region influence initiation of regular smoking in current smokers.

METHODS

A total of 3637 individuals from three independent studies (two of adults and one of adolescents) were analysed in relation to the age of first regular smoking (AFRS). Haplotypes and genotypes were obtained for the polymorphisms TH01, IGF2 ApaI, INS HphI and DRD4 VNTR (48 bp)n. Association between IGF2-INS-TH haplotypes and AFRS was tested by a regression model. A genotype-based genetic model assuming additivity was followed in order to estimate the effect of individual loci.

RESULTS

Overall, no significant associations were found after correcting for multiple tests. However, an IGF2-INS-TH haplotype (*5) was found to be nominally associated with AFRS at younger ages in adult smokers. Analyses of individual loci points to TH01 as a possible candidate influencing initiation of regular smoking. An AFRS-lowering trend nominally associated with allele 9 in a dosage-dependent manner was identified in both adult cohorts. TH01 did not show association or trend with age of initiation (first puff) either in adolescents or in the adolescents smoking regularly at age 18.

CONCLUSION

This study adds to the genetic evaluation of the associations of TH01 with smoking predisposition. Differences between historical and prospective surveys, different biological pathways and possible functional roles of this microsatellite in smoking initiation are discussed.

Host microsatellite alleles in malaria predisposition?

BACKGROUND

Malaria is a serious, sometimes fatal, disease caused by Plasmodium infection of human red blood cells. The host-parasite co-evolutionary processes are well understood by the association of coding variations such as G6PD, Duffy blood group receptor, HLA, and beta-globin gene variants with malaria resistance. The profound genetic diversity in host is attributed to polymorphic microsatellites loci. The microsatellite alleles in bacterial species are known to have aided their survival in fatal environmental conditions. The fascinating question is whether microsatellites are genomic cushion in the human genome to combat disease stress and has cause-effect relationships with infections.

PRESENTATION OF THE HYPOTHESIS

It is hypothesized that repeat units or alleles of microsatellites TH01 and D5S818, located in close proximity to beta-globin gene and immune regulatory region in human play a role in malaria predisposition. Association of alleles at aforesaid microsatellites with malaria infection was analysed. To overrule the false association in unrecognized population stratification, structure analysis and AMOVA were performed among the sampled groups.

TESTING OF HYPOTHESIS

Associations of microsatellite alleles with malaria infection were verified using recombination rate, Chi-square, and powerful likelihood tests. Further investigation of population genetic structure, and AMOVA was done to rule out the confounding effects of population stratification in interpretation of association studies.

IMPLICATION OF THE HYPOTHESIS

Lower recombination rate (theta) between microsatellites and genes implicated in host fitness; positive association between alleles-13 (D5S818), 9 (TH01) and strong susceptibility to Plasmodium falciparum; and alleles-12 (D5S818) and 6 (TH01) rendering resistance to human host were evident. The interesting fact emerging from the study was that while predisposition to malaria was a prehistoric attribute, among TH01 alleles; evolution of resistant allele-6 was a recent phenomenon, which could conceivably be driven by infection related selective forces. The host's microsatellite allelic associations with malaria infection were valid in the light of low genetic variance between sampled groups and no population stratification.

Association of tumor necrosis factor-alpha polymorphisms and ozone-induced change in lung function

Ozone is a major air pollutant with adverse health effects which exhibit marked inter-individual variability. In mice, regions of genetic linkage with ozone-induced lung injury include the tumor necrosis factor-alpha (TNF), lymphotoxin-alpha (LTA), Toll-like receptor 4 (TLR4), superoxide dismutase (SOD2), and glutathione peroxidase (GPX1) genes. We genotyped polymorphisms in these genes in 51 individuals who had undergone ozone challenge. Mean change in FEV1 with ozone challenge, as a percentage of baseline, was -3% in TNF -308G/A or A/A individuals, compared with -9% in G/G individuals (p = 0.024). When considering TNF haplotypes, the smallest change in FEV1 with ozone exposure was associated with the TNF haplotype comprising LTA +252G/TNF -1031T/TNF -308A/TNF -238G. This association remained statistically significant after correction for age, sex, disease, and ozone concentration (p = 0.047). SOD2 or GPX1 genotypes were not associated with lung function, and the TLR4 polymorphism was too infrequent to analyze. The results of this study support TNF as a genetic factor for susceptibility to ozone-induced changes in lung function in humans, and has potential implications for stratifying health risks of air pollution.

Genetic variability of the marine mussel Mytilus galloprovincialis assessed using two-dimensional electrophoresis

Two-dimensional electrophoresis (2-DE) has been used to measure the degree of genetic variability of the marine mussel Mytilus galloprovincialis. Genetic polymorphisms were detected in 33 of a total of 86 polypeptides scored among the most abundant proteins from foot samples in 38 individuals. Estimates of average heterozygosity were 0.101+/-0.018 and 0.114+/-0.021 in a natural and a cultured population, respectively, from the NW of the Iberian Peninsula. These are the highest estimates of average heterozygosity reported by 2-DE in an animal species to date. We consider that these data throw open the question of the level of genetic variability detectable by two-dimensional electrophoresis. Multilocus genotype data were used to infer haplotypic frequencies by means of the EM algorithm in order to detect linkage disequilibrium between loci coding abundant proteins. Significant associations were found in 22.7% of the 406 two-locus pairs analysed. Also, clusters of loci in which all pairwise combinations exhibit statistically significant associations were detected and physical linkage between some of these loci is postulated from the linkage disequilibrium data.

Distribution of nonrandom associations between pairs of protein loci along the third chromosome of Drosophila melanogaster

The within-chromosome distribution of gametic disequilibrium (GD) between protein loci, and the underlying evolutionary factors of this distribution, are still largely unknown. Here, we report a detailed study of GD between a large number of protein loci (15) spanning 87% of the total length of the third chromosome of Drosophila melanogaster in a large sample of haplotypes (600) drawn from a single natural population. We used a sign-based GD estimation method recently developed for multiallelic systems, which considerably increases both the statistical power and the accuracy of estimation of the intensity of GD. We found that strong GD between pairs of protein loci was widespread throughout the chromosome. In total, 22% of both the pairs of alleles and pairs of loci were in significant GD, with mean intensities (as measured by D' coefficients) of 0.43 and 0.31, respectively. In addition, strong GD often occurs between loci that are far apart. By way of illustration, 32% of the allele pairs in significant GD occurred within pairs of loci separated by effective frequencies of recombination (EFRs) of 15-20 cM, the mean D' value being 0.49. These observations are in sharp contrast with previous studies showing that GD between protein loci is rarely found in natural populations of outcrossing species, even between very closely linked loci. Interestingly, we found that most instances of significant interallelic GD (68%) involved functionally related protein loci. Specifically, GD was markedly more frequent between protein loci related by the functions of hormonal control, molybdenum control, antioxidant defense system, and reproduction than between loci without known functional relationship, which is indicative of epistatic selection. Furthermore, long-distance GD between functionally related loci (mean EFR 9 cM) suggests that epistatic interactions must be very strong along the chromosome. This evidence is hardly compatible with the neutral theory and has far-reaching implications for understanding the multilocus architecture of the functional genome. Our findings also suggest that GD may be a useful tool for discovering networks of functionally interacting proteins.

Typing dinucleotide repeats under nondenaturing conditions with single-base resolution and high sizing precision

Dinucleotide repeats are genetic markers that are useful for many purposes, including genetic epidemiology, population genetics, and genetic diagnostics. The accuracy of analyses based on dinucleotide repeat polymorphisms is highly dependent on the success achieved in minimizing genotyping errors. Genotyping errors in dinucleotide repeat typing may arise for various reasons, including polymerase chain reaction (PCR) processing errors and the use of unsuitable electrophoretic conditions for resolving amplification products (i.e., lack of single-base resolution and inadequate precision in allele sizing). We have recently described a nondenaturing electrophoretic system useful for detecting PCR processing errors that lead to misidentification of heterozygotes as homozygotes in (AC)n repeat typing. Here, we show that this system also allows resolution of (AC)n repeats in native conditions with single-base resolution and high sizing precision, on the basis of an analysis of seven human (AC)n repeats ranging in size from 72 to 217 bp. This PAGE system is thus also useful for reducing the likelihood both of allele misidentification due to the absence of single-base resolution and of inaccuracies in allele sizing due to anomalous electrophoretic migrations among the alleles within an (AC)n repeat.