Positional cloning of disease genes: advantages of genetic isolates

Whole-exome sequencing of individuals from an isolated population implicates rare risk variants in bipolar disorder

Bipolar disorder affects about 1% of the world's population, and its estimated heritability is about 75%. Only few whole genome or whole-exome sequencing studies in bipolar disorder have been reported, and no rare coding variants have yet been robustly identified. The use of isolated populations might help finding variants with a recent origin, more likely to have drifted to higher frequency by chance. Following this approach, we investigated 28 bipolar cases and 214 controls from the Faroe Islands by whole exome sequencing, and the results were followed-up in a British sample of 2025 cases and 1358 controls. Seventeen variants in 16 genes in the single-variant analysis, and 3 genes in the gene-based statistics surpassed exome-wide significance in the discovery phase. The discovery findings were supported by enrichment analysis of common variants from genome-wide association studies (GWAS) data and interrogation of protein-protein interaction networks. The replication in the British sample confirmed the association with NOS1 (missense variant rs79487279) and NCL (gene-based test). A number of variants from the discovery set were not present in the replication sample, including a novel PITPNM2 missense variant, which is located in a highly significant schizophrenia GWAS locus. Likewise, PIK3C2A identified in the gene-based analysis is located in a combined bipolar and schizophrenia GWAS locus. Our results show support both for existing findings in the literature, as well as for new risk genes, and identify rare variants that might provide additional information on the underlying biology of bipolar disorder.

PRIMAL: Fast and accurate pedigree-based imputation from sequence data in a founder population

Founder populations and large pedigrees offer many well-known advantages for genetic mapping studies, including cost-efficient study designs. Here, we describe PRIMAL (PedigRee IMputation ALgorithm), a fast and accurate pedigree-based phasing and imputation algorithm for founder populations. PRIMAL incorporates both existing and original ideas, such as a novel indexing strategy of Identity-By-Descent (IBD) segments based on clique graphs. We were able to impute the genomes of 1,317 South Dakota Hutterites, who had genome-wide genotypes for ~300,000 common single nucleotide variants (SNVs), from 98 whole genome sequences. Using a combination of pedigree-based and LD-based imputation, we were able to assign 87% of genotypes with >99% accuracy over the full range of allele frequencies. Using the IBD cliques we were also able to infer the parental origin of 83% of alleles, and genotypes of deceased recent ancestors for whom no genotype information was available. This imputed data set will enable us to better study the relative contribution of rare and common variants on human phenotypes, as well as parental origin effect of disease risk alleles in >1,000 individuals at minimal cost.

Systems biology of the vervet monkey

Nonhuman primates (NHP) provide crucial biomedical model systems intermediate between rodents and humans. The vervet monkey (also called the African green monkey) is a widely used NHP model that has unique value for genetic and genomic investigations of traits relevant to human diseases. This article describes the phylogeny and population history of the vervet monkey and summarizes the use of both captive and wild vervet monkeys in biomedical research. It also discusses the effort of an international collaboration to develop the vervet monkey as the most comprehensively phenotypically and genomically characterized NHP, a process that will enable the scientific community to employ this model for systems biology investigations.

Insight into genetic determinants of resting heart rate

BACKGROUND

Recent studies suggested that resting heart rate (RHR) might be an independent predictor of cardiovascular mortality and morbidity. Nonetheless, the interrelation between RHR and cardiovascular diseases is not clear. In order to resolve this puzzle, the importance of genetic determinants of RHR has been recently suggested, but it needs to be further investigated.

OBJECTIVE

The aim of this study was to estimate the contribution of common genetic variations on RHR using Genome Wide Association Study.

METHODS

We performed a Genome Wide Association Study in an isolated population cohort of 1737 individuals, the Italian Network on Genetic Isolates - Friuli Venezia Giulia (INGI-FVG). Moreover, a haplotype analysis was performed. A regression tree analysis was run to highlight the effect of each haplotype combination on the phenotype.

RESULTS

A significant level of association (p<5 × 10(-8)) was detected for Single Nucleotide Polymorphisms (SNPs) in two genes expressed in the heart: MAML1 and CANX. Founding that the three different variants of the haplotype, which encompass both genes, yielded a phenotypic correlation. Indeed, a haplotype in homozygosity is significantly associated with the lower quartile of RHR (RHR ≤ 58 bpm). Moreover no significant association was found between cardiovascular risk factors and the different haplotype combinations.

CONCLUSION

Mastermind-like 1 and Calnexin were found to be associated with RHR. We demonstrated a relation between a haplotype and the lower quartile of RHR in our populations. Our findings highlight that genetic determinants of RHR may be implicated in determining cardiovascular diseases and could allow a better risk stratification.

Assessing the genetic landscape of a contact zone: the case of European hare in northeastern Greece

The European hare populations of the Balkan Peninsula comprise two divergent phylogenetic lineages with discrete geographical distribution slightly overlapping in the area of northeastern Greece and Bulgaria. Here we elucidate their contact zone, by defining the spatial distributional pattern of the two highly divergent groups, detecting individuals of hybrid origin, and identifying genetic barriers present in the area of their co-existence. Specimens from northeastern Greece were assayed for lineage assignment and population genetic inference based on a 511 bp fragment of mitochondrial DNA control region and allelic data from 10 microsatellite loci. Bayesian analyses on original and simulated genotypes were performed allowing for the contact zone delineation. Our results indicate high genetic diversity in both nuclear and mitochondrial DNA, strong population structure and non random spatial distribution of the differentiated gene pools. The information provided by the two types of molecular markers yielded consistent results. This study comprises a fine scale analysis of the contact zone between the two evolutionary lineages of European brown hares in northeastern Greece. Specific questions on the spatial patterns where addressed for the first time. Furthermore, hypotheses regarding the presence of hybrids were also tested. As a result, interpretive power to the diversity patterns observed today in the Balkans was added and previously overlooked aspects of the species biology were highlighted.

Custom genotyping for substance addiction susceptibility genes in Jordanians of Arab descent

Background

Both environmental and genetic factors contribute to individual susceptibility to initiation of substance use and vulnerability to addiction. Determining genetic risk factors can make an important contribution to understanding the processes leading to addiction. In order to identify gene(s) and mechanisms associated with substance addiction, a custom platform array search for a genetic association in a case/control of homogenous Jordanian Arab population was undertaken. Patients meeting the DSM-VI criteria for substance dependence (n = 220) and entering eight week treatment program at two Jordanian Drug Rehabilitation Centres were genotyped. In addition, 240 healthy controls were also genotyped. The sequenom MassARRAY system (iPLEX GOLD) was used to genotype 49 single nucleotide polymorphisms (SNPs) within 8 genes (DRD1, DRD2, DRD3, DRD4, DRD5, BDNF, SLC6A3 and COMT).

Results

