Search

Deaminase base editing has emerged as a tool to install or correct point mutations in the genomes of living cells in a wide range of organisms. However, the genome-wide off-target effects introduced by base editors in the mammalian genome have been examined in only one study. Here, we have investigated the fidelity of cytosine base editor 4 (BE4) and adenine base editors (ABE) in mouse embryos using unbiased whole-genome sequencing of a family-based trio cohort. The same sgRNA was used for BE4 and ABE. We demonstrate that BE4-edited mice carry an excess of single-nucleotide variants and deletions compared to ABE-edited mice and controls. Therefore, an optimization of cytosine base editors is required to improve its fidelity. While the remarkable fidelity of ABE has implications for a wide range of applications, the occurrence of rare aberrant C-to-T conversions at specific target sites needs to be addressed., Hye Kyung Lee, Harold E. Smith et al. examined the fidelity of cytosine base editor 4 (BE4) and adenine base editors (ABEs) in mouse embryos using whole-genome sequencing of a family-based trio cohort. They show that BE4-edited mice carry more single-nucleotide variants and deletions than ABE-edited mice.

The C. elegans dauer is an alternative third stage larva induced by dense population and adverse environmental conditions. Genes whose mutants caused dauer formation constitutive (Daf-c) and dauer formation defective (Daf-d) phenotypes were ordered via epistasis into a signaling network, with upstream DAF-7/TGF-beta and DAF-11/receptor guanylyl cyclase defining sensory branches and downstream DAF-2/Insulin receptor and DAF-12/nuclear hormone receptor executing the dauer decision. Mutations in the Scd genes were defined as incompletely penetrant suppressors of the constitutive dauer phenotype conferred by mutation of the DAF-7/TGF-beta signaling axis. SCD-2 was previously shown to be an ortholog of mammalian ALK (Anaplastic Lymphoma Kinase), a receptor tyrosine kinase. Mutations disrupting the HEN-1/Jeb ligand, SOC-1/DOS/GAB adaptor protein and SMA-5/ERK5 atypical MAP Kinase caused Scd phenotypes similar to that of mutant SCD-2. This group regulated expression from a TGF-beta-responsive GFP reporter. Here we find that a strain harboring a mutation in the uncharacterized SCD-4 is mutant for MLK-1, the C. elegans ortholog of mammalian Mixed Lineage Kinase and Drosophila slipper (slpr), a MAP3 kinase. We validated this finding by showing that a previously characterized deletion in MLK-1 caused a Scd phenotype similar to that of mutant SCD-4 and altered expression from the TGF-beta-responsive GFP reporter, suggesting that SCD-4 and MLK-1 are the same protein. Based on shared phenotypes and molecular identities, we hypothesize that MLK-1 functions as a MAP3K in the SCD-2/ALK cascade that signals through SMA-5/ERK5 MAP Kinase to modulate the output of the TGF-beta cascade controlling dauer formation in response to environmental cues.

BACKGROUND AND OBJECTIVES Next generation sequencing (NGS) has promising applications in transfusion medicine. Exome sequencing (ES) is increasingly used in the clinical setting, and blood group interpretation is an additional value that could be extracted from existing data sets. We provide the first release of an open-source software tailored for this purpose and describe its validation with three blood group systems. MATERIALS AND METHODS The DTM-Tools algorithm was designed and used to analyse 1018 ES NGS files from the ClinSeq® cohort. Predictions were correlated with serology for 5 antigens in a subset of 108 blood samples. Discrepancies were investigated with alternative phenotyping and genotyping methods, including a long-read NGS platform. RESULTS Of 116 genomic variants queried, those corresponding to 18 known KEL, FY and JK alleles were identified in this cohort. 596 additional exonic variants were identified KEL, ACKR1 and SLC14A1, including 58 predicted frameshifts. Software predictions were validated by serology in 108 participants; one case in the FY blood group and three cases in the JK blood group were discrepant. Investigation revealed that these discrepancies resulted from (1) clerical error, (2) serologic failure to detect weak antigenic expression and (3) a frameshift variant absent in blood group databases. CONCLUSION DTM-Tools can be employed for rapid Kell, Duffy and Kidd blood group antigen prediction from existing ES data sets; for discrepancies detected in the validation data set, software predictions proved accurate. DTM-Tools is open-source and in continuous development.

