A global metagenomic map of urban microbiomes and antimicrobial resistance.

We present a global atlas of 4,728 metagenomic samples from mass-transit systems in 60 cities over 3 years, representing the first systematic, worldwide catalog of the urban microbial ecosystem. This atlas provides an annotated, geospatial profile of microbial strains, functional characteristics, antimicrobial resistance (AMR) markers, and genetic elements, including 10,928 viruses, 1,302 bacteria, 2 archaea, and 838,532 CRISPR arrays not found in reference databases. We identified 4,246 known species of urban microorganisms and a consistent set of 31 species found in 97% of samples that were distinct from human commensal organisms. Profiles of AMR genes varied widely in type and density across cities. Cities showed distinct microbial taxonomic signatures that were driven by climate and geographic differences. These results constitute a high-resolution global metagenomic atlas that enables discovery of organisms and genes, highlights potential public health and forensic applications, and provides a culture-independent view of AMR burden in cities.

Disruptions of topological chromatin domains cause pathogenic rewiring of gene-enhancer interactions.

Mammalian genomes are organized into megabase-scale topologically associated domains (TADs). We demonstrate that disruption of TADs can rewire long-range regulatory architecture and result in pathogenic phenotypes. We show that distinct human limb malformations are caused by deletions, inversions, or duplications altering the structure of the TAD-spanning WNT6/IHH/EPHA4/PAX3 locus. Using CRISPR/Cas genome editing, we generated mice with corresponding rearrangements. Both in mouse limb tissue and patient-derived fibroblasts, disease-relevant structural changes cause ectopic interactions between promoters and non-coding DNA, and a cluster of limb enhancers normally associated with Epha4 is misplaced relative to TAD boundaries and drives ectopic limb expression of another gene in the locus. This rewiring occurred only if the variant disrupted a CTCF-associated boundary domain. Our results demonstrate the functional importance of TADs for orchestrating gene expression via genome architecture and indicate criteria for predicting the pathogenicity of human structural variants, particularly in non-coding regions of the human genome.

In vitro and in vivo evidence that the switch from calcineurin to mTOR inhibitors may be a strategy for immunosuppression in Epstein-Barr virus-associated post-transplant lymphoproliferative disorder.

Post-transplant lymphoproliferative disorder is a life-threatening complication of immunosuppression following transplantation mediated by failure of T cells to control Epstein-Barr virus (EBV)-infected and transformed B cells. Typically, a modification or reduction of immunosuppression is recommended, but insufficiently defined thus far. In order to help delineate this, we characterized EBV-antigen-specific T cells and lymphoblastoid cell lines from healthy donors and in patients with a kidney transplant in the absence or presence of the standard immunosuppressants tacrolimus, cyclosporin A, prednisolone, rapamycin, and mycophenolic acid. Phenotypes of lymphoblastoid cell-lines and T cells, T cell-receptor-repertoire diversity, and T-cell reactivity upon co-culture with autologous lymphoblastoid cell lines were analyzed. Rapamycin and mycophenolic acid inhibited lymphoblastoid cell-line proliferation. T cells treated with prednisolone and rapamycin showed nearly normal cytokine production. Proliferation and the viability of T cells were decreased by mycophenolic acid, while tacrolimus and cyclosporin A were strong suppressors of T-cell function including their killing activity. Overall, our study provides a basis for the clinical decision for the modification and reduction of immunosuppression and adds information to the complex balance of maintaining anti-viral immunity while preventing acute rejection. Thus, an immunosuppressive regime based on mTOR inhibition and reduced or withdrawn calcineurin inhibitors could be a promising strategy for patients with increased risk of or manifested EBV-associated post-transplant lymphoproliferative disorder.

Alternative splicing of BUD13 determines the severity of a developmental disorder with lipodystrophy and progeroid features.

