Somatic Point Mutation Calling in Low Cellularity Tumors

Somatic mutation calling from next-generation sequencing data remains a challenge due to the difficulties of distinguishing true somatic events from artifacts arising from PCR, sequencing errors or mis-mapping. Tumor cellularity or purity, sub-clonality and copy number changes also confound the identification of true somatic events against a background of germline variants. We have developed a heuristic strategy and software (http://www.qcmg.org/bioinformatics/qsnp/) for somatic mutation calling in samples with low tumor content and we show the superior sensitivity and precision of our approach using a previously sequenced cell line, a series of tumor/normal admixtures, and 3,253 putative somatic SNVs verified on an orthogonal platform.     

Remarkably Divergent Regions Punctuate the Genome Assembly of the Caenorhabditis elegans Hawaiian Strain CB4856

The Hawaiian strain (CB4856) of Caenorhabditis elegans is one of the most divergent from the canonical laboratory strain N2 and has been widely used in developmental, population, and evolutionary studies. To enhance the utility of the strain, we have generated a draft sequence of the CB4856 genome, exploiting a variety of resources and strategies. When compared against the N2 reference, the CB4856 genome has 327,050 single nucleotide variants (SNVs) and 79,529 insertion–deletion events that result in a total of 3.3 Mb of N2 sequence missing from CB4856 and 1.4 Mb of sequence present in CB4856 but not present in N2. As previously reported, the density of SNVs varies along the chromosomes, with the arms of chromosomes showing greater average variation than the centers. In addition, we find 61 regions totaling 2.8 Mb, distributed across all six chromosomes, which have a greatly elevated SNV density, ranging from 2 to 16% SNVs. A survey of other wild isolates show that the two alternative haplotypes for each region are widely distributed, suggesting they have been maintained by balancing selection over long evolutionary times. These divergent regions contain an abundance of genes from large rapidly evolving families encoding F-box, MATH, BATH, seven-transmembrane G-coupled receptors, and nuclear hormone receptors, suggesting that they provide selective advantages in natural environments. The draft sequence makes available a comprehensive catalog of sequence differences between the CB4856 and N2 strains that will facilitate the molecular dissection of their phenotypic differences. Our work also emphasizes the importance of going beyond simple alignment of reads to a reference genome when assessing differences between genomes.     

SomatiCA: Identifying,Characterizing and Quantifying Somatic Copy Number Aberrations from Cancer Genome Sequencing Data

Whole genome sequencing of matched tumor-normal sample pairs is becoming routine in cancer research. However, analysis of somatic copy-number changes from sequencing data is still challenging because of insufficient sequencing coverage, unknown tumor sample purity and subclonal heterogeneity. Here we describe a computational framework, named SomatiCA, which explicitly accounts for tumor purity and subclonality in the analysis of somatic copy-number profiles. Taking read depths (RD) and lesser allele frequencies (LAF) as input, SomatiCA will output 1) admixture rate for each tumor sample, 2) somatic allelic copy-number for each genomic segment, 3) fraction of tumor cells with subclonal change in each somatic copy number aberration (SCNA), and 4) a list of substantial genomic aberration events including gain, loss and LOH. SomatiCA is available as a Bioconductor R package at http://www.bioconductor.org/packages/2.13/bioc/html/SomatiCA.html.     

The DNA Damage Response and Checkpoint Adaptation in Saccharomyces cerevisiae: Distinct Roles for the Replication Protein A2 (Rfa2) N-Terminus

In response to DNA damage, two general but fundamental processes occur in the cell: (1) a DNA lesion is recognized and repaired, and (2) concomitantly, the cell halts the cell cycle to provide a window of opportunity for repair to occur. An essential factor for a proper DNA-damage response is the heterotrimeric protein complex Replication Protein A (RPA). Of particular interest is hyperphosphorylation of the 32-kDa subunit, called RPA2, on its serine/threonine-rich amino (N) terminus following DNA damage in human cells. The unstructured N-terminus is often referred to as the phosphorylation domain and is conserved among eukaryotic RPA2 subunits, including Rfa2 in Saccharomyces cerevisiae. An aspartic acid/alanine-scanning and genetic interaction approach was utilized to delineate the importance of this domain in budding yeast. It was determined that the Rfa2 N-terminus is important for a proper DNA-damage response in yeast, although its phosphorylation is not required. Subregions of the Rfa2 N-terminus important for the DNA-damage response were also identified. Finally, an Rfa2 N-terminal hyperphosphorylation-mimetic mutant behaves similarly to another Rfa1 mutant (rfa1-t11) with respect to genetic interactions, DNA-damage sensitivity, and checkpoint adaptation. Our data indicate that post-translational modification of the Rfa2 N-terminus is not required for cells to deal with “repairable” DNA damage; however, post-translational modification of this domain might influence whether cells proceed into M-phase in the continued presence of unrepaired DNA lesions as a “last-resort” mechanism for cell survival.     

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.     

