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

Title: Full-length transcriptome assembly from RNA-Seq data without a reference genome. Content: 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.

Title: Effects of GC Bias in Next-Generation-Sequencing Data on De Novo Genome Assembly Content: 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.

Title: RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome Content: 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.

Title: Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Content: High-throughput mRNA sequencing (RNA-Seq) promises simultaneous transcript discovery and abundance estimation. However, this would require algorithms that are not restricted by prior gene annotations and that account for alternative transcription and splicing. Here we introduce such algorithms in an open-source software program called Cufflinks. To test Cufflinks, we sequenced and analyzed >430 million paired 75-bp RNA-Seq reads from a mouse myoblast cell line over a differentiation time series. We detected 13,692 known transcripts and 3,724 previously unannotated ones, 62% of which are supported by independent expression data or by homologous genes in other species. Over the time series, 330 genes showed complete switches in the dominant transcription start site (TSS) or splice isoform, and we observed more subtle shifts in 1,304 other genes. These results suggest that Cufflinks can illuminate the substantial regulatory flexibility and complexity in even this well-studied model of muscle development and that it can improve transcriptome-based genome annotation.

Title: De novo assembly of a PML nuclear subcompartment occurs through multiple pathways and induces telomere elongation. Content: 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.

Title: CENP-B Controls Centromere Formation Depending on the Chromatin Context Content: 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.

Title: Whole-genome sequencing of Oryza brachyantha reveals mechanisms underlying Oryza genome evolution Content: The wild species of the genus Oryza contain a largely untapped reservoir of agronomically important genes for rice improvement. Here we report the 261-Mb de novo assembled genome sequence of Oryza brachyantha. Low activity of long-terminal repeat retrotransposons and massive internal deletions of ancient long-terminal repeat elements lead to the compact genome of Oryza brachyantha. We model 32,038 protein-coding genes in the Oryza brachyantha genome, of which only 70% are located in collinear positions in comparison with the rice genome. Analysing breakpoints of non-collinear genes suggests that double-strand break repair through non-homologous end joining has an important role in gene movement and erosion of collinearity in the Oryza genomes. Transition of euchromatin to heterochromatin in the rice genome is accompanied by segmental and tandem duplications, further expanded by transposable element insertions. The high-quality reference genome sequence of Oryza brachyantha provides an important resource for functional and evolutionary studies in the genus Oryza.

Title: Human artificial chromosomes with alpha satellite-based de novo centromeres show increased frequency of nondisjunction and anaphase lag. Content: 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.

Title: Junctional actin assembly is mediated by Formin-like 2 downstream of Rac1 Content: Epithelial integrity is vitally important, and its deregulation causes early stage cancer. De novo formation of an adherens junction (AJ) between single epithelial cells requires coordinated, spatial actin dynamics, but the mechanisms steering nascent actin polymerization for cell-cell adhesion initiation are not well understood. Here we investigated real-time actin assembly during daughter cell-cell adhesion formation in human breast epithelial cells in 3D environments. We identify formin-like 2 (FMNL2) as being specifically required for actin assembly and turnover at newly formed cell-cell contacts as well as for human epithelial lumen formation. FMNL2 associates with components of the AJ complex involving Rac1 activity and the FMNL2 C terminus. Optogenetic control of Rac1 in living cells rapidly drove FMNL2 to epithelial cell-cell contact zones. Furthermore, Rac1-induced actin assembly and subsequent AJ formation critically depends on FMNL2. These data uncover FMNL2 as a driver for human epithelial AJ formation downstream of Rac1.

Title: A culture-independent sequence-based metagenomics approach to the investigation of an outbreak of Shiga-toxigenic Escherichia coli O104:H4. Content: IMPORTANCE Identification of the bacterium responsible for an outbreak can aid in disease management. However, traditional culture-based diagnosis can be difficult, particularly if no specific diagnostic test is available for an outbreak strain. OBJECTIVE To explore the potential of metagenomics, which is the direct sequencing of DNA extracted from microbiologically complex samples, as an open-ended clinical discovery platform capable of identifying and characterizing bacterial strains from an outbreak without laboratory culture. DESIGN, SETTING, AND PATIENTS In a retrospective investigation, 45 samples were selected from fecal specimens obtained from patients with diarrhea during the 2011 outbreak of Shiga-toxigenic Escherichia coli (STEC) O104:H4 in Germany. Samples were subjected to high-throughput sequencing (August-September 2012), followed by a 3-phase analysis (November 2012-February 2013). In phase 1, a de novo assembly approach was developed to obtain a draft genome of the outbreak strain. In phase 2, the depth of coverage of the outbreak strain genome was determined in each sample. In phase 3, sequences from each sample were compared with sequences from known bacteria to identify pathogens other than the outbreak strain. MAIN OUTCOMES AND MEASURES The recovery of genome sequence data for the purposes of identification and characterization of the outbreak strain and other pathogens from fecal samples. RESULTS During phase 1, a draft genome of the STEC outbreak strain was obtained. During phase 2, the outbreak strain genome was recovered from 10 samples at greater than 10-fold coverage and from 26 samples at greater than 1-fold coverage. Sequences from the Shiga-toxin genes were detected in 27 of 40 STEC-positive samples (67%). In phase 3, sequences from Clostridium difficile, Campylobacter jejuni, Campylobacter concisus, and Salmonella enterica were recovered. CONCLUSIONS AND RELEVANCE These results suggest the potential of metagenomics as a culture-independent approach for the identification of bacterial pathogens during an outbreak of diarrheal disease. Challenges include improving diagnostic sensitivity, speeding up and simplifying workflows, and reducing costs.

Title: The miniature genome of a carnivorous plant Genlisea aurea contains a low number of genes and short non-coding sequences Content: 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.

Title: DNA methylation status predicts cell type-specific enhancer activity. Content: 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.

Title: Quantitative analysis of tumor-derived methylated p16INK4a sequences in plasma, serum, and blood cells of hepatocellular carcinoma patients. Content: PURPOSE AND EXPERIMENTAL DESIGN Using real-time quantitative methylation-specific PCR (RTQ-MSP), we quantified methylated p16INK4a sequences and determined the fractional concentrations of circulating tumor DNA in plasma, serum, and peripheral blood cells collected preoperatively, intraoperatively, and postoperatively from 49 patients with hepatocellular carcinoma (HCC). RESULTS RTQ-MSP was sufficiently sensitive to detect down to 10 genome-equivalents of methylated p16INK4a sequences. Quantitative MSP data were expressed in terms of the methylation index, which was the percentage of bisulfite converted unmethylated and methylated p16INK4a sequences that consisted of methylated p16INK4a sequences. Quantities of methylated p16INK4a sequences were detected in peripheral circulation of 80% (23 of 29) of HCC patients. No significant difference was seen in the detectability and concentrations of methylated p16INK4a sequences (range: 10-4046 genome-equivalents/ml) between preoperative plasma and serum samples from HCC patients. Preoperatively, the p16INK4a methylation indices ranged from 0.2 to 100% and from 0.012 to 0.075% in the patients' plasma and buffy coat samples, respectively. After surgical resection, the median p16INK4a methylation indices in plasma and buffy coat concordantly decreased 12- and 15-fold, respectively. These results demonstrated the clinical usefulness and effectiveness of peripheral blood RTQ-MSP for detecting and monitoring HCC after treatment. Furthermore, none of the intraoperative plasma samples and only two of the intraoperative buffy coat samples were p16INK4a methylation positive. CONCLUSIONS Quantification of epigenetic changes in peripheral blood by RTQ-MSP is useful for the detection and monitoring of HCC.

Title: STAR: ultrafast universal RNA-seq aligner. Content: 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/.

Title: Production of de novo cardiomyocytes: human pluripotent stem cell differentiation and direct reprogramming. Content: Cardiovascular disease is a leading cause of death worldwide. The limited capability of heart tissue to regenerate has prompted methodological developments for creating de novo cardiomyocytes, both in vitro and in vivo. Beyond uses in cell replacement therapy, patient-specific cardiomyocytes may find applications in drug testing, drug discovery, and disease modeling. Recently, approaches for generating cardiomyocytes have expanded to encompass three major sources of starting cells: human pluripotent stem cells (hPSCs), adult heart-derived cardiac progenitor cells (CPCs), and reprogrammed fibroblasts. We discuss state-of-the-art methods for generating de novo cardiomyocytes from hPSCs and reprogrammed fibroblasts, highlighting potential applications and future challenges.

Title: Nucleation of nuclear bodies by RNA Content: 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.

Title: Expression profiles during dedifferentiation in newt lens regeneration revealed by expressed sequence tags Content: 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).

Title: IgBLAST: an immunoglobulin variable domain sequence analysis tool Content: 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.

Title: Chaperoning 5S RNA assembly. Content: In eukaryotes, three of the four ribosomal RNAs (rRNAs)—the 5.8S, 18S, and 25S/28S rRNAs—are processed from a single pre-rRNA transcript and assembled into ribosomes. The fourth rRNA, the 5S rRNA, is transcribed by RNA polymerase III and is assembled into the 5S ribonucleoprotein particle (RNP), containing ribosomal proteins Rpl5/uL18 and Rpl11/uL5, prior to its incorporation into preribosomes. In mammals, the 5S RNP is also a central regulator of the homeostasis of the tumor suppressor p53. The nucleolar localization of the 5S RNP and its assembly into preribosomes are performed by a specialized complex composed of Rpf2 and Rrs1 in yeast or Bxdc1 and hRrs1 in humans. Here we report the structural and functional characterization of the Rpf2-Rrs1 complex alone, in complex with the 5S RNA, and within pre-60S ribosomes. We show that the Rpf2-Rrs1 complex contains a specialized 5S RNA E-loop-binding module, contacts the Rpl5 protein, and also contacts the ribosome assembly factor Rsa4 and the 25S RNA. We propose that the Rpf2-Rrs1 complex establishes a network of interactions that guide the incorporation of the 5S RNP in preribosomes in the initial conformation prior to its rotation to form the central protuberance found in the mature large ribosomal subunit.

