REMiner-II: A tool for rapid identification and configuration of repetitive element arrays from large mammalian chromosomes as a single query

Genes occupy ~3% of the human and mouse genomes whereas repetitive elements (REs), whose biologic functions are largely uncharacterized, constitute greater than 50%. A heterogeneous population of RE arrays (arrangement structures) is formed by combinations of various REs in mammalian genomes. In this study, REMiner-II was refined from the original REMiner for a more efficient identification and configuration of RE arrays from large queries (e.g., human chromosomes) using an unbiased self-alignment protocol. Chromosome-wide RE array profiles for the entire sets of human and mouse chromosomes were obtained using REMiner-II on a personal computer. REMiner-II provides 10 adjustable parameters and three data output modes to accommodate different experimental settings and/or goals. Examination of the human and mouse chromosome data using the REMiner-II viewer revealed species-specific libraries of complexly organized RE arrays. In conclusion, REMiner-II is an efficient tool for chromosome-wide identification and characterization of RE arrays from mammalian genomes.     

Unique profile of ordered arrangements of repetitive elements in the C57BL/6J mouse genome implicating their functional roles

The entirety of all protein coding sequences is reported to represent a small fraction (~2%) of the mouse and human genomes; the vast majority of the rest of the genome is presumed to be repetitive elements (REs). In this study, the C57BL/6J mouse reference genome was subjected to an unbiased RE mining to establish a whole-genome profile of RE occurrence and arrangement. The C57BL/6J mouse genome was fragmented into an initial set of 5,321 units of 0.5 Mb, and surveyed for REs using unbiased self-alignment and dot-matrix protocols. The survey revealed that individual chromosomes had unique profiles of RE arrangement structures, named RE arrays. The RE populations in certain genomic regions were arranged into various forms of complexly organized structures using combinations of direct and/or inverse repeats. Some of these RE arrays spanned stretches of over 2 Mb, which may contribute to the structural configuration of the respective genomic regions. There were substantial differences in RE density among the 21 chromosomes, with chromosome Y being the most densely populated. In addition, the RE array population in the mouse chromosomes X and Y was substantially different from those of the reference human chromosomes. Conversion of the dot-matrix data pertaining to a tandem 13-repeat structure within the Ch7.032 genome unit into a line map of known REs revealed a repeat unit of ~11.3 Kb as a mosaic of six different RE types. The data obtained from this study allowed for a comprehensive RE profiling, including the establishment of a library of RE arrays, of the reference mouse genome. Some of these RE arrays may participate in a spectrum of normal and disease biology that are specific for mice.     

Genome-wide targeted search for human specific and polymorphic L1 integrations

Retroelements (REs) occupy up to 40% of the human genome. Newly integrated REs can change the pattern of expression of pre-existing host genes and therefore might play a significant role in evolution. In particular, human- and primate-specific REs could affect the divergence of the Hominoidea superfamily. A comparative genome-wide analysis of RE sites of integration, neighboring genes, and their regulatory interplay in human and ape genomes would be of help in understanding the impact of REs on evolution and genome regulation. We have developed a technique for the genome-wide comparison of the integrations of transposable elements in genomic DNAs of closely related species. The technique called targeted genome differences analysis (TGDA) is also useful for the detection of deletion/insertion polymorphisms of REs. The technique is based on an enhanced version of subtractive hybridization and does not require preliminary knowledge of the genome sequences under comparison. In this report, we describe its application to the detection and analysis of human specific L1 integrations and their polymorphisms. We obtained a library highly enriched in human-specific L1 insertions and identified 24 such new insertions. Many of these insertions are polymorphic in human populations. The total number of human-specific L1 inserts was estimated to be approximately 4000. The results suggest that TGDA is a universal method that can be successfully used for the detection of evolutionary and polymorphic markers in any closely related genomes.     

Whole-genome experimental identification of insertion/deletion polymorphisms of interspersed repeats by a new general approach

A new experimental technique for genome-wide detection of integration sites of polymorphic retroelements (REs) is described. The technique allows one to reveal the absence of a retroelement in an individual genome provided that this retroelement is present in at least one of several other genomes under comparison. Since quite a number of genomes are compared simultaneously, the search for polymorphic REs insertions is very efficient. The technique includes two whole-genome selective PCR amplifications of sequences flanking REs: one for a particular genome and another one for a mixture of ten different genomes. A subsequent subtractive hybridization of the obtained amplicons with DNA of a particular genome as driver results in isolation of polymorphic insertions. The technique was successfully applied for identification of 41 new polymorphic human AluYa5/Ya8 insertions. Among them, 18 individual Alu elements first sequenced in this work were not found in the available human genome databases. This result suggests that significant part of polymorphic REs were not identified during genome sequencing and remain to be detected and characterized. The proposed method does not depend on preliminary knowledge of evolutionary history of retroelements and can be applied for identification of insertion/deletion polymorphic markers in genomes of different species.     

