Genome-wide mapping of brain phenotypes in extended pedigrees with strong genetic loading for bipolar disorder.

Bipolar disorder is a highly heritable illness, associated with alterations of brain structure. As such, identification of genes influencing inter-individual differences in brain morphology may help elucidate the underlying pathophysiology of bipolar disorder (BP). To identify quantitative trait loci (QTL) that contribute to phenotypic variance of brain structure, structural neuroimages were acquired from family members (n = 527) of extended pedigrees heavily loaded for bipolar disorder ascertained from genetically isolated populations in Latin America. Genome-wide linkage and association analysis were conducted on the subset of heritable brain traits that showed significant evidence of association with bipolar disorder (n = 24) to map QTL influencing regional measures of brain volume and cortical thickness. Two chromosomal regions showed significant evidence of linkage; a QTL on chromosome 1p influencing corpus callosum volume and a region on chromosome 7p linked to cortical volume. Association analysis within the two QTLs identified three SNPs correlated with the brain measures.

Author Correction: The evolutionary history of lethal metastatic prostate cancer.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Contribution of common and rare variants to bipolar disorder susceptibility in extended pedigrees from population isolates.

Current evidence from case/control studies indicates that genetic risk for psychiatric disorders derives primarily from numerous common variants, each with a small phenotypic impact. The literature describing apparent segregation of bipolar disorder (BP) in numerous multigenerational pedigrees suggests that, in such families, large-effect inherited variants might play a greater role. To identify roles of rare and common variants on BP, we conducted genetic analyses in 26 Colombia and Costa Rica pedigrees ascertained for bipolar disorder 1 (BP1), the most severe and heritable form of BP. In these pedigrees, we performed microarray SNP genotyping of 838 individuals and high-coverage whole-genome sequencing of 449 individuals. We compared polygenic risk scores (PRS), estimated using the latest BP1 genome-wide association study (GWAS) summary statistics, between BP1 individuals and related controls. We also evaluated whether BP1 individuals had a higher burden of rare deleterious single-nucleotide variants (SNVs) and rare copy number variants (CNVs) in a set of genes related to BP1. We found that compared with unaffected relatives, BP1 individuals had higher PRS estimated from BP1 GWAS statistics (P = 0.001 ~ 0.007) and displayed modest increase in burdens of rare deleterious SNVs (P = 0.047) and rare CNVs (P = 0.002 ~ 0.033) in genes related to BP1. We did not observe rare variants segregating in the pedigrees. These results suggest that small-to-moderate effect rare and common variants are more likely to contribute to BP1 risk in these extended pedigrees than a few large-effect rare variants.

Sequencing of prostate cancers identifies new cancer genes, routes of progression and drug targets.

Prostate cancer represents a substantial clinical challenge because it is difficult to predict outcome and advanced disease is often fatal. We sequenced the whole genomes of 112 primary and metastatic prostate cancer samples. From joint analysis of these cancers with those from previous studies (930 cancers in total), we found evidence for 22 previously unidentified putative driver genes harboring coding mutations, as well as evidence for NEAT1 and FOXA1 acting as drivers through noncoding mutations. Through the temporal dissection of aberrations, we identified driver mutations specifically associated with steps in the progression of prostate cancer, establishing, for example, loss of CHD1 and BRCA2 as early events in cancer development of ETS fusion-negative cancers. Computational chemogenomic (canSAR) analysis of prostate cancer mutations identified 11 targets of approved drugs, 7 targets of investigational drugs, and 62 targets of compounds that may be active and should be considered candidates for future clinical trials.

PHIP - a novel candidate breast cancer susceptibility locus on 6q14.1.

