CRISPR Activation Reverses Haploinsufficiency and Functional Deficits Caused by TTN Truncation Variants

BACKGROUND

TTN truncation variants (TTNtvs) are the most common genetic lesion identified in individuals with dilated cardiomyopathy, a disease with high morbidity and mortality rates. TTNtvs reduce normal TTN (titin) protein levels, produce truncated proteins, and impair sarcomere content and function. Therapeutics targeting TTNtvs have been elusive because of the immense size of TTN, the rarity of specific TTNtvs, and incomplete knowledge of TTNtv pathogenicity.

METHODS

We adapted CRISPR activation using dCas9-VPR to functionally interrogate TTNtv pathogenicity and develop a therapeutic in human cardiomyocytes and 3-dimensional cardiac microtissues engineered from induced pluripotent stem cell models harboring a dilated cardiomyopathy-associated TTNtv. We performed guide RNA screening with custom TTN reporter assays, agarose gel electrophoresis to quantify TTN protein levels and isoforms, and RNA sequencing to identify molecular consequences of TTN activation. Cardiomyocyte epigenetic assays were also used to nominate DNA regulatory elements to enable cardiomyocyte-specific TTN activation.

RESULTS

CRISPR activation of TTN using single guide RNAs targeting either the TTN promoter or regulatory elements in spatial proximity to the TTN promoter through 3-dimensional chromatin interactions rescued TTN protein deficits disturbed by TTNtvs. Increasing TTN protein levels normalized sarcomere content and contractile function despite increasing truncated TTN protein. In addition to TTN transcripts, CRISPR activation also increased levels of myofibril assembly-related and sarcomere-related transcripts.

CONCLUSIONS

TTN CRISPR activation rescued TTNtv-related functional deficits despite increasing truncated TTN levels, which provides evidence to support haploinsufficiency as a relevant genetic mechanism underlying heterozygous TTNtvs. CRISPR activation could be developed as a therapeutic to treat a large proportion of TTNtvs.

Multifaceted roles of cohesin in regulating transcriptional loops

Cohesin is required for chromatin loop formation. However, its precise role in regulating gene transcription remains largely unknown. We investigated the relationship between cohesin and RNA Polymerase II (RNAPII) using single-molecule mapping and live-cell imaging methods in human cells. Cohesin-mediated transcriptional loops were highly correlated with those of RNAPII and followed the direction of gene transcription. Depleting RAD21, a subunit of cohesin, resulted in the loss of long-range (>100 kb) loops between distal (super-)enhancers and promoters of cell-type-specific genes. By contrast, the short-range (<50 kb) loops were insensitive to RAD21 depletion and connected genes that are mostly housekeeping. This result explains why only a small fraction of genes are affected by the loss of long-range chromatin interactions due to cohesin depletion. Remarkably, RAD21 depletion appeared to up-regulate genes located in early initiation zones (EIZ) of DNA replication, and the EIZ signals were amplified drastically without RAD21. Our results revealed new mechanistic insights of cohesin's multifaceted roles in establishing transcriptional loops, preserving long-range chromatin interactions for cell-specific genes, and maintaining timely order of DNA replication.

ChIATAC is an efficient strategy for multi-omics mapping of 3D epigenomes from low-cell inputs

Connecting genes to their cis-regulatory elements has been enabled by genome-wide mapping of chromatin interactions using proximity ligation in ChIA-PET, Hi-C, and their derivatives. However, these methods require millions of input cells for high-quality data and thus are unsuitable for many studies when only limited cells are available. Conversely, epigenomic profiling via transposase digestion in ATAC-seq requires only hundreds to thousands of cells to robustly map open chromatin associated with transcription activity, but it cannot directly connect active genes to their distal enhancers. Here, we combine proximity ligation in ChIA-PET and transposase accessibility in ATAC-seq into ChIATAC to efficiently map interactions between open chromatin loci in low numbers of input cells. We validate ChIATAC in Drosophila cells and optimize it for mapping 3D epigenomes in human cells robustly. Applying ChIATAC to primary human T cells, we reveal mechanisms that topologically regulate transcriptional programs during T cell activation.

Epigenetic activation of the FLT3 gene by ZNF384 fusion confers a therapeutic susceptibility in acute lymphoblastic leukemia

FLT3 is an attractive therapeutic target in acute lymphoblastic leukemia (ALL) but the mechanism for its activation in this cancer is incompletely understood. Profiling global gene expression in large ALL cohorts, we identify over-expression of FLT3 in ZNF384-rearranged ALL, consistently across cases harboring different fusion partners with ZNF384. Mechanistically, we discover an intergenic enhancer element at the FLT3 locus that is exclusively activated in ZNF384-rearranged ALL, with the enhancer-promoter looping directly mediated by the fusion protein. There is also a global enrichment of active enhancers within ZNF384 binding sites across the genome in ZNF384-rearranged ALL cells. Downregulation of ZNF384 blunts FLT3 activation and decreases ALL cell sensitivity to FLT3 inhibitor gilteritinib in vitro. In patient-derived xenograft models of ZNF384-rearranged ALL, gilteritinib exhibits significant anti-leukemia efficacy as a monotherapy in vivo. Collectively, our results provide insights into FLT3 regulation in ALL and point to potential genomics-guided targeted therapy for this patient population.

Guilt by association: EcDNA as a mobile transactivator in cancer

Extrachromosomal DNA (ecDNA), first described in the 1960s, is emerging as a prevalent but poorly characterized oncogenic alteration in cancer. ecDNA is a reservoir for oncogene amplification and is associated with an aggressive tumor phenotype and poor patient outcome. Despite the long-held knowledge of its existence, little is known about how ecDNA affects tumor cell behavior. Recent data reveal that ecDNA hubs are mobile transcriptional enhancers which can transactivate gene expression through chromatin interactions. Given its prevalence, structural complexity, and unequal segregation into daughter cells, ecDNA can offer selective growth advantages, contribute to intratumor heterogeneity (ITH), and accelerate tumor evolution. Future technology development is expected to transform the current paradigm for studying ecDNA and lead to therapeutic strategies targeting ecDNA vulnerabilities.

Discovery of photosynthesis genes through whole-genome sequencing of acetate-requiring mutants of Chlamydomonas reinhardtii

Large-scale mutant libraries have been indispensable for genetic studies, and the development of next-generation genome sequencing technologies has greatly advanced efforts to analyze mutants. In this work, we sequenced the genomes of 660 Chlamydomonas reinhardtii acetate-requiring mutants, part of a larger photosynthesis mutant collection previously generated by insertional mutagenesis with a linearized plasmid. We identified 554 insertion events from 509 mutants by mapping the plasmid insertion sites through paired-end sequences, in which one end aligned to the plasmid and the other to a chromosomal location. Nearly all (96%) of the events were associated with deletions, duplications, or more complex rearrangements of genomic DNA at the sites of plasmid insertion, and together with deletions that were unassociated with a plasmid insertion, 1470 genes were identified to be affected. Functional annotations of these genes were enriched in those related to photosynthesis, signaling, and tetrapyrrole synthesis as would be expected from a library enriched for photosynthesis mutants. Systematic manual analysis of the disrupted genes for each mutant generated a list of 253 higher-confidence candidate photosynthesis genes, and we experimentally validated two genes that are essential for photoautotrophic growth, CrLPA3 and CrPSBP4. The inventory of candidate genes includes 53 genes from a phylogenomically defined set of conserved genes in green algae and plants. Altogether, 70 candidate genes encode proteins with previously characterized functions in photosynthesis in Chlamydomonas, land plants, and/or cyanobacteria; 14 genes encode proteins previously shown to have functions unrelated to photosynthesis. Among the remaining 169 uncharacterized genes, 38 genes encode proteins without any functional annotation, signifying that our results connect a function related to photosynthesis to these previously unknown proteins. This mutant library, with genome sequences that reveal the molecular extent of the chromosomal lesions and resulting higher-confidence candidate genes, will aid in advancing gene discovery and protein functional analysis in photosynthesis.

