lociNGS: A Lightweight Alternative for Assessing Suitability of Next-Generation Loci for Evolutionary Analysis

Genomic enrichment methods and next-generation sequencing produce uneven coverage for the portions of the genome (the loci) they target; this information is essential for ascertaining the suitability of each locus for further analysis. lociNGS is a user-friendly accessory program that takes multi-FASTA formatted loci, next-generation sequence alignments and demographic data as input and collates, displays and outputs information about the data. Summary information includes the parameters coverage per locus, coverage per individual and number of polymorphic sites, among others. The program can output the raw sequences used to call loci from next-generation sequencing data. lociNGS also reformats subsets of loci in three commonly used formats for multi-locus phylogeographic and population genetics analyses – NEXUS, IMa2 and Migrate. lociNGS is available at https://github.com/SHird/lociNGS and is dependent on installation of MongoDB (freely available at http://www.mongodb.org/downloads). lociNGS is written in Python and is supported on MacOSX and Unix; it is distributed under a GNU General Public License.     

MutAid: Sanger and NGS Based Integrated Pipeline for Mutation Identification,Validation and Annotation in Human Molecular Genetics

Traditional Sanger sequencing as well as Next-Generation Sequencing have been used for the identification of disease causing mutations in human molecular research. The majority of currently available tools are developed for research and explorative purposes and often do not provide a complete, efficient, one-stop solution. As the focus of currently developed tools is mainly on NGS data analysis, no integrative solution for the analysis of Sanger data is provided and consequently a one-stop solution to analyze reads from both sequencing platforms is not available. We have therefore developed a new pipeline called MutAid to analyze and interpret raw sequencing data produced by Sanger or several NGS sequencing platforms. It performs format conversion, base calling, quality trimming, filtering, read mapping, variant calling, variant annotation and analysis of Sanger and NGS data under a single platform. It is capable of analyzing reads from multiple patients in a single run to create a list of potential disease causing base substitutions as well as insertions and deletions. MutAid has been developed for expert and non-expert users and supports four sequencing platforms including Sanger, Illumina, 454 and Ion Torrent. Furthermore, for NGS data analysis, five read mappers including BWA, TMAP, Bowtie, Bowtie2 and GSNAP and four variant callers including GATK-HaplotypeCaller, SAMTOOLS, Freebayes and VarScan2 pipelines are supported. MutAid is freely available at https://sourceforge.net/projects/mutaid.     

Personalized Oncology Suite: integrating next-generation sequencing data and whole-slide bioimages

Background

Cancer immunotherapy has recently entered a remarkable renaissance phase with the approval of several agents for treatment. Cancer treatment platforms have demonstrated profound tumor regressions including complete cure in patients with metastatic cancer. Moreover, technological advances in next-generation sequencing (NGS) as well as the development of devices for scanning whole-slide bioimages from tissue sections and image analysis software for quantitation of tumor-infiltrating lymphocytes (TILs) allow, for the first time, the development of personalized cancer immunotherapies that target patient specific mutations. However, there is currently no bioinformatics solution that supports the integration of these heterogeneous datasets.

Results

We have developed a bioinformatics platform – Personalized Oncology Suite (POS) – that integrates clinical data, NGS data and whole-slide bioimages from tissue sections. POS is a web-based platform that is scalable, flexible and expandable. The underlying database is based on a data warehouse schema, which is used to integrate information from different sources. POS stores clinical data, genomic data (SNPs and INDELs identified from NGS analysis), and scanned whole-slide images. It features a genome browser as well as access to several instances of the bioimage management application Bisque. POS provides different visualization techniques and offers sophisticated upload and download possibilities. The modular architecture of POS allows the community to easily modify and extend the application.

Conclusions

The web-based integration of clinical, NGS, and imaging data represents a valuable resource for clinical researchers and future application in medical oncology. POS can be used not only in the context of cancer immunology but also in other studies in which NGS data and images of tissue sections are generated. The application is open-source and can be downloaded at http://www.icbi.at/POS.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2105-15-306) contains supplementary material, which is available to authorized users.     