Title: Understanding mechanisms underlying human gene expression variation with RNA sequencing Content: Understanding the genetic mechanisms underlying natural variation in gene expression is a central goal of both medical and evolutionary genetics, and studies of expression quantitative trait loci (eQTLs) have become an important tool for achieving this goal. Although all eQTL studies so far have assayed messenger RNA levels using expression microarrays, recent advances in RNA sequencing enable the analysis of transcript variation at unprecedented resolution. We sequenced RNA from 69 lymphoblastoid cell lines derived from unrelated Nigerian individuals that have been extensively genotyped by the International HapMap Project. By pooling data from all individuals, we generated a map of the transcriptional landscape of these cells, identifying extensive use of unannotated untranslated regions and more than 100 new putative protein-coding exons. Using the genotypes from the HapMap project, we identified more than a thousand genes at which genetic variation influences overall expression levels or splicing. We demonstrate that eQTLs near genes generally act by a mechanism involving allele-specific expression, and that variation that influences the inclusion of an exon is enriched within and near the consensus splice sites. Our results illustrate the power of high-throughput sequencing for the joint analysis of variation in transcription, splicing and allele-specific expression across individuals.

Title: MG53 nucleates assembly of cell membrane repair machinery Content: Dynamic membrane repair and remodelling is an elemental process that maintains cell integrity and mediates efficient cellular function. Here we report that MG53, a muscle-specific tripartite motif family protein (TRIM72), is a component of the sarcolemmal membrane-repair machinery. MG53 interacts with phosphatidylserine to associate with intracellular vesicles that traffic to and fuse with sarcolemmal membranes. Mice null for MG53 show progressive myopathy and reduced exercise capability, associated with defective membrane-repair capacity. Injury of the sarcolemmal membrane leads to entry of the extracellular oxidative environment and MG53 oligomerization, resulting in recruitment of MG53-containing vesicles to the injury site. After vesicle translocation, entry of extracellular Ca2+ facilitates vesicle fusion to reseal the membrane. Our data indicate that intracellular vesicle translocation and Ca2+-dependent membrane fusion are distinct steps involved in the repair of membrane damage and that MG53 may initiate the assembly of the membrane repair machinery in an oxidation-dependent manner.

Title: High-throughput sequencing reveals extensive variation in human-specific L1 content in individual human genomes. Content: Using high-throughput sequencing, we devised a technique to determine the insertion sites of virtually all members of the human-specific L1 retrotransposon family in any human genome. Using diagnostic nucleotides, we were able to locate the approximately 800 L1Hs copies corresponding specifically to the pre-Ta, Ta-0, and Ta-1 L1Hs subfamilies, with over 90% of sequenced reads corresponding to human-specific elements. We find that any two individual genomes differ at an average of 285 sites with respect to L1 insertion presence or absence. In total, we assayed 25 individuals, 15 of which are unrelated, at 1139 sites, including 772 shared with the reference genome and 367 nonreference L1 insertions. We show that L1Hs profiles recapitulate genetic ancestry, and determine the chromosomal distribution of these elements. Using these data, we estimate that the rate of L1 retrotransposition in humans is between 1/95 and 1/270 births, and the number of dimorphic L1 elements in the human population with gene frequencies greater than 0.05 is between 3000 and 10,000.

Title: Sequential domain assembly of ribosomal protein S3 drives 40S subunit maturation. Content: 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.

Title: siRNA specificity searching incorporating mismatch tolerance data. Content: UNLABELLED Artificially synthesized short interfering RNAs (siRNAs) are widely used in functional genomics to knock down specific target genes. One ongoing challenge is to guarantee that the siRNA does not elicit off-target effects. Initial reports suggested that siRNAs were highly sequence-specific; however, subsequent data indicates that this is not necessarily the case. It is still uncertain what level of similarity and other rules are required for an off-target effect to be observed, and scoring schemes have not been developed to look beyond simple measures such as the number of mismatches or the number of consecutive matching bases present. We created design rules for predicting the likelihood of a non-specific effect and present a web server that allows the user to check the specificity of a given siRNA in a flexible manner using a combination of methods. The server finds potential off-target matches in the corresponding RefSeq database and ranks them according to a scoring system based on experimental studies of specificity. AVAILABILITY The server is available at http://informatics-eskitis.griffith.edu.au/SpecificityServer.

Title: Integrative Genomics Viewer (IGV): high-performance genomics data visualization and exploration Content: Data visualization is an essential component of genomic data analysis. However, the size and diversity of the data sets produced by today's sequencing and array-based profiling methods present major challenges to visualization tools. The Integrative Genomics Viewer (IGV) is a high-performance viewer that efficiently handles large heterogeneous data sets, while providing a smooth and intuitive user experience at all levels of genome resolution. A key characteristic of IGV is its focus on the integrative nature of genomic studies, with support for both array-based and next-generation sequencing data, and the integration of clinical and phenotypic data. Although IGV is often used to view genomic data from public sources, its primary emphasis is to support researchers who wish to visualize and explore their own data sets or those from colleagues. To that end, IGV supports flexible loading of local and remote data sets, and is optimized to provide high-performance data visualization and exploration on standard desktop systems. IGV is freely available for download from http://www.broadinstitute.org/igv, under a GNU LGPL open-source license.

Title: De novo expression of Trpm4 initiates secondary hemorrhage in spinal cord injury Content: The role of transient receptor potential M4 (Trpm4), an unusual member of the Trp family of ion channels, is poorly understood. Using rodent models of spinal cord injury, we studied involvement of Trpm4 in the progressive expansion of secondary hemorrhage associated with capillary fragmentation, the most destructive mechanism of secondary injury in the central nervous system. Trpm4 mRNA and protein were abundantly upregulated in capillaries preceding their fragmentation and formation of petechial hemorrhages. Trpm4 expression in vitro rendered COS-7 cells highly susceptible to oncotic swelling and oncotic death following ATP depletion. After spinal cord injury, in vivo gene suppression in rats treated with Trpm4 antisense or in Trpm4−/− mice preserved capillary structural integrity, eliminated secondary hemorrhage, yielded a threefold to fivefold reduction in lesion volume and produced a substantial improvement in neurological function. To our knowledge, this is the first example of a Trp channel that must undergo de novo expression for manifestation of central nervous system pathology.

Title: BLAST+: architecture and applications Content: 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.

Title: Improved variant discovery through local re-alignment of short-read next-generation sequencing data using SRMA Content: 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.

Title: Hierarchical phylogenetic models for analyzing multipartite sequence data. Content: Debate exists over how to incorporate information from multipartite sequence data in phylogenetic analyses. Strict combined-data approaches argue for concatenation of all partitions and estimation of one evolutionary history, maximizing the explanatory power of the data. Consensus/independence approaches endorse a two-step procedure where partitions are analyzed independently and then a consensus is determined from the multiple results. Mixtures across the model space of a strict combined-data approach and a priori independent parameters are popular methods to integrate these methods. We propose an alternative middle ground by constructing a Bayesian hierarchical phylogenetic model. Our hierarchical framework enables researchers to pool information across data partitions to improve estimate precision in individual partitions while permitting estimation and testing of tendencies in across-partition quantities. Such across-partition quantities include the distribution from which individual topologies relating the sequences within a partition are drawn. We propose standard hierarchical priors on continuous evolutionary parameters across partitions, while the structure on topologies varies depending on the research problem. We illustrate our model with three examples. We first explore the evolutionary history of the guinea pig (Cavia porcellus) using alignments of 13 mitochondrial genes. The hierarchical model returns substantially more precise continuous parameter estimates than an independent parameter approach without losing the salient features of the data. Second, we analyze the frequency of horizontal gene transfer using 50 prokaryotic genes. We assume an unknown species-level topology and allow individual gene topologies to differ from this with a small estimable probability. Simultaneously inferring the species and individual gene topologies returns a transfer frequency of 17%. We also examine HIV sequences longitudinally sampled from HIV+ patients. We ask whether posttreatment development of CCR5 coreceptor virus represents concerted evolution from middisease CXCR4 virus or reemergence of initial infecting CCR5 virus. The hierarchical model pools partitions from multiple unrelated patients by assuming that the topology for each patient is drawn from a multinomial distribution with unknown probabilities. Preliminary results suggest evolution and not reemergence.

Title: PML induces compaction, TRF2 depletion and DNA damage signaling at telomeres and promotes their alternative lengthening. Content: The alternative lengthening of telomeres (ALT) mechanism allows cancer cells to escape senescence and apoptosis in the absence of active telomerase. A characteristic feature of this pathway is the assembly of ALT-associated promyelocytic leukemia (PML) nuclear bodies (APBs) at telomeres. Here, we dissected the role of APBs in a human ALT cell line by performing an RNA interference screen using an automated 3D fluorescence microscopy platform and advanced 3D image analysis. We identified 29 proteins that affected APB formation, which included proteins involved in telomere and chromatin organization, protein sumoylation and DNA repair. By integrating and extending these findings, we found that APB formation induced clustering of telomere repeats, telomere compaction and concomitant depletion of the shelterin protein TRF2 (also known as TERF2). These APB-dependent changes correlated with the induction of a DNA damage response at telomeres in APBs as evident by a strong enrichment of the phosphorylated form of the ataxia telangiectasia mutated (ATM) kinase. Accordingly, we propose that APBs promote telomere maintenance by inducing a DNA damage response in ALT-positive tumor cells through changing the telomeric chromatin state to trigger ATM phosphorylation.