Sox2 controls neural stem cell self-renewal through a Fos-centered gene regulatory network

The Sox2 transcription factor is necessary for the long-term self-renewal of neural stem cells (NSCs). Its mechanism of action is still poorly defined. To identify molecules regulated by Sox2, and acting in mouse NSC maintenance, we transduced, into Sox2-deleted NSC, genes whose expression is strongly downregulated following Sox2 loss (Fos, Jun, Egr2), individually or in combination. Fos alone rescued long-term proliferation, as shown by in vitro cell growth and clonal analysis. Furthermore, pharmacological inhibition by T-5224 of FOS/JUN AP1 complex binding to its targets decreased cell proliferation and expression of the putative target Suppressor of cytokine signaling 3 (Socs3). Additionally, Fos requirement for efficient long-term proliferation was demonstrated by the reduction of NSC clones capable of long-term expansion following CRISPR/Cas9-mediated Fos inactivation. Previous work showed that the Socs3 gene is strongly downregulated following Sox2 deletion, and its re-expression by lentiviral transduction rescues long-term NSC proliferation. Fos appears to be an upstream regulator of Socs3, possibly together with Jun and Egr2; indeed, Sox2 re-expression in Sox2-deleted NSC progressively activates both Fos and Socs3 expression; in turn, Fos transduction activates Socs3 expression. Based on available SOX2 ChIPseq and ChIA-PET data, we propose a model whereby Sox2 is a direct activator of both Socs3 and Fos, as well as possibly Jun and Egr2; furthermore, we provide direct evidence for FOS and JUN binding on Socs3 promoter, suggesting direct transcriptional regulation. These results provide the basis for developing a model of a network of interactions, regulating critical effectors of NSC proliferation and long-term maintenance.

Oncogenic extrachromosomal DNA functions as mobile enhancers to globally amplify chromosomal transcription

Extrachromosomal, circular DNA (ecDNA) is emerging as a prevalent yet less characterized oncogenic alteration in cancer genomes. We leverage ChIA-PET and ChIA-Drop chromatin interaction assays to characterize genome-wide ecDNA-mediated chromatin contacts that impact transcriptional programs in cancers. ecDNAs in glioblastoma patient-derived neurosphere and prostate cancer cell cultures are marked by widespread intra-ecDNA and genome-wide chromosomal interactions. ecDNA-chromatin contact foci are characterized by broad and high-level H3K27ac signals converging predominantly on chromosomal genes of increased expression levels. Prostate cancer cells harboring synthetic ecDNA circles composed of characterized enhancers result in the genome-wide activation of chromosomal gene transcription. Deciphering the chromosomal targets of ecDNAs at single-molecule resolution reveals an association with actively expressed oncogenes spatially clustered within ecDNA-directed interaction networks. Our results suggest that ecDNA can function as mobile transcriptional enhancers to promote tumor progression and manifest a potential synthetic aneuploidy mechanism of transcription control in cancer.

[Evaluation of intervention effect in the occupational protection of glass fiber workers by occupational health education]

Objective: To understand the mastery of occupational hygienic knowledge and the implementing of occupational health protection measures in the group which were exposed to the procedure of manufacture and use in glass fiber company, and to explore the feasibility of the prevention of the skin injury by occupational health education in glass fiber workers. Methods: We selected 257 on-the-job employees as the research object in a ceramic enterprise in Nanjing from June 2018 to August 2019, with the method of cluster random sampling. According to Solomon's design, the intervention group in which we took measures with health education was divided into RG(1) (O(1)XO(2)) and RG(3) (XO(5)) group, and the control group where we didn't take any intervention was divided into RG(2) (O(3)-O(4)) and RG(4) (-O(6)) group. The intervention effect of health education on the occupational protection of glass fiber workers was evaluated by the results of questionnaire. Results: After training, the average score of occupational health knowledge in the intervention group was 27.34 points higher than that before training, the intervention index was 1.42, 23.62-27.73 points higher than the control glass fiber workers and 33.62-35.52 points higher than the control glass non-glass fiber workers; Compared with the control group, the positive attitude rate of fiber glass workers in the intervention group increased by 13.28%, 13.51%, 11.68% and 11.48%, and the intervention indexes were 1.18, 1.17, 1.14 and 1.15, which was corresponding to using protective cream, wearing gloves, wearing working clothes, washing hands and bathing after work, respectively; Compared with the control group, the implementation rate of occupational protection measures which were represented by wearing gloves、washing hands and bathing for glass fiber workers in the intervention group increased by 29.25% and 7.27% respectively, and the intervention indexes were 1.43 and 1.08 respectively; The skin injury rate of fiberglass workers in the intervention group was reduced by 11.43% comparing to the control group, the intervention index was 1.67. Conclusion: According to the occupational health education of fiberglass workers, it improves the mastery of occupational health knowledge, positive atti-tude rate and the implementation rate of occupational protection measures, meanwhile, it reduces the skin injury rate of the intervention objects to a certain extent.

ChIA-PIPE: A fully automated pipeline for comprehensive ChIA-PET data analysis and visualization

ChIA-PET (chromatin interaction analysis with paired-end tags) enables genome-wide discovery of chromatin interactions involving specific protein factors, with base pair resolution. Interpretation of ChIA-PET data requires a robust analytic pipeline. Here, we introduce ChIA-PIPE, a fully automated pipeline for ChIA-PET data processing, quality assessment, visualization, and analysis. ChIA-PIPE performs linker filtering, read mapping, peak calling, and loop calling and automates quality control assessment for each dataset. To enable visualization, ChIA-PIPE generates input files for two-dimensional contact map viewing with Juicebox and HiGlass and provides a new dockerized visualization tool for high-resolution, browser-based exploration of peaks and loops. To enable structural interpretation, ChIA-PIPE calls chromatin contact domains, resolves allele-specific peaks and loops, and annotates enhancer-promoter loops. ChIA-PIPE also supports the analysis of other related chromatin-mapping data types.

Chromatin topology reorganization and transcription repression by PML-RARα in acute promyeloid leukemia

BACKGROUND

Acute promyeloid leukemia (APL) is characterized by the oncogenic fusion protein PML-RARα, a major etiological agent in APL. However, the molecular mechanisms underlying the role of PML-RARα in leukemogenesis remain largely unknown.

RESULTS

Using an inducible system, we comprehensively analyze the 3D genome organization in myeloid cells and its reorganization after PML-RARα induction and perform additional analyses in patient-derived APL cells with native PML-RARα. We discover that PML-RARα mediates extensive chromatin interactions genome-wide. Globally, it redefines the chromatin topology of the myeloid genome toward a more condensed configuration in APL cells; locally, it intrudes RNAPII-associated interaction domains, interrupts myeloid-specific transcription factors binding at enhancers and super-enhancers, and leads to transcriptional repression of genes critical for myeloid differentiation and maturation.

CONCLUSIONS

Our results not only provide novel topological insights for the roles of PML-RARα in transforming myeloid cells into leukemia cells, but further uncover a topological framework of a molecular mechanism for oncogenic fusion proteins in cancers.

Mapping the Global Chromatin Connectivity Network for Sox2 Function in Neural Stem Cell Maintenance

