The lipidome in major depressive disorder: Shared genetic influence for ether-phosphatidylcholines, a plasma-based phenotype related to inflammation, and disease risk.

BACKGROUND: The lipidome is rapidly garnering interest in the field of psychiatry. Recent studies have implicated lipidomic changes across numerous psychiatric disorders. In particular, there is growing evidence that the concentrations of several classes of lipids are altered in those diagnosed with MDD. However, for lipidomic abnormalities to be considered potential treatment targets for MDD (rather than secondary manifestations of the disease), a shared etiology between lipid concentrations and MDD should be demonstrated. METHODS: In a sample of 567 individuals from 37 extended pedigrees (average size 13.57 people, range=3-80), we used mass spectrometry lipidomic measures to evaluate the genetic overlap between twenty-three biologically distinct lipid classes and a dimensional scale of MDD. RESULTS: We found that the lipid class with the largest endophenotype ranking value (ERV, a standardized parametric measure of pleiotropy) were ether-phosphodatidylcholines (alkylphosphatidylcholine, PC(O) and alkenylphosphatidylcholine, PC(P) subclasses). Furthermore, we examined the cluster structure of the twenty-five species within the top-ranked lipid class, and the relationship of those clusters with MDD. This analysis revealed that species containing arachidonic acid generally exhibited the greatest degree of genetic overlap with MDD. CONCLUSIONS: This study is the first to demonstrate a shared genetic etiology between MDD and ether-phosphatidylcholine species containing arachidonic acid, an omega-6 fatty acid that is a precursor to inflammatory mediators, such as prostaglandins. The study highlights the potential utility of the well-characterized linoleic/arachidonic acid inflammation pathway as a diagnostic marker and/or treatment target for MDD.

Transcriptome sequencing study implicates immune-related genes differentially expressed in schizophrenia: new data and a meta-analysis.

We undertook an RNA sequencing (RNAseq)-based transcriptomic profiling study on lymphoblastoid cell lines of a European ancestry sample of 529 schizophrenia cases and 660 controls, and found 1058 genes to be differentially expressed by affection status. These differentially expressed genes were enriched for involvement in immunity, especially the 697 genes with higher expression in cases. Comparing the current RNAseq transcriptomic profiling to our previous findings in an array-based study of 268 schizophrenia cases and 446 controls showed a highly significant positive correlation over all genes. Fifteen (18%) of the 84 genes with significant (false discovery rate<0.05) expression differences between cases and controls in the previous study and analyzed here again were differentially expressed by affection status here at a genome-wide significance level (Bonferroni P<0.05 adjusted for 8141 analyzed genes in total, or P<~6.1 × 10-6), all with the same direction of effect, thus providing corroborative evidence despite each sample of fully independent subjects being studied by different technological approaches. Meta-analysis of the RNAseq and array data sets (797 cases and 1106 controls) showed 169 additional genes (besides those found in the primary RNAseq-based analysis) to be differentially expressed, and provided further evidence of immune gene enrichment. In addition to strengthening our previous array-based gene expression differences in schizophrenia cases versus controls and providing transcriptomic support for some genes implicated by other approaches for schizophrenia, our study detected new genes differentially expressed in schizophrenia. We highlight RNAseq-based differential expression of various genes involved in neurodevelopment and/or neuronal function, and discuss caveats of the approach.

Genome-wide significant loci for addiction and anxiety.

