ExoPLOT: Representation of alternative splicing in human tissues and developmental stages with transposed element (TE) involvement.

Advances in RNA high-throughput sequencing and large-scale functional assays yield new insights into the multifaceted activities of transposed elements (TE) and many other previously undiscovered sequence elements. Currently, no tool for easy access, analysis, quantification, and visualization of alternatively spliced exons across multiple tissues or developmental stages is available. Also, analysis pipelines demand computational skills or hardware requirements, which often are hard to meet by wet-lab scientists. We developed ExoPLOT to enable simplified access to massive RNA high throughput sequencing datasets to facilitate the analysis of alternative splicing across many biological samples. To demonstrate the functonality of ExoPLOT, we analyzed the contributon of exonized TEs to human coding sequences (CDS). mRNA splice variants containing the TE-derived exon were quantified and compared to expression levels of TE-free splice variants. For analysis, we utilized 313 human cerebrum, cerebellum, heart, kidney, liver, ovary, and testis transcriptomes, representing various pre- and postnatal developmental stages. ExoPLOT visualizes the relative expression levels of alternative transcripts, e.g., caused by the insertion of new TE-derived exons, across different developmental stages of and among multiple tissues. This tool also provides a unique link between evolution and function during exonization (gain of a new exon) and exaptation (recruitment/co-optation) of a new exon. As input for analysis, we derived a database of 1151 repeat-masked, exonized TEs, representing all prominent families of transposons in the human genome and the collection of human consensus coding sequences (CCDS). ExoPLOT screened preprocessed RNA high-throughput sequencing datasets from seven human tissues to quantify and visualize the dynamics in RNA splicing for these 1151 TE-derived exons during the entire human organ development. In addition, we successfully mapped and analyzed 993 recently described exonized sequences from the human frontal cortex onto these 313 transcriptome libraries. ExoPLOT's approach to preprocessing RNA deep sequencing datasets facilitates alternative splicing analysis and significantly reduces processing times. In addition, ExoPLOT's design allows studying alternative RNA isoforms other than TE-derived in a customized - coordinate-based manner and is available at http://retrogenomics3.uni-muenster.de:3838/exz-plot-d/.

Enhancer occlusion transcripts regulate the activity of human enhancer domains via transcriptional interference: a computational perspective.

Analysis of ENCODE long RNA-Seq and ChIP-seq (Chromatin Immunoprecipitation Sequencing) datasets for HepG2 and HeLa cell lines uncovered 1647 and 1958 transcripts that interfere with transcription factor binding to human enhancer domains. TFBSs (Transcription Factor Binding Sites) intersected by these 'Enhancer Occlusion Transcripts' (EOTrs) displayed significantly lower relative transcription factor (TF) binding affinities compared to TFBSs for the same TF devoid of EOTrs. Expression of most EOTrs was regulated in a cell line specific manner; analysis for the same TFBSs across cell lines, i.e. in the absence or presence of EOTrs, yielded consistently higher relative TF/DNA-binding affinities for TFBSs devoid of EOTrs. Lower activities of EOTr-associated enhancer domains coincided with reduced occupancy levels for histone tail modifications H3K27ac and H3K9ac. Similarly, the analysis of EOTrs with allele-specific expression identified lower activities for alleles associated with EOTrs. ChIA-PET (Chromatin Interaction Analysis by Paired-End Tag Sequencing) and 5C (Carbon Copy Chromosome Conformation Capture) uncovered that enhancer domains associated with EOTrs preferentially interacted with poised gene promoters. Analysis of EOTr regions with GRO-seq (Global run-on) data established the correlation of RNA polymerase pausing and occlusion of TF-binding. Our results implied that EOTr expression regulates human enhancer domains via transcriptional interference.

Transcriptional interference by small transcripts in proximal promoter regions.

