De novo assembly of sequence data has more specific contigs than unassembled sequence data.

Effects of GC Bias in Next-Generation-Sequencing Data on De Novo Genome Assembly

Next-generation-sequencing (NGS) has revolutionized the field of genome assembly because of its much higher data throughput and much lower cost compared with traditional Sanger sequencing. However, NGS poses new computational challenges to de novo genome assembly. Among the challenges, GC bias in NGS data is known to aggravate genome assembly. However, it is not clear to what extent GC bias affects genome assembly in general. In this work, we conduct a systematic analysis on the effects of GC bias on genome assembly. Our analyses reveal that GC bias only lowers assembly completeness when the degree of GC bias is above a threshold. At a strong GC bias, the assembly fragmentation due to GC bias can be explained by the low coverage of reads in the GC-poor or GC-rich regions of a genome. This effect is observed for all the assemblers under study. Increasing the total amount of NGS data thus rescues the assembly fragmentation because of GC bias. However, the amount of data needed for a full rescue depends on the distribution of GC contents. Both low and high coverage depths due to GC bias lower the accuracy of assembly. These pieces of information provide guidance toward a better de novo genome assembly in the presence of GC bias.

BLAST+: architecture and applications

BACKGROUND Sequence similarity searching is a very important bioinformatics task. While Basic Local Alignment Search Tool (BLAST) outperforms exact methods through its use of heuristics, the speed of the current BLAST software is suboptimal for very long queries or database sequences. There are also some shortcomings in the user-interface of the current command-line applications. RESULTS We describe features and improvements of rewritten BLAST software and introduce new command-line applications. Long query sequences are broken into chunks for processing, in some cases leading to dramatically shorter run times. For long database sequences, it is possible to retrieve only the relevant parts of the sequence, reducing CPU time and memory usage for searches of short queries against databases of contigs or chromosomes. The program can now retrieve masking information for database sequences from the BLAST databases. A new modular software library can now access subject sequence data from arbitrary data sources. We introduce several new features, including strategy files that allow a user to save and reuse their favorite set of options. The strategy files can be uploaded to and downloaded from the NCBI BLAST web site. CONCLUSION The new BLAST command-line applications, compared to the current BLAST tools, demonstrate substantial speed improvements for long queries as well as chromosome length database sequences. We have also improved the user interface of the command-line applications.

RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome

BACKGROUND RNA-Seq is revolutionizing the way transcript abundances are measured. A key challenge in transcript quantification from RNA-Seq data is the handling of reads that map to multiple genes or isoforms. This issue is particularly important for quantification with de novo transcriptome assemblies in the absence of sequenced genomes, as it is difficult to determine which transcripts are isoforms of the same gene. A second significant issue is the design of RNA-Seq experiments, in terms of the number of reads, read length, and whether reads come from one or both ends of cDNA fragments. RESULTS We present RSEM, an user-friendly software package for quantifying gene and isoform abundances from single-end or paired-end RNA-Seq data. RSEM outputs abundance estimates, 95% credibility intervals, and visualization files and can also simulate RNA-Seq data. In contrast to other existing tools, the software does not require a reference genome. Thus, in combination with a de novo transcriptome assembler, RSEM enables accurate transcript quantification for species without sequenced genomes. On simulated and real data sets, RSEM has superior or comparable performance to quantification methods that rely on a reference genome. Taking advantage of RSEM's ability to effectively use ambiguously-mapping reads, we show that accurate gene-level abundance estimates are best obtained with large numbers of short single-end reads. On the other hand, estimates of the relative frequencies of isoforms within single genes may be improved through the use of paired-end reads, depending on the number of possible splice forms for each gene. CONCLUSIONS RSEM is an accurate and user-friendly software tool for quantifying transcript abundances from RNA-Seq data. As it does not rely on the existence of a reference genome, it is particularly useful for quantification with de novo transcriptome assemblies. In addition, RSEM has enabled valuable guidance for cost-efficient design of quantification experiments with RNA-Seq, which is currently relatively expensive.

Full-length transcriptome assembly from RNA-Seq data without a reference genome.

Massively parallel sequencing of cDNA has enabled deep and efficient probing of transcriptomes. Current approaches for transcript reconstruction from such data often rely on aligning reads to a reference genome, and are thus unsuitable for samples with a partial or missing reference genome. Here we present the Trinity method for de novo assembly of full-length transcripts and evaluate it on samples from fission yeast, mouse and whitefly, whose reference genome is not yet available. By efficiently constructing and analyzing sets of de Bruijn graphs, Trinity fully reconstructs a large fraction of transcripts, including alternatively spliced isoforms and transcripts from recently duplicated genes. Compared with other de novo transcriptome assemblers, Trinity recovers more full-length transcripts across a broad range of expression levels, with a sensitivity similar to methods that rely on genome alignments. Our approach provides a unified solution for transcriptome reconstruction in any sample, especially in the absence of a reference genome.

IgBLAST: an immunoglobulin variable domain sequence analysis tool

The variable domain of an immunoglobulin (IG) sequence is encoded by multiple genes, including the variable (V) gene, the diversity (D) gene and the joining (J) gene. Analysis of IG sequences typically requires identification of each gene, as well as a comparison of sequence variations in the context of defined regions. General purpose tools, such as the BLAST program, have only limited use for such tasks, as the rearranged nature of an IG sequence and the variable length of each gene requires multiple rounds of BLAST searches for a single IG sequence. Additionally, manual assembly of different genes is difficult and error-prone. To address these issues and to facilitate other common tasks in analysing IG sequences, we have developed the sequence analysis tool IgBLAST (http://www.ncbi.nlm.nih.gov/igblast/). With this tool, users can view the matches to the germline V, D and J genes, details at rearrangement junctions, the delineation of IG V domain framework regions and complementarity determining regions. IgBLAST has the capability to analyse nucleotide and protein sequences and can process sequences in batches. Furthermore, IgBLAST allows searches against the germline gene databases and other sequence databases simultaneously to minimize the chance of missing possibly the best matching germline V gene.

De novo assembly of a PML nuclear subcompartment occurs through multiple pathways and induces telomere elongation.

Telomerase-negative tumor cells use an alternative lengthening of telomeres (ALT) pathway that involves DNA recombination and repair to maintain their proliferative potential. The cytological hallmark of this process is the accumulation of promyelocytic leukemia (PML) nuclear protein at telomeric DNA to form ALT-associated PML bodies (APBs). Here, the de novo formation of a telomeric PML nuclear subcompartment was investigated by recruiting APB protein components. We show that functionally distinct proteins were able to initiate the formation of bona fide APBs with high efficiency in a self-organizing and self-propagating manner. These included: (1) PML and Sp100 as the constituting components of PML nuclear bodies, (2) telomere repeat binding factors 1 and 2 (TRF1 and TRF2, respectively), (3) the DNA repair protein NBS1 and (4) the SUMO E3 ligase MMS21, as well as the isolated SUMO1 domain, through an interacting domain of another protein factor. By contrast, the repair factors Rad9, Rad17 and Rad51 were less efficient in APB nucleation but were recruited to preassembled APBs. The artificially created APBs induced telomeric extension through a DNA repair mechanism, as inferred from their colocalization with sites of non-replicative DNA synthesis and histone H2A.X phosphorylation, and an increase of the telomere repeat length. These activities were absent after recruitment of the APB factors to a pericentric locus and establish APBs as functional intermediates of the ALT pathway.

position-specific gap penalties

The sensitivity of the commonly used progressive multiple sequence alignment method has been greatly improved for the alignment of divergent protein sequences. Firstly, individual weights are assigned to each sequence in a partial alignment in order to down-weight near-duplicate sequences and up-weight the most divergent ones. Secondly, amino acid substitution matrices are varied at different alignment stages according to the divergence of the sequences to be aligned. Thirdly, residue-specific gap penalties and locally reduced gap penalties in hydrophilic regions encourage new gaps in potential loop regions rather than regular secondary structure. Fourthly, positions in early alignments where gaps have been opened receive locally reduced gap penalties to encourage the opening up of new gaps at these positions. These modifications are incorporated into a new program, CLUSTAL W which is freely available.

Human artificial chromosomes with alpha satellite-based de novo centromeres show increased frequency of nondisjunction and anaphase lag.

Human artificial chromosomes have been used to model requirements for human chromosome segregation and to explore the nature of sequences competent for centromere function. Normal human centromeres require specialized chromatin that consists of alpha satellite DNA complexed with epigenetically modified histones and centromere-specific proteins. While several types of alpha satellite DNA have been used to assemble de novo centromeres in artificial chromosome assays, the extent to which they fully recapitulate normal centromere function has not been explored. Here, we have used two kinds of alpha satellite DNA, DXZ1 (from the X chromosome) and D17Z1 (from chromosome 17), to generate human artificial chromosomes. Although artificial chromosomes are mitotically stable over many months in culture, when we examined their segregation in individual cell divisions using an anaphase assay, artificial chromosomes exhibited more segregation errors than natural human chromosomes (P < 0.001). Naturally occurring, but abnormal small ring chromosomes derived from chromosome 17 and the X chromosome also missegregate more than normal chromosomes, implicating overall chromosome size and/or structure in the fidelity of chromosome segregation. As different artificial chromosomes missegregate over a fivefold range, the data suggest that variable centromeric DNA content and/or epigenetic assembly can influence the mitotic behavior of artificial chromosomes.

Strand-specific mismatch correction in nuclear extracts of human and Drosophila melanogaster cell lines.

Nuclear extracts derived from HeLa and Drosophila melanogaster KC cell lines have been found to correct single base-base mispairs within open circular DNA heteroduplexes containing a strand-specific, site-specific incision located 808 base pairs from the mismatch. Correction in both extract systems is strand specific, being highly biased to the incised DNA strand. Different mispairs within a homologous set of heteroduplexes were processed with different efficiencies (G.T greater than G.G approximately equal to A.C greater than C.C), and correction was accompanied by mismatch-dependent DNA synthesis localized to the region spanning the mispair and the strand break, thus demonstrating that mismatch recognition is associated with the repair reaction. Correction of each of these heteroduplexes was abolished by aphidicolin but was relatively insensitive to the presence of high concentrations of ddTTP, indicating probable involvement of alpha and/or delta class DNA polymerase(s). These findings suggest that higher eukaryotic cells possess a general, strand-specific mismatch repair system analogous to the Escherichia coli mutHLS and the Streptococcus pneumoniae hexAB pathways, systems that contribute in a major way to the genetic stability of these bacterial species.

Monitoring spatiotemporal biogenesis of macromolecular assemblies by pulse-chase epitope labeling.

Many cellular proteins perform their roles within macromolecular assemblies. Hence, an understanding of how these multiprotein complexes form is a fundamental question in cell biology. We developed a translation-controlled pulse-chase system that allows time-resolved isolation of newly forming multiprotein complexes in chemical quantities suitable for biochemical and cell biological analysis. The "pulse" is triggered by an unnatural amino acid, which induces immediate translation of an amber stop codon repressed mRNA encoding the protein of interest with a built-in tag for detection and purification. The "chase" is elicited by stopping translation of this bait via a riboswitch in the respective mRNA. Over the course of validating our method, we discovered a distinct time-resolved assembly step during NPC biogenesis and could directly monitor the spatiotemporal maturation of preribosomes via immunofluorescence detection and purification of a pulse-labeled ribosomal protein. Thus, we provide an innovative strategy to study dynamic protein assembly within cellular networks.

DNA methylation status predicts cell type-specific enhancer activity.

Cell-selective glucocorticoid receptor (GR) binding to distal regulatory elements is associated with cell type-specific regions of locally accessible chromatin. These regions can either pre-exist in chromatin (pre-programmed) or be induced by the receptor (de novo). Mechanisms that create and maintain these sites are not well understood. We observe a global enrichment of CpG density for pre-programmed elements, and implicate their demethylated state in the maintenance of open chromatin in a tissue-specific manner. In contrast, sites that are actively opened by GR (de novo) are characterized by low CpG density, and form a unique class of enhancers devoid of suppressive effect of agglomerated methyl-cytosines. Furthermore, treatment with glucocorticoids induces rapid changes in methylation levels at selected CpGs within de novo sites. Finally, we identify GR-binding elements with CpGs at critical positions, and show that methylation can affect GR-DNA interactions in vitro. The findings present a unique link between tissue-specific chromatin accessibility, DNA methylation and transcription factor binding and show that DNA methylation can be an integral component of gene regulation by nuclear receptors.

N- and E-cadherins in Xenopus are specifically required in the neural and non-neural ectoderm, respectively, for F-actin assembly and morphogenetic movements.

Transmembrane cadherins are calcium-dependent intercellular adhesion molecules. Recently, they have also been shown to be sites of actin assembly during adhesive contact formation. However, the roles of actin assembly on transmembrane cadherins during development are not fully understood. We show here, using the developing ectoderm of the Xenopus embryo as a model, that F-actin assembly is a primary function of both N-cadherin in the neural ectoderm and E-cadherin in the non-neural (epidermal) ectoderm, and that each cadherin is essential for the characteristic morphogenetic movements of these two tissues. However, depletion of N-cadherin and E-cadherin did not cause dissociation in these tissues at the neurula stage, probably owing to the expression of C-cadherin in each tissue. Depletion of each of these cadherins is not rescued by the other, nor by the expression of C-cadherin, which is expressed in both tissues. One possible reason for this is that each cadherin is expressed in a different domain of the cell membrane. These data indicate the combinatorial nature of cadherin function, the fact that N- and E-cadherin play primary roles in F-actin assembly in addition to roles in cell adhesion, and that this function is specific to individual cadherins. They also show how cell adhesion and motility can be combined in morphogenetic tissue movements that generate the form and shape of the embryonic organs.

Nucleation of nuclear bodies by RNA

The biogenesis of the many functional compartments contained in the mammalian cell nucleus is poorly understood. More specifically, little is known regarding the initial nucleation step required for nuclear body formation. Here we show that RNA can function as a structural element and a nucleator of nuclear bodies. We find that several types of coding and noncoding RNAs are sufficient to de novo assemble, and are physiologically enriched in, histone locus bodies (with associated Cajal bodies), nuclear speckles, paraspeckles and nuclear stress bodies. Formation of nuclear bodies occurs through recruitment and accumulation of proteins resident in the nuclear bodies by nucleating RNA. These results demonstrate that transcription is a driving force in nuclear body formation and RNA transcripts can function as a scaffold in the formation of major nuclear bodies. Together, these data suggest that RNA-primed biogenesis of nuclear bodies is a general principle of nuclear organization.

The miniature genome of a carnivorous plant Genlisea aurea contains a low number of genes and short non-coding sequences

BACKGROUND Genlisea aurea (Lentibulariaceae) is a carnivorous plant with unusually small genome size - 63.6 Mb - one of the smallest known among higher plants. Data on the genome sizes and the phylogeny of Genlisea suggest that this is a derived state within the genus. Thus, G. aurea is an excellent model organism for studying evolutionary mechanisms of genome contraction. RESULTS Here we report sequencing and de novo draft assembly of G. aurea genome. The assembly consists of 10,687 contigs of the total length of 43.4 Mb and includes 17,755 complete and partial protein-coding genes. Its comparison with the genome of Mimulus guttatus, another representative of higher core Lamiales clade, reveals striking differences in gene content and length of non-coding regions. CONCLUSIONS Genome contraction was a complex process, which involved gene loss and reduction of lengths of introns and intergenic regions, but not intron loss. The gene loss is more frequent for the genes that belong to multigenic families indicating that genetic redundancy is an important prerequisite for genome size reduction.

CENP-B Controls Centromere Formation Depending on the Chromatin Context

The centromere is a chromatin region that serves as the spindle attachment point and directs accurate inheritance of eukaryotic chromosomes during cell divisions. However, the mechanism by which the centromere assembles and stabilizes at a specific genomic region is not clear. The de novo formation of a human/mammalian artificial chromosome (HAC/MAC) with a functional centromere assembly requires the presence of alpha-satellite DNA containing binding motifs for the centromeric CENP-B protein. We demonstrate here that de novo centromere assembly on HAC/MAC is dependent on CENP-B. In contrast, centromere formation is suppressed in cells expressing CENP-B when alpha-satellite DNA was integrated into a chromosomal site. Remarkably, on those integration sites CENP-B enhances histone H3-K9 trimethylation and DNA methylation, thereby stimulating heterochromatin formation. Thus, we propose that CENP-B plays a dual role in centromere formation, ensuring de novo formation on DNA lacking a functional centromere but preventing the formation of excess centromeres on chromosomes.

Improved tools for biological sequence comparison.

We have developed three computer programs for comparisons of protein and DNA sequences. They can be used to search sequence data bases, evaluate similarity scores, and identify periodic structures based on local sequence similarity. The FASTA program is a more sensitive derivative of the FASTP program, which can be used to search protein or DNA sequence data bases and can compare a protein sequence to a DNA sequence data base by translating the DNA data base as it is searched. FASTA includes an additional step in the calculation of the initial pairwise similarity score that allows multiple regions of similarity to be joined to increase the score of related sequences. The RDF2 program can be used to evaluate the significance of similarity scores using a shuffling method that preserves local sequence composition. The LFASTA program can display all the regions of local similarity between two sequences with scores greater than a threshold, using the same scoring parameters and a similar alignment algorithm; these local similarities can be displayed as a "graphic matrix" plot or as individual alignments. In addition, these programs have been generalized to allow comparison of DNA or protein sequences based on a variety of alternative scoring matrices.

Characterization of constitutive CTCF/cohesin loci: a possible role in establishing topological domains in mammalian genomes

BACKGROUND Recent studies suggested that human/mammalian genomes are divided into large, discrete domains that are units of chromosome organization. CTCF, a CCCTC binding factor, has a diverse role in genome regulation including transcriptional regulation, chromosome-boundary insulation, DNA replication, and chromatin packaging. It remains unclear whether a subset of CTCF binding sites plays a functional role in establishing/maintaining chromatin topological domains. RESULTS We systematically analysed the genomic, transcriptomic and epigenetic profiles of the CTCF binding sites in 56 human cell lines from ENCODE. We identified ~24,000 CTCF sites (referred to as constitutive sites) that were bound in more than 90% of the cell lines. Our analysis revealed: 1) constitutive CTCF loci were located in constitutive open chromatin and often co-localized with constitutive cohesin loci; 2) most constitutive CTCF loci were distant from transcription start sites and lacked CpG islands but were enriched with the full-spectrum CTCF motifs: a recently reported 33/34-mer and two other potentially novel (22/26-mer); 3) more importantly, most constitutive CTCF loci were present in CTCF-mediated chromatin interactions detected by ChIA-PET and these pair-wise interactions occurred predominantly within, but not between, topological domains identified by Hi-C.   CONCLUSIONS Our results suggest that the constitutive CTCF sites may play a role in organizing/maintaining the recently identified topological domains that are common across most human cells.