Most non-BRCA1/2 breast cancer families have no identified genetic cause. We used linkage and haplotype analyses in familial and sporadic breast cancer cases to identify a susceptibility locus on chromosome 6q. Two independent genome-wide linkage analysis studies suggested a 3 Mb locus on chromosome 6q and two unrelated Swedish families with a LOD >2 together seemed to share a haplotype in 6q14.1. We hypothesized that this region harbored a rare high-risk founder allele contributing to breast cancer in these two families. Sequencing of DNA and RNA from the two families did not detect any pathogenic mutations. Finally, 29 SNPs in the region were analyzed in 44,214 cases and 43,532 controls from BCAC, and the original haplotypes in the two families were suggested as low-risk alleles for European and Swedish women specifically. There was also some support for one additional independent moderate-risk allele in Swedish familial samples. The results were consistent with our previous findings in familial breast cancer and supported a breast cancer susceptibility locus at 6q14.1 around the PHIP gene.

Appraising the relevance of DNA copy number loss and gain in prostate cancer using whole genome DNA sequence data.

A variety of models have been proposed to explain regions of recurrent somatic copy number alteration (SCNA) in human cancer. Our study employs Whole Genome DNA Sequence (WGS) data from tumor samples (n = 103) to comprehensively assess the role of the Knudson two hit genetic model in SCNA generation in prostate cancer. 64 recurrent regions of loss and gain were detected, of which 28 were novel, including regions of loss with more than 15% frequency at Chr4p15.2-p15.1 (15.53%), Chr6q27 (16.50%) and Chr18q12.3 (17.48%). Comprehensive mutation screens of genes, lincRNA encoding sequences, control regions and conserved domains within SCNAs demonstrated that a two-hit genetic model was supported in only a minor proportion of recurrent SCNA losses examined (15/40). We found that recurrent breakpoints and regions of inversion often occur within Knudson model SCNAs, leading to the identification of ZNF292 as a target gene for the deletion at 6q14.3-q15 and NKX3.1 as a two-hit target at 8p21.3-p21.2. The importance of alterations of lincRNA sequences was illustrated by the identification of a novel mutational hotspot at the KCCAT42, FENDRR, CAT1886 and STCAT2 loci at the 16q23.1-q24.3 loss. Our data confirm that the burden of SCNAs is predictive of biochemical recurrence, define nine individual regions that are associated with relapse, and highlight the possible importance of ion channel and G-protein coupled-receptor (GPCR) pathways in cancer development. We concluded that a two-hit genetic model accounts for about one third of SCNA indicating that mechanisms, such haploinsufficiency and epigenetic inactivation, account for the remaining SCNA losses.

Genetic contributions to circadian activity rhythm and sleep pattern phenotypes in pedigrees segregating for severe bipolar disorder.

Abnormalities in sleep and circadian rhythms are central features of bipolar disorder (BP), often persisting between episodes. We report here, to our knowledge, the first systematic analysis of circadian rhythm activity in pedigrees segregating severe BP (BP-I). By analyzing actigraphy data obtained from members of 26 Costa Rican and Colombian pedigrees [136 euthymic (i.e., interepisode) BP-I individuals and 422 non-BP-I relatives], we delineated 73 phenotypes, of which 49 demonstrated significant heritability and 13 showed significant trait-like association with BP-I. All BP-I-associated traits related to activity level, with BP-I individuals consistently demonstrating lower activity levels than their non-BP-I relatives. We analyzed all 49 heritable phenotypes using genetic linkage analysis, with special emphasis on phenotypes judged to have the strongest impact on the biology underlying BP. We identified a locus for interdaily stability of activity, at a threshold exceeding genome-wide significance, on chromosome 12pter, a region that also showed pleiotropic linkage to two additional activity phenotypes.

The evolutionary history of lethal metastatic prostate cancer.