The SOX2 transcription factor is critical for neural stem cell (NSC) maintenance and brain development. Through chromatin immunoprecipitation (ChIP) and chromatin interaction analysis (ChIA-PET), we determined genome-wide SOX2-bound regions and Pol II-mediated long-range chromatin interactions in brain-derived NSCs. SOX2-bound DNA was highly enriched in distal chromatin regions interacting with promoters and carrying epigenetic enhancer marks. Sox2 deletion caused widespread reduction of Pol II-mediated long-range interactions and decreased gene expression. Genes showing reduced expression in Sox2-deleted cells were significantly enriched in interactions between promoters and SOX2-bound distal enhancers. Expression of one such gene, Suppressor of Cytokine Signaling 3 (Socs3), rescued the self-renewal defect of Sox2-ablated NSCs. Our work identifies SOX2 as a major regulator of gene expression through connections to the enhancer network in NSCs. Through the definition of such a connectivity network, our study shows the way to the identification of genes and enhancers involved in NSC maintenance and neurodevelopmental disorders.

•

Sox2-bound enhancers are enriched within long-range interactions in neural stem cells

•

SOX2 loss decreases chromatin interactivity genome-wide

•

Sox2-bound enhancers from interactions activate reporter genes in zebrafish forebrain

•

Socs3, a gene downregulated in Sox2 mutant NSCs, rescues their self-renewal

Sox2-Dependent 3D Chromatin Interactomes in Transcription, Neural Stem Cell Proliferation and Neurodevelopmental Diseases

In our article, we asked whether Sox2, a transcription factor important in brain development and disease, is involved in gene regulation through its action on long-range interactions between promoters and distant enhancers. Our findings highlight that Sox2 shapes a genome-wide network of promoter-enhancer interactions, acting by direct binding to these elements. Sox2 loss affects the three-dimensional (3D) genome and decreases the activity of a subset of genes involved in Sox2-bound interactions. At least one of such downregulated genes, Socs3, is critical for long-term neural stem cell maintenance. These results point to the possibility of identifying a transcriptional network downstream to Sox2, and involved in neural stem cell maintenance. In addition, interacting Sox2-bound enhancers are often connected to genes which are relevant, in man, to neurodevelopmental disease; this may facilitate the detection of functionally relevant mutations in regulatory elements in man, contributing to neural disease.

Ultra-long Read Sequencing for Whole Genomic DNA Analysis

Third generation single-molecule DNA sequencing technologies offer significantly longer read length that can facilitate the assembly of complex genomes and analysis of complex structural variants. Nanopore platforms perform single-molecule sequencing by directly measuring the current changes mediated by DNA passage through the pores and can generate hundreds of kilobase (kb) reads with minimal capital cost. This platform has been adopted by many researchers for a variety of applications. Achieving longer sequencing read lengths is the most critical factor to leverage the value of nanopore sequencing platforms. To generate ultra-long reads, special consideration is required to avoid DNA breakages and gain efficiency to generate productive sequencing templates. Here, we provide the detailed protocol of ultra-long DNA sequencing including high molecular weight (HMW) DNA extraction from fresh or frozen cells, library construction by mechanical shearing or transposase fragmentation, and sequencing on a nanopore device. From 20-25 µg of HMW DNA, the method can achieve N50 read length of 50-70 kb with mechanical shearing and N50 of 90-100 kb read length with transposase mediated fragmentation. The protocol can be applied to DNA extracted from mammalian cells to perform whole genome sequencing for the detection of structural variants and genome assembly. Additional improvements on the DNA extraction and enzymatic reactions will further increase the read length and expand its utility.

Complex Microbiome in Brain Abscess Revealed by Whole-Genome Culture-Independent and Culture-Based Sequencing

Brain abscess is a severe infectious disease with high mortality and mobility. Although culture-based techniques have been widely used for the investigation of microbial composition of brain abscess, these approaches are inherent biased. Recent studies using 16S ribosomal sequencing approaches revealed high complexity of the bacterial community involved in brain abscess but fail to detect fungal and viral composition. In the study, both culture-independent nanopore metagenomic sequencing and culture-based whole-genome sequencing using both the Illumina and the Nanopore platforms were conducted to investigate the microbial composition and genomic characterization in brain abscess. Culture-independent metagenomic sequencing revealed not only a larger taxonomic diversity of bacteria but also the presence of fungi and virus communities. The culture-based whole-genome sequencing identified a novel species in Prevotella and reconstructs a Streptococcus constellatus with a high GC-skew genome. Antibiotic-resistance genes CfxA and ErmF associated with resistance to penicillin and clindamycin were also identified in culture-based and culture-free sequencing. This study implies current understanding of brain abscess need to consider the broader diversity of microorganisms.

Multiplex chromatin interactions with single-molecule precision

The genomes of multicellular organisms are extensively folded into 3D chromosome territories within the nucleus¹. Advanced 3D genome-mapping methods that combine proximity ligation and high-throughput sequencing (such as chromosome conformation capture, Hi-C)², and chromatin immunoprecipitation techniques (such as chromatin interaction analysis by paired-end tag sequencing, ChIA-PET)³, have revealed topologically associating domains⁴ with frequent chromatin contacts, and have identified chromatin loops mediated by specific protein factors for insulation and regulation of transcription^5-7. However, these methods rely on pairwise proximity ligation and reflect population-level views, and thus cannot reveal the detailed nature of chromatin interactions. Although single-cell Hi-C⁸ potentially overcomes this issue, this method may be limited by the sparsity of data that is inherent to current single-cell assays. Recent advances in microfluidics have opened opportunities for droplet-based genomic analysis⁹ but this approach has not yet been adapted for chromatin interaction analysis. Here we describe a strategy for multiplex chromatin-interaction analysis via droplet-based and barcode-linked sequencing, which we name ChIA-Drop. We demonstrate the robustness of ChIA-Drop in capturing complex chromatin interactions with single-molecule precision, which has not been possible using methods based on population-level pairwise contacts. By applying ChIA-Drop to Drosophila cells, we show that chromatin topological structures predominantly consist of multiplex chromatin interactions with high heterogeneity; ChIA-Drop also reveals promoter-centred multivalent interactions, which provide topological insights into transcription.

Transcriptional responses of the marine green alga Micromonas pusilla and an infecting prasinovirus under different phosphate conditions

Prasinophytes are widespread marine algae for which responses to nutrient limitation and viral infection are not well understood. We studied the picoprasinophyte, Micromonas pusilla, grown under phosphate-replete (0.65 ± 0.07 d^-1 ) and 10-fold lower (low)-phosphate (0.11 ± 0.04 d^-1 ) conditions, and infected by the phycodnavirus MpV-SP1. Expression of 17% of Micromonas genes in uninfected cells differed by >1.5-fold (q < 0.01) between nutrient conditions, with genes for P-metabolism and the uniquely-enriched Sel1-like repeat (SLR) family having higher relative transcript abundances, while phospholipid-synthesis genes were lower in low-P than P-replete. Approximately 70% (P-replete) and 30% (low-P) of cells were lysed 24 h post-infection, and expression of ≤5.8% of host genes changed relative to uninfected treatments. Host genes for CAZymes and glycolysis were activated by infection, supporting importance in viral production, which was significantly lower in slower growing (low-P) hosts. All MpV-SP1 genes were expressed, and our analyses suggest responses to differing host-phosphate backgrounds involve few viral genes, while the temporal program of infection involves many more, and is largely independent of host-phosphate background. Our study (i) identifies genes previously unassociated with nutrient acclimation or viral infection, (ii) provides insights into the temporal program of prasinovirus gene expression by hosts and (iii) establishes cell biological aspects of an ecologically important host-prasinovirus system that differ from other marine algal-virus systems.

Specialized proteomic responses and an ancient photoprotection mechanism sustain marine green algal growth during phosphate limitation