lumi: a pipeline for processing Illumina microarray.

UNLABELLED Illumina microarray is becoming a popular microarray platform. The BeadArray technology from Illumina makes its preprocessing and quality control different from other microarray technologies. Unfortunately, most other analyses have not taken advantage of the unique properties of the BeadArray system, and have just incorporated preprocessing methods originally designed for Affymetrix microarrays. lumi is a Bioconductor package especially designed to process the Illumina microarray data. It includes data input, quality control, variance stabilization, normalization and gene annotation portions. In specific, the lumi package includes a variance-stabilizing transformation (VST) algorithm that takes advantage of the technical replicates available on every Illumina microarray. Different normalization method options and multiple quality control plots are provided in the package. To better annotate the Illumina data, a vendor independent nucleotide universal identifier (nuID) was devised to identify the probes of Illumina microarray. The nuID annotation packages and output of lumi processed results can be easily integrated with other Bioconductor packages to construct a statistical data analysis pipeline for Illumina data. AVAILABILITY The lumi Bioconductor package, www.bioconductor.org

Identification and characterization of a common set of complex I assembly intermediates in mitochondria from patients with complex I deficiency.

Deficiencies in the activity of complex I (NADH: ubiquinone oxidoreductase) are an important cause of human mitochondrial disease. Complex I is composed of at least 46 structural subunits that are encoded in both nuclear and mitochondrial DNA. Enzyme deficiency can result from either impaired catalytic efficiency or an inability to assemble the holoenzyme complex; however, the assembly process remains poorly understood. We have used two-dimensional Blue-Native/SDS gel electrophoresis and a panel of 11 antibodies directed against structural subunits of the enzyme to investigate complex I assembly in the muscle mitochondria from four patients with complex I deficiency caused by either mitochondrial or nuclear gene defects. Immunoblot analyses of second dimension denaturing gels identified seven distinct complex I subcomplexes in the patients studied, five of which could also be detected in nondenaturing gels in the first dimension. Although the abundance of these intermediates varied among the different patients, a common constellation of subcomplexes was observed in all cases. A similar profile of subcomplexes was present in a human/mouse hybrid fibroblast cell line with a severe complex I deficiency due to an almost complete lack of assembly of the holoenzyme complex. The finding that diverse causes of complex I deficiency produce a similar pattern of complex I subcomplexes suggests that these are intermediates in the assembly of the holoenzyme complex. We propose a possible assembly pathway for the complex, which differs significantly from that proposed for Neurospora, the current model for complex I assembly.

Frozen robust multiarray analysis (fRMA).

Robust multiarray analysis (RMA) is the most widely used preprocessing algorithm for Affymetrix and Nimblegen gene expression microarrays. RMA performs background correction, normalization, and summarization in a modular way. The last 2 steps require multiple arrays to be analyzed simultaneously. The ability to borrow information across samples provides RMA various advantages. For example, the summarization step fits a parametric model that accounts for probe effects, assumed to be fixed across arrays, and improves outlier detection. Residuals, obtained from the fitted model, permit the creation of useful quality metrics. However, the dependence on multiple arrays has 2 drawbacks: (1) RMA cannot be used in clinical settings where samples must be processed individually or in small batches and (2) data sets preprocessed separately are not comparable. We propose a preprocessing algorithm, frozen RMA (fRMA), which allows one to analyze microarrays individually or in small batches and then combine the data for analysis. This is accomplished by utilizing information from the large publicly available microarray databases. In particular, estimates of probe-specific effects and variances are precomputed and frozen. Then, with new data sets, these are used in concert with information from the new arrays to normalize and summarize the data. We find that fRMA is comparable to RMA when the data are analyzed as a single batch and outperforms RMA when analyzing multiple batches. The methods described here are implemented in the R package fRMA and are currently available for download from the software section of http://rafalab.jhsph.edu.

Structural features of the interaction of the 3'-untranslated region of mRNA containing exosomal RNA-specific motifs with YB-1, a potential mediator of mRNA sorting.

We have previously shown that YB-1 is the only protein of the HEK293 cell cytoplasmic (S100) extract that specifically interacts with RNA hairpins each containing one of the motifs ACCAGCCU (1), CAGUGAGC (2) and UAAUCCCA (3), which had been identified as often found in exosomal RNA and proposed as potential cis-acting elements targeting RNAs into exosomes. Here we explored the interactions of YB-1 with a fragment of the 3'-untranslated region (UTR) of septin 14 mRNA (SEPT14 RNA), which contains all three motifs. We demonstrated the occurrence of YB-1 among proteins pulled down from the HEK293 S100 extract using biotinylated SEPT14 RNA. With recombinant YB-1, it was found that SEPT14 RNA can bind up to 5 moles of protein per mole of RNA in a cooperative manner, which was shown to be mainly facilitated by the presence of the above motifs. RNA hairpins with motifs 1 and 2 competed with SEPT14 RNA for binding to the protein, whereas that with motif 3 was less competitive, in accordance with the affinity of YB-1 for these RNA hairpins. With YB-1-bound RNA, nucleotides protected from attack by hydroxyl radicals were revealed in all three motifs, although hairpins with motif 2 and especially with motif 1 contained many protected nucleotides outside the motifs, suggesting that the specific environments of these motifs contribute significantly to the YB-1 binding. An analysis of the environments of motifs 1-3 in the HEK293 cell mRNA 3' UTRs gained from RNA-seq data led us to conclude that the primary binding sites of YB-1 in the 3' UTRs are hairpins containing some part of the motif along with its specific surroundings; the consensus sequences of these hairpins were derived. Thus, our findings provide a new understanding of the structural basis of the interactions between YB-1 and mRNAs carrying the aforementioned motifs.

Specific HDV RNA-templated transcription by pol II in vitro.

RNA polymerase II is implicated in the RNA-templated RNA synthesis during replication of viroids and Hepatitis Delta Virus (HDV); however, neither the RNA template nor protein factor requirements for this process are well defined. We have developed an in vitro transcription system based on HeLa cell nuclear extract (NE), in which a segment of antigenomic RNA corresponding to the left-hand tip region of the HDV rod-like structure serves as a template for efficient and highly specific RNA synthesis. Accumulation of the unique RNA product is highly sensitive to alpha-amanitin in HeLa NE and only partially sensitive to this drug in NE from PMG cells that contain an allele of the alpha-amanitin-resistant subunit of pol II, strongly suggesting pol II involvement in this reaction. Detailed analysis of the RNA product revealed that it represents a chimeric molecule composed of a newly synthesized transcript covalently attached to the 5' half of the RNA template. Selection of the start site for transcription is remarkably specific and depends on the secondary structure of the RNA template, rather than on its primary sequence. Some features of this reaction resemble the RNA cleavage-extension process observed for pol II-arrested complexes in vitro. A possible involvement of the described reaction in HDV replication is discussed.

Mitotic chromosomes fold by condensin-dependent helical winding of chromatin loop arrays

During mitosis, chromosomes fold into compacted rod shaped structures. We combined imaging and Hi-C of synchronous DT40 cell cultures with polymer simulations to determine how interphase chromosomes are converted into compressed arrays of loops characteristic of mitotic chromosomes. We found that the interphase organization is disassembled within minutes of prophase entry and by late prophase chromosomes are already folded as arrays of consecutive loops. During prometaphase, this array reorganizes to form a helical arrangement of nested loops. Polymer simulations reveal that Hi-C data are inconsistent with solenoidal coiling of the entire chromatid, but instead suggest a centrally located helically twisted axis from which consecutive loops emanate as in a spiral staircase. Chromosomes subsequently shorten through progressive helical winding, with the numbers of loops per turn increasing so that the size of a helical turn grows from around 3 Mb (~40 loops) to ~12 Mb (~150 loops) in fully condensed metaphase chromosomes. Condensin is essential to disassemble the interphase chromatin conformation. Analysis of mutants revealed differing roles for condensin I and II during these processes. Either condensin can mediate formation of loop arrays. However, condensin II was required for helical winding during prometaphase, whereas condensin I modulated the size and arrangement of loops inside the helical turns. These observations identify a mitotic chromosome morphogenesis pathway in which folding of linear loop arrays produces long thin chromosomes during prophase that then shorten by progressive growth of loops and helical winding during prometaphase.

Dnmt1-Independent CG Methylation Contributes to Nucleosome Positioning in Diverse Eukaryotes

Dnmt1 epigenetically propagates symmetrical CG methylation in many eukaryotes. Their genomes are typically depleted of CG dinucleotides because of imperfect repair of deaminated methylcytosines. Here, we extensively survey diverse species lacking Dnmt1 and show that, surprisingly, symmetrical CG methylation is nonetheless frequently present and catalyzed by a different DNA methyltransferase family, Dnmt5. Numerous Dnmt5-containing organisms that diverged more than a billion years ago exhibit clustered methylation, specifically in nucleosome linkers. Clustered methylation occurs at unprecedented densities and directly disfavors nucleosomes, contributing to nucleosome positioning between clusters. Dense methylation is enabled by a regime of genomic sequence evolution that enriches CG dinucleotides and drives the highest CG frequencies known. Species with linker methylation have small, transcriptionally active nuclei that approach the physical limits of chromatin compaction. These features constitute a previously unappreciated genome architecture, in which dense methylation influences nucleosome positions, likely facilitating nuclear processes under extreme spatial constraints.