Cancers emerge from an ongoing Darwinian evolutionary process, often leading to multiple competing subclones within a single primary tumour. This evolutionary process culminates in the formation of metastases, which is the cause of 90% of cancer-related deaths. However, despite its clinical importance, little is known about the principles governing the dissemination of cancer cells to distant organs. Although the hypothesis that each metastasis originates from a single tumour cell is generally supported, recent studies using mouse models of cancer demonstrated the existence of polyclonal seeding from and interclonal cooperation between multiple subclones. Here we sought definitive evidence for the existence of polyclonal seeding in human malignancy and to establish the clonal relationship among different metastases in the context of androgen-deprived metastatic prostate cancer. Using whole-genome sequencing, we characterized multiple metastases arising from prostate tumours in ten patients. Integrated analyses of subclonal architecture revealed the patterns of metastatic spread in unprecedented detail. Metastasis-to-metastasis spread was found to be common, either through de novo monoclonal seeding of daughter metastases or, in five cases, through the transfer of multiple tumour clones between metastatic sites. Lesions affecting tumour suppressor genes usually occur as single events, whereas mutations in genes involved in androgen receptor signalling commonly involve multiple, convergent events in different metastases. Our results elucidate in detail the complex patterns of metastatic spread and further our understanding of the development of resistance to androgen-deprivation therapy in prostate cancer.

Analysis of the genetic phylogeny of multifocal prostate cancer identifies multiple independent clonal expansions in neoplastic and morphologically normal prostate tissue.

Genome-wide DNA sequencing was used to decrypt the phylogeny of multiple samples from distinct areas of cancer and morphologically normal tissue taken from the prostates of three men. Mutations were present at high levels in morphologically normal tissue distant from the cancer, reflecting clonal expansions, and the underlying mutational processes at work in morphologically normal tissue were also at work in cancer. Our observations demonstrate the existence of ongoing abnormal mutational processes, consistent with field effects, underlying carcinogenesis. This mechanism gives rise to extensive branching evolution and cancer clone mixing, as exemplified by the coexistence of multiple cancer lineages harboring distinct ERG fusions within a single cancer nodule. Subsets of mutations were shared either by morphologically normal and malignant tissues or between different ERG lineages, indicating earlier or separate clonal cell expansions. Our observations inform on the origin of multifocal disease and have implications for prostate cancer therapy in individual cases.

Multisystem component phenotypes of bipolar disorder for genetic investigations of extended pedigrees.

IMPORTANCE: Genetic factors contribute to risk for bipolar disorder (BP), but its pathogenesis remains poorly understood. A focus on measuring multisystem quantitative traits that may be components of BP psychopathology may enable genetic dissection of this complex disorder, and investigation of extended pedigrees from genetically isolated populations may facilitate the detection of specific genetic variants that affect BP as well as its component phenotypes. OBJECTIVE: To identify quantitative neurocognitive, temperament-related, and neuroanatomical phenotypes that appear heritable and associated with severe BP (bipolar I disorder [BP-I]) and therefore suitable for genetic linkage and association studies aimed at identifying variants contributing to BP-I risk. DESIGN, SETTING, AND PARTICIPANTS: Multigenerational pedigree study in 2 closely related, genetically isolated populations: the Central Valley of Costa Rica and Antioquia, Colombia. A total of 738 individuals, all from Central Valley of Costa Rica and Antioquia pedigrees, participated; among them, 181 have BP-I. MAIN OUTCOMES AND MEASURES: Familial aggregation (heritability) and association with BP-I of 169 quantitative neurocognitive, temperament, magnetic resonance imaging, and diffusion tensor imaging phenotypes. RESULTS: Of 169 phenotypes investigated, 126 (75%) were significantly heritable and 53 (31%) were associated with BP-I. About one-quarter of the phenotypes, including measures from each phenotype domain, were both heritable and associated with BP-I. Neuroimaging phenotypes, particularly cortical thickness in prefrontal and temporal regions as well as volume and microstructural integrity of the corpus callosum, represented the most promising candidate traits for genetic mapping related to BP based on strong heritability and association with disease. Analyses of phenotypic and genetic covariation identified substantial correlations among the traits, at least some of which share a common underlying genetic architecture. CONCLUSIONS AND RELEVANCE: To our knowledge, this is the most extensive investigation of BP-relevant component phenotypes to date. Our results identify brain and behavioral quantitative traits that appear to be genetically influenced and show a pattern of BP-I association within families that is consistent with expectations from case-control studies. Together, these phenotypes provide a basis for identifying loci contributing to BP-I risk and for genetic dissection of the disorder.