Marine algae perform approximately half of global carbon fixation, but their growth is often limited by the availability of phosphate or other nutrients^1,2. As oceans warm, the area of phosphate-limited surface waters is predicted to increase, resulting in ocean desertification^3,4. Understanding the responses of key eukaryotic phytoplankton to nutrient limitation is therefore critical^5,6. We used advanced photo-bioreactors to investigate how the widespread marine green alga Micromonas commoda grows under transitions from replete nutrients to chronic phosphate limitation and subsequent relief, analysing photosystem changes and broad cellular responses using proteomics, transcriptomics and biophysical measurements. We find that physiological and protein expression responses previously attributed to stress are critical to supporting stable exponential growth when phosphate is limiting. Unexpectedly, the abundance of most proteins involved in light harvesting does not change, but an ancient light-harvesting-related protein, LHCSR, is induced and dissipates damaging excess absorbed light as heat throughout phosphate limitation. Concurrently, a suite of uncharacterized proteins with narrow phylogenetic distributions increase multifold. Notably, of the proteins that exhibit significant changes, 70% are not differentially expressed at the mRNA transcript level, highlighting the importance of post-transcriptional processes in microbial eukaryotes. Nevertheless, transcript-protein pairs with concordant changes were identified that will enable more robust interpretation of eukaryotic phytoplankton responses in the field from metatranscriptomic studies. Our results show that P-limited Micromonas responds quickly to a fresh pulse of phosphate by rapidly increasing replication, and that the protein network associated with this ability is composed of both conserved and phylogenetically recent proteome systems that promote dynamic phosphate homeostasis. That an ancient mechanism for mitigating light stress is central to sustaining growth during extended phosphate limitation highlights the possibility of interactive effects arising from combined stressors under ocean change, which could reduce the efficacy of algal strategies for optimizing marine photosynthesis.

Picky comprehensively detects high-resolution structural variants in nanopore long reads