PURPOSE: In this study we aimed to identify the molecular genetic cause of a progressive multisystem disease with prominent lipodystrophy. METHODS: In total, 5 affected individuals were investigated using exome sequencing. Dermal fibroblasts were characterized using RNA sequencing, proteomics, immunoblotting, immunostaining, and electron microscopy. Subcellular localization and rescue studies were performed. RESULTS: We identified a lipodystrophy phenotype with a typical facial appearance, corneal clouding, achalasia, progressive hearing loss, and variable severity. Although 3 individuals showed stunted growth, intellectual disability, and died within the first decade of life (A1, A2, and A3), 2 are adults with normal intellectual development (A4 and A5). All individuals harbored an identical homozygous nonsense variant affecting the retention and splicing complex component BUD13. The nucleotide substitution caused alternative splicing of BUD13 leading to a stable truncated protein whose expression positively correlated with disease expression and life expectancy. In dermal fibroblasts, we found elevated intron retention, a global reduction of spliceosomal proteins, and nuclei with multiple invaginations, which were more pronounced in A1, A2, and A3. Overexpression of both BUD13 isoforms normalized the nuclear morphology. CONCLUSION: Our results define a hitherto unknown syndrome and show that the alternative splice product converts a loss-of-function into a hypomorphic allele, thereby probably determining the severity of the disease and the survival of affected individuals.

SARS-CoV-2 Seroprevalence and Antibody Kinetics Among Health Care Workers in a Spanish Hospital After 3 Months of Follow-up.

BACKGROUND: At the COVID-19 spring 2020 pandemic peak in Spain, prevalence of SARS-CoV-2 infection in a cohort of 578 randomly selected health care workers (HCWs) from Hospital Clínic de Barcelona was 11.2%. METHODS: A follow-up survey 1 month later (April-May 2020) measured infection by rRT-PCR and IgM, IgA, and IgG to the receptor-binding domain of the spike protein by Luminex. Antibody kinetics, including IgG subclasses, was assessed until month 3. RESULTS: At month 1, the prevalence of infection measured by rRT-PCR and serology was 14.9% (84/565) and seroprevalence 14.5% (82/565). We found 25 (5%) new infections in 501 participants without previous evidence of infection. IgM, IgG, and IgA levels declined in 3 months (antibody decay rates 0.15 [95% CI, .11-.19], 0.66 [95% CI, .54-.82], and 0.12 [95% CI, .09-.16], respectively), and 68.33% of HCWs had seroreverted for IgM, 3.08% for IgG, and 24.29% for IgA. The most frequent subclass responses were IgG1 (highest levels) and IgG2, followed by IgG3, and only IgA1 but no IgA2 was detected. CONCLUSIONS: Continuous and improved surveillance of SARS-CoV-2 infections in HCWs remains critical, particularly in high-risk groups. The observed fast decay of IgA and IgM levels has implications for seroprevalence studies using these isotypes.

Genome-wide strategies identify downstream target genes of chick connective tissue-associated transcription factors.

Connective tissues support organs and play crucial roles in development, homeostasis and fibrosis, yet our understanding of their formation is still limited. To gain insight into the molecular mechanisms of connective tissue specification, we selected five zinc-finger transcription factors - OSR1, OSR2, EGR1, KLF2 and KLF4 - based on their expression patterns and/or known involvement in connective tissue subtype differentiation. RNA-seq and ChIP-seq profiling of chick limb micromass cultures revealed a set of common genes regulated by all five transcription factors, which we describe as a connective tissue core expression set. This common core was enriched with genes associated with axon guidance and myofibroblast signature, including fibrosis-related genes. In addition, each transcription factor regulated a specific set of signalling molecules and extracellular matrix components. This suggests a concept whereby local molecular niches can be created by the expression of specific transcription factors impinging on the specification of local microenvironments. The regulatory network established here identifies common and distinct molecular signatures of limb connective tissue subtypes, provides novel insight into the signalling pathways governing connective tissue specification, and serves as a resource for connective tissue development.

De Novo Mutations in SLC25A24 Cause a Craniosynostosis Syndrome with Hypertrichosis, Progeroid Appearance, and Mitochondrial Dysfunction.