Title: Metagenomic analysis of microbial consortium from natural crude oil that seeps into the marine ecosystem offshore Southern California Content: Crude oils can be major contaminants of the marine ecosystem and microorganisms play a significant role in the degradation of its main constituents. To increase our understanding of the microbial hydrocarbon degradation process in the marine ecosystem, we collected crude oil from an active seep area located in the Santa Barbara Channel (SBC) and generated a total of about 52 Gb of raw metagenomic sequence data. The assembled data comprised ~500 Mb, representing ~1.1 million genes derived primarily from chemolithoautotrophic bacteria. Members of Oceanospirillales, a bacterial order belonging to the Deltaproteobacteria, recruited less than 2% of the assembled genes within the SBC metagenome. In contrast, the microbial community associated with the oil plume that developed in the aftermath of the Deepwater Horizon (DWH) blowout in 2010, was dominated by Oceanospirillales, which comprised more than 60% of the metagenomic data generated from the DWH oil plume. This suggests that Oceanospirillales might play a less significant role in the microbially mediated hydrocarbon conversion within the SBC seep oil compared to the DWH plume oil. We hypothesize that this difference results from the SBC oil seep being mostly anaerobic, while the DWH oil plume is aerobic. Within the Archaea, the phylum Euryarchaeota, recruited more than 95% of the assembled archaeal sequences from the SBC oil seep metagenome, with more than 50% of the sequences assigned to members of the orders Methanomicrobiales and Methanosarcinales. These orders contain organisms capable of anaerobic methanogenesis and methane oxidation (AOM) and we hypothesize that these orders - and their metabolic capabilities - may be fundamental to the ecology of the SBC oil seep.

Title: Human artificial chromosomes containing chromosome 17 alphoid DNA maintain an active centromere in murine cells but are not stable. Content: 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.

Title: The binding of maize DHN1 to lipid vesicles. Gain of structure and lipid specificity. Content: 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.

Title: Ubiquitin-dependent intramembrane rhomboid protease promotes ERAD of membrane proteins. Content: The ER-associated degradation (ERAD) pathway serves as an important cellular safeguard by directing incorrectly folded and unassembled proteins from the ER to the proteasome. Still, however, little is known about the components mediating ERAD of membrane proteins. Here we show that the evolutionary conserved rhomboid family protein RHBDL4 is a ubiquitin-dependent ER-resident intramembrane protease that is upregulated upon ER stress. RHBDL4 cleaves single-spanning and polytopic membrane proteins with unstable transmembrane helices, leading to their degradation by the canonical ERAD machinery. RHBDL4 specifically binds the AAA+-ATPase p97, suggesting that proteolytic processing and dislocation into the cytosol are functionally linked. The phylogenetic relationship between rhomboids and the ERAD factor derlin suggests that substrates for intramembrane proteolysis and protein dislocation are recruited by a shared mechanism.

Title: The Ndc80 Kinetochore Complex Forms Load-Bearing Attachments to Dynamic Microtubule Tips via Biased Diffusion Content: 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.

Title: Long non-coding RNA profiling of human lymphoid progenitor cells reveals transcriptional divergence of B cell and T cell lineages Content: To elucidate the transcriptional 'landscape' that regulates human lymphoid commitment during postnatal life, we used RNA sequencing to assemble the long non-coding transcriptome across human bone marrow and thymic progenitor cells spanning the earliest stages of B lymphoid and T lymphoid specification. Over 3,000 genes encoding previously unknown long non-coding RNAs (lncRNAs) were revealed through the analysis of these rare populations. Lymphoid commitment was characterized by lncRNA expression patterns that were highly stage specific and were more lineage specific than those of protein-coding genes. Protein-coding genes co-expressed with neighboring lncRNA genes showed enrichment for ontologies related to lymphoid differentiation. The exquisite cell-type specificity of global lncRNA expression patterns independently revealed new developmental relationships among the earliest progenitor cells in the human bone marrow and thymus.

Title: Suppression of intestinal neoplasia by deletion of Dnmt3b Content: Aberrant gene silencing accompanied by DNA methylation is associated with neoplastic progression in many tumors that also show global loss of DNA methylation. Using conditional inactivation of de novo methyltransferase Dnmt3b in Apc(Min/+) mice, we demonstrate that the loss of Dnmt3b has no impact on microadenoma formation, which is considered the earliest stage of intestinal tumor formation. Nevertheless, we observed a significant decrease in the formation of macroscopic colonic adenomas. Interestingly, many large adenomas showed regions with Dnmt3b inactivation, indicating that Dnmt3b is required for initial outgrowth of macroscopic adenomas but is not required for their maintenance. These results support a role for Dnmt3b in the transition stage between microadenoma formation and macroscopic colonic tumor growth and further suggest that Dnmt3b, and by extension de novo methylation, is not required for maintaining tumor growth after this transition stage has occurred.

Title: D3 phosphoinositides and outside-in integrin signaling by glycoprotein IIb-IIIa mediate platelet actin assembly and filopodial extension induced by phorbol 12-myristate 13-acetate. Content: Phorbol 12-myristate 13-acetate (PMA) uncaps a small number of the fast-growing (barbed) ends of actin filaments, thereby eliciting slow actin assembly and extension of filopodia in human blood platelets. These reactions, which also occur in response to immunologic perturbation of the integrin glycoprotein (GP) IIb-IIIa, are sensitive to the phosphoinositide 3-kinase inhibitor wortmannin. Platelets deficient in GPIIb-IIIa integrins or with GPIIb-IIIa function inhibited by calcium chelation or the peptide RGDS have diminished PMA responsiveness. The effects of PMA contrast with thrombin receptor stimulation by >/=5 microM thrombin receptor-activating peptide (TRAP), which causes rapid and massive wortmannin-insensitive actin assembly and lamellar and filopodial extension. However, we show here that wortmannin can inhibit filopod formation if the thrombin receptor is ligated using suboptimal doses (<1 microM) of TRAP. Phosphatidylinositol 3,4-bisphosphate inhibits actin filament severing and capping by human gelsolin in vitro. The findings implicate D3 polyphosphoinositides and integrin signaling in PMA-mediated platelet stimulation and implicate D3 containing phosphoinositides generated in response to protein kinase C activation and GPIIb-IIIa signaling as late-acting intermediates leading to filopodial actin assembly.

Title: Acentrosomal Microtubule Assembly in Mitosis: The Where, When, and How. Content: 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.

Title: Haplotype estimation using sequencing reads. Content: 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.

Title: Functional coupling of RNAP II transcription to spliceosome assembly. Content: 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.

Title: PEP-FOLD: an updated de novo structure prediction server for both linear and disulfide bonded cyclic peptides Content: In the context of the renewed interest of peptides as therapeutics, it is important to have an on-line resource for 3D structure prediction of peptides with well-defined structures in aqueous solution. We present an updated version of PEP-FOLD allowing the treatment of both linear and disulphide bonded cyclic peptides with 9-36 amino acids. The server makes possible to define disulphide bonds and any residue-residue proximity under the guidance of the biologists. Using a benchmark of 34 cyclic peptides with one, two and three disulphide bonds, the best PEP-FOLD models deviate by an average RMS of 2.75 Å from the full NMR structures. Using a benchmark of 37 linear peptides, PEP-FOLD locates lowest-energy conformations deviating by 3 Å RMS from the NMR rigid cores. The evolution of PEP-FOLD comes as a new on-line service to supersede the previous server. The server is available at: http://bioserv.rpbs.univ-paris-diderot.fr/PEP-FOLD.

Title: Acetylation of Histone H3 Lysine 56 Regulates Replication-Coupled Nucleosome Assembly Content: 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.

Title: The Sequence Alignment/Map format and SAMtools Content: 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.

Title: Targeted proteomics reveals compositional dynamics of 60S pre-ribosomes after nuclear export Content: 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.

Title: Methylomic trajectories across human fetal brain development. Content: Epigenetic processes play a key role in orchestrating transcriptional regulation during development. The importance of DNA methylation in fetal brain development is highlighted by the dynamic expression of de novo DNA methyltransferases during the perinatal period and neurodevelopmental deficits associated with mutations in the methyl-CpG binding protein 2 (MECP2) gene. However, our knowledge about the temporal changes to the epigenome during fetal brain development has, to date, been limited. We quantified genome-wide patterns of DNA methylation at ∼ 400,000 sites in 179 human fetal brain samples (100 male, 79 female) spanning 23 to 184 d post-conception. We identified highly significant changes in DNA methylation across fetal brain development at >7% of sites, with an enrichment of loci becoming hypomethylated with fetal age. Sites associated with developmental changes in DNA methylation during fetal brain development were significantly underrepresented in promoter regulatory regions but significantly overrepresented in regions flanking CpG islands (shores and shelves) and gene bodies. Highly significant differences in DNA methylation were observed between males and females at a number of autosomal sites, with a small number of regions showing sex-specific DNA methylation trajectories across brain development. Weighted gene comethylation network analysis (WGCNA) revealed discrete modules of comethylated loci associated with fetal age that are significantly enriched for genes involved in neurodevelopmental processes. This is, to our knowledge, the most extensive study of DNA methylation across human fetal brain development to date, confirming the prenatal period as a time of considerable epigenomic plasticity.

Title: Cortactin regulates cofilin and N-WASp activities to control the stages of invadopodium assembly and maturation Content: Invadopodia are matrix-degrading membrane protrusions in invasive carcinoma cells. The mechanisms regulating invadopodium assembly and maturation are not understood. We have dissected the stages of invadopodium assembly and maturation and show that invadopodia use cortactin phosphorylation as a master switch during these processes. In particular, cortactin phosphorylation was found to regulate cofilin and Arp2/3 complex-dependent actin polymerization. Cortactin directly binds cofilin and inhibits its severing activity. Cortactin phosphorylation is required to release this inhibition so cofilin can sever actin filaments to create barbed ends at invadopodia to support Arp2/3-dependent actin polymerization. After barbed end formation, cortactin is dephosphorylated, which blocks cofilin severing activity thereby stabilizing invadopodia. These findings identify novel mechanisms for actin polymerization in the invadopodia of metastatic carcinoma cells and define four distinct stages of invadopodium assembly and maturation consisting of invadopodium precursor formation, actin polymerization, stabilization, and matrix degradation.