The Importance of cGMP Signaling in Sensory Cilia for Body Size Regulation in Caenorhabditis elegans

The body size of Caenorhabditis elegans is thought to be controlled by sensory inputs because many mutants with sensory cilium structure defects exhibit small body size. The EGL-4 cGMP-dependent protein kinase acts in sensory neurons to reduce body size when animals fail to perceive sensory signals. In addition to body size control, EGL-4 regulates various other behavioral and developmental pathways, including those involved in the regulation of egg laying and chemotaxis behavior. Here we have identified gcy-12, which encodes a receptor-type guanylyl cyclase, as a gene involved in the sensory regulation of body size. Analyses with GFP fusion constructs showed that gcy-12 is expressed in several sensory neurons and localizes to sensory cilia. Genetic analyses indicated that GCY-12 acts upstream of EGL-4 in body size control but does not affect other EGL-4 functions. Our studies indicate that the function of the GCY-12 guanylyl cyclase is to provide cGMP to the EGL-4 cGMP-dependent kinase only for limited tasks including body size regulation. We also found that the PDE-2 cyclic nucleotide phosphodiesterase negatively regulates EGL-4 in controlling body size. Thus, the cGMP level is precisely controlled by GCY-12 and PDE-2 to determine body size through EGL-4, and the defects in the sensory cilium structure may disturb the balanced control of the cGMP level. The large number of guanylyl cyclases encoded in the C. elegans genome suggests that EGL-4 exerts pleiotropic effects by partnering with different guanylyl cyclases for different downstream functions.     

Evidence for a Nonendosomal Function of the Saccharomyces cerevisiae ESCRT-III-Like Protein Chm7

Endosomal sorting complex required for transport (ESCRT) proteins are involved in a number of cellular processes, such as endosomal protein sorting, HIV budding, cytokinesis, plasma membrane repair, and resealing of the nuclear envelope during mitosis. Here we explored the function of a noncanonical member of the ESCRT-III protein family, the Saccharomyces cerevisiae ortholog of human CHMP7. Very little is known about this protein. In silico analysis predicted that Chm7 (yeast ORF YJL049w) is a fusion of an ESCRT-II and ESCRT-III-like domain, which would suggest a role in endosomal protein sorting. However, our data argue against a role of Chm7 in endosomal protein sorting. The turnover of the endocytic cargo protein Ste6 and the vacuolar protein sorting of carboxypeptidase S (CPS) were not affected by CHM7 deletion, and Chm7 also responded very differently to a loss in Vps4 function compared to a canonical ESCRT-III protein. Our data indicate that the Chm7 function could be connected to the endoplasmic reticulum (ER). In line with a function at the ER, we observed a strong negative genetic interaction between the deletion of a gene function (APQ12) implicated in nuclear pore complex assembly and messenger RNA (mRNA) export and the CHM7 deletion. The patterns of genetic interactions between the APQ12 deletion and deletions of ESCRT-III genes, two-hybrid interactions, and the specific localization of mCherry fusion proteins are consistent with the notion that Chm7 performs a novel function at the ER as part of an alternative ESCRT-III complex.     

Neuronal Development in Caenorhabditis elegans Is Regulated by Inhibition of an MLK MAP Kinase Pathway

We show that loss-of-function mutations in kinases of the MLK-1 pathway (mlk-1, mek-1, and kgb-1/jnk) function cell-autonomously in neurons to suppress defects in synapse formation and axon termination caused by rpm-1 loss of function. Our genetic analysis also suggests that the phosphatase PPM-1, like RPM-1, is a potential inhibitor of kinases in the MLK-1 pathway.     

