Pseudogenes in the ENCODE regions: consensus annotation, analysis of transcription, and evolution.

Arising from either retrotransposition or genomic duplication of functional genes, pseudogenes are "genomic fossils" valuable for exploring the dynamics and evolution of genes and genomes. Pseudogene identification is an important problem in computational genomics, and is also critical for obtaining an accurate picture of a genome's structure and function. However, no consensus computational scheme for defining and detecting pseudogenes has been developed thus far. As part of the ENCyclopedia Of DNA Elements (ENCODE) project, we have compared several distinct pseudogene annotation strategies and found that different approaches and parameters often resulted in rather distinct sets of pseudogenes. We subsequently developed a consensus approach for annotating pseudogenes (derived from protein coding genes) in the ENCODE regions, resulting in 201 pseudogenes, two-thirds of which originated from retrotransposition. A survey of orthologs for these pseudogenes in 28 vertebrate genomes showed that a significant fraction ( approximately 80%) of the processed pseudogenes are primate-specific sequences, highlighting the increasing retrotransposition activity in primates. Analysis of sequence conservation and variation also demonstrated that most pseudogenes evolve neutrally, and processed pseudogenes appear to have lost their coding potential immediately or soon after their emergence. In order to explore the functional implication of pseudogene prevalence, we have extensively examined the transcriptional activity of the ENCODE pseudogenes. We performed systematic series of pseudogene-specific RACE analyses. These, together with complementary evidence derived from tiling microarrays and high throughput sequencing, demonstrated that at least a fifth of the 201 pseudogenes are transcribed in one or more cell lines or tissues.

The estrogen-responsive element as an inducible enhancer: DNA sequence requirements and conversion to a glucocorticoid-responsive element.

The estrogen-responsive element (ERE) present in the 5'-flanking region of the Xenopus laevis vitellogenin (vit) gene B1 has been characterized by transient expression analysis of chimeric vit-tk-CAT (chloramphenicol acetyltransferase) gene constructs transfected into the human estrogen-responsive MCF-7 cell line. The vit B1 ERE behaves like an inducible enhancer, since it is able to confer estrogen inducibility to the heterologous HSV thymidine kinase (tk) promoter in a relative position- and orientation-independent manner. In this assay, the minimal B1 ERE is 33 bp long and consists of two 13 bp imperfect palindromic elements both of which are required for the enhancer activity. A third imperfect palindromic element is present further upstream within the 5'-flanking region of the gene but is unable to confer hormone responsiveness by itself. Similarly, neither element forming the B1 ERE can alone confer estrogen inducibility to the tk promoter. However, in combinations of two, all three imperfect palindromes can act cooperatively to form a functional ERE. In contrast a single 13 bp perfect palindromic element, GGTCACTGTGACC, such as the one found upstream of the vit gene A2, is itself sufficient to act as a fully active ERE. Single point mutations within this element abolish estrogen inducibility, while a defined combination of two mutations converts this ERE into a glucocorticoid-responsive element.

Sequence homologies in the region preceding the transcription initiation site of the liver estrogen-responsive vitellogenin and apo-VLDLII genes.

In the liver of oviparous vertebrates vitellogenin gene expression is controlled by estrogen. The nucleotide sequence of the 5' flanking region of the Xenopus laevis vitellogenin genes A1, A2, B1 and B2 has been determined. These sequences have been compared to each other and to the equivalent region of the chicken vitellogenin II and apo-VLDLII genes which are also expressed in the liver in response to estrogen. The homology between the 5' flanking region of the Xenopus genes B1 and B2 is higher than between the corresponding regions of the other closely related genes A1 and A2. Four short blocks of sequence homology which are present at equivalent positions in the vitellogenin genes of both Xenopus laevis and chicken are characterized. A short sequence with two-fold rotational symmetry (GGTCANNNTGACC) was found at similar positions upstream of the five vitellogenin genes and is also present in two copies close to the 5' end of the chicken apo-VLDLII gene. The possible functional significance of this sequence, common to liver estrogen-responsive genes, is discussed.

The evolution of lncRNA repertoires and expression patterns in tetrapods.

Only a very small fraction of long noncoding RNAs (lncRNAs) are well characterized. The evolutionary history of lncRNAs can provide insights into their functionality, but the absence of lncRNA annotations in non-model organisms has precluded comparative analyses. Here we present a large-scale evolutionary study of lncRNA repertoires and expression patterns, in 11 tetrapod species. We identify approximately 11,000 primate-specific lncRNAs and 2,500 highly conserved lncRNAs, including approximately 400 genes that are likely to have originated more than 300 million years ago. We find that lncRNAs, in particular ancient ones, are in general actively regulated and may function predominantly in embryonic development. Most lncRNAs evolve rapidly in terms of sequence and expression levels, but tissue specificities are often conserved. We compared expression patterns of homologous lncRNA and protein-coding families across tetrapods to reconstruct an evolutionarily conserved co-expression network. This network suggests potential functions for lncRNAs in fundamental processes such as spermatogenesis and synaptic transmission, but also in more specific mechanisms such as placenta development through microRNA production.