BACKGROUND: Psychiatric comorbidity is common among individuals with addictive disorders, with patients frequently suffering from anxiety disorders. While the genetic architecture of comorbid addictive and anxiety disorders remains unclear, elucidating the genes involved could provide important insights into the underlying etiology. METHODS: Here we examine a sample of 1284 Mexican-Americans from randomly selected extended pedigrees. Variance decomposition methods were used to examine the role of genetics in addiction phenotypes (lifetime history of alcohol dependence, drug dependence or chronic smoking) and various forms of clinically relevant anxiety. Genome-wide univariate and bivariate linkage scans were conducted to localize the chromosomal regions influencing these traits. RESULTS: Addiction phenotypes and anxiety were shown to be heritable and univariate genome-wide linkage scans revealed significant quantitative trait loci for drug dependence (14q13.2-q21.2, LOD=3.322) and a broad anxiety phenotype (12q24.32-q24.33, LOD=2.918). Significant positive genetic correlations were observed between anxiety and each of the addiction subtypes (ρg=0.550-0.655) and further investigation with bivariate linkage analyses identified significant pleiotropic signals for alcohol dependence-anxiety (9q33.1-q33.2, LOD=3.054) and drug dependence-anxiety (18p11.23-p11.22, LOD=3.425). CONCLUSIONS: This study confirms the shared genetic underpinnings of addiction and anxiety and identifies genomic loci involved in the etiology of these comorbid disorders. The linkage signal for anxiety on 12q24 spans the location of TMEM132D, an emerging gene of interest from previous GWAS of anxiety traits, whilst the bivariate linkage signal identified for anxiety-alcohol on 9q33 peak coincides with a region where rare CNVs have been associated with psychiatric disorders. Other signals identified implicate novel regions of the genome in addiction genetics.

Detection of a quantitative trait locus associated with resistance to infection with Trichuris suis in pigs.

Whipworms (Trichuris spp.) infect a variety of hosts, including domestic animals and humans. Of considerable interest is the porcine whipworm, T. suis, which is particularly prevalent in outdoor production systems. High infection levels may cause growth retardation, anaemia and haemorrhagic diarrhoea. A significant proportion of the variation in Trichuris faecal egg count (FEC) has been attributed to the host's genetic make-up. The aim of the present study was to identify genetic loci associated with resistance to T. suis in pigs. We used single nucleotide polymorphism (SNP) markers to perform a whole-genome scan of an F1 resource population (n = 195) trickle-infected with T. suis. A measured genotype analysis revealed a putative quantitative trait locus (QTL) for T. suis FEC on chromosome 13 covering ∼ 4.5 Mbp, although none of the SNPs reached genome-wide significance. We tested the hypothesis that this region of SSC13 harboured genes with effects on T. suis burden by genotyping three SNPs within the putative QTL in unrelated pigs exposed to either experimental or natural T. suis infections and from which we had FEC (n = 113) or worm counts (n = 178). In these studies, two of the SNPs (rs55618716, ST) were associated with FEC (P < 0.01), thus confirming our initial findings. However, we did not find any of the SNPs to be associated with T. suis worm burden. In conclusion, our study demonstrates that genetic markers for resistance to T. suis as indicated by low FEC can be identified in pigs.

Evidence that duplications of 22q11.2 protect against schizophrenia.

A number of large, rare copy number variants (CNVs) are deleterious for neurodevelopmental disorders, but large, rare, protective CNVs have not been reported for such phenotypes. Here we show in a CNV analysis of 47 005 individuals, the largest CNV analysis of schizophrenia to date, that large duplications (1.5-3.0 Mb) at 22q11.2--the reciprocal of the well-known, risk-inducing deletion of this locus--are substantially less common in schizophrenia cases than in the general population (0.014% vs 0.085%, OR=0.17, P=0.00086). 22q11.2 duplications represent the first putative protective mutation for schizophrenia.

Epidemiology and genetic determinants of progressive deterioration of glycaemia in American Indians: the Strong Heart Family Study.