Ciliary microtubules are subject to post-translational modifications that act as a “Tubulin Code” to regulate motor traffic, binding proteins and stability. In humans, loss of CCP1, a cytosolic carboxypeptidase and tubulin deglutamylating enzyme, causes infantile-onset neurodegeneration. In C. elegans, mutations in ccpp-1, the homolog of CCP1, result in progressive degeneration of neuronal cilia and loss of neuronal function. To identify genes that regulate microtubule glutamylation and ciliary integrity, we performed a forward genetic screen for suppressors of ciliary degeneration in ccpp-1 mutants. We isolated the ttll-5(my38) suppressor, a mutation in the tubulin tyrosine ligase-like glutamylase gene. We show that mutation in ttll-4, ttll-5, or ttll-11 gene suppressed the hyperglutamylation-induced loss of microtubules and kinesin-2 mislocalization in ccpp-1 cilia. We also identified the nekl-4(my31) suppressor, an allele affecting the NIMA (Never in Mitosis A)-related kinase NEKL-4/NEK10. In humans, NEK10 mutation causes bronchiectasis, an airway and mucociliary transport disorder caused by defective motile cilia. C. elegans NEKL-4 does not localize to cilia yet plays a role in regulating axonemal microtubule stability. This work defines a pathway in which glutamylation, a component of the Tubulin Code, is written by TTLL-4, TTLL-5, and TTLL-11; is erased by CCPP-1; is read by ciliary kinesins; and its downstream effects are modulated by NEKL-4 activity. Identification of regulators of microtubule glutamylation in diverse cellular contexts is important to the development of effective therapies for disorders characterized by changes in microtubule glutamylation. By identifying C. elegans genes important for neuronal and ciliary stability, our work may inform research into human ciliopathies and neurodegenerative diseases.

Given their copy number differences and unique modes of inheritance, the evolved gene content and expression of sex chromosomes is unusual. In many organisms the X and Y chromosomes are inactivated in spermatocytes, possibly as a defense mechanism against insertions into unpaired chromatin. In addition to current sex chromosomes, Drosophila has a small gene-poor X-chromosome relic (4th) that re-acquired autosomal status. Here we use single cell RNA-Seq on fly larvae to demonstrate that the single X and pair of 4th chromosomes are specifically inactivated in primary spermatocytes, based on measuring all genes or a set of broadly expressed genes in testis we identified. In contrast, genes on the single Y chromosome become maximally active in primary spermatocytes. Reduced X transcript levels are due to failed activation of RNA-Polymerase-II by phosphorylation of Serine 2 and 5., Sex chromosome gene content and expression is unusual. Here the authors use single cell RNA-Seq on Drosophila larvae to demonstrate that the single X and pair of 4th chromosomes are specifically inactivated in primary spermatocytes, while genes on the single Y chromosome become maximally active in primary spermatocytes.

Oxidative damage contributes to human diseases of aging including diabetes, cancer, and cardiovascular disorders. Reactive oxygen species resulting from xenobiotic and endogenous metabolites are sensed by a poorly understood process, triggering a cascade of regulatory factors and leading to the activation of the transcription factor Nrf2 (Nuclear factor-erythroid-related factor 2, SKN-1 in Caenorhabditis elegans). Nrf2/SKN-1 activation promotes the induction of the phase II detoxification system that serves to limit oxidative stress. We have extended a previous C. elegans genetic approach to explore the mechanisms by which a phase II enzyme is induced by endogenous and exogenous oxidants. The xrep (xenobiotics response pathway) mutants were isolated as defective in their ability to properly regulate the induction of a glutathione S-transferase (GST) reporter. The xrep-1 gene was previously identified as wdr-23, which encodes a C. elegans homolog of the mammalian β-propeller repeat-containing protein WDR-23. Here, we identify and confirm the mutations in xrep-2, xrep-3, and xrep-4. The xrep-2 gene is alh-6, an ortholog of a human gene mutated in familial hyperprolinemia. The xrep-3 mutation is a gain-of-function allele of skn-1. The xrep-4 gene is F46F11.6, which encodes a F-box-containing protein. We demonstrate that xrep-4 alters the stability of WDR-23 (xrep-1), a key regulator of SKN-1 (xrep-3). Epistatic relationships among the xrep mutants and their interacting partners allow us to propose an ordered genetic pathway by which endogenous and exogenous stressors induce the phase II detoxification response.

The zinc finger protein CTCF has been invoked in establishing boundaries between genes, thereby controlling spatial and temporal enhancer activities. However, there is limited genetic evidence to support the concept that these boundaries restrict the search space of enhancers. We have addressed this question in the casein locus containing five mammary and two non-mammary genes under the control of at least seven putative enhancers. We have identified two CTCF binding sites flanking the locus and two associated with a super-enhancer. Individual deletion of these sites from the mouse genome did not alter expression of any of the genes. However, deletion of the border CTCF site separating the Csn1s1 mammary enhancer from neighboring genes resulted in the activation of Sult1d1 at a distance of more than 95 kb but not the more proximal and silent Sult1e1 gene. Loss of this CTCF site led to de novo interactions between the Sult1d1 promoter and several enhancers in the casein locus. Our study demonstrates that only one out of the four CTCF sites in the casein locus had a measurable in vivo activity. Studies on additional loci are needed to determine the biological role of CTCF sites associated with enhancers.