Gorlin-Chaudhry-Moss syndrome (GCMS) is a dysmorphic syndrome characterized by coronal craniosynostosis and severe midface hypoplasia, body and facial hypertrichosis, microphthalmia, short stature, and short distal phalanges. Variable lipoatrophy and cutis laxa are the basis for a progeroid appearance. Using exome and genome sequencing, we identified the recurrent de novo mutations c.650G>A (p.Arg217His) and c.649C>T (p.Arg217Cys) in SLC25A24 in five unrelated girls diagnosed with GCMS. Two of the girls had pronounced neonatal progeroid features and were initially diagnosed with Wiedemann-Rautenstrauch syndrome. SLC25A24 encodes a mitochondrial inner membrane ATP-Mg/Pi carrier. In fibroblasts from affected individuals, the mutated SLC25A24 showed normal stability. In contrast to control cells, the probands' cells showed mitochondrial swelling, which was exacerbated upon treatment with hydrogen peroxide (H2O2). The same effect was observed after overexpression of the mutant cDNA. Under normal culture conditions, the mitochondrial membrane potential of the probands' fibroblasts was intact, whereas ATP content in the mitochondrial matrix was lower than that in control cells. However, upon H2O2 exposure, the membrane potential was significantly elevated in cells harboring the mutated SLC25A24. No reduction of mitochondrial DNA copy number was observed. These findings demonstrate that mitochondrial dysfunction with increased sensitivity to oxidative stress is due to the SLC25A24 mutations. Our results suggest that the SLC25A24 mutations induce a gain of pathological function and link mitochondrial ATP-Mg/Pi transport to the development of skeletal and connective tissue.

Genome-Wide Binding of Posterior HOXA/D Transcription Factors Reveals Subgrouping and Association with CTCF.

Homeotic genes code for key transcription factors (HOX-TFs) that pattern the animal body plan. During embryonic development, Hox genes are expressed in overlapping patterns and function in a partially redundant manner. In vitro biochemical screens probing the HOX-TF sequence specificity revealed largely overlapping sequence preferences, indicating that co-factors might modulate the biological function of HOX-TFs. However, due to their overlapping expression pattern, high protein homology, and insufficiently specific antibodies, little is known about their genome-wide binding preferences. In order to overcome this problem, we virally expressed tagged versions of limb-expressed posterior HOX genes (HOXA9-13, and HOXD9-13) in primary chicken mesenchymal limb progenitor cells (micromass). We determined the effect of each HOX-TF on cellular differentiation (chondrogenesis) and gene expression and found that groups of HOX-TFs induce distinct regulatory programs. We used ChIP-seq to determine their individual genome-wide binding profiles and identified between 12,721 and 28,572 binding sites for each of the nine HOX-TFs. Principal Component Analysis (PCA) of binding profiles revealed that the HOX-TFs are clustered in two subgroups (Group 1: HOXA/D9, HOXA/D10, HOXD12, and HOXA13 and Group 2: HOXA/D11 and HOXD13), which are characterized by differences in their sequence specificity and by the presence of cofactor motifs. Specifically, we identified CTCF binding sites in Group 1, indicating that this subgroup of HOX-proteins cooperates with CTCF. We confirmed this interaction by an independent biological assay (Proximity Ligation Assay) and demonstrated that CTCF is a novel HOX cofactor that specifically associates with Group 1 HOX-TFs, pointing towards a possible interplay between HOX-TFs and chromatin architecture.

GOPHER: Generator Of Probes for capture Hi-C Experiments at high Resolution.