This study revealed six new associations involving SNPs within DRD2 gene on chromosome 11. These six SNPs within the DRD2 were found to be most strongly associated with substance addiction in the Jordanian Arabic sample. The strongest statistical evidence for these new association signals were from rs1799732 in the C/−C promoter and rs1125394 in A/G intron 1 regions of DRD2, with the overall estimate of effects returning an odds ratio of 3.37 (χ2 (2, N = 460) = 21, p-value = 0.000026) and 1.78 (χ2 (2, N = 460) = 8, p-value = 0.001), respectively. It has been suggested that DRD2, dopamine receptor D2, plays an important role in dopamine secretion and the signal pathways of dopaminergic reward and drug addiction.

Conclusion

This study is the first to show a genetic link to substance addiction in a Jordanian population of Arab descent. These findings may contribute to our understanding of drug addiction mechanisms in Middle Eastern populations and how to manage or dictate therapy for individuals. Comparative analysis with different ethnic groups could assist further improving our understanding of these mechanisms.

Copy number variation across European populations

Genome analysis provides a powerful approach to test for evidence of genetic variation within and between geographical regions and local populations. Copy number variants which comprise insertions, deletions and duplications of genomic sequence provide one such convenient and informative source. Here, we investigate copy number variants from genome wide scans of single nucleotide polymorphisms in three European population isolates, the island of Vis in Croatia, the islands of Orkney in Scotland and the South Tyrol in Italy. We show that whereas the overall copy number variant frequencies are similar between populations, their distribution is highly specific to the population of origin, a finding which is supported by evidence for increased kinship correlation for specific copy number variants within populations.

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.

Familial transmission of schizophrenia in Palau: A 20-year genetic epidemiological study in three generations

Our genetic epidemiological studies of schizophrenia and other psychotic disorders (SCZ) in the isolated population of Palau have been ongoing for 20 years. Results from the first decade showed that Palau has an elevated prevalence of SCZ and that cases cluster in extended multigenerational pedigrees interconnected via complex genetic relationships after centuries of endogamous, but not consanguineous, marriages. The aim of our second decade of research, which extended data collection into a third generation of young, high-risk (HR) Palauans, was to identify significant predictors of intergenerational transmission of illness. Our findings revealed that degree of familial loading and gender effects on reproductive fitness are important modifiers of risk for transmission of SCZ. Among 45 distinct multiplex families, we identified 10 high-density (HD) Palauan families, each with 7-29 SCZ cases, which contain half of Palau's 260 SCZ cases and 80% of the 113 SCZ cases with one or more affected first-degree relatives, indicating that familial loading is a major risk factor for SCZ in Palau. Cases that belong to multiply affected sibships are more common than cases with an affected parent. Furthermore, only 6/38 multiply affected sibships have an affected parent, strong evidence that many unaffected parents are obligate carriers of susceptibility genes. Although reproductive fitness is dramatically reduced in affected males, the 30% minority who do become fathers are twice as likely as affected mothers to transmit SCZ to an offspring. As they evolve, these HD families can help to elucidate the genetic mechanisms that predict intergenerational transmission of SCZ.

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.

Case-control association testing in the presence of unknown relationships

Genome-wide association studies result in inflated false-positive results when unrecognized cryptic relatedness exists. A number of methods have been proposed for testing association between markers and disease with a correction for known pedigree-based relationships. However, in most case-control studies, relationships are generally unknown, yet the design is predicated on the assumption of at least ancestral relatedness among cases. Here, we focus on adjusting cryptic relatedness when the genealogy of the sample is unknown, particularly in the context of samples from isolated populations where cryptic relatedness may be problematic. We estimate cryptic relatedness using maximum-likelihood methods and use a corrected chi(2) test with estimated kinship coefficients for testing in the context of unknown cryptic relatedness. Estimated kinship coefficients characterize precisely the relatedness between truly related people, but are biased for unrelated pairs. The proposed test substantially reduces spurious positive results, producing a uniform null distribution of P-values. Especially with missing pedigree information, estimated kinship coefficients can still be used to correct non-independence among individuals. The corrected test was applied to real data sets from genetic isolates and created a distribution of P-value that was close to uniform. Thus, the proposed test corrects the non-uniform distribution of P-values obtained with the uncorrected test and illustrates the advantage of the approach on real data.

Comparing population structure as inferred from genealogical versus genetic information

Algorithms for inferring population structure from genetic data (ie, population assignment methods) have shown to effectively recognize genetic clusters in human populations. However, their performance in identifying groups of genealogically related individuals, especially in scanty-differentiated populations, has not been tested empirically thus far. For this study, we had access to both genealogical and genetic data from two closely related, isolated villages in southern Italy. We found that nearly all living individuals were included in a single pedigree, with multiple inbreeding loops. Despite F(st) between villages being a low 0.008, genetic clustering analysis identified two clusters roughly corresponding to the two villages. Average kinship between individuals (estimated from genealogies) increased at increasing values of group membership (estimated from the genetic data), showing that the observed genetic clusters represent individuals who are more closely related to each other than to random members of the population. Further, average kinship within clusters and F(st) between clusters increases with increasingly stringent membership threshold requirements. We conclude that a limited number of genetic markers is sufficient to detect structuring, and that the results of genetic analyses faithfully mirror the structuring inferred from detailed analyses of population genealogies, even when F(st) values are low, as in the case of the two villages. We then estimate the impact of observed levels of population structure on association studies using simulated data.

A novel locus for arterial hypertension on chromosome 1p36 maps to a metabolic syndrome trait cluster in the Sorbs, a Slavic population isolate in Germany

OBJECTIVE

Genome-wide linkage studies and genome-wide association studies have not as yet identified major genes contributing to primary hypertension in the general population. This state-of-affairs suggests considerable heterogeneity with small contributing effects for primary hypertension, or other complex genetic traits, in outbred populations. Isolated populations, as recent data from Iceland and French Canada suggest, could offer a solution to this problem.

METHODS

We studied a Slavic isolate in Germany, the Sorbs, and genotyped 1040 polymorphic microsatellite markers in 87 multigeneration families.

RESULTS

Our genome-wide linkage scan revealed a locus on chromosome 1p36.13 at D1S3669-D1S2826 (40.95 cM Marshfield coordinates; logarithm of the odds = 3.45, nominal P = 0.00003) that reached genome-wide significance (P = 0.004), indicating the increased power in isolated populations. The chromosome 1 locus maps to a region in which traits such as diabetes, hyperlipidemia, obesity and BMI cluster.

CONCLUSION

Our results suggest that this locus contributes to the metabolic syndrome, and that further attention in this and other populations is warranted.

A genome-wide association scan of RR and QT interval duration in 3 European genetically isolated populations: the EUROSPAN project

BACKGROUND

We set out to identify common genetic determinants of the length of the RR and QT intervals in 2325 individuals from isolated European populations.

METHODS AND RESULTS