AIMS/HYPOTHESIS: Type 2 diabetes is a chronic, heterogeneous disease and a major risk factor for cardiovascular diseases. The underlying mechanisms leading to progression to type 2 diabetes are not fully understood and genetic tools may help to identify important pathways of glycaemic deterioration. METHODS: Using prospective data on American Indians from the Strong Heart Family Study, we identified 373 individuals defined as progressors (diabetes incident cases), 566 individuals with transitory impaired fasting glucose (IFG) and 1,011 controls (normal fasting glycaemia at all visits). We estimated the heritability (h(2)) of the traits and the evidence for association with 16 known variants identified in type 2 diabetes genome-wide association studies. RESULTS: We noted high h(2) for diabetes progression (h(2) = 0.65 ± 0.16, p = 2.7 × 10(-6)) but little contribution of genetic factors to transitory IFG (h(2) = 0.09 ± 0.10, p = 0.19) for models adjusted for multiple risk factors. At least three variants (in WFS1, TSPAN8 and THADA) were nominally associated with diabetes progression in age- and sex-adjusted analyses with estimates showing the same direction of effects as reported in the discovery European ancestry studies. CONCLUSIONS/INTERPRETATION: Our findings do not exclude these loci for diabetes susceptibility in American Indians and suggest phenotypic heterogeneity of the IFG trait, which may have implications for genetic studies when diagnosis is based on a single time-point measure.

Transcriptomics of cortical gray matter thickness decline during normal aging.

INTRODUCTION: We performed a whole-transcriptome correlation analysis, followed by the pathway enrichment and testing of innate immune response pathway analyses to evaluate the hypothesis that transcriptional activity can predict cortical gray matter thickness (GMT) variability during normal cerebral aging. METHODS: Transcriptome and GMT data were available for 379 individuals (age range=28-85) community-dwelling members of large extended Mexican American families. Collection of transcriptome data preceded that of neuroimaging data by 17 years. Genome-wide gene transcriptome data consisted of 20,413 heritable lymphocytes-based transcripts. GMT measurements were performed from high-resolution (isotropic 800 µm) T1-weighted MRI. Transcriptome-wide and pathway enrichment analysis was used to classify genes correlated with GMT. Transcripts for sixty genes from seven innate immune pathways were tested as specific predictors of GMT variability. RESULTS: Transcripts for eight genes (IGFBP3, LRRN3, CRIP2, SCD, IDS, TCF4, GATA3, and HN1) passed the transcriptome-wide significance threshold. Four orthogonal factors extracted from this set predicted 31.9% of the variability in the whole-brain and between 23.4 and 35% of regional GMT measurements. Pathway enrichment analysis identified six functional categories including cellular proliferation, aggregation, differentiation, viral infection, and metabolism. The integrin signaling pathway was significantly (p<10(-6)) enriched with GMT. Finally, three innate immune pathways (complement signaling, toll-receptors and scavenger and immunoglobulins) were significantly associated with GMT. CONCLUSION: Expression activity for the genes that regulate cellular proliferation, adhesion, differentiation and inflammation can explain a significant proportion of individual variability in cortical GMT. Our findings suggest that normal cerebral aging is the product of a progressive decline in regenerative capacity and increased neuroinflammation.

QTL-based association analyses reveal novel genes influencing pleiotropy of metabolic syndrome (MetS).

OBJECTIVE: Metabolic Syndrome (MetS) is a phenotype cluster predisposing to type 2 diabetes and cardiovascular disease. We conducted a study to elucidate the genetic basis underlying linkage signals for multiple representative traits of MetS that we had previously identified at two significant QTLs on chromosomes 3q27 and 17p12. DESIGN AND METHODS: We performed QTL-specific genomic and transcriptomic analyses in 1,137 individuals from 85 extended families that contributed to the original linkage. We tested in SOLAR association of MetS phenotypes with QTL-specific haplotype-tagging SNPs as well as transcriptional profiles of peripheral blood mononuclear cells (PBMCs). RESULTS: SNPs significantly associated with MetS phenotypes under the prior hypothesis of linkage mapped to seven genes at 3q27 and seven at 17p12. Prioritization based on biologic relevance, SNP association, and expression analyses identified two genes: insulin-like growth factor 2 mRNA-binding protein 2 (IGF2BP2) at 3q27 and tumor necrosis factor receptor 13B (TNFRSF13B) at 17p12. Prioritized genes could influence cell-cell adhesion and adipocyte differentiation, insulin/glucose responsiveness, cytokine effectiveness, plasma lipid levels, and lipoprotein densities. CONCLUSIONS: Using an approach combining genomic, transcriptomic, and bioinformatic data we identified novel candidate genes for MetS.