Acquired genomic structural variants (SVs) are major hallmarks of cancer genomes, but they are challenging to reconstruct from short-read sequencing data. Here we exploited the long reads of the nanopore platform using our customized pipeline, Picky ( https://github.com/TheJacksonLaboratory/Picky ), to reveal SVs of diverse architecture in a breast cancer model. We identified the full spectrum of SVs with superior specificity and sensitivity relative to short-read analyses, and uncovered repetitive DNA as the major source of variation. Examination of genome-wide breakpoints at nucleotide resolution uncovered micro-insertions as the common structural features associated with SVs. Breakpoint density across the genome is associated with the propensity for interchromosomal connectivity and was found to be enriched in promoters and transcribed regions of the genome. Furthermore, we observed an over-representation of reciprocal translocations from chromosomal double-crossovers through phased SVs. We demonstrate that Picky analysis is an effective tool for comprehensive detection of SVs in cancer genomes from long-read data.

Sox2 conditional mutation in mouse causes ataxic symptoms, cerebellar vermis hypoplasia, and postnatal defects of Bergmann glia

Sox2 is a transcription factor active in the nervous system, within different cell types, ranging from radial glia neural stem cells to a few specific types of differentiated glia and neurons. Mutations in the human SOX2 transcription factor gene cause various central nervous system (CNS) abnormalities, involving hippocampus and eye defects, as well as ataxia. Conditional Sox2 mutation in mouse, with different Cre transgenes, previously recapitulated different essential features of the disease, such as hippocampus and eye defects. In the cerebellum, Sox2 is active from early embryogenesis in the neural progenitors of the cerebellar primordium; Sox2 expression is maintained, postnatally, within Bergmann glia (BG), a differentiated cell type essential for Purkinje neurons functionality and correct motor control. By performing Sox2 Cre-mediated ablation in the developing and postnatal mouse cerebellum, we reproduced ataxia features. Embryonic Sox2 deletion (with Wnt1Cre) leads to reduction of the cerebellar vermis, known to be commonly related to ataxia, preceded by deregulation of Otx2 and Gbx2, critical regulators of vermis development. Postnatally, BG is progressively disorganized, mislocalized, and reduced in mutants. Sox2 postnatal deletion, specifically induced in glia (with GLAST-CreERT2), reproduces the BG defect, and causes (milder) ataxic features. Our results define a role for Sox2 in cerebellar function and development, and identify a functional requirement for Sox2 within postnatal BG, of potential relevance for ataxia in mouse mutants, and in human patients.

An improved DNA marker technique for genetic characterization using RAMP-PCR with high-GC primers

Random amplified polymorphic DNA (RAPD) is a widely used molecular marker technique. As traditional RAPD has poor reproducibility and productivity, we previously developed an improved RAPD method (termed RAMP-PCR), which increased the reproducibility, number of bands, and efficiency of studies on polymorphism. To further develop the efficiency of this method, we used high-GC content primers for improved RAMP-PCR with DNA samples from Lonicera japonica. Comparison of amplification profiles obtained by standard RAPD primers with those obtained by regular PCR and RAMP-PCR, and high-GC primers with regular PCR and RAMP-PCR showed that the average number of bands and polymorphisms per primer gradually and significantly increased (from 6.4 to 15.0 and from 4.6 to 10.2, respectively). Cluster dendrograms showed similar results, indicating that this new method is consistent and reproducible. A total of 22 samples from different species, including plants, animals, and humans, were used for RAMP-PCR with high-GC primers. Multiple bands were successfully amplified from all samples, demonstrating that this method is a reliable technique with consistent results and may be of general interest in studies on different genera and species. We developed highly effective DNA markers, which can provide a more effective and potentially valuable approach than traditional RAPD for the genetic identification of various organisms, particularly of medicinal plants.

Development and significance of SCAR marker QG12-5 for Canarium album (Lour.) Raeusch by molecular cloning from improved RAPD amplification

Sequence-characterized amplified region (SCAR) is a valuable molecular marker for the genetic identification of any species. This marker is mainly derived from molecular cloning of random amplified polymorphic DNA (RAPD). We have previously reported the use of an improved RAPD technique for the genetic characterization of different samples of Canarium album (Lour.) Raeusch (C. album). In this study, DNA fragments were amplified using improved RAPD amplified from different samples of C. album. The amplified DNA fragment was excised, purified from an agarose gel and cloned into a pGM-T vector; subsequently, a positive clone, called QG12-5 was identified by PCR amplification and enzymatic digestion and sequenced by Sanger di-deoxy sequencing method. This clone was revealed consisting of 510 nucleotides of C. album. The SCAR marker QG12-5 was developed using specifically designed PCR primers and optimized PCR conditions. This SCAR marker expressed seven continuous "TATG" [(TATG)n] tandem repeats, which was found to characterize C. album. Subsequently, this novel SCAR marker was deposited in GenBank with accession No. KT359568. Therefore, we successfully developed a C. album-specific SCAR marker for the identification and authentication of different C. album species in this study.

Development of two novel specific SCAR markers by cloning improved RAPD fragments from the medicinal mushroom Ganoderma lucidium (Leysser: Fr) Karst

Development of sequence-characterized amplified region (SCAR) markers from random-amplified polymorphic DNA (RAPD) fragments is a valuable molecular approach for the genetic identification of different species. By using SCAR markers, molecular analysis is reduced to a simple polymerase chain reaction (PCR) analysis using primers designed from the amplicon sequence of RAPD. In this study, the DNA fragments from an improved RAPD amplification of Ganoderma species were cloned into a pGM-T vector; positive clones were identified by PCR amplification and enzymatic digestion, and finally, DNA fragments were sequenced using the Sanger sequencing method for developing the SCAR markers. Two SCAR markers, named LZ4-1 with 534 nucleotides, and LZ5-2 with 337 nucleotides were identified, which are specific to Ganoderma lucidium (Leysser: Fr) Karst species. BLAST of these two nucleotide sequences in the GenBank database showed no identity to other species. We deposited these sequences into the GenBank database (LZ4-1 accession No. KM391933, LZ5-2 accession No. KM391934). PCR assays confirmed them as novel molecular markers for G. lucidium (Leysser: Fr) Karst, which might be used for genetic authentication of adulterant samples. Thus, our study developed two specific SCAR markers for identifying and distinguishing the medicinal mushroom G. lucidium (Leysser: Fr) Karst from other Ganoderma species.

Identifying Aspects of the Post-Transcriptional Program Governing the Proteome of the Green Alga Micromonas pusilla

Micromonas is a unicellular motile alga within the Prasinophyceae, a green algal group that is related to land plants. This picoeukaryote (<2 μm diameter) is widespread in the marine environment but is not well understood at the cellular level. Here, we examine shifts in mRNA and protein expression over the course of the day-night cycle using triplicated mid-exponential, nutrient replete cultures of Micromonas pusilla CCMP1545. Samples were collected at key transition points during the diel cycle for evaluation using high-throughput LC-MS proteomics. In conjunction, matched mRNA samples from the same time points were sequenced using pair-ended directional Illumina RNA-Seq to investigate the dynamics and relationship between the mRNA and protein expression programs of M. pusilla. Similar to a prior study of the marine cyanobacterium Prochlorococcus, we found significant divergence in the mRNA and proteomics expression dynamics in response to the light:dark cycle. Additionally, expressional responses of genes and the proteins they encoded could also be variable within the same metabolic pathway, such as we observed in the oxygenic photosynthesis pathway. A regression framework was used to predict protein levels from both mRNA expression and gene-specific sequence-based features. Several features in the genome sequence were found to influence protein abundance including codon usage as well as 3’ UTR length and structure. Collectively, our studies provide insights into the regulation of the proteome over a diel cycle as well as the relationships between transcriptional and translational programs in the widespread marine green alga Micromonas.

[Retrospective cohort study for the impact on readmission of patients with ischemic stroke after treatment of aspirin plus clopidogrel or aspirin mono-therapy]

OBJECTIVE

To see the influence of different antiplatelet therapies on stroke patients' readmission by performing a deep data-mining into Beijing Healthcare Insuring Database, based on a large sample size.

METHODS

Aretrospective cohort study, was adopted to extract patients primarily diagnosed as ischemic stroke from healthcare database. The first hospital records were considered as the patient's baseline in this study, who were divided into MAPT (aspirin) and DAPT (aspirin and clopidogrel) according to the patient's baseline medications. A follow-up was conducted to see whether the patients would have rehospitalization record because of major result events after medication. The major result events, included: (1) recurrence of ischemic stroke; (2) hemorrhagic transformation of ischemic stroke; (3) myocardial infarction; (4) the digestive hemorrhage. The Kaplan-Meier figure was used to compare the survival situations between these two groups, the log-rank test was used to test the difference of the survival curve, and 1:1 propensity score matching was calculated from the patients' baseline data. Cox proportional hazards model was used to calculate the hazard ratio (HR).

RESULTS

A total of 27 695 patients From January 2010 to September 2013 were included, 4 047 with DAPT, and 23 648 with MAPT. Because the baseline characteristics of the patients was disequilibrium, so we used 1:1 propensity score matching, after which, the number of the two groups was 4 046 each. Adjusted for the general demographic characteristics such as age, sex, nationality, complication and drug combination, no statistical significance was observed between the survival curves of the two groups (P=0.06). HR value of major result events between the groups was 0.91 (0.82-1.01, P=0.07), which was not statistically significant. The covariate gender HR=1.36 (1.20-1.55, P<0.05), accompanied by diabetes HR= 1.36 (1.20-1.54, P<0.05), dyslipidemia HR=1.13 (1.00-1.27, P=1.13), heart disease HR=1.39 (1.22-1.58, P<0.05) was statistically significant. Drug combination with other antiplatelet agents HR=1.05 (0.95-1.17, P>1.05) did not increase the risk of readmission.

CONCLUSION

There was no difference in prevention of readmission between patients with DAPT and MAPT. Patients with complications should actively treat the complications at the same time as they prevent recurrence after first attack.

Development of novel SCAR markers for genetic characterization of Lonicera japonica from high GC-RAMP-PCR and DNA cloning

Sequence-characterized amplified region (SCAR) markers were further developed from high-GC primer RAMP-PCR-amplified fragments from Lonicera japonica DNA by molecular cloning. The four DNA fragments from three high-GC primers (FY-27, FY-28, and FY-29) were successfully cloned into a pGM-T vector. The positive clones were sequenced; their names, sizes, and GenBank numbers were JYHGC1-1, 345 bp, KJ620024; YJHGC2-1, 388 bp, KJ620025; JYHGC7-2, 1036 bp, KJ620026; and JYHGC6-2, 715 bp, KJ620027, respectively. Four novel SCAR markers were developed by designing specific primers, optimizing conditions, and PCR validation. The developed SCAR markers were used for the genetic authentication of L. japonica from its substitutes. This technique provides another means of developing DNA markers for the characterization and authentication of various organisms including medicinal plants and their substitutes.

Evidence-based green algal genomics reveals marine diversity and ancestral characteristics of land plants

Background

Prasinophytes are widespread marine green algae that are related to plants. Cellular abundance of the prasinophyte Micromonas has reportedly increased in the Arctic due to climate-induced changes. Thus, studies of these unicellular eukaryotes are important for marine ecology and for understanding Viridiplantae evolution and diversification.

Results

We generated evidence-based Micromonas gene models using proteomics and RNA-Seq to improve prasinophyte genomic resources. First, sequences of four chromosomes in the 22 Mb Micromonas pusilla (CCMP1545) genome were finished. Comparison with the finished 21 Mb genome of Micromonas commoda (RCC299; named herein) shows they share ≤8,141 of ~10,000 protein-encoding genes, depending on the analysis method. Unlike RCC299 and other sequenced eukaryotes, CCMP1545 has two abundant repetitive intron types and a high percent (26 %) GC splice donors. Micromonas has more genus-specific protein families (19 %) than other genome sequenced prasinophytes (11 %). Comparative analyses using predicted proteomes from other prasinophytes reveal proteins likely related to scale formation and ancestral photosynthesis. Our studies also indicate that peptidoglycan (PG) biosynthesis enzymes have been lost in multiple independent events in select prasinophytes and plants. However, CCMP1545, polar Micromonas CCMP2099 and prasinophytes from other classes retain the entire PG pathway, like moss and glaucophyte algae. Surprisingly, multiple vascular plants also have the PG pathway, except the Penicillin-Binding Protein, and share a unique bi-domain protein potentially associated with the pathway. Alongside Micromonas experiments using antibiotics that halt bacterial PG biosynthesis, the findings highlight unrecognized phylogenetic complexity in PG-pathway retention and implicate a role in chloroplast structure or division in several extant Viridiplantae lineages.

Conclusions

Extensive differences in gene loss and architecture between related prasinophytes underscore their divergence. PG biosynthesis genes from the cyanobacterial endosymbiont that became the plastid, have been selectively retained in multiple plants and algae, implying a biological function. Our studies provide robust genomic resources for emerging model algae, advancing knowledge of marine phytoplankton and plant evolution.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2585-6) contains supplementary material, which is available to authorized users.

Isolation and analysis of cell-free fetal DNA from maternal peripheral blood in Chinese women

Non-invasive prenatal diagnosis is used to detect the genetic material of the fetus by isolating the cell-free fetal DNA (cffDNA) from maternal peripheral blood. In order to establish an isolation method for cffDNA from maternal peripheral blood in Chinese women, the cffDNA was acquired with a two-step centrifugation using a QlAamp DNA Blood mini kit. The SRY gene of plasma DNA was amplified by polymerase chain reaction (PCR). Real-time quantitative PCR was used to measure the concentration of cffDNA in maternal peripheral blood in different pregnant women. The results of the SRY gene amplification of plasma DNA from pregnant women was the same as that of the amniocyte DNA. The average concentration of cffDNA in maternal peripheral blood of pregnant women in different gestational stages was 0.98 ng/mL (0.26-1.49 ng/mL), 1.43 ng/mL (0.46- 2.34 ng/mL), and 1.95 ng/mL (0.65-6.81 ng/mL) from early, middle, and late gestational stages, respectively. The mean of cffDNA from total DNA in plasma in different stages of gestation was 22.28% (9.86-27.81%). The lowest concentration of DNA amplified by nested-PCR in our research was 10-4-10-3 ng/μL. The isolation method for cffDNA from maternal peripheral blood was successfully established and further research into its applications will be conducted.