Deaminase base editing has emerged as a tool to install or correct point mutations in the genomes of living cells in a wide range of organisms and its ultimate success therapeutically depends on its accuracy. Here we have investigated the fidelity of cytosine base editor 4 (BE4) and adenine base editor (ABE) in mouse embryos using unbiased whole genome sequencing of a family-based trio cohort. We demonstrate that BE4-edited mice carry an excess of single-nucleotide variants and deletions compared to ABE-edited mice and controls.

A particular challenge in genome engineering has been the simultaneous introduction of mutations into linked (located on the same chromosome) loci. Although CRISPR/Cas9 has been widely used to mutate individual sites, its application in simultaneously targeting of linked loci is limited as multiple nearby double-stranded DNA breaks created by Cas9 routinely result in the deletion of sequences between the cleavage sites. Base editing is a newer form of genome editing that directly converts C∙G-to-T∙A, or A∙T-to-G∙C, base pairs without introducing double-stranded breaks, thus opening the possibility to generate linked mutations without disrupting the entire locus. Through the co-injection of two base editors and two sgRNAs into mouse zygotes, we introduced C∙G-to-T∙A transitions into two cytokine-sensing transcription factor binding sites separated by 9 kb. We determined that one enhancer activates the two flanking genes in mammary tissue during pregnancy and lactation. The ability to introduce linked mutations simultaneously in one step into the mammalian germline has implications for a wide range of applications, including the functional analysis of linked cis-elements creating disease models and correcting pathogenic mutations.

The production of viable embryos requires the coordination of many cellular processes, including protein synthesis, cytoskeletal reorganization, establishment of polarity, cell migration, cell division, and in Caenorhabditis elegans, eggshell formation. Defects in any of these processes can lead to embryonic lethality. We examined six temperature-sensitive mutants as well as one nonconditional mutant that were previously identified in genetic screens as either embryonic lethal (maternal-effect or zygotic lethal) or eggshell defective. The responsible molecular lesion for each had never been determined. After confirmation of temperature sensitivity and lethality, we performed whole-genome sequencing using a single-nucleotide polymorphism mapping strategy to pinpoint the molecular lesions. Gene candidates were confirmed by RNA interference phenocopy and/or complementation tests and one mutant was further validated by CRISPR (Clustered Regularly Interspaced Short Palidromic Repeats)/Cas9 gene editing. This approach identified new alleles of several genes that had only been previously studied by RNA interference depletion. Our identification of temperature-sensitive alleles for all of these essential genes provides an extremely useful tool for further investigation for the C. elegans community, such as the ability to address mutant phenotypes at various developmental stages and the ability to carry out suppressor/enhancer screens to identify other genes that function in a specific cellular process.

Oxidative damage contributes to human diseases of aging including diabetes, cancer, and cardiovascular disorders. Reactive oxygen species resulting from xenobiotic and endogenous metabolites are sensed by a poorly understood process, triggering a cascade of regulatory factors and leading to the activation of the transcription factor Nrf2 (Nuclear factor-erythroid-related factor 2, SKN-1 in Caenorhabditis elegans). Nrf2/SKN-1 activation promotes the induction of the phase II detoxification system that serves to limit oxidative stress. We have extended a previous C. elegans genetic approach to explore the mechanisms by which a phase II enzyme is induced by endogenous and exogenous oxidants. The xrep (xenobiotics response pathway) mutants were isolated as defective in their ability to properly regulate the induction of a glutathione S-transferase (GST) reporter. The xrep-1 gene was previously identified as wdr-23, which encodes a C. elegans homolog of the mammalian β-propeller repeat-containing protein WDR-23. Here, we identify and confirm the mutations in xrep-2, xrep-3, and xrep-4. The xrep-2 gene is alh-6, an ortholog of a human gene mutated in familial hyperprolinemia. The xrep-3 mutation is a gain-of-function allele of skn-1. The xrep-4 gene is F46F11.6, which encodes a F-box-containing protein. We demonstrate that xrep-4 alters the stability of WDR-23 (xrep-1), a key regulator of SKN-1 (xrep-3). Epistatic relationships among the xrep mutants and their interacting partners allow us to propose an ordered genetic pathway by which endogenous and exogenous stressors induce the phase II detoxification response.