We analyzed the heart rate at rest, measured as the RR interval, and the length of the corrected QT interval for association with 318 237 single-nucleotide polymorphisms. The RR interval was associated with common variants within GPR133, a G-protein-coupled receptor (rs885389, P=3.9 x 10(-8)). The QT interval was associated with the earlier reported NOS1AP gene (rs2880058, P=2.00 x 10(-10)) and with a region on chromosome 13 (rs2478333, P=4.34 x 10(-8)), which is 100 kb from the closest known transcript LOC730174 and has previously not been associated with the length of the QT interval.

CONCLUSIONS

Our results suggested an association between the RR interval and GPR133 and confirmed an association between the QT interval and NOS1AP.

Bipolar disorder in the Bulgarian Gypsies: genetic heterogeneity in a young founder population

We report the results of follow-up analyses of 12 genomic regions showing evidence of linkage in a genome-wide scan (GWS) of Gypsy families with bipolar affective disorder (BPAD). The Gypsies are a young founder population comprising multiple genetically differentiated sub-isolates with strong founder effect and limited genetic diversity. The BPAD families belong to a single sub-isolate and are connected by numerous inter-marriages, resulting in a super-pedigree with 181 members. We aimed to re-assess the positive GWS findings and search for evidence of a founder susceptibility allele after the addition of newly recruited subjects, some changes in diagnostic assignment, and the use of denser genetic maps. Linkage analysis was conducted with SimWalk2, accommodating the full complexity of pedigree structure and using a conservative narrow phenotype definition (BPAD only). Six regions were rejected, while 1p36, 13q31, 17p11, 17q21, 6q24, and 4q31 produced nominally significant results in both the individual families and the super-pedigree. Haplotypes were reconstructed and joint tests for linkage and association were done for the most promising regions. No common ancestral haplotype was identified by sequencing a strong positional and functional candidate gene (GRM1) and additional STR genotyping in the top GWS region, 6q24. The best supported region was a 12 cM interval on 4q31, also implicated in previous studies, where we obtained significant results in the super-pedigree using both SimWalk2 (P = 0.004) and joint Pseudomarker analysis of linkage and linkage disequilibrium (P = 0.000056). The size of the region and the characteristics of the Gypsy population make it suitable for LD mapping.

Linkage mapping of CVD risk traits in the isolated Norfolk Island population

To understand the underlying genetic architecture of cardiovascular disease (CVD) risk traits, we undertook a genome-wide linkage scan to identify CVD quantitative trait loci (QTLs) in 377 individuals from the Norfolk Island population. The central aim of this research focused on the utilization of a genetically and geographically isolated population of individuals from Norfolk Island for the purposes of variance component linkage analysis to identify QTLs involved in CVD risk traits. Substantial evidence supports the involvement of traits such as systolic and diastolic blood pressures, high-density lipoprotein-cholesterol, low-density lipoprotein-cholesterol, body mass index and triglycerides as important risk factors for CVD pathogenesis. In addition to the environmental influences of poor diet, reduced physical activity, increasing age, cigarette smoking and alcohol consumption, many studies have illustrated a strong involvement of genetic components in the CVD phenotype through family and twin studies. We undertook a genome scan using 400 markers spaced approximately 10 cM in 600 individuals from Norfolk Island. Genotype data was analyzed using the variance components methods of SOLAR. Our results gave a peak LOD score of 2.01 localizing to chromosome 1p36 for systolic blood pressure and replicated previously implicated loci for other CVD relevant QTLs.

Genetic landscape of the people of India: a canvas for disease gene exploration

Analyses of frequency profiles of markers on disease or drug-response related genes in diverse populations are important for the dissection of common diseases. We report the results of analyses of data on 405 SNPs from 75 such genes and a 5.2 Mb chromosome, 22 genomic region in 1871 individuals from diverse 55 endogamous Indian populations. These include 32 large (>10 million individuals) and 23 isolated populations, representing a large fraction of the people of India. We observe high levels of genetic divergence between groups of populations that cluster largely on the basis of ethnicity and language. Indian populations not only overlap with the diversity of HapMap populations, but also contain population groups that are genetically distinct. These data and results are useful for addressing stratification and study design issues in complex traits especially for heterogeneous populations.

Association of IL23R, TNFRSF1A, and HLA-DRB1*0103 allele variants with inflammatory bowel disease phenotypes in the Finnish population

BACKGROUND

Crohn's disease (CD) and ulcerative colitis (UC), 2 major forms of inflammatory bowel disease (IBD), are complex disorders with significant genetic predisposition. The first CD-associated gene, CARD15/NOD2, was recently identified and since then several reports on novel IBD candidate genes have emerged. We investigated disease phenotype association to genetic variations in IL23R, ATG16L1, DLG5, ABCB1/MDR1, TLR4, TNFRSF1A, chromosome 5 risk haplotype including SLC22A4 and SLC22A5, and HLA-DRB1*0103 allele among Finnish IBD patients.

METHODS

A total of 699 IBD patients were genotyped for disease-associated variants by polymerase chain reaction (PCR) and restriction enzyme digestion or Sequenom iPLEX method.

RESULTS

Five markers spanning the IL23R gene were associated with CD. The SNP (single nucleotide polymorphism) rs2201841 gave the strongest association (P = 0.002). The rare HLA-DRB1*0103 allele was found to associate with UC (P = 0.008), and the TNFRSF1A A36G variant was associated with familial UC (P = 0.007). Upon phenotypic analysis we detected association between familial UC and rare TNFRSF1A alleles 36G and IVS6+10G (P = 0.001 and P = 0.042, respectively). In addition, IL23R markers were associated with stricturing CD (P = 0.010-0.017), and ileocolonic CD was more prevalent in the carriers of the same 2 TNFRSF1A variants (P = 0.021 and P = 0.028, respectively). Less significant genotype-phenotype associations were observed for the TLR4 and HLA variants.

CONCLUSIONS

We were able to replicate the association of the IL23R variants with CD as well as HLA-DRB1*0103 with UC; confirmation of TNFRSF1A association with UC needs additional studies. Our findings also suggest that polymorphisms at IL23R and TNFRSF1A, and possibly HLA and TLR4, loci may account for phenotypic variation in IBD.

Are the determinants of vertebral endplate changes and severe disc degeneration in the lumbar spine the same? A magnetic resonance imaging study in middle-aged male workers

BACKGROUND

Modic changes are bone marrow lesions visible in magnetic resonance imaging (MRI), and they are assumed to be associated with symptomatic intervertebral disc disease, especially changes located at L5-S1. Only limited information exists about the determinants of Modic changes. The objective of this study was to evaluate the determinants of vertebral endplate (Modic) changes, and whether they are similar for Modic changes and severe disc degeneration focusing on L5-S1 level.

METHODS

