Novel linkage mapping approach using DNA pooling in human and animal genetics. I. Detection of complex disease loci

DNA pooling is a potential methodology for genetic loci with small effect contributing to complex diseases and quantitative traits. This is accomplished by the rapid preliminary screening of the genome for the allelic association with the most common class of polymorphic short tandem repeat markers. The methodology assumes as a common founder for the linked disease locus of interest and searches for a region of a chromosome shared between affected individuals. The general theory of DNA pooling basically relies on the observed differences in the allelic distribution between pools from affected and unaffected individuals, including a reduction in the number of alleles in the affected pool, which indicate the sharing of a chromosomal region. The power of statistic for associated linkage mapping can be determined using two recently developed strategies, firstly, by measuring the differences of allelic image patterns produced by two DNA pools of extreme character and secondly, by measuring total allele content differences by comparing between two pools containing large numbers of DNA samples. These strategies have effectively been utilized to identify the shared chromosomal regions for linkage studies and to investigate the candidate disease loci for fine structure gene mapping using allelic association. This paper outlines the utilization of DNA pooling as a potential tool to locate the complex disease loci, statistical methods for accurate estimates of allelic frequencies from DNA pools, its advantages, drawbacks and significance in associate linkage mapping using pooled DNA samples.     

VNTR loci reveal differentiation between and structure within populations of the eastern barred bandicoot Perameles gunnii

The Eastern Barred Bandicoot Perameles gunnii has declined in abundance within mainland south-eastern Australia, to a relict wild population of less than 100 individuals in Hamilton, Victoria. It is more common, but is also declining in Tasmania. Genomic DN A variability was compared within and between surviving populations of P. gunnii using variable number of tandem repeat (VNTR) markers in one of two ways. First, average percentage differences (APDs) were determined between profiles for two VNTR probe—endonuclease combinations. Secondly, because one of these combinations revealed two multiallelic VNTR loci, genotypes were assigned and analysed for homogeneity of allele frequencies among subpopulations, for deviation of heterozygosity from Hardy-Weinberg equilibrium within populations and for genetic structuring among individuals from different subpopulations. The results of both the APD and defined locus approaches showed consistent trends within and between populations. Genetic variability was higher among mainland P. gunnii than in Tasmanian populations (higher APDs, number of alleles, and heterozygosity at one locus), despite the known decline and subdivision of the Hamilton population. Eleven per cent of the variability detected in Hamilton was attributed to genetic differentiation between east and west subdivisions of the population. Departure from random mating indicating local inbreeding within collecting localities was evident for one locus in both north and south Tasmania, particularly at one locality. AH alleles at both loci were unique to either Hamilton or Tasmanian P. gunnii. The initial captive colony contains high heterozygosity for these loci. It is concluded that VNTR markers can be of benefit for use in studies of population differentiation and for conservation management.     

Disease swamps molecular signatures of genetic-environmental associations to abiotic factors in Tasmanian devil (Sarcophilus harrisii) populations

Landscape genomics studies focus on identifying candidate genes under selection via spatial variation in abiotic environmental variables, but rarely by biotic factors (i.e., disease). The Tasmanian devil (Sarcophilus harrisii) is found only on the environmentally heterogeneous island of Tasmania and is threatened with extinction by a transmissible cancer, devil facial tumor disease (DFTD). Devils persist in regions of long-term infection despite epidemiological model predictions of species’ extinction, suggesting possible adaptation to DFTD. Here, we test the extent to which spatial variation and genetic diversity are associated with the abiotic environment (i.e., climatic variables, elevation, vegetation cover) and/or DFTD. We employ genetic-environment association analyses using 6886 SNPs from 3287 individuals sampled pre- and post-disease arrival across the devil's geographic range. Pre-disease, we find significant correlations of allele frequencies with environmental variables, including 365 unique loci linked to 71 genes, suggesting local adaptation to abiotic environment. The majority of candidate loci detected pre-DFTD are not detected post-DFTD arrival. Several post-DFTD candidate loci are associated with disease prevalence and were in linkage disequilibrium with genes involved in tumor suppression and immune response. Loss of apparent signal of abiotic local adaptation post-disease suggests swamping by strong selection resulting from the rapid onset of DFTD.     

The Human Obesity Gene Map: The 1997 Update