CNV analysis in Tourette syndrome implicates large genomic rearrangements in COL8A1 and NRXN1.

Tourette syndrome (TS) is a neuropsychiatric disorder with a strong genetic component. However, the genetic architecture of TS remains uncertain. Copy number variation (CNV) has been shown to contribute to the genetic make-up of several neurodevelopmental conditions, including schizophrenia and autism. Here we describe CNV calls using SNP chip genotype data from an initial sample of 210 TS cases and 285 controls ascertained in two Latin American populations. After extensive quality control, we found that cases (N = 179) have a significant excess (P = 0.006) of large CNV (>500 kb) calls compared to controls (N = 234). Amongst 24 large CNVs seen only in the cases, we observed four duplications of the COL8A1 gene region. We also found two cases with ∼400 kb deletions involving NRXN1, a gene previously implicated in neurodevelopmental disorders, including TS. Follow-up using multiplex ligation-dependent probe amplification (and including 53 more TS cases) validated the CNV calls and identified additional patients with rearrangements in COL8A1 and NRXN1, but none in controls. Examination of available parents indicates that two out of three NRXN1 deletions detected in the TS cases are de-novo mutations. Our results are consistent with the proposal that rare CNVs play a role in TS aetiology and suggest a possible role for rearrangements in the COL8A1 and NRXN1 gene regions.

Asociación de esquizofrenia y sus dimensiones clínicas con el gen NOS1AP en población colombiana / Association of schizophrenia and its clinical implications with the nos1ap gene in the colombian population

Introducción: El gen NOS1AP codifica para la proteína adaptadora de óxido nítrico sintasa neuronal 1, que posiblemente está implicada en la etiopatogénesis de la esquizofrenia. Objetivos: Determinar si existe asociación de variantes en el gen NOS1AP con esquizofrenia y si estas variantes tienen relación con las dimensiones clínicas del trastorno en población colombiana. Metodología: Es un estudio de casos y controles, con 255 sujetos por grupo. Se tipificaron marcadores dentro del gen NOS1AP y otros informativos de origen genético, con el fin de ajustar por estratificación de la población. Se hizo un análisis factorial de componentes principales de cada uno de los ítems de las escalas de evaluación de síntomas negativos (SANS) y de síntomas positivos (SAPS) para determinar las dimensiones clínicas. Posteriormente, se evaluó la asociación de las variantes genéticas con la esquizofrenia y con cada una de las dimensiones. Resultados: Se encontró asociación entre el genotipo C/C del marcador rs945713 con esquizofrenia (OR = 1,79, IC95%: 1,13-2,84). El genotipo C/C de rs945713 se asoció con puntuaciones más altas en la dimensión “aplanamiento afectivo y alogia” y el genotipo A/A del marcador rs4657181 se relacionó con puntuaciones más bajas en esa misma dimensión. Conclusiones: Se encontró asociación significativa de marcadores dentro de NOS1AP con esquizofrenia y la dimensión clínica “aplanamiento afectivo y alogia”. Estos resultados son consistentes con estudios previos y apoyan la posibilidad de que NOS1AP influya en la susceptibilidad a esquizofrenia y que sea un modificador de sus características clínicas…