228 middle-aged male workers (159 train engineers and 69 sedentary factory workers) from northern Finland underwent sagittal T1- and T2-weighted MRI. Modic changes and disc degeneration were analyzed from the scans. The participants responded to a questionnaire including items of occupational history and lifestyle factors. Logistic regression analysis was used to evaluate the associations between selected determinants (age, lifetime exercise, weight-related factors, fat percentage, smoking, alcohol use, lifetime whole-body vibration) and Modic type I and II changes, and severe disc degeneration (= grade V on Pfirrmann's classification).

RESULTS

The prevalences of the Modic changes and severe disc degeneration were similar in the occupational groups. Age was significantly associated with all degenerative changes. In the age-adjusted analyses, only weight-related determinants (BMI, waist circumference) were associated with type II changes. Exposure to whole-body vibration, besides age, was the only significant determinant for severe disc degeneration. In the multivariate model, BMI was associated with type II changes at L5-S1 (OR 2.75 per one SD = 3 unit increment in BMI), and vibration exposure with severe disc degeneration at L5-S1 (OR 1.08 per one SD = 11-year increment in vibration exposure).

CONCLUSION

Besides age, weight-related factors seem important in the pathogenesis of Modic changes, whereas whole-body vibration was the only significant determinant of severe disc degeneration.

Identification of a chromosome 8p locus for early-onset coronary heart disease in a French Canadian population

Susceptibility to coronary heart disease (CHD) has long been known to exhibit familial aggregation, with heritability estimated to be greater than 50%. The French Canadian population of the Saguenay-Lac Saint-Jean region of Quebec, Canada is descended from a founder population that settled this region 300-400 years ago and this may provide increased power to detect genes contributing to complex traits such as CHD. Probands with early-onset CHD, defined by angiographically determined coronary stenosis, and their relatives were recruited from this population (average sibship size of 6.4). Linkage analysis was performed following a genome-wide microsatellite marker scan on 42 families with 284 individuals. Nonparametric linkage (NPL) analysis provided suggestive evidence for a CHD susceptibility locus on chromosome 8 with an NPL score of 3.14 (P=0.001) at D8S1106. Linkage to this locus was verified by fine mapping in an enlarged sample of 50 families with 320 individuals. This analysis provided evidence of linkage at D8S552 (NPL score=3.53, P=0.0003), a marker that maps to the same location as D8S1106. Candidate genes in this region, including macrophage scavenger receptor 1, farnesyl-diphosphate farnesyltransferase 1, fibrinogen-like 1, and GATA-binding protein 4, were resequenced in all coding exons in both affected and unaffected individuals. Association studies with variants in these and five other genes did not identify a disease-associated mutation. In conclusion, a genome-wide scan and additional fine mapping provide evidence for a locus on chromosome 8 that contributes to CHD in a French Canadian population.

The genetics of IgA nephropathy

IgA nephropathy is the most common form of primary glomerulonephritis. Variations in clinical manifestations indicate that a diagnosis of IgA nephropathy encompasses multiple disease subsets that cannot be distinguished on the basis of renal pathology or clinical variables alone. Familial forms of the disease have been reported throughout the world, but are probably under-recognized because associated urinary abnormalities are often intermittent in affected family members. IgA nephropathy has complex determination, with different genes probably causing disease in different patient subgroups. Of the many pathogenic mechanisms reported, defects in IgA1 glycosylation that lead to formation of immune complexes have been consistently implicated. Here, we present the evidence for genetic contributions to the disease, review clinical patterns of familial disease, and summarize some of the most promising genetic studies conducted to date. Linkage-based approaches to the study of familial forms of the disease have identified significant or suggestive loci on chromosomes 6q22-23, 2q36, 4q26-31, 17q12-22 and 3p24-23, but no causal gene has yet been identified. Many interesting, but poorly replicated, genetic association studies have also been reported. We discuss recent developments in analytic tools that should enable genetic studies of sporadic forms of disease by the genome-wide association approach.

Detecting population structure using STRUCTURE software: effect of background linkage disequilibrium

STRUCTURE is the most widely used clustering software to detect population genetic structure. The last version of this software (STRUCTURE 2.1) has been enhanced recently to take into account the occurrence of linkage disequilibrium (LD) caused by admixture between populations. This last version, however, still does not consider the effects of strong background LD caused by genetic drift, and which may cause spurious results. STRUCTURE authors have, therefore, suggested a rough threshold value of the distance (1.0 cM) between two loci below which the pair of loci should not be used. Because of the sensitiveness of LD to demographic events, the distance between loci is not always a good indicator of the strength of LD. In this study, we examine the link between genomic distance and the strength of the correlation between loci (r(LD)) in a free-ranging population of mouflon (Ovis aries), and we present an empirical test of effect of r(LD) on the clustering results provided by the linkage model in STRUCTURE. We showed that a high r(LD) value increases the probability of detecting spurious clustering. We propose to use r(LD) as an index to base a decision on whether or not to use a pair of loci in a clustering analysis.

Linkage to chromosome 11p12 in two Maltese families with a highly penetrant form of osteoporosis

Osteoporosis is a metabolic bone disease with a strong genetic component. Family-based linkage studies were performed by a number of investigators to try to identify loci that might contain genes responsible for an increased susceptibility to osteoporosis. A whole-genome linkage scan using 400 microsatellite markers was performed in 27 members from two Maltese families with a highly penetrant form of osteoporosis. The phenotype was defined by lumbar and femoral z-scores calculated after measurement of bone mineral density by DEXA. Both males and females were among the affected individuals. Multipoint parametric and non-parametric linkage analyses were performed by EasyLinkage v4.01 using GENEHUNTER v2.1, assuming dominant and recessive modes of inheritance with variable penetrance. Evidence of linkage was observed to a marker at 11p12 where a non-parametric LOD score of 5.77 (P=0.0006) was obtained. A maximum heterogeneity LOD score of 2.55 for this region was obtained for the dominant mode of inheritance with 90% penetrance and a phenocopy rate of 1%. Following fine mapping, the critical interval was narrowed to a region that is 52.94 cM from 11p-telomere. In this region, the gene for tumour necrosis factor receptor-associated factor 6 (TRAF6) is located approximately 1 cM away from the indicated marker. Sequencing of the promoter region and exons of the TRAF6 gene revealed three sequence variants, one of which was found in three affected members within one family.

Power and precision of alternate methods for linkage disequilibrium mapping of quantitative trait loci

Linkage disequilibrium (LD) analysis in outbred populations uses historical recombinations to detect and fine map quantitative trait loci (QTL). Our objective was to evaluate the effect of various factors on power and precision of QTL detection and to compare LD mapping methods on the basis of regression and identity by descent (IBD) in populations of limited effective population size (N(e)). An 11-cM region with 6-38 segregating single-nucleotide polymorphisms (SNPs) and a central QTL was simulated. After 100 generations of random mating with N(e) of 50, 100, or 200, SNP genotypes and phenotypes were generated on 200, 500, or 1000 individuals with the QTL explaining 2 or 5% of phenotypic variance. To detect and map the QTL, phenotypes were regressed on genotypes or (assumed known) haplotypes, in comparison with the IBD method. Power and precision to detect QTL increased with sample size, marker density, and QTL effect. Power decreased with N(e), but precision was affected little by N(e). Single-marker regression had similar or greater power and precision than other regression models, and was comparable to the IBD method. Thus, for rapid initial screening of samples of adequate size in populations in which drift is the primary force that has created LD, QTL can be detected and mapped by regression on SNP genotypes without recovering haplotypes.

Recurrence risk to offspring in extended multiplex schizophrenia pedigrees from a Pacific Island isolate

Genetic transmission plays a major role in the pathogenesis of schizophrenia. Family, twin, and adoption studies have consistently shown that risks in relatives are many times greater than the general population risk of approximately 1%. McGue, Gottesman, and Rao (1983; Am J Hum Genet 35:1161-1178) calculated risk estimates of 12.8% for offspring and 3.5% for nieces/nephews of schizophrenia patients based on a large data set of Western European families. The present study evaluated corresponding risk levels in Palau, an isolated population in Micronesia where the prevalence of narrowly (broadly) defined schizophrenia is 1.99% (2.67%) and cases cluster in extended pedigrees, 20 of which contain 80% of affected individuals. We hypothesized that offspring in these extended families would have a higher risk for schizophrenia than offspring in smaller schizophrenia pedigrees from more genetically heterogeneous populations. RDC diagnostic data based on complete ascertainment of cases and their families covering the past two generations were used to quantify empirical recurrence risks in the offspring and nieces/nephews of Palauan schizophrenia patients. Risks to 1st- and 2nd-degree offspring were approximately double the rates found in the smaller Western European families: 23.4% in the offspring of an affected parent, 6.4% in offspring with one affected aunt/uncle, and 15.0% in offspring with two or more affected aunts/uncles. Recurrence rates in offspring of an affected parent were 1.6 times higher in males (27.9%) than in females (17.7%). The high risk levels we found in Palauan offspring reflect the elevated population prevalence, strong familial aggregation, and multi-lineal transmission pattern of schizophrenia in Palau.

Association mapping of complex trait loci with context-dependent effects and unknown context variable

A novel method for Bayesian analysis of genetic heterogeneity and multilocus association in random population samples is presented. The method is valid for quantitative and binary traits as well as for multiallelic markers. In the method, individuals are stochastically assigned into two etiological groups that can have both their own, and possibly different, subsets of trait-associated (disease-predisposing) loci or alleles. The method is favorable especially in situations when etiological models are stratified by the factors that are unknown or went unmeasured, that is, if genetic heterogeneity is due to, for example, unknown genes x environment or genes x gene interactions. Additionally, a heterogeneity structure for the phenotype does not need to follow the structure of the general population; it can have a distinct selection history. The performance of the method is illustrated with simulated example of genes x environment interaction (quantitative trait with loosely linked markers) and compared to the results of single-group analysis in the presence of missing data. Additionally, example analyses with previously analyzed cystic fibrosis and type 2 diabetes data sets (binary traits with closely linked markers) are presented. The implementation (written in WinBUGS) is freely available for research purposes from http://www.rni.helsinki.fi/ approximately mjs/.

Linkage disequilibrium mapping in domestic dog breeds narrows the progressive rod-cone degeneration interval and identifies ancestral disease-transmitting chromosome

Canine progressive rod-cone degeneration (prcd) is a retinal disease previously mapped to a broad, gene-rich centromeric region of canine chromosome 9. As allelic disorders are present in multiple breeds, we used linkage disequilibrium (LD) to narrow the approximately 6.4-Mb interval candidate region. Multiple dog breeds, each representing genetically isolated populations, were typed for SNPs and other polymorphisms identified from BACs. The candidate region was initially localized to a 1.5-Mb zero recombination interval between growth factor receptor-bound protein 2 (GRB2) and SEC14-like 1 (SEC14L). A fine-scale haplotype of the region was developed, which reduced the LD interval to 106 kb and identified a conserved haplotype of 98 polymorphisms present in all prcd-affected chromosomes from 14 different dog breeds. The findings strongly suggest that a common ancestor transmitted the prcd disease allele to many of the modern dog breeds and demonstrate the power of the LD approach in the canine model.

Familial aggregation of primary glomerulonephritis in an Italian population isolate: Valtrompia study

Hereditary factors are suspected to contribute to the pathogenesis of sporadic primary glomerulonephritis, but their contribution is difficult to delineate in the general population. We studied the prevalence of primary glomerulonephritis in an isolated population from the extreme northern Valtrompia valley, Northern Italy. Investigation of medical records, community urinary screening program and molecular characterization of the population's ancestry were performed; genealogies of affected individuals were researched. Forty-three patients with primary glomerulonephritis were identified: 25 had biopsy-proven disease (11 immunoglobulin A (IgA) nephropathy; eight mesangial proliferative glomerulonephritis without IgA deposits; four focal segmental glomerular sclerosis; two membranous nephropathy), and 18 had clinical glomerulonephritis. All 43 patients originated from three mountain villages (Collio, San Colombano, and Bovegno). In contrast, we found only four cases of primary glomerulonephritis in two nearby villages (Pezzaze and Tavernole) that shared similar population histories and lifestyles, demonstrating heterogeneity of risk factors for glomerulonephritis (P=3 x 10(-5)). All 43 affected individuals could be traced back to common ancestors (XVI-XVII centuries), enabling the construction of three large pedigree including three parent-child affected pairs and five affected siblings pairs. Molecular data showed lower genetic diversity and increased inbreeding in the Valtrompia population compared to the control population. Molecular and genealogical evidence of limited set of founders and the absence of shared nephritogenic environmental factors suggest that our patients share a common genetic susceptibility to the development of primary glomerulonephritis. Further molecular study of our families will offer the possibility to shed light on the genetic background underlying these glomerular disorders.

Complex trait mapping in isolated populations: Are specific statistical methods required?

In this paper, we review the statistical methods that can be used in isolated populations to map genes involved in complex diseases. Our intention is to highlight the fact that if the features of population isolates may help in the identification of susceptibility factors for complex traits, the choice and design of methods for statistical analysis in these populations deserve particular care. We show that methods designed for outbred samples are generally not appropriate for isolated populations and could lead to false conclusions.

A major locus for hereditary prostate cancer in Finland: localization by linkage disequilibrium of a haplotype in the HPCX region

BACKGROUND

Prostate cancer (PRCA) is the most common cancer in males in the western world. In Finland PRCA has an age-adjusted incidence of 81.5 per 100,000. We previously reported that in Finland, the late-onset cases in families with "no-male-to-male" (NMM) transmission of PRCA accounted for most of the linkage to the HPCX region (Xq27-28). The aim of the present study was to test for linkage disequilibrium (LD) and haplotype-sharing around marker DXS1205 between cases from hereditary prostate cancer (HPC) families and population controls. The initial allelic association was performed between 108 PRCA cases and 257 population controls genotyped for 23 markers in the Xq26-28 region. This resulted in a highly significant nominal one-sided Fisher's exact P-value of 0.0003 for allele ''180'' of marker DXS1205. Subsequently, a similar level of significance was observed for the same allele for marker DXS1205 (P=0.0002) when comparing 60 NMM cases and 257 controls. These results were still significant after Bonferroni correction for multiple testing. Fine mapping efforts included the genotyping of four additional markers D3S2390, bG82i1.9, bG82i1.1, bG82i1.0 and four single nucleotide polymorphisms (SNPs) to augment the original markers around DXS1205.

RESULTS

Our major finding is that markers extending from ''D3S2390'' to ''bG82i1.0'' flank the critical locus, about 150 kb. Levin and Bertell's LD measure (delta), a guide to localization of a possible variant, was 0.42 and 0.41 for alleles of markers bG82i1.9 and DXS1205, respectively.

CONCLUSIONS

In this study, the most significant haplotype comprised the three tightly linked, contiguous markers: ''cen-bG82i1.9-SNP-Hap B-bG82i1.1-tel'' [''197-2-234''] among several possible haplotypes (nominal Fisher's one-sided P=0.003). The two transcription units mapping within this interval are the LDOC1 and SPANXC genes. Positional cloning of the HPCX gene(s) is being facilitated by this exploration of the Xq26-28 region. This study represents the first report identifying a haplotype in the Xq27-28 region for an association between HPCX and X-linked prostate cancer with no-male-to-male transmission in the Finnish population.

Replication in genetic studies of complex traits

Disappointments in replicating initial findings in gene mapping for complex traits are often attributed to small sample sizes and inadequate techniques to determine the threshold value. This is clearly not the whole truth. More fundamental reasons lie in the inherent heterogeneity related to disease, including genetic heterogeneity, differences in allele frequencies, and context-dependency in genetic architecture. There are also other reasons related to the data collection and analysis. Replication may remain a source of frustration unless more emphasis is put on controlling these sources of heterogeneity between studies.

Mapping genes of complex psychiatric diseases in Daghestan genetic isolates

Genetic isolates, which provide outstanding opportunities for identification of susceptibility genes for complex diseases, can be classified as primary (having an ancient demographic history in a stable environment) or secondary (having a younger demographic history) Neel [1992: Minority populations: Genetics, demography, and health, pp. 1-13]. Daghestan contains 26 out of 50 indigenous Caucasus ethnicities that have been in existence for hundreds of generations in the same highland region. The ethnic groups are subdivided into numerous primary isolates. The founder effect and gene drift in these primary isolates may have caused aggregation of specific haplotypes with limited numbers of pathogenic alleles and loci in some isolates relative to others. These are expressed as inter-population differences in lifetime prevalence and features of certain complex clinical phenotypes and in patterns of genetic linkage and linkage disequilibrium (LD). Stable highland and ethnic-cultural environments have led to increased penetrance and a reduced number of phenocopies, which typically hamper the identification of any susceptibility genes for complex diseases. Owing to these characteristics of the primary isolates, a comparative linkage study in the primary isolates allows us to define the number of susceptibility genes for any complex disease and to identify the source of variability and non-replication of linkage analysis results. As part of an ongoing study, seven extended schizophrenia and one nonspecific mental retardation kindreds have been ascertained from Daghestan isolates. Lifetime morbid risk for schizophrenia in the isolates varied from 0 to 5%. A genome scan with markers spaced 10 cM apart was carried out on these pedigrees and linkage analysis was performed using descent graph methods, as implemented in Simwalk2. To identify regions containing susceptibility genes within these kindreds, we followed up those regions with non-parametric and parametric linkage analyses, with the choice of genetic model guided by the results obtained in the NPL. While the analyses are ongoing, the most positive findings were made in different isolated pedigrees on chromosomes 17p11, 3q24, and 22q for schizophrenia and on chromosome 12q for nonspecific mental retardation.

Genetic studies of neuropsychiatric disorders in Costa Rica: a model for the use of isolated populations

The importance of genetics in understanding the etiology of mental illness has become increasingly clear in recent years, as more evidence has mounted that almost all neuropsychiatric disorders have a genetic component. It has also become clear, however, that these disorders are etiologically complex, and multiple genetic and environmental factors contribute to their makeup. So far, traditional linkage mapping studies have not definitively identified specific disease genes for neuropsychiatric disorders, although some potential candidates have been identified via these methods (e.g. the dysbindin gene in schizophrenia; Straub et al., 2002; Schwab et al., 2003). For this reason, alternative approaches are being attempted, including studies in genetically isolated populations. Because isolated populations have a high degree of genetic homogeneity, their use may simplify the process of identifying disease genes in disorders where multiple genes may play a role. Several areas of Latin America contain genetically isolated populations that are well suited for the study of neuropsychiatric disorders. Genetic studies of several major psychiatric illnesses, including bipolar disorder, major depression, schizophrenia, Tourette Syndrome, alcohol dependence, attention deficit hyperactivity disorder, and obsessive-compulsive disorder, are currently underway in these regions. In this paper we highlight the studies currently being conducted by our groups in the Central Valley of Costa Rica to illustrate the potential advantages of this population for genetic studies.

Mutation analysis of INSL3 and GREAT/LGR8 genes in familial cryptorchidism

OBJECTIVES

Male mice deficient in insulin-like 3 hormone (Insl3) or its receptor, Great/Lgr8, exhibit cryptorchidism. Recently, sequence analysis of the human INSL3 and GREAT genes identified several allelic variants. These include polymorphisms without apparent functional consequence and a few alleles encoding products with compromised function. However, loss-of-function alleles appear to be rare in human cryptorchidism. Most patients studied to date are presumed to have had sporadic cryptorchidism. We postulated that any genotypic variants predisposing to cryptorchidism would be more prevalent among patients with familial cryptorchidism.

METHODS

We isolated genomic DNA from 13 individuals with personal and family histories of cryptorchidism and used polymerase chain reaction to amplify all exons of both INSL3 and GREAT, as well as INSL3 proximal promoter sequence, including a putative SF-1 transcription factor binding site. We directly sequenced all 20 amplicons and compared them with the wild-type alleles.

RESULTS

We detected two silent substitutions and one missense (A60T) substitution in exon 1 of INSL3 and two silent substitutions in exon 12 and one missense (I604V) substitution in exon 17 of GREAT, all previously described. We found that in vitro the I604V GREAT variant receptor responds to INSL3 stimulation similarly to the wild-type receptor.

CONCLUSIONS

We found polymorphic alleles of INSL3 and GREAT, but no deleterious mutations among individuals with familial cryptorchidism. Thus, mutations in these two genes are responsible only for a small proportion of familial cryptorchidism.

Linkage disequilibrium mapping in the Newfoundland population: a re-evaluation of the refinement of the Bardet?Biedl syndrome 1 critical interval

Genetically isolated populations, such as Newfoundland, have contributed greatly to the identification of disease-causing genes. A linkage disequilibrium (LD) study involving six Newfoundland families predicted a critical interval for Bardet-Biedl syndrome 1 (BBS1) (Young et al. in Am J Hum Genet 65:1680-1687, 1999), but the subsequent identification of BBS1 revealed that it lies outside this region. This suggested that either there is another gene responsible for BBS in these families or the Newfoundland population may not be ideal for LD studies. We screened these six Newfoundland families for mutations in BBS1 and found that affected individuals in five of them were homozygous for the same M390R mutation. There was no evidence for any BBS1 mutation in the affected individual in the sixth family. Therefore, one of the criteria for LD mapping was not met; namely, there should be a single disease-causing allele in the population. Haplotype analysis of unaffected individuals from south-west Newfoundland and English BBS1 patients homozygous for M390R, revealed that a second criterion for LD mapping was violated. The M390R mutation occurred in a common haplotype and both of these chromosomes, the ancestral wild-type and disease-causing haplotypes, were introduced to Newfoundland and spread by a founder effect. Moreover, it was found that disease-associated alleles occurred at relatively high frequencies in normal haplotypes and this probably accounted for the incorrect prediction in the previous LD study. Knowing the amount of genetic variation and its distribution in the Newfoundland population would be useful to maximize its potential for mapping hereditary disorders.

Genome-wide scan in Portuguese Island families implicates multiple loci in bipolar disorder: fine mapping adds support on chromosomes 6 and 11

As part of an extensive study in the Portuguese Island population of families with multiple patients suffering from bipolar disorder and schizophrenia, we performed an initial genome-wide scan of 16 extended families with bipolar disorder that identified three regions on chromosomes 2, 11, and 19 with genome-wide suggestive linkage and several other regions, including chromosome 6q, also approached suggestive levels of significance. Dick et al. [2003: Am J Hum Genet 73:107-114] recently reported in a study of 250 families with bipolar disorder a maxLOD score of 3.61 near marker D6S1021 on chromosome 6q. This study replicates this finding having detected a peak NPL = 2.02 (P = 0.025) with the same marker D6S1021(104.7 Mb). Higher-density mapping provided additional support for loci on chromosome 6 including marker D6S1021 with an NPL = 2.59 (P = 0.0068) and peaking at marker D6S1639 (125 Mb) with an NPL = 3.06 (P = 0.0019). A similar pattern was detected with higher-density mapping of chromosome 11 with an NPL = 3.15 (P = 0.0014) at marker D11S1883 (63.1 Mb). Simulations at the density of our fine mapping data indicate that less than 1 scan out of 10 would find two such scores genome-wide in the same scan by chance. Our findings provide additional support for a susceptibility locus for bipolar disorder on 6q, as well as, suggesting the importance of denser scans. Published 2004 Wiley-Liss, Inc.

Genomewide linkage analysis of bipolar disorder by use of a high-density single-nucleotide-polymorphism (SNP) genotyping assay: a comparison with microsatellite marker assays and finding of significant linkage to chromosome 6q22

We performed a linkage analysis on 25 extended multiplex Portuguese families segregating for bipolar disorder, by use of a high-density single-nucleotide-polymorphism (SNP) genotyping assay, the GeneChip Human Mapping 10K Array (HMA10K). Of these families, 12 were used for a direct comparison of the HMA10K with the traditional 10-cM microsatellite marker set and the more dense 4-cM marker set. This comparative analysis indicated the presence of significant linkage peaks in the SNP assay in chromosomal regions characterized by poor coverage and low information content on the microsatellite assays. The HMA10K provided consistently high information and enhanced coverage throughout these regions. Across the entire genome, the HMA10K had an average information content of 0.842 with 0.21-Mb intermarker spacing. In the 12-family set, the HMA10K-based analysis detected two chromosomal regions with genomewide significant linkage on chromosomes 6q22 and 11p11; both regions had failed to meet this strict threshold with the microsatellite assays. The full 25-family collection further strengthened the findings on chromosome 6q22, achieving genomewide significance with a maximum nonparametric linkage (NPL) score of 4.20 and a maximum LOD score of 3.56 at position 125.8 Mb. In addition to this highly significant finding, several other regions of suggestive linkage have also been identified in the 25-family data set, including two regions on chromosome 2 (57 Mb, NPL = 2.98; 145 Mb, NPL = 3.09), as well as regions on chromosomes 4 (91 Mb, NPL = 2.97), 16 (20 Mb, NPL = 2.89), and 20 (60 Mb, NPL = 2.99). We conclude that at least some of the linkage peaks we have identified may have been largely undetected in previous whole-genome scans for bipolar disorder because of insufficient coverage or information content, particularly on chromosomes 6q22 and 11p11.

Are common disease susceptibility alleles the same in outbred and founder populations?

Founder populations have been the subjects of complex disease studies because of their decreased genetic heterogeneity, increased linkage disequilibrium and more homogeneous environmental exposures. However, it is possible that disease alleles identified in founder populations may not contribute significantly to susceptibility in outbred populations. In this study we examine the Hutterites, a founder population of European descent, for 103 polymorphisms in 66 genes that are candidates for cardiovascular or inflammatory diseases. We compare the frequencies of alleles at these loci in the Hutterites to their frequencies in outbred European-American populations and test for associations with cardiovascular disease-associated phenotypes in the Hutterites. We show that alleles at these loci are found at similar frequencies in the Hutterites and in outbred populations. In addition, we report associations between 39 alleles or haplotypes and cardiovascular disease phenotypes (P<0.05), with five loci remaining significant after adjusting for multiple comparisons. These data indicate that this founder population is informative and offers considerable advantages for genetic studies of common complex diseases.

Genome-wide scan in Portuguese Island families identifies 5q31-5q35 as a susceptibility locus for schizophrenia and psychosis

Schizophrenia is a common psychiatric disorder with a complex genetic etiology. To understand the genetic basis of this syndrome in Portuguese Island populations, we performed a genome-wide scan of 29 families with schizophrenia, which identified a single region on 5q31-5q35 with strong linkage (NPL=3.09, P=0.0012 at D5S820). Empirical simulations set a genome-wide threshold of NPL=3.10 for significant linkage. Additional support for this locus in schizophrenia comes from higher-density mapping and mapping of 11 additional families. The combined set of 40 families had a peak NPL=3.28 (P=0.00066) at markers D5S2112-D5S820. These data and previous linkage findings from other investigators provide strong and consistent evidence for this genomic region as a susceptibility locus for schizophrenia. Exploratory analyses of a novel phenotype, psychosis, in families with schizophrenia and bipolar disorder detected evidence for linkage to the same markers as found in schizophrenia (peak NPL=3.03, P=0.0012 at D5S820), suggesting that this locus may be responsible for the psychotic symptoms observed in both diseases. Molecular Psychiatry (2004) 9, 213-218. doi:10.1038/sj.mp.4001418 Published online 30 December 2003

The SLC6A14 gene shows evidence of association with obesity

In our previous genome-wide scan of Finnish nuclear families, obesity was linked to chromosome Xq24. Here we analyzed this 15-Mb region by genotyping 9 microsatellite markers and 36 single nucleotide polymorphisms (SNPs) for 11 positional and functional candidate genes in an extended sample of 218 obese Finnish sibling pairs (sibpairs) (BMI > 30 kg/m2). Evidence of linkage emerged mainly from the obese male sibpairs, suggesting a gender-specific effect for the underlying gene. By constructing haplotypes among the obese male sibpairs, we restricted the region from 15 Mb to 4 Mb, between markers DXS8088 and DXS8067. Regional functional candidate genes were tested for association in an initial sample of 117 cases and 182 controls. Significant evidence was observed for association for an SNP in the 3'-untranslated region of the solute carrier family 6 member 14 (SLC6A14) gene (P = 0.0002) and for SNP haplotypes of the SLC6A14 gene (P = 0.0007-0.006). Furthermore, an independent replication study sample of 837 cases and 968 controls from Finland and Sweden also showed significant differences in allele frequencies between obese and non-obese individuals (P = 0.003). The SLC6A14 gene is an interesting novel candidate for obesity because it encodes an amino acid transporter, which potentially regulates tryptophan availability for serotonin synthesis and thus possibly affects appetite control.

Genetic Diversity and Background Linkage Disequilibrium in the North American Holstein Cattle Population

The objectives of this study were to 1) identify highly heterozygous Holstein bulls that are as unrelated as possible and widely used in the US dairy industry; 2) quantify the level of genetic diversity in US Holsteins; and 3) determine the extent of background linkage disequilibrium (BLD) and disease trait associated linkage disequilibrium (DLD) in the US Holstein population. Twenty-three Holstein bulls that are not closely related but were widely used in the US dairy industry were genotyped for 54 microsatellite loci. The genotyping was performed on automated DNA sequencers (PE Applied Biosystems, CA), following polymerase chain reaction amplification with fluorescent dye-labeled primers. The heterozygosity for the sampled population ranged from 0.43 to 0.80. This wide range of heterozygosity allows selection of the most heterozygous bulls to develop informative families for gene mapping studies. The degree of genetic diversity in this population is significant and allows selection for traits of economic importance. As expected, there is extensive linkage disequilibrium (LD) in the US Holstein population. About half of the syntenic marker pairs presented a typical pattern of LD produced by DLD. Most of the nonsyntenic marker pairs had a typical pattern of LD arising from BLD. These results suggest that the observed LD is not purely due to genetic drift and migration and that a portion might be due to DLD. This raises our hopes of successful fine-localization of genes for complex traits using LD mapping.