CTCF-Mediated Human 3D Genome Architecture Reveals Chromatin Topology for Transcription

Spatial genome organization and its effect on transcription remains a fundamental question. We applied an advanced chromatin interaction analysis by paired-end tag sequencing (ChIA-PET) strategy to comprehensively map higher-order chromosome folding and specific chromatin interactions mediated by CCCTC-binding factor (CTCF) and RNA polymerase II (RNAPII) with haplotype specificity and nucleotide resolution in different human cell lineages. We find that CTCF/cohesin-mediated interaction anchors serve as structural foci for spatial organization of constitutive genes concordant with CTCF-motif orientation, whereas RNAPII interacts within these structures by selectively drawing cell-type-specific genes toward CTCF foci for coordinated transcription. Furthermore, we show that haplotype variants and allelic interactions have differential effects on chromosome configuration, influencing gene expression, and may provide mechanistic insights into functions associated with disease susceptibility. 3D genome simulation suggests a model of chromatin folding around chromosomal axes, where CTCF is involved in defining the interface between condensed and open compartments for structural regulation. Our 3D genome strategy thus provides unique insights in the topological mechanism of human variations and diseases.

Identification of the origin of marker chromosomes by two-color fluorescence in situ hybridization and polymerase chain reaction in azoospermic patients

Y chromosomal microdeletions at the azoospermia factor locus and chromosome abnormalities have been implicated as the major causes of idiopathic male infertility. A marker chromosome is a structurally abnormal chromosome in which no part can be identified by cytogenetics. In this study, to identify the origin of the marker chromosomes and to perform a genetic diagnosis of patients with azoospermia, two-color fluorescence in situ hybridization (FISH) and polymerase chain reaction (PCR) techniques were carried out. The marker chromosomes for the two patients with azoospermia originated in the Y chromosome; it was ascertained that the karyotype of both patients was 46,X, ish del(Y)(q11)(DYZ3+, DXZ1-). The combination of two-color FISH and PCR techniques is an important method for the identification of the origin of marker chromosomes. Thus, genetic counseling and a clear genetic diagnosis of patients with azoospermia before intracytoplasmic sperm injection or other clinical managements are important.

Diversity and population structure of northern switchgrass as revealed through exome capture sequencing

Panicum virgatum L. (switchgrass) is a polyploid, perennial grass species that is native to North America, and is being developed as a future biofuel feedstock crop. Switchgrass is present primarily in two ecotypes: a northern upland ecotype, composed of tetraploid and octoploid accessions, and a southern lowland ecotype, composed of primarily tetraploid accessions. We employed high-coverage exome capture sequencing (~2.4 Tb) to genotype 537 individuals from 45 upland and 21 lowland populations. From these data, we identified ~27 million single-nucleotide polymorphisms (SNPs), of which 1 590 653 high-confidence SNPs were used in downstream analyses of diversity within and between the populations. From the 66 populations, we identified five primary population groups within the upland and lowland ecotypes, a result that was further supported through genetic distance analysis. We identified conserved, ecotype-restricted, non-synonymous SNPs that are predicted to affect the protein function of CONSTANS (CO) and EARLY HEADING DATE 1 (EHD1), key genes involved in flowering, which may contribute to the phenotypic differences between the two ecotypes. We also identified, relative to the near-reference Kanlow population, 17 228 genes present in more copies than in the reference genome (up-CNVs), 112 630 genes present in fewer copies than in the reference genome (down-CNVs) and 14 430 presence/absence variants (PAVs), affecting a total of 9979 genes, including two upland-specific CNV clusters. In total, 45 719 genes were affected by an SNP, CNV, or PAV across the panel, providing a firm foundation to identify functional variation associated with phenotypic traits of interest for biofuel feedstock production.

Lineage-specific chromatin signatures reveal a regulator of lipid metabolism in microalgae

Alga-derived lipids represent an attractive potential source of biofuels. However, lipid accumulation in algae is a stress response tightly coupled to growth arrest, thereby imposing a major limitation on productivity. To identify transcriptional regulators of lipid accumulation, we performed an integrative chromatin signature and transcriptomic analysis to decipher the regulation of lipid biosynthesis in the alga Chlamydomonas reinhardtii. Genome-wide histone modification profiling revealed remarkable differences in functional chromatin states between the algae and higher eukaryotes and uncovered regulatory components at the core of lipid accumulation pathways. We identified the transcription factor, PSR1, as a pivotal switch that triggers cytosolic lipid accumulation. Dissection of the PSR1-induced lipid profiles corroborates its role in coordinating multiple lipid-inducing stress responses. The comprehensive maps of functional chromatin signatures in a major clade of eukaryotic life and the discovery of a transcriptional regulator of algal lipid metabolism will facilitate targeted engineering strategies to mediate high lipid production in microalgae.

A genome-wide map of adeno-associated virus-mediated human gene targeting

To determine which genomic features promote homologous recombination, we created a genome-wide map of gene targeting sites. We used an adeno-associated virus vector to target identical loci introduced as transcriptionally active retroviral vectors. A comparison of ~2,000 targeted and untargeted sites showed that targeting occurred throughout the human genome and was not influenced by the presence of nearby CpG islands, sequence repeats or DNase I-hypersensitive sites. Targeted sites were preferentially located within transcription units, especially when the target loci were transcribed in the opposite orientation to their surrounding chromosomal genes. We determined the impact of DNA replication by mapping replication forks, which revealed a preference for recombination at target loci transcribed toward an incoming fork. Our results constitute the first genome-wide screen of gene targeting in mammalian cells and demonstrate a strong recombinogenic effect of colliding polymerases.

Activation and repression by oncogenic MYC shape tumour-specific gene expression profiles

In mammalian cells, the MYC oncoprotein binds to thousands of promoters. During mitogenic stimulation of primary lymphocytes, MYC promotes an increase in the expression of virtually all genes. In contrast, MYC-driven tumour cells differ from normal cells in the expression of specific sets of up- and downregulated genes that have considerable prognostic value. To understand this discrepancy, we studied the consequences of inducible expression and depletion of MYC in human cells and murine tumour models. Changes in MYC levels activate and repress specific sets of direct target genes that are characteristic of MYC-transformed tumour cells. Three factors account for this specificity. First, the magnitude of response parallels the change in occupancy by MYC at each promoter. Functionally distinct classes of target genes differ in the E-box sequence bound by MYC, suggesting that different cellular responses to physiological and oncogenic MYC levels are controlled by promoter affinity. Second, MYC both positively and negatively affects transcription initiation independent of its effect on transcriptional elongation. Third, complex formation with MIZ1 (also known as ZBTB17) mediates repression of multiple target genes by MYC and the ratio of MYC and MIZ1 bound to each promoter correlates with the direction of response.

Genome sequence of the Listia angolensis microsymbiont Microvirga lotononidis strain WSM3557(T.)

