Affymetrix SNP arrays have been widely used for single-nucleotide polymorphism (SNP) genotype calling and DNA copy number variation inference. Although numerous methods have achieved high accuracy in these fields, most studies have paid little attention to the modeling of hybridization of probes to off-target allele sequences, which can affect the accuracy greatly. In this study, we address this issue and demonstrate that hybridization with mismatch nucleotides (HWMMN) occurs in all SNP probe-sets and has a critical effect on the estimation of allelic concentrations (ACs). We study sequence binding through binding free energy and then binding affinity, and develop a probe intensity composite representation (PICR) model. The PICR model allows the estimation of ACs at a given SNP through statistical regression. Furthermore, we demonstrate with cell-line data of known true copy numbers that the PICR model can achieve reasonable accuracy in copy number estimation at a single SNP locus, by using the ratio of the estimated AC of each sample to that of the reference sample, and can reveal subtle genotype structure of SNPs at abnormal loci. We also demonstrate with HapMap data that the PICR model yields accurate SNP genotype calls consistently across samples, laboratories and even across array platforms.

Chromatin remodelers are essential for establishing and maintaining the placement of nucleosomes along genomic DNA. Yet how chromatin remodelers recognize and respond to distinct chromatin environments surrounding nucleosomes is poorly understood. Here, we use Lac repressor as a tool to probe how a DNA-bound factor influences action of the Chd1 remodeler. We show that Chd1 preferentially shifts nucleosomes away from Lac repressor, demonstrating that a DNA-bound factor defines a barrier for nucleosome positioning. Rather than an absolute block in sliding, the barrier effect was achieved by altered rates of nucleosome sliding that biased redistribution of nucleosomes away from the bound Lac repressor site. Remarkably, in addition to slower sliding toward the LacO site, the presence of Lac repressor also stimulated sliding in the opposite direction. These experiments therefore demonstrate that Chd1 responds to the presence of a bound protein on both entry and exit sides of the nucleosome. This sensitivity to both sides of the nucleosome allows for a faster and sharper response than would be possible by responding to only the entry side, and we speculate that dual entry/exit sensitivity is also important for regularly spaced nucleosome arrays generated by Chd1 and the related ISWI remodelers.

Prokaryotes memorize invader information by incorporating alien DNA as spacers into CRISPR arrays. Although the spacer size has been suggested to be predefined by the architecture of the acquisition complex, there is usually an unexpected heterogeneity. Here, we explored the causes of this heterogeneity in Haloarcula hispanica I-B CRISPR. High-throughput sequencing following adaptation assays demonstrated significant size variation among 37 957 new spacers, which appeared to be sequence-dependent. Consistently, the third nucleotide at the spacer 3΄-end (PAM-distal end) showed an evident bias for cytosine and mutating this cytosine in the protospacer sequence could change the final spacer size. In addition, slippage of the 5΄-end (PAM-end), which contributed to most of the observed PAM (protospacer adjacent motif) inaccuracy, also tended to change the spacer size. We propose that both ends of the PAM-protospacer sequence should exhibit nucleotide selectivity (with different stringencies), which fine-tunes the structural ruler, to a certain extent, to specify the spacer size.

The hexametric T7 helicase (gp4) adopts a spiral lock-washer form and encircles a coil-like DNA (tracking) strand with two nucleotides bound to each subunit. However, the chemo-mechanical coupling mechanism in unwinding has yet to be elucidated. Here, we utilized nanotensioner-enhanced Förster resonance energy transfer with one nucleotide precision to investigate gp4-induced unwinding of DNA that contains an abasic lesion. We observed that the DNA unwinding activity of gp4 is hindered but not completely blocked by abasic lesions. Gp4 moves back and forth repeatedly when it encounters an abasic lesion, whereas it steps back only occasionally when it unwinds normal DNA. We further observed that gp4 translocates on the tracking strand in step sizes of one to four nucleotides. We propose that a hypothetical intermediate conformation of the gp4-DNA complex during DNA unwinding can help explain how gp4 molecules pass lesions, providing insights into the unwinding dynamics of gp4.

