James R. Iben

High incidence of unrecognized visceral/neurological late-onset Niemann-Pick disease, type C1, predicted by analysis of massively parallel sequencing data sets

Purpose:Niemann-Pick disease type C (NPC) is a recessive, neurodegenerative, lysosomal storage disease caused by mutations in either NPC1 or NPC2. The diagnosis is difficult and frequently delayed. Ascertainment is likely incomplete because of both these factors and because the full phenotypic spectrum may not have been fully delineated. Given the recent development of a blood-based diagnostic test and the development of potential therapies, understanding the incidence of NPC and defining at-risk patient populations are important.Method:We evaluated data from four large, massively parallel exome sequencing data sets. Variant sequences were identified and classified as pathogenic or nonpathogenic based on a combination of literature review and bioinformatic analysis. This methodology provided an unbiased approach to determining the allele frequency.Results:Our data suggest an incidence rate for NPC1 and NPC2 of 1/92,104 and 1/2,858,998, respectively. Evaluation of common NPC1 variants, however, suggests that there may be a late-onset NPC1 phenotype with a markedly higher incidence, on the order of 1/19,000-1/36,000.Conclusion:We determined a combined incidence of classical NPC of 1/89,229, or 1.12 affected patients per 100,000 conceptions, but predict incomplete ascertainment of a late-onset phenotype of NPC1. This finding strongly supports the need for increased screening of potential patients.Genet Med 18 1, 41–48.

Lack of tRNA Modification Isopentenyl-A37 Alters mRNA Decoding and Causes Metabolic Deficiencies in Fission Yeast

ABSTRACT tRNA isopentenyltransferases (Tit1) modify tRNA position 37, adjacent to the anticodon, to N6-isopentenyladenosine (i6A37) in all cells, yet the tRNA subsets selected for modification vary among species, and their relevance to phenotypes is unknown. We examined the function of i6A37 in Schizosaccharomyces pombe tit1+ and tit1-Δ cells by using a β-galactosidase codon-swap reporter whose catalytic activity is sensitive to accurate decoding of codon 503. i6A37 increased the activity of tRNACys at a cognate codon and that of tRNATyr at a near-cognate codon, suggesting that i6A37 promotes decoding activity generally and increases fidelity at cognate codons while decreasing fidelity at noncognate codons. S. pombe cells lacking tit1+ exhibit slow growth in glycerol or rapamycin. While existing data link wobble base U34 modifications to translation of functionally related mRNAs, whether this might extend to the anticodon-adjacent position 37 was unknown. Indeed, we found a biased presence of i6A37-cognate codons in high-abundance mRNAs for ribosome subunits and energy metabolism, congruent with the observed phenotypes and the idea that i6A37 promotes translational efficiency. Polysome profiles confirmed the decreased translational efficiency of mRNAs in tit1-Δ cells. Because subsets of i6A37-tRNAs differ among species, as do their cognate codon-sensitive mRNAs, these genomic variables may underlie associated phenotypic differences.

Global ‘bootprinting’ reveals the elastic architecture of the yeast TFIIIB–TFIIIC transcription complex in vivo

TFIIIB and TFIIIC are multi-subunit factors required for transcription by RNA polymerase III. We present a genome-wide high-resolution footprint map of TFIIIB–TFIIIC complexes in Saccharomyces cerevisiae, obtained by paired-end sequencing of micrococcal nuclease-resistant DNA. On tRNA genes, TFIIIB and TFIIIC form stable complexes with the same distinctive occupancy pattern but in mirror image, termed ‘bootprints’. Global analysis reveals that the TFIIIB–TFIIIC transcription complex exhibits remarkable structural elasticity: tRNA genes vary significantly in length but remain protected by TFIIIC. Introns, when present, are markedly less protected. The RNA polymerase III transcription terminator is flexibly accommodated within the transcription complex and, unexpectedly, plays a major structural role by delimiting its 3′-boundary. The ETC sites, where TFIIIC binds without TFIIIB, exhibit different bootprints, suggesting that TFIIIC forms complexes involving other factors. We confirm six ETC sites and report a new site (ETC10). Surprisingly, TFIIIC, but not TFIIIB, interacts with some centromeric nucleosomes, suggesting that interactions between TFIIIC and the centromere may be important in the 3D organization of the nucleus.

tRNAomics: tRNA gene copy number variation and codon use provide bioinformatic evidence of a new anticodon:codon wobble pair in a eukaryote