Using the yeast gene deletion collection to customize gene expression.

Substitute Teacher Despite the relative ease of genome manipulation in S. cerevisiae, researchers are always looking to learn still more convenient and rapid methods for substituting yeast promoters. Replacing a gene's native promoter with a heterolo-gous promoter of choice allows regulated expression and simplifies the task of discerning functional relevance. Although a host of clever chromosomal insertion strategies have been described over the years, the advent of the S. cerevisiae Genome Deletion Project provides an incredible resource for a further streamlined workflow. The strategy, explained by Liko et al. on p. 728 is appealingly simple. The genome deletion project resulted in a collection of strains in which a single ORF is replaced with a kanamycin resistance module. Although the purpose of the collection is to have a comprehensive resource of essentially all possible knockouts, the authors point out that for almost any given yeast promoter of interest there will be a strain in which the ORF imm...

A conserved mechanism for centromeric nucleosome recognition by centromere protein CENP-C.

Chromosome segregation during mitosis requires assembly of the kinetochore complex at the centromere. Kinetochore assembly depends on specific recognition of the histone variant CENP-A in the centromeric nucleosome by centromere protein C (CENP-C). We have defined the determinants of this recognition mechanism and discovered that CENP-C binds a hydrophobic region in the CENP-A tail and docks onto the acidic patch of histone H2A and H2B. We further found that the more broadly conserved CENP-C motif uses the same mechanism for CENP-A nucleosome recognition. Our findings reveal a conserved mechanism for protein recruitment to centromeres and a histone recognition mode whereby a disordered peptide binds the histone tail through hydrophobic interactions facilitated by nucleosome docking.

In vivo microtubules are copolymers of available beta-tubulin isotypes: localization of each of six vertebrate beta-tubulin isotypes using polyclonal antibodies elicited by synthetic peptide antigens

beta-Tubulin is encoded in the genomes of higher animals by a small multigene family comprising approximately seven functional genes. These genes produce a family of closely related, but distinct polypeptide isotypes that are distinguished principally by sequences within the approximately 15 carboxy-terminal amino acid residues. By immunizing rabbits with chemically synthesized peptides corresponding to these variable domain sequences, we have now prepared polyclonal antibodies specific for each of six distinct isotypes. Specificity of each antiserum has been demonstrated unambiguously by antibody binding to bacterially produced, cloned proteins representing each isotype and by the inhibition of such binding by preincubation of each antiserum only with the immunizing peptide and not with heterologous peptides. Protein blotting of known amounts of cloned, isotypically pure polypeptides has permitted accurate quantitative measurement of the amount of each beta-tubulin isotype present in the soluble and polymer forms in various cells, but has not revealed a bias for preferential assembly of any isotype. Localization of each isotype in three different cell types using indirect immunofluorescence has demonstrated that in vivo each class of microtubules distinguishable by light microscopy is assembled as copolymers of all isotypes expressed in a single cell.

Interplay between microtubule bundling and sorting factors ensures acentriolar spindle stability during C. elegans oocyte meiosis

In many species, oocyte meiosis is carried out in the absence of centrioles. As a result, microtubule organization, spindle assembly, and chromosome segregation proceed by unique mechanisms. Here, we report insights into the principles underlying this specialized form of cell division, through studies of C. elegans KLP-15 and KLP-16, two highly homologous members of the kinesin-14 family of minus-end-directed kinesins. These proteins localize to the acentriolar oocyte spindle and promote microtubule bundling during spindle assembly; following KLP-15/16 depletion, microtubule bundles form but then collapse into a disorganized array. Surprisingly, despite this defect we found that during anaphase, microtubules are able to reorganize into a bundled array that facilitates chromosome segregation. This phenotype therefore enabled us to identify factors promoting microtubule organization during anaphase, whose contributions are normally undetectable in wild-type worms; we found that SPD-1 (PRC1) bundles microtubules and KLP-18 (kinesin-12) likely sorts those bundles into a functional orientation capable of mediating chromosome segregation. Therefore, our studies have revealed an interplay between distinct mechanisms that together promote spindle formation and chromosome segregation in the absence of structural cues such as centrioles.

Extended-connectivity fingerprints.

Extended-connectivity fingerprints (ECFPs) are a novel class of topological fingerprints for molecular characterization. Historically, topological fingerprints were developed for substructure and similarity searching. ECFPs were developed specifically for structure-activity modeling. ECFPs are circular fingerprints with a number of useful qualities: they can be very rapidly calculated; they are not predefined and can represent an essentially infinite number of different molecular features (including stereochemical information); their features represent the presence of particular substructures, allowing easier interpretation of analysis results; and the ECFP algorithm can be tailored to generate different types of circular fingerprints, optimized for different uses. While the use of ECFPs has been widely adopted and validated, a description of their implementation has not previously been presented in the literature.

Functional coupling of RNAP II transcription to spliceosome assembly.

The pathway of gene expression in higher eukaryotes involves a highly complex network of physical and functional interactions among the different machines involved in each step of the pathway. Here we established an efficient in vitro system to determine how RNA polymerase II (RNAP II) transcription is functionally coupled to pre-mRNA splicing. Strikingly, our data show that nascent pre-messenger RNA (pre-mRNA) synthesized by RNAP II is immediately and quantitatively directed into the spliceosome assembly pathway. In contrast, nascent pre-mRNA synthesized by T7 RNA polymerase is quantitatively assembled into the nonspecific H complex, which consists of heterogeneous nuclear ribonucleoprotein (hnRNP) proteins and is inhibitory for spliceosome assembly. Consequently, RNAP II transcription results in a dramatic increase in both the kinetics of splicing and overall yield of spliced mRNA relative to that observed for T7 transcription. We conclude that RNAP II mediates the functional coupling of transcription to splicing by directing the nascent pre-mRNA into spliceosome assembly, thereby bypassing interaction of the pre-mRNA with the inhibitory hnRNP proteins.

Inferring nucleosome positions with their histone mark annotation from ChIP data

MOTIVATION The nucleosome is the basic repeating unit of chromatin. It contains two copies each of the four core histones H2A, H2B, H3 and H4 and about 147 bp of DNA. The residues of the histone proteins are subject to numerous post-translational modifications, such as methylation or acetylation. Chromatin immunoprecipitiation followed by sequencing (ChIP-seq) is a technique that provides genome-wide occupancy data of these modified histone proteins, and it requires appropriate computational methods. RESULTS We present NucHunter, an algorithm that uses the data from ChIP-seq experiments directed against many histone modifications to infer positioned nucleosomes. NucHunter annotates each of these nucleosomes with the intensities of the histone modifications. We demonstrate that these annotations can be used to infer nucleosomal states with distinct correlations to underlying genomic features and chromatin-related processes, such as transcriptional start sites, enhancers, elongation by RNA polymerase II and chromatin-mediated repression. Thus, NucHunter is a versatile tool that can be used to predict positioned nucleosomes from a panel of histone modification ChIP-seq experiments and infer distinct histone modification patterns associated to different chromatin states. AVAILABILITY The software is available at http://epigen.molgen.mpg.de/nuchunter/.

The non-catalytic function of XPG protein during dual incision in human nucleotide excision repair.

XPG is a member of the FEN-1 structure-specific endonuclease family. It has 3'-junction cutting activity on bubble substrates and makes the 3'-incision in the human dual incision (excision nuclease) repair system. To investigate the precise role of XPG in nucleotide excision repair, we mutagenized two amino acid residues thought to be involved in DNA binding and catalysis, overproduced the mutant proteins using a baculovirus/insect cell system, and purified and characterized the mutant proteins. The mutation D77A had a modest effect on junction cutting and excision activity and gave rise to uncoupled 5'-incision by mammalian cell-free extracts. The D812A mutation completely abolished the junction cutting and 3'-incision activities of XPG, but the excision nuclease reconstituted with XPG (D812A) carried out normal 5'-incision at the 23rd-24th phosphodiester bonds 5' to a (6-4) photoproduct without producing any 3'-incision. It is concluded that Asp-812 is an active site residue of XPG and that in addition to making the 3'-incision, the physical presence of XPG in the protein-DNA complex is required non-catalytically for subsequent 5'-incision by XPF-ERCC1.

A 3D Map of the Yeast Kinetochore Reveals the Presence of Core and Accessory Centromere-Specific Histone

The budding yeast kinetochore is ~68 nm in length with a diameter slightly larger than a 25 nm microtubule. The kinetochores from the 16 chromosomes are organized in a stereotypic cluster encircling central spindle microtubules. Quantitative analysis of the inner kinetochore cluster (Cse4, COMA) reveals structural features not apparent in singly attached kinetochores. The cluster of Cse4-containing kinetochores is physically larger perpendicular to the spindle axis relative to the cluster of Ndc80 molecules. If there was a single Cse4 (molecule or nucleosome) at the kinetochore attached to each microtubule plus end, the cluster of Cse4 would appear geometrically identical to Ndc80. Thus, the structure of the inner kinetochore at the surface of the chromosomes remains unsolved. We have used point fluorescence microscopy and statistical probability maps to deduce the two-dimensional mean position of representative components of the yeast kinetochore relative to the mitotic spindle in metaphase. Comparison of the experimental images to three-dimensional architectures from convolution of mathematical models reveals a pool of Cse4 radially displaced from Cse4 at the kinetochore and kinetochore microtubule plus ends. The pool of displaced Cse4 can be experimentally depleted in mRNA processing pat1Δ or xrn1Δ mutants. The peripheral Cse4 molecules do not template outer kinetochore components. This study suggests an inner kinetochore plate at the centromere-microtubule interface in budding yeast and yields information on the number of Ndc80 molecules at the microtubule attachment site.

P446L-importin-beta inhibits nuclear envelope assembly by sequestering nuclear envelope assembly factors to the microtubules.

The P446L mutant Drosophila importin-beta (P446L-imp-beta) has been reported to prohibit--in dominant negative fashion--nuclear envelope (NE) assembly. Along elucidating the mode of action of P446L-imp-beta we studied in vitro NE assembly on Sepharose beads. While Drosophila embryo extracts support NE assembly over Sepharose beads coated with Ran, NE assembly does not take place in extracts supplied with exogenous P446L-imp-beta. A NE also forms over importin-beta-coated beads. Surprisingly, when immobilized to Sepharose beads P446L-imp-beta as efficiently recruits NE vesicles as normal importin-beta. The discrepancy in behavior of cytoplasmic and bead-bound P446L-imp-beta appears to be related to icreased--as compared to normal importin-beta--microtubule (MT) binding ability of P446L-imp-beta. While wild-type importin-beta is able to bind MTs and the binding decreases upon RanGTP interaction, P446L-imp-beta cannot be removed from the MTs by RanGTP. P446L-imp-beta, like normal importin-beta, binds some types of the nucleoporins that have been known to be required for NE assembly at the end of mitosis. It appears that the inhibitory effect of P446L-imp-beta on NE assembly is caused by sequestering some of the nucleoporins required for NE assembly to the MTs.

Expression profiles during dedifferentiation in newt lens regeneration revealed by expressed sequence tags

PURPOSE The adult newt can regenerate lens from pigmented epithelial cells (PECs) of the dorsal iris via dedifferentiation. The purpose of this research is to obtain sequence resources for a newt lens regeneration study and to obtain insights of dedifferentiation at the molecular level. METHODS mRNA was purified from iris during dedifferentiation and its cDNA library was constructed. From the cDNA library 10,449 clones were sequenced and analyzed. RESULTS From 10,449 reads, 780 contigs and 1,666 singlets were annotated. The presence of several cancer- and apoptosis-related genes during newt dedifferentiation was revealed. Moreover, several candidate genes, which might participate in reprogramming during dedifferentiation, were also found. CONCLUSIONS The expression of cancer- and apoptosis-related genes could be hallmarks during dedifferentiation. The expression sequence tag (EST) resource is useful for the future study of newt dedifferentiation, and the sequence information is available in GenBank (accession numbers; FS290155-FS300559).

The Ndc80 Kinetochore Complex Forms Load-Bearing Attachments to Dynamic Microtubule Tips via Biased Diffusion

Kinetochores couple chromosomes to the assembling and disassembling tips of microtubules, a dynamic behavior that is fundamental to mitosis in all eukaryotes but poorly understood. Genetic, biochemical, and structural studies implicate the Ndc80 complex as a direct point of contact between kinetochores and microtubules, but these approaches provide only a static view. Here, using techniques for manipulating and tracking individual molecules in vitro, we demonstrate that the Ndc80 complex is capable of forming the dynamic, load-bearing attachments to assembling and disassembling tips required for coupling in vivo. We also establish that Ndc80-based coupling likely occurs through a biased diffusion mechanism and that this activity is conserved from yeast to humans. Our findings demonstrate how an ensemble of Ndc80 complexes may provide the combination of plasticity and strength that allows kinetochores to maintain load-bearing tip attachments during both microtubule assembly and disassembly.

JTK_CYCLE: an efficient nonparametric algorithm for detecting rhythmic components in genome-scale data sets.