Introduction: The nitric oxide synthase 1 adaptor protein (NOS1AP) gene is possibly implicated in schizophrenia etiopathogenesis. Objective: To determine the association of NOS1AP gene variants with schizophrenia and the relationship of variants with the clinical dimensions of the disorder in the Colombian population. Methodology: It is a case-control study with 255 subjects per group. Markers within the NOS1AP gene were typified as well as other informative material of genetic origin so as to adjust by population stratification. A factorial analysis of the main components for each item in the Scales for Evaluating Negative Symptoms (SENS) together with the Scales for Evaluating Positive Symptoms (SEPS) to determine clinical dimensions. Results: Association between the C/C genotype of the rs945713 marker with schizophrenia (OR = 1.79, 95% CI: 1.13 – 2.84) was found. The C/C genotype of the rs945713 was related to higher scores in the “affective flattening and alogia” dimension; and the A/A genotype of the rs4657181 marker was associated to lower scores in the same dimension. Conclusions: Significant associations of markers inside the NOS1AP gene with schizophrenia and the “affective flattening and alogia” clinical dimension were found. These results are consistent with previous studies and support the possibility that NOS1AP influences schizophrenia susceptibility. Furthermore, NOS1AP might be a modifier of schizophrenia clinical characteristics…

[Association of Schizophrenia and its Clinical Implications with the NOS1AP Gene in the Colombian Population]. / Asociación de esquizofrenia y sus dimensiones clínicas con el gen NOS1AP en población colombiana.

INTRODUCTION: The nitric oxide synthase 1 adaptor protein (NOS1AP) gene is possibly implicated in schizophrenia etiopathogenesis. OBJECTIVE: To determine the association of NOS1AP gene variants with schizophrenia and the relationship of variants with the clinical dimensions of the disorder in the Colombian population. METHODOLOGY: It is a case-control study with 255 subjects per group. Markers within the NOS1AP gene were typified as well as other informative material of genetic origin so as to adjust by population stratification. A factorial analysis of the main components for each item in the Scales for Evaluating Negative Symptoms (SENS) together with the Scales for Evaluating Positive Symptoms (SEPS) to determine clinical dimensions. RESULTS: Association between the C/C genotype of the rs945713 marker with schizophrenia (OR = 1.79, 95% CI: 1.13 - 2.84) was found. The C/C genotype of the rs945713 was related to higher scores in the "affective flattening and alogia" dimension; and the A/A genotype of the rs4657181 marker was associated to lower scores in the same dimension. CONCLUSIONS: Significant associations of markers inside the NOS1AP gene with schizophrenia and the "affective flattening and alogia" clinical dimension were found. These results are consistent with previous studies and support the possibility that NOS1AP influences schizophrenia susceptibility. Furthermore, NOS1AP might be a modifier of schizophrenia clinical characteristics.

A gain-of-function mutation in TRPA1 causes familial episodic pain syndrome.

Human monogenic pain syndromes have provided important insights into the molecular mechanisms that underlie normal and pathological pain states. We describe an autosomal-dominant familial episodic pain syndrome characterized by episodes of debilitating upper body pain, triggered by fasting and physical stress. Linkage and haplotype analysis mapped this phenotype to a 25 cM region on chromosome 8q12-8q13. Candidate gene sequencing identified a point mutation (N855S) in the S4 transmembrane segment of TRPA1, a key sensor for environmental irritants. The mutant channel showed a normal pharmacological profile but altered biophysical properties, with a 5-fold increase in inward current on activation at normal resting potentials. Quantitative sensory testing demonstrated normal baseline sensory thresholds but an enhanced secondary hyperalgesia to punctate stimuli on treatment with mustard oil. TRPA1 antagonists inhibit the mutant channel, promising a useful therapy for this disorder. Our findings provide evidence that variation in the TRPA1 gene can alter pain perception in humans.

Association of DRD2 variants and Gilles de la Tourette syndrome in a family-based sample from a South American population isolate.

Gilles de la Tourette Syndrome (GTS) is a chronic neuropsychiatric disorder characterized by motor and vocal tics. Epidemiological evidence supports the importance of genetic factors in disease susceptibility, whereas pharmacological and neuroimaging studies have suggested a defect in the dopamine system. The dopamine receptor D2 gene (DRD2) has been reported to be associated with GTS and related phenotypes. Here, we evaluate genetic association between DRD2 and GTS in a sample from a South American population isolate (Antioquia, Colombia). We genotyped nine single nucleotide polymorphisms (SNPs) across the DRD2 gene region in 69 GTS patients and their nuclear families and carried out both SNP and haplotype-based transmission distortion analysis. Evidence for association was found for three SNPs (rs6279, rs1079597 and rs4648318) and a five marker-haplotype comprising both rs6279 and rs1079597. Our findings replicate the association of DRD2 and GTS, and are consistent with the proposed connection between the dopamine system and this complex neuropsychiatric disease.