tRNA genes are interspersed throughout eukaryotic DNA, contributing to genome architecture and evolution in addition to translation of the transcriptome. Codon use correlates with tRNA gene copy number in noncomplex organisms including yeasts. Synonymous codons impact translation with various outcomes, dependent on relative tRNA abundances. Availability of whole-genome sequences allowed us to examine tRNA gene copy number variation (tgCNV) and codon use in four Schizosaccharomyces species and Saccharomyces cerevisiae. tRNA gene numbers vary from 171 to 322 in the four Schizosaccharomyces despite very high similarity in other features of their genomes. In addition, we performed whole-genome sequencing of several related laboratory strains of Schizosaccharomyces pombe and found tgCNV at a cluster of tRNA genes. We examined for the first time effects of wobble rules on correlation of tRNA gene number and codon use and showed improvement for S. cerevisiae and three of the Schizosaccharomyces species. In contrast, correlation in Schizosaccharomyces japonicus is poor due to markedly divergent tRNA gene content, and much worsened by the wobble rules. In japonicus, some tRNA iso-acceptor genes are absent and others are greatly reduced relative to the other yeasts, while genes for synonymous wobble iso-acceptors are amplified, indicating wobble use not apparent in any other eukaryote. We identified a subset of japonicus-specific wobbles that improves correlation of codon use and tRNA gene content in japonicus. We conclude that tgCNV is high among Schizo species and occurs in related laboratory strains of S. pombe (and expectedly other species), and tRNAome-codon analyses can provide insight into species-specific wobble decoding.

RSC-dependent constructive and destructive interference between opposing arrays of phased nucleosomes in yeast

RSC and SWI/SNF are related ATP-dependent chromatin remodeling machines that move nucleosomes, regulating access to DNA. We addressed their roles in nucleosome phasing relative to transcription start sites in yeast. SWI/SNF has no effect on phasing at the global level. In contrast, RSC depletion results in global nucleosome repositioning: Both upstream and downstream nucleosomal arrays shift toward the nucleosome-depleted region (NDR), with no change in spacing, resulting in a narrower and partly filled NDR. The global picture of RSC-depleted chromatin represents the average of a range of chromatin structures, with most genes showing a shift of the +1 or the -1 nucleosome into the NDR. Using RSC ChIP data reported by others, we show that RSC occupancy is highest on the coding regions of heavily transcribed genes, though not at their NDRs. We propose that RSC has a role in restoring chromatin structure after transcription. Analysis of gene pairs in different orientations demonstrates that phasing patterns reflect competition between phasing signals emanating from neighboring NDRs. These signals may be in phase, resulting in constructive interference and a regular array, or out of phase, resulting in destructive interference and fuzzy positioning. We propose a modified barrier model, in which a stable complex located at the NDR acts as a bidirectional phasing barrier. In RSC-depleted cells, this barrier has a smaller footprint, resulting in narrower NDRs. Thus, RSC plays a critical role in organizing yeast chromatin.

RNA Polymerase III Output Is Functionally Linked to tRNA Dimethyl-G26 Modification.

Control of the differential abundance or activity of tRNAs can be important determinants of gene regulation. RNA polymerase (RNAP) III synthesizes all tRNAs in eukaryotes and it derepression is associated with cancer. Maf1 is a conserved general repressor of RNAP III under the control of the target of rapamycin (TOR) that acts to integrate transcriptional output and protein synthetic demand toward metabolic economy. Studies in budding yeast have indicated that the global tRNA gene activation that occurs with derepression of RNAP III via maf1-deletion is accompanied by a paradoxical loss of tRNA-mediated nonsense suppressor activity, manifested as an antisuppression phenotype, by an unknown mechanism. We show that maf1-antisuppression also occurs in the fission yeast S. pombe amidst general activation of RNAP III. We used tRNA-HydroSeq to document that little changes occurred in the relative levels of different tRNAs in maf1Δ cells. By contrast, the efficiency of N2,N2-dimethyl G26 (m2 2G26) modification on certain tRNAs was decreased in response to maf1-deletion and associated with antisuppression, and was validated by other methods. Over-expression of Trm1, which produces m2 2G26, reversed maf1-antisuppression. A model that emerges is that competition by increased tRNA levels in maf1Δ cells leads to m2 2G26 hypomodification due to limiting Trm1, reducing the activity of suppressor-tRNASerUCA and accounting for antisuppression. Consistent with this, we show that RNAP III mutations associated with hypomyelinating leukodystrophy decrease tRNA transcription, increase m2 2G26 efficiency and reverse antisuppression. Extending this more broadly, we show that a decrease in tRNA synthesis by treatment with rapamycin leads to increased m2 2G26 modification and that this response is conserved among highly divergent yeasts and human cells.

LEDGF/p75 interacts with mRNA splicing factors and targets HIV-1 integration to highly spliced genes

