Sol Katzman

An RNA gene expressed during cortical development evolved rapidly in humans

The developmental and evolutionary mechanisms behind the emergence of human-specific brain features remain largely unknown. However, the recent ability to compare our genome to that of our closest relative, the chimpanzee, provides new avenues to link genetic and phenotypic changes in the evolution of the human brain. We devised a ranking of regions in the human genome that show significant evolutionary acceleration. Here we report that the most dramatic of these ‘human accelerated regions’, HAR1, is part of a novel RNA gene (HAR1F) that is expressed specifically in Cajal–Retzius neurons in the developing human neocortex from 7 to 19 gestational weeks, a crucial period for cortical neuron specification and migration. HAR1F is co-expressed with reelin, a product of Cajal–Retzius neurons that is of fundamental importance in specifying the six-layer structure of the human cortex. HAR1 and the other human accelerated regions provide new candidates in the search for uniquely human biology.

FragSeq: transcriptome-wide RNA structure probing using high-throughput sequencing

Classical approaches to determine structures of noncoding RNA (ncRNA) probed only one RNA at a time with enzymes and chemicals, using gel electrophoresis to identify reactive positions. To accelerate RNA structure inference, we developed fragmentation sequencing (FragSeq), a high-throughput RNA structure probing method that uses high-throughput RNA sequencing of fragments generated by digestion with nuclease P1, which specifically cleaves single-stranded nucleic acids. In experiments probing the entire mouse nuclear transcriptome, we accurately and simultaneously mapped single-stranded RNA regions in multiple ncRNAs with known structure. We probed in two cell types to verify reproducibility. We also identified and experimentally validated structured regions in ncRNAs with, to our knowledge, no previously reported probing data.

An evolutionary arms race between KRAB zinc-finger genes ZNF91/93 and SVA/L1 retrotransposons

Throughout evolution primate genomes have been modified by waves of retrotransposon insertions. For each wave, the host eventually finds a way to repress retrotransposon transcription and prevent further insertions. In mouse embryonic stem cells, transcriptional silencing of retrotransposons requires KAP1 (also known as TRIM28) and its repressive complex, which can be recruited to target sites by KRAB zinc-finger (KZNF) proteins such as murine-specific ZFP809 which binds to integrated murine leukaemia virus DNA elements and recruits KAP1 to repress them. KZNF genes are one of the fastest growing gene families in primates and this expansion is hypothesized to enable primates to respond to newly emerged retrotransposons. However, the identity of KZNF genes battling retrotransposons currently active in the human genome, such as SINE-VNTR-Alu (SVA) and long interspersed nuclear element 1 (L1), is unknown. Here we show that two primate-specific KZNF genes rapidly evolved to repress these two distinct retrotransposon families shortly after they began to spread in our ancestral genome. ZNF91 underwent a series of structural changes 8–12 million years ago that enabled it to repress SVA elements. ZNF93 evolved earlier to repress the primate L1 lineage until ∼12.5 million years ago when the L1PA3-subfamily of retrotransposons escaped ZNF93’s restriction through the removal of the ZNF93-binding site. Our data support a model where KZNF gene expansion limits the activity of newly emerged retrotransposon classes, and this is followed by mutations in these retrotransposons to evade repression, a cycle of events that could explain the rapid expansion of lineage-specific KZNF genes.

SAM‐T04: What is new in protein–structure prediction for CASP6

The SAM‐T04 method for predicting protein structures uses a single protocol across the entire range of targets, from comparative modeling to new folds. This protocol is similar to the SAM‐T02 protocol used in CASP5, but has improvements in the iterative search for similar sequences in finding and aligning templates, in creating fragment libraries, in generating protein conformations, and in scoring the conformations. The automatic procedure made some improvements over simply selecting an alignment to the highest‐scoring template, and human intervention made substantial improvements over the automatic procedure. The main improvements made by human intervention were from adding constraints to build (or retain) β‐sheets and from splitting multidomain proteins into separate domains. The uniform protocol was moderately successful across the entire range of target difficulty, but was somewhat less successful than other approaches in CASP6 on the comparative modeling targets. Proteins 2005;Suppl 7:135–142.

Human tRNA genes function as chromatin insulators

Insulators help separate active chromatin domains from silenced ones. In yeast, gene promoters act as insulators to block the spread of Sir and HP1 mediated silencing while in metazoans most insulators are multipartite autonomous entities. tDNAs are repetitive sequences dispersed throughout the human genome and we now show that some of these tDNAs can function as insulators in human cells. Using computational methods, we identified putative human tDNA insulators. Using silencer blocking, transgene protection and repressor blocking assays we show that some of these tDNA‐containing fragments can function as barrier insulators in human cells. We find that these elements also have the ability to block enhancers from activating RNA pol II transcribed promoters. Characterization of a putative tDNA insulator in human cells reveals that the site possesses chromatin signatures similar to those observed at other better‐characterized eukaryotic insulators. Enhanced 4C analysis demonstrates that the tDNA insulator makes long‐range chromatin contacts with other tDNAs and ETC sites but not with intervening or flanking RNA pol II transcribed genes.

Competition between pre-mRNAs for the splicing machinery drives global regulation of splicing

During meiosis in yeast, global splicing efficiency increases and then decreases. Here we provide evidence that splicing improves due to reduced competition for the splicing machinery. The timing of this regulation corresponds to repression and reactivation of ribosomal protein genes (RPGs) during meiosis. In vegetative cells, RPG repression by rapamycin treatment also increases splicing efficiency. Downregulation of the RPG-dedicated transcription factor gene IFH1 genetically suppresses two spliceosome mutations, prp11-1 and prp4-1, and globally restores splicing efficiency in prp4-1 cells. We conclude that the splicing apparatus is limiting and that pre-messenger RNAs compete. Splicing efficiency of a pre-mRNA therefore depends not just on its own concentration and affinity for limiting splicing factor(s), but also on those of competing pre-mRNAs. Competition between RNAs for limiting processing factors appears to be a general condition in eukaryotes for a variety of posttranscriptional control mechanisms including microRNA (miRNA) repression, polyadenylation, and splicing.