BACKGROUND: Target enrichment combined with chromosome conformation capturing methodologies such as capture Hi-C (CHC) can be used to investigate spatial layouts of genomic regions with high resolution and at scalable costs. A common application of CHC is the investigation of regulatory elements that are in contact with promoters, but CHC can be used for a range of other applications. Therefore, probe design for CHC needs to be adapted to experimental needs, but no flexible tool is currently available for this purpose. RESULTS: We present a Java desktop application called GOPHER (Generator Of Probes for capture Hi-C Experiments at high Resolution) that implements three strategies for CHC probe design. GOPHER's simple approach is similar to the probe design of previous approaches that employ CHC to investigate all promoters, with one probe being placed at each margin of a single digest that overlaps the transcription start site (TSS) of each promoter. GOPHER's simple-patched approach extends this methodology with a heuristic that improves coverage of viewpoints in which the TSS is located near to one of the boundaries of the digest. GOPHER's extended approach is intended mainly for focused investigations of smaller gene sets. GOPHER can also be used to design probes for regions other than TSS such as GWAS hits or large blocks of genomic sequence. GOPHER additionally provides a number of features that allow users to visualize and edit viewpoints, and outputs a range of files useful for documentation, ordering probes, and downstream analysis. CONCLUSION: GOPHER is an easy-to-use and robust desktop application for CHC probe design. Source code and a precompiled executable can be downloaded from the GOPHER GitHub page at https://github.com/TheJacksonLaboratory/Gopher .

Saturation analysis of ChIP-seq data for reproducible identification of binding peaks.

Chromatin immunoprecipitation coupled with next-generation sequencing (ChIP-seq) is a powerful technology to identify the genome-wide locations of transcription factors and other DNA binding proteins. Computational ChIP-seq peak calling infers the location of protein-DNA interactions based on various measures of enrichment of sequence reads. In this work, we introduce an algorithm, Q, that uses an assessment of the quadratic enrichment of reads to center candidate peaks followed by statistical analysis of saturation of candidate peaks by 5' ends of reads. We show that our method not only is substantially faster than several competing methods but also demonstrates statistically significant advantages with respect to reproducibility of results and in its ability to identify peaks with reproducible binding site motifs. We show that Q has superior performance in the delineation of double RNAPII and H3K4me3 peaks surrounding transcription start sites related to a better ability to resolve individual peaks. The method is implemented in C++ and is freely available under an open source license.

Molecular mechanism of CHRDL1-mediated X-linked megalocornea in humans and in Xenopus model.

Chordin-Like 1 (CHRDL1) mutations cause non-syndromic X-linked megalocornea (XMC) characterized by enlarged anterior eye segments. Mosaic corneal degeneration, presenile cataract and secondary glaucoma are associated with XMC. Beside that CHRDL1 encodes Ventroptin, a secreted bone morphogenetic protein (BMP) antagonist, the molecular mechanism of XMC is not well understood yet. In a family with broad phenotypic variability of XMC, we identified the novel CHRDL1 frameshift mutation c.807_808delTC [p.H270Wfs*22] presumably causing CHRDL1 loss of function. Using Xenopus laevis as model organism, we demonstrate that chrdl1 is specifically expressed in the ocular tissue at late developmental stages. The chrdl1 knockdown directly resembles the human XMC phenotype and confirms CHRDL1 deficiency to cause XMC. Interestingly, secondary to this bmp4 is down-regulated in the Xenopus eyes. Moreover, phospho-SMAD1/5 is altered and BMP receptor 1A is reduced in a XMC patient. Together, we classify these observations as negative-feedback regulation due to the deficient BMP antagonism in XMC. As CHRDL1 is preferentially expressed in the limbal stem cell niche of adult human cornea, we assume that CHRDL1 plays a key role in cornea homeostasis. In conclusion, we provide novel insights into the molecular mechanism of XMC as well as into the specific role of CHRDL1 during cornea organogenesis, among others by the establishment of the first XMC in vivo model. We show that unravelling monogenic cornea disorders like XMC-with presumably disturbed cornea growth and differentiation-contribute to the identification of potential limbal stem cell niche factors that are promising targets for regenerative therapies of corneal injuries.

Mutations in PGAP3 impair GPI-anchor maturation, causing a subtype of hyperphosphatasia with mental retardation.