Feasibility of RNA studies on illegitimate transcription for molecular characterization of splicing mutations in the ATP7B gene: a case report.

Approximately 520 Wilson disease-causing mutations in the ATP7B gene have been described to date. In this study we report DNA and RNA analyses carried out for molecular characterization of a consensus sequence splicing mutation found in homozygosity in a Swiss Wilson disease patient. RNA analysis of 1946 +6 T→C in both the peripheral lymphoblasts and liver resulted in the production in the propositus of only an alternative transcript lacking exons 6, 7, and 8 resulting most likely in alterations of cell biochemistry and disease. The patient presents an early form of severe hepatic disease characterized by hepatosplenomegaly, reduced hepatic function, anemia and thrombocytopenia indicating that 1946 +6 T→C is a severe mutation. Since identical results were obtained from both peripheral lymphoblasts and liver they also suggest that RNA studies of illegitimate transcripts can be safely used for molecular characterization of ATP7B splicing mutations, thus improving genetic counseling and diagnosis of Wilson disease. Moreover these studies, contribute to reveal the exact molecular mechanisms producing Wilson disease.

Taxonomic status of Hylomys parvus and Hylomys suillus (Insectivora: Erinaceidae): biochemical and morphological analyses

The genus Hylomys was thought to be represented by a single widespread species. Biochemical and morphometric analyses of several Southeast Asian populations reveal that Sumatra is inhabited by two distinct species, the dwarf gymnure (H. parvus) and the lesser gymnure (H. suillus). The absence of interbreeding between these two groups along with their relatively ancient common origins are documented by several diagnostic loci and a large Nei's genetic distance (D = 0.353 +/- 0.035). The dwarf gymnure has been reported only from the slopes of the Mt. Kerinci volcano in Sumatra, where the species lives at higher elevations than its potential competitor, the lesser gymnure. Other populations of Hylomys from Java, Borneo, and Malaysia are more closely related to the Sumatran sample of H. suillus, but they exhibit strong interpopulational genetic differentiation (D = 0.165 +/- 0.040) that may be accounted for by their isolated montane habitat. In addition, a principal-components analysis based on 16 measurements of the skull clearly separates adult specimens of both species. There is little overlap in the measurements between H. suillus (which is larger) and H. parvus. On Sumatra where both species may be sympatric, the notched space between premaxillary tips, soft texture of the fur, and more delicate skull and dentition are diagnostic of H. parvus.

Genetical and environmental analyses of yield in six biparental soybean crosses

The yearly genetic progress obtained by breeding for increased soybean yield has been considered acceptable worldwide. It is common sense, however, that this progress can be improved further if refined breeding techniques, developed from the knowledge of the genetic mechanisms controlling soybean yield, are used. In this paper, data from four cultivars and/or lines and their derived sets of F2, F3, F7, F8, F9 and F10 generations assayed in 17 environments were analyzed to allow an insight of the genetic control of soybean yield under different environmental conditions. The general picture was of a complex polygene system controlling yield in soybeans. Additive genetic effects predominated although dominance was often found to be significant. Complications such as epistasis, linkage and macro and micro genotype x environment (G x E) interactions were also commonly detected. The overall heritability was 0.29. The relative magnitude of the additive effects and the complicating factors allowed the inference that the latter are not a serious problem to the breeder. The low heritability values and the considerable magnitude of G x E interactions for yield, however, indicated that careful evaluation through experiments designed to allow for the presence of these effects is necessary for successful selection.

Phylogeography and Pleistocene refugia of the adder (Vipera berus) as inferred from mitochondrial DNA sequence data.

In order to contribute to the debate about southern glacial refugia used by temperate species and more northern refugia used by boreal or cold-temperate species, we examined the phylogeography of a widespread snake species (Vipera berus) inhabiting Europe up to the Arctic Circle. The analysis of the mitochondrial DNA (mtDNA) sequence variation in 1043 bp of the cytochrome b gene and in 918 bp of the noncoding control region was performed with phylogenetic approaches. Our results suggest that both the duplicated control region and cytochrome b evolve at a similar rate in this species. Phylogenetic analysis showed that V. berus is divided into three major mitochondrial lineages, probably resulting from an Italian, a Balkan and a Northern (from France to Russia) refugial area in Eastern Europe, near the Carpathian Mountains. In addition, the Northern clade presents an important substructure, suggesting two sequential colonization events in Europe. First, the continent was colonized from the three main refugial areas mentioned above during the Lower-Mid Pleistocene. Second, recolonization of most of Europe most likely originated from several refugia located outside of the Mediterranean peninsulas (Carpathian region, east of the Carpathians, France and possibly Hungary) during the Mid-Late Pleistocene, while populations within the Italian and Balkan Peninsulas fluctuated only slightly in distribution range, with larger lowland populations during glacial times and with refugial mountain populations during interglacials, as in the present time. The phylogeographical structure revealed in our study suggests complex recolonization dynamics of the European continent by V. berus, characterized by latitudinal as well as altitudinal range shifts, driven by both climatic changes and competition with related species.