An update of the human obesity gene map incorporating published results up to October 1997 is presented. Evidence from Mendelian disorders exhibiting obesity as a clinical feature; single-gene mutation rodent models; quantitative trait loci uncovered in human genome-wide scans and in crossbreeding experiments with mouse, rat, and pig models; association and case-control studies with candidate genes; and linkage studies with genes and other markers is reviewed. All chromosomal locations of the animal loci are converted into human genome locations based on syntenic relationships between the genomes. A complete listing of all of these loci reveals that all but chromosome Y of the 24 human chromosomes are represented. Some chromosomes show at least three putative loci related to obesity on both arms (1, 2, 6, 8, 11, and 20) and several on one chromosome arm only (3p, 4q, 5q, 7q, 12q, 13q, 15q, 15p, 22q, and Xq). Studies reporting negative association and linkage results are also listed, with the exception of the unlinked markers from genome-wide scans.     

Mapping the gene for hereditary hyperparathyroidism and prolactinoma (MEN1Burin) to chromosome 11q: evidence for a founder effect in patients from Newfoundland.

An autosomal dominant syndrome of prolactinomas, carcinoids, and hyperparathyroidism was described in four Newfoundland kindreds in 1980 and in one kindred from the Pacific Northwest in 1983. Because this syndrome shares many features with multiple endocrine neoplasia type 1, the gene for which maps to proximal chromosome 11q, we performed linkage studies with chromosome 11 markers in prolactinoma families to determine whether the two genes map to the same location. All proximal chromosome 11q markers gave positive LOD scores, and no recombinants were seen with PYGM (LOD score 15.25, recombination fraction .0). All affected individuals from Newfoundland shared the same PYGM allele, providing evidence for a founder effect. The disease in the Pacific Northwest kindred cosegregated with a different PYGM allele.     

Localization of a recessive gene for North American Indian childhood cirrhosis to chromosome region 16q22-and identification of a shared haplotype

North American Indian childhood cirrhosis (NAIC, or CIRH1A) is an isolated nonsyndromic form of familial cholestasis reported in Ojibway-Cree children and young adults in northwestern Quebec. The pattern of transmission is consistent with an autosomal recessive mode of inheritance. To map the NAIC locus, we performed a genomewide scan on three DNA pools of samples from 13 patients, 16 unaffected siblings, and 22 parents from five families. Analysis of 333 highly polymorphic markers revealed 3 markers with apparent excess allele sharing among affected individuals. Additional mapping identified a chromosome 16q segment shared by all affected individuals. When the program FASTLINK/LINKAGE was used and a completely penetrant autosomal recessive mode of inheritance was assumed, a maximum LOD score of 4.44 was observed for a recombination fraction of 0, with marker D16S3067. A five-marker haplotype (D16S3067, D16S752, D16S2624, D16S3025, and D16S3106) spanning 4.9 cM was shared by all patients. These results provide significant evidence of linkage for a candidate gene on chromosome 16q22.     

Mapping of the variegate porphyria (VP) gene: Contradictory evidence for linkage between VP and microsatellite markers at chromosome 14q32

The gene for variegate porphyria (VP), an autosomal dominant disease with a high prevalence in South Africa, evidently due to a founder effect, was previously mapped to chromosome 14q32. In the current study this localization was evaluated by linkage and haplotype analyses using microsatellite markers spanning a region of more than 20 cM on chromosome 14q32. In many recent studies linkage disequilibrium between disease and marker loci has been utilized to map genes in founder populations, but we could not find any association between VP and the markers used in this study. Our data suggest that the allocation of VP to chromosome 14q32 may be incorrect. Received: 1 September 1995 / Revised: 1 November 1995     

Prostate cancer susceptibility Loci identified on chromosome 12 in African Americans

Prostate cancer (PCa) is a complex disease that disproportionately affects African Americans and other individuals of African descent. A number of regions across the genome have been associated to PCa, most of them with moderate effects. A few studies have reported chromosomal changes on 12p and 12q that occur during the onset and development of PCa but to date no consistent association of the disease with chromosome 12 polymorphic variation has been identified. In order to unravel genetic risk factors that underlie PCa health disparities we investigated chromosome 12 using ancestry informative markers (AIMs), which allow us to distinguish genomic regions of European or West African origin, and tested them for association with PCa. Additional SNPs were genotyped in those areas where significant signals of association were detected. The strongest signal was discovered at the SNP rs12827748, located upstream of the PAWR gene, a tumor suppressor, which is amply expressed in the prostate. The most frequent allele in Europeans was the risk allele among African Americans. We also examined vitamin D related genes, VDR and CYP27B1, and found a significant association of PCa with the TaqI polymorphism (rs731236) in the former. Although our results warrant further investigation we have uncovered a genetic susceptibility factor for PCa in a likely candidate by means of an approach that takes advantage of the differential contribution of parental groups to an admixed population.     