Proximal promoter regions (PPR) are heavily transcribed yielding different types of small RNAs. The act of transcription within PPRs might regulate downstream gene expression via transcriptional interference (TI). For analysis, we investigated capped and polyadenylated small RNA transcripts within PPRs of human RefSeq genes in eight different cell lines. Transcripts of our datasets overlapped with experimentally determined transcription factor binding sites (TFBS). For TFBSs intersected by these small RNA transcripts, we established negative correlation of sRNA expression levels and transcription factor (TF) DNA binding affinities; suggesting that the transcripts acted via TI. Accordingly, datasets were designated as TFbiTrs (TF-binding interfering transcripts). Expression of most TFbiTrs was restricted to certain cell lines. This facilitated the analysis of effects related to TFbiTr expression for the same RefSeq genes across cell lines. We consistently uncovered higher relative TF/DNA binding affinities and concomitantly higher expression levels for RefSeq genes in the absence of TFbiTrs. Analysis of corresponding chromatin landscapes supported these results. ChIA-PET revealed the participation of distal enhancers in TFbiTr transcription. Enhancers regulating TFbiTrs, in effect, act as repressors for corresponding downstream RefSeq genes. We demonstrate the significant impact of TI on gene expression using selected small RNA datasets.

Evidence for a Novel Mechanism of Influenza Virus-Induced Type I Interferon Expression by a Defective RNA-Encoded Protein.

Influenza A virus (IAV) defective RNAs are generated as byproducts of error-prone viral RNA replication. They are commonly derived from the larger segments of the viral genome and harbor deletions of various sizes resulting in the generation of replication incompatible viral particles. Furthermore, small subgenomic RNAs are known to be strong inducers of pattern recognition receptor RIG-I-dependent type I interferon (IFN) responses. The present study identifies a novel IAV-induced defective RNA derived from the PB2 segment of A/Thailand/1(KAN-1)/2004 (H5N1). It encodes a 10 kDa protein (PB2∆) sharing the N-terminal amino acid sequence of the parental PB2 protein followed by frame shift after internal deletion. PB2∆ induces the expression of IFNß and IFN-stimulated genes by direct interaction with the cellular adapter protein MAVS, thereby reducing viral replication of IFN-sensitive viruses such as IAV or vesicular stomatitis virus. This induction of IFN is completely independent of the defective RNA itself that usually serves as pathogen-associated pattern and thus does not require the cytoplasmic sensor RIG-I. These data suggest that not only defective RNAs, but also some defective RNA-encoded proteins can act immunostimulatory. In this particular case, the KAN-1-induced defective RNA-encoded protein PB2∆ enhances the overwhelming immune response characteristic for highly pathogenic H5N1 viruses, leading to a more severe phenotype in vivo.

GPAC-genome presence/absence compiler: a web application to comparatively visualize multiple genome-level changes.