The Newfoundland population: a unique resource for genetic investigation of complex diseases

The population of the province of Newfoundland and Labrador is genetically isolated. This isolation is evidenced by an overabundance of several monogenic disorders. The Newfoundland population, like that of other isolates, is now the focus of interest for identification of genes implicated in common diseases. However, the utility of such populations for this purpose remains unproven. In this paper, we review the current genetic architecture of the province, with respect to geographic isolation, homogeneity, founder effect, genetic drift and extended linkage disequilibrium. Based on these factors, we propose that the population of Newfoundland offers many advantages for genetic mapping of common diseases, compared with admixed populations, and even compared with other isolates.

Genetic linkage and association between chromosome 1q and working memory function in schizophrenia

Schizophrenia is substantially heritable, but specific susceptibility genes remain to be identified. Progress in this endeavor has been hindered by non-Mendelian transmission patterns, probable genetic heterogeneity, and an inability to detect premorbid and nonpenetrant carriers of predisposing genes. To circumvent these complexities, this study employed quantitative measures of liability, or "endophenotypes," within a sample of twins discordant for schizophrenia, drawn from the relatively genetically isolated population of Finland. A region on the distal portion of chromosome 1 has shown evidence for linkage to schizophrenia in two prior studies using Finnish patient samples. To elucidate further the nature and location of this potential susceptibility gene, linkage and association analyses were carried out across the chromosome 1 region of interest using quantitative neuropsychological measures of liability. Analyses with a composite measure of liability yielded suggestive evidence for linkage at marker D1S2833 (P = 0.04). Follow-up analyses of the individual trait measures showed that the Visual Span subtest of the Wechsler Memory Scale (WMS), an indicator of spatial working memory function, was uniquely sensitive to marker D1S2833 (P = 0.007). Association analysis confirmed that allelic variation in D1S2833 is associated with variation in spatial working memory performance as measured by the Visual Span subtest (P = 0.003), the significance of which was confirmed in an analysis of 10,000 Monte Carlo permutations. These data support the utility of this approach and provide evidence for a gene affecting spatial working memory function in schizophrenia patients and their unaffected co-twins.