Primary sex-determination "switches" evolve rapidly, but Doublesex (DSX)-related transcription factors (DMRTs) act downstream of these switches to control sexual development in most animal species. Drosophila dsx encodes female- and male-specific isoforms (DSX(F) and DSX(M)), but little is known about how dsx controls sexual development, whether DSX(F) and DSX(M) bind different targets, or how DSX proteins direct different outcomes in diverse tissues. We undertook genome-wide analyses to identify DSX targets using in vivo occupancy, binding site prediction, and evolutionary conservation. We find that DSX(F) and DSX(M) bind thousands of the same targets in multiple tissues in both sexes, yet these targets have sex- and tissue-specific functions. Interestingly, DSX targets show considerable overlap with targets identified for mouse DMRT1. DSX targets include transcription factors and signaling pathway components providing for direct and indirect regulation of sex-biased expression.

Homology-directed repair (HDR) of double-strand DNA breaks is a promising method for genome editing, but is thought to be less efficient than error-prone nonhomologous end joining in most cell types. We have investigated HDR of double-strand breaks induced by CRISPR-associated protein 9 (Cas9) in Caenorhabditis elegans. We find that HDR is very robust in the C. elegans germline. Linear repair templates with short (∼30–60 bases) homology arms support the integration of base and gene-sized edits with high efficiency, bypassing the need for selection. Based on these findings, we developed a systematic method to mutate, tag, or delete any gene in the C. elegans genome without the use of co-integrated markers or long homology arms. We generated 23 unique edits at 11 genes, including premature stops, whole-gene deletions, and protein fusions to antigenic peptides and GFP. Whole-genome sequencing of five edited strains revealed the presence of passenger variants, but no mutations at predicted off-target sites. The method is scalable for multi-gene editing projects and could be applied to other animals with an accessible germline.

Form follows function, and this maxim holds particularly true for the nematode sperm cell. Motility is essential for fertilization, and the process of spermatogenesis culminates in the production of a crawling spermatozoon with an extended pseudopod. However, the morphological similarity to amoeboid cells of other organisms is not conserved at the molecular level. Instead of utilizing the actin cytoskeleton and motor proteins, the pseudopod moves via the regulated assembly and disassembly of filaments composed of the major sperm protein (MSP). The current work reviews the structure and dynamics of MSP filament formation, the critical role of pH in MSP assembly, and the components that regulate this process. The combination of cytological, biochemical, and genetic approaches in this relatively simple system make nematode sperm an attractive model for investigating the mechanics of amoeboid cell motility.

Screening for suppressor mutations is a powerful method to isolate genes that function in a common pathway or process. Because suppressor mutations often do not have phenotypes on their own, cloning of suppressor loci can be challenging. A method combining whole-genome sequencing (WGS) and single nucleotide polymorphism (SNP) mapping (WGS/SNP mapping) was developed to identify mutations with visible phenotypes in C. elegans. We show here that WGS/SNP mapping is an efficient method to map suppressor mutations without the need for previous phenotypic characterization. Using RNA-mediated interference to test candidate loci identified by WGS/SNP mapping, we identified 10 extragenic and six intragenic suppressors of mbk-2, a DYRK family kinase required for the transition from oocyte to zygote. Remarkably, seven suppressors are mutations in cell-cycle regulators that extend the timing of the oocyte-to-zygote transition.

Background—Although in vitro studies have identified numerous possible targets, the molecules that mediate the in vivo effects of volatile anesthetics remain largely unknown. The mammalian ryanodine receptor (Ryr) is a known halothane target, and we hypothesized that it has a central role in anesthesia. Methods—Gene function of the Drosophila Ryr (dRyr) was manipulated in the whole body or in specific tissues using a collection of mutants and transgenes, and responses to halothane were measured with a reactive climbing assay. Cellular responses to halothane were studied using Ca 2+ imaging and patch clamp electrophysiology. Results—Halothane potency strongly correlates with dRyr gene copy number, and missense mutations in regions known to be functionally important in mammalian Ryrs gene cause dominant hypersensitivity. Tissue-specific manipulation of dRyr shows that expression in neurons and glia, but not muscle, mediates halothane sensitivity. In cultured cells, halothane-induced Ca 2+ efflux is strictly dRyr-dependent, suggesting a close interaction between halothane and dRyr. Ca 2+ imaging and electrophysiology of Drosophila central neurons reveal halothane-induced Ca 2+ flux that is altered in dRyr mutants and correlates with strong hyperpolarization. Conclusions—In Drosophila, neurally-expressed dRyr mediates a substantial proportion of halothane's anesthetic effects in vivo, is potently activated by halothane in vitro, and activates an inhibitory conductance. Our results provide support for Ryr as an important mediator of immobilization by volatile anesthetics.