Efficient In Silico Identification of a Common Insertion in the MAK Gene which Causes Retinitis Pigmentosa

Background

Next generation sequencing (NGS) offers a rapid and comprehensive method of screening for mutations associated with retinitis pigmentosa and related disorders. However, certain sequence alterations such as large insertions or deletions may remain undetected using standard NGS pipelines. One such mutation is a recently-identified Alu insertion into the Male Germ Cell-Associated Kinase (MAK) gene, which is missed by standard NGS-based variant callers. Here, we developed an in silico method of searching NGS raw sequence reads to detect this mutation, without the need to recalculate sequence alignments or to screen every sample by PCR.

Methods

The Linux program grep was used to search for a 23 bp “probe” sequence containing the known junction sequence of the insert. A corresponding search was performed with the wildtype sequence. The matching reads were counted and further compared to the known sequences of the full wildtype and mutant genomic loci. (See https://github.com/MEEIBioinformaticsCenter/grepsearch.)

Results

In a test sample set consisting of eleven previously published homozygous mutants, detection of the MAK-Alu insertion was validated with 100% sensitivity and specificity. As a discovery cohort, raw NGS reads from 1,847 samples (including custom and whole exome selective capture) were searched in ~1 hour on a local computer cluster, yielding an additional five samples with MAK-Alu insertions and solving two previously unsolved pedigrees. Of these, one patient was homozygous for the insertion, one compound heterozygous with a missense change on the other allele (c. 46G>A; p.Gly16Arg), and three were heterozygous carriers.

Conclusions

Using the MAK-Alu grep program proved to be a rapid and effective method of finding a known, disease-causing Alu insertion in a large cohort of patients with NGS data. This simple approach avoids wet-lab assays or computationally expensive algorithms, and could also be used for other known disease-causing insertions and deletions.     

QC-Chain: Fast and Holistic Quality Control Method for Next-Generation Sequencing Data

Next-generation sequencing (NGS) technologies have been widely used in life sciences. However, several kinds of sequencing artifacts, including low-quality reads and contaminating reads, were found to be quite common in raw sequencing data, which compromise downstream analysis. Therefore, quality control (QC) is essential for raw NGS data. However, although a few NGS data quality control tools are publicly available, there are two limitations: First, the processing speed could not cope with the rapid increase of large data volume. Second, with respect to removing the contaminating reads, none of them could identify contaminating sources de novo, and they rely heavily on prior information of the contaminating species, which is usually not available in advance. Here we report QC-Chain, a fast, accurate and holistic NGS data quality-control method. The tool synergeticly comprised of user-friendly tools for (1) quality assessment and trimming of raw reads using Parallel-QC, a fast read processing tool; (2) identification, quantification and filtration of unknown contamination to get high-quality clean reads. It was optimized based on parallel computation, so the processing speed is significantly higher than other QC methods. Experiments on simulated and real NGS data have shown that reads with low sequencing quality could be identified and filtered. Possible contaminating sources could be identified and quantified de novo, accurately and quickly. Comparison between raw reads and processed reads also showed that subsequent analyses (genome assembly, gene prediction, gene annotation, etc.) results based on processed reads improved significantly in completeness and accuracy. As regard to processing speed, QC-Chain achieves 7–8 time speed-up based on parallel computation as compared to traditional methods. Therefore, QC-Chain is a fast and useful quality control tool for read quality process and de novo contamination filtration of NGS reads, which could significantly facilitate downstream analysis. QC-Chain is publicly available at: http://www.computationalbioenergy.org/qc-chain.html.     

Genome Destabilizing Mutator Alleles Drive Specific Mutational Trajectories in Saccharomyces cerevisiae

In addition to environmental factors and intrinsic variations in base substitution rates, specific genome-destabilizing mutations can shape the mutational trajectory of genomes. How specific alleles influence the nature and position of accumulated mutations in a genomic context is largely unknown. Understanding the impact of genome-destabilizing alleles is particularly relevant to cancer genomes where biased mutational signatures are identifiable. We first created a more complete picture of cellular pathways that impact mutation rate using a primary screen to identify essential Saccharomyces cerevisiae gene mutations that cause mutator phenotypes. Drawing primarily on new alleles identified in this resource, we measure the impact of diverse mutator alleles on mutation patterns directly by whole-genome sequencing of 68 mutation-accumulation strains derived from wild-type and 11 parental mutator genotypes. The accumulated mutations differ across mutator strains, displaying base-substitution biases, allele-specific mutation hotspots, and break-associated mutation clustering. For example, in mutants of POLα and the Cdc13–Stn1–Ten1 complex, we find a distinct subtelomeric bias for mutations that we show is independent of the target sequence. Together our data suggest that specific genome-instability mutations are sufficient to drive discrete mutational signatures, some of which share properties with mutation patterns seen in tumors. Thus, in a population of cells, genome-instability mutations could influence clonal evolution by establishing discrete mutational trajectories for genomes.     

Dominant mutations in familial lethal and severe osteogenesis imperfecta

Summary Four families presenting with familial osteogenesis imperfecta (OI) have been studied: 2 with the lethal type II and 2 with the severe type III form. Fibroblasts of the patients, all issue from non-consanguineous parents, produced normal and abnormal (I) chains. These heterozygous mutations differentiate the recurrent forms from homozygous mutations characteristic of autosomal recessive forms. Although the identity of the mutations could not be determined, such recurrence of autosomal dominant OI is probably the result of germinal mosaicism in one of the parents. Biochemical results were consistent with a somatic mosaicism in the father's fibroblasts in one family. Moreover, our studies show that not only OI type II but also severe OI type III can arise from gonadal mosaicism. We discuss the importance of such a phenomenon for genetic counseling.     

Molecular Pathology of Rare Bleeding Disorders (RBDs) in India: A Systematic Review

Background

Though rare in occurrence, patients with rare bleeding disorders (RBDs) are highly heterogeneous and may manifest with severe bleeding diathesis. Due to the high rate of consanguinity in many caste groups, these autosomal recessive bleeding disorders which are of rare occurrence in populations across the world, may not be as rare in India.

Objectives

To comprehensively analyze the frequency and nature of mutations in Indian patients with RBDs.

Methods

Pubmed search was used (www.pubmed.com) to explore the published literature from India on RBDs using the key words “rare bleeding disorders”, “mutations”, “India”, “fibrinogen”, “afibrinogenemia”, “factor II deficiency”, “prothrombin” “factor VII deficiency”, “factor V deficiency”, “factor X deficiency”, “factor XI deficiency”, “combined factor V and VIII deficiency”, “factor XIII deficiency”, “Bernard Soulier syndrome” and “Glanzmanns thrombasthenia” in different combinations. A total of 60 relevant articles could be retrieved. The distribution of mutations from India was compared with that of the world literature by referring to the Human Gene Mutation Database (HGMD) (www.hgmd.org).

Results

Taken together, 181 mutations in 270 patients with different RBDs have been reported from India. Though the types of mutations reported from India and their percentage distribution with respect to the world data are largely similar, yet much higher percentage of small deletions, duplication mutations, insertions, indels were observed in this analysis. Besides the identification of novel mutations and polymorphisms, several common mutations have also been reported, which will allow to develop a strategy for mutation screening in Indian patients with RBDs.

Conclusion

There is a need for a consortium of Institutions working on the molecular pathology of RBDs in India. This will facilitate a quicker and cheaper diagnosis of RBDs besides its utility in first trimester prenatal diagnosis of the affected families.     

Too Much of a Good Thing: The Unique and Repeated Paths Toward Copper Adaptation

Copper is a micronutrient essential for growth due to its role as a cofactor in enzymes involved in respiration, defense against oxidative damage, and iron uptake. Yet too much of a good thing can be lethal, and yeast cells typically do not have tolerance to copper levels much beyond the concentration in their ancestral environment. Here, we report a short-term evolutionary study of Saccharomyces cerevisiae exposed to levels of copper sulfate that are inhibitory to the initial strain. We isolated and identified adaptive mutations soon after they arose, reducing the number of neutral mutations, to determine the first genetic steps that yeast take when adapting to copper. We analyzed 34 such strains through whole-genome sequencing and by assaying fitness within different environments; we also isolated a subset of mutations through tetrad analysis of four lines. We identified a multilayered evolutionary response. In total, 57 single base-pair mutations were identified across the 34 lines. In addition, gene amplification of the copper metallothionein protein, CUP1-1, was rampant, as was chromosomal aneuploidy. Four other genes received multiple, independent mutations in different lines (the vacuolar transporter genes VTC1 and VTC4; the plasma membrane H+-ATPase PMA1; and MAM3, a protein required for normal mitochondrial morphology). Analyses indicated that mutations in all four genes, as well as CUP1-1 copy number, contributed significantly to explaining variation in copper tolerance. Our study thus finds that evolution takes both common and less trodden pathways toward evolving tolerance to an essential, but highly toxic, micronutrient.     

The Abundant Histone Chaperones Spt6 and FACT Collaborate to Assemble,Inspect, and Maintain Chromatin Structure in Saccharomyces cerevisiae

Saccharomyces cerevisiae Spt6 protein is a conserved chromatin factor with several distinct functional domains, including a natively unstructured 30-residue N-terminal region that binds competitively with Spn1 or nucleosomes. To uncover physiological roles of these interactions, we isolated histone mutations that suppress defects caused by weakening Spt6:Spn1 binding with the spt6-F249K mutation. The strongest suppressor was H2A-N39K, which perturbs the point of contact between the two H2A-H2B dimers in an assembled nucleosome. Substantial suppression also was observed when the H2A-H2B interface with H3-H4 was altered, and many members of this class of mutations also suppressed a defect in another essential histone chaperone, FACT. Spt6 is best known as an H3-H4 chaperone, but we found that it binds with similar affinity to H2A-H2B or H3-H4. Like FACT, Spt6 is therefore capable of binding each of the individual components of a nucleosome, but unlike FACT, Spt6 did not produce endonuclease-sensitive reorganized nucleosomes and did not displace H2A-H2B dimers from nucleosomes. Spt6 and FACT therefore have distinct activities, but defects can be suppressed by overlapping histone mutations. We also found that Spt6 and FACT together are nearly as abundant as nucleosomes, with ∼24,000 Spt6 molecules, ∼42,000 FACT molecules, and ∼75,000 nucleosomes per cell. Histone mutations that destabilize interfaces within nucleosomes therefore reveal multiple spatial regions that have both common and distinct roles in the functions of these two essential and abundant histone chaperones. We discuss these observations in terms of different potential roles for chaperones in both promoting the assembly of nucleosomes and monitoring their quality.     

The Mechanism of Nucleotide Excision Repair-Mediated UV-Induced Mutagenesis in Nonproliferating Cells

Following the irradiation of nondividing yeast cells with ultraviolet (UV) light, most induced mutations are inherited by both daughter cells, indicating that complementary changes are introduced into both strands of duplex DNA prior to replication. Early analyses demonstrated that such two-strand mutations depend on functional nucleotide excision repair (NER), but the molecular mechanism of this unique type of mutagenesis has not been further explored. In the experiments reported here, an ade2 adeX colony-color system was used to examine the genetic control of UV-induced mutagenesis in nondividing cultures of Saccharomyces cerevisiae. We confirmed a strong suppression of two-strand mutagenesis in NER-deficient backgrounds and demonstrated that neither mismatch repair nor interstrand crosslink repair affects the production of these mutations. By contrast, proteins involved in the error-prone bypass of DNA damage (Rev3, Rev1, PCNA, Rad18, Pol32, and Rad5) and in the early steps of the DNA-damage checkpoint response (Rad17, Mec3, Ddc1, Mec1, and Rad9) were required for the production of two-strand mutations. There was no involvement, however, for the Pol η translesion synthesis DNA polymerase, the Mms2-Ubc13 postreplication repair complex, downstream DNA-damage checkpoint factors (Rad53, Chk1, and Dun1), or the Exo1 exonuclease. Our data support models in which UV-induced mutagenesis in nondividing cells occurs during the Pol ζ-dependent filling of lesion-containing, NER-generated gaps. The requirement for specific DNA-damage checkpoint proteins suggests roles in recruiting and/or activating factors required to fill such gaps.     

Resection Activity of the Sgs1 Helicase Alters the Affinity of DNA Ends for Homologous Recombination Proteins in Saccharomyces cerevisiae

The RecQ helicase family is critical during DNA damage repair, and mutations in these proteins are associated with Bloom, Werner, or Rothmund-Thompson syndromes in humans, leading to cancer predisposition and/or premature aging. In the budding yeast Saccharomyces cerevisiae, mutations in the RecQ homolog, SGS1, phenocopy many of the defects observed in the human syndromes. One challenge to studying RecQ helicases is that their disruption leads to a pleiotropic phenotype. Using yeast, we show that the separation-of-function allele of SGS1, sgs1-D664Δ, has impaired activity at DNA ends, resulting in a resection processivity defect. Compromising Sgs1 resection function in the absence of the Sae2 nuclease causes slow growth, which is alleviated by making the DNA ends accessible to Exo1 nuclease. Furthermore, fluorescent microscopy studies reveal that, when Sgs1 resection activity is compromised in sae2Δ cells, Mre11 repair foci persist. We suggest a model where the role of Sgs1 in end resection along with Sae2 is important for removing Mre11 from DNA ends during repair.     

Novel compound heterozygous mutations in SLC5A2 are responsible for autosomal recessive renal glucosuria

Familial renal glucosuria is an inherited renal tubular disorder. A homozygous nonsense mutation in the SLC5A2 gene, encoding the sodium/glucose co-transporter SGLT2, has recently been identified in an affected child of consanguineous parents. We now report novel compound heterozygous mutations in the son of non-consanguineous parents. One allele has a p.Q167fsX186 mutation, which is expected to produce a truncated protein, and the other a p.N654S mutation involving a highly conserved residue. These findings confirm that mutations in the SLC5A2 gene are responsible for recessive renal glucosuria.     

A database of natural products and chemical entities from marine habitat

Marine compound database consists of marine natural products and chemical entities, collected from various literature sources, which are known to possess bioactivity against human diseases. The database is constructed using html code. The 12 categories of 182 compounds are provided with the source, compound name, 2-dimensional structure, bioactivity and clinical trial information. The database is freely available online and can be accessed at http://www.progenebio.in/mcdb/index.htm     

DNA Replication Error-Induced Extinction of Diploid Yeast

Genetic defects in DNA polymerase accuracy, proofreading, or mismatch repair (MMR) induce mutator phenotypes that accelerate adaptation of microbes and tumor cells. Certain combinations of mutator alleles synergistically increase mutation rates to levels that drive extinction of haploid cells. The maximum tolerated mutation rate of diploid cells is unknown. Here, we define the threshold for replication error-induced extinction (EEX) of diploid Saccharomyces cerevisiae. Double-mutant pol3 alleles that carry mutations for defective DNA polymerase-δ proofreading (pol3-01) and accuracy (pol3-L612M or pol3-L612G) induce strong mutator phenotypes in heterozygous diploids (POL3/pol3-01,L612M or POL3/pol3-01,L612G). Both pol3-01,L612M and pol3-01,L612G alleles are lethal in the homozygous state; cells with pol3-01,L612M divide up to 10 times before arresting at random stages in the cell cycle. Antimutator eex mutations in the pol3 alleles suppress this lethality (pol3-01,L612M,eex or pol3-01,L612G,eex). MMR defects synergize with pol3-01,L612M,eex and pol3-01,L612G,eex alleles, increasing mutation rates and impairing growth. Conversely, inactivation of the Dun1 S-phase checkpoint kinase suppresses strong pol3-01,L612M,eex and pol3-01,L612G,eex mutator phenotypes as well as the lethal pol3-01,L612M phenotype. Our results reveal that the lethal error threshold in diploids is 10 times higher than in haploids and likely determined by homozygous inactivation of essential genes. Pronounced loss of fitness occurs at mutation rates well below the lethal threshold, suggesting that mutator-driven cancers may be susceptible to drugs that exacerbate replication errors.