Interactions of Chromatin Context,Binding Site Sequence Content,and Sequence Evolution in Stress-Induced p53 Occupancy and Transactivation

Cellular stresses activate the tumor suppressor p53 protein leading to selective binding to DNA response elements (REs) and gene transactivation from a large pool of potential p53 REs (p53REs). To elucidate how p53RE sequences and local chromatin context interact to affect p53 binding and gene transactivation, we mapped genome-wide binding localizations of p53 and H3K4me3 in untreated and doxorubicin (DXR)-treated human lymphoblastoid cells. We examined the relationships among p53 occupancy, gene expression, H3K4me3, chromatin accessibility (DNase 1 hypersensitivity, DHS), ENCODE chromatin states, p53RE sequence, and evolutionary conservation. We observed that the inducible expression of p53-regulated genes was associated with the steady-state chromatin status of the cell. Most highly inducible p53-regulated genes were suppressed at baseline and marked by repressive histone modifications or displayed CTCF binding. Comparison of p53RE sequences residing in different chromatin contexts demonstrated that weaker p53REs resided in open promoters, while stronger p53REs were located within enhancers and repressed chromatin. p53 occupancy was strongly correlated with similarity of the target DNA sequences to the p53RE consensus, but surprisingly, inversely correlated with pre-existing nucleosome accessibility (DHS) and evolutionary conservation at the p53RE. Occupancy by p53 of REs that overlapped transposable element (TE) repeats was significantly higher (p<10⁻⁷) and correlated with stronger p53RE sequences (p<10⁻¹¹⁰) relative to nonTE-associated p53REs, particularly for MLT1H, LTR10B, and Mer61 TEs. However, binding at these elements was generally not associated with transactivation of adjacent genes. Occupied p53REs located in L2-like TEs were unique in displaying highly negative PhyloP scores (predicted fast-evolving) and being associated with altered H3K4me3 and DHS levels. These results underscore the systematic interaction between chromatin status and p53RE context in the induced transactivation response. This p53 regulated response appears to have been tuned via evolutionary processes that may have led to repression and/or utilization of p53REs originating from primate-specific transposon elements.     

Effects of varying the position of thyroid hormone response elements within the rat growth hormone promoter: implications for positive and negative regulation by 3,5,3'-triiodothyronine

The thyroid hormone response element (T3RE) of the rat GH (rGH) promoter is located at -188 to -165 relative to the mRNA start site (TSS). Similar sites have been identified in other genes regulated by T3. We have investigated some of these T3REs in positions within the rGH promoter to assess the relative influences of DNA-binding site and position on positive and negative regulation by T3. Synthetic oligonucleotides were used with sequences from the rGH T3RE and proposed negative T3REs (nT3RE) from the rat and human alpha-subunit and rat beta TSH genes. The nT3REs were placed in the background of the wild-type rGH promoter in two positions, at -55 and down-stream of the TSS, with up- and down-mutations of the rGH T3RE. Rat GH T3RE elements were placed 700 basepairs up-stream of a basal rGH promoter and some also at the -55 and TSS positions. Constructions were tested in a transient transfection assay in rat pituitary tumor cells. Two copies of the rGHPAL (palindromic T3RE) placed 700 basepairs up-stream of the rGH promoter conferred 10-fold T3 induction. In the -55 position, the rGHPAL increased T3 induction compared to that in controls, whereas a fragment from the rat and human alpha-subunit gene in the same position reduced induction. Negative T3REs from rat beta TSH and human alpha-subunit reduced T3 induction 50% when placed at the TSS position of a rGH promoter containing an up-mutant T3RE. The T3REPAL placed at the same site increased T3 induction.(ABSTRACT TRUNCATED AT 250 WORDS)     

Genome structure analysis of molluscs revealed whole genome duplication and lineage specific repeat variation

Comparative genome structure analysis allows us to identify novel genes, repetitive sequences and gene duplications. To explore lineage-specific genomic changes of the molluscs that is good model for development of nervous system in invertebrate, we conducted comparative genome structure analyses of three molluscs, pygmy squid, nautilus and scallops using partial genome shotgun sequencing. Most effective elements on the genome structural changes are repetitive elements (REs) causing expansion of genome size and whole genome duplication producing large amount of novel functional genes. Therefore, we investigated variation and proportion of REs and whole genome duplication. We, first, identified variations of REs in the three molluscan genomes by homology-based and de novo RE detection. Proportion of REs were 9.2%, 4.0%, and 3.8% in the pygmy squid, nautilus and scallop, respectively. We, then, estimated genome size of the species as 2.1, 4.2 and 1.8 Gb, respectively, with 2× coverage frequency and DNA sequencing theory. We also performed a gene duplication assay based on coding genes, and found that large-scale duplication events occurred after divergence from the limpet Lottia, an out-group of the three molluscan species. Comparison of all the results suggested that RE expansion did not relate to the increase in genome size of nautilus. Despite close relationships to nautilus, the squid has the largest portion of REs and smaller genome size than nautilus. We also identified lineage-specific RE and gene-family expansions, possibly relate to acquisition of the most complicated eye and brain systems in the three species.     