Linkage analysis in psychiatric disorders: the emerging picture

Gene finding in genetically complex diseases has been difficult as a result of many factors that have diagnostic and methodologic considerations. For bipolar disorder and schizophrenia, numerous family, twin, and adoption studies have identified a strong genetic component to these behavioral psychiatric disorders. Despite difficulties that include diagnostic differences between sample populations and the lack of statistical significance in many individual studies, several promising patterns have emerged, suggesting that true susceptibility loci for schizophrenia and bipolar disorder may have been identified. In this review, the genetic epidemiology of these disorders is covered as well as linkage findings on chromosomes 4, 12, 13, 18, 21, and 22 in bipolar disorder and on chromosomes 1, 6, 8, 10, 13, 15, and 22 in schizophrenia. The sequencing of the human genome and identification of numerous single nucleotide polymorphisms (SNP) should substantially enhance the ability of investigators to identify disease-causing genes in these areas of the genome.

Manic-depression genes and the new millennium: poised for discovery

Manic-depressive illness is a common psychiatric disorder with complex etiology that likely involves multiple genes and non-genetic influences. The uncertain path to gene discovery has spurred considerable debate over genetic findings and gene-finding strategies. In this article, I review the main findings, with a focus on: (1) putative linked loci on chromosomes 1q31-32, 4p16, 6pter-p24, 10p14, 10q21-26, 12q23-24, 13q31-32, 18p11, 18q21-23, 21q22, 22q11-13, and Xq24-28; and (2) association studies with candidate genes, dynamic mutations, mitochondrial mutations, and chromosomal aberrations. Although no gene has been identified, promising findings are emerging. I then discuss the challenges and opportunities ahead, with special emphasis on gene-finding methods-in particular, questions pertaining to phenotype definition, linkage and association mapping, gene markers, sampling, study population, multigene systems, lessons from other disorders, animal models, and bioinformatics. The progress to date, together with rapid advances in genomics, analytical and computational methods, and bioinformatics, holds promise for new insights into the genetics of manic-depression, in the new millennium.