BCFtools/liftover: an accurate and comprehensive tool to convert genetic variants across genome assemblies.

MOTIVATION: Many genetics studies report results tied to genomic coordinates of a legacy genome assembly. However, as assemblies are updated and improved, researchers are faced with either realigning raw sequence data using the updated coordinate system or converting legacy datasets to the updated coordinate system to be able to combine results with newer datasets. Currently available tools to perform the conversion of genetic variants have numerous shortcomings, including poor support for indels and multi-allelic variants, that lead to a higher rate of variants being dropped or incorrectly converted. As a result, many researchers continue to work with and publish using legacy genomic coordinates. RESULTS: Here we present BCFtools/liftover, a tool to convert genomic coordinates across genome assemblies for variants encoded in the variant call format with improved support for indels represented by different reference alleles across genome assemblies and full support for multi-allelic variants. It further supports variant annotation fields updates whenever the reference allele changes across genome assemblies. The tool has the lowest rate of variants being dropped with an order of magnitude less indels dropped or incorrectly converted and is an order of magnitude faster than other tools typically used for the same task. It is particularly suited for converting variant callsets from large cohorts to novel telomere-to-telomere assemblies as well as summary statistics from genome-wide association studies tied to legacy genome assemblies. AVAILABILITY AND IMPLEMENTATION: The tool is written in C and freely available under the MIT open source license as a BCFtools plugin available at http://github.com/freeseek/score.

Connexin 50 gene on human chromosome 1q21 is associated with schizophrenia in matched case control and family-based studies.

BACKGROUND: The gap junction subunit connexin permits direct intercellular exchange of ions and molecules including glutamate, and plays an important role in the central nervous system. The connexin 40 (Cx40) and connexin 50 (Cx50) genes are located on chromosome 1q21.1, a region strongly linked with schizophrenia. These lines of evidence suggest that Cx40 and Cx50 may play a role in schizophrenia. METHODS: Using an allele-specific PCR assay, four polymorphisms each were genotyped for Cx40 and Cx50 in 190 Caucasian patients with schizophrenia and 190 controls matched for sex, age and ethnicity. Following up, Cx50 rs989192 and rs4950495 were investigated in 99 Canadian and 163 Portuguese trios and nuclear families with schizophrenia probands. Hardy-Weinberg equilibrium and linkage disequilibrium (LD) block identification was carried out with HaploView, and association analysis for alleles and haplotypes with a permutation test of 10 000 simulations was carried out using the UNPHASED software program. RESULTS: Distributions of genotype frequencies of all markers were in Hardy-Weinberg equilibrium in Caucasian patients, controls and families. One rs989192-rs4950495 LD block was found in patients but not in controls. We found a significant association between the Cx50 rs989192-rs4950495 haplotype and schizophreniay (chi(2) = 29.55, p<0.01). The A-C haplotype had a higher frequency in patients (chi(2) = 7.153, p<0.01). Family studies also showed that the A-C haplotype was transmitted more often to patients with schizophrenia (chi(2) = 8.43, p<0.01). No association of Cx40 with schizophrenia was found for allele, genotype or haplotype analyses. CONCLUSIONS: Our matched case-control and family study indicate that Cx50, but not Cx40, may play a role in the genetic susceptibility to schizophrenia.

Human p53 tumor suppressor gene (TP53) and schizophrenia: case-control and family studies.

The human p53 tumor suppressor gene (TP53) is considered as a candidate susceptibility gene for schizophrenia because of its functions in neurodevelopment. To test for an association between TP53 and schizophrenia, both the case-control study and the transmission disequilibrium test (TDT) were performed on genotype data from eight polymorphisms in TP53. Our samples included 286 Toronto schizophrenia cases and 264 controls, and 163 Portuguese nuclear families. In the Toronto case-control study significant differences of allele frequencies of the CAA Ins/Del (p=0.027) and the 16bp Ins/Del (p=0.022) were detected. In TDT analysis we found significant differences for transmission of the CAA Ins/Del (p=0.017) in Portuguese schizophrenia families. Haplotype analysis also showed a significant association between TP53 and schizophrenia. These results provide further evidence that TP53 may play a role in the pathogenesis of schizophrenia.