Impact of DISC1 variation on neuroanatomical and neurocognitive phenotypes.

Although disrupted in schizophrenia 1 (DISC1) has been implicated in many psychiatric disorders, including schizophrenia, bipolar disorder, schizoaffective disorder and major depression, its biological role in these disorders is unclear. To better understand this gene and its role in psychiatric disease, we conducted transcriptional profiling and genome-wide association analysis in 1232 pedigreed Mexican-American individuals for whom we have neuroanatomic images, neurocognitive assessments and neuropsychiatric diagnoses. SOLAR was used to determine heritability, identify gene expression patterns and perform association analyses on 188 quantitative brain-related phenotypes. We found that the DISC1 transcript is highly heritable (h(2)=0.50; P=1.97 × 10(-22)), and that gene expression is strongly cis-regulated (cis-LOD=3.89) but is also influenced by trans-effects. We identified several DISC1 polymorphisms that were associated with cortical gray matter thickness within the parietal, temporal and frontal lobes. Associated regions affiliated with memory included the entorhinal cortex (rs821639, P=4.11 × 10(-5); rs2356606, P=4.71 × 10(-4)), cingulate cortex (rs16856322, P=2.88 × 10(-4)) and parahippocampal gyrus (rs821639, P=4.95 × 10(-4)); those affiliated with executive and other cognitive processing included the transverse temporal gyrus (rs9661837, P=5.21 × 10(-4); rs17773946, P=6.23 × 10(-4)), anterior cingulate cortex (rs2487453, P=4.79 × 10(-4); rs3738401, P=5.43 × 10(-4)) and medial orbitofrontal cortex (rs9661837; P=7.40 × 10(-4)). Cognitive measures of working memory (rs2793094, P=3.38 × 10(-4)), as well as lifetime history of depression (rs4658966, P=4.33 × 10(-4); rs12137417, P=4.93 × 10(-4)) and panic (rs12137417, P=7.41 × 10(-4)) were associated with DISC1 sequence variation. DISC1 has well-defined genetic regulation and clearly influences important phenotypes related to psychiatric disease.

Genome-wide linkage scan for plasma high density lipoprotein cholesterol, apolipoprotein A-1 and triglyceride variation among American Indian populations: the Strong Heart Family Study.

BACKGROUND: Recent studies have identified chromosomal regions linked to variation in high density lipoprotein cholesterol (HDL-C), apolipoprotein A-1 (apo A-1) and triglyceride (TG), although results have been inconsistent and previous studies of American Indian populations are limited. OBJECTIVE: In an attempt to localise quantitative trait loci (QTLs) influencing HDL-C, apo A-1 and TG, we conducted genome-wide linkage scans of subjects of the Strong Heart Family Study. METHODS: We implemented analyses in 3484 men and women aged 18 years or older, at three study centres. RESULTS: With adjustment for age, sex and centre, we detected a QTL influencing both HDL-C (logarithm of odds (LOD) = 4.4, genome-wide p = 0.001) and apo A-1 (LOD = 3.2, genome-wide p = 0.020) nearest marker D6S289 at 6p23 in the Arizona sample. Another QTL influencing apo A-1 was found nearest marker D9S287 at 9q22.2 (LOD = 3.0, genome-wide p = 0.033) in the North and South Dakotas. We detected a QTL influencing TG nearest marker D15S153 at 15q22.31 (LOD = 4.5 in the overall sample and LOD = 3.8 in the Dakotas sample, genome-wide p = 0.0044) and when additionally adjusted for waist, current smoking, current alcohol, current oestrogen, lipid treatment, impaired fasting glucose, and diabetes, nearest marker D10S217 at 10q26.2 (LOD = 3.7, genome-wide p = 0.0058) in the Arizona population. CONCLUSIONS: The replication of QTLs in regions of the genome that harbour well known candidate genes suggest that chromosomes 6p, 9q and 15q warrant further investigation with fine mapping for causative polymorphisms in American Indians.