Ongoing GC-biased evolution is widespread in the human genome and enriched near recombination hot spots.

Fast evolving regions of many metazoan genomes show a bias toward substitutions that change weak (A,T) into strong (G,C) base pairs. Single-nucleotide polymorphisms (SNPs) do not share this pattern, suggesting that it results from biased fixation rather than biased mutation. Supporting this hypothesis, analyses of polymorphism in specific regions of the human genome have identified a positive correlation between weak to strong (W→S) SNPs and derived allele frequency (DAF), suggesting that SNPs become increasingly GC biased over time, especially in regions of high recombination. Using polymorphism data generated by the 1000 Genomes Project from 179 individuals from 4 human populations, we evaluated the extent and distribution of ongoing GC-biased evolution in the human genome. We quantified GC fixation bias by comparing the DAFs of W→S mutations and S→W mutations using a Mann–Whitney U test. Genome-wide, W→S SNPs have significantly higher DAFs than S→W SNPs. This pattern is widespread across the human genome but varies in magnitude along the chromosomes. We found extreme GC-biased evolution in neighborhoods of recombination hot spots, a significant correlation between GC bias and recombination rate, and an inverse correlation between GC bias and chromosome arm length. These findings demonstrate the presence of ongoing fixation bias favoring G and C alleles throughout the human genome and suggest that the bias is caused by a recombination-associated process, such as GC-biased gene conversion.

Double Minute Chromosomes in Glioblastoma Multiforme Are Revealed by Precise Reconstruction of Oncogenic Amplicons

DNA sequencing offers a powerful tool in oncology based on the precise definition of structural rearrangements and copy number in tumor genomes. Here, we describe the development of methods to compute copy number and detect structural variants to locally reconstruct highly rearranged regions of the tumor genome with high precision from standard, short-read, paired-end sequencing datasets. We find that circular assemblies are the most parsimonious explanation for a set of highly amplified tumor regions in a subset of glioblastoma multiforme samples sequenced by The Cancer Genome Atlas (TCGA) consortium, revealing evidence for double minute chromosomes in these tumors. Further, we find that some samples harbor multiple circular amplicons and, in some cases, further rearrangements occurred after the initial amplicon-generating event. Fluorescence in situ hybridization analysis offered an initial confirmation of the presence of double minute chromosomes. Gene content in these assemblies helps identify likely driver oncogenes for these amplicons. RNA-seq data available for one double minute chromosome offered additional support for our local tumor genome assemblies, and identified the birth of a novel exon made possible through rearranged sequences present in the double minute chromosomes. Our method was also useful for analysis of a larger set of glioblastoma multiforme tumors for which exome sequencing data are available, finding evidence for oncogenic double minute chromosomes in more than 20% of clinical specimens examined, a frequency consistent with previous estimates.

Predict-2nd

MOTIVATION Predictions of protein local structure, derived from sequence alignment information alone, provide visualization tools for biologists to evaluate the importance of amino acid residue positions of interest in the absence of X-ray crystal/NMR structures or homology models. They are also useful as inputs to sequence analysis and modeling tools, such as hidden Markov models (HMMs), which can be used to search for homology in databases of known protein structure. In addition, local structure predictions can be used as a component of cost functions in genetic algorithms that predict protein tertiary structure. We have developed a program (predict-2nd) that trains multilayer neural networks and have applied it to numerous local structure alphabets, tuning network parameters such as the number of layers, the number of units in each layer and the window sizes of each layer. We have had the most success with four-layer networks, with gradually increasing window sizes at each layer. RESULTS Because the four-layer neural nets occasionally get trapped in poor local optima, our training protocol now uses many different random starts, with short training runs, followed by more training on the best performing networks from the short runs. One recent addition to the program is the option to add a guide sequence to the profile inputs, increasing the number of inputs per position by 20. We find that use of a guide sequence provides a small but consistent improvement in the predictions for several different local-structure alphabets. AVAILABILITY Local structure prediction with the methods described here is available for use online at http://www.soe.ucsc.edu/compbio/SAM_T08/T08-query.html. The source code and example networks for PREDICT-2ND are available at http://www.soe.ucsc.edu/~karplus/predict-2nd/ A required C++ library is available at http://www.soe.ucsc.edu/~karplus/ultimate/

Foxg1 Coordinates the Switch from Nonradially to Radially Migrating Glutamatergic Subtypes in the Neocortex through Spatiotemporal Repression

The specification of neuronal subtypes in the cerebral cortex proceeds in a temporal manner; however, the regulation of the transitions between the sequentially generated subtypes is poorly understood. Here, we report that the forkhead box transcription factor Foxg1 coordinates the production of neocortical projection neurons through the global repression of a default gene program. The delayed activation of Foxg1 was necessary and sufficient to induce deep-layer neurogenesis, followed by a sequential wave of upper-layer neurogenesis. A genome-wide analysis revealed that Foxg1 binds to mammalian-specific noncoding sequences to repress over 12 transcription factors expressed in early progenitors, including Ebf2/3, Dmrt3, Dmrta1, and Eya2. These findings reveal an unexpected prolonged competence of progenitors to initiate corticogenesis at a progressed stage during development and identify Foxg1 as a critical initiator of neocorticogenesis through spatiotemporal repression, a system that balances the production of nonradially and radially migrating glutamatergic subtypes during mammalian cortical expansion.