Microvirga lotononidis is a recently described species of root-nodule bacteria that is an effective nitrogen- (N2) fixing microsymbiont of the symbiotically specific African legume Listia angolensis (Welw. ex Bak.) B.-E. van Wyk & Boatwr. M. lotononidis possesses several properties that are unusual in root-nodule bacteria, including pigmentation and the ability to grow at temperatures of up to 45°C. Strain WSM3557(T) is an aerobic, motile, Gram-negative, non-spore-forming rod isolated from a L. angolensis root nodule collected in Chipata, Zambia in 1963. This is the first report of a complete genome sequence for the genus Microvirga. Here we describe the features of Microvirga lotononidis strain WSM3557(T), together with genome sequence information and annotation. The 7,082,538 high-quality-draft genome is arranged in 18 scaffolds of 104 contigs, contains 6,956 protein-coding genes and 84 RNA-only encoding genes, and is one of 20 rhizobial genomes sequenced as part of the DOE Joint Genome Institute 2010 Community Sequencing Program.

A near complete snapshot of the Zea mays seedling transcriptome revealed from ultra-deep sequencing

RNA-sequencing (RNA-seq) enables in-depth exploration of transcriptomes, but typical sequencing depth often limits its comprehensiveness. In this study, we generated nearly 3 billion RNA-Seq reads, totaling 341 Gb of sequence, from a Zea mays seedling sample. At this depth, a near complete snapshot of the transcriptome was observed consisting of over 90% of the annotated transcripts, including lowly expressed transcription factors. A novel hybrid strategy combining de novo and reference-based assemblies yielded a transcriptome consisting of 126,708 transcripts with 88% of expressed known genes assembled to full-length. We improved current annotations by adding 4,842 previously unannotated transcript variants and many new features, including 212 maize transcripts, 201 genes, 10 genes with undocumented potential roles in seedlings as well as maize lineage specific gene fusion events. We demonstrated the power of deep sequencing for large transcriptome studies by generating a high quality transcriptome, which provides a rich resource for the research community.

Interactome maps of mouse gene regulatory domains reveal basic principles of transcriptional regulation

A key finding of the ENCODE project is that the enhancer landscape of mammalian cells undergoes marked alterations during ontogeny. However, the nature and extent of these changes are unclear. As part of the NIH Mouse Regulome Project, we here combined DNaseI hypersensitivity, ChIP-seq, and ChIA-PET technologies to map the promoter-enhancer interactomes of pluripotent ES cells and differentiated B lymphocytes. We confirm that enhancer usage varies widely across tissues. Unexpectedly, we find that this feature extends to broadly transcribed genes, including Myc and Pim1 cell-cycle regulators, which associate with an entirely different set of enhancers in ES and B cells. By means of high-resolution CpG methylomes, genome editing, and digital footprinting, we show that these enhancers recruit lineage-determining factors. Furthermore, we demonstrate that the turning on and off of enhancers during development correlates with promoter activity. We propose that organisms rely on a dynamic enhancer landscape to control basic cellular functions in a tissue-specific manner.

Genome sequence of the Ornithopus/Lupinus-nodulating Bradyrhizobium sp. strain WSM471

Bradyrhizobium sp. strain WSM471 is an aerobic, motile, Gram-negative, non-spore-forming rod that was isolated from an effective nitrogen- (N₂) fixing root nodule formed on the annual legume Ornithopus pinnatus (Miller) Druce growing at Oyster Harbour, Albany district, Western Australia in 1982. This strain is in commercial production as an inoculant for Lupinus and Ornithopus. Here we describe the features of Bradyrhizobium sp. strain WSM471, together with genome sequence information and annotation. The 7,784,016 bp high-quality-draft genome is arranged in 1 scaffold of 2 contigs, contains 7,372 protein-coding genes and 58 RNA-only encoding genes, and is one of 20 rhizobial genomes sequenced as part of the DOE Joint Genome Institute 2010 Community Sequencing Program.

Genome sequence of the Lebeckia ambigua-nodulating "Burkholderia sprentiae" strain WSM5005(T.)

"Burkholderia sprentiae" strain WSM5005(T) is an aerobic, motile, Gram-negative, non-spore-forming rod that was isolated in Australia from an effective N2-fixing root nodule of Lebeckia ambigua collected in Klawer, Western Cape of South Africa, in October 2007. Here we describe the features of "Burkholderia sprentiae" strain WSM5005(T), together with the genome sequence and its annotation. The 7,761,063 bp high-quality-draft genome is arranged in 8 scaffolds of 236 contigs, contains 7,147 protein-coding genes and 76 RNA-only encoding genes, and is one of 20 rhizobial genomes sequenced as part of the DOE Joint Genome Institute 2010 Community Sequencing Program.

Genome sequence of the lupin-nodulating Bradyrhizobium sp. strain WSM1417

Bradyrhizobium sp. strain WSM1417 is an aerobic, motile, Gram-negative, non-spore-forming rod that was isolated from an effective nitrogen (N₂) fixing root nodule of Lupinus sp. collected in Papudo, Chile, in 1995. However, this microsymbiont is a poorly effective N₂ fixer with the legume host Lupinus angustifolius L.; a lupin species of considerable economic importance in both Chile and Australia. The symbiosis formed with L. angustifolius produces less than half of the dry matter achieved by the symbioses with commercial inoculant strains such as Bradyrhizobium sp. strain WSM471. Therefore, WSM1417 is an important candidate strain with which to investigate the genetics of effective N₂ fixation in the lupin-bradyrhizobia symbioses. Here we describe the features of Bradyrhizobium sp. strain WSM1417, together with genome sequence information and annotation. The 8,048,963 bp high-quality-draft genome is arranged in a single scaffold of 2 contigs, contains 7,695 protein-coding genes and 77 RNA-only encoding genes, and is one of 20 rhizobial genomes sequenced as part of the DOE Joint Genome Institute 2010 Community Sequencing Program.

Genome sequence of the Trifolium rueppellianum -nodulating Rhizobium leguminosarum bv. trifolii strain WSM2012

Rhizobium leguminosarum bv. trifolii WSM2012 (syn. MAR1468) is an aerobic, motile, Gram-negative, non-spore-forming rod that was isolated from an ineffective root nodule recovered from the roots of the annual clover Trifolium rueppellianum Fresen growing in Ethiopia. WSM2012 has a narrow, specialized host range for N2-fixation. Here we describe the features of R. leguminosarum bv. trifolii strain WSM2012, together with genome sequence information and annotation. The 7,180,565 bp high-quality-draft genome is arranged into 6 scaffolds of 68 contigs, contains 7,080 protein-coding genes and 86 RNA-only encoding genes, and is one of 20 rhizobial genomes sequenced as part of the DOE Joint Genome Institute 2010 Community Sequencing Program.

Genome sequence of the clover-nodulating Rhizobium leguminosarum bv. trifolii strain TA1

Rhizobium leguminosarum bv. trifolii strain TA1 is an aerobic, motile, Gram-negative, non-spore-forming rod that is an effective nitrogen fixing microsymbiont on the perennial clovers originating from Europe and the Mediterranean basin. TA1 however is ineffective with many annual and perennial clovers originating from Africa and America. Here we describe the features of R. leguminosarum bv. trifolii strain TA1, together with genome sequence information and annotation. The 8,618,824 bp high-quality-draft genome is arranged in a 6 scaffold of 32 contigs, contains 8,493 protein-coding genes and 83 RNA-only encoding genes, and is one of 20 rhizobial genomes sequenced as part of the DOE Joint Genome Institute 2010 Community Sequencing Program.

Genome sequence of the South American clover-nodulating Rhizobium leguminosarum bv. trifolii strain WSM597

Rhizobium leguminosarum bv. trifolii strain WSM597 is an aerobic, motile, Gram-negative, non-spore-forming rod isolated from a root nodule of the annual clover Trifolium pallidum L. growing at Glencoe Research Station near Tacuarembó, Uruguay. This strain is generally ineffective for nitrogen (N₂) fixation with clovers of Mediterranean, North American and African origin, but is effective on the South American perennial clover T. polymorphum Poir. Here we describe the features of R. leguminosarum bv. trifolii strain WSM597, together with genome sequence information and annotation. The 7,634,384 bp high-quality-draft genome is arranged in 2 scaffolds of 53 contigs, contains 7,394 protein-coding genes and 87 RNA-only encoding genes, and is one of 20 rhizobial genomes sequenced as part of the DOE Joint Genome Institute 2010 Community Sequencing Program.