Genome-wide scan identifies loci associated with classical BSE occurrence

Classical bovine spongiform encephalopathy (BSE) is an acquired prion disease that is invariably fatal in cattle and has been implicated as a significant human health risk. Sequence variations in the coding region of the prion gene (PRNP) have been associated with acquired transmissible spongiform encephalopathy (TSE) susceptibility in mammals; however, this is not the case in cattle. It has been hypothesized that genes, in addition to the prion gene, contribute to genetic susceptibility of acquired TSEs. Accordingly, genetic studies of classical BSE in cattle identified loci other than PRNP that are associated with disease incidence. The objective of this study was to utilize a genome-wide association study to test for genetic loci associated with classical BSE. The samples include 143 BSE affected (case) and 173 unaffected half sib (control) animals collected in the mid 1990s in Southern England. The data analysis identifies loci on two different chromosomes associated with BSE disease occurrence. Most notable is a single nucleotide polymorphism on chromosome 1 at 29.15 Mb that is associated with BSE disease (p = 3.09E-05). Additionally, a locus on chromosome 14, within a cluster of SNPs showed a trend toward significance (p = 5.24E-05). It is worth noting that in a human vCJD study markers on human chromosome 8, a region with shared synteny to the region identified on cattle chromosome 14, were associated with disease. Further, our candidate genes appear to have plausible biological relevance with the known etiology of TSE disease. One of the candidate genes is hypothetical gene LOC521010, similar to FK506 binding protein 2 located on chromosome 1 at 29.32 Mb. This gene encodes a protein that is a member of the immunophilin protein family and is involved in basic cellular processes including protein folding. The chromosomal regions identified in this study and candidate genes within these regions merit further investigation.     

Evidence for shared ancestral polymorphism rather than recurrent gene flow at microsatellite loci differentiating two hybridizing oaks (Quercus spp.)

Quercus petraea and Quercus robur are two closely related oak species, considered to hybridize. Genetic markers, however, indicate that despite sharing most alleles, the two species remain separate genetic units. Analysis of 20 microsatellite loci in multiple populations from both species suggested a genome-wide differentiation. Thus, the allele sharing between both species could be explained either by low rates of gene flow or shared ancestral variation. We performed further analyses of population differentiation in a biogeographical setting and an admixture analysis in mixed oak stands to distinguish between both hypotheses. Based on our results we propose that the low genetic differentiation among these species results from shared ancestry rather than high rates of gene flow.     

The Human Obesity Gene Map: The 1996 Update

An update of the human obesity gene map up to October 1996 is presented. Evidence from Mendelian disorders exhibiting obesity as a clinical feature, single-gene mutation rodent models, quantitative trait loci uncovered in crossbreeding experiments with mouse, rat, and pig models, association and case-control studies with candidate genes, and linkage studies with genes and other markers is reviewed. All chromosomal locations of the animal loci are converted into human genome locations based on syntenic relationships between the genomes. A complete listing of all these loci reveals that only 4 of the 24 human chromosomes are not yet represented, i.e., 9, 18, 21, and Y. Several chromosome arms are characterized by the presence of several putative loci. The following arms include at least three such loci: 1p, 1q, 3p, 4q, 6p, 7q, 8p, 8q, 11p, 11q, 15q, 20q, and Xq. Studies with negative association and linkage results are also reviewed.     

Linkage of Bipolar Affective Disorder to Chromosome 18 Markers in a New Pedigree Series

Several groups have reported evidence suggesting linkage of bipolar affective disorder (BPAD) to chromosome 18. We have reported data from 28 pedigrees that showed linkage to marker loci on 18p and to loci 40 cM distant on 18q. Most of the linkage evidence derived from families with affected phenotypes in only the paternal lineage and from marker alleles transmitted on the paternal chromosome. We now report results from a series of 30 new pedigrees (259 individuals) genotyped for 13 polymorphic markers spanning chromosome 18. Subjects were interviewed by a psychiatrist and were diagnosed by highly reliable methods. Genotypes were generated with automated technology and were scored blind to phenotype. Affected sib pairs showed excess allele sharing at the 18q markers D18S541 and D18S38. A parent-of-origin effect was observed, but it was not consistently paternal. No robust evidence of linkage was detected for markers elsewhere on chromosome 18. Multipoint nonparametric linkage analysis in the new sample combined with the original sample of families supports linkage on chromosome 18q, but the susceptibility gene is not well localized.     