Title: The EMBL nucleotide sequence database. Content: The EMBL Nucleotide Sequence Database (http://www.ebi.ac.uk/embl. html ) constitutes Europe's primary nucleotide sequence resource. DNA and RNA sequences are directly submitted from researchers and genome sequencing groups and collected from the scientific literature and patent applications (Fig. 1). In collaboration with DDBJ and GenBank the database is produced, maintained and distributed at the European Bioinformatics Institute. Database releases are produced quarterly and are distributed on CD-ROM. EBI's network services allow access to the most up-to-date data collection via Internet and World Wide Web interface, providing database searching and sequence similarity facilities plus access to a large number of additional databases.

Title: Spliceosome assembly pathways for different types of alternative splicing converge during commitment to splice site pairing in the A complex. Content: Differential splice site pairing establishes alternative splicing patterns resulting in the generation of multiple mRNA isoforms. This process is carried out by the spliceosome, which is activated by a series of sequential structural rearrangements of its five core snRNPs. To determine when splice sites become functionally paired, we carried out a series of kinetic trap experiments using pre-mRNAs that undergo alternative 5' splice site selection or alternative exon inclusion. We show that commitment to splice site pairing in both cases occurs in the A complex, which is characterized by the ATP-dependent association of the U2 snRNP with the branch point. Interestingly, the timing of splice site pairing is independent of the intron or exon definition modes of splice site recognition. Using the ATP analog ATPgammaS, we showed that ATP hydrolysis is required for splice site pairing independent from U2 snRNP binding to the pre-mRNA. These results identify the A complex as the spliceosomal assembly step dedicated to splice site pairing and suggest that ATP hydrolysis locks splice sites into a splicing pattern after stable U2 snRNP association to the branch point.

Title: A specific amyloid-β protein assembly in the brain impairs memory Content: Memory function often declines with age, and is believed to deteriorate initially because of changes in synaptic function rather than loss of neurons. Some individuals then go on to develop Alzheimer's disease with neurodegeneration. Here we use Tg2576 mice, which express a human amyloid-β precursor protein (APP) variant linked to Alzheimer's disease, to investigate the cause of memory decline in the absence of neurodegeneration or amyloid-β protein amyloidosis. Young Tg2576 mice (< 6 months old) have normal memory and lack neuropathology, middle-aged mice (6–14 months old) develop memory deficits without neuronal loss, and old mice (> 14 months old) form abundant neuritic plaques containing amyloid-β (refs 3–6). We found that memory deficits in middle-aged Tg2576 mice are caused by the extracellular accumulation of a 56-kDa soluble amyloid-β assembly, which we term Aβ*56 (Aβ star 56). Aβ*56 purified from the brains of impaired Tg2576 mice disrupts memory when administered to young rats. We propose that Aβ*56 impairs memory independently of plaques or neuronal loss, and may contribute to cognitive deficits associated with Alzheimer's disease.

Title: Differential assembly of inwardly rectifying K+ channel subunits, Kir4.1 and Kir5.1, in brain astrocytes. Content: The inwardly rectifying K+ channel subunit Kir5.1 is expressed abundantly in the brain, but its precise distribution and function are still largely unknown. Because Kir5.1 is co-expressed with Kir4.1 in retinal glial Muller cells, we have compared the biochemical and immunological properties of Kir5.1 and Kir4.1 in the mouse brain. Immunoprecipitation experiments suggested that brain expressed at least two subsets of Kir channels, heteromeric Kir4.1/5.1 and homomeric Kir4.1. Immunolabeling using specific antibodies showed that channels comprising Kir4.1 and Kir5.1 subunits were assembled in a region-specific fashion. Heteromeric Kir4.1/5.1 was identified in the neocortex and in the glomeruli of the olfactory bulb. Homomeric Kir4.1 was confined to the hippocampus and the thalamus. Homomeric Kir5.1 was not identified. Kir4.1/5.1 and Kir4.1 expression appeared to occur only in astrocytes, specifically in the membrane domains facing the pia mater and blood vessels or in the processes surrounding synapses. Both Kir4.1/5.1 and Kir4.1 could be associated with PDZ domain-containing syntrophins, which might be involved in the subcellular targeting of these astrocyte Kir channels. Because heteromeric Kir4.1/5.1 and homomeric Kir4.1 have distinct ion channel properties (Tanemoto, M., Kittaka, N., Inanobe, A., and Kurachi, Y. (2000) J. Physiol. (Lond.) 525, 587-592 and Tucker, S. J., Imbrici, P., Salvatore, L., D'Adamo, M. C., and Pessia, M. (2000) J. Biol. Chem. 275, 16404-16407), it is plausible that these channels play differential physiological roles in the K+ -buffering action of brain astrocytes in a region-specific manner.

Title: CLASPs link focal adhesion-associated microtubule capture to localized exocytosis and adhesion site turnover Content: Turnover of integrin-based focal adhesions (FAs) with the extracellular matrix (ECM) is essential for coordinated cell movement. In collectively migrating human keratinocytes, FAs assemble near the leading edge, grow and mature as a result of contractile forces and disassemble underneath the advancing cell body. We report that clustering of microtubule-associated CLASP1 and CLASP2 proteins around FAs temporally correlates with FA turnover. CLASPs and LL5β (also known as PHLDB2), which recruits CLASPs to FAs, facilitate FA disassembly. CLASPs are further required for FA-associated ECM degradation, and matrix metalloprotease inhibition slows FA disassembly similarly to CLASP or PHLDB2 (LL5β) depletion. Finally, CLASP-mediated microtubule tethering at FAs establishes an FA-directed transport pathway for delivery, docking and localized fusion of exocytic vesicles near FAs. We propose that CLASPs couple microtubule organization, vesicle transport and cell interactions with the ECM, establishing a local secretion pathway that facilitates FA turnover by severing cell-matrix connections.

Title: Exploring the temporal structure of heterochronous sequences using TempEst (formerly Path-O-Gen) Content: Gene sequences sampled at different points in time can be used to infer molecular phylogenies on a natural timescale of months or years, provided that the sequences in question undergo measurable amounts of evolutionary change between sampling times. Data sets with this property are termed heterochronous and have become increasingly common in several fields of biology, most notably the molecular epidemiology of rapidly evolving viruses. Here we introduce the cross-platform software tool, TempEst (formerly known as Path-O-Gen), for the visualization and analysis of temporally sampled sequence data. Given a molecular phylogeny and the dates of sampling for each sequence, TempEst uses an interactive regression approach to explore the association between genetic divergence through time and sampling dates. TempEst can be used to (1) assess whether there is sufficient temporal signal in the data to proceed with phylogenetic molecular clock analysis, and (2) identify sequences whose genetic divergence and sampling date are incongruent. Examination of the latter can help identify data quality problems, including errors in data annotation, sample contamination, sequence recombination, or alignment error. We recommend that all users of the molecular clock models implemented in BEAST first check their data using TempEst prior to analysis.

Title: Automated genome sequence analysis and annotation. Content: MOTIVATION Large-scale genome projects generate a rapidly increasing number of sequences, most of them biochemically uncharacterized. Research in bioinformatics contributes to the development of methods for the computational characterization of these sequences. However, the installation and application of these methods require experience and are time consuming. RESULTS We present here an automatic system for preliminary functional annotation of protein sequences that has been applied to the analysis of sets of sequences from complete genomes, both to refine overall performance and to make new discoveries comparable to those made by human experts. The GeneQuiz system includes a Web-based browser that allows examination of the evidence leading to an automatic annotation and offers additional information, views of the results, and links to biological databases that complement the automatic analysis. System structure and operating principles concerning the use of multiple sequence databases, underlying sequence analysis tools, lexical analyses of database annotations and decision criteria for functional assignments are detailed. The system makes automatic quality assessments of results based on prior experience with the underlying sequence analysis tools; overall error rates in functional assignment are estimated at 2.5-5% for cases annotated with highest reliability ('clear' cases). Sources of over-interpretation of results are discussed with proposals for improvement. A conservative definition for reporting 'new findings' that takes account of database maturity is presented along with examples of possible kinds of discoveries (new function, family and superfamily) made by the system. System performance in relation to sequence database coverage, database dynamics and database search methods is analysed, demonstrating the inherent advantages of an integrated automatic approach using multiple databases and search methods applied in an objective and repeatable manner. AVAILABILITY The GeneQuiz system is publicly available for analysis of protein sequences through a Web server at http://www.sander.ebi.ac. uk/gqsrv/submit

Title: Cutadapt removes adapter sequences from high-throughput sequencing reads Content: 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/

Title: Construction of plant bacterial artificial chromosome (BAC) libraries: an illustrated guide. Content: J. Agric. Genomics, 5 ABSTRACT Bacterial artificial chromosome (BAC) libraries have become invaluable tools in plant genetic research. However, it is difficult for new practitioners to create plant BAC libraries de novo because published protocols are not particularly detailed, and plant cells possess features that make isolation of clean, high molecular weight DNA troublesome. In this document we present an illustrated, step-by-step protocol for constructing plant BAC libraries. This protocol is sufficiently detailed to be of use to both new and experienced investigators. We hope that by reducing the obstacles to BAC cloning in plants, we will foster new and accelerated progress in plant genomics.

Title: Dislocation of a type I membrane protein requires interactions between membrane-spanning segments within the lipid bilayer. Content: The human cytomegalovirus gene product US11 causes rapid degradation of class I major histocompatibility complex (MHCI) heavy chains by inducing their dislocation from the endoplasmic reticulum (ER) and subsequent degradation by the proteasome. This set of reactions resembles the endogenous cellular quality control pathway that removes misfolded or unassembled proteins from the ER. We show that the transmembrane domain (TMD) of US11 is essential for MHCI heavy chain dislocation, but dispensable for MHCI binding. A Gln residue at position 192 in the US11 TMD is crucial for the ubiquitination and degradation of MHCI heavy chains. Cells that express US11 TMD mutants allow formation of MHCI-beta2m complexes, but their rate of egress from the ER is significantly impaired. Further mutagenesis data are consistent with the presence of an alpha-helical structure in the US11 TMD essential for MHCI heavy chain dislocation. The failure of US11 TMD mutants to catalyze dislocation is a unique instance in which a polar residue in the TMD of a type I membrane protein is required for that protein's function. Targeting of MHCI heavy chains for dislocation by US11 thus requires the formation of interhelical hydrogen bonds within the ER membrane.

Title: DNA Deamination Mediates Innate Immunity to Retroviral Infection Content: CEM15/APOBEC3G is a cellular protein required for resistance to infection by virion infectivity factor (Vif)-deficient human immunodeficiency virus (HIV). Here, using a murine leukemia virus (MLV)-based system, we provide evidence that CEM15/APOBEC3G is a DNA deaminase that is incorporated into virions during viral production and subsequently triggers massive deamination of deoxycytidine to deoxyuridine within the retroviral minus (first)-strand cDNA, thus providing a probable trigger for viral destruction. Furthermore, HIV Vif can protect MLV from this CEM15/APOBEC3G-dependent restriction. These findings imply that targeted DNA deamination is a major strategy of innate immunity to retroviruses and likely also contributes to the sequence variation observed in many viruses (including HIV).

Title: Model-based Analysis of ChIP-Seq (MACS) Content: We present Model-based Analysis of ChIP-Seq data, MACS, which analyzes data generated by short read sequencers such as Solexa's Genome Analyzer. MACS empirically models the shift size of ChIP-Seq tags, and uses it to improve the spatial resolution of predicted binding sites. MACS also uses a dynamic Poisson distribution to effectively capture local biases in the genome, allowing for more robust predictions. MACS compares favorably to existing ChIP-Seq peak-finding algorithms, and is freely available.

Title: position-specific gap penalties Content: 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.

Title: Identification and characterization of a common set of complex I assembly intermediates in mitochondria from patients with complex I deficiency. Content: 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.

Title: Trimmomatic: a flexible trimmer for Illumina sequence data Content: 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.

Title: Evolution of metastasis revealed by mutational landscapes of chemically induced skin cancers Content: Human tumors show a high level of genetic heterogeneity, but the processes that influence the timing and route of metastatic dissemination of the subclones are unknown. Here we have used whole-exome sequencing of 103 matched benign, malignant and metastatic skin tumors from genetically heterogeneous mice to demonstrate that most metastases disseminate synchronously from the primary tumor, supporting parallel rather than linear evolution as the predominant model of metastasis. Shared mutations between primary carcinomas and their matched metastases have the distinct A-to-T signature of the initiating carcinogen dimethylbenzanthracene, but non-shared mutations are primarily G-to-T, a signature associated with oxidative stress. The existence of carcinomas that either did or did not metastasize in the same host animal suggests that there are tumor-intrinsic factors that influence metastatic seeding. We also demonstrate the importance of germline polymorphisms in determining allele-specific mutations, and we identify somatic genetic alterations that are specifically related to initiation of carcinogenesis by Hras or Kras mutations. Mouse tumors that mimic the genetic heterogeneity of human cancers can aid our understanding of the clonal evolution of metastasis and provide a realistic model for the testing of novel therapies.

Title: Specificity Residues Determine Binding Affinity for Two-Component Signal Transduction Systems Content: UNLABELLED Two-component systems (TCS) comprise histidine kinases and their cognate response regulators and allow bacteria to sense and respond to a wide variety of signals. Histidine kinases (HKs) phosphorylate and dephosphorylate their cognate response regulators (RRs) in response to stimuli. In general, these reactions appear to be highly specific and require an appropriate association between the HK and RR proteins. The Myxococcus xanthus genome encodes one of the largest repertoires of signaling proteins in bacteria (685 open reading frames [ORFs]), including at least 127 HKs and at least 143 RRs. Of these, 27 are bona fide NtrC-family response regulators, 21 of which are encoded adjacent to their predicted cognate kinases. Using system-wide profiling methods, we determined that the HK-NtrC RR pairs display a kinetic preference during both phosphotransfer and phosphatase functions, thereby defining cognate signaling systems in M. xanthus. Isothermal titration calorimetry measurements indicated that cognate HK-RR pairs interact with dissociation constants (Kd) of approximately 1 µM, while noncognate pairs had no measurable binding. Lastly, a chimera generated between the histidine kinase, CrdS, and HK1190 revealed that residues conferring phosphotransfer and phosphatase specificity dictate binding affinity, thereby establishing discrete protein-protein interactions which prevent cross talk. The data indicate that binding affinity is a critical parameter governing system-wide signaling fidelity for bacterial signal transduction proteins. IMPORTANCE Using in vitro phosphotransfer and phosphatase profiling assays and isothermal titration calorimetry, we have taken a system-wide approach to demonstrate specificity for a family of two-component signaling proteins in Myxococcus xanthus. Our results demonstrate that previously identified specificity residues dictate binding affinity and that phosphatase specificity follows phosphotransfer specificity for cognate HK-RR pairs. The data indicate that preferential binding affinity is the basis for signaling fidelity in bacterial two-component systems.

Title: Folding and Self-Assembly of the TatA Translocation Pore Based on a Charge Zipper Mechanism Content: We propose a concept for the folding and self-assembly of the pore-forming TatA complex from the Twin-arginine translocase and of other membrane proteins based on electrostatic "charge zippers. " Each subunit of TatA consists of a transmembrane segment, an amphiphilic helix (APH), and a C-terminal densely charged region (DCR). The sequence of charges in the DCR is complementary to the charge pattern on the APH, suggesting that the protein can be "zipped up" by a ladder of seven salt bridges. The length of the resulting hairpin matches the lipid bilayer thickness, hence a transmembrane pore could self-assemble via intra- and intermolecular salt bridges. The steric feasibility was rationalized by molecular dynamics simulations, and experimental evidence was obtained by monitoring the monomer-oligomer equilibrium of specific charge mutants. Similar "charge zippers" are proposed for other membrane-associated proteins, e.g., the biofilm-inducing peptide TisB, the human antimicrobial peptide dermcidin, and the pestiviral E(RNS) protein.

Title: N- and E-cadherins in Xenopus are specifically required in the neural and non-neural ectoderm, respectively, for F-actin assembly and morphogenetic movements. Content: 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.

Title: An RNA-sequencing transcriptome and splicing database of glia, neurons, and vascular cells of the cerebral cortex. Content: The major cell classes of the brain differ in their developmental processes, metabolism, signaling, and function. To better understand the functions and interactions of the cell types that comprise these classes, we acutely purified representative populations of neurons, astrocytes, oligodendrocyte precursor cells, newly formed oligodendrocytes, myelinating oligodendrocytes, microglia, endothelial cells, and pericytes from mouse cerebral cortex. We generated a transcriptome database for these eight cell types by RNA sequencing and used a sensitive algorithm to detect alternative splicing events in each cell type. Bioinformatic analyses identified thousands of new cell type-enriched genes and splicing isoforms that will provide novel markers for cell identification, tools for genetic manipulation, and insights into the biology of the brain. For example, our data provide clues as to how neurons and astrocytes differ in their ability to dynamically regulate glycolytic flux and lactate generation attributable to unique splicing of PKM2, the gene encoding the glycolytic enzyme pyruvate kinase. This dataset will provide a powerful new resource for understanding the development and function of the brain. To ensure the widespread distribution of these datasets, we have created a user-friendly website (http://web.stanford.edu/group/barres_lab/brain_rnaseq.html) that provides a platform for analyzing and comparing transciption and alternative splicing profiles for various cell classes in the brain.

Title: NCBI GEO: archive for functional genomics data sets—update Content: The Gene Expression Omnibus (GEO, http://www.ncbi.nlm.nih.gov/geo/) is an international public repository for high-throughput microarray and next-generation sequence functional genomic data sets submitted by the research community. The resource supports archiving of raw data, processed data and metadata which are indexed, cross-linked and searchable. All data are freely available for download in a variety of formats. GEO also provides several web-based tools and strategies to assist users to query, analyse and visualize data. This article reports current status and recent database developments, including the release of GEO2R, an R-based web application that helps users analyse GEO data.

Title: Eic1 links Mis18 with the CCAN/Mis6/Ctf19 complex to promote CENP-A assembly Content: 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.

Title: Exome sequencing identifies somatic gain-of-function PPM1D mutations in brainstem gliomas Content: Gliomas arising in the brainstem and thalamus are devastating tumors that are difficult to surgically resect. To determine the genetic and epigenetic landscape of these tumors, we performed exomic sequencing of 14 brainstem gliomas (BSGs) and 12 thalamic gliomas. We also performed targeted mutational analysis of an additional 24 such tumors and genome-wide methylation profiling of 45 gliomas. This study led to the discovery of tumor-specific mutations in PPM1D, encoding wild-type p53-induced protein phosphatase 1D (WIP1), in 37.5% of the BSGs that harbored hallmark H3F3A mutations encoding p. Lys27Met substitutions. PPM1D mutations were mutually exclusive with TP53 mutations in BSG and attenuated p53 activation in vitro. PPM1D mutations were truncating alterations in exon 6 that enhanced the ability of PPM1D to suppress the activation of the DNA damage response checkpoint protein CHK2. These results define PPM1D as a frequent target of somatic mutation and as a potential therapeutic target in brainstem gliomas.

Title: Integrated molecular genetic profiling of pediatric high-grade gliomas reveals key differences with the adult disease. Content: PURPOSE To define copy number alterations and gene expression signatures underlying pediatric high-grade glioma (HGG). PATIENTS AND METHODS We conducted a high-resolution analysis of genomic imbalances in 78 de novo pediatric HGGs, including seven diffuse intrinsic pontine gliomas, and 10 HGGs arising in children who received cranial irradiation for a previous cancer using single nucleotide polymorphism microarray analysis. Gene expression was analyzed with gene expression microarrays for 53 tumors. Results were compared with publicly available data from adult tumors. RESULTS Significant differences in copy number alterations distinguish childhood and adult glioblastoma. PDGFRA was the predominant target of focal amplification in childhood HGG, including diffuse intrinsic pontine gliomas, and gene expression analyses supported an important role for deregulated PDGFRalpha signaling in pediatric HGG. No IDH1 hotspot mutations were found in pediatric tumors, highlighting molecular differences with adult secondary glioblastoma. Pediatric and adult glioblastomas were clearly distinguished by frequent gain of chromosome 1q (30% v 9%, respectively) and lower frequency of chromosome 7 gain (13% v 74%, respectively) and 10q loss (35% v 80%, respectively). PDGFRA amplification and 1q gain occurred at significantly higher frequency in irradiation-induced tumors, suggesting that these are initiating events in childhood gliomagenesis. A subset of pediatric HGGs showed minimal copy number changes. CONCLUSION Integrated molecular profiling showed substantial differences in the molecular features underlying pediatric and adult HGG, indicating that findings in adult tumors cannot be simply extrapolated to younger patients. PDGFRalpha may be a useful target for pediatric HGG, including diffuse pontine gliomas.

Title: Brain-Specific Phosphorylation of MeCP2 Regulates Activity-Dependent Bdnf Transcription, Dendritic Growth, and Spine Maturation Content: Mutations or duplications in MECP2 cause Rett and Rett-like syndromes, neurodevelopmental disorders characterized by mental retardation, motor dysfunction, and autistic behaviors. MeCP2 is expressed in many mammalian tissues and functions as a global repressor of transcription; however, the molecular mechanisms by which MeCP2 dysfunction leads to the neural-specific phenotypes of RTT remain poorly understood. Here, we show that neuronal activity and subsequent calcium influx trigger the de novo phosphorylation of MeCP2 at serine 421 (S421) by a CaMKII-dependent mechanism. MeCP2 S421 phosphorylation is induced selectively in the brain in response to physiological stimuli. Significantly, we find that S421 phosphorylation controls the ability of MeCP2 to regulate dendritic patterning, spine morphogenesis, and the activity-dependent induction of Bdnf transcription. These findings suggest that, by triggering MeCP2 phosphorylation, neuronal activity regulates a program of gene expression that mediates nervous system maturation and that disruption of this process in individuals with mutations in MeCP2 may underlie the neural-specific pathology of RTT.

Title: A library of yeast transcription factor motifs reveals a widespread function for Rsc3 in targeting nucleosome exclusion at promoters. Content: The sequence specificity of DNA-binding proteins is the primary mechanism by which the cell recognizes genomic features. Here, we describe systematic determination of yeast transcription factor DNA-binding specificities. We obtained binding specificities for 112 DNA-binding proteins representing 19 distinct structural classes. One-third of the binding specificities have not been previously reported. Several binding sequences have striking genomic distributions relative to transcription start sites, supporting their biological relevance and suggesting a role in promoter architecture. Among these are Rsc3 binding sequences, containing the core CGCG, which are found preferentially approximately 100 bp upstream of transcription start sites. Mutation of RSC3 results in a dramatic increase in nucleosome occupancy in hundreds of proximal promoters containing a Rsc3 binding element, but has little impact on promoters lacking Rsc3 binding sequences, indicating that Rsc3 plays a broad role in targeting nucleosome exclusion at yeast promoters.

Title: Resequencing studies of nonmodel organisms using closely related reference genomes: optimal experimental designs and bioinformatics approaches for population genomics. Content: Decreasing costs of next-generation sequencing (NGS) experiments have made a wide range of genomic questions open for study with nonmodel organisms. However, experimental designs and analysis of NGS data from less well-known species are challenging because of the lack of genomic resources. In this work, we investigate the performance of alternative experimental designs and bioinformatics approaches in estimating variability and neutrality tests based on the site-frequency-spectrum (SFS) from individual resequencing data. We pay particular attention to challenges faced in the study of nonmodel organisms, in particular the absence of a species-specific reference genome, although phylogenetically close genomes are assumed to be available. We compare the performance of three alternative bioinformatics approaches – genotype calling, genotype–haplotype calling and direct estimation without calling genotypes. We find that relying on genotype calls provides biased estimates of population genetic statistics at low to moderate read depth (2–8X). Genotype–haplotype calling returns more accurate estimates irrespective of the divergence to the reference genome, but requires moderate depth (8–20X). Direct estimation without calling genotypes returns the most accurate estimates of variability and of most SFS tests investigated, including at low read depth (2–4X). Studies without species-specific reference genome should thus aim for low read depth and avoid genotype calling whenever individual genotypes are not essential. Otherwise, aiming for moderate to high depth at the expense of number of individuals, and using genotype–haplotype calling, is recommended.

Title: LncRNA NEAT1/let-7a-5p axis regulates the cisplatin resistance in nasopharyngeal carcinoma by targeting Rsf-1 and modulating the Ras-MAPK pathway. Content: The long non-coding RNA nuclear paraspeckle assembly transcript 1 (NEAT1) was reported to be upregulated and be involved in oncogenic growth and drug resistance in nasopharyngeal carcinoma (NPC). However, the exact roles of NEAT1 and its underlying mechanisms in the drug resistance of NPC remain largely unclear. In this study, the expressions of NEAT1, let-72-5p and Rsf-1 mRNA were detected by reverse transcription-quantitative polymerase chain reaction (RT-qPCR). The effects of NEAT1 and let-72-5p on cell proliferation and cisplatin resistance of NPC cells were investigated by 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyl tetrazolium bromide (MTT) assay and 5-ethynyl-20-deoxyuridine (EdU) assay. Western blot analysis was performed to detect the protein levels of Rsf-1, Ras, p-Raf1, Raf1, p-MEK1, MEK1, p-ERK1/2 and ERK1/2. Xenograft tumor assay was done to elucidate the role of NEAT1 involved in NPC tumor growth in vivo. We found that NEAT1 was upregulated and let-7a-5p was downregulated in NPC tissues, as well as NPC cell lines. Inhibition of NEAT1 markedly repressed the cisplatin resistance of NPC cells. NEAT1 was demonstrated to interact with let-7a-5p. Besides, a negative correlation between NEAT1 and let-7a-5p expression was observed in NPC tissues. Rsf-1 was confirmed as a target of let-7a-5p. NEAT1 remarkably reversed the inhibitory effect of let-7q-5p on the cisplatin resistance of NPC cells in vitro. Additionally, NEAT1 knockdown inhibited the Ras-MAPK pathway in NPC cells. NEAT1 knockdown suppressed tumor growth in the presence of cisplatin in vivo. Overall, these findings suggest that NEAT1/let-7a-5p axis regulates the cisplatin resistance in NPC by targeting Rsf-1 and modulating the Ras-MAPK signaling pathway.

Title: Frozen robust multiarray analysis (fRMA). Content: 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.

Title: The DNA Methylome of Human Peripheral Blood Mononuclear Cells Content: DNA methylation plays an important role in biological processes in human health and disease. Recent technological advances allow unbiased whole-genome DNA methylation (methylome) analysis to be carried out on human cells. Using whole-genome bisulfite sequencing at 24.7-fold coverage (12.3-fold per strand), we report a comprehensive (92.62%) methylome and analysis of the unique sequences in human peripheral blood mononuclear cells (PBMC) from the same Asian individual whose genome was deciphered in the YH project. PBMC constitute an important source for clinical blood tests world-wide. We found that 68.4% of CpG sites and <0.2% of non-CpG sites were methylated, demonstrating that non-CpG cytosine methylation is minor in human PBMC. Analysis of the PBMC methylome revealed a rich epigenomic landscape for 20 distinct genomic features, including regulatory, protein-coding, non-coding, RNA-coding, and repeat sequences. Integration of our methylome data with the YH genome sequence enabled a first comprehensive assessment of allele-specific methylation (ASM) between the two haploid methylomes of any individual and allowed the identification of 599 haploid differentially methylated regions (hDMRs) covering 287 genes. Of these, 76 genes had hDMRs within 2 kb of their transcriptional start sites of which >80% displayed allele-specific expression (ASE). These data demonstrate that ASM is a recurrent phenomenon and is highly correlated with ASE in human PBMCs. Together with recently reported similar studies, our study provides a comprehensive resource for future epigenomic research and confirms new sequencing technology as a paradigm for large-scale epigenomics studies.

Title: Circulating tumor DNA as a non-invasive substitute to metastasis biopsy for tumor genotyping and personalized medicine in a prospective trial across all tumor types. Content: Cell-free tumor DNA (ctDNA) has the potential to enable non-invasive diagnostic tests for personalized medicine in providing similar molecular information as that derived from invasive tumor biopsies. The histology-independent phase II SHIVA trial matches patients with targeted therapeutics based on previous screening of multiple somatic mutations using metastatic biopsies. To evaluate the utility of ctDNA in this trial, as an ancillary study we performed de novo detection of somatic mutations using plasma DNA compared to metastasis biopsies in 34 patients covering 18 different tumor types, scanning 46 genes and more than 6800 COSMIC mutations with a multiplexed next-generation sequencing panel. In 27 patients, 28 of 29 mutations identified in metastasis biopsies (97%) were detected in matched ctDNA. Among these 27 patients, one additional mutation was found in ctDNA only. In the seven other patients, mutation detection from metastasis biopsy failed due to inadequate biopsy material, but was successful in all plasma DNA samples providing three more potential actionable mutations. These results suggest that ctDNA analysis is a potential alternative and/or replacement to analyses using costly, harmful and lengthy tissue biopsies of metastasis, irrespective of cancer type and metastatic site, for multiplexed mutation detection in selecting personalized therapies based on the patient's tumor genetic content.

Title: A Protein Inventory of Human Ribosome Biogenesis Reveals an Essential Function of Exportin 5 in 60S Subunit Export Content: The assembly of ribosomal subunits in eukaryotes is a complex, multistep process so far mostly studied in yeast. In S. cerevisiae, more than 200 factors including ribosomal proteins and trans-acting factors are required for the ordered assembly of 40S and 60S ribosomal subunits. To date, only few human homologs of these yeast ribosome synthesis factors have been characterized. Here, we used a systematic RNA interference (RNAi) approach to analyze the contribution of 464 candidate factors to ribosomal subunit biogenesis in human cells. The screen was based on visual readouts, using inducible, fluorescent ribosomal proteins as reporters. By performing computer-based image analysis utilizing supervised machine-learning techniques, we obtained evidence for a functional link of 153 human proteins to ribosome synthesis. Our data show that core features of ribosome assembly are conserved from yeast to human, but differences exist for instance with respect to 60S subunit export. Unexpectedly, our RNAi screen uncovered a requirement for the export receptor Exportin 5 (Exp5) in nuclear export of 60S subunits in human cells. We show that Exp5, like the known 60S exportin Crm1, binds to pre-60S particles in a RanGTP-dependent manner. Interference with either Exp5 or Crm1 function blocks 60S export in both human cells and frog oocytes, whereas 40S export is compromised only upon inhibition of Crm1. Thus, 60S subunit export is dependent on at least two RanGTP-binding exportins in vertebrate cells.

Title: Assessing sequence comparison methods with reliable structurally identified distant evolutionary relationships. Content: 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.

Title: Unifying Viral Genetics and Human Transportation Data to Predict the Global Transmission Dynamics of Human Influenza H3N2 Content: Information on global human movement patterns is central to spatial epidemiological models used to predict the behavior of influenza and other infectious diseases. Yet it remains difficult to test which modes of dispersal drive pathogen spread at various geographic scales using standard epidemiological data alone. Evolutionary analyses of pathogen genome sequences increasingly provide insights into the spatial dynamics of influenza viruses, but to date they have largely neglected the wealth of information on human mobility, mainly because no statistical framework exists within which viral gene sequences and empirical data on host movement can be combined. Here, we address this problem by applying a phylogeographic approach to elucidate the global spread of human influenza subtype H3N2 and assess its ability to predict the spatial spread of human influenza A viruses worldwide. Using a framework that estimates the migration history of human influenza while simultaneously testing and quantifying a range of potential predictive variables of spatial spread, we show that the global dynamics of influenza H3N2 are driven by air passenger flows, whereas at more local scales spread is also determined by processes that correlate with geographic distance. Our analyses further confirm a central role for mainland China and Southeast Asia in maintaining a source population for global influenza diversity. By comparing model output with the known pandemic expansion of H1N1 during 2009, we demonstrate that predictions of influenza spatial spread are most accurate when data on human mobility and viral evolution are integrated. In conclusion, the global dynamics of influenza viruses are best explained by combining human mobility data with the spatial information inherent in sampled viral genomes. The integrated approach introduced here offers great potential for epidemiological surveillance through phylogeographic reconstructions and for improving predictive models of disease control.

Title: A conditionally lethal yeast mutant blocked at the first step in glycosyl phosphatidylinositol anchor synthesis. Content: Glycosyl phosphatidylinositols (GPIs) anchor many proteins to the surface of eukaryotic cells and may also serve as sorting signals on proteins and participate in signal transduction. We have isolated a Saccharomyces cerevisiae GPI anchoring mutant, gpi1, using a colony screen for cells blocked in [3H]inositol incorporation into protein. The gpi1 mutant is defective in vitro in the synthesis of N-acetylglucosaminyl phosphatidylinositol, the first intermediate in GPI synthesis, and is also temperature-sensitive for growth. Completion of the first step in GPI assembly is therefore required for growth of the unicellular eukaryote S. cerevisiae. GPI synthesis could therefore be exploited as a target for antifungal or antiparasitic agents.

Title: The First Myriapod Genome Sequence Reveals Conservative Arthropod Gene Content and Genome Organisation in the Centipede Strigamia maritima Content: Myriapods (e.g., centipedes and millipedes) display a simple homonomous body plan relative to other arthropods. All members of the class are terrestrial, but they attained terrestriality independently of insects. Myriapoda is the only arthropod class not represented by a sequenced genome. We present an analysis of the genome of the centipede Strigamia maritima. It retains a compact genome that has undergone less gene loss and shuffling than previously sequenced arthropods, and many orthologues of genes conserved from the bilaterian ancestor that have been lost in insects. Our analysis locates many genes in conserved macro-synteny contexts, and many small-scale examples of gene clustering. We describe several examples where S. maritima shows different solutions from insects to similar problems. The insect olfactory receptor gene family is absent from S. maritima, and olfaction in air is likely effected by expansion of other receptor gene families. For some genes S. maritima has evolved paralogues to generate coding sequence diversity, where insects use alternate splicing. This is most striking for the Dscam gene, which in Drosophila generates more than 100,000 alternate splice forms, but in S. maritima is encoded by over 100 paralogues. We see an intriguing linkage between the absence of any known photosensory proteins in a blind organism and the additional absence of canonical circadian clock genes. The phylogenetic position of myriapods allows us to identify where in arthropod phylogeny several particular molecular mechanisms and traits emerged. For example, we conclude that juvenile hormone signalling evolved with the emergence of the exoskeleton in the arthropods and that RR-1 containing cuticle proteins evolved in the lineage leading to Mandibulata. We also identify when various gene expansions and losses occurred. The genome of S. maritima offers us a unique glimpse into the ancestral arthropod genome, while also displaying many adaptations to its specific life history.

Title: Sequence-based species delimitation for the DNA taxonomy of undescribed insects. Content: Cataloging the very large number of undescribed species of insects could be greatly accelerated by automated DNA based approaches, but procedures for large-scale species discovery from sequence data are currently lacking. Here, we use mitochondrial DNA variation to delimit species in a poorly known beetle radiation in the genus Rivacindela from arid Australia. Among 468 individuals sampled from 65 sites and multiple morphologically distinguishable types, sequence variation in three mtDNA genes (cytochrome oxidase subunit 1, cytochrome b, 16S ribosomal RNA) was strongly partitioned between 46 or 47 putative species identified with quantitative methods of species recognition based on fixed unique ("diagnostic") characters. The boundaries between groups were also recognizable from a striking increase in branching rate in clock-constrained calibrated trees. Models of stochastic lineage growth (Yule models) were combined with coalescence theory to develop a new likelihood method that determines the point of transition from species-level (speciation and extinction) to population-level (coalescence) evolutionary processes. Fitting the location of the switches from speciation to coalescent nodes on the ultrametric tree of Rivacindela produced a transition in branching rate occurring at 0.43 Mya, leading to an estimate of 48 putative species (confidence interval for the threshold ranging from 47 to 51 clusters within 2 logL units). Entities delimited in this way exhibited biological properties of traditionally defined species, showing coherence of geographic ranges, broad congruence with morphologically recognized species, and levels of sequence divergence typical for closely related species of insects. The finding of discontinuous evolutionary groupings that are readily apparent in patterns of sequence variation permits largely automated species delineation from DNA surveys of local communities as a scaffold for taxonomy in this poorly known insect group.

Title: COHCAP: an integrative genomic pipeline for single-nucleotide resolution DNA methylation analysis Content: COHCAP (City of Hope CpG Island Analysis Pipeline) is an algorithm to analyze single-nucleotide resolution DNA methylation data produced by either an Illumina methylation array or targeted bisulfite sequencing. The goal of the COHCAP algorithm is to identify CpG islands that show a consistent pattern of methylation among CpG sites. COHCAP is currently the only DNA methylation package that provides integration with gene expression data to identify a subset of CpG islands that are most likely to regulate downstream gene expression, and it can generate lists of differentially methylated CpG islands with ∼50% concordance with gene expression from both cell line data and heterogeneous patient data. For example, this article describes known breast cancer biomarkers (such as estrogen receptor) with a negative correlation between DNA methylation and gene expression. COHCAP also provides visualization for quality control metrics, regions of differential methylation and correlation between methylation and gene expression. This software is freely available at https://sourceforge.net/projects/cohcap/.

Title: Improved tools for biological sequence comparison. Content: 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.

Title: Haplotype-based variant detection from short-read sequencing Content: The direct detection of haplotypes from short-read DNA sequencing data requires changes to existing small-variant detection methods. Here, we develop a Bayesian statistical framework which is capable of modeling multiallelic loci in sets of individuals with non-uniform copy number. We then describe our implementation of this framework in a haplotype-based variant detector, FreeBayes.

Title: Histone chaperones and nucleosome assembly. Content: Recent structures of the nucleosome core particle reveal details of histone-histone and histone-DNA interactions. These structures have now set the stage for understanding chromatin assembly and dynamics during replication and transcription. Histone chaperones and chromatin remodeling complexes are important in both of these processes. The nucleosome and its protein core, the histone octamer, have twofold symmetry, which histone chaperones may use to bind core histones. Recent studies suggest that the nucleoplasmin pentamer may mediate histone storage, sperm chromatin decondensation and nucleosome assembly, by dimerizing to form a decamer. In this model, histone binding on the lateral surface of the chaperone involves stereospecific interactions and a shared twofold axis.

Title: Plant Nuclear RNA Polymerase IV Mediates siRNA and DNA Methylation-Dependent Heterochromatin Formation Content: All eukaryotes have three nuclear DNA-dependent RNA polymerases, namely, Pol I, II, and III. Interestingly, plants have catalytic subunits for a fourth nuclear polymerase, Pol IV. Genetic and biochemical evidence indicates that Pol IV does not functionally overlap with Pol I, II, or III and is nonessential for viability. However, disruption of the Pol IV catalytic subunit genes NRPD1 or NRPD2 inhibits heterochromatin association into chromocenters, coincident with losses in cytosine methylation at pericentromeric 5S gene clusters and AtSN1 retroelements. Loss of CG, CNG, and CNN methylation in Pol IV mutants implicates a partnership between Pol IV and the methyltransferase responsible for RNA-directed de novo methylation. Consistent with this hypothesis, 5S gene and AtSN1 siRNAs are essentially eliminated in Pol IV mutants. The data suggest that Pol IV helps produce siRNAs that target de novo cytosine methylation events required for facultative heterochromatin formation and higher-order heterochromatin associations.

Title: Monitoring spatiotemporal biogenesis of macromolecular assemblies by pulse-chase epitope labeling. Content: 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.

Title: YAP/TAZ-Dependent Reprogramming of Colonic Epithelium Links ECM Remodeling to Tissue Regeneration Content: Tissue regeneration requires dynamic cellular adaptation to the wound environment. It is currently unclear how this is orchestrated at the cellular level and how cell fate is affected by severe tissue damage. Here we dissect cell fate transitions during colonic regeneration in a mouse dextran sulfate sodium (DSS) colitis model, and we demonstrate that the epithelium is transiently reprogrammed into a primitive state. This is characterized by de novo expression of fetal markers as well as suppression of markers for adult stem and differentiated cells. The fate change is orchestrated by remodeling the extracellular matrix (ECM), increased FAK/Src signaling, and ultimately YAP/TAZ activation. In a defined cell culture system recapitulating the extracellular matrix remodeling observed in vivo, we show that a collagen 3D matrix supplemented with Wnt ligands is sufficient to sustain endogenous YAP/TAZ and induce conversion of cell fate. This provides a simple model for tissue regeneration, implicating cellular reprogramming as an essential element.

Title: The expanding problem of adipose depot remodeling and postnatal adipocyte progenitor recruitment. Content: The rising incidence of obesity and associated metabolic diseases has increased the urgency in understanding all aspects of adipose tissue biology. This includes the function of adipocytes, how adipose tissue expands in obesity, and how expanded adipose tissues in adults can impact physiology. Here, we highlight the growing appreciation for the importance of de novo adipocyte differentiation to adipose tissue expansion in adult humans and animals. We detail recent efforts to identify adipose precursor populations that contribute to the physiological postnatal recruitment of white, brown, and beige adipocytes in mice, and summarize new data that reveal the complexity of adipose tissue development in vivo.

Title: A novel ATAC-seq approach reveals lineage-specific reinforcement of the open chromatin landscape via cooperation between BAF and p63 Content: BACKGROUND Open chromatin regions are correlated with active regulatory elements in development and are dysregulated in diseases. The BAF (SWI/SNF) complex is essential for development, and has been demonstrated to remodel reconstituted chromatin in vitro and to control the accessibility of a few individual regions in vivo. However, it remains unclear where and how BAF controls the open chromatin landscape to regulate developmental processes, such as human epidermal differentiation. RESULTS Using a novel "on-plate" ATAC-sequencing approach for profiling open chromatin landscapes with a low number of adherent cells, we demonstrate that the BAF complex is essential for maintaining 11.6 % of open chromatin regions in epidermal differentiation. These BAF-dependent open chromatin regions are highly cell-type-specific and are strongly enriched for binding sites for p63, a master epidermal transcription factor. The DNA sequences of p63 binding sites intrinsically favor nucleosome formation and are inaccessible in other cell types without p63 to prevent ectopic activation. In epidermal cells, BAF and p63 mutually recruit each other to maintain 14,853 open chromatin regions. We further demonstrate that BAF and p63 cooperatively position nucleosomes away from p63 binding sites and recruit transcriptional machinery to control tissue differentiation. CONCLUSIONS BAF displays high specificity in controlling the open chromatin landscape during epidermal differentiation by cooperating with the master transcription factor p63 to maintain lineage-specific open chromatin regions.

Title: Four newly located pseudouridylate residues in Escherichia coli 23S ribosomal RNA are all at the peptidyltransferase center: analysis by the application of a new sequencing technique. Content: A new technique has been developed for the facile location of pseudouridylate (psi) residues in any RNA molecule. The method uses two known modification procedures which in combination uniquely identify U residues which have been converted into psi. The first procedure involves reaction of all U-like and G-like residues with N-cyclohexyl-N'-beta-(4-methylmorpholinium)ethylcarbodiimide p-tosylate (CMC), followed by alkaline removal of all CMC groups except those linked to the N3 of psi. This stops reverse transcription, resulting in a gel band which identifies the U residue. The second procedure is uridine-specific hydrazinolysis which cleaves the RNA chain at all U residues and produces a gel band upon reverse transcription. psi residues, being resistant to hydrazinolysis, are not cleaved and do not stop reverse transcription. This leads to the absence of a band at psi residues. The combined method can also distinguish psi from 5-methyluridine, 4-thiouridine, uridine-5-oxyacetic acid, and 2-thio-5-methylaminomethyluridine as shown by treating rRNA and tRNA species known to contain these modified bases at defined sites. By this procedure, four new sites for psi in Escherichia coli 23S RNA were discovered, and one was disproven. The four new sites are at positions 2457, 2504, 2580, and 2605. The erroneous site is at position 2555. These four new psi residues, which are all in or within 2-3 residues of the peptidyltransferase ring, are thus in a position to play a functional and/or structural role at the peptidyltransferase center.(ABSTRACT TRUNCATED AT 250 WORDS)

Title: Chronophin mediates an ATP-sensing mechanism for cofilin dephosphorylation and neuronal cofilin-actin rod formation. Content: Actin and its key regulatory component, cofilin, are found together in large rod-shaped assemblies in neurons subjected to energy stress. Such inclusions are also enriched in Alzheimer's disease brain, and appear in transgenic models of neurodegeneration. Neuronal insults, such as energy loss and/or oxidative stress, result in rapid dephosphorylation of the cellular cofilin pool prior to its assembly into rod-shaped inclusions. Although these events implicate a role for phosphatases in cofilin rod formation, a mechanism linking energy stress, phosphocofilin turnover, and subsequent rod assembly has been elusive. We demonstrate the ATP-sensitive interaction of the cofilin phosphatase chronophin (CIN) with the chaperone hsp90 to form a biosensor that mediates cofilin/actin rod formation. Our results suggest a model whereby attenuated interactions between CIN and hsp90 during ATP depletion enhance CIN-dependent cofilin dephosphorylation and consequent rod assembly, thereby providing a mechanism for the formation of pathological actin/cofilin aggregates during neurodegenerative energy flux.

Title: Membrane lipids in invadopodia and podosomes: Key structures for cancer invasion and metastasis Content: Invadopodia are extracellular matrix (ECM)-degrading protrusions formed by invasive cancer cells. Podosomes are structures functionally similar to invadopodia that are found in oncogene-transformed fibroblasts and monocyte-derived cells, including macrophages and osteoclasts. These structures are thought to play important roles in the pericellular remodeling of ECM during cancer invasion and metastasis. Much effort has been directed toward identification of the molecular components and regulators of invadopodia/podosomes, which could be therapeutic targets in the treatment of malignant cancers. However, it remains largely unknown how these components are assembled into invadopodia/podosomes and how the assembly process is spatially and temporally regulated. This review will summarize recent progress on the molecular mechanisms of invadopodia/podosome formation, with strong emphasis on the roles of lipid rafts and phosphoinositides.

Title: Canonical Wnt signaling regulates organ-specific assembly and differentiation of CNS vasculature. Content: Every organ depends on blood vessels for oxygen and nutrients, but the vasculature associated with individual organs can be structurally and molecularly diverse. The central nervous system (CNS) vasculature consists of a tightly sealed endothelium that forms the blood-brain barrier, whereas blood vessels of other organs are more porous. Wnt7a and Wnt7b encode two Wnt ligands produced by the neuroepithelium of the developing CNS coincident with vascular invasion. Using genetic mouse models, we found that these ligands directly target the vascular endothelium and that the CNS uses the canonical Wnt signaling pathway to promote formation and CNS-specific differentiation of the organ's vasculature.

Title: Memory and Modularity in Cell-Fate Decision Making Content: Genetically identical cells sharing an environment can display markedly different phenotypes. It is often unclear how much of this variation derives from chance, external signals, or attempts by individual cells to exert autonomous phenotypic programs. By observing thousands of cells for hundreds of consecutive generations under constant conditions, we dissect the stochastic decision between a solitary, motile state and a chained, sessile state in Bacillus subtilis. We show that the motile state is 'memoryless', exhibiting no autonomous control over the time spent in the state. In contrast, the time spent as connected chains of cells is tightly controlled, enforcing coordination among related cells in the multicellular state. We show that the three-protein regulatory circuit governing the decision is modular, as initiation and maintenance of chaining are genetically separable functions. As stimulation of the same initiating pathway triggers biofilm formation, we argue that autonomous timing allows a trial commitment to multicellularity that external signals could extend.

Title: Cerebral cortex assembly: generating and reprogramming projection neuron diversity Content: The mammalian cerebral cortex is responsible for the highest levels of associative, cognitive and motor functions. In the central nervous system (CNS) the cortex stands as a prime example of extreme neuronal diversity, broadly classified into excitatory projection neurons (PNs) and inhibitory interneurons (INs). We review here recent progress made in understanding the strategies and mechanisms that shape PN diversity during embryogenesis, and discuss how PN classes may be maintained, postnatally, for the life of the organism. In addition, we consider the intriguing possibility that PNs may be amenable to directed reprogramming of their class-specific features to allow enhanced cortical plasticity in the adult.

Title: A framework for improving microRNA prediction in non-human genomes Content: 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.