The permeation of antibiotics through bacterial membranes to their target site is a crucial determinant of drug activity but in many cases remains poorly understood. During screening efforts to discover new broad-spectrum antibiotic compounds from marine sponge samples, we identified a new analog of the peptidyl nucleoside antibiotic blasticidin S that exhibited up to 16-fold-improved potency against a range of laboratory and clinical bacterial strains which we named P10. Whole-genome sequencing of laboratory-evolved strains of Staphylococcus aureus resistant to blasticidin S and P10, combined with genome-wide assessment of the fitness of barcoded Escherichia coli knockout strains in the presence of the antibiotics, revealed that restriction of cellular access was a key feature in the development of resistance to this class of drug. In particular, the gene encoding the well-characterized multidrug efflux pump NorA was found to be mutated in 69% of all S. aureus isolates resistant to blasticidin S or P10. Unexpectedly, resistance was associated with inactivation of norA , suggesting that the NorA transporter facilitates cellular entry of peptidyl nucleosides in addition to its known role in the efflux of diverse compounds, including fluoroquinolone antibiotics.

The development of the single cell layer skin or hypodermis of Caenorhabditis elegans is an excellent model for understanding cell fate specification and differentiation. Early in C. elegans embryogenesis, six rows of hypodermal cells adopt dorsal, lateral or ventral fates that go on to display distinct behaviors during larval life. Several transcription factors are known that function in specifying these major hypodermal cell fates, but our knowledge of the specification of these cell types is sparse, particularly in the case of the ventral hypodermal cells, which become Vulval Precursor Cells and form the vulval opening in response to extracellular signals. Previously, the gene pvl-4 was identified in a screen for mutants with defects in vulval development. We found by whole genome sequencing that pvl-4 is the Paired-box gene pax-3, which encodes the sole PAX-3 transcription factor homolog in C. elegans. pax-3 mutants show embryonic and larval lethality, and body morphology abnormalities indicative of hypodermal cell defects. We report that pax-3 is expressed in ventral P cells and their descendants during embryogenesis and early larval stages, and that in pax-3 reduction-of-function animals the ventral P cells undergo a cell fate transformation and express several markers of the lateral seam cell fate. Furthermore, forced expression of pax-3 in the lateral hypodermal cells causes them to lose expression of seam cell markers. We propose that pax-3 functions in the ventral hypodermal cells to prevent these cells from adopting the lateral seam cell fate. pax-3 represents the first gene required for specification solely of the ventral hypodermal fate in C. elegans providing insights into cell type diversification.

SummaryEvolutionary transitions between sexual modes could be potent forces in genome evolution [1, 2]. Several Caenorhabditis nematode species have evolved self-fertile hermaphrodites from the obligately outcrossing females of their ancestors [3]. We explored the relationship between sexual mode and global gene expression by comparing two selfing species, C. elegans and C. briggsae, with three phylogenetically informative outcrossing relatives, C. remanei, C. brenneri, and C. japonica. Adult transcriptome assemblies from the selfing species are consistently and strikingly smaller than those of the outcrossing species. Against this background of overall simplification, genes conserved in multiple outcrossing species with strong sex-biased expression are even more likely to be missing from the genomes of the selfing species. In addition, the sexual regulation of remaining transcripts has diverged markedly from the ancestral pattern in both selfing lineages, though in distinct ways. Thus, both the complexity and the sexual specialization of transciptomes are rapidly altered in response to the evolution of self-fertility. These changes may result from the combination of relaxed sexual selection and a recently reported genetic mechanism favoring genome shrinkage [4] in partial selfers.

Signal Transducers and Activators of Transcription (STAT) 5A/B regulate cytokine-inducible genes upon binding to GAS motifs. It is not known what percentage of genes with GAS motifs bind to and are regulated by STAT5. Moreover, it is not clear whether genome-wide STAT5 binding is modulated by its concentration. To clarify these issues we established genome-wide STAT5 binding upon growth hormone (GH) stimulation of wild-type (WT) mouse embryonic fibroblasts (MEFs) and MEFs overexpressing STAT5A more than 20-fold. Upon GH stimulation, 23 827 and 111 939 STAT5A binding sites were detected in WT and STAT5A overexpressing MEFs, respectively. 13 278 and 71 561 peaks contained at least one GAS motif. 1586 and 8613 binding sites were located within 2.5 kb of promoter sequences, respectively. Stringent filtering revealed 78 genes in which the promoter/upstream region (-10 kb to +0.5 kb) was recognized by STAT5 both in WT and STAT5 overexpressing MEFs and 347 genes that bound STAT5 only in overexpressing cells. Genome-wide expression analyses identified that the majority of STAT5-bound genes was not under GH control. Up to 40% of STAT5-bound genes were not expressed. For the first time we demonstrate the magnitude of opportunistic genomic STAT5 binding that does not translate into transcriptional activation of neighboring genes.