Evidence for a novel glaucoma locus at chromosome 3p21-22

Primary open-angle glaucoma (POAG) is one of the leading causes of blindness in the world. It is a clinically variable group of diseases with the majority of cases presenting as the late onset adult type. Several chromosomal loci have been implicated in disease aetiology, but causal mutations have only been identified in a small proportion of glaucoma. We have previously described a large six-generation Tasmanian family with POAG exhibiting genetic heterogeneity. In this family, approximately one third of affected individuals presented with a glutamine-368-STOP (Q368STOP) mutation in the myocilin gene. We now use a Markov Chain Monte Carlo (MCMC) method to identify a second disease region in this family on the short arm of chromosome 3. This disease locus was initially mapped to the marker D3S1298 and a subsequent minimum disease region of 9 cM between markers D3S1298 and D3S1289 was identified through additional mapping. The region did not overlap with any previously described locus for POAG. Using a multiplicative relative risk model, we identified a positive association between this region and the Q368STOP mutation of myocilin on chromosome 1 in affected individuals. These findings provide evidence of a new autosomal dominant glaucoma locus on the short arm of chromosome 3.     

Confirmation of linkage of Van der Woude syndrome to chromosome 1q32: evidence of association with STR alleles suggests possible unique origin of the disease mutation

Van der Woude syndrome (VWS) is an autosomal dominant craniofacial disorder with high penetrance and variable expression. Its clinical features are variably expressed, but include cleft lip and/or cleft palate, lip pits and hypodontia. All VWS families studied to date map the disease gene to a < 2 cM region of chromosome 1q32, with no evidence of locus heterogeneity. The aim of this study is to refine the localization of the VWS gene and to further assess possible heterogeneity. We analyzed four multiplex VWS families. All available members were clinically assessed and genotyped for 19 short tandem repeat markers on chromosome 1 in the VWS candidate gene region. We performed two-point and multipoint limit of detection (LOD) score analyses using a high penetrance autosomal dominant model. All families showed positive LOD scores without any recombination in the candidate region. The largest two-point LOD score was 5.87. Our assay method for short tandem repeat (STR) markers provided highly accurate size estimation of marker allele fragment sizes, and therefore enabled us to determine the specific alleles segregating with the VWS gene in each of our four families. We observed a striking pattern of STR allele sharing at several closely linked loci among our four Caucasian VWS families recruited at three different locations in the US. These results suggest the possibility of a unique origin for a mutation responsible for many or most cases of VWS.     

Fine localization of a major disease-susceptibility locus for diffuse panbronchiolitis

Diffuse panbronchiolitis affecting East Asians is strongly associated with the class I human leukocyte antigen (HLA) alleles. Recent observations suggest that a major disease-susceptibility gene may be located between the HLA-B and HLA-A loci in the class I region of the major histocompatibility complex on chromosome 6. To test this possibility, we analyzed 14 polymorphic markers in 92 Japanese patients and 93 healthy controls. Of these, seven marker alleles, including HLA-B54 and HLA-A11, were significantly associated with the disease. Maximum-likelihood haplotype analysis and subsequent direct determination of individual haplotypes identified a group of disease-associated haplotypes, one of which contained all seven disease-associated marker alleles. Another haplotype, containing HLA-B*5504, was also associated with the disease. All these haplotypes seem to have diverged from a common ancestral haplotype in East Asians and share a specific segment containing three consecutive markers between the S and TFIIH loci in the class I region. Furthermore, one of the markers within the candidate region showed the highest delta value, indicating the strongest association. Of 20 Korean patients with diffuse panbronchiolitis, 17 also shared the combination of the disease-associated marker alleles within the candidate region. These results indicate that an HLA-associated major susceptibility gene for diffuse panbronchiolitis is probably located within the 200 kb in the class I region 300 kb telomeric of the HLA-B locus on the chromosome 6p21.3.     

Biogeography in a continental island: population structure of the relict endemic centipede Craterostigmus tasmanianus (Chilopoda, Craterostigmomorpha) in Tasmania using 16S rRNA and COI