A new polymorphic retroelement database (PRED) for the human genome

Comparison of primate genomic sequences has demonstrated that the intra-and interspecific genetic variation is provided by retroelements (REs). The human genome contains many thousands of polymorphic RE copies, which are regarded as a promising source of new generation molecular genetic markers. However, the absence of systematized data on the RE number, distribution, genomic context, and abundance in various human populations limits the use of RE insertion polymorphism. We designed the first bilingual (Russian/English) web resource on the known polymorphic REs discovered both by our team and other researchers. The database contains the information about the genomic location of each RE, its position relative to known and predicted genes, abundance in human populations, and other data. Our web portal () allows a search of the database with user-specified parameters. The database makes it possible to most comprehensively analyze the RE distribution in the human genome and to design molecular genetic markers for studies of human genome diversity and biomedical applications.     

Lysosomal acid lipase is the major acid retinyl ester hydrolase in cultured human hepatic stellate cells but not essential for retinyl ester degradation

Vitamin A is stored as retinyl esters (REs) in lipid droplets of hepatic stellate cells (HSCs). To date, two different pathways are known to facilitate the breakdown of REs: (i) Hydrolysis of REs by neutral lipases, and (ii) whole lipid droplet degradation in autolysosomes by acid hydrolysis.In this study, we evaluated the contribution of neutral and acid RE hydrolases to the breakdown of REs in human HSCs. (R)-Bromoenol lactone (R-BEL), inhibitor of adipose triglyceride lipase (ATGL) and patatin-like phospholipase domain-containing 3 (PNPLA3), the hormone-sensitive lipase (HSL) inhibitor 76-0079, as well as the serine-hydrolase inhibitor Orlistat reduced neutral RE hydrolase activity of LX-2 cell-lysates between 20 and 50%. Interestingly, in pulse-chase experiments, R-BEL, 76-0079, as well as Orlistat exerted little to no effect on cellular RE breakdown of LX-2 cells as well as primary human HSCs. In contrast, Lalistat2, a specific lysosomal acid lipase (LAL) inhibitor, virtually blunted acid in vitro RE hydrolase activity of LX-2 cells. Accordingly, HSCs isolated from LAL-deficient mice showed RE accumulation and were virtually devoid of acidic RE hydrolase activity. In pulse-chase experiments however, LAL-deficient HSCs, similar to LX-2 cells and primary human HSCs, were not defective in degrading REs.In summary, results demonstrate that ATGL, PNPLA3, and HSL contribute to neutral RE hydrolysis of human HSCs. LAL is the major acid RE hydrolase in HSCs. Yet, LAL is not limiting for RE degradation under serum-starvation. Together, results suggest that RE breakdown of HSCs is facilitated by (a) so far unknown, non-Orlistat inhibitable RE-hydrolase(s).     

单纯疱疹病毒基因组的限制性核酸内切酶分析

为了获得单纯疱疹病毒(HSV)基因组的限制性核酸内切酶(RE)分析资料和选择合适的RE去研究HSV感染,用11种常用的RE对HSV两个型别的实验室标准毒株的基因组分别作了分析。比较研究的结果表明,BamHI、HpaI和PstI等RE较适于HSV的研究和分型。本研究所采用的小量提取HSV DNA的方法具有快速、简便和实用的特点,值得在HSV感染的诊断、分型和HSV分子流行病学诸研究中推广使用。     

KIAA1363 affects retinyl ester turnover in cultured murine and human hepatic stellate cells

Large quantities of vitamin A are stored as retinyl esters (REs) in specialized liver cells, the hepatic stellate cells (HSCs). To date, the enzymes controlling RE degradation in HSCs are poorly understood. In this study, we identified KIAA1363 (also annotated as arylacetamide deacetylase 1 or neutral cholesterol ester hydrolase 1) as a novel RE hydrolase. We show that KIAA1363 is expressed in the liver, mainly in HSCs, and exhibits RE hydrolase activity at neutral pH. Accordingly, addition of the KIAA1363-specific inhibitor JW480 largely reduced RE hydrolase activity in lysates of cultured murine and human HSCs. Furthermore, cell fractionation experiments and confocal microscopy studies showed that KIAA1363 localizes to the endoplasmic reticulum. We demonstrate that overexpression of KIAA1363 in cells led to lower cellular RE content after a retinol loading period. Conversely, pharmacological inhibition or shRNA-mediated silencing of KIAA1363 expression in cultured murine and human HSCs attenuated RE degradation. Together, our data suggest that KIAA1363 affects vitamin A metabolism of HSCs by hydrolyzing REs at the endoplasmic reticulum, thereby counteracting retinol esterification and RE storage in lipid droplets.     