Glycosylphophatidylinositol (GPI)-anchored proteins play important roles in many biological processes, and mutations affecting proteins involved in the synthesis of the GPI anchor are reported to cause a wide spectrum of intellectual disabilities (IDs) with characteristic additional phenotypic features. Here, we describe a total of five individuals (from three unrelated families) in whom we identified mutations in PGAP3, encoding a protein that is involved in GPI-anchor maturation. Three siblings in a consanguineous Pakistani family presented with profound developmental delay, severe ID, no speech, psychomotor delay, and postnatal microcephaly. A combination of autozygosity mapping and exome sequencing identified a 13.8 Mb region harboring a homozygous c.275G>A (p.Gly92Asp) variant in PGAP3 region 17q11.2-q21.32. Subsequent testing showed elevated serum alkaline phosphatase (ALP), a GPI-anchored enzyme, in all three affected children. In two unrelated individuals in a cohort with developmental delay, ID, and elevated ALP, we identified compound-heterozygous variants c.439dupC (p.Leu147Profs(∗)16) and c.914A>G (p.Asp305Gly) and homozygous variant c.314C>G (p.Pro105Arg). The 1 bp duplication causes a frameshift and nonsense-mediated decay. Further evidence supporting pathogenicity of the missense mutations c.275G>A, c.314C>G, and c.914A>G was provided by the absence of the variants from ethnically matched controls, phylogenetic conservation, and functional studies on Chinese hamster ovary cell lines. Taken together with recent data on PGAP2, these results confirm the importance of the later GPI-anchor remodelling steps for normal neuronal development. Impairment of PGAP3 causes a subtype of hyperphosphatasia with ID, a congenital disorder of glycosylation that is also referred to as Mabry syndrome.

Homozygous and compound-heterozygous mutations in TGDS cause Catel-Manzke syndrome.

Catel-Manzke syndrome is characterized by Pierre Robin sequence and a unique form of bilateral hyperphalangy causing a clinodactyly of the index finger. We describe the identification of homozygous and compound heterozygous mutations in TGDS in seven unrelated individuals with typical Catel-Manzke syndrome by exome sequencing. Six different TGDS mutations were detected: c.892A>G (p.Asn298Asp), c.270_271del (p.Lys91Asnfs(∗)22), c.298G>T (p.Ala100Ser), c.294T>G (p.Phe98Leu), c.269A>G (p.Glu90Gly), and c.700T>C (p.Tyr234His), all predicted to be disease causing. By using haplotype reconstruction we showed that the mutation c.298G>T is probably a founder mutation. Due to the spectrum of the amino acid changes, we suggest that loss of function in TGDS is the underlying mechanism of Catel-Manzke syndrome. TGDS (dTDP-D-glucose 4,6-dehydrogenase) is a conserved protein belonging to the SDR family and probably plays a role in nucleotide sugar metabolism.

The germline variants in DNA repair genes in pediatric medulloblastoma: a challenge for current therapeutic strategies.

BACKGROUND: The defects in DNA repair genes are potentially linked to development and response to therapy in medulloblastoma. Therefore the purpose of this study was to establish the spectrum and frequency of germline variants in selected DNA repair genes and their impact on response to chemotherapy in medulloblastoma patients. METHODS: The following genes were investigated in 102 paediatric patients: MSH2 and RAD50 using targeted gene panel sequencing and NBN variants (p.I171V and p.K219fs*19) by Sanger sequencing. In three patients with presence of rare life-threatening adverse events (AE) and no detected variants in the analyzed genes, whole exome sequencing was performed. Based on combination of molecular and immunohistochemical evaluations tumors were divided into molecular subgroups. Presence of variants was tested for potential association with the occurrence of rare life-threatening AE and other clinical features. RESULTS: We have identified altogether six new potentially pathogenic variants in MSH2 (p.A733T and p.V606I), RAD50 (p.R1093*), FANCM (p.L694*), ERCC2 (p.R695C) and EXO1 (p.V738L), in addition to two known NBN variants. Five out of twelve patients with defects in either of MSH2, RAD50 and NBN genes suffered from rare life-threatening AE, more frequently than in control group (p = 0.0005). When all detected variants were taken into account, the majority of patients (8 out of 15) suffered from life-threatening toxicity during chemotherapy. CONCLUSION: Our results, based on the largest systematic study performed in a clinical setting, provide preliminary evidence for a link between defects in DNA repair genes and treatment related toxicity in children with medulloblastoma. The data suggest that patients with DNA repair gene variants could need special vigilance during and after courses of chemotherapy.