Complete genome sequence of Dehalobacter restrictus PER-K23(T.)

Dehalobacter restrictus strain PER-K23 (DSM 9455) is the type strain of the species Dehalobacter restrictus. D. restrictus strain PER-K23 grows by organohalide respiration, coupling the oxidation of H₂ to the reductive dechlorination of tetra- or trichloroethene. Growth has not been observed with any other electron donor or acceptor, nor has fermentative growth been shown. Here we introduce the first full genome of a pure culture within the genus Dehalobacter. The 2,943,336 bp long genome contains 2,826 protein coding and 82 RNA genes, including 5 16S rRNA genes. Interestingly, the genome contains 25 predicted reductive dehalogenase genes, the majority of which appear to be full length. The reductive dehalogenase genes are mainly located in two clusters, suggesting a much larger potential for organohalide respiration than previously anticipated.

Draft Genome Sequence of Frankia sp. Strain BCU110501, a Nitrogen-Fixing Actinobacterium Isolated from Nodules of Discaria trinevis

Frankia forms a nitrogen-fixing symbiosis with actinorhizal plants. We report a draft genome sequence for Frankia sp. strain BCU110501, a nitrogen-fixing actinobacterium isolated from nodules of Discaria trinevis grown in the Patagonia region of Argentina.

Thermus oshimai JL-2 and T. thermophilus JL-18 genome analysis illuminates pathways for carbon, nitrogen, and sulfur cycling

The complete genomes of Thermus oshimai JL-2 and T. thermophilus JL-18 each consist of a circular chromosome, 2.07 Mb and 1.9 Mb, respectively, and two plasmids ranging from 0.27 Mb to 57.2 kb. Comparison of the T. thermophilus JL-18 chromosome with those from other strains of T. thermophilus revealed a high degree of synteny, whereas the megaplasmids from the same strains were highly plastic. The T. oshimai JL-2 chromosome and megaplasmids shared little or no synteny with other sequenced Thermus strains. Phylogenomic analyses using a concatenated set of conserved proteins confirmed the phylogenetic and taxonomic assignments based on 16S rRNA phylogenetics. Both chromosomes encode a complete glycolysis, tricarboxylic acid (TCA) cycle, and pentose phosphate pathway plus glucosidases, glycosidases, proteases, and peptidases, highlighting highly versatile heterotrophic capabilities. Megaplasmids of both strains contained a gene cluster encoding enzymes predicted to catalyze the sequential reduction of nitrate to nitrous oxide; however, the nitrous oxide reductase required for the terminal step in denitrification was absent, consistent with their incomplete denitrification phenotypes. A sox gene cluster was identified in both chromosomes, suggesting a mode of chemolithotrophy. In addition, nrf and psr gene clusters in T. oshmai JL-2 suggest respiratory nitrite ammonification and polysulfide reduction as possible modes of anaerobic respiration.

BatMeth: improved mapper for bisulfite sequencing reads on DNA methylation

DNA methylation plays a crucial role in higher organisms. Coupling bisulfite treatment with next generation sequencing enables the interrogation of 5-methylcytosine sites in the genome. However, bisulfite conversion introduces mismatches between the reads and the reference genome, which makes mapping of Illumina and SOLiD reads slow and inaccurate. BatMeth is an algorithm that integrates novel Mismatch Counting, List Filtering, Mismatch Stage Filtering and Fast Mapping onto Two Indexes components to improve unique mapping rate, speed and precision. Experimental results show that BatMeth is faster and more accurate than existing tools. BatMeth is freely available at http://code.google.com/p/batmeth/.

Long span DNA paired-end-tag (DNA-PET) sequencing strategy for the interrogation of genomic structural mutations and fusion-point-guided reconstruction of amplicons

Structural variations (SVs) contribute significantly to the variability of the human genome and extensive genomic rearrangements are a hallmark of cancer. While genomic DNA paired-end-tag (DNA-PET) sequencing is an attractive approach to identify genomic SVs, the current application of PET sequencing with short insert size DNA can be insufficient for the comprehensive mapping of SVs in low complexity and repeat-rich genomic regions. We employed a recently developed procedure to generate PET sequencing data using large DNA inserts of 10–20 kb and compared their characteristics with short insert (1 kb) libraries for their ability to identify SVs. Our results suggest that although short insert libraries bear an advantage in identifying small deletions, they do not provide significantly better breakpoint resolution. In contrast, large inserts are superior to short inserts in providing higher physical genome coverage for the same sequencing cost and achieve greater sensitivity, in practice, for the identification of several classes of SVs, such as copy number neutral and complex events. Furthermore, our results confirm that large insert libraries allow for the identification of SVs within repetitive sequences, which cannot be spanned by short inserts. This provides a key advantage in studying rearrangements in cancer, and we show how it can be used in a fusion-point-guided-concatenation algorithm to study focally amplified regions in cancer.

Target gene analysis by microarrays and chromatin immunoprecipitation identifies HEY proteins as highly redundant bHLH repressors

HEY bHLH transcription factors have been shown to regulate multiple key steps in cardiovascular development. They can be induced by activated NOTCH receptors, but other upstream stimuli mediated by TGFß and BMP receptors may elicit a similar response. While the basic and helix-loop-helix domains exhibit strong similarity, large parts of the proteins are still unique and may serve divergent functions. The striking overlap of cardiac defects in HEY2 and combined HEY1/HEYL knockout mice suggested that all three HEY genes fulfill overlapping function in target cells. We therefore sought to identify target genes for HEY proteins by microarray expression and ChIPseq analyses in HEK293 cells, cardiomyocytes, and murine hearts. HEY proteins were found to modulate expression of their target gene to a rather limited extent, but with striking functional interchangeability between HEY factors. Chromatin immunoprecipitation revealed a much greater number of potential binding sites that again largely overlap between HEY factors. Binding sites are clustered in the proximal promoter region especially of transcriptional regulators or developmental control genes. Multiple lines of evidence suggest that HEY proteins primarily act as direct transcriptional repressors, while gene activation seems to be due to secondary or indirect effects. Mutagenesis of putative DNA binding residues supports the notion of direct DNA binding. While class B E-box sequences (CACGYG) clearly represent preferred target sequences, there must be additional and more loosely defined modes of DNA binding since many of the target promoters that are efficiently bound by HEY proteins do not contain an E-box motif. These data clearly establish the three HEY bHLH factors as highly redundant transcriptional repressors in vitro and in vivo, which explains the combinatorial action observed in different tissues with overlapping expression.

NOTCH signaling is a central developmental pathway that influences a multitude of cell fate decisions and differentiation steps as well as later tissue homeostasis and regeneration. The three HEY genes encode basic helix-loop-helix transcription factors that are critical effectors to convey signaling by NOTCH receptors and similar signaling systems. This is underscored by the multitude of developmental defects observed in HEY single- and double-mutant mice. The mode of action of HEY proteins remained largely unexplored, however. By gene expression analysis and chromatin immunoprecipitation we have now identified a large set of HEY target genes. While only 500–2,000 mRNAs are regulated by HEY1 or HEY2, there are around 10,000 binding sites in the genome. HEY proteins act as transcriptional repressors that bind close to transcriptional start sites in all cases tested. In contrast, gene activation seems to be mediated by indirect/secondary mechanisms. The extent of regulation is rather limited, implicating HEY genes in modulating rather than switching on or off target gene expression. All our in vitro and in vivo data point to a high degree of redundancy between the three HEY genes, suggesting that tissue specific patterns and expression levels determine the final outcome of HEY induced cellular responses.