Our understanding of genome-wide and comparative sequence information has been broadened considerably by the databases available from the University of California Santa Cruz (UCSC) Genome Bioinformatics Department. In particular, the identification and visualization of genomic sequences, present in some species but absent in others, led to fundamental insights into gene and genome evolution. However, the UCSC tools currently enable one to visualize orthologous genomic loci for a range of species in only a single locus. For large-scale comparative analyses of such presence/absence patterns a multilocus view would be more desirable. Such a tool would enable us to compare thousands of relevant loci simultaneously and to resolve many different questions about, for example, phylogeny, specific aspects of genome and gene evolution, such as the gain or loss of exons and introns, the emergence of novel transposed elements, nonprotein-coding RNAs, and viral genomic particles. Here, we present the first tool to facilitate the parallel analysis of thousands of genomic loci for cross-species presence/absence patterns based on multiway genome alignments. This genome presence/absence compiler uses annotated or other compilations of coordinates of genomic locations and compiles all presence/absence patterns in a flexible, color-coded table linked to the individual UCSC Genome Browser alignments. We provide examples of the versatile information content of such a screening system especially for 7SL-derived transposed elements, nuclear mitochondrial DNA, DNA transposons, and miRNAs in primates (http://www.bioinformatics.uni-muenster.de/tools/gpac, last accessed October 1, 2014).

Exploring Massive Incomplete Lineage Sorting in Arctoids (Laurasiatheria, Carnivora).

Freed from the competition of large raptors, Paleocene carnivores could expand their newly acquired habitats in search of prey. Such changing conditions might have led to their successful distribution and rapid radiation. Today, molecular evolutionary biologists are faced, however, with the consequences of such accelerated adaptive radiations, because they led to sequential speciation more rapidly than phylogenetic markers could be fixed. The repercussions being that current genealogies based on such markers are incongruent with species trees.Our aim was to explore such conflicting phylogenetic zones of evolution during the early arctoid radiation, especially to distinguish diagnostic from misleading phylogenetic signals, and to examine other carnivore-related speciation events. We applied a combination of high-throughput computational strategies to screen carnivore and related genomes in silico for randomly inserted retroposed elements that we then used to identify inconsistent phylogenetic patterns in the Arctoidea group, which is well known for phylogenetic discordances.Our combined retrophylogenomic and in vitro wet lab approach detected hundreds of carnivore-specific insertions, many of them confirming well-established splits or identifying and solving conflicting species distributions. Our systematic genome-wide screens for Long INterspersed Elements detected homoplasy-free markers with insertion-specific truncation points that we used to distinguish phylogenetically informative markers from conflicting signals. The results were independently confirmed by phylogenetic diagnostic Short INterspersed Elements. As statistical analysis ruled out ancestral hybridization, these doubly verified but still conflicting patterns were statistically determined to be genomic remnants from a time of ancestral incomplete lineage sorting that especially accompanied large parts of Arctoidea evolution.

What is an RNA? A top layer for RNA classification.

Every ribonucleic acid begins its cellular life as a transcript. If the transcript or its processing product has a function it should be regarded an RNA. Nonfunctional transcripts, by-products from processing, degradation intermediates, even those originating from (functional) RNAs, and non-functional products of transcriptional gene regulation accomplished via the act of transcription, as well as stochastic (co)transcripts could simply be addressed as transcripts (class 0). The copious functional RNAs (class I), often maturing after one or more processing steps, already are systematized into ever expanding sub-classifications ranging from micro RNAs to rRNAs. Established sub-classifications addressing a wide functional diversity remain unaffected. mRNAs (class II) are distinct from any other RNA by virtue of their potential to be translated into (poly)peptide(s) on ribosomes. We are not proposing a novel RNA classification, but wish to add a basic concept with existing terminology (transcript, RNA, and mRNA) that should serve as an additional framework for carefully delineating RNA function from an avalanche of RNA sequencing data. At the same time, this top level hierarchical model should illuminate important principles of RNA evolution and biology thus heightening our awareness that in biology boundaries and categorizations are typically fuzzy.

Ancestry of the Australian termitivorous numbat.

The Australian numbat, Myrmecobius fasciatus, is the only marsupial that feeds almost exclusively on termites and that has a life following the diurnally restricted and dynamic geographical distribution of termites. The millions of years of this adaptation led to unique morphological and anatomical features, especially basicranial and dental characteristics, that make it difficult to identify a clear phylogenetic affiliation to other marsupials. From DNA sequence analyses, the family Myrmecobiidae is placed within the dasyuromorph marsupials, but the exact position varies from study to study, and support values are mostly rather modest. Here, we report the recovery and analysis of approximately 110,000 quasifossilized traces of mobile element insertions into the genome of a dasyurid marsupial (Tasmanian devil), 25 of which are phylogenetically informative for early dasyuromorphial evolution. Fourteen of these ancient retroposon insertions are shared by the 16 Dasyuromorphia species analyzed, including the numbat, but are absent in the outgroups. An additional 11 other insertions are present in all Dasyuridae but are absent in the numbat. These findings place numbats as the sister group to all living Dasyuridae and show that the investigated Dasyuromorphia, including the Myrmecobiidae, constitutes a monophyletic group that is separated from Peramelemorphia, Notoryctemorphia, and other marsupials.

A universal method for the study of CR1 retroposons in nonmodel bird genomes.

Presence/absence patterns of retroposon insertions at orthologous genomic loci constitute straightforward markers for phylogenetic or population genetic studies. In birds, the convenient identification and utility of these markers has so far been mainly restricted to the lineages leading to model birds (i.e., chicken and zebra finch). We present an easy-to-use, rapid, and cost-effective method for the experimental isolation of chicken repeat 1 (CR1) insertions from virtually any bird genome and potentially nonavian genomes. The application of our method to the little grebe genome yielded insertions belonging to new CR1 subfamilies that are scattered all across the phylogenetic tree of avian CR1s. Furthermore, presence/absence analysis of these insertions provides the first retroposon evidence grouping flamingos + grebes as Mirandornithes and several markers for all subsequent branching events within grebes (Podicipediformes). Five markers appear to be species-specific insertions, including the hitherto first evidence in birds for biallelic CR1 insertions that could be useful in future population genetic studies.

Retroposon insertion patterns of neoavian birds: strong evidence for an extensive incomplete lineage sorting era.

More than 150 Ma, the avian lineage separated from that of other dinosaurs and later diversified into the more than 10,000 species extant today. The early neoavian bird radiations most likely occurred in the late Cretaceous (more than 65 Ma) but left behind few if any molecular signals of their archaic evolutionary past. Retroposed elements, once established in an ancestral population, are highly valuable, virtually homoplasy-free markers of species evolution; after applying stringent orthology criteria, their phylogenetically informative presence/absence patterns are free of random noise and independent of evolutionary rate or nucleotide composition effects. We screened for early neoavian orthologous retroposon insertions and identified six markers with conflicting presence/absence patterns, whereas six additional retroposons established before or after the presumed major neoavian radiation show consistent phylogenetic patterns. The exceptionally frequent conflicting retroposon presence/absence patterns of neoavian orders are strong indicators of an extensive incomplete lineage sorting era, potentially induced by an early rapid successive speciation of ancestral Neoaves.

Tracking marsupial evolution using archaic genomic retroposon insertions.

The Australasian and South American marsupial mammals, such as kangaroos and opossums, are the closest living relatives to placental mammals, having shared a common ancestor around 130 million years ago. The evolutionary relationships among the seven marsupial orders have, however, so far eluded resolution. In particular, the relationships between the four Australasian and three South American marsupial orders have been intensively debated since the South American order Microbiotheria was taxonomically moved into the group Australidelphia. Australidelphia is significantly supported by both molecular and morphological data and comprises the four Australasian marsupial orders and the South American order Microbiotheria, indicating a complex, ancient, biogeographic history of marsupials. However, the exact phylogenetic position of Microbiotheria within Australidelphia has yet to be resolved using either sequence or morphological data analysis. Here, we provide evidence from newly established and virtually homoplasy-free retroposon insertion markers for the basal relationships among marsupial orders. Fifty-three phylogenetically informative markers were retrieved after in silico and experimental screening of approximately 217,000 retroposon-containing loci from opossum and kangaroo. The four Australasian orders share a single origin with Microbiotheria as their closest sister group, supporting a clear divergence between South American and Australasian marsupials. In addition, the new data place the South American opossums (Didelphimorphia) as the first branch of the marsupial tree. The exhaustive computational and experimental evidence provides important insight into the evolution of retroposable elements in the marsupial genome. Placing the retroposon insertion pattern in a paleobiogeographic context indicates a single marsupial migration from South America to Australia. The now firmly established phylogeny can be used to determine the direction of genomic changes and morphological transitions within marsupials.

Retroposon insertions and the chronology of avian sex chromosome evolution.

The vast majority of extant birds possess highly differentiated Z and W sex chromosomes. Nucleotide sequence data from gametologs (homologs on opposite sex chromosomes) suggest that this divergence occurred throughout early bird evolution via stepwise cessation of recombination between identical sex chromosomal regions. Here, we investigated avian sex chromosome differentiation from a novel perspective, using retroposon insertions and random insertions/deletions for the reconstruction of gametologous gene trees. Our data confirm that the CHD1Z/CHD1W genes differentiated in the ancestor of the neognaths, whereas the NIPBLZ/NIPBLW genes diverged in the neoavian ancestor and independently within Galloanserae. The divergence of the ATP5A1Z/ATP5A1W genes in galloanserans occurred independently in the chicken, the screamer, and the ancestor of duck-related birds. In Neoaves, this gene pair differentiated in each of the six sampled representatives, respectively. Additionally, three of our investigated loci can be utilized as universal, easy-to-use independent tools for molecular sexing of Neoaves or Neognathae.

Application of housekeeping npcRNAs for quantitative expression analysis of human transcriptome by real-time PCR.

In recent years the improvements in high-throughput gene expression analysis have led to the discovery of numerous non-protein-coding RNA (npcRNA) molecules. They form an abundant class of untranslated RNAs that have shown to play a crucial role in different biochemical pathways in the cell. Reverse transcription quantitative real-time polymerase chain reaction (RT-qPCR) is an efficient tool to measure RNA abundance and gene expression levels in tiny amounts of material. Despite its sensitivity, the lack of appropriate internal controls necessary for accurate data analysis is a limiting factor for its application in npcRNA research. Common internal controls applied are protein-coding reference genes, also termed "housekeeping" genes (HKGs). However, their expression levels reportedly vary among tissues and different experimental conditions. Moreover, application of HKGs as reference in npcRNA expression analyses is questionable, due to the differences in biogenesis. To address the issue of optimal RT-qPCR normalizers in npcRNA analysis, we performed a systematic evaluation of 18 npcRNAs along with four common HKGs in 20 different human tissues. To determine the most suitable internal control with least expression variance, four evaluation strategies, geNORM, NormFinder, BestKeeper, and the comparative delta C(q) method, were applied. Our data strongly suggest that five npcRNAs, which we term housekeeping RNAs (HKRs), exhibit significantly better constitutive expression levels in 20 different human tissues than common HKGs. Determined HKRs are ideal candidates for RT-qPCR data normalization in human transcriptome analysis, and might also be used as reference genes irrespective of the nature of the genes under investigation.

Rodent evolution: back to the root.

Some 70 Ma, rodents arose along a branch of our own mammalian lineage. Today, about 40% of all mammalian species are rodents and are found in vast numbers on almost every continent. Not only is their proliferation extensive but also the rates of DNA evolution vary significantly among lineages, which has hindered attempts to reconstruct, especially the root of, their evolutionary history. The presence or absence of rare genomic changes, such as short interspersed elements (SINEs), are, however, independent of high molecular substitution rates and provide a powerful, virtually homoplasy-free source for solving such phylogenetic problems. We screened 12 Gb of rodent genomic information using whole-genome three-way alignments, multiple lineage-specific sequences, high-throughput polymerase chain reaction amplifications, and sequencing to reveal 65 phylogenetically informative SINE insertions dispersed over 23 rodent phylogenetic nodes. Eight SINEs and six indels provide significant support for an early association of the Mouse-related and Ctenohystrica (guinea pig and relatives) clades, the Squirrel-related clade being the sister group. This early speciation scenario was also evident in the genomewide distribution pattern of B1-related retroposons, as mouse and guinea pig genomes share six such retroposon subfamilies, containing hundreds of thousands of elements that are clearly absent in the ground squirrel genome. Interestingly, however, two SINE insertions and one diagnostic indel support an association of Ctenohystrica with the Squirrel-related clade. Lineage sorting or a more complex evolutionary scenario that includes an early divergence of the Squirrel-related ancestor and a subsequent hybridization of the latter and the Ctenohystrica lineage best explains such apparently contradictory insertions.

Retroposon insertions provide insights into deep lagomorph evolution.

The homogenous mammalian order Lagomorpha comprises about 80 species in two families, Ochotonidae (pikas) and Leporidae (rabbits and hares). However, the phylogenetic relationships among leporids are controversial. Molecular data, particularly from mitochondrial sequences, give highly homoplasious signals. To resolve the controversy between mitochondrial and nuclear data, we analyzed genomic orthologous retroposon insertion sites, a virtually homoplasy-free marker system. From a differential screen of rabbit genomic data for intronic retroposon insertions of CSINE elements, we polymerase chain reaction-amplified and sequenced 11 retroposons in eight representative lagomorphs. We found three retroposons shared among all lagomorphs but absent in outgroups, four confirmed the monophyly of leporids, and three significantly supported Pronolagus as the sister group to all other leporids. One retroposon supported the monophyly of Lepus. The position of Pronolagus outside of the remaining leporids supports the sequence-based signals of nuclear genes and clearly refutes the misleading signals of mitochondrial genes.

Spatial codes in dendritic BC1 RNA.

BC1 RNA is a dendritic untranslated RNA that has been implicated in local translational control mechanisms in neurons. Prerequisite for a functional role of the RNA in synaptodendritic domains is its targeted delivery along the dendritic extent. We report here that the targeting-competent 5' BC1 domain carries two dendritic targeting codes. One code, specifying somatic export, is located in the medial-basal region of the 5' BC1 stem-loop structure. It is defined by an export-determinant stem-bulge motif. The second code, specifying long-range dendritic delivery, is located in the apical part of the 5' stem-loop domain. This element features a GA kink-turn (KT) motif that is indispensable for distal targeting. It specifically interacts with heterogeneous nuclear ribonucleoprotein A2, a trans-acting targeting factor that has previously been implicated in the transport of MBP mRNA in oligodendrocytes and neurons. Our work suggests that a BC1 KT motif encodes distal targeting via the A2 pathway and that architectural RNA elements, such as KT motifs, may function as spatial codes in neural cells.

On BC1 RNA and the fragile X mental retardation protein.

The fragile X mental retardation protein (FMRP), the functional absence of which causes fragile X syndrome, is an RNA-binding protein that has been implicated in the regulation of local protein synthesis at the synapse. The mechanism of FMRP's interaction with its target mRNAs, however, has remained controversial. In one model, it has been proposed that BC1 RNA, a small non-protein-coding RNA that localizes to synaptodendritic domains, operates as a requisite adaptor by specifically binding to both FMRP and, via direct base-pairing, to FMRP target mRNAs. Other models posit that FMRP interacts with its target mRNAs directly, i.e., in a BC1-independent manner. Here five laboratories independently set out to test the BC1-FMRP model. We report that specific BC1-FMRP interactions could be documented neither in vitro nor in vivo. Interactions between BC1 RNA and FMRP target mRNAs were determined to be of a nonspecific nature. Significantly, the association of FMRP with bona fide target mRNAs was independent of the presence of BC1 RNA in vivo. The combined experimental evidence is discordant with a proposed scenario in which BC1 RNA acts as a bridge between FMRP and its target mRNAs and rather supports a model in which BC1 RNA and FMRP are translational repressors that operate independently.

A novel web-based TinT application and the chronology of the Primate Alu retroposon activity.

BACKGROUND: DNA sequences afford access to the evolutionary pathways of life. Particularly mobile elements that constantly co-evolve in genomes encrypt recent and ancient information of their host's history. In mammals there is an extraordinarily abundant activity of mobile elements that occurs in a dynamic succession of active families, subfamilies, types, and subtypes of retroposed elements. The high frequency of retroposons in mammals implies that, by chance, such elements also insert into each other. While inactive elements are no longer able to retropose, active elements retropose by chance into other active and inactive elements. Thousands of such directional, element-in-element insertions are found in present-day genomes. To help analyze these events, we developed a computational algorithm (Transpositions in Transpositions, or TinT) that examines the different frequencies of nested transpositions and reconstructs the chronological order of retroposon activities. RESULTS: By examining the different frequencies of such nested transpositions, the TinT application reconstructs the chronological order of retroposon activities. We use such activity patterns as a comparative tool to (1) delineate the historical rise and fall of retroposons and their relations to each other, (2) understand the retroposon-induced complexity of recent genomes, and (3) find selective informative homoplasy-free markers of phylogeny. The efficiency of the new application is demonstrated by applying it to dimeric Alu Short INterspersed Elements (SINE) to derive a complete chronology of such elements in primates. CONCLUSION: The user-friendly, web-based TinT interface presented here affords an easy, automated screening for nested transpositions from genome assemblies or trace data, assembles them in a frequency-matrix, and schematically displays their chronological activity history.

Evolutionary history of 7SL RNA-derived SINEs in Supraprimates.

The evolutionary relationships of 7SL RNA-derived SINEs such as the primate Alu or the rodent B1 elements have hitherto been obscure. We established an unambiguous phylogenetic tree for Supraprimates, and derived intraordinal relationships of the 7SL RNA-derived SINEs. As well as new elements in Tupaia and primates, we also found that the purported ancestral fossil Alu monomer was restricted to Primates, and provide here the first description of a potential chimeric promoter box region in SINEs.

Deletion of the MBII-85 snoRNA gene cluster in mice results in postnatal growth retardation.

Prader-Willi syndrome (PWS [MIM 176270]) is a neurogenetic disorder characterized by decreased fetal activity, muscular hypotonia, failure to thrive, short stature, obesity, mental retardation, and hypogonadotropic hypogonadism. It is caused by the loss of function of one or more imprinted, paternally expressed genes on the proximal long arm of chromosome 15. Several potential PWS mouse models involving the orthologous region on chromosome 7C exist. Based on the analysis of deletions in the mouse and gene expression in PWS patients with chromosomal translocations, a critical region (PWScr) for neonatal lethality, failure to thrive, and growth retardation was narrowed to the locus containing a cluster of neuronally expressed MBII-85 small nucleolar RNA (snoRNA) genes. Here, we report the deletion of PWScr. Mice carrying the maternally inherited allele (PWScr(m-/p+)) are indistinguishable from wild-type littermates. All those with the paternally inherited allele (PWScr(m+/p-)) consistently display postnatal growth retardation, with about 15% postnatal lethality in C57BL/6, but not FVB/N crosses. This is the first example in a multicellular organism of genetic deletion of a C/D box snoRNA gene resulting in a pronounced phenotype.