We used 16S ribosomal RNA (rRNA) and cytochrome c oxidase subunit I (COI) sequence data to investigate the population structure in the centipede Craterostigmus tasmanianus Pocock, 1902 (Chilopoda: Craterostigmomorpha: Craterostigmidae) and to look for possible barriers to gene flow on the island of Tasmania, where C. tasmanianus is a widespread endemic. We first confirmed a molecular diagnostic character in 28S rRNA separating Tasmanian Craterostigmus from its sister species Craterostigmus crabilli (Edgecombe and Giribet 2008) in New Zealand and found no shared polymorphism in this marker for the 2 species. In Tasmania, analysis of molecular variance analysis showed little variation at the 16S rRNA and COI loci within populations (6% and 13%, respectively), but substantial variation (56% and 48%, respectively) among populations divided geographically into groups. We found no clear evidence of isolation by distance using a Mantel test. Bayesian clustering and gene network analysis both group the C. tasmanianus populations in patterns which are broadly concordant with previously known biogeographical divisions within Tasmania, but we did not find that genetic distance varied in a simple way across cluster boundaries. The coarse-scale geographical sampling on which this study was based should be followed in the future by sampling at a finer spatial scale and to investigate genetic structure within clusters and across cluster boundaries.     

A genomewide screen for autism: strong evidence for linkage to chromosomes 2q, 7q, and 16p

Autism is characterized by impairments in reciprocal communication and social interaction and by repetitive and stereotyped patterns of activities and interests. Evidence for a strong underlying genetic predisposition comes from twin and family studies, although susceptibility genes have not yet been identified. A whole-genome screen for linkage, using 83 sib pairs with autism, has been completed, and 119 markers have been genotyped in 13 candidate regions in a further 69 sib pairs. The addition of new families and markers provides further support for previous reports of linkages on chromosomes 7q and 16p. Two new regions of linkage have also been identified on chromosomes 2q and 17q. The most significant finding was a multipoint maximum LOD score (MLS) of 3.74 at marker D2S2188 on chromosome 2; this MLS increased to 4.80 when only sib pairs fulfilling strict diagnostic criteria were included. The susceptibility region on chromosome 7 was the next most significant, generating a multipoint MLS of 3.20 at marker D7S477. Chromosome 16 generated a multipoint MLS of 2.93 at D16S3102, whereas chromosome 17 generated a multipoint MLS of 2.34 at HTTINT2. With the addition of new families, there was no increased allele sharing at a number of other loci originally showing some evidence of linkage. These results support the continuing collection of multiplex sib-pair families to identify autism-susceptibility genes.     

汉族群体23个Y染色体双等位基因标记遗传多态性研究

为了筛选在汉族群体中具有多态性的Y染色体双等位基因标记并获取其群体遗传学数据。采用片段长度差异等位基因特异性PCR和PAGE技术对武汉地区160名男性汉族无关个体的23个Y染色体双等位基因标记（M7,M9,M50,M88,M89,M95,M111,M117,M119,M121,M122,M134,M159,M164,M175,M214,LINE1,MSY2,RPS4Y711,SRY＋465,IMS-JST164520,IMS-JST021354和IMS-JST003305）进行分型。除M50、M159和M164外,其余20个标记在武汉汉族群体中均具有遗传多态性,其基因多样性（GD）范围为0.0126～0.4855,共检出35种不同单体群组合（Hg1～35）,单体群多样性（HD）为0.9471。表明20个Y染色体双等位基因标记组成的单体群具有较高的遗传多样性,在法医学应用和群体进化研究中具有较高的实用价值。     

Random search for shared chromosomal regions in four affected individuals: the assignment of a new hereditary ataxia locus.

Infantile-onset spinocerebellar ataxia (IOSCA) is an autosomal recessively inherited progressive neurological disorder of unknown etiology. This ataxia, identified so far only in the genetically isolated Finnish population, does not share gene locus with any of the previously identified hereditary ataxias, and a random mapping approach was adopted to assign the IOSCA locus. Based on the assumption of one founder mutation, a primary screening of the genome was performed using samples from just four affected individuals in two consanguineous pedigrees. The identification of a shared chromosomal region in these four patients provided the first evidence that the IOSCA gene locus is on chromosome 10q23.3-q24.1, which was confirmed by conventional linkage analysis in the complete family material. Strong linkage disequilibrium observed between IOSCA and the linked markers was utilized to define accurately the critical chromosomal region. The results showed the power of linkage disequilibrium in the locus assignment of diseases with very limited family materials.