The Principal Role of Ku in Telomere Length Maintenance Is Promotion of Est1 Association with Telomeres

Telomere length is tightly regulated in cells that express telomerase. The Saccharomyces cerevisiae Ku heterodimer, a DNA end-binding complex, positively regulates telomere length in a telomerase-dependent manner. Ku associates with the telomerase RNA subunit TLC1, and this association is required for TLC1 nuclear retention. Ku–TLC1 interaction also impacts the cell-cycle-regulated association of the telomerase catalytic subunit Est2 to telomeres. The promotion of TLC1 nuclear localization and Est2 recruitment have been proposed to be the principal role of Ku in telomere length maintenance, but neither model has been directly tested. Here we study the impact of forced recruitment of Est2 to telomeres on telomere length in the absence of Ku’s ability to bind TLC1 or DNA ends. We show that tethering Est2 to telomeres does not promote efficient telomere elongation in the absence of Ku–TLC1 interaction or DNA end binding. Moreover, restoration of TLC1 nuclear localization, even when combined with Est2 recruitment, does not bypass the role of Ku. In contrast, forced recruitment of Est1, which has roles in telomerase recruitment and activation, to telomeres promotes efficient and progressive telomere elongation in the absence of Ku–TLC1 interaction, Ku DNA end binding, or Ku altogether. Ku associates with Est1 and Est2 in a TLC1-dependent manner and enhances Est1 recruitment to telomeres independently of Est2. Together, our results unexpectedly demonstrate that the principal role of Ku in telomere length maintenance is to promote the association of Est1 with telomeres, which may in turn allow for efficient recruitment and activation of the telomerase holoenzyme.     

New Insights into the Post-Translational Regulation of DNA Damage Response and Double-Strand Break Repair in Caenorhabditis elegans

Although a growing number of studies have reported the importance of SUMOylation in genome maintenance and DNA double-strand break repair (DSBR), relevant target proteins and how this modification regulates their functions are yet to be clarified. Here, we analyzed SUMOylation of ZTF-8, the homolog of mammalian RHINO, to test the functional significance of this protein modification in the DSBR and DNA damage response (DDR) pathways in the Caenorhabditis elegans germline. We found that ZTF-8 is a direct target for SUMOylation in vivo and that its modification is required for DNA damage checkpoint induced apoptosis and DSBR. Non-SUMOylatable mutants of ZTF-8 mimic the phenotypes observed in ztf-8 null mutants, including reduced fertility, impaired DNA damage repair, and defective DNA damage checkpoint activation. However, while mutants for components acting in the SUMOylation pathway fail to properly localize ZTF-8, its localization is not altered in the ZTF-8 non-SUMOylatable mutants. Taken together, these data show that direct SUMOylation of ZTF-8 is required for its function in DSBR as well as DDR but not its localization. ZTF-8’s human ortholog is enriched in the germline, but its meiotic role as well as its post-translational modification has never been explored. Therefore, our discovery may assist in understanding the regulatory mechanism of this protein in DSBR and DDR in the germline.     

Analysis of PHA-1 Reveals a Limited Role in Pharyngeal Development and Novel Functions in Other Tissues

PHA-1 encodes a cytoplasmic protein that is required for embryonic morphogenesis and attachment of the foregut (pharynx) to the mouth (buccal capsule). Previous reports have in some cases suggested that PHA-1 is essential for the differentiation of most or all pharyngeal cell types. By performing mosaic analysis with a recently acquired pha-1 null mutation (tm3671), we found that PHA-1 is not required within most or all pharyngeal cells for their proper specification, differentiation, or function. Rather, our evidence suggests that PHA-1 acts in the arcade or anterior epithelial cells of the pharynx to promote attachment of the pharynx to the future buccal capsule. In addition, PHA-1 appears to be required in the epidermis for embryonic morphogenesis, in the excretory system for osmoregulation, and in the somatic gonad for normal ovulation and fertility. PHA-1 activity is also required within at least a subset of intestinal cells for viability. To better understand the role of PHA-1 in the epidermis, we analyzed several apical junction markers in pha-1(tm3671) homozygous embryos. PHA-1 regulates the expression of several components of two apical junction complexes including AJM-1–DLG-1/discs large complex and the classical cadherin–catenin complex, which may account for the role of PHA-1 in embryonic morphogenesis.     