High heritability for Ascaris and Trichuris infection levels in pigs.

Aggregated distributions of macroparasites within their host populations are characteristic of most natural and experimental infections. We designed this study to measure the amount of variation that is attributable to host genetic factors in a pig-helminth system. In total, 195 piglets were produced after artificial insemination of 19 sows (Danish Landrace-Yorkshire crossbreds) with semen selected from 13 individual Duroc boars (1 or 2 sows per boar; mean litter size: 10.3; 5-14 piglets per litter). Starting at 10 weeks of age, piglets were repeatedly infected with the gastrointestinal helminths Trichuris suis and Ascaris suum by administering eggs in the feed for 14 weeks until necropsy. Faecal egg counts (FECs) were estimated regularly and A. suum worm burden was obtained at necropsy. Heritability calculations for log (FEC+1) at weeks 7-10 post-infection (p.i.) showed that 0.32-0.73 of the phenotypic variation for T. suis could be attributed to genetic factors. For A. suum, heritabilities of 0.29-0.31 were estimated for log (FEC+1) at weeks 7-14 p.i., whereas the heritability of log worm counts was 0.45. Strong positive genetic correlations (0.75-0.89) between T. suis and A. suum FECs suggest that resistance to both infections involves regulation by overlapping genes. Our data demonstrate that there is a strong genetic component in resistance to A. suum and T. suis infections in pigs. Identification of responsible genes would enhance our understanding of the host immune response to these common nematodes and for the closely related species (T. trichiura and A. lumbricoides) in man infecting more than a billion people.

Cross-species replication of a resistin mRNA QTL, but not QTLs for circulating levels of resistin, in human and baboon.

Resistin has been associated with inflammation and risk for cardiovascular disease. We previously reported evidence of a QTL on chromosome 19p13 affecting the abundance of resistin (RETN) mRNA in the omental adipose tissue of baboons (L0D score 3.8). In this study, whole genome transcription levels were assessed in human lymphocyte samples from 1240 adults participating in the San Antonio Family Heart Study, using the Sentrix Human-6 Expression Beadchip. Lymphocytes were surveyed, as it has been proposed that their expression levels may reflect those in harder to ascertain tissues, such as adipose tissue, that are thought to be more directly relevant to disease procesn was conducted to detect loci affecting RETN mRNA levels. We obtained significant evidence for a QTL influencing the RETN expression (LOD score 10.7) on chromosome 19p. This region is orthologous/homologous to the region previously localized on baboon chromosome 19. The strongest positional candidate gene in this region is the structural gene for resistin, itself. We also found evidence for a QTL influencing resistin protein levels (LOD score 5.3) on chromosome 14q. This differs from our previously reported QTL on chromosome 18 in baboons. The different QTLs for circulating protein suggests that post-translational processing and turnover may be influenced by different or multiple genes in baboons and humans. The parallel findings of a cis-eQTL for RETN mRNA in baboon omental tissue and human lymphocytes lends support to the strategy of using lymphocyte gene expression levels as a surrogate for gene expression levels in other tissues.

Genome-wide linkage scan for loci of musical aptitude in Finnish families: evidence for a major locus at 4q22.