Bacterial production of recombinant proteins offers several advantages over alternative expression methods and remains the system of choice for many structural genomics projects. However, a large percentage of targets accumulate as insoluble inclusion bodies rather than soluble protein, creating a significant bottleneck in the protein production pipeline. Numerous strategies have been reported that can improve in vivo protein solubility, but most do not scale easily for high-throughput expression screening. To understand better the host cell response to the accumulation of insoluble protein, we determined genome-wide changes in bacterial gene expression upon induction of either soluble or insoluble target proteins. By comparing transcriptional profiles for multiple examples from the soluble or insoluble class, we identified a pattern of gene expression that correlates strongly with protein solubility. Direct targets of the sigma32 heat shock sigma factor, which includes genes involved in protein folding and degradation, were highly expressed in response to induction of insoluble protein. This same group of genes was also upregulated by insoluble protein accumulation under a different growth regime, indicating that sigma32-mediated gene expression is a general response to protein insolubility. This knowledge provides a starting point for the rational design of growth parameters and host strains with improved protein solubility characteristics. Summary Problems with protein solubility are frequently encountered when recombinant proteins are expressed in E. coli. The bacterial host responds to this problem by increasing expression of the protein folding machinery via the heat shock sigma factor sigma32. Manipulation of the sigma32 regulon might provide a general mechanism for improving recombinant protein solubility.

Next-generation sequencing provides a rapid and powerful method for mutation identification. Herein is described a workflow for sample preparation to allow the simultaneous mapping and identification of candidate mutations by whole-genome sequencing in Caenorhabditis elegans. The protocol is designed for small numbers of worms to accommodate classes of mutations, such as lethal and sterile alleles, that are difficult to identify by traditional means.

This protocols is from: Alexandre Paix, et al. (2014) Scalable and Versatile Genome Editing Using Linear DNAs with Microhomology to Cas9 Sites in Caenorhabditis elegans.Genetics198:1347-1356;doi:10.1534/genetics.114.170423 Please see the full manuscript or additional details

Bone morphogenetic proteins (BMPs) belong to the transforming growth factor β (TGFβ) superfamily of secreted molecules. BMPs play essential roles in multiple developmental and homeostatic processes in metazoans. Malfunction of the BMP pathway can cause a variety of diseases in humans, including cancer, skeletal disorders and cardiovascular diseases. Identification of factors that ensure proper spatiotemporal control of BMP signaling is critical for understanding how this pathway is regulated. We have used a unique and sensitive genetic screen to identify the plasma membrane-localized tetraspanin TSP-21 as a key new factor in the C. elegans BMP-like “Sma/Mab” signaling pathway that controls body size and postembryonic M lineage development. We showed that TSP-21 acts in the signal-receiving cells and genetically functions at the ligand-receptor level. We further showed that TSP-21 can associate with itself and with two additional tetraspanins, TSP-12 and TSP-14, which also promote Sma/Mab signaling. TSP-12 and TSP-14 can also associate with SMA-6, the type I receptor of the Sma/Mab pathway. Finally, we found that glycosphingolipids, major components of the tetraspanin-enriched microdomains, are required for Sma/Mab signaling. Our findings suggest that the tetraspanin-enriched membrane microdomains are important for proper BMP signaling. As tetraspanins have emerged as diagnostic and prognostic markers for tumor progression, and TSP-21, TSP-12 and TSP-14 are all conserved in humans, we speculate that abnormal BMP signaling due to altered expression or function of certain tetraspanins may be a contributing factor to cancer development., Author Summary The bone morphogenetic protein (BMP) signaling pathway is required for multiple developmental processes during metazoan development. Various diseases, including cancer, can result from mis-regulation of the BMP pathway. Thus, it is critical to identify factors that ensure proper regulation of BMP signaling. Using the nematode C. elegans, we have devised a highly specific and sensitive genetic screen to identify new modulators in the BMP pathway. Through this screen, we identified three conserved tetraspanin molecules as novel factors that function to promote BMP signaling in a living organism. We further showed that these three tetraspanins likely form a complex and function together with glycosphingolipids to promote BMP signaling. Recent studies have implicated several tetraspanins in cancer initiation, progression and metastasis in mammals. Our findings suggest that the involvement of tetraspanins in cancer may partially be due to their function in modulating the activity of BMP signaling.