The advent of single-cell sequencing opened a new era in transcriptomic and genomic research. To understand cell composition using single-cell studies, a variety of cell markers have been widely used to label individual cell types. However, the specific database of cell markers for use by the plant research community remains very limited. To overcome this problem, we developed the Plant Cell Marker DataBase (PCMDB, http://www.tobaccodb.org/pcmdb/), which is based on a uniform annotation pipeline. By manually curating over 130 000 research publications, we collected a total of 81 117 cell marker genes of 263 cell types in 22 tissues across six plant species. Tissue- and cell-specific expression patterns can be visualized using multiple tools: eFP Browser, Bar, and UMAP/TSNE graph. The PCMDB also supports several analysis tools, including SCSA and SingleR, which allows for user annotation of cell types. To provide information about plant species currently unsupported in PCMDB, potential marker genes for other plant species can be searched based on homology with the supported species. PCMDB is a user-friendly hierarchical platform that contains five built-in search engines. We believe PCMDB will constitute a useful resource for researchers working on cell type annotation and the prediction of the biological function of individual cells.

Glycogen synthase kinase 3 beta (GSK3β) is a critical protein kinase that phosphorylates numerous proteins in cells and thereby impacts multiple pathways including the β-Catenin/TCF/LEF-1 pathway. MicroRNAs (miRs) are a class of noncoding small RNAs of ∼22 nucleotides in length. Both GSK3β and miR play myriad roles in cell functions including stem cell development, apoptosis, embryogenesis and tumorigenesis. Here we show that GSK3β inhibits the expression of miR-96, miR-182 and miR-183 through the β-Catenin/TCF/LEF-1 pathway. Knockout of GSK3β in mouse embryonic fibroblast cells increases expression of miR-96, miR-182 and miR-183, coinciding with increases in the protein level and nuclear translocation of β-Catenin. In addition, overexpression of β-Catenin enhances the expression of miR-96, miR-182 and miR-183 in human gastric cancer AGS cells. GSK3β protein levels are decreased in human gastric cancer tissue compared with surrounding normal gastric tissue, coinciding with increases of β-Catenin protein, miR-96, miR-182, miR-183 and primary miR-183-96-182 cluster (pri-miR-183). Furthermore, suppression of miR-183-96-182 cluster with miRCURY LNA miR inhibitors decreases the proliferation and migration of AGS cells. Knockdown of GSK3β with siRNA increases the proliferation of AGS cells. Mechanistically, we show that β-Catenin/TCF/LEF-1 binds to the promoter of miR-183-96-182 cluster gene and thereby activates the transcription of the cluster. In summary, our findings identify a novel role for GSK3β in the regulation of miR-183-96-182 biogenesis through β-Catenin/TCF/LEF-1 pathway in gastric cancer cells.

To determine early somatic changes in high-grade serous ovarian cancer (HGSOC), we performed whole genome sequencing on a rare collection of 16 low stage HGSOCs. The majority showed extensive structural alterations (one had an ultramutated profile), exhibited high levels of p53 immunoreactivity, and harboured a TP53 mutation, deletion or inactivation. BRCA1 and BRCA2 mutations were observed in two tumors, with nine showing evidence of a homologous recombination (HR) defect. Combined Analysis with The Cancer Genome Atlas (TCGA) indicated that low and late stage HGSOCs have similar mutation and copy number profiles. We also found evidence that deleterious TP53 mutations are the earliest events, followed by deletions or loss of heterozygosity (LOH) of chromosomes carrying TP53, BRCA1 or BRCA2. Inactivation of HR appears to be an early event, as 62.5% of tumours showed a LOH pattern suggestive of HR defects. Three tumours with the highest ploidy had little genome-wide LOH, yet one of these had a homozygous somatic frame-shift BRCA2 mutation, suggesting that some carcinomas begin as tetraploid then descend into diploidy accompanied by genome-wide LOH. Lastly, we found evidence that structural variants (SV) cluster in HGSOC, but are absent in one ultramutated tumor, providing insights into the pathogenesis of low stage HGSOC.

Alternative DNA structures that deviate from B-form double-stranded DNA such as G-quadruplex (G4) DNA can be formed by G-rich sequences that are widely distributed throughout the human genome. We have previously shown that Pif1p not only unfolds G4, but also unwinds the downstream duplex DNA in a G4-stimulated manner. In the present study, we further characterized the G4-stimulated duplex DNA unwinding phenomenon by means of single-molecule fluorescence resonance energy transfer. It was found that Pif1p did not unwind the partial duplex DNA immediately after unfolding the upstream G4 structure, but rather, it would dwell at the ss/dsDNA junction with a 'waiting time'. Further studies revealed that the waiting time was in fact related to a protein dimerization process that was sensitive to ssDNA sequence and would become rapid if the sequence is G-rich. Furthermore, we identified that the G-rich sequence, as the G4 structure, equally stimulates duplex DNA unwinding. The present work sheds new light on the molecular mechanism by which G4-unwinding helicase Pif1p resolves physiological G4/duplex DNA structures in cells.

The histone chaperone FACT (FAcilitates Chromatin Transcription) plays an essential role in transcription and DNA replication by its dual functions on nucleosome assembly to maintain chromatin integrity and nucleosome disassembly to destabilize nucleosome and facilitate its accessibility simultaneously. Mono-ubiquitination at Lysine 119 of H2A (ubH2A) has been suggested to repress transcription by preventing the recruitment of FACT at early elongation process. However, up to date, how ubH2A directly affects FACT on nucleosome assembly and disassembly remains elusive. In this study, we demonstrated that the dual functions of FACT are differently regulated by ubH2A. The H2A ubiquitination does not affect FACT's chaperone function in nucleosome assembly and FACT can deposit ubH2A-H2B dimer on tetrasome to form intact nucleosome. However, ubH2A greatly restricts FACT binding on nucleosome and inhibits its activity of nucleosome disassembly. Interestingly, deubiquitination of ubH2A rescues the nucleosome disassembly function of FACT to activate gene transcription. Our findings provide mechanistic insights of how H2A ubiquitination affects FACT in breaking nucleosome and maintaining its integrity, which sheds light on the biological function of ubH2A and various FACT's activity under different chromatin states.

Recent studies on CRISPR adaptation revealed that priming is a major pathway of spacer acquisition, at least for the most prevalent type I systems. Priming is guided by a CRISPR RNA which fully/partially matches the invader DNA, but the plasticity of this RNA guide has not yet been characterized. In this study, we extensively modified the two conserved handles of a priming crRNA in Haloarcula hispanica, and altered the size of its central spacer part. Interestingly, priming is insusceptible to the full deletion of 3' handle, which seriously impaired crRNA stability and interference effects. With 3' handle deletion, further truncation of 5' handle revealed that its spacer-proximal 6 nucleotides could provide the least conserved sequence required for priming. Subsequent scanning mutation further identified critical nucleotides within 5' handle. Besides, priming was also shown to tolerate a wider size variation of the spacer part, compared to interference. These data collectively illustrate the high tolerance of priming to extensive structural/size variations of the crRNA guide, which highlights the structural flexibility of the crRNA-effector ribonucleoprotein complex. The observed high priming effectiveness suggests that primed adaptation promotes clearance of the fast-replicating and ever-evolving viral DNA, by rapidly and persistently multiplexing the interference pathway.

Recent discovery of ectopic repeats (outside CRISPR arrays) provided unprecedented insights into the nondefense roles of CRISPR-Cas. A striking example is the addiction module CreTA (CRISPR-regulated toxin-antitoxins), where one or two (in most cases) ectopic repeats produce CRISPR-resembling antitoxic (CreA) RNAs that direct the CRISPR effector Cascade to transcriptionally repress a toxic RNA (CreT). Here, we demonstrated that CreTA repeats are extensively degenerated in sequence, with the first repeat (ψR1) being more diverged than the second one (ψR2). As a result, such addiction modules become highly specific to their physically-linked CRISPR-Cas loci, and in most cases, CreA could not harness a heterologous CRISPR-Cas to suppress its cognate toxin. We further disclosed that this specificity primarily derives from the degeneration of ψR1, and could generally be altered by modifying this repeat element. We also showed that the degenerated repeats of CreTA were insusceptible to recombination and thus more stable compared to a typical CRISPR array, which could be exploited to develop highly stable CRISPR-based tools. These data illustrated that repeat degeneration (a common feature of ectopic repeats) improves the stability and specificity of CreTA in protecting CRISPR-Cas, which could have contributed to the widespread occurrence and deep diversification of CRISPR systems.

Pif1 is a prototypical member of the 5' to 3' DNA helicase family conserved from bacteria to human. It has a high binding affinity for DNA, but unwinds double-stranded DNA (dsDNA) with a low processivity. Efficient DNA unwinding has been observed only at high protein concentrations that favor dimerization of Pif1. In this research, we used single-molecule fluorescence resonance energy transfer (smFRET) and magnetic tweezers (MT) to study the DNA unwinding activity of Saccharomyces cerevisiae Pif1 (Pif1) under different forces exerted on the tails of a forked dsDNA. We found that Pif1 can unwind the forked DNA repetitively for many unwinding-rezipping cycles at zero force. However, Pif1 was found to have a very limited processivity in each cycle because it loosened its strong association with the tracking strand readily, which explains why Pif1 cannot be observed to unwind DNA efficiently in bulk assays at low protein concentrations. The force enhanced the unwinding rate and the total unwinding length of Pif1 significantly. With a force of 9 pN, the rate and length were enhanced by more than 3- and 20-fold, respectively. Our results imply that the DNA unwinding activity of Pif1 can be regulated by force. The relevance of this characteristic of Pif1 to its cellular functions is discussed.

Mutations in the RecQ DNA helicase gene BLM give rise to Bloom's syndrome, which is a rare autosomal recessive disorder characterized by genetic instability and cancer predisposition. BLM helicase is highly active in binding and unwinding G-quadruplexes (G4s), which are physiological targets for BLM, as revealed by genome-wide characterizations of gene expression of cells from BS patients. With smFRET assays, we studied the molecular mechanism of BLM-catalyzed G4 unfolding and showed that ATP is required for G4 unfolding. Surprisingly, depending on the molecular environments of G4, BLM unfolds G4 through different mechanisms: unfolding G4 harboring a 3'-ssDNA tail in three discrete steps with unidirectional translocation, and unfolding G4 connected to dsDNA by ssDNA in a repetitive manner in which BLM remains anchored at the ss/dsDNA junction, and G4 was unfolded by reeling in ssDNA. This indicates that one BLM molecule may unfold G4s in different molecular environments through different mechanisms.

Pif1 helicases are ubiquitous members of the SF1B family and are essential for maintaining genome stability. It was speculated that Pif1-specific motifs may fold in specific structures, conferring distinct activities upon it. Here, we report the crystal structures of the Pif1 helicase from Bacteroides spp with and without adenosine triphosphate (ATP) analog/ssDNA. BsPif1 shares structural similarities with RecD2 and Dda helicases but has specific features in the 1B and 2B domains. The highly conserved Pif1 family specific sequence motif interacts with and constraints a putative pin-loop in domain 1B in a precise conformation. More importantly, we found that the 2B domain which contains a specific extended hairpin undergoes a significant rotation and/or movement upon ATP and DNA binding, which is absolutely required for DNA unwinding. We therefore propose a mechanism for DNA unwinding in which the 2B domain plays a predominant role. The fact that the conformational change regulates Pif1 activity may provide insight into the puzzling observation that Pif1 becomes highly processive during break-induced replication in association with Polδ, while the isolated Pif1 has low processivity.

We have developed a novel method for estimating the parameters of hidden Markov models for gene finding in newly sequenced species. Our approach does not rely on curated training data sets, but instead uses extrinsic evidence (including paired-end ditags that have not been used in gene finding previously) and iterative training. This new method is particularly suitable for annotation of species with large evolutionary distance to the closest annotated species. We have used our approach to produce an initial annotation of more than 16,000 genes in the newly sequenced Schistosoma japonicum draft genome. We established the high quality of our predictions by comparison to full-length cDNAs (withdrawn from the extrinsic evidence) and to CEGMA core genes. We also evaluated the effectiveness of the new training procedure on Caenorhabditis elegans genome. ExonHunter and the newest parametric files for S. japonicum genome are available for download at www.bioinformatics.uwaterloo.ca/downloads/exonhunter.

Archaeal ribosomes have many domain-specific features; however, our understanding of these structures is limited. We present 10 cryo-electron microscopy (cryo-EM) structures of the archaeal ribosome from crenarchaeota Sulfolobus acidocaldarius (Sac) at 2.7-5.7 Å resolution. We observed unstable conformations of H68 and h44 of ribosomal RNA (rRNA) in the subunit structures, which may interfere with subunit association. These subunit structures provided models for 12 rRNA expansion segments and 3 novel r-proteins. Furthermore, the 50S-aRF1 complex structure showed the unique domain orientation of aRF1, possibly explaining P-site transfer RNA (tRNA) release after translation termination. Sac 70S complexes were captured in seven distinct steps of the tRNA translocation reaction, confirming conserved structural features during archaeal ribosome translocation. In aEF2-engaged 70S ribosome complexes, 3D classification of cryo-EM data based on 30S head domain identified two new translocation intermediates with 30S head domain tilted 5-6° enabling its disengagement from the translocated tRNA and its release post-translocation. Additionally, we observed conformational changes to aEF2 during ribosome binding and switching from three different states. Our structural and biochemical data provide new insights into archaeal translation and ribosome translocation.

Understanding the role of a given transcription factor (TF) in regulating gene expression requires precise mapping of its binding sites in the genome. Chromatin immunoprecipitation-exo, an emerging technique using λ exonuclease to digest TF unbound DNA after ChIP, is designed to reveal transcription factor binding site (TFBS) boundaries with near-single nucleotide resolution. Although ChIP-exo promises deeper insights into transcription regulation, no dedicated bioinformatics tool exists to leverage its advantages. Most ChIP-seq and ChIP-chip analytic methods are not tailored for ChIP-exo, and thus cannot take full advantage of high-resolution ChIP-exo data. Here we describe a novel analysis framework, termed MACE (model-based analysis of ChIP-exo) dedicated to ChIP-exo data analysis. The MACE workflow consists of four steps: (i) sequencing data normalization and bias correction; (ii) signal consolidation and noise reduction; (iii) single-nucleotide resolution border peak detection using the Chebyshev Inequality and (iv) border matching using the Gale-Shapley stable matching algorithm. When applied to published human CTCF, yeast Reb1 and our own mouse ONECUT1/HNF6 ChIP-exo data, MACE is able to define TFBSs with high sensitivity, specificity and spatial resolution, as evidenced by multiple criteria including motif enrichment, sequence conservation, direct sequence pileup, nucleosome positioning and open chromatin states. In addition, we show that the fundamental advance of MACE is the identification of two boundaries of a TFBS with high resolution, whereas other methods only report a single location of the same event. The two boundaries help elucidate the in vivo binding structure of a given TF, e.g. whether the TF may bind as dimers or in a complex with other co-factors.

Dr.VIS is a database of human disease-related viral integration sites (VIS). The number of VIS has grown rapidly since Dr.VIS was first released in 2011, and there is growing recognition of the important role that viral integration plays in the development of malignancies. The updated database version, Dr.VIS v2.0 (http://www.bioinfo.org/drvis or bminfor.tongji.edu.cn/drvis_v2), represents 25 diseases, covers 3340 integration sites of eight oncogenic viruses in human chromosomes and provides more accurate information about VIS from high-throughput deep sequencing results obtained mainly after 2012. Data of VISes for three newly identified oncogenic viruses for 14 related diseases have been added to this 2015 update, which has a 5-fold increase of VISes compared to Dr.VIS v1.0. Dr.VIS v2.0 has 2244 precise integration sites, 867 integration regions and 551 junction sequences. A total of 2295 integration sites are located near 1730 involved genes. Of the VISes, 1153 are detected in the exons or introns of genes, with 294 located up to 5 kb and a further 112 located up to 10 kb away. As viral integration may alter chromosome stability and gene expression levels, characterizing VISes will contribute toward the discovery of novel oncogenes, tumor suppressor genes and tumor-associated pathways.

The clustered regularly interspaced short palindromic repeat (CRISPR)-Cas system mediates adaptive immunity against foreign nucleic acids in prokaryotes. However, efficient adaptation of a native CRISPR to purified viruses has only been observed for the type II-A system from a Streptococcus thermophilus industry strain, and rarely reported for laboratory strains. Here, we provide a second native system showing efficient adaptation. Infected by a newly isolated virus HHPV-2, Haloarcula hispanica type I-B CRISPR system acquired spacers discriminatively from viral sequences. Unexpectedly, in addition to Cas1, Cas2 and Cas4, this process also requires Cas3 and at least partial Cascade proteins, which are involved in interference and/or CRISPR RNA maturation. Intriguingly, a preexisting spacer partially matching a viral sequence is also required, and spacer acquisition from upstream and downstream sequences of its target sequence (i.e. priming protospacer) shows different strand bias. These evidences strongly indicate that adaptation in this system strictly requires a priming process. This requirement, if validated also true for other CRISPR systems as implied by our bioinformatic analysis, may help to explain failures to observe efficient adaptation to purified viruses in many laboratory strains, and the discrimination mechanism at the adaptation level that has confused scientists for years.

Clustered Regularly Interspaced Short Palindromic Repeats (CRISPRs) acquire new spacers to generate adaptive immunity in prokaryotes. During spacer integration, the leader-preceded repeat is always accurately duplicated, leading to speculations of a repeat-length ruler. Here in Haloarcula hispanica, we demonstrate that the accurate duplication of its 30-bp repeat requires two conserved mid-repeat motifs, AACCC and GTGGG. The AACCC motif was essential and needed to be ∼10 bp downstream from the leader-repeat junction site, where duplication consistently started. Interestingly, repeat duplication terminated sequence-independently and usually with a specific distance from the GTGGG motif, which seemingly served as an anchor site for a molecular ruler. Accordingly, altering the spacing between the two motifs led to an aberrant duplication size (29, 31, 32 or 33 bp). We propose the adaptation complex may recognize these mid-repeat elements to enable measuring the repeat DNA for spacer integration.