The host chromatin-binding factor LEDGF/p75 interacts with HIV-1 integrase and directs integration to active transcription units. To understand how LEDGF/p75 recognizes transcription units, we sequenced 1 million HIV-1 integration sites isolated from cultured HEK293T cells. Analysis of integration sites showed that cancer genes were preferentially targeted, raising concerns about using lentivirus vectors for gene therapy. Additional analysis led to the discovery that introns and alternative splicing contributed significantly to integration site selection. These correlations were independent of transcription levels, size of transcription units, and length of the introns. Multivariate analysis with five parameters previously found to predict integration sites showed that intron density is the strongest predictor of integration density in transcription units. Analysis of previously published HIV-1 integration site data showed that integration density in transcription units in mouse embryonic fibroblasts also correlated strongly with intron number, and this correlation was absent in cells lacking LEDGF. Affinity purification showed that LEDGF/p75 is associated with a number of splicing factors, and RNA sequencing (RNA-seq) analysis of HEK293T cells lacking LEDGF/p75 or the LEDGF/p75 integrase-binding domain (IBD) showed that LEDGF/p75 contributes to splicing patterns in half of the transcription units that have alternative isoforms. Thus, LEDGF/p75 interacts with splicing factors, contributes to exon choice, and directs HIV-1 integration to transcription units that are highly spliced.

tRNA gene copy number variation in humans

The human tRNAome consists of more than 500 interspersed tRNA genes comprising 51 anticodon families of largely unequal copy number. We examined tRNA gene copy number variation (tgCNV) in six individuals; two kindreds of two parents and a child, using high coverage whole genome sequence data. Such differences may be important because translation of some mRNAs is sensitive to the relative amounts of tRNAs and because tRNA competition determines translational efficiency vs. fidelity and production of native vs. misfolded proteins. We identified several tRNA gene clusters with CNV, which in some cases were part of larger iterations. In addition there was an isolated tRNALysCUU gene that was absent as a homozygous deletion in one of the parents. When assessed by semiquantitative PCR in 98 DNA samples representing a wide variety of ethnicities, this allele was found deleted in hetero- or homozygosity in all groups at ~50% frequency. This is the first report of copy number variation of human tRNA genes. We conclude that tgCNV exists at significant levels among individual humans and discuss the results in terms of genetic diversity and prior genome wide association studies (GWAS) that suggest the importance of the ratio of tRNALys isoacceptors in Type-2 diabetes.

Determination of the allelic frequency in Smith–Lemli–Opitz syndrome by analysis of massively parallel sequencing data sets

Data from massively parallel sequencing or ‘Next Generation Sequencing’ of the human exome has reached a critical mass in both public and private databases, in that these collections now allow researchers to critically evaluate population genetics in a manner that was not feasible a decade ago. The ability to determine pathogenic allele frequencies by evaluation of the full coding sequences and not merely a single nucleotide polymorphism (SNP) or series of SNPs will lead to more accurate estimations of incidence. For demonstrative purposes, we analyzed the causative gene for the disorder Smith–Lemli–Opitz Syndrome (SLOS), the 7‐dehydrocholesterol reductase (DHCR7) gene and determined both the carrier frequency for DHCR7 mutations, and predicted an expected incidence of the disorder. Estimations of the incidence of SLOS have ranged widely from 1:10,000 to 1:70,000 while the carrier frequency has been reported as high as 1 in 30. Using four exome data sets with a total of 17,836 chromosomes, we ascertained a carrier frequency of pathogenic DHRC7 mutations of 1.01%, and predict a SLOS disease incidence of 1/39,215 conceptions. This approach highlights yet another valuable aspect of the exome sequencing databases, to inform clinical and health policy decisions related to genetic counseling, prenatal testing and newborn screening.

Heavy transcription of yeast genes correlates with differential loss of histone H2B relative to H4 and queued RNA polymerases

Eukaryotic chromatin is composed of nucleosomes, which contain nearly two coils of DNA wrapped around a central histone octamer. The octamer contains an H3-H4 tetramer and two H2A-H2B dimers. Gene activation is associated with chromatin disruption: a wider nucleosome-depleted region (NDR) at the promoter and reduced nucleosome occupancy over the coding region. Here, we examine the nature of disrupted chromatin after induction, using MNase-seq to map nucleosomes and subnucleosomes, and a refined high-resolution ChIP-seq method to map H4, H2B and RNA polymerase II (Pol II) genome-wide. Over coding regions, induced genes show a differential loss of H2B relative to H4, which correlates with Pol II density and the appearance of subnucleosomes. After induction, Pol II is surprisingly low at the promoter, but accumulates on the gene and downstream of the termination site, implying that dissociation is very slow. Thus, induction-dependent chromatin disruption reflects both eviction of H2A-H2B dimers and the presence of queued Pol II elongation complexes. We propose that slow Pol II dissociation after transcription is a major factor in chromatin disruption and that it may be of critical importance in gene regulation.