New stopping criteria for segmenting DNA sequences

We propose a solution on the stopping criterion in segmenting inhomogeneous DNA sequences with complex statistical patterns. This new stopping criterion is based on Bayesian Information Criterion (BIC) in the model selection framework. When this stopping criterion is applied to a left telomere sequence of yeast Saccharomyces cerevisiae and the complete genome sequence of bacterium Escherichia coli, borders of biologically meaningful units were identified (e.g. subtelomeric units, replication origin, and replication terminus), and a more reasonable number of domains was obtained. We also introduce a measure called segmentation strength which can be used to control the delineation of large domains. The relationship between the average domain size and the threshold of segmentation strength is determined for several genome sequences.Comment: 4 pages, 4 figures, Physical Review Letters, to appea

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

Simplifying the mosaic description of DNA sequences

By using the Jensen-Shannon divergence, genomic DNA can be divided into compositionally distinct domains through a standard recursive segmentation procedure. Each domain, while significantly different from its neighbours, may however share compositional similarity with one or more distant (non--neighbouring) domains. We thus obtain a coarse--grained description of the given DNA string in terms of a smaller set of distinct domain labels. This yields a minimal domain description of a given DNA sequence, significantly reducing its organizational complexity. This procedure gives a new means of evaluating genomic complexity as one examines organisms ranging from bacteria to human. The mosaic organization of DNA sequences could have originated from the insertion of fragments of one genome (the parasite) inside another (the host), and we present numerical experiments that are suggestive of this scenario.Comment: 16 pages, 1 figure, Accepted for publication in Phys. Rev.

A Yeast Recombination Assay to Characterize Human BRCA1 Missense Variants of Unknown Pathological Significance

The BRCA1 tumor suppressor gene is found mutated in familial breast cancer. Although many of the mutations are clearly pathological because they give rise to truncated proteins, several missense variants of uncertain pathological consequences have been identified. A novel functional assay to screen for BRCA1 missense variants in a simple genetic system could be very useful for the identification of potentially deleterious mutations. By using two prediction computer programs, Sorting Intolerant from Tolerant (SIFT) and Polymorphism Phenotyping (PolyPhen), seven nonsynonymous missense BRCA1 variants likely disrupting the gene function were selected as potentially deleterious. The budding yeast Saccharomyces cerevisiae (S. cerevisiae) was used to test these cancer-related missense mutations for their ability to affect cell growth and homologous recombination (HR) at the HIS3 and ADE2 loci. The variants localized in the BRCA1 C-Terminus (BRCT) domain did not show any growth inhibition when overexpressed in agreement with previous results. Overexpression of either wild-type BRCA1 or two neutral missense variants did not increase yeast HR but when cancer-related variants were overexpressed a significant increase in recombination was observed. Results clearly showed that this genetic system can be useful to discriminate between neutral and deleterious BRCA1 missense variants