Q-nexus: a comprehensive and efficient analysis pipeline designed for ChIP-nexus.

BACKGROUND: ChIP-nexus, an extension of the ChIP-exo protocol, can be used to map the borders of protein-bound DNA sequences at nucleotide resolution, requires less input DNA and enables selective PCR duplicate removal using random barcodes. However, the use of random barcodes requires additional preprocessing of the mapping data, which complicates the computational analysis. To date, only a very limited number of software packages are available for the analysis of ChIP-exo data, which have not yet been systematically tested and compared on ChIP-nexus data. RESULTS: Here, we present a comprehensive software package for ChIP-nexus data that exploits the random barcodes for selective removal of PCR duplicates and for quality control. Furthermore, we developed bespoke methods to estimate the width of the protected region resulting from protein-DNA binding and to infer binding positions from ChIP-nexus data. Finally, we applied our peak calling method as well as the two other methods MACE and MACS2 to the available ChIP-nexus data. CONCLUSIONS: The Q-nexus software is efficient and easy to use. Novel statistics about duplication rates in consideration of random barcodes are calculated. Our method for the estimation of the width of the protected region yields unbiased signatures that are highly reproducible for biological replicates and at the same time very specific for the respective factors analyzed. As judged by the irreproducible discovery rate (IDR), our peak calling algorithm shows a substantially better reproducibility. An implementation of Q-nexus is available at http://charite.github.io/Q/ .

PGAP2 mutations, affecting the GPI-anchor-synthesis pathway, cause hyperphosphatasia with mental retardation syndrome.

Recently, mutations in genes involved in the biosynthesis of the glycosylphosphatidylinositol (GPI) anchor have been identified in a new subclass of congenital disorders of glycosylation (CDGs) with a distinct spectrum of clinical features. To date, mutations have been identified in six genes (PIGA, PIGL, PIGM, PIGN, PIGO, and PIGV) encoding proteins in the GPI-anchor-synthesis pathway in individuals with severe neurological features, including seizures, muscular hypotonia, and intellectual disability. We developed a diagnostic gene panel for targeting all known genes encoding proteins in the GPI-anchor-synthesis pathway to screen individuals matching these features, and we detected three missense mutations in PGAP2, c.46C>T, c.380T>C, and c.479C>T, in two unrelated individuals with hyperphosphatasia with mental retardation syndrome (HPMRS). The mutations cosegregated in the investigated families. PGAP2 is involved in fatty-acid GPI-anchor remodeling, which occurs in the Golgi apparatus and is required for stable association between GPI-anchored proteins and the cell-surface membrane rafts. Transfection of the altered protein constructs, p.Arg16Trp (NP_001243169.1), p.Leu127Ser, and p.Thr160Ile, into PGAP2-null cells showed only partial restoration of GPI-anchored marker proteins, CD55 and CD59, on the cell surface. In this work, we show that an impairment of GPI-anchor remodeling also causes HPMRS and conclude that targeted sequencing of the genes encoding proteins in the GPI-anchor-synthesis pathway is an effective diagnostic approach for this subclass of CDGs.

Mutations in WNT1 cause different forms of bone fragility.

We report that hypofunctional alleles of WNT1 cause autosomal-recessive osteogenesis imperfecta, a congenital disorder characterized by reduced bone mass and recurrent fractures. In consanguineous families, we identified five homozygous mutations in WNT1: one frameshift mutation, two missense mutations, one splice-site mutation, and one nonsense mutation. In addition, in a family affected by dominantly inherited early-onset osteoporosis, a heterozygous WNT1 missense mutation was identified in affected individuals. Initial functional analysis revealed that altered WNT1 proteins fail to activate canonical LRP5-mediated WNT-regulated ß-catenin signaling. Furthermore, osteoblasts cultured in vitro showed enhanced Wnt1 expression with advancing differentiation, indicating a role of WNT1 in osteoblast function and bone development. Our finding that homozygous and heterozygous variants in WNT1 predispose to low-bone-mass phenotypes might advance the development of more effective therapeutic strategies for congenital forms of bone fragility, as well as for common forms of age-related osteoporosis.