Evidence for a role of the NOS1AP (CAPON) gene in schizophrenia and its clinical dimensions: an association study in a South American population isolate.

BACKGROUND/AIMS: Recent studies have implicated a region on chromosome 1q21-23, including the NOS1AP gene, in susceptibility to schizophrenia. However, replication studies have been inconsistent, a fact that could partly relate to the marked psychopathological heterogeneity of schizophrenia. The aim of this study is to evaluate association of polymorphisms in the NOS1AP gene region to schizophrenia, in patients from a South American population isolate, and to assess if these variants are associated with specific clinical dimensions of the disorder. METHODS: We genotyped 24 densely spaced SNPs in the NOS1AP gene region in a schizophrenia trio sample. The transmission disequilibrium test (TDT) was applied to single marker and haplotype data. Association to clinical dimensions (identified by factor analysis) was evaluated using a quantitative transmission disequilibrium test (QTDT). RESULTS: We found significant association between eight SNPs in the NOS1AP gene region to schizophrenia (minimum p value = 0.004). The QTDT analysis of clinical dimensions revealed an association to a dimension consisting mainly of negative symptoms (minimum p value 0.001). CONCLUSIONS: Our findings are consistent with a role for NOS1AP in susceptibility to schizophrenia, especially for the 'negative syndrome' of the disorder.

Convergent linkage evidence from two Latin-American population isolates supports the presence of a susceptibility locus for bipolar disorder in 5q31-34.

We performed a whole genome microsatellite marker scan in six multiplex families with bipolar (BP) mood disorder ascertained in Antioquia, a historically isolated population from North West Colombia. These families were characterized clinically using the approach employed in independent ongoing studies of BP in the closely related population of the Central Valley of Costa Rica. The most consistent linkage results from parametric and non-parametric analyses of the Colombian scan involved markers on 5q31-33, a region implicated by the previous studies of BP in Costa Rica. Because of these concordant results, a follow-up study with additional markers was undertaken in an expanded set of Colombian and Costa Rican families; this provided a genome-wide significant evidence of linkage of BPI to a candidate region of approximately 10 cM in 5q31-33 (maximum non-parametric linkage score=4.395, P<0.00004). Interestingly, this region has been implicated in several previous genetic studies of schizophrenia and psychosis, including disease association with variants of the enthoprotin and gamma-aminobutyric acid receptor genes.

Transmission distortion of BDNF variants to bipolar disorder type I patients from a South American population isolate.

Recent reports have implicated polymorphisms in the brain derived neurotrophic factor (BDNF) gene region in the etiology of several psychiatric phenotypes, including bipolar disorder. Significant disease association has been reported for the G allele at SNP rs6265, which encodes for Valine at position 66 of BDNF (Val66Met), an apparently functional variant of this key BDNF. Here we examined a sample of 224 bipolar type I patients and available parents (comprising a total of 212 nuclear families) ascertained in a South American population isolate (Antioquia, Colombia). We tested for transmission distortion to bipolar patients of alleles at the rs6265 polymorphism and at a microsatellite marker 1.3 kb away from this SNP. Significant excess transmission of the rs6265 G allele to cases was observed (chi(2) = 10.77, d.f. = 1, P = 0.001). Two-locus haplotype analysis showed a significant global transmission distortion (chi(2) = 16.059, d.f. = 7, P = 0.025) with an excess transmission of a haplotype comprising the rs6265 G allele and microsatellite allele 227. These results are consistent with previous studies pointing to a role for BDNF in susceptibility to mood disorders.