Introduction: Accurate typing of patient and donor red blood cell (RBC) antigens is critical for safe transfusion practice. Although blood typing is traditionally accomplished by serology, genotyping methods to predict RBC antigens have proven valuable in a growing number of situations such as recently-transfused patients, scarcity of typing reagents, and indeterminate serologic results. However current RBC genotyping assays address a limited number of blood group genes and associated variants, and may not detect novel genetic changes and certain rare but clinically-significant variants. Next generation sequencing (NGS) technology provides an appealing alternative technology, allowing the user to examine a patient's entire genome or exome in a high-throughput manner. Whereas efforts are underway in multiple fields to apply exome sequencing (ES) for diagnostic, prognostic, and treatment purposes, Transfusion Medicine, with its extensive clinical genomic database, should find ready application from this approach. We describe here the creation of an algorithm to interpret NGS into a predicted extended RBC phenotype, and its application to analyze ES data from 245 participants of the ClinSeq® sequencing cohort. Methods: RyLAN (Red Cell and Lymphocyte Antigen prediction from NGS) was created as an open-source Python application that takes an NGS sorted binary alignment matrix (.bam) file and index as input. The software interacts with a non-relational database that encodes genomic blood group coordinates and phenotype interpretation rules, and yields a predicted extended RBC phenotype and quality parameters. Hard filters for mapping quality, depth, vcf QUAL, and fraction of alternate allele can be modified per individual genomic coordinate. The output is provided as a MongoDB document to facilitate advanced bulk queries and statistical analysis. We employed RyLAN to analyze 245 ES NGS files from the ClinSeq® cohort, using a database of 176 known antigenic, null, and weak blood group single nucleotide variants in 27 blood group genes as input. Results: The cohort consisted of 115 females and 130 males; 89% of participants self-described as white race, non- Hispanic ethnicity. Three percent of participants self-described as Hispanic or Latino, 4% as Asian, 2% with African ancestry, and the remaining as mixed or unknown race. From the total 176 genomic positions analyzed, 160 were not addressed by current commercially-available RBC genotyping platforms. The average read depth for the positions of interest was 78.2, and the average vcf QUAL value was 968. The highest variant nucleotide frequency was observed at the Fya/Fyb and Jka/Jkb loci (275 and 223 total haplotype variant calls, respectively). Among other phenotypes, RyLAN predicted 4 instances of heterozygosity for the KEL*02N.17 allele, 5 heterozygous individuals for the weak FY* X allele, 32 total heterozygous samples for various weak Kidd alleles, 2 homozygous individuals for weak Kidd expression, 1 heterozygosity for Lu6/Lu9, 1 SC:1,2 case, 1 Co(a-b+) predicted phenotype, and a total of 19 RHAG*01.04 and 47 KLF1*BGM12 alleles. Limited areas of the BCAM, KLF1, KEL, FUT7, ERMAP and CR1 genes failed quality filters repeatedly, and careful review indicated that these regions were not captured in the ES libraries. The ACKR1 promoter GATA-binding site variant was present in every sample and predicted all cases of self-reported African ancestry. Conclusions: We describe a new, open-source informatics tool to translate NGS data into a predicted extended RBC phenotype, and demonstrate its application through the analysis of 245 ClinSeq® ES files. Most predicted antigen frequencies were as expected for the ethnic composition of our cohort. We detected a higher frequency of the RHAG p.V270I and KLF1 p.S102P variants than expected, findings that are in agreement with the 1000 Genomes Project and warrant further study. Our analysis also corroborates the relative frequency of the JK*01W.01 allele, and the presence of the JK*01W.03 and JK*01W.04 alleles in the Caucasian population, which can lead to serologic discrepancies in other genotyping platforms. Serologic confirmation of these findings is being conducted. Further study of genomic data across multiple ethnic groups can help refine knowledge of blood group gene polymorphisms and their clinical association. Disclosures No relevant conflicts of interest to declare.