Distinct global shifts in genomic binding profiles of limb malformation-associated HOXD13 mutations.

Gene regulation by transcription factors (TFs) determines developmental programs and cell identity. Consequently, mutations in TFs can lead to dramatic phenotypes in humans by disrupting gene regulation. To date, the molecular mechanisms that actually cause these phenotypes have been difficult to address experimentally. ChIP-seq, which couples chromatin immunoprecipitation with high-throughput sequencing, allows TF function to be investigated on a genome-wide scale, enabling new approaches for the investigation of gene regulation. Here, we present the application of ChIP-seq to explore the effect of missense mutations in TFs on their genome-wide binding profile. Using a retroviral expression system in chicken mesenchymal stem cells, we elucidated the mechanism underlying a novel missense mutation in HOXD13 (Q317K) associated with a complex hand and foot malformation phenotype. The mutated glutamine (Q) is conserved in most homeodomains, a notable exception being bicoid-type homeodomains that have lysine (K) at this position. Our results show that the mutation results in a shift in the binding profile of the mutant toward a bicoid/PITX1 motif. Gene expression analysis and functional assays using in vivo overexpression studies confirm that the mutation results in a partial conversion of HOXD13 into a TF with bicoid/PITX1 properties. A similar shift was not observed with another mutation, Q317R, which is associated with brachysyndactyly, suggesting that the bicoid/PITX1-shift observed for Q317K might be related to the severe clinical phenotype. The methodology described can be used to investigate a wide spectrum of TFs and mutations that have not previously been amenable to ChIP-seq experiments.

Efficient CRISPR/Cas9 genome editing with low off-target effects in zebrafish.

Gene modifications in animal models have been greatly facilitated through the application of targeted genome editing tools. The prokaryotic CRISPR/Cas9 type II genome editing system has recently been applied in cell lines and vertebrates. However, we still have very limited information about the efficiency of mutagenesis, germline transmission rates and off-target effects in genomes of model organisms. We now demonstrate that CRISPR/Cas9 mutagenesis in zebrafish is highly efficient, reaching up to 86.0%, and is heritable. The efficiency of the CRISPR/Cas9 system further facilitated the targeted knock-in of a protein tag provided by a donor oligonucleotide with knock-in efficiencies of 3.5-15.6%. Mutation rates at potential off-target sites are only 1.1-2.5%, demonstrating the specificity of the CRISPR/Cas9 system. The ease and efficiency of the CRISPR/Cas9 system with limited off-target effects make it a powerful genome engineering tool for in vivo studies.

Strategies to improve the performance of rare variant association studies by optimizing the selection of controls.

MOTIVATION: When analyzing a case group of patients with ultra-rare disorders the ethnicities are often diverse and the data quality might vary. The population substructure in the case group as well as the heterogeneous data quality can cause substantial inflation of test statistics and result in spurious associations in case-control studies if not properly adjusted for. Existing techniques to correct for confounding effects were especially developed for common variants and are not applicable to rare variants. RESULTS: We analyzed strategies to select suitable controls for cases that are based on similarity metrics that vary in their weighting schemes. We simulated different disease entities on real exome data and show that a similarity-based selection scheme can help to reduce false positive associations and to optimize the performance of the statistical tests. Especially when data quality as well as ethnicities vary a lot in the case group, a matching approach that puts more weight on rare variants shows the best performance. We reanalyzed collections of unrelated patients with Kabuki make-up syndrome, Hyperphosphatasia with Mental Retardation syndrome and Catel-Manzke syndrome for which the disease genes were recently described. We show that rare variant association tests are more sensitive and specific in identifying the disease gene than intersection filters and should thus be considered as a favorable approach in analyzing even small patient cohorts. AVAILABILITY AND IMPLEMENTATION: Datasets used in our analysis are available at ftp://ftp.1000genomes.ebi.ac.uk./vol1/ftp/ CONTACT: : peter.krawitz@charite.de SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.