BACKGROUND: Music perception and performance are comprehensive human cognitive functions and thus provide an excellent model system for studying human behaviour and brain function. However, the molecules involved in mediating music perception and performance are so far uncharacterised. OBJECTIVE: To unravel the biological background of music perception, using molecular and statistical genetic approaches. METHODS: 15 Finnish multigenerational families (with a total of 234 family members) were recruited via a nationwide search. The phenotype of all family members was determined using three tests used in defining musical aptitude: a test for auditory structuring ability (Karma Music test; KMT) commonly used in Finland, and the Seashore pitch and time discrimination subtests (SP and ST respectively) used internationally. We calculated heritabilities and performed a genome-wide variance components-based linkage scan using genotype data for 1113 microsatellite markers. RESULTS: The heritability estimates were 42% for KMT, 57% for SP, 21% for ST and 48% for the combined music test scores. Significant evidence of linkage was obtained on chromosome 4q22 (LOD 3.33) and suggestive evidence of linkage at 8q13-21 (LOD 2.29) with the combined music test scores, using variance component linkage analyses. The major contribution of the 4q22 locus was obtained for the KMT (LOD 2.91). Interestingly, a positive LOD score of 1.69 was shown at 18q, a region previously linked to dyslexia (DYX6) using combined music test scores. CONCLUSION: Our results show that there is a genetic contribution to musical aptitude that is likely to be regulated by several predisposing genes or variants.

Linkage to CFA29 detected in a genome-wide linkage screen of a canine pedigree segregating Sry-negative XX sex reversal.

Canine Sry-negative XX sex reversal is a disorder of gonadal development wherein individuals having a female karyotype develop testes or ovotestes. In this study, linkage mapping was undertaken in a pedigree derived from one proven carrier American cocker spaniel founder male and beagle females. All affected dogs in the analysis were XX true hermaphrodites and confirmed to be Sry negative by polymerase chain reaction. A genome-wide linkage screen conducted using 245 microsatellite markers revealed highest LOD score of 3.4 (marker CPH9) on CFA29. Fine mapping with additional microsatellites in the region containing CPH9 localized the Sry-negative XX sex reversal locus to a 5.4-Mb candidate region between markers CPH9 and FH3003 (LOD score 3.15). Insignificant LOD scores were found at genome-wide screen or fine mapping markers that were within 10 Mb of 45 potential candidate genes reported to have a role in mammalian sex determination or differentiation. Together, these results suggest that a novel locus on CFA29 may be responsible for sex reversal in this pedigree.

Linkage analysis of circulating levels of adiponectin in Hispanic children.

INTRODUCTION: Adiponectin, a hormone produced exclusively by adipose tissue, is inversely associated with insulin resistance and proinflammatory conditions. The aim of this study was to find quantitative trait loci (QTLs) that affect circulating levels of adiponectin in Hispanic children participating in the VIVA LA FAMILIA Study by use of a systematic genome scan. METHODS: The present study included extended families with at least one overweight child between 4 and 19 years old. Overweight was defined as body mass index (BMI) 95th percentile. Fasting blood was collected from 466 children from 127 families. Adiponectin was assayed by radioimmunoassay (RIA) technique in fasting serum. A genome-wide scan on circulating levels of adiponectin as a quantitative phenotype was conducted using the variance decomposition approach. RESULTS: The highest logarithm of odds (LOD) score (4.2) was found on chromosome 11q23.2-11q24.2, and a second significant signal (LOD score=3.0) was found on chromosome 8q12.1-8q21.3. In addition, a signal suggestive of linkage (LOD score=2.5) was found between 18q21.3 and 18q22.3. After adjustment for BMI-Z score, the LOD score on chromosome 11 remained unchanged, but the signals on chromosomes 8 and 18 dropped to 1.6 and 1.7, respectively. Two other signals suggestive of linkage were found on chromosome 3 (LOD score=2.1) and 10 (LOD score=2.5). Although the region on chromosome 11 has been associated with obesity and diabetes-related traits in adult populations, this is the first observation of linkage in this region for adiponectin levels. Our suggestive linkages on chromosomes 10 and 3 replicate results for adiponectin seen in other populations. The influence of loci on chromosomes 18 and 8 on circulating adiponectin seemed to be mediated by BMI in the present study. CONCLUSION: Our genome scan in children has identified a novel QTL and replicated QTLs in chromosomal regions previously shown to be linked with obesity and type 2 diabetes (T2D)-related phenotypes in adults. The genetic contribution of loci to adiponectin levels may vary across different populations and age groups. The strong linkage signal on chromosome 11 is most likely underlain by a gene(s) that may contribute to the high susceptibility of these Hispanic children to obesity and T2D.