Accurate gene model annotation of reference genomes is critical for making them useful. The modENCODE project has improved the D. melanogaster genome annotation by using deep and diverse high-throughput data. Since transcriptional activity that has been evolutionarily conserved is likely to have an advantageous function, we have performed large-scale interspecific comparisons to increase confidence in predicted annotations. To support comparative genomics, we filled in divergence gaps in the Drosophila phylogeny by generating draft genomes for eight new species. For comparative transcriptome analysis, we generated mRNA expression profiles on 81 samples from multiple tissues and developmental stages of 15 Drosophila species, and we performed cap analysis of gene expression in D. melanogaster and D. pseudoobscura. We also describe conservation of four distinct core promoter structures composed of combinations of elements at three positions. Overall, each type of genomic feature shows a characteristic divergence rate relative to neutral models, highlighting the value of multispecies alignment in annotating a target genome that should prove useful in the annotation of other high priority genomes, especially human and other mammalian genomes that are rich in noncoding sequences. We report that the vast majority of elements in the annotation are evolutionarily conserved, indicating that the annotation will be an important springboard for functional genetic testing by the Drosophila community., National Institutes of Health (U.S.) (NIDDK (DK015600-18)), National Institutes of Health (U.S.) (Extramural program (1ROIGM082843)), National Institutes of Health (U.S.) (Extramural program (U01HB004271))

The production of multiple transcript isoforms from one gene is a major source of transcriptome complexity. RNA-Seq experiments, in which transcripts are converted to cDNA and sequenced, allow the resolution and quantification of alternative transcript isoforms. However, methods to analyze splicing are underdeveloped and errors resulting in incorrect splicing calls occur in every experiment. We used RNA-Seq data to develop sequencing and aligner error models. By applying these error models to known input from simulations, we found that errors result from false alignment to minor splice motifs and antisense stands, shifted junction positions, paralog joining, and repeat induced gaps. By using a series of quantitative and qualitative filters, we eliminated diagnosed errors in the simulation, and applied this to RNA-Seq data from Drosophila melanogaster heads. We used high-confidence junction detections to specifically interrogate local splicing differences between transcripts. This method out-performed commonly used RNA-seq methods to identify known alternative splicing events in the Drosophila sex determination pathway. We describe a flexible software package to perform these tasks called Splicing Analysis Kit (Spanki), available at http://www.cbcb.umd.edu/software/spanki . Splice-junction centric analysis of RNA-Seq data provides advantages in specificity for detection of alternative splicing. Our software provides tools to better understand error profiles in RNA-Seq data and improve inference from this new technology. The splice-junction centric approach that this software enables will provide more accurate estimates of differentially regulated splicing than current tools.

IME1 is required in yeast for meiosis and for expression of IME2 and other early meiotic genes. IME1 is a 360-amino acid polypeptide with central and C-terminal tyrosine-rich regions. We report here that a fusion protein composed of the lexA DNA-binding domain and IME1 activates transcription in vivo of a reporter gene containing upstream lexA binding sites. Activation by the fusion protein shares several features with natural IME1 activity: both are dependent on the RIM11 gene product; both are impaired by the same ime1 missense mutations; both are restored by intragenic suppressors. The central tyrosine-rich region is sufficient to activate transcription when fused to lexA. Deletion of this putative activation domain results in a defective IME1 derivative. Function of the deletion derivative is restored by fusion to the acidic Herpesvirus VP16 activation domain. The C-terminal tyrosine-rich region is dispensable for transcriptional activation; rather it renders activation dependent upon starvation and RIM11. Immunofluorescence studies indicate that an IME1-lacZ fusion protein is concentrated in the nucleus. These observations are consistent with a model in which IME1 normally stimulates IME2 expression by providing a transcriptional activation domain at the IME2 5' regulatory region.

Small regulatory RNAs are key regulators of gene expression. One class of small regulatory RNAs, termed the endogenous small interfering RNAs (endo siRNAs), is thought to negatively regulate cellular transcripts via an RNA interference (RNAi)-like mechanism termed endogenous RNAi (endo RNAi). A complex of proteins composed of ERI-1/3/5, RRF-3, and DICER (the ERI/DICER complex) mediates endo RNAi processes in Caenorhabditis elegans. We conducted a genetic screen to identify additional components of the endo RNAi machinery. Our screen recovered alleles of eri-9, which encodes a novel DICER-interacting protein, and a missense mutation within the helicase domain of DICER [DCR-1(G492R)]. ERI-9(−) and DCR-1(G492) animals exhibit defects in endo siRNA expression and a concomitant failure to regulate mRNAs that exhibit sequence homology to these endo siRNAs, indicating that ERI-9 and the DCR-1 helicase domain function in the C. elegans endo RNAi pathway. We define a subset of Eri mutant animals (including eri-1, rrf-3, eri-3, and dcr-1, but not eri-9 or ergo-1) that exhibit temperature-sensitive, sperm-specific sterility and defects in X chromosome segregation. Among these mutants we find multiple aberrations in sperm development beginning with cytokinesis and extending through terminal differentiation. These results identify novel components of the endo RNAi machinery, demonstrate differential requirements for the Eri factors in the sperm-producing germline, and begin to delineate the functional requirement for the ERI/DICER complex in sperm development.