Circadian rhythms are oscillations of physiology, behavior, and metabolism that have period lengths near 24 hours. In several model organisms and humans, circadian clock genes have been characterized and found to be transcription factors. Because of this, researchers have used microarrays to characterize global regulation of gene expression and algorithmic approaches to detect cycling. This article presents a new algorithm, JTK_CYCLE, designed to efficiently identify and characterize cycling variables in large data sets. Compared with COSOPT and the Fisher's G test, two commonly used methods for detecting cycling transcripts, JTK_CYCLE distinguishes between rhythmic and nonrhythmic transcripts more reliably and efficiently. JTK_CYCLE's increased resistance to outliers results in considerably greater sensitivity and specificity. Moreover, JTK_CYCLE accurately measures the period, phase, and amplitude of cycling transcripts, facilitating downstream analyses. Finally, JTK_CYCLE is several orders of magnitude faster than COSOPT, making it ideal for large-scale data sets. JTK_CYCLE was used to analyze legacy data sets including NIH3T3 cells, which have comparatively low amplitude oscillations. JTK_CYCLE's improved power led to the identification of a novel cluster of RNA-interacting genes whose abundance is under clear circadian regulation. These data suggest that JTK_CYCLE is an ideal tool for identifying and characterizing oscillations in genome-scale data sets.

Histone chaperones: 30 years from isolation to elucidation of the mechanisms of nucleosome assembly and disassembly

Some three decades have passed since the discovery of nucleosomes in 1974 and the first isolation of a histone chaperone in 1978. While various types of histone chaperones have been isolated and functionally analyzed, the elementary processes of nucleosome assembly and disassembly have been less well characterized. Recently, the tertiary structure of a hetero-trimeric complex composed of the histone chaperone CIA/ASF1 and the histone H3-H4 dimer was determined, and this complex was proposed to be an intermediate in nucleosome assembly and disassembly reactions. In addition, CIA alone was biochemically shown to dissociate the histone (H3-H4)2 tetramer into two histone H3-H4 dimers. This activity suggested that CIA regulates the semi-conservative replication of nucleosomes. Here, we provide an overview of prominent histone chaperones with the goal of elucidating the mechanisms that preserve and modify epigenetic information. We also discuss the reactions involved in nucleosome assembly and disassembly.

Increased dendritic complexity and axonal length in cultured mouse cortical neurons overexpressing methyl-CpG-binding protein MeCP2

Rett syndrome is caused by loss-of-function mutations in the gene encoding the methyl DNA-binding factor MeCP2. As brain mass and neuronal complexity tend to be diminished in Rett patients, we tested whether MeCP2 directly influences the morphological complexity of developing neurons. Our results show that cultured mouse neurons overexpressing MeCP2beta (MECP2A) develop more complex morphologies, having longer axonal and dendritic processes, and an increased number of axonal and dendritic terminal endings. We then tested whether overexpressing a mutant form of MeCP2beta lacking its carboxyl terminus would elicit the same effects. Interestingly, while neurons overexpressing this mutant failed to enhance axonal and dendritic process elongation, the complexity of their axonal and dendritic processes remained significantly elevated. Taken together, these data support the hypothesis that MeCP2 directly regulates neuronal maturation and/or synaptogenesis, and provides evidence that MeCP2 may influence neuritic elongation and process branching through different mechanisms.

Haplotype estimation using sequencing reads.

High-throughput sequencing technologies produce short sequence reads that can contain phase information if they span two or more heterozygote genotypes. This information is not routinely used by current methods that infer haplotypes from genotype data. We have extended the SHAPEIT2 method to use phase-informative sequencing reads to improve phasing accuracy. Our model incorporates the read information in a probabilistic model through base quality scores within each read. The method is primarily designed for high-coverage sequence data or data sets that already have genotypes called. One important application is phasing of single samples sequenced at high coverage for use in medical sequencing and studies of rare diseases. Our method can also use existing panels of reference haplotypes. We tested the method by using a mother-father-child trio sequenced at high-coverage by Illumina together with the low-coverage sequence data from the 1000 Genomes Project (1000GP). We found that use of phase-informative reads increases the mean distance between switch errors by 22% from 274.4 kb to 328.6 kb. We also used male chromosome X haplotypes from the 1000GP samples to simulate sequencing reads with varying insert size, read length, and base error rate. When using short 100 bp paired-end reads, we found that using mixtures of insert sizes produced the best results. When using longer reads with high error rates (5-20 kb read with 4%-15% error per base), phasing performance was substantially improved.

Trimmomatic: a flexible trimmer for Illumina sequence data

MOTIVATION Although many next-generation sequencing (NGS) read preprocessing tools already existed, we could not find any tool or combination of tools that met our requirements in terms of flexibility, correct handling of paired-end data and high performance. We have developed Trimmomatic as a more flexible and efficient preprocessing tool, which could correctly handle paired-end data. RESULTS The value of NGS read preprocessing is demonstrated for both reference-based and reference-free tasks. Trimmomatic is shown to produce output that is at least competitive with, and in many cases superior to, that produced by other tools, in all scenarios tested. AVAILABILITY AND IMPLEMENTATION Trimmomatic is licensed under GPL V3. It is cross-platform (Java 1.5+ required) and available at http://www.usadellab.org/cms/index.php?page=trimmomatic   CONTACT usadel@bio1.rwth-aachen.de   SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.

Cutadapt removes adapter sequences from high-throughput sequencing reads

When small RNA is sequenced on current sequencing machines, the resulting reads are usually longer than the RNA and therefore contain parts of the 3' adapter. That adapter must be found and removed error-tolerantly from each read before read mapping. Previous solutions are either hard to use or do not offer required features, in particular support for color space data. As an easy to use alternative, we developed the command-line tool cutadapt, which supports 454, Illumina and SOLiD (color space) data, offers two adapter trimming algorithms, and has other useful features. Cutadapt, including its MIT-licensed source code, is available for download at http://code.google.com/p/cutadapt/

The binding of maize DHN1 to lipid vesicles. Gain of structure and lipid specificity.

Dehydrins (DHNs; late embryogenesis abundant D-11) are a family of plant proteins induced in response to abiotic stresses such as drought, low temperature, and salinity or during the late stages of embryogenesis. Spectral and thermal properties of these proteins in purified form suggest that they are "intrinsically unstructured. " However, DHNs contain at least one copy of a consensus 15-amino acid sequence, the "K segment," which resembles a class A2 amphipathic alpha-helical, lipid-binding domain found in other proteins such as apolipoproteins and alpha-synuclein. The presence of the K segment raises the question of whether DHNs bind lipids, bilayers, or phospholipid vesicles. Here, we show that maize (Zea mays) DHN DHN1 can bind to lipid vesicles that contain acidic phospholipids. We also observe that DHN1 binds more favorably to vesicles of smaller diameter than to larger vesicles, and that the association of DHN1 with vesicles results in an apparent increase of alpha-helicity of the protein. Therefore, DHNs, and presumably somewhat similar plant stress proteins in the late embryogenesis abundant and cold-regulated classes may undergo function-related conformational changes at the water/membrane interface, perhaps related to the stabilization of vesicles or other endomembrane structures under stress conditions.

Eic1 links Mis18 with the CCAN/Mis6/Ctf19 complex to promote CENP-A assembly

CENP-A chromatin forms the foundation for kinetochore assembly. Replication-independent incorporation of CENP-A at centromeres depends on its chaperone HJURP(Scm3), and Mis18 in vertebrates and fission yeast. The recruitment of Mis18 and HJURP(Scm3) to centromeres is cell cycle regulated. Vertebrate Mis18 associates with Mis18BP1(KNL2), which is critical for the recruitment of Mis18 and HJURP(Scm3). We identify two novel fission yeast Mis18-interacting proteins (Eic1 and Eic2), components of the Mis18 complex. Eic1 is essential to maintain Cnp1(CENP-A) at centromeres and is crucial for kinetochore integrity; Eic2 is dispensable. Eic1 also associates with Fta7(CENP-Q/Okp1), Cnl2(Nkp2) and Mal2(CENP-O/Mcm21), components of the constitutive CCAN/Mis6/Ctf19 complex. No Mis18BP1(KNL2) orthologue has been identified in fission yeast, consequently it remains unknown how the key Cnp1(CENP-A) loading factor Mis18 is recruited. Our findings suggest that Eic1 serves a function analogous to that of Mis18BP1(KNL2), thus representing the functional counterpart of Mis18BP1(KNL2) in fission yeast that connects with a module within the CCAN/Mis6/Ctf19 complex to allow the temporally regulated recruitment of the Mis18/Scm3(HJURP) Cnp1(CENP-A) loading factors. The novel interactions identified between CENP-A loading factors and the CCAN/Mis6/Ctf19 complex are likely to also contribute to CENP-A maintenance in other organisms.

Targeted proteomics reveals compositional dynamics of 60S pre-ribosomes after nuclear export

Construction and intracellular targeting of eukaryotic pre-ribosomal particles involve a multitude of diverse transiently associating trans-acting assembly factors, energy-consuming enzymes, and transport factors. The ability to rapidly and reliably measure co-enrichment of multiple factors with maturing pre-ribosomal particles presents a major biochemical bottleneck towards revealing their function and the precise contribution of >50 energy-consuming steps that drive ribosome assembly. Here, we devised a workflow that combines genetic trapping, affinity-capture, and selected reaction monitoring mass spectrometry (SRM-MS), to overcome this deficiency. We exploited this approach to interrogate the dynamic proteome of pre-60S particles after nuclear export. We uncovered assembly factors that travel with pre-60S particles to the cytoplasm, where they are released before initiating translation. Notably, we identified a novel shuttling factor that facilitates nuclear export of pre-60S particles. Capturing and quantitating protein interaction networks of trapped intermediates of macromolecular complexes by our workflow is a reliable discovery tool to unveil dynamic processes that contribute to their in vivo assembly and transport.

Ribosome Profiling Provides Evidence that Large Noncoding RNAs Do Not Encode Proteins

Large noncoding RNAs are emerging as an important component in cellular regulation. Considerable evidence indicates that these transcripts act directly as functional RNAs rather than through an encoded protein product. However, a recent study of ribosome occupancy reported that many large intergenic ncRNAs (lincRNAs) are bound by ribosomes, raising the possibility that they are translated into proteins. Here, we show that classical noncoding RNAs and 5' UTRs show the same ribosome occupancy as lincRNAs, demonstrating that ribosome occupancy alone is not sufficient to classify transcripts as coding or noncoding. Instead, we define a metric based on the known property of translation whereby translating ribosomes are released upon encountering a bona fide stop codon. We show that this metric accurately discriminates between protein-coding transcripts and all classes of known noncoding transcripts, including lincRNAs. Taken together, these results argue that the large majority of lincRNAs do not function through encoded proteins.

Specificity of the dRP/AP lyase of Ku promotes nonhomologous end joining (NHEJ) fidelity at damaged ends.