Linkage analysis of LDL cholesterol in American Indian populations: the Strong Heart Family Study.

Previous studies have demonstrated that low density lipoprotein cholesterol (LDL-C) concentration is influenced by both genes and environment. Although rare genetic variants associated with Mendelian causes of increased LDL-C are known, only one common genetic variant has been identified, the apolipoprotein E gene (APOE). In an attempt to localize quantitative trait loci (QTLs) influencing LDL-C, we conducted a genome-wide linkage scan of LDL-C in participants of the Strong Heart Family Study (SHFS). Nine hundred eighty men and women, age 18 years or older, in 32 extended families at three centers (in Arizona, Oklahoma, and North and South Dakota) were phenotyped for LDL-C concentration and other risk factors. Using a variance component approach and the program SOLAR, and after accounting for the effects of covariates, we detected a QTL influencing LDL-C on chromosome 19, nearest marker D19S888 at 19q13.41 [logarithm of odds (LOD) = 4.3] in the sample from the Dakotas. This region on chromosome 19 includes many possible candidate genes, including the APOE/C1/C4/C2 gene cluster. In follow-up association analyses, no significant evidence for an association was detected with the APOE*2 and APOE*4 alleles (P = 0.76 and P = 0.53, respectively). Suggestive evidence of linkage to LDL-C was detected on chromosomes 3q, 4q, 7p, 9q, 10p, 14q, and 17q. These linkage signals overlap positive findings for lipid-related traits and harbor plausible candidate genes for LDL-C.

The role of sequence variations within the genes encoding collagen II, IX and XI in non-syndromic, early-onset osteoarthritis.

OBJECTIVE: We sought to determine whether sequence variations in cartilage collagen genes are associated with primary, early-onset osteoarthritis (OA). METHODS: The cartilage collagen genes, COL2A1, COL9A1, COL9A2, COL9A3, COL11A1 and COL11A2, were screened for sequence variations in 72 Finnish probands and one US family with primary early-onset hip and/or knee OA. In addition, allelic association studies were performed using six to 12 common polymorphisms from each gene by genotyping 72 OA patients and 103 controls. RESULTS: Altogether 239 sequence variations were found, of which 16 were not present in the controls. Seven of the unique variations, four in COL11A1, two in COL11A2 and one in COL2A1, were studied further, because they resulted in the substitution of conserved amino acids or were predicted to affect mRNA splicing. Co-segregation of a sequence variation and the phenotype was found in all four families available for study. Association analysis failed to identify any common predisposing alleles. CONCLUSIONS: Early-onset OA demonstrates locus and allelic heterogeneity since the identified variations were in three different collagen genes and each of the six probands had a different mutation. It is also possible that some OA cases represent the mild end of the chondrodysplasia phenotypic spectrum. The major susceptibility alleles in this form of OA, however, remain to be identified.

Fine mapping and SNP analysis of positional candidates at the preeclampsia susceptibility locus (PREG1) on chromosome 2.

Genome scans in Icelandic, Australian and New Zealand, and Finnish families have localized putative susceptibility loci for preeclampsia/ eclampsia to chromosome 2. The locus mapped in the Australian and New Zealand study (designated PREG1) was thought to be the same locus as that identified in the Icelandic study. In both these studies, two distinct quantitative trait locus (QTL) regions were evident on chromosome 2. Here, we describe our fine mapping of the PREG1 locus and a genetic analysis of two positional candidate genes. Twenty-five additional microsatellite markers were genotyped within the 74-cM linkage region defined by the combined Icelandic and Australian and New Zealand genome scans. The overall position and shape of the localization evidence obtained using nonparametric multipoint analysis did not change from that seen previously in our 10-cM resolution genome scan; two peaks were displayed, one on chromosome 2p at marker D2S388 (107.46 cM) and the other on chromosome 2q at 151.5 cM at marker D2S2313. Using the robust two-point linkage analysis implemented in the Analyze program, all 25 markers gave positive LOD scores with significant evidence of linkage being seen at marker D2S2313 (151.5 cM), achieving a LOD score of 3.37 under a strict diagnostic model. Suggestive evidence of linkage was seen at marker D2S388 (107.46 cM) with a LOD score of 2.22 under the general diagnostic model. Two candidate genes beneath the peak on chromosome 2p were selected for further analysis using public single nucleotide polymorphisms (SNPs) within these genes. Maximum LOD scores were obtained for an SNP in TACR1 (LOD = 3.5) and for an SNP in TCF7L1 (LOD = 3.33), both achieving genome-wide significance. However, no evidence of association was seen with any of the markers tested. These data strongly support the presence of a susceptibility gene on chromosome 2p11-12 and substantiate the possibility of a second locus on chromosome 2q23.

