LoFreq: a sequence-quality aware, ultra-sensitive variant caller for uncovering cell-population heterogeneity from high-throughput sequencing datasets.

The study of cell-population heterogeneity in a range of biological systems, from viruses to bacterial isolates to tumor samples, has been transformed by recent advances in sequencing throughput. While the high-coverage afforded can be used, in principle, to identify very rare variants in a population, existing ad hoc approaches frequently fail to distinguish true variants from sequencing errors. We report a method (LoFreq) that models sequencing run-specific error rates to accurately call variants occurring in <0.05% of a population. Using simulated and real datasets (viral, bacterial and human), we show that LoFreq has near-perfect specificity, with significantly improved sensitivity compared with existing methods and can efficiently analyze deep Illumina sequencing datasets without resorting to approximations or heuristics. We also present experimental validation for LoFreq on two different platforms (Fluidigm and Sequenom) and its application to call rare somatic variants from exome sequencing datasets for gastric cancer. Source code and executables for LoFreq are freely available at http://sourceforge.net/projects/lofreq/.

Genetic variability of RyR2 and CASQ2 genes in an Asian population.

We analyzed the coding regions of the cardiac calcium-handling genes, ryanodine receptor 2 (RyR2) and calsequestrin 2 (CASQ2) for genetic variants in a healthy Chinese population (n=95) and in a cohort of 28 sudden unexplained death victims. Mutations in RyR2 and CASQ2 have been shown to alter calcium homeostasis during excitation-contraction coupling and predispose individuals to fatal cardiac arrhythmias. The genetic screening was accomplished by denaturing high-performance liquid chromatography and DNA sequencing methods. Genetic analysis revealed the following non-synonymous genetic variations: two reported RyR2 polymorphisms; 5654G>A (G1885E) and 5656G>A (G1886S), two reported CASQ2 polymorphisms; 196A>G (T66A) and 226G>A (V76M) and one novel CASQ2 mutation; 529G>C (E177Q). The functional significance of the novel CASQ2 mutation has not been evaluated and characterized. This study shows that multiple genetic variations of the RyR2 and CASQ2 genes exist in the two study populations. The inter-individual genetic variability may underlie the different susceptibility of individuals to developing ventricular tachycardia. The research results will be valuable for which future work involving clinical and forensic samples can be based upon to distinguish potential disease-associated mutations from common polymorphisms.

Induction of human sulfotransferase 1A3 (SULT1A3) by glucocorticoids.

Sulfotransferases (SULTs) play an important role in the detoxification and bioactivation of endogenous compounds and xenobiotics. Studies on rat sulfotransferases had shown that SULT genes, like cytochrome P450 genes, can be regulated by ligands that bind nuclear receptors. For human SULT genes, the regulation of human SULT2A1 expression is currently the best characterized. In this study, we systematically examined the regulation of human SULT1A genes by glucocorticoids. Treatment of the human hepatocellular carcinoma derived HepG2 cells with 10(-7) M dexamethasone did not affect the SULT1A1 activity toward p-nitrophenol. In contrast, SULT1A3 activity toward dopamine was significantly induced. Transient transfection of the SULT1A3 5'-flanking region/luciferase reporter construct showed that SULT1A3 was responsive to dexamethasone and prednisolone in a concentration-dependent manner with maximal induction at 10(-7) M dexamethasone or 1 microM prednisolone. In addition, induction by dexamethasone was dependent on the level of expression of the glucocorticoid receptor. Analysis of the 5'-flanking region led to the identification of a putative glucocorticoid response element at position (-1211 to -1193) upstream of the transcription start site and deletion or mutation of this element resulted in a loss of response. In summary, the data from this study shows that the human SULT1A3 gene is inducible by glucocorticoids through a glucocorticoid receptor-mediated mechanism and the glucocorticoid response element at position (-1211 to -1193) is necessary for this induction.