Nonhomologous end joining (NHEJ) is essential for efficient repair of chromosome breaks. However, the NHEJ ligation step is often obstructed by break-associated nucleotide damage, including base loss (abasic site or 5'-dRP/AP sites). Ku, a 5'-dRP/AP lyase, can excise such damage at ends in preparation for the ligation step. We show here that this activity is greatest if the abasic site is within a short 5' overhang, when this activity is necessary and sufficient to prepare such termini for ligation. In contrast, Ku is less active near 3' strand termini, where excision would leave a ligation-blocking α,β-unsaturated aldehyde. The Ku AP lyase activity is also strongly suppressed by as little as two paired bases 5' of the abasic site. Importantly, in vitro end joining experiments show that abasic sites significantly embedded in double-stranded DNA do not block the NHEJ ligation step. Suppression of the excision activity of Ku in this context therefore is not essential for ligation and further helps NHEJ retain terminal sequence in junctions. We show that the DNA between the 5' terminus and the abasic site can also be retained in junctions formed by cellular NHEJ, indicating that these sites are at least partly resistant to other abasic site-cleaving activities as well. High levels of the 5'-dRP/AP lyase activity of Ku are thus restricted to substrates where excision of an abasic site is required for ligation, a degree of specificity that promotes more accurate joining.

Insights into G-quadruplex specific recognition by the DEAH-box helicase RHAU: Solution structure of a peptide-quadruplex complex.

Four-stranded nucleic acid structures called G-quadruplexes have been associated with important cellular processes, which should require G-quadruplex-protein interaction. However, the structural basis for specific G-quadruplex recognition by proteins has not been understood. The DEAH (Asp-Glu-Ala-His) box RNA helicase associated with AU-rich element (RHAU) (also named DHX36 or G4R1) specifically binds to and resolves parallel-stranded G-quadruplexes. Here we identified an 18-amino acid G-quadruplex-binding domain of RHAU and determined the structure of this peptide bound to a parallel DNA G-quadruplex. Our structure explains how RHAU specifically recognizes parallel G-quadruplexes. The peptide covers a terminal guanine base tetrad (G-tetrad), and clamps the G-quadruplex using three-anchor-point electrostatic interactions between three positively charged amino acids and negatively charged phosphate groups. This binding mode is strikingly similar to that of most ligands selected for specific G-quadruplex targeting. Binding to an exposed G-tetrad represents a simple and efficient way to specifically target G-quadruplex structures.

Assessing sequence comparison methods with reliable structurally identified distant evolutionary relationships.

Pairwise sequence comparison methods have been assessed using proteins whose relationships are known reliably from their structures and functions, as described in the SCOP database [Murzin, A. G., Brenner, S. E., Hubbard, T. & Chothia C. (1995) J. Mol. Biol. 247, 536-540]. The evaluation tested the programs BLAST [Altschul, S. F., Gish, W., Miller, W., Myers, E. W. & Lipman, D. J. (1990). J. Mol. Biol. 215, 403-410], WU-BLAST2 [Altschul, S. F. & Gish, W. (1996) Methods Enzymol. 266, 460-480], FASTA [Pearson, W. R. & Lipman, D. J. (1988) Proc. Natl. Acad. Sci. USA 85, 2444-2448], and SSEARCH [Smith, T. F. & Waterman, M. S. (1981) J. Mol. Biol. 147, 195-197] and their scoring schemes. The error rate of all algorithms is greatly reduced by using statistical scores to evaluate matches rather than percentage identity or raw scores. The E-value statistical scores of SSEARCH and FASTA are reliable: the number of false positives found in our tests agrees well with the scores reported. However, the P-values reported by BLAST and WU-BLAST2 exaggerate significance by orders of magnitude. SSEARCH, FASTA ktup = 1, and WU-BLAST2 perform best, and they are capable of detecting almost all relationships between proteins whose sequence identities are >30%. For more distantly related proteins, they do much less well; only one-half of the relationships between proteins with 20-30% identity are found. Because many homologs have low sequence similarity, most distant relationships cannot be detected by any pairwise comparison method; however, those which are identified may be used with confidence.

The docking domain of histone H2A is required for H1 binding and RSC-mediated nucleosome remodeling

Histone variants within the H2A family show high divergences in their C-terminal regions. In this work, we have studied how these divergences and in particular, how a part of the H2A COOH-terminus, the docking domain, is implicated in both structural and functional properties of the nucleosome. Using biochemical methods in combination with Atomic Force Microscopy and Electron Cryo-Microscopy, we show that the H2A-docking domain is a key structural feature within the nucleosome. Deletion of this domain or replacement with the incomplete docking domain from the variant H2A.Bbd results in significant structural alterations in the nucleosome, including an increase in overall accessibility to nucleases, un-wrapping of ∼10 bp of DNA from each end of the nucleosome and associated changes in the entry/exit angle of DNA ends. These structural alterations are associated with a reduced ability of the chromatin remodeler RSC to both remodel and mobilize the nucleosomes. Linker histone H1 binding is also abrogated in nucleosomes containing the incomplete docking domain of H2A.Bbd. Our data illustrate the unique role of the H2A-docking domain in coordinating the structural-functional aspects of the nucleosome properties. Moreover, our data suggest that incorporation of a 'defective' docking domain may be a primary structural role of H2A.Bbd in chromatin.

Neurite sprouting and synapse deterioration in the aging Caenorhabditis elegans nervous system.

Caenorhabditis elegans is a powerful model for analysis of the conserved mechanisms that modulate healthy aging. In the aging nematode nervous system, neuronal death and/or detectable loss of processes are not readily apparent, but because dendrite restructuring and loss of synaptic integrity are hypothesized to contribute to human brain decline and dysfunction, we combined fluorescence microscopy and electron microscopy (EM) to screen at high resolution for nervous system changes. We report two major components of morphological change in the aging C. elegans nervous system: (1) accumulation of novel outgrowths from specific neurons, and (2) physical decline in synaptic integrity. Novel outgrowth phenotypes, including branching from the main dendrite or new growth from somata, appear at a high frequency in some aging neurons, but not all. Mitochondria are often associated with age-associated branch sites. Lowered insulin signaling confers some maintenance of ALM and PLM neuron structural integrity into old age, and both DAF-16/FOXO and heat shock factor transcription factor HSF-1 exert neuroprotective functions. hsf-1 can act cell autonomously in this capacity. EM evaluation in synapse-rich regions reveals a striking decline in synaptic vesicle numbers and a diminution of presynaptic density size. Interestingly, old animals that maintain locomotory prowess exhibit less synaptic decline than same-age decrepit animals, suggesting that synaptic integrity correlates with locomotory healthspan. Our data reveal similarities between the aging C. elegans nervous system and mammalian brain, suggesting conserved neuronal responses to age. Dissection of neuronal aging mechanisms in C. elegans may thus influence the development of brain healthspan-extending therapies.

RecQ helicases: caretakers of the genome

RecQ helicases are highly conserved from bacteria to man. Germline mutations in three of the five known family members in humans give rise to debilitating disorders that are characterized by, amongst other things, a predisposition to the development of cancer. One of these disorders — Bloom's syndrome — is uniquely associated with a predisposition to cancers of all types. So how do RecQ helicases protect against cancer? They seem to maintain genomic stability by functioning at the interface between DNA replication and DNA repair.

Shaping the landscape of the Escherichia coli chromosome: replication-transcription encounters in cells with an ectopic replication origin

Each cell division requires the unwinding of millions of DNA base pairs to allow chromosome duplication and gene transcription. As DNA replication and transcription share the same template, conflicts between both processes are unavoidable and head-on collisions are thought to be particularly problematic. Surprisingly, a recent study reported unperturbed cell cycle progression in Escherichia coli cells with an ectopic replication origin in which highly transcribed rrn operons were forced to be replicated opposite to normal. In this study we have re-generated a similar strain and found the doubling time to be twice that of normal cells. Replication profiles of this background revealed significant deviations in comparison to wild-type profiles, particularly in highly transcribed regions and the termination area. These deviations were alleviated by mutations that either inactivate the termination area or destabilise RNA polymerase complexes and allow their easier displacement by replication forks. Our data demonstrate that head-on replication-transcription conflicts are highly problematic. Indeed, analysis of the replication profile of the previously published E. coli construct revealed a chromosomal rearrangement that alleviates replication-transcription conflicts in an intriguingly simple way. Our data support the idea that avoiding head-on collisions has significantly contributed to shaping the distinct architecture of bacterial chromosomes.

Mfold web server for nucleic acid folding and hybridization prediction

The abbreviated name, 'mfold web server', describes a number of closely related software applications available on the World Wide Web (WWW) for the prediction of the secondary structure of single stranded nucleic acids. The objective of this web server is to provide easy access to RNA and DNA folding and hybridization software to the scientific community at large. By making use of universally available web GUIs (Graphical User Interfaces), the server circumvents the problem of portability of this software. Detailed output, in the form of structure plots with or without reliability information, single strand frequency plots and 'energy dot plots', are available for the folding of single sequences. A variety of 'bulk' servers give less information, but in a shorter time and for up to hundreds of sequences at once. The portal for the mfold web server is http://www.bioinfo.rpi.edu/applications/mfold. This URL will be referred to as 'MFOLDROOT'.

STAR: ultrafast universal RNA-seq aligner.

MOTIVATION Accurate alignment of high-throughput RNA-seq data is a challenging and yet unsolved problem because of the non-contiguous transcript structure, relatively short read lengths and constantly increasing throughput of the sequencing technologies. Currently available RNA-seq aligners suffer from high mapping error rates, low mapping speed, read length limitation and mapping biases. RESULTS To align our large (>80 billon reads) ENCODE Transcriptome RNA-seq dataset, we developed the Spliced Transcripts Alignment to a Reference (STAR) software based on a previously undescribed RNA-seq alignment algorithm that uses sequential maximum mappable seed search in uncompressed suffix arrays followed by seed clustering and stitching procedure. STAR outperforms other aligners by a factor of >50 in mapping speed, aligning to the human genome 550 million 2 × 76 bp paired-end reads per hour on a modest 12-core server, while at the same time improving alignment sensitivity and precision. In addition to unbiased de novo detection of canonical junctions, STAR can discover non-canonical splices and chimeric (fusion) transcripts, and is also capable of mapping full-length RNA sequences. Using Roche 454 sequencing of reverse transcription polymerase chain reaction amplicons, we experimentally validated 1960 novel intergenic splice junctions with an 80-90% success rate, corroborating the high precision of the STAR mapping strategy. AVAILABILITY AND IMPLEMENTATION STAR is implemented as a standalone C++ code. STAR is free open source software distributed under GPLv3 license and can be downloaded from http://code.google.com/p/rna-star/.

Spatial genome organization: contrasting views from chromosome conformation capture and fluorescence in situ hybridization.

Although important for gene regulation, most studies of genome organization use either fluorescence in situ hybridization (FISH) or chromosome conformation capture (3C) methods. FISH directly visualizes the spatial relationship of sequences but is usually applied to a few loci at a time. The frequency at which sequences are ligated together by formaldehyde cross-linking can be measured genome-wide by 3C methods, with higher frequencies thought to reflect shorter distances. FISH and 3C should therefore give the same views of genome organization, but this has not been tested extensively. We investigated the murine HoxD locus with 3C carbon copy (5C) and FISH in different developmental and activity states and in the presence or absence of epigenetic regulators. We identified situations in which the two data sets are concordant but found other conditions under which chromatin topographies extrapolated from 5C or FISH data are not compatible. We suggest that products captured by 3C do not always reflect spatial proximity, with ligation occurring between sequences located hundreds of nanometers apart, influenced by nuclear environment and chromatin composition. We conclude that results obtained at high resolution with either 3C methods or FISH alone must be interpreted with caution and that views about genome organization should be validated by independent methods.

Data analysis methods for detection of differential protein expression in two-dimensional gel electrophoresis.

The recent development of microarray technology has led statisticians and bioinformaticians to develop new statistical methodologies for comparing different biological samples. The objective is to identify a small number of differentially expressed genes from among thousands. In quantitative proteomics, analysis of protein expression using two-dimensional gel electrophoresis shows some similarities with transcriptomic studies. Thus, the goal of this study was to evaluate different data analysis methodologies widely used in array analysis using different proteomic data sets of hundreds of proteins. Even with few replications, the significance analysis of microarrays method appeared to be more powerful than the Student's t test in truly declaring differentially expressed proteins. This procedure will avoid wasting time due to false positives and losing information with false negatives.

Direct-coupling analysis of residue co-evolution captures native contacts across many protein families

The similarity in the three-dimensional structures of homologous proteins imposes strong constraints on their sequence variability. It has long been suggested that the resulting correlations among amino acid compositions at different sequence positions can be exploited to infer spatial contacts within the tertiary protein structure. Crucial to this inference is the ability to disentangle direct and indirect correlations, as accomplished by the recently introduced Direct Coupling Analysis (DCA) (Weigt et al. (2009) Proc Natl Acad Sci 106:67). Here we develop a computationally efficient implementation of DCA, which allows us to evaluate the accuracy of contact prediction by DCA for a large number of protein domains, based purely on sequence information. DCA is shown to yield a large number of correctly predicted contacts, recapitulating the global structure of the contact map for the majority of the protein domains examined. Furthermore, our analysis captures clear signals beyond intra- domain residue contacts, arising, e.g., from alternative protein conformations, ligand- mediated residue couplings, and inter-domain interactions in protein oligomers. Our findings suggest that contacts predicted by DCA can be used as a reliable guide to facilitate computational predictions of alternative protein conformations, protein complex formation, and even the de novo prediction of protein domain structures, provided the existence of a large number of homologous sequences which are being rapidly made available due to advances in genome sequencing.

Multiplex genome engineering using CRISPR/Cas systems.

Functional elucidation of causal genetic variants and elements requires precise genome editing technologies. The type II prokaryotic CRISPR (clustered regularly interspaced short palindromic repeats)/Cas adaptive immune system has been shown to facilitate RNA-guided site-specific DNA cleavage. We engineered two different type II CRISPR/Cas systems and demonstrate that Cas9 nucleases can be directed by short RNAs to induce precise cleavage at endogenous genomic loci in human and mouse cells. Cas9 can also be converted into a nicking enzyme to facilitate homology-directed repair with minimal mutagenic activity. Lastly, multiple guide sequences can be encoded into a single CRISPR array to enable simultaneous editing of several sites within the mammalian genome, demonstrating easy programmability and wide applicability of the RNA-guided nuclease technology.

Improved variant discovery through local re-alignment of short-read next-generation sequencing data using SRMA

A primary component of next-generation sequencing analysis is to align short reads to a reference genome, with each read aligned independently. However, reads that observe the same non-reference DNA sequence are highly correlated and can be used to better model the true variation in the target genome. A novel short-read micro realigner, SRMA, that leverages this correlation to better resolve a consensus of the underlying DNA sequence of the targeted genome is described here.

Global Reorganization of Replication Domains During Embryonic Stem Cell Differentiation

DNA replication in mammals is regulated via the coordinate firing of clusters of replicons that duplicate megabase-sized chromosome segments at specific times during S-phase. Cytogenetic studies show that these "replicon clusters" coalesce as subchromosomal units that persist through multiple cell generations, but the molecular boundaries of such units have remained elusive. Moreover, the extent to which changes in replication timing occur during differentiation and their relationship to transcription changes has not been rigorously investigated. We have constructed high-resolution replication-timing profiles in mouse embryonic stem cells (mESCs) before and after differentiation to neural precursor cells. We demonstrate that chromosomes can be segmented into multimegabase domains of coordinate replication, which we call "replication domains," separated by transition regions whose replication kinetics are consistent with large originless segments. The molecular boundaries of replication domains are remarkably well conserved between distantly related ESC lines and induced pluripotent stem cells. Unexpectedly, ESC differentiation was accompanied by the consolidation of smaller differentially replicating domains into larger coordinately replicated units whose replication time was more aligned to isochore GC content and the density of LINE-1 transposable elements, but not gene density. Replication-timing changes were coordinated with transcription changes for weak promoters more than strong promoters, and were accompanied by rearrangements in subnuclear position. We conclude that replication profiles are cell-type specific, and changes in these profiles reveal chromosome segments that undergo large changes in organization during differentiation. Moreover, smaller replication domains and a higher density of timing transition regions that interrupt isochore replication timing define a novel characteristic of the pluripotent state.

The Sequence Alignment/Map format and SAMtools

SUMMARY The Sequence Alignment/Map (SAM) format is a generic alignment format for storing read alignments against reference sequences, supporting short and long reads (up to 128 Mbp) produced by different sequencing platforms. It is flexible in style, compact in size, efficient in random access and is the format in which alignments from the 1000 Genomes Project are released. SAMtools implements various utilities for post-processing alignments in the SAM format, such as indexing, variant caller and alignment viewer, and thus provides universal tools for processing read alignments. AVAILABILITY http://samtools.sourceforge.net.

H3K9 demethylase KDM4E is an epigenetic regulator for bovine embryonic development and a defective factor for nuclear reprogramming.

Aberrant epigenetic reprogramming often results in developmental defects in somatic cell nuclear transfer (SCNT) embryos during embryonic genome activation (EGA). Bovine eight-cell SCNT embryos exhibit global hypermethylation of histone H3 lysine 9 tri- and di-methylation (H3K9me3/2), but the intrinsic reason for this remains elusive. Here, we provide evidence that two H3K9 demethylase genes, lysine-specific demethylase 4D (KDM4D) and 4E (KDM4E), are related to active H3K9me3/2 demethylation in in vitro fertilized (IVF) embryos and are deficiently expressed in cloned embryos at the time of EGA. Moreover, KDM4E plays a more crucial role in IVF and SCNT embryonic development, and overexpression of KDM4E can restore the global transcriptome, improve blastocyst formation and increase the cloning efficiency of SCNT embryos. Our results thereby indicate that KDM4E can function as a crucial epigenetic regulator of EGA and as an internal defective factor responsible for persistent H3K9me3/2 barriers to SCNT-mediated reprogramming. Furthermore, we show that interactions between RNA and KDM4E are essential for H3K9 demethylation during EGA. These observations advance the understanding of incomplete nuclear reprogramming and are of great importance for transgenic cattle procreation.

Optimized CRISPR/Cas tools for efficient germline and somatic genome engineering in Drosophila.

The type II clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated (Cas) system has emerged recently as a powerful method to manipulate the genomes of various organisms. Here, we report a toolbox for high-efficiency genome engineering of Drosophila melanogaster consisting of transgenic Cas9 lines and versatile guide RNA (gRNA) expression plasmids. Systematic evaluation reveals Cas9 lines with ubiquitous or germ-line-restricted patterns of activity. We also demonstrate differential activity of the same gRNA expressed from different U6 snRNA promoters, with the previously untested U6:3 promoter giving the most potent effect. An appropriate combination of Cas9 and gRNA allows targeting of essential and nonessential genes with transmission rates ranging from 25-100%. We also demonstrate that our optimized CRISPR/Cas tools can be used for offset nicking-based mutagenesis. Furthermore, in combination with oligonucleotide or long double-stranded donor templates, our reagents allow precise genome editing by homology-directed repair with rates that make selection markers unnecessary. Last, we demonstrate a novel application of CRISPR/Cas-mediated technology in revealing loss-of-function phenotypes in somatic cells following efficient biallelic targeting by Cas9 expressed in a ubiquitous or tissue-restricted manner. Our CRISPR/Cas tools will facilitate the rapid evaluation of mutant phenotypes of specific genes and the precise modification of the genome with single-nucleotide precision. Our results also pave the way for high-throughput genetic screening with CRISPR/Cas.

A kinetochore-independent mechanism drives anaphase chromosome separation during acentrosomal meiosis

Although assembly of acentrosomal meiotic spindles has been extensively studied, little is known about the segregation of chromosomes on these spindles. Here, we show in Caenorhabditis elegans oocytes that the kinetochore protein, KNL-1, directs assembly of meiotic kinetochores that orient chromosomes. However, in contrast to mitosis, chromosome separation during meiotic anaphase is kinetochore-independent. Before anaphase, meiotic kinetochores and spindle poles disassemble along with the microtubules on the poleward side of chromosomes. During anaphase, microtubules then form between the separating chromosomes. Functional analysis implicated a set of proteins that localize to a ring-shaped domain between kinetochores during pre-anaphase spindle assembly and anaphase separation. These proteins are localized by the chromosomal passenger complex, which regulates the loss of meiotic chromosome cohesion. Thus, meiotic segregation in C. elegans is a two-stage process, where kinetochores orient chromosomes, but are then dispensable for their separation. We suggest that separation is controlled by a meiosis-specific chromosomal domain to coordinate cohesin removal and chromosome segregation.

Expansion and deletion of CTG repeats from human disease genes are determined by the direction of replication in E. coli

Several human hereditary neurological and neurodegenerative disease genes are associated with the expansion of CTG repeats. Here we show that the frequency of genetic expansions or deletions in Escherichia coli depends on the direction of replication. Large expansions occur predominantly when the CTGs are in the leading strand template rather than the lagging strand. However, deletions are more prominant when the CTGs are in the opposite orientation. Most deletions generated products of defined size classes. Strand slippage coupled with non–classical DMA structures may account for these observations and relate to expansion–deletion mechanisms in eukaryotic chromosomes for disease genes.

Functional elements within the dynein microtubule-binding domain

Dynein interacts with microtubules through an ATP-sensitive linkage mapped to a structurally complex region of the heavy chain following the fourth P-loop motif. Virtually nothing is known regarding how binding affinity is achieved and modulated during ATP hydrolysis. We have performed a detailed dissection of the microtubule contact site, using fragment expression, alanine substitution, and peptide competition. Our work identifies three clusters of amino acids important for the physical contact with microtubules; two of these fall within a region sharing sequence homology with MAP1B, the third in a region just downstream. Amino acid substitutions within any one of these regions can eliminate or weaken microtubule binding (KK3379, 80, E3385, K3387, K3397, KK3410,11, W3414, RKK3418-20, F3426, R3464, S3466, and K3467), suggesting that their activities are highly coordinated. A peptide that actively displaces MAP1B from microtubules perturbs dynein binding, supporting previous evidence for similar sites of interaction. We have also identified four amino acids whose substitutions affect release of the motor from the microtubule (E3413, R3444, E3460, and C3469). These suggest that nucleotide-sensitive affinity may be locally controlled at the site of contact. Our work is the first detailed description of dynein-tubulin interactions and provides a framework for understanding how affinity is achieved and modulated.

CRISPR Interference Limits Horizontal Gene Transfer in Staphylococci by Targeting DNA

Horizontal gene transfer (HGT) in bacteria and archaea occurs through phage transduction, transformation, or conjugation, and the latter is particularly important for the spread of antibiotic resistance. Clustered, regularly interspaced, short palindromic repeat (CRISPR) loci confer sequence-directed immunity against phages. A clinical isolate of Staphylococcus epidermidis harbors a CRISPR spacer that matches the nickase gene present in nearly all staphylococcal conjugative plasmids. Here we show that CRISPR interference prevents conjugation and plasmid transformation in S. epidermidis. Insertion of a self-splicing intron into nickase blocks interference despite the reconstitution of the target sequence in the spliced mRNA, which indicates that the interference machinery targets DNA directly. We conclude that CRISPR loci counteract multiple routes of HGT and can limit the spread of antibiotic resistance in pathogenic bacteria.

A key role for Ctf4 in coupling the MCM2-7 helicase to DNA polymerase alpha within the eukaryotic replisome.

The eukaryotic replisome is a crucial determinant of genome stability, but its structure is still poorly understood. We found previously that many regulatory proteins assemble around the MCM2-7 helicase at yeast replication forks to form the replisome progression complex (RPC), which might link MCM2-7 to other replisome components. Here, we show that the RPC associates with DNA polymerase alpha that primes each Okazaki fragment during lagging strand synthesis. Our data indicate that a complex of the GINS and Ctf4 components of the RPC is crucial to couple MCM2-7 to DNA polymerase alpha. Others have found recently that the Mrc1 subunit of RPCs binds DNA polymerase epsilon, which synthesises the leading strand at DNA replication forks. We show that cells lacking both Ctf4 and Mrc1 experience chronic activation of the DNA damage checkpoint during chromosome replication and do not complete the cell cycle. These findings indicate that coupling MCM2-7 to replicative polymerases is an important feature of the regulation of chromosome replication in eukaryotes, and highlight a key role for Ctf4 in this process.

Sequential domain assembly of ribosomal protein S3 drives 40S subunit maturation.

Eukaryotic ribosomes assemble by association of ribosomal RNA with ribosomal proteins into nuclear precursor particles, which undergo a complex maturation pathway coordinated by non-ribosomal assembly factors. Here, we provide functional insights into how successive structural re-arrangements in ribosomal protein S3 promote maturation of the 40S ribosomal subunit. We show that S3 dimerizes and is imported into the nucleus with its N-domain in a rotated conformation and associated with the chaperone Yar1. Initial assembly of S3 with 40S precursors occurs via its C-domain, while the N-domain protrudes from the 40S surface. Yar1 is replaced by the assembly factor Ltv1, thereby fixing the S3 N-domain in the rotated orientation and preventing its 40S association. Finally, Ltv1 release, triggered by phosphorylation, and flipping of the S3 N-domain into its final position results in the stable integration of S3. Such a stepwise assembly may represent a new paradigm for the incorporation of ribosomal proteins.

Combinatorial control of exon recognition.

Pre-mRNA splicing is a fundamental process required for the expression of most metazoan genes. It is carried out by the spliceosome, which catalyzes the removal of noncoding intronic sequences to assemble exons into mature mRNAs prior to export and translation. Given the complexity of higher eukaryotic genes and the relatively low level of splice site conservation, the precision of the splicing machinery in recognizing and pairing splice sites is impressive. Introns ranging in size from <100 up to 100,000 bases are removed efficiently. At the same time, a large number of alternative splicing events are observed between different cell types, during development, or during other biological processes. This extensive alternative splicing implies a significant flexibility of the spliceosome to identify and process exons within a given pre-mRNA. To reach this flexibility, splice site selection in higher eukaryotes has evolved to depend on multiple parameters such as splice site strength, the presence or absence of splicing regulators, RNA secondary structures, the exon/intron architecture, and the process of pre-mRNA synthesis itself. The relative contributions of each of these parameters control how efficiently splice sites are recognized and flanking introns are removed.

Direct isolation and identification of promoters in the human genome.

Transcriptional regulatory elements play essential roles in gene expression during animal development and cellular response to environmental signals, but our knowledge of these regions in the human genome is limited despite the availability of the complete genome sequence. Promoters mark the start of every transcript and are an important class of regulatory elements. A large, complex protein structure known as the pre-initiation complex (PIC) is assembled on all active promoters, and the presence of these proteins distinguishes promoters from other sequences in the genome. Using components of the PIC as tags, we isolated promoters directly from human cells as protein-DNA complexes and identified the resulting DNA sequences using genomic tiling microarrays. Our experiments in four human cell lines uncovered 252 PIC-binding sites in 44 semirandomly selected human genomic regions comprising 1% (30 megabase pairs) of the human genome. Nearly 72% of the identified fragments overlap or immediately flank 5' ends of known cDNA sequences, while the remainder is found in other genomic regions that likely harbor putative promoters of unannotated transcripts. Indeed, molecular analysis of the RNA isolated from one cell line uncovered transcripts initiated from over half of the putative promoter fragments, and transient transfection assays revealed promoter activity for a significant proportion of fragments when they were fused to a luciferase reporter gene. These results demonstrate the specificity of a genome-wide analysis method for mapping transcriptional regulatory elements and also indicate that a small, yet significant number of human genes remains to be discovered.

Genetic and physical analysis of double-strand break repair and recombination in Saccharomyces cerevisiae.

We have investigated HO endonuclease-induced double-strand break (DSB) recombination and repair in a LACZ duplication plasmid in yeast. A 117-bp MATa fragment, embedded in one copy of LACZ, served as a site for initiation of a DSB when HO endonuclease was expressed. The DSB could be repaired using wild-type sequences located on a second, promoterless, copy of LACZ on the same plasmid. In contrast to normal mating-type switching, crossing-over associated with gene conversion occurred at least 50% of the time. The proportion of conversion events accompanied by exchange was greater when the two copies of LACZ were in direct orientation (80%), than when inverted (50%). In addition, the fraction of plasmids lost was significantly greater in the inverted orientation. The kinetics of appearance of intermediates and final products were also monitored. The repair of the DSB is slow, requiring at least an hour from the detection of the HO-cut fragments to completion of repair. Surprisingly, the appearance of the two reciprocal products of crossing over did not occur with the same kinetics. For example, when the two LACZ sequences were in the direct orientation, the HO-induced formation of a large circular deletion product was not accompanied by the appearance of a small circular reciprocal product. We suggest that these differences may reflect two kinetically separable processes, one involving only one cut end and the other resulting from the concerted participation of both ends of the DSB.

ModLoop: automated modeling of loops in protein structures.

SUMMARY ModLoop is a web server for automated modeling of loops in protein structures. The input is the atomic coordinates of the protein structure in the Protein Data Bank format, and the specification of the starting and ending residues of one or more segments to be modeled, containing no more than 20 residues in total. The output is the coordinates of the non-hydrogen atoms in the modeled segments. A user provides the input to the server via a simple web interface, and receives the output by e-mail. The server relies on the loop modeling routine in MODELLER that predicts the loop conformations by satisfaction of spatial restraints, without relying on a database of known protein structures. For a rapid response, ModLoop runs on a cluster of Linux PC computers. AVAILABILITY The server is freely accessible to academic users at http://salilab.org/modloop

hnRNP A1 associates with telomere ends and stimulates telomerase activity.

Telomerase is a ribonucleoprotein enzyme complex that reverse-transcribes an integral RNA template to add short DNA repeats to the 3'-ends of telomeres. G-quadruplex structure in a DNA substrate can block its extension by telomerase. We have found that hnRNP A1--which was previously implicated in telomere length regulation--binds to both single-stranded and structured human telomeric repeats, and in the latter case, it disrupts their higher-order structure. Using an in vitro telomerase assay, we observed that depletion of hnRNP A/B proteins from 293 human embryonic kidney cell extracts dramatically reduced telomerase activity, which was fully recovered upon addition of purified recombinant hnRNP A1. This finding suggests that hnRNP A1 functions as an auxiliary, if not essential, factor of telomerase holoenzyme. We further show, using chromatin immunoprecipitation, that hnRNP A1 associates with human telomeres in vivo. We propose that hnRNP A1 stimulates telomere elongation through unwinding of a G-quadruplex or G-G hairpin structure formed at each translocation step.

Crystal structure of a nucleosome core particle containing the variant histone H2A.Z

Activation of transcription within chromatin has been correlated with the incorporation of the essential histone variant H2A.Z into nucleosomes. H2A.Z and other histone variants may establish structurally distinct chromosomal domains; however, the molecular mechanism by which they function is largely unknown. Here we report the 2.6 Å crystal structure of a nucleosome core particle containing the histone variant H2A.Z. The overall structure is similar to that of the previously reported 2.8 Å nucleosome structure containing major histone proteins. However, distinct localized changes result in the subtle destabilization of the interaction between the (H2A.Z–H2B) dimer and the (H3–H4)2 tetramer. Moreover, H2A.Z nucleosomes have an altered surface that includes a metal ion. This altered surface may lead to changes in higher order structure, and/or could result in the association of specific nuclear proteins with H2A.Z. Finally, incorporation of H2A.Z and H2A within the same nucleosome is unlikely, due to significant changes in the interface between the two H2A.Z–H2B dimers.

US-SOMO HPLC-SAXS module: dealing with capillary fouling and extraction of pure component patterns from poorly resolved SEC-SAXS data

Size-exclusion chromatography coupled with SAXS (small-angle X-ray scattering), often performed using a flow-through capillary, should allow direct collection of monodisperse sample data. However, capillary fouling issues and non-baseline-resolved peaks can hamper its efficacy. The UltraScan solution modeler (US-SOMO) HPLC-SAXS (high-performance liquid chromatography coupled with SAXS) module provides a comprehensive framework to analyze such data, starting with a simple linear baseline correction and symmetrical Gaussian decomposition tools [Brookes, Pérez, Cardinali, Profumo, Vachette & Rocco (2013 ▸). J. Appl. Cryst.46, 1823-1833]. In addition to several new features, substantial improvements to both routines have now been implemented, comprising the evaluation of outcomes by advanced statistical tools. The novel integral baseline-correction procedure is based on the more sound assumption that the effect of capillary fouling on scattering increases monotonically with the intensity scattered by the material within the X-ray beam. Overlapping peaks, often skewed because of sample interaction with the column matrix, can now be accurately decomposed using non-symmetrical modified Gaussian functions. As an example, the case of a polydisperse solution of aldolase is analyzed: from heavily convoluted peaks, individual SAXS profiles of tetramers, octamers and dodecamers are extracted and reliably modeled.

Acentrosomal Microtubule Assembly in Mitosis: The Where, When, and How.

In mitosis the cell assembles the bipolar spindle, a microtubule (MT)-based apparatus that segregates the duplicated chromosomes into two daughter cells. Most animal cells enter mitosis with duplicated centrosomes that provide an active source of dynamic MTs. However, it is now established that spindle assembly relies on the nucleation of acentrosomal MTs occurring around the chromosomes after nuclear envelope breakdown, and on pre-existing microtubules. Where chromosome-dependent MT nucleation occurs, when MT amplification takes place and how the two pathways function are still key questions that generate some controversies. We reconcile the data and present an integrated model accounting for acentrosomal microtubule assembly in the dividing cell.

Efficient Generation of Fully Reprogrammed Human iPS Cells via Polycistronic Retroviral Vector and a New Cocktail of Chemical Compounds

Direct reprogramming of human somatic cells into induced pluripotent stem (iPS) cells by defined transcription factors (TFs) provides great potential for regenerative medicine and biomedical research. This procedure has many challenges, including low reprogramming efficiency, many partially reprogrammed colonies, somatic coding mutations in the genome, etc. Here, we describe a simple approach for generating fully reprogrammed human iPS cells by using a single polycistronic retroviral vector expressing four human TFs in a single open reading frame (ORF), combined with a cocktail containing three small molecules (Sodium butyrate, SB431542, and PD0325901). Our results demonstrate that human iPS cells generated by this approach express human ES cells markers and exhibit pluripotency demonstrated by their abilities to differentiate into the three germ layers in vitro and in vivo. Notably, this approach not only provides a much faster reprogramming process but also significantly diminishes partially reprogrammed iPS cell colonies, thus facilitating efficient isolation of desired fully reprogrammed iPS cell colonies.

The Protein Data Bank

The Protein Data Bank [PDB; Berman, Westbrook et al. (2000), Nucleic Acids Res. 28, 235-242; http://www.pdb.org/] is the single worldwide archive of primary structural data of biological macromolecules. Many secondary sources of information are derived from PDB data. It is the starting point for studies in structural bioinformatics. This article describes the goals of the PDB, the systems in place for data deposition and access, how to obtain further information and plans for the future development of the resource. The reader should come away with an understanding of the scope of the PDB and what is provided by the resource.

GSDS 2.0: an upgraded gene feature visualization server

UNLABELLED : Visualizing genes' structure and annotated features helps biologists to investigate their function and evolution intuitively. The Gene Structure Display Server (GSDS) has been widely used by more than 60 000 users since its first publication in 2007. Here, we reported the upgraded GSDS 2.0 with a newly designed interface, supports for more types of annotation features and formats, as well as an integrated visual editor for editing the generated figure. Moreover, a user-specified phylogenetic tree can be added to facilitate further evolutionary analysis. The full source code is also available for downloading. AVAILABILITY AND IMPLEMENTATION Web server and source code are freely available at http://gsds.cbi.pku.edu.cn. CONTACT gaog@mail.cbi.pku.edu.cn or gsds@mail.cbi.pku.edu.cn   SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.

CLIP-170 interacts with dynactin complex and the APC-binding protein EB1 by different mechanisms.

CLIP-170 is a "cytoplasmic linker protein" implicated in endosome-microtubule interactions and in control of microtubule dynamics. CLIP-170 localizes dynamically to growing microtubule plus ends, colocalizing with the dynein activator dynactin and the APC-binding protein EB1. This shared "plus-end tracking" behavior suggests that CLIP-170 might interact with dynactin and/or EB1. We have used site-specific mutagenesis of CLIP-170 and a transfection/colocalization assay to address this question in mammalian tissue culture cells. Our results indicate that CLIP-170 interacts, directly or indirectly, with both dynactin and EB1. We find that the CLIP-170/dynactin interaction is mediated by the second metal binding motif of the CLIP-170 tail. In contrast, the CLIP-170/EB1 interaction requires neither metal binding motif. In addition, our experiments suggest that the CLIP-170/dynactin interaction occurs via the shoulder/sidearm subcomplex of dynactin and can occur in the cytosol (i.e., it does not require microtubule binding). These results have implications for the targeting of both dynactin and EB1 to microtubule plus ends. Our data suggest that the CLIP-170/dynactin interaction can target dynactin complex to microtubule plus ends, although dynactin likely also targets MT plus ends directly via the microtubule binding motif of the p150(Glued) subunit. We find that CLIP-170 mutants alter p150(Glued) localization without affecting EB1, indicating that EB1 can target microtubule plus ends independently of dynactin.

A framework for improving microRNA prediction in non-human genomes

The prediction of novel pre-microRNA (miRNA) from genomic sequence has received considerable attention recently. However, the majority of studies have focused on the human genome. Previous studies have demonstrated that sensitivity (correctly detecting true miRNA) is sustained when human-trained methods are applied to other species, however they have failed to report the dramatic drop in specificity (the ability to correctly reject non-miRNA sequences) in non-human genomes. Considering the ratio of true miRNA sequences to pseudo-miRNA sequences is on the order of 1:1000, such low specificity prevents the application of most existing tools to non-human genomes, as the number of false positives overwhelms the true predictions. We here introduce a framework (SMIRP) for creating species-specific miRNA prediction systems, leveraging sequence conservation and phylogenetic distance information. Substantial improvements in specificity and precision are obtained for four non-human test species when our framework is applied to three different prediction systems representing two types of classifiers (support vector machine and Random Forest), based on three different feature sets, with both human-specific and taxon-wide training data. The SMIRP framework is potentially applicable to all miRNA prediction systems and we expect substantial improvement in precision and specificity, while sustaining sensitivity, independent of the machine learning technique chosen.

Acetylation of Histone H3 Lysine 56 Regulates Replication-Coupled Nucleosome Assembly

Chromatin assembly factor 1 (CAF-1) and Rtt106 participate in the deposition of newly synthesized histones onto replicating DNA to form nucleosomes. This process is critical for the maintenance of genome stability and inheritance of functionally specialized chromatin structures in proliferating cells. However, the molecular functions of the acetylation of newly synthesized histones in this DNA replication-coupled nucleosome assembly pathway remain enigmatic. Here we show that histone H3 acetylated at lysine 56 (H3K56Ac) is incorporated onto replicating DNA and, by increasing the binding affinity of CAF-1 and Rtt106 for histone H3, H3K56Ac enhances the ability of these histone chaperones to assemble DNA into nucleosomes. Genetic analysis indicates that H3K56Ac acts in a nonredundant manner with the acetylation of the N-terminal residues of H3 and H4 in nucleosome assembly. These results reveal a mechanism by which H3K56Ac regulates replication-coupled nucleosome assembly mediated by CAF-1 and Rtt106.

Human artificial chromosomes containing chromosome 17 alphoid DNA maintain an active centromere in murine cells but are not stable.

Human artificial chromosomes (HACs) are autonomous molecules that can function and segregate as normal chromosomes in human cells. De novo HACs have successfully been used as gene expression vectors to complement genetic deficiencies in human cultured cells. HACs now offer the possibility of studying the regulation and expression of large genes in a variety of cell types from different tissues and correcting gene deficiencies caused by human inherited diseases. Complementary gene expression studies in mice, especially in mouse models of human genetic diseases, are also important in determining if large human transgenes can be expressed appropriately from artificial chromosomes. Toward this aim we are establishing artificial chromosomes in murine cells as novel gene expression vectors. Initially we transferred HAC vectors into murine cells, but were unable to generate de novo HACs at a reasonable frequency. We then transferred HACs previously established in human HT1080 cells to three different murine cell types by microcell fusion, followed by positive selection. We observed that the HACs in murine cells bound centromere protein C (CENP-C), a marker of active centromeres, and were detected under selection but rapidly lost when selection was removed. These results suggest that the HACs maintain at least a partially functional centromere complex in murine cells, but other factors are required for stability and segregation. Artificial chromosomes containing mouse centromeric sequences may be required for better stability and maintenance in murine cells.

The nucleosome family: dynamic and growing.

Ever since the discovery of the nucleosome in 1974, scientists have stumbled upon discrete particles in which DNA is wrapped around histone complexes of different stoichiometries: octasomes, hexasomes, tetrasomes, "split" half-nucleosomes, and, recently, bona fide hemisomes. Do all these particles exist in vivo? Under what conditions? What is their physiological significance in the complex DNA transactions in the eukaryotic nucleus? What are their dynamics? This review summarizes research spanning more than three decades and provides a new meaning to the term "nucleosome. " The nucleosome can no longer be viewed as a single static entity: rather, it is a family of particles differing in their structural and dynamic properties, leading to different functionalities.

The structure of DdrB from Deinococcus: a new fold for single-stranded DNA binding proteins

Deinococcus spp. are renowned for their amazing ability to recover rapidly from severe genomic fragmentation as a result of exposure to extreme levels of ionizing radiation or desiccation. Despite having been originally characterized over 50 years ago, the mechanism underlying this remarkable repair process is still poorly understood. Here, we report the 2.8 A structure of DdrB, a single-stranded DNA (ssDNA) binding protein unique to Deinococcus spp. that is crucial for recovery following DNA damage. DdrB forms a pentameric ring capable of binding single-stranded but not double-stranded DNA. Unexpectedly, the crystal structure reveals that DdrB comprises a novel fold that is structurally and topologically distinct from all other single-stranded binding (SSB) proteins characterized to date. The need for a unique ssDNA binding function in response to severe damage, suggests a distinct role for DdrB which may encompass not only standard SSB protein function in protection of ssDNA, but also more specialized roles in protein recruitment or DNA architecture maintenance. Possible mechanisms of DdrB action in damage recovery are discussed.

Circos: an information aesthetic for comparative genomics.

We created a visualization tool called Circos to facilitate the identification and analysis of similarities and differences arising from comparisons of genomes. Our tool is effective in displaying variation in genome structure and, generally, any other kind of positional relationships between genomic intervals. Such data are routinely produced by sequence alignments, hybridization arrays, genome mapping, and genotyping studies. Circos uses a circular ideogram layout to facilitate the display of relationships between pairs of positions by the use of ribbons, which encode the position, size, and orientation of related genomic elements. Circos is capable of displaying data as scatter, line, and histogram plots, heat maps, tiles, connectors, and text. Bitmap or vector images can be created from GFF-style data inputs and hierarchical configuration files, which can be easily generated by automated tools, making Circos suitable for rapid deployment in data analysis and reporting pipelines.

In vitro splicing pathways of pre-mRNAs containing multiple intervening sequences?

We analyzed the in vitro splicing pathways of three multi-intervening-sequence (IVS) pre-mRNAs: human beta-globin, which contains two IVSs (K. M. Lang, V. L. van Santen, and R. A. Spritz, EMBO J. 4:1991-1996, 1985); rat alpha-lactalbumin, which contains three IVSs; and murine interleukin-3, which contains four IVSs. We found that there are highly preferred pathways of IVS removal from these multi-IVS pre-mRNAs in vitro. The three IVSs of rat alpha-lactalbumin pre-mRNA were excised sequentially from 5' to 3'; in most molecules, IVS1 was removed first, followed by IVS2 and finally by IVS3. The splicing pathway of interleukin-3 pre-mRNA in vitro was more complex. The four IVSs were excised in a highly preferred temporal order, but the order was not strictly sequential or directional. In most molecules, IVS1 and IVS4 were removed first, either simultaneously or in rapid succession. Subsequently, IVS2 was excised, followed by IVS3. The observed splicing pathways apparently resulted from differences in lag times and maximum excision rates of the different IVSs. We detected no exon skipping during splicing of these transcripts in vitro. These observations have implication for proposed models of splice site selection.

Rapid and Sensitive RT-QuIC Detection of Human Creutzfeldt-Jakob Disease Using Cerebrospinal Fluid

UNLABELLED Fast, definitive diagnosis of Creutzfeldt-Jakob disease (CJD) is important in assessing patient care options and transmission risks. Real-time quaking-induced conversion (RT-QuIC) assays of cerebrospinal fluid (CSF) and nasal-brushing specimens are valuable in distinguishing CJD from non-CJD conditions but have required 2.5 to 5 days. Here, an improved RT-QuIC assay is described which identified positive CSF samples within 4 to 14 h with better analytical sensitivity. Moreover, analysis of 11 CJD patients demonstrated that while 7 were RT-QuIC positive using the previous conditions, 10 were positive using the new assay. In these and further analyses, a total of 46 of 48 CSF samples from sporadic CJD patients were positive, while all 39 non-CJD patients were negative, giving 95.8% diagnostic sensitivity and 100% specificity. This second-generation RT-QuIC assay markedly improved the speed and sensitivity of detecting prion seeds in CSF specimens from CJD patients. This should enhance prospects for rapid and accurate ante mortem CJD diagnosis. IMPORTANCE A long-standing problem in dealing with various neurodegenerative protein misfolding diseases is early and accurate diagnosis. This issue is particularly important with human prion diseases, such as CJD, because prions are deadly, transmissible, and unusually resistant to decontamination. The recently developed RT-QuIC test allows for highly sensitive and specific detection of CJD in human cerebrospinal fluid and is being broadly implemented as a key diagnostic tool. However, as currently applied, RT-QuIC takes 2.5 to 5 days and misses 11 to 23% of CJD cases. Now, we have markedly improved RT-QuIC analysis of human CSF such that CJD and non-CJD patients can be discriminated in a matter of hours rather than days with enhanced sensitivity. These improvements should allow for much faster, more accurate, and practical testing for CJD. In broader terms, our study provides a prototype for tests for misfolded protein aggregates that cause many important amyloid diseases, such as Alzheimer's, Parkinson's, and tauopathies.

Reprogramming factor expression initiates widespread targeted chromatin remodeling.

Despite rapid progress in characterizing transcription factor-driven reprogramming of somatic cells to an induced pluripotent stem cell (iPSC) state, many mechanistic questions still remain. To gain insight into the earliest events in the reprogramming process, we systematically analyzed the transcriptional and epigenetic changes that occur during early factor induction after discrete numbers of divisions. We observed rapid, genome-wide changes in the euchromatic histone modification, H3K4me2, at more than a thousand loci including large subsets of pluripotency-related or developmentally regulated gene promoters and enhancers. In contrast, patterns of the repressive H3K27me3 modification remained largely unchanged except for focused depletion specifically at positions where H3K4 methylation is gained. These chromatin regulatory events precede transcriptional changes within the corresponding loci. Our data provide evidence for an early, organized, and population-wide epigenetic response to ectopic reprogramming factors that clarify the temporal order through which somatic identity is reset during reprogramming.

The evolution of RNAs with multiple functions.

Increasing numbers of transcripts have been reported to transmit both protein-coding and regulatory information. Apart from challenging our conception of the gene, this observation raises the question as to what extent this phenomenon occurs across the genome and how and why such dual encoding of function has evolved in the eukaryotic genome. To address this question, we consider the evolutionary path of genes in the earliest forms of life on Earth, where it is generally regarded that proteins evolved from a cellular machinery based entirely within RNA. This led to the domination of protein-coding genes in the genomes of microorganisms, although it is likely that RNA never lost its other capacities and functionalities, as evidenced by cis-acting riboswitches and UTRs. On the basis that the subsequent evolution of a more sophisticated regulatory architecture to provide higher levels of epigenetic control and accurate spatiotemporal expression in developmentally complex organisms is a complicated task, we hypothesize: (i) that mRNAs have been and remain subject to secondary selection to provide trans-acting regulatory capability in parallel with protein-coding functions; (ii) that some and perhaps many protein-coding loci, possibly as a consequence of gene duplication, have lost protein-coding functions en route to acquiring more sophisticated trans-regulatory functions; (iii) that many transcripts have become subject to secondary processing to release different products; and (iv) that novel proteins have emerged within loci that previously evolved functionality as regulatory RNAs. In support of the idea that there is a dynamic flux between different types of informational RNAs in both evolutionary and real time, we review recent observations that have arisen from transcriptomic surveys of complex eukaryotes and reconsider how these observations impact on the notion that apparently discrete loci may express transcripts with more than one function. In conclusion, we posit that many eukaryotic loci have evolved the capacity to transact a multitude of overlapping and potentially independent functions as both regulatory and protein-coding RNAs.

Hotspots of aberrant epigenomic reprogramming in human induced pluripotent stem cells

Induced pluripotent stem cells (iPSCs) offer immense potential for regenerative medicine and studies of disease and development. Somatic cell reprogramming involves epigenomic reconfiguration, conferring iPSCs with characteristics similar to embryonic stem (ES) cells. However, it remains unknown how complete the reestablishment of ES-cell-like DNA methylation patterns is throughout the genome. Here we report the first whole-genome profiles of DNA methylation at single-base resolution in five human iPSC lines, along with methylomes of ES cells, somatic cells, and differentiated iPSCs and ES cells. iPSCs show significant reprogramming variability, including somatic memory and aberrant reprogramming of DNA methylation. iPSCs share megabase-scale differentially methylated regions proximal to centromeres and telomeres that display incomplete reprogramming of non-CG methylation, and differences in CG methylation and histone modifications. Lastly, differentiation of iPSCs into trophoblast cells revealed that errors in reprogramming CG methylation are transmitted at a high frequency, providing an iPSC reprogramming signature that is maintained after differentiation.

Generation and characterization of methyl-lysine histone antibodies.

Publisher Summary This chapter summarizes the development and characterization of rabbit polyclonal antibodies named histone that are directed against the methylated H3-K9 position. It provides protocols for peptide design, rabbit immunizations, and quality controls of methyl-lysine histone antibodies, followed by their in vivo characterization using indirect IF of inter-and metaphase chromatin in wild-type (wt) and mutant mouse cells that are deficient for the Suv39h histone methyltransferases (HMTases). Histone amino-termini (tails) protrude from the nucleosome core and are subject to a variety of post-translational modifications, including acetylation (on lysine residues), phosphorylation (on serine and threonine residues), methylation (on lysine and arginine residues), ubiquitination (on lysine residues), and ADP-ribosylation (on glutamic acid residues). In addition to their structural roles, histones play important functions in the control of gene expression by regulating access to the underlying nucleosomal template. It is without doubt that the development of high-quality, position-specific methyl-lysine histone antibodies can provide important tools for the further decoding of the epigenetic information, which is in part, indexed by distinct methylation states of selective lysine residues in the histone amino-termini. A comparative analysis indicates significant discrepancies in the specificity and avidity of the available methyl-lysine histone antibodies and highlights the need for extensive quality controls, such that experimental data can be correctly interpreted despite the exquisite complexity of histone lysine methylation.

Nucleosome stability mediated by histone variants H3.3 and H2A.Z.

Nucleosomes containing the histone variant H3.3 tend to be clustered in vivo in the neighborhood of transcriptionally active genes and over regulatory elements. It has not been clear, however, whether H3.3-containing nucleosomes possess unique properties that would affect transcription. We report here that H3.3 nucleosomes isolated from vertebrates, regardless of whether they are partnered with H2A or H2A.Z, are unusually sensitive to salt-dependent disruption, losing H2A/H2B or H2A.Z/H2B dimers. Immunoprecipitation studies of nucleosome core particles (NCPs) show that NCPs that contain both H3.3 and H2A.Z are even less stable than NCPs containing H3.3 and H2A. Intriguingly, NCPs containing H3 and H2A.Z are at least as stable as H3/H2A NCPs. These results establish an hierarchy of stabilities for native nucleosomes carrying different complements of variants, and suggest how H2A.Z could play different roles depending on its partners within the NCP. They also are consistent with the idea that H3.3 plays an active role in maintaining accessible chromatin structures in enhancer regions and transcribed regions. Consistent with this idea, promoters and enhancers at transcriptionally active genes and coding regions at highly expressed genes have nucleosomes that simultaneously carry both H3.3 and H2A.Z, and should therefore be extremely sensitive to disruption.

A Role for NuSAP in Linking Microtubules to Mitotic Chromosomes

The spindle apparatus is a microtubule (MT)-based machinery that attaches to and segregates the chromosomes during mitosis and meiosis. Self-organization of the spindle around chromatin involves the assembly of MTs, their attachment to the chromosomes, and their organization into a bipolar array. One regulator of spindle self-organization is RanGTP. RanGTP is generated at chromatin and activates a set of soluble, Ran-regulated spindle factors such as TPX2, NuMA, and NuSAP . How the spindle factors direct and attach MTs to the chromosomes are key open questions. Nucleolar and Spindle-Associated Protein (NuSAP) was recently identified as an essential MT-stabilizing and bundling protein that is enriched at the central part of the spindle . Here, we show by biochemical reconstitution that NuSAP efficiently adsorbs to isolated chromatin and DNA and that it can directly produce and retain high concentrations of MTs in the immediate vicinity of chromatin or DNA. Moreover, our data reveal that NuSAP-chromatin interaction is subject to Ran regulation and can be suppressed by Importin alpha (Impalpha) and Imp7. We propose that the presence of MT binding agents such as NuSAP, which can be directly immobilized on chromatin, are critical for targeting MT production to vertebrate chromosomes during spindle self-organization.

Genomic Nucleosome Organization Reconstituted with Pure Proteins

Chromatin remodelers regulate genes by organizing nucleosomes around promoters, but their individual contributions are obfuscated by the complex in vivo milieu of factor redundancy and indirect effects. Genome-wide reconstitution of promoter nucleosome organization with purified proteins resolves this problem and is therefore a critical goal. Here, we reconstitute four stages of nucleosome architecture using purified components: yeast genomic DNA, histones, sequence-specific Abf1/Reb1, and remodelers RSC, ISW2, INO80, and ISW1a. We identify direct, specific, and sufficient contributions that in vivo observations validate. First, RSC clears promoters by translating poly(dA:dT) into directional nucleosome removal. Second, partial redundancy is recapitulated where INO80 alone, or ISW2 at Abf1/Reb1sites, positions +1 nucleosomes. Third, INO80 and ISW2 each align downstream nucleosomal arrays. Fourth, ISW1a tightens the spacing to canonical repeat lengths. Such a minimal set of rules and proteins establishes core mechanisms by which promoter chromatin architecture arises through a blend of redundancy and specialization.

Osteoclast nuclei of myeloma patients show chromosome translocations specific for the myeloma cell clone: a new type of cancer-host partnership?

A major clinical manifestation of bone cancers is bone destruction. It is widely accepted that this destruction is not caused by the malignant cells themselves, but by osteoclasts, multinucleated cells of monocytic origin that are considered to be the only cells able to degrade bone. The present study demonstrates that bone-resorbing osteoclasts from myeloma patients contain nuclei with translocated chromosomes of myeloma B-cell clone origin, in addition to nuclei without these translocations, by using combined FISH and immunohistochemistry on bone sections. These nuclei of malignant origin are transcriptionally active and appear fully integrated amongst the other nuclei. The contribution of malignant nuclei to the osteoclast population analysed in this study was greater than 30%. Osteoclast-myeloma clone hybrids contained more nuclei than normal osteoclasts and their occurrence correlated with the proximity of myeloma cells. Similar hybrid cells were generated in myeloma cell-osteoclast co-cultures, as revealed by tracing myeloma nuclei using translocations, bromo-deoxyuridine, or the Y chromosome of male myeloma cells in female osteoclasts. These observations indicate that hybrid cells can originate through fusion between myeloma cells and osteoclasts. In conclusion, malignant cells contribute significantly to the formation of bone-resorbing osteoclasts in multiple myeloma. Osteoclast-myeloma clone hybrids reflect a previously unrecognized mechanism of bone destruction in which malignant cells participate directly. The possibility that malignant cells corrupt host cells by the transfer of malignant DNA may have been underestimated to date in cancer research.

Consensus structural features of purified bacterial TatABC complexes.

The twin-arginine translocation (Tat) system transports folded proteins across bacterial plasma membranes and the chloroplast thylakoid membrane. Here, we investigate the composition and structural organization of three different purified Tat complexes from Escherichia coli, Salmonella typhimurium and Agrobacterium tumefaciens. First, we demonstrate the functional activity of these Tat systems in vivo, since expression of the tatABC operons from S.typhimurium or A.tumefaciens in an E.coli tat null mutant strain resulted in efficient Tat-dependent export of an E.coli cofactor-containing substrate, TMAO reductase. The three isolated, affinity-tagged Tat complexes comprised TatA, TatB and TatC in each case, demonstrating a strong interaction between these three subunits. Single-particle electron microscopy studies of all three complexes revealed approximately oval-shaped, asymmetric particles with maximal dimensions up to 13 nm. A common feature is a number of stain-excluding densities surrounding more or less central pools of stain, suggesting protein-lined pores or cavities. The characteristics of size variation among the particles suggest a modular form of assembly and/or the recruitment of varying numbers of TatBC/TatA units. Despite low levels of sequence homology, the combined data indicate structural and functional conservation in the Tat systems of these three bacterial species.

Sperm cross-over activity in regions of the human genome showing extreme breakdown of marker association.

Population diversity data have recently provided profound, albeit inferential, insights into meiotic recombination across the human genome, revealing a landscape dominated by thousands of cross-over hotspots. However, very few of these putative hotspots have been directly analyzed for cross-over activity. We now describe a search for very active hotspots, by using extreme breakdown of marker association as a guide for high-resolution sperm cross-over analysis. This strategy has led to the isolation of the most active cross-over hotspots yet described. Their morphology, sequence attributes, and cross-over processes are very similar to those seen at less active hotspots, but their activity in sperm is poorly predicted from population diversity information. Several of these hotspots showed evidence for biased gene conversion accompanying cross-over, in some cases associated with variation between men in cross-over activity and with two hotspots showing complete presence/absence polymorphism in different men. Hotspot polymorphism is very common at less active hotspots but curiously was not seen at any of the most active hotspots. This contrasts with the prediction that extreme hotspots should be the most vulnerable to attenuation by meiotic drive in favor of mutations that suppress recombination and should therefore show rapid rate evolution and thus variation in activity between men. Finally, these very intense hotspots provide a valuable resource for dissecting meiotic recombination processes and pathways in humans.