Topographic ERP analyses: a step-by-step tutorial review

In this tutorial review, we detail both the rationale for as well as the implementation of a set of analyses of surface-recorded event-related potentials (ERPs) that uses the reference-free spatial (i.e. topographic) information available from high-density electrode montages to render statistical information concerning modulations in response strength, latency, and topography both between and within experimental conditions. In these and other ways these topographic analysis methods allow the experimenter to glean additional information and neurophysiologic interpretability beyond what is available from canonical waveform analyses. In this tutorial we present the example of somatosensory evoked potentials (SEPs) in response to stimulation of each hand to illustrate these points. For each step of these analyses, we provide the reader with both a conceptual and mathematical description of how the analysis is carried out, what it yields, and how to interpret its statistical outcome. We show that these topographic analysis methods are intuitive and easy-to-use approaches that can remove much of the guesswork often confronting ERP researchers and also assist in identifying the information contained within high-density ERP datasets

Islands of euchromatin-like sequence and expressed polymorphic sequences within the short arm of human chromosome 21.

The goals of the human genome project did not include sequencing of the heterochromatic regions. We describe here an initial sequence of 1.1 Mb of the short arm of human chromosome 21 (HSA21p), estimated to be 10% of 21p. This region contains extensive euchromatic-like sequence and includes on average one transcript every 100 kb. These transcripts show multiple inter- and intrachromosomal copies, and extensive copy number and sequence variability. The sequencing of the "heterochromatic" regions of the human genome is likely to reveal many additional functional elements and provide important evolutionary information.

Cost-plus-loss analyses of forest inventory strategies based on kNN-assigned reference sample plot data

Abstract

S100A1 deficiency results in prolonged ventricular repolarization in response to sympathetic activation.

S100A1 is a Ca(2+)-binding protein and predominantly expressed in the heart. We have generated a mouse line of S100A1 deficiency by gene trap mutagenesis to investigate the impact of S100A1 ablation on heart function. Electrocardiogram recordings revealed that after beta-adrenergic stimulation S100A1-deficient mice had prolonged QT, QTc and ST intervals and intraventricular conduction disturbances reminiscent of 2 : 1 bundle branch block. In order to identify genes affected by the loss of S100A1, we profiled the mutant and wild type cardiac transcriptomes by gene array analysis. The expression of several genes functioning to the electrical activity of the heart were found to be significantly altered. Although the default prediction would be that mRNA and protein levels are highly correlated, comprehensive immunoblot analyses of salient up- or down-regulated candidate genes of any cellular network revealed no significant changes on protein level. Taken together, we found that S100A1 deficiency results in cardiac repolarization delay and alternating ventricular conduction defects in response to sympathetic activation accompanied by a significantly different transcriptional regulation.

Diffusion tensor echo planar imaging using surface coil transceiver with a semiadiabatic RF pulse sequence at 14.1T.

Diffusion magnetic resonance studies of the brain are typically performed using volume coils. Although in human brain this leads to a near optimal filling factor, studies of rodent brain must contend with the fact that only a fraction of the head volume can be ascribed to the brain. The use of surface coil as transceiver increases Signal-to-Noise Ratio (SNR), reduces radiofrequency power requirements and opens the possibility of parallel transmit schemes, likely to allow efficient acquisition schemes, of critical importance for reducing the long scan times implicated in diffusion tensor imaging. This study demonstrates the implementation of a semiadiabatic echo planar imaging sequence (echo time=40 ms, four interleaves) at 14.1T using a quadrature surface coil as transceiver. It resulted in artifact free images with excellent SNR throughout the brain. Diffusion tensor derived parameters obtained within the rat brain were in excellent agreement with reported values.

Letter to the Editor : on the value of probability for evaluating results of comparative pattern analyses

This letter to the Editor comments on the article Practical relevance of pattern uniqueness in forensic science by P.T. Jayaprakash (Forensic Science International, in press).

Pliocene and Pleistocene diversification and multiple refugia in a Eurasian shrew (Crocidura suaveolens group).

We sequenced 998 base pairs (bp) of mitochondrial DNA cytochrome b and 799 bp of nuclear gene BRCA1 in the Lesser white-toothed shrew (Crocidura suaveolens group) over its geographic range from Portugal to Japan. The aims of the study were to identify the main clades within the group and respective refugia resulting from Pleistocene glaciations. Analyses revealed the Asian lesser white-toothed shrew (C. shantungensis) as the basal clade, followed by a major branch of C. suaveolens, subdivided sensu stricto into six clades, which split-up in the Upper Pliocene and Lower Pleistocene (1.9-0.9 Myr). The largest clade, occurring over a huge range from east Europe to Mongolia, shows evidence of population expansion after a bottleneck. West European clades originated from Iberian and Italo-Balkanic refugia. In the Near East, three clades evolved in an apparent hotspot of refugia (west Turkey, south-west and south-east of the Caucasus). Most clades include specimens of different morphotypes and the validity of many taxa in the C. suaveolens group has to be re-evaluated.