Chromosomal aberrations induced by the restriction endonucleases EcoR I, Pst I, Sal I and Bam HI in CHO cells

4 widely used cohesive end-producing restriction endonucleases (REs), EcoR I, Pst I, Sal I and Bam HI were tested in CHO cells for their aberration-inducing effects. It was demonstrated that all these REs significantly increased the frequencies of aberrant cells, the aberration frequencies per cell and the aberration frequencies per chromosome. The effects of REs on chromosomal aberrations are similar to ionizing radiation, but more minutes and interchange figures are observed. Polyploid cells are more susceptible to RE treatment, an interesting finding which may be explained by the mechanisms leading to the formation of polyploid cells.     

Retrotransposon distribution and copy number variation in gymnosperm genomes

Retrotransposable elements (REs) and related sequences form a large proportion of conifer genomes. During genome evolution, some RE sequences are degraded or eliminated, but some are evolutionarily stable, and can be identified even in distantly related species. Use of genome sequence information from loblolly pine (Pinus taeda) enables investigation of divergent non-coding RE sequences in other pine and conifer species, including Scots pine (Pinus sylvestris). Non-specific inter-retrotransposon amplified polymorphism technique (IRAP) as well as the amplification polymorphism of 12 RE families were investigated in 80 gymnosperm species. The obtained results were compared with phylogenetic relationships among gymnosperms. Investigation of distantly related gymnosperm species reveals persistent RE sequences, such as IFG and Pineywoods, distributed among a wide range of plant lineages. RE sequence divergence was observed, reflecting periods of inactivity and degradation during speciation of pine lineages, as demonstrated by the delineation of the main pine subgenera. Intraspecific variation of 10 RE copy numbers (CN) between Scots pine individuals ranged from 8.9 to 26.6% of the overall mean estimates. CN analyses were performed in 16 additional gymnosperm species. The analysed pine species contained a similar complement of RE families; however, CN and genome occupation proportions differ. A decrease in RE CN estimates can reflect sequence divergence, associated with independent transposition events. Transposition of some REs can be induced by stress conditions; therefore, even distantly related species inhabiting extreme environments could have similar patterns or distribution of these elements.     

Recycling endosomes can serve as intermediates during transport from the Golgi to the plasma membrane of MDCK cells

The AP-1B clathrin adaptor complex is responsible for the polarized transport of many basolateral membrane proteins in epithelial cells. Localization of AP-1B to recycling endosomes (REs) along with other components (exocyst subunits and Rab8) involved in AP-1B-dependent transport suggested that RE might be an intermediate between the Golgi and the plasma membrane. Although the involvement of endosomes in the secretory pathway has long been suspected, we now present direct evidence using four independent methods that REs play a role in basolateral transport in MDCK cells. Newly synthesized AP-1B-dependent cargo, vesicular stomatitis virus glycoprotein G (VSV-G), was found by video microscopy, immunoelectron microscopy, and cell fractionation to enter transferrin-positive REs within a few minutes after exit from the trans-Golgi network. Although transient, RE entry appears essential because enzymatic inactivation of REs blocked VSV-G delivery to the cell surface. Because an apically targeted VSV-G mutant behaved similarly, these results suggest that REs not only serve as an intermediate but also as a common site for polarized sorting on the endocytic and secretory pathways.     

Hidden heterochromatin: Characterization in the Rodentia species Cricetus cricetus, Peromyscus eremicus (Cricetidae) and Praomys tullbergi (Muridae)

The use of in situ restriction endonuclease (RE) (which cleaves DNA at specific sequences) digestion has proven to be a useful technique in improving the dissection of constitutive heterochromatin (CH), and in the understanding of the CH evolution in different genomes. In the present work we describe in detail the CH of the three Rodentia species, Cricetus cricetus, Peromyscus eremicus (family Cricetidae) and Praomys tullbergi (family Muridae) using a panel of seven REs followed by C-banding. Comparison of the amount, distribution and molecular nature of C-positive heterochromatin revealed molecular heterogeneity in the heterochromatin of the three species. The large number of subclasses of CH identified in Praomys tullbergi chromosomes indicated that the karyotype of this species is the more derived when compared with the other two genomes analyzed, probably originated by a great number of complex chromosomal rearrangements. The high level of sequence heterogeneity identified in the CH of the three genomes suggests the coexistence of different satellite DNA families, or variants of these families in these genomes.