Large upward bias in estimation of locus-specific effects from genomewide scans.

The primary goal of a genomewide scan is to estimate the genomic locations of genes influencing a trait of interest. It is sometimes said that a secondary goal is to estimate the phenotypic effects of each identified locus. Here, it is shown that these two objectives cannot be met reliably by use of a single data set of a currently realistic size. Simulation and analytical results, based on variance-components linkage analysis as an example, demonstrate that estimates of locus-specific effect size at genomewide LOD score peaks tend to be grossly inflated and can even be virtually independent of the true effect size, even for studies on large samples when the true effect size is small. However, the bias diminishes asymptotically. The explanation for the bias is that the LOD score is a function of the locus-specific effect-size estimate, such that there is a high correlation between the observed statistical significance and the effect-size estimate. When the LOD score is maximized over the many pointwise tests being conducted throughout the genome, the locus-specific effect-size estimate is therefore effectively maximized as well. We argue that attempts at bias correction give unsatisfactory results, and that pointwise estimation in an independent data set may be the only way of obtaining reliable estimates of locus-specific effect-and then only if one does not condition on statistical significance being obtained. We further show that the same factors causing this bias are responsible for frequent failures to replicate initial claims of linkage or association for complex traits, even when the initial localization is, in fact, correct. The findings of this study have wide-ranging implications, as they apply to all statistical methods of gene localization. It is hoped that, by keeping this bias in mind, we will more realistically interpret and extrapolate from the results of genomewide scans.

A spectrum of ABCC6 mutations is responsible for pseudoxanthoma elasticum.

To better understand the pathogenetics of pseudoxanthoma elasticum (PXE), we performed a mutational analysis of ATP-binding cassette subfamily C member 6 (ABCC6) in 122 unrelated patients with PXE, the largest cohort of patients yet studied. Thirty-six mutations were characterized, and, among these, 28 were novel variants (for a total of 43 PXE mutations known to date). Twenty-one alleles were missense variants, six were small insertions or deletions, five were nonsense, two were alleles likely to result in aberrant mRNA splicing, and two were large deletions involving ABCC6. Although most mutations appeared to be unique variants, two disease-causing alleles occurred frequently in apparently unrelated individuals. R1141X was found in our patient cohort at a frequency of 18.8% and was preponderant in European patients. ABCC6del23-29 occurred at a frequency of 12.9% and was prevalent in patients from the United States. These results suggested that R1141X and ABCC6del23-29 might have been derived regionally from founder alleles. Putative disease-causing mutations were identified in approximately 64% of the 244 chromosomes studied, and 85.2% of the 122 patients were found to have at least one disease-causing allele. Our results suggest that a fraction of the undetected mutant alleles could be either genomic rearrangements or mutations occurring in noncoding regions of the ABCC6 gene. The distribution pattern of ABCC6 mutations revealed a cluster of disease-causing variants within exons encoding a large C-terminal cytoplasmic loop and in the C-terminal nucleotide-binding domain (NBD2). We discuss the potential structural and functional significance of this mutation pattern within the context of the complex relationship between the PXE phenotype and the function of ABCC6.