Linking Gene Expression in the Intestine to Production of Gametes Through the Phosphate Transporter PITR-1 in Caenorhabditis elegans

Inorganic phosphate is an essential mineral for both prokaryotic and eukaryotic cell metabolism and structure. Its uptake into the cell is mediated by membrane-bound transporters and coupled to Na⁺ transport. Mammalian sodium-dependent Pi cotransporters have been grouped into three families NaPi-I, NaPi-II, and NaPi-III. Despite being discovered more than two decades ago, very little is known about requirements for NaPi-III transporters in vivo, in the context of intact animal models. Here we find that impaired function of the Caenorhabditis elegans NaPi-III transporter, pitr-1, results in decreased brood size and dramatically increased expression of vitellogenin by the worm intestine. Unexpectedly, we found that the effects of pitr-1 mutation on vitellogenin expression in the intestine could only be rescued by expression of pitr-1 in the germline, and not by expression of pitr-1 in the intestine itself. Our results indicate the existence of a signal from the germline that regulates gene expression in the intestine, perhaps linking nutrient export from the intestine to production of gametes by the germline.     

The Role of Dbf4-Dependent Protein Kinase in DNA Polymerase ζ-Dependent Mutagenesis in Saccharomyces cerevisiae

The yeast Dbf4-dependent kinase (DDK) (composed of Dbf4 and Cdc7 subunits) is an essential, conserved Ser/Thr protein kinase that regulates multiple processes in the cell, including DNA replication, recombination and induced mutagenesis. Only DDK substrates important for replication and recombination have been identified. Consequently, the mechanism by which DDK regulates mutagenesis is unknown. The yeast mcm5-bob1 mutation that bypasses DDK’s essential role in DNA replication was used here to examine whether loss of DDK affects spontaneous as well as induced mutagenesis. Using the sensitive lys2ΔA746 frameshift reversion assay, we show DDK is required to generate “complex” spontaneous mutations, which are a hallmark of the Polζ translesion synthesis DNA polymerase. DDK co-immunoprecipitated with the Rev7 regulatory, but not with the Rev3 polymerase subunit of Polζ. Conversely, Rev7 bound mainly to the Cdc7 kinase subunit and not to Dbf4. The Rev7 subunit of Polζ may be regulated by DDK phosphorylation as immunoprecipitates of yeast Cdc7 and also recombinant Xenopus DDK phosphorylated GST-Rev7 in vitro. In addition to promoting Polζ-dependent mutagenesis, DDK was also important for generating Polζ-independent large deletions that revert the lys2ΔA746 allele. The decrease in large deletions observed in the absence of DDK likely results from an increase in the rate of replication fork restart after an encounter with spontaneous DNA damage. Finally, nonepistatic, additive/synergistic UV sensitivity was observed in cdc7Δ pol32Δ and cdc7Δ pol30-K127R,K164R double mutants, suggesting that DDK may regulate Rev7 protein during postreplication “gap filling” rather than during “polymerase switching” by ubiquitinated and sumoylated modified Pol30 (PCNA) and Pol32.     

DNA Replication Checkpoint Signaling Depends on a Rad53–Dbf4 N-Terminal Interaction in Saccharomyces cerevisiae

Dbf4-dependent kinase (DDK) and cyclin-dependent kinase (CDK) are essential to initiate DNA replication at individual origins. During replication stress, the S-phase checkpoint inhibits the DDK- and CDK-dependent activation of late replication origins. Rad53 kinase is a central effector of the replication checkpoint and both binds to and phosphorylates Dbf4 to prevent late-origin firing. The molecular basis for the Rad53–Dbf4 physical interaction is not clear but occurs through the Dbf4 N terminus. Here we found that both Rad53 FHA1 and FHA2 domains, which specifically recognize phospho-threonine (pT), interacted with Dbf4 through an N-terminal sequence and an adjacent BRCT domain. Purified Rad53 FHA1 domain (but not FHA2) bound to a pT Dbf4 peptide in vitro, suggesting a possible phospho-threonine-dependent interaction between FHA1 and Dbf4. The Dbf4–Rad53 interaction is governed by multiple contacts that are separable from the Cdc5- and Msa1-binding sites in the Dbf4 N terminus. Importantly, abrogation of the Rad53–Dbf4 physical interaction blocked Dbf4 phosphorylation and allowed late-origin firing during replication checkpoint activation. This indicated that Rad53 must stably bind to Dbf4 to regulate its activity.     

Chromatin Remodeling Factors Isw2 and Ino80 Regulate Checkpoint Activity and Chromatin Structure in S Phase

When cells undergo replication stress, proper checkpoint activation and deactivation are critical for genomic stability and cell survival and therefore must be highly regulated. Although mechanisms of checkpoint activation are well studied, mechanisms of checkpoint deactivation are far less understood. Previously, we reported that chromatin remodeling factors Isw2 and Ino80 attenuate the S-phase checkpoint activity in Saccharomyces cerevisiae, especially during recovery from hydroxyurea. In this study, we found that Isw2 and Ino80 have a more pronounced role in attenuating checkpoint activity during late S phase in the presence of methyl methanesulfonate (MMS). We therefore screened for checkpoint factors required for Isw2 and Ino80 checkpoint attenuation in the presence of MMS. Here we demonstrate that Isw2 and Ino80 antagonize checkpoint activators and attenuate checkpoint activity in S phase in MMS either through a currently unknown pathway or through RPA. Unexpectedly, we found that Isw2 and Ino80 increase chromatin accessibility around replicating regions in the presence of MMS through a novel mechanism. Furthermore, through growth assays, we provide additional evidence that Isw2 and Ino80 partially counteract checkpoint activators specifically in the presence of MMS. Based on these results, we propose that Isw2 and Ino80 attenuate S-phase checkpoint activity through a novel mechanism.     

A γ-Secretase Independent Role for Presenilin in Calcium Homeostasis Impacts Mitochondrial Function and Morphology in Caenorhabditis elegans

Mutations in the presenilin (PSEN) encoding genes (PSEN1 and PSEN2) occur in most early onset familial Alzheimer’s Disease. Despite the identification of the involvement of PSEN in Alzheimer’s Disease (AD) ∼20 years ago, the underlying role of PSEN in AD is not fully understood. To gain insight into the biological function of PSEN, we investigated the role of the PSEN homolog SEL-12 in Caenorhabditis elegans. Using genetic, cell biological, and pharmacological approaches, we demonstrate that mutations in sel-12 result in defects in calcium homeostasis, leading to mitochondrial dysfunction. Moreover, consistent with mammalian PSEN, we provide evidence that SEL-12 has a critical role in mediating endoplasmic reticulum (ER) calcium release. Furthermore, we found that in SEL-12-deficient animals, calcium transfer from the ER to the mitochondria leads to fragmentation of the mitochondria and mitochondrial dysfunction. Additionally, we show that the impact that SEL-12 has on mitochondrial function is independent of its role in Notch signaling, γ-secretase proteolytic activity, and amyloid plaques. Our results reveal a critical role for PSEN in mediating mitochondrial function by regulating calcium transfer from the ER to the mitochondria.     

Structural Integrity of Centromeric Chromatin and Faithful Chromosome Segregation Requires Pat1

The kinetochore (centromeric DNA and associated protein complex) is essential for faithful chromosome segregation and maintenance of genome stability. Here we report that an evolutionarily conserved protein Pat1 is a structural component of Saccharomyces cerevisiae kinetochore and associates with centromeres in a NDC10-dependent manner. Consistent with a role for Pat1 in kinetochore structure and function, a deletion of PAT1 results in delay in sister chromatid separation, errors in chromosome segregation, and defects in structural integrity of centromeric chromatin. Pat1 is involved in topological regulation of minichromosomes as altered patterns of DNA supercoiling were observed in pat1Δ cells. Studies with pat1 alleles uncovered an evolutionarily conserved region within the central domain of Pat1 that is required for its association with centromeres, sister chromatid separation, and faithful chromosome segregation. Taken together, our data have uncovered a novel role for Pat1 in maintaining the structural integrity of centromeric chromatin to facilitate faithful chromosome segregation and proper kinetochore function.