Regulation of mammalian gene expression by retroelements and non-coding tandem repeats

Comparative analysis of the diversity of trinucleotide repeats in bacterial genomes

The human gut is the most favorable niche for microbial populations, and few studies have explored the possibilities of horizontal gene transfer between host and pathogen. Trinucleotide repeat (TNR) expansion in humans can cause more than 40 neurodegenerative diseases. Further, TNRs are a type of microsatellite that resides on coding regions can contribute to the synthesis of homopolymeric amino acids. Hence, the present study aims to estimate the occurrence and diversity of TNRs in bacterial genomes available in the NCBI Genome database. Genome-wide analyses revealed that several bacterial genomes contain different types of uninterrupted TNRs. It was found that TNRs are abundant in the genomes of Alcaligenes faecalis, Mycoplasma gallisepticum, Mycoplasma genitalium, Sorangium cellulosum, and Thermus thermophilus. Interestingly, the genome of Bacillus thuringiensis strain YBT-1518 contained 169 uninterrupted ATT repeats. The genome of Leclercia adecarboxylata had 46 uninterrupted CAG repeats, which potentially translate into polyglutamine. In some instances, the TNRs were present in genes that potentially encode essential functions. Similar occurrences in human genes are known to cause genetic disorders. Further analysis of the occurrence of TNRs in bacterial genomes is likely to provide a better understanding of mismatch repair, genetic disorders, host-pathogen interaction, and homopolymeric amino acids.

Drosophila Satellite Repeats at the Intersection of Chromatin, Gene Regulation and Evolution

Satellite repeats make up a large fraction of the genomes of many higher eukaryotes. Until recently these sequences were viewed as molecular parasites with few functions. Drosophila melanogaster and related species have a wealth of diverse satellite repeats. Comparative studies of Drosophilids have been instrumental in understanding how these rapidly evolving sequences change and move. Remarkably, satellite repeats have been found to modulate gene expression and mediate genetic conflicts between chromosomes and between closely related fly species. This suggests that satellites play a key role in speciation. We have taken advantage of the depth of research on satellite repeats in flies to review the known functions of these sequences and consider their central role in evolution and gene expression.

The wide distribution and horizontal transfers of beta satellite DNA in eukaryotes

Beta satellite DNA (satDNA), also known as Sau3A sequences, are repetitive DNA sequences reported in human and primate genomes. It is previously thought that beta satDNAs originated in old world monkeys and bursted in great apes. In this study, we searched 7821 genome assemblies of 3767 eukaryotic species and found that beta satDNAs are widely distributed across eukaryotes. The four major branches of eukaryotes, animals, fungi, plants and Harosa/SAR, all have multiple clades containing beta satDNAs. These results were also confirmed by searching whole genome sequencing data (SRA) and PCR assay. Beta satDNA sequences were found in all the primate clades, as well as in Dermoptera and Scandentia, indicating that the beta satDNAs in primates might originate in the common ancestor of Primatomorpha or Euarchonta. In contrast, the widely patchy distribution of beta satDNAs across eukaryotes presents a typical scenario of multiple horizontal transfers.

Tissue and cancer-specific expression of DIEXF is epigenetically mediated by an Alu repeat

Alu repeats constitute a major fraction of human genome and for a small subset of them a role in gene regulation has been described. The number of studies focused on the functional characterization of particular Alu elements is very limited. Most Alu elements are DNA methylated and then assumed to lie in repressed chromatin domains. We hypothesize that Alu elements with low or variable DNA methylation are candidates for a functional role. In a genome-wide study in normal and cancer tissues, we pinpointed an Alu repeat (AluSq2) with differential methylation located upstream of the promoter region of the DIEXF gene. DIEXF encodes a highly conserved factor essential for the development of zebrafish digestive tract. To characterize the contribution of the Alu element to the regulation of DIEXF we analysed the epigenetic landscapes of the gene promoter and flanking regions in different cell types and cancers. Alternate epigenetic profiles (DNA methylation and histone modifications) of the AluSq2 element were associated with DIEXF transcript diversity as well as protein levels, while the epigenetic profile of the CpG island associated with the DIEXF promoter remained unchanged. These results suggest that AluSq2 might directly contribute to the regulation of DIEXF transcription and protein expression. Moreover, AluSq2 was DNA hypomethylated in different cancer types, pointing out its putative contribution to DIEXF alteration in cancer and its potential as tumoural biomarker.

Identification of self-regulatory network motifs in reverse engineering gene regulatory networks using microarray gene expression data

Gene Regulatory Networks (GRNs) are reconstructed from the microarray gene expression data through diversified computational approaches. This process ensues in symmetric and diagonal interaction of gene pairs that cannot be modelled as direct activation, inhibition, and self-regulatory interactions. The values of gene co-expressions could help in identifying co-regulations among them. The proposed approach aims at computing the differences in variances of co-expressed genes rather than computing differences in values of mean expressions across experimental conditions. It adopts multivariate co-variances using principal component analysis (PCA) to predict an asymmetric and non-diagonal gene interaction matrix, to select only those gene pair interactions that exhibit the maximum variances in gene regulatory expressions. The asymmetric gene regulatory interactions help in identifying the controlling regulatory agents, thus lowering the false positive rate by minimizing the connections between previously unlinked network components. The experimental results on real as well as in silico datasets including time-series RTX therapy, Arabidopsis thaliana, DREAM-3, and DREAM-8 datasets, in comparison with existing state-of-the-art approaches demonstrated the enhanced performance of the proposed approach for predicting positive and negative feedback loops and self-regulatory interactions. The generated GRNs hold the potential in determining the real nature of gene pair regulatory interactions.

Regular Higher Order Repeat Structures in Beetle Tribolium castaneum Genome

Higher order repeats (HORs) containing tandems of primary and secondary repeat units (head-to-tail “tandem within tandem pattern”), referred to as regular HORs, are typical for primate alpha satellite DNAs and most pronounced in human genome. Regular HORs are known to be a result of recent evolutionary processes. In non-primate genomes mostly so called complex HORs have been found, without head to tail tandem of primary repeat units. In beetle Tribolium castaneum, considered as a model case for genome studies, large tandem repeats have been identified, but no HORs have been reported. Here, using our novel robust repeat finding algorithm Global Repeat Map, we discover two regular and six complex HORs in T. castaneum. In organizational pattern, the integrity and homogeneity of regular HORs in T. castaneum resemble human regular HORs (with T. castaneum monomers different from human alpha satellite monomers), involving a wider range of monomer lengths than in human HORs. Similar regular higher order repeat structures have previously not been found in insects. Some of these novel HORs in T. castaneum appear as most regular among known HORs in non-primate genomes, although with substantial riddling. This is intriguing, in particular from the point of view of role of non-coding repeats in modulation of gene expression.

Dispersion Profiles and Gene Associations of Repetitive DNAs in the Euchromatin of the Beetle Tribolium castaneum

Satellite DNAs are tandemly repeated sequences clustered within heterochromatin. However, in some cases, such as the major TCAST1 satellite DNA from the beetle Tribolium castaneum, they are found partially dispersed within euchromatin. Such organization together with transcriptional activity enables TCAST1 to modulate the activity of neighboring genes. In order to explore if other T. castaneum repetitive families have features that could provide them with a possible gene-modulatory role, we compare here the structure, organization, dispersion profiles, and transcription activity of 10 distinct TCAST repetitive families including TCAST1. The genome organization of TCAST families exhibit either satellite-like or transposon-like characteristics. In addition to heterochromatin localization, bioinformatic searches of the assembled genome have revealed dispersion of all families within euchromatin, preferentially in the form of single repeats. Dispersed TCAST repeats are mutually correlated in distribution and are grouped in distinct regions of euchromatin. The repeats are associated with genes, are enriched in introns relative to intergenic regions, and very rarely overlap exons. In spite of the different mechanisms of repeat proliferation, such as transposition and homologous recombination, all TCAST families share a similar frequency of spreading as well as dispersion and gene association profiles. Additionally, TCAST families are transcribed and their transcription is significantly activated by heat stress. A possibility that such common features of TCAST families might be related to their potential gene-modulatory role is discussed.

Genomic landscape of CpG rich elements in human

Background

The studies on CpG islands (CGI) and Alu elements functions, evolution, and distribution in the genome started since the discovery in nineteen eighties (1981, 1986, correspondingly). Their highly skewed genome wide distribution implies the non-random retrotransposition pattern. Besides CGIs in gene promoters, CGIs clusters were observed in the homeobox gene regions and in the macrosatellites, but the whole picture of their distribution specifics was not grasped. Attempts to identify any causative features upon their (genome wide) distribution, such as the DNA context mediated preferred insertion sites of Alu repeats, have been made to ascribe their clusters location.

Methods

Recent emergence of high resolution 3D map of human genome allowed segregating the genome into the large scale chromatin domains of naturally observable nuclear subcompartments, or Topologically Associated Domains (TADs), designated by spatial chromatin distribution. We utilized the chromatin map to elucidate relations between large scale chromatin state and CpG rich elements landscape.

In the course of analysis it was confirmed that genes, Alu and CGI clusters maintain obvious, albeit different in strength, preference for open chromatin. For the first time it was clearly shown that the clusters density of the Alu and CGIs monotonically depend on the chromatin accessibility rate. In particular, the highest density of these elements is found in A1 euchromatin regions characterized by a high density of small length genes replicating in the early S-phase. It implies that these elements mediate (CGIs) or are a side element (Alus) of chromatin accessibility.

Results

We elucidated that both methylated and non-methylated CGIs display the affinity to chromatin accessibility. As a part of comparative genomics section, we elucidated that the dog’s genome non-canonical structure, outstanding in mammals for its high CGIs abundance compared to gene number, is explained by the presence of dense tandem CGI extended hotspots (500 kb on average) in subtelomeric and pericentromeric regions with highly skewed CG content, and not by CGIs global distribution pattern shift.

Conclusions

The study underlines the close association of CG-rich elements distribution with the newly introduced large scale chromatin state map, proposing a refined standpoint on interrelation of aforementioned genome elements and the chromatin state. To our expertise, the TAD-associated partition model employed in the study is likely the most substantial one regarding CpG rich clusters distribution among the whole genome chromatin/isochores maps available.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0864-0) contains supplementary material, which is available to authorized users.

SBMDb: first whole genome putative microsatellite DNA marker database of sugarbeet for bioenergy and industrial applications

DNA marker plays important role as valuable tools to increase crop productivity by finding plausible answers to genetic variations and linking the Quantitative Trait Loci (QTL) of beneficial trait. Prior approaches in development of Short Tandem Repeats (STR) markers were time consuming and inefficient. Recent methods invoking the development of STR markers using whole genomic or transcriptomics data has gained wide importance with immense potential in developing breeding and cultivator improvement approaches. Availability of whole genome sequences and in silico approaches has revolutionized bulk marker discovery. We report world's first sugarbeet whole genome marker discovery having 145 K markers along with 5 K functional domain markers unified in common platform using MySQL, Apache and PHP in SBMDb. Embedded markers and corresponding location information can be selected for desired chromosome, location/interval and primers can be generated using Primer3 core, integrated at backend. Our analyses revealed abundance of 'mono' repeat (76.82%) over 'di' repeats (13.68%). Highest density (671.05 markers/Mb) was found in chromosome 1 and lowest density (341.27 markers/Mb) in chromosome 6. Current investigation of sugarbeet genome marker density has direct implications in increasing mapping marker density. This will enable present linkage map having marker distance of ∼2 cM, i.e. from 200 to 2.6 Kb, thus facilitating QTL/gene mapping. We also report e-PCR-based detection of 2027 polymorphic markers in panel of five genotypes. These markers can be used for DUS test of variety identification and MAS/GAS in variety improvement program. The present database presents wide source of potential markers for developing and implementing new approaches for molecular breeding required to accelerate industrious use of this crop, especially for sugar, health care products, medicines and color dye. Identified markers will also help in improvement of bioenergy trait of bioethanol and biogas production along with reaping advantage of crop efficiency in terms of low water and carbon footprint especially in era of climate change. Database URL: http://webapp.cabgrid.res.in/sbmdb/.

Transcriptional regulation of gilthead seabream bone morphogenetic protein (BMP) 2 gene by bone- and cartilage-related transcription factors

Bone morphogenetic protein (BMP) 2 belongs to the transforming growth factor β (TGFβ) superfamily of cytokines and growth factors. While it plays important roles in embryo morphogenesis and organogenesis, BMP2 is also critical to bone and cartilage formation. Protein structure and function have been remarkably conserved throughout evolution and BMP2 transcription has been proposed to be tightly regulated, although few data is available. In this work we report the cloning and functional analysis of gilthead seabream BMP2 promoter. As in other vertebrates, seabream BMP2 gene has a 5′ non-coding exon, a feature already present in DPP gene, the fruit fly ortholog of vertebrate BMP2 gene, and maintained throughout evolution. In silico analysis of seabream BMP2 promoter revealed several binding sites for bone and cartilage related transcription factors (TFs) and their functionality was evaluated using promoter-luciferase constructions and TF-expressing vectors. Runt-related transcription factor 3 (RUNX3) was shown to negatively regulate BMP2 transcription and combination with the core binding factor β (CBFβ) further reduced transcriptional activity of the promoter. Although to a lesser extent, myocyte enhancer factor 2C (MEF2C) had also a negative effect on the regulation of BMP2 gene transcription, when associated with SRY (sex determining region Y)-box 9 (SOX9b). Finally, v-ets avian erythroblastosis virus E26 oncogene homolog 1 (ETS1) was able to slightly enhance BMP2 transcription. Data reported here provides new insights toward the better understanding of the transcriptional regulation of BMP2 gene in a bone and cartilage context.

Satellite DNA Modulates Gene Expression in the Beetle Tribolium castaneum after Heat Stress

Non-coding repetitive DNAs have been proposed to perform a gene regulatory role, however for tandemly repeated satellite DNA no such role was defined until now. Here we provide the first evidence for a role of satellite DNA in the modulation of gene expression under specific environmental conditions. The major satellite DNA TCAST1 in the beetle Tribolium castaneum is preferentially located within pericentromeric heterochromatin but is also dispersed as single repeats or short arrays in the vicinity of protein-coding genes within euchromatin. Our results show enhanced suppression of activity of TCAST1-associated genes and slower recovery of their activity after long-term heat stress relative to the same genes without associated TCAST1 satellite DNA elements. The level of gene suppression is not influenced by the distance of TCAST1 elements from the associated genes up to 40 kb from the genes' transcription start sites, but it does depend on the copy number of TCAST1 repeats within an element, being stronger for the higher number of copies. The enhanced gene suppression correlates with the enrichment of the repressive histone marks H3K9me2/3 at dispersed TCAST1 elements and their flanking regions as well as with increased expression of TCAST1 satellite DNA. The results reveal transient, RNAi based heterochromatin formation at dispersed TCAST1 repeats and their proximal regions as a mechanism responsible for enhanced silencing of TCAST1-associated genes. Differences in the pattern of distribution of TCAST1 elements contribute to gene expression diversity among T. castaneum strains after long-term heat stress and might have an impact on adaptation to different environmental conditions.

Differential regulation of genes by retrotransposons in rice promoters

Rice genome harbors genes and promoters with retrotransposon insertions. There is very little information about their function. The effect of retrotransposon insertions in four rice promoter regions on gene regulation, was investigated using promoter-reporter gene constructs with and without retrotransposons. Differences in expression levels of gus and egfp reporter genes in forward orientation and rfp in reverse orientation were evaluated in rice plants with transient expression employing quantitative RT-PCR analysis, histochemical GUS staining, and eGFP and RFP fluorescent microscopy. The presence of SINE in the promoter 1 (P1) resulted in higher expression levels of the reporter genes, whereas the presence of LINE in P2 or gypsy LTR retrotransposon in P3 reduced expression of the reporter genes. Furthermore, the SINE in P1 acts as an enhancer in contrast with the LINE in P2 and the gypsy LTR retrotransposon in P3 which act as silencers. CTAA and CGG motifs in these retrotransposons are the likely candidates for the downregulation compared to TCTT motif (SINE) which is a candidate for the upregulation of gene expression. The effect of retrotransposons on gene regulation correlated with the earlier investigation of conservation patterns of these four retrotransposon insertions in several rice accessions implying their evolutionary significance.

Astrocytic abnormalities and global DNA methylation patterns in depression and suicide

Astrocytes are glial cells specific to the central nervous system and involved in numerous brain functions, including regulation of synaptic transmission and of immune reactions. There is mounting evidence suggesting astrocytic dysfunction in psychopathologies such as major depression, however, little is known about the underlying etiological mechanisms. Here we report a two-stage study investigating genome-wide DNA methylation associated with astrocytic markers in depressive psychopathology. We first characterized prefrontal cortex samples from 121 individuals (76 who died during a depressive episode and 45 healthy controls) for the astrocytic markers GFAP, ALDH1L1, SOX9, GLUL, SCL1A3, GJA1 and GJB6. A subset of 22 cases with consistently downregulated astrocytic markers was then compared with 17 matched controls using methylation binding domain-2 (MBD2) sequencing followed by validation with high-resolution melting and bisulfite Sanger sequencing. With these data, we generated a genome-wide methylation map unique to altered astrocyte-associated depressive psychopathology. The map revealed differentially methylated regions (DMRs) between cases and controls, the majority of which displayed reduced methylation levels in cases. Among intragenic DMRs, those found in GRIK2 (glutamate receptor, ionotropic kainate 2) and BEGAIN (brain-enriched guanylate kinase-associated protein) were most significant and also showed significant correlations with gene expression. Cell-sorted fractions were investigated and demonstrated an important non-neuronal contribution of methylation status in BEGAIN. Functional cell assays revealed promoter and enhancer-like properties in this region that were markedly decreased by methylation. Furthermore, a large number of our DMRs overlapped known Encyclopedia of DNA elements (ENCODE)-identified regulatory elements. Taken together, our data indicate significant differences in the methylation patterns specific to astrocytic dysfunction associated with depressive psychopathology, providing a potential framework for better understanding this disease phenotype.

Identification, cross-taxon transferability and application of full-length cDNA SSR markers in Phyllostachys pubescens

Current databases of Phyllostachys pubescens full-length cDNAs (FL-cDNAs) provide a rich source of sequences for the development of potential FL-cDNA simple sequence repeat (SSR) markers. We screened 10,608 P. pubescens cDNAs, discovering 1614 SSRs in 1382 SSR-containing FL-cDNAs. The SSRs were more abundant within transposable elements (TEs) than expressed sequence tags (ESTs) and genome survey sequences (GSSs), and specific dinucleotide repeats tended to associate with particular TE families: (TA)n with En/Spm and (CT)n with Mutator. A selected panel of 100 FL-cDNAs containing type I SSRs yielded 68 functional SSR markers with an average polymorphism information content (PIC) value of 0.12, among which 22 loci contained polymorphisms. These markers became less transferrable (83.1% → 69.9% → 49.3%) but more polymorphic (79.4% → 92.3% → 92.8%) with increasing phylogenetic distance (intra-genus → intra-subtribe → intra-family). Transferability and polymorphism also depended on the location of the marker, with those located in the coding region being more transferrable (69.1%) and less polymorphic (89.4%) than those in the 5'-UTR (63.4% transferable, 90.7% polymorphic) and the 3'-UTR (61.8% transferable, 91.4% polymorphic). As proof of principle, we were able to use our FL-cDNA SSR markers to identify the parental stocks in interspecific hybrids of bamboo within and beyond P. pubescens, and estimate the outcrossing rate for P. pubescens. Our research should facilitate molecular breeding in bamboo species where original genetic markers are scarce.

Sequence variability of the MspI satellite DNA family of the pinewood nematode Bursaphelenchus xylophilus at different geographic scales

Tandemly repeated sequences known as satellite DNA (satDNA) generally exhibit complex evolutionary patterns of concerted evolution in which mutations are homogenized and fixed in a stochastic process of molecular drive. Here, the nucleotidic variability of the MspI satDNA family of the pinewood nematode Bursaphelenchus xylophilus is analyzed in order to understand the evolutionary dynamics of satDNA at the intraspecific level. A total of 425 MspI monomer units, either PCR-amplified from isolates of local (Peninsula of Setúbal, Portugal) or worldwide origin, or retrieved from the B. xylophilus genome sequence, were characterized and compared. Whatever their origin, sliding window analysis of sequence variability patterns among monomers revealed low, moderate and highly variant domains, indicating that variable levels of evolutionary constraint may act upon the entire monomers. The phylogenetic inference based on the different sets of MspI satDNA family for this species shows a broad polymorphism of the individual monomers, which were distributed into four main clusters. However, such clustering appeared independent from the geographic origin of the nematodes, and could not discriminate isolates or groups of geographically close isolates. Rather, the formation of different phylogenetic groups within this satDNA family suggests an a priori embodying of a set of diverging repeats from a common ancestor satDNA library, which have been differently amplified along the evolutionary pathway of this species. The present work improves knowledge on the evolutionary dynamics of satDNA at the intraspecific level, and provides new information on satDNA sequence variability among natural populations sampled at a local geographic scale.

Transposable element 'roo' attaches to nuclear matrix of the Drosophila melanogaster

The genome of eukaryotes is organized into structural units of chromatin loops. This higher order organization is supported by a nuclear skeleton called the nuclear matrix. The genomic DNA associated with the nuclear matrix is called the matrix associated region (MAR). Only a few genome-wide screens have been attempted, although many studies have characterized locusspecific MAR DNA sequences. In this study, a MAR DNA library was prepared from the Drosophila melanogaster Meigen (Diptera: Drosophilidae) genome. One of the sequences identified as a MAR was from a long terminal repeat region of 'roo' retrotransposon (roo MAR). Sequence analysis of roo MAR showed its distribution across the D. melanogaster genome. roo MAR also showed high sequence similarity with a previously identified MAR in Drosophila, namely the 'gypsy' retrotransposon. Analysis of the genes flanking roo MAR insertions in the Drosophila genome showed that genes were co-ordinately expressed. The results from the present study in D. melanogaster suggest this sequence plays an important role in genome organization and function. The findings point to an evolutionary role of retrotransposons in shaping the genomic architecture of eukaryotes.

Transposable elements: an abundant and natural source of regulatory sequences for host genes

The fact that transposable elements (TEs) can influence host gene expression was first recognized more than 50 years ago. However, since that time, TEs have been widely regarded as harmful genetic parasites-selfish elements that are rarely co-opted by the genome to serve a beneficial role. Here, we survey recent findings that relate to TE impact on host genes and remind the reader that TEs, in contrast to other noncoding parts of the genome, are uniquely suited to gene regulatory functions. We review recent studies that demonstrate the role of TEs in establishing and rewiring gene regulatory networks and discuss the overall ubiquity of exaptation. We suggest that although individuals within a population can be harmed by the deleterious effects of new TE insertions, the presence of TE sequences in a genome is of overall benefit to the population.

Evolutionary history of linked D4Z4 and Beta satellite clusters at the FSHD locus (4q35)

We performed a detailed genomic investigation of the chimpanzee locus syntenic to human chromosome 4q35.2, associated to the facioscapulohumeral dystrophy. Two contigs of approximately 150 kb and 200 kb were derived from PTR chromosomes 4q35 and 3p12, respectively: both regions showed a very similar sequence organization, including D4Z4 and Beta satellite linked clusters. Starting from these findings, we derived a hypothetical evolutionary history of human 4q35, 10q26 and 3p12 chromosome regions focusing on the D4Z4–Beta satellite linked organization. The D4Z4 unit showed an open reading frame (DUX4) at both PTR 4q35 and 3p12 regions; furthermore some subregions of the Beta satellite unit showed a high degree of conservation between chimpanzee and humans. In conclusion, this paper provides evidence that at the 4q subtelomere the linkage between D4Z4 and Beta satellite arrays is a feature that appeared late during evolution and is conserved between chimpanzee and humans.

Large-scale transcriptome analysis of retroelements in the migratory locust, Locusta migratoria

Background

Retroelements can successfully colonize eukaryotic genome through RNA-mediated transposition, and are considered to be some of the major mediators of genome size. The migratory locust Locusta migratoria is an insect with a large genome size, and its genome is probably subject to the proliferation of retroelements. An analysis of deep-sequencing transcriptome data will elucidate the structure, diversity and expression characteristics of retroelements.

Results

We performed a de novo assembly from deep sequencing RNA-seq data and identified 105 retroelements in the locust transcriptome. Phylogenetic analysis of reverse transcriptase sequences revealed 1 copia, 1 BEL, 8 gypsy and 23 non-long terminal repeat (LTR) retroelements in the locust transcriptome. A novel approach was developed to identify full-length LTR retroelements. A total of 5 full-length LTR retroelements and 2 full-length non-LTR retroelements that contained complete structures for retrotransposition were identified. Structural analysis indicated that all these retroelements may have been activated or deprived of retrotransposition activities very recently. Expression profiling analysis revealed that the retroelements exhibited a unique expression pattern at the egg stage and showed differential expression profiles between the solitarious and gregarious phases at the fifth instar and adult stage.

Conclusion

We hereby present the first de novo transcriptome analysis of retroelements in a species whose genome is not available. This work contributes to a comprehensive understanding of the landscape of retroelements in the locust transcriptome. More importantly, the results reveal that non-LTR retroelements are abundant and diverse in the locust transcriptome.

Introduction to genome biology: features, processes, and structures

Genomic analyses increasingly make use of sophisticated statistical and computational approaches in investigations of genomic function and evolution. Scientists implementing and developing these approaches are often computational scientists, physicists, or mathematicians. This article aims to provide a compact overview of genome biology for these scientists. Thus, the article focuses on providing biological context to the genomic features, processes, and structures analysed by these approaches. Topics covered include (1) differences between eukaryotic and prokaryotic cells; (2) the physical structure of genomes and chromatin; (3) different categories of genomic regions, including those serving as templates for RNA and protein synthesis, regulatory regions, repetitive regions, and "architectural" or "organisational" regions, such as centromeres and telomeres; (4) the cell cycle; (5) an overview of transcription, translation, and protein structure; and (6) a glossary of relevant terms.

Characterization and potential functional significance of human-chimpanzee large INDEL variation

BACKGROUND

Although humans and chimpanzees have accumulated significant differences in a number of phenotypic traits since diverging from a common ancestor about six million years ago, their genomes are more than 98.5% identical at protein-coding loci. This modest degree of nucleotide divergence is not sufficient to explain the extensive phenotypic differences between the two species. It has been hypothesized that the genetic basis of the phenotypic differences lies at the level of gene regulation and is associated with the extensive insertion and deletion (INDEL) variation between the two species. To test the hypothesis that large INDELs (80 to 12,000 bp) may have contributed significantly to differences in gene regulation between the two species, we categorized human-chimpanzee INDEL variation mapping in or around genes and determined whether this variation is significantly correlated with previously determined differences in gene expression.

RESULTS

Extensive, large INDEL variation exists between the human and chimpanzee genomes. This variation is primarily attributable to retrotransposon insertions within the human lineage. There is a significant correlation between differences in gene expression and large human-chimpanzee INDEL variation mapping in genes or in proximity to them.

CONCLUSIONS

The results presented herein are consistent with the hypothesis that large INDELs, particularly those associated with retrotransposons, have played a significant role in human-chimpanzee regulatory evolution.

HLA-DRB1 may be antagonistically regulated by the coordinately evolved promoter and 3'-UTR under stabilizing selection

HLA-DRB1 is the most polymorphic MHC (major histocompatibility complex) class II gene in human, and plays a crucial role in the development and function of the immune system. Extensive polymorphisms exist in the promoter and 3'-UTR of HLA-DRB1, especially a LTR (Long terminal repeat) element in the promoter, which may be involved in the expression regulation. However, it remains unknown how the polymorphisms in the whole promoter region and 3'-UTR to regulate the gene expression. In this study, we investigated the extensive polymorphisms in the HLA-DRB1 promoter and 3'-UTR, and how these polymorphisms affect the gene expression in both independent and jointly manners. It was observed that most of the haplotypes in the DRB1 promoter and 3'-UTR were clustered into 4 conserved lineages (H1, H2, H3 and H4), and showed high linkage disequilibrium. Compared with H1 and H2 lineage, a LTR element in the promoter of H3 and H4 lineage significantly suppressed the promoter activity, whereas the activity of the linked 3'-UTR increased, leading to no apparent difference in the final expression product between H1/H2 and H3/H4 lineage. Nevertheless, compared with the plasmid with a promoter and 3'-UTR from the same lineage, the recombinant plasmid with a promoter from H2 and a 3'-UTR from H3 showed about double fold increased luciferase activity, Conversely, the recombinant plasmid with a promoter from H3 and a 3'-UTR from H2 resulted in about 2-fold decreased luciferase activity. These results indicate that the promoter and 3'-UTR of HLA-DRB1 may antagonistically regulate the gene expression, which may be subjected to stabilizing selection. These findings may provide a novel insight into the mechanisms of the diseases associated with HLA-DRB1 genes.

Genomic relationship between SINE retrotransposons, Pol III-Pol II transcription, and chromatin organization: the journey from junk to jewel

A typical eukaryotic genome harbors a rich variety of repetitive elements. The most abundant are retrotransposons, mobile retroelements that utilize reverse transcriptase and an RNA intermediate to relocate to a new location within the cellular genomes. A vast majority of the repetitive mammalian genome content has originated from the retrotransposition of SINE (100-300 bp short interspersed nuclear elements that are derived from the structural 7SL RNA or tRNA), LINE (7kb long interspersed nuclear element), and LTR (2-3 kb long terminal repeats) transposable element superfamilies. Broadly labeled as "evolutionary junkyard" or "fossils", this enigmatic "dark matter" of the genome possesses many yet to be discovered properties.

Transcriptional repression of repeat-derived transcripts correlates with histone hypoacetylation at repetitive DNA elements in aged mice brain

In order to better characterize epigenetic alterations at repetitive DNA elements with aging, DNA methylation and histone marks at various repeat classes were investigated. Repetitive DNA elements were hypermethylated in the brains of old mice. Histone hypoacetylation and altered histone trimethylation at repetitive sequences were detected in brain tissues during aging. The expression of repeat-derived transcripts (RDTs) was then measured to explore any correlations with the observed epigenetic alterations. Large numbers of RDTs investigated were down-regulated along with age. Bisulfite sequencing revealed that CpG dinucleotide methylation patterns at the repeats of the RDT promoter region were mostly well maintained during aging. ChIP assay showed that histones were deacetylated at the promoter region of RDTs in aged mice brain. The observations indicate that the transcriptional repression of RDTs appears to be related to histone hypoacetylation, but not to DNA hypermethylation at repeat DNA elements in the brains of aged mice.

The 1.688 repetitive DNA of Drosophila: concerted evolution at different genomic scales and association with genes

Concerted evolution leading to homogenization of tandemly repeated DNA arrays is widespread and important for genome evolution. We investigated the range and nature of the process at chromosomal and array levels using the 1.688 tandem repeats of Drosophila melanogaster where large arrays are present in the heterochromatin of chromosomes 2, 3, and X, and short arrays are found in the euchromatin of the same chromosomes. Analysis of 326 euchromatic and heterochromatic repeats from 52 arrays showed that the homogenization of 1.688 repeats occurred differentially for distinct genomic regions, from euchromatin to heterochromatin and from local arrays to chromosomes. We further found that most euchromatic arrays are either close to, or are within introns of, genes. The short size of euchromatic arrays (one to five repeats) could be selectively constrained by their role as gene regulators, a situation similar to the so-called "tuning knobs."

A reverse transcriptase-dependent mechanism is essential for murine preimplantation development

LINE-1 (Long Interspersed Nuclear elements) and HERVs (Human Endogenous Retroviruses) are two families of retrotransposons which together account for about 28% of the human genome. Genes harbored within LINE-1 and HERV retrotransposons, particularly that encoding the reverse transcriptase (RT) enzyme, are generally expressed at low levels in differentiated cells, but their expression is up-regulated in embryonic tissues and transformed cells. Here we review evidence indicating that the LINE-1-encoded RT plays regulatory roles in early embryonic development. Indeed, antisense-mediated inhibition of expression of a highly expressed LINE-1 family in mouse zygotes caused developmental arrest at the two- or four-cell embryo stages. Development is also arrested when the embryo endogenous RT activity is pharmacologically inhibited by nevirapine, an RT inhibitor currently employed in AIDS treatment. The arrest of embryonic development is irreversible even after RT inhibition is removed and it is associated with subverted gene expression profiles. These data indicate an early requirement for LINE-1-encoded RT to support early developmental progression. Consistent with this, recent findings indicate that a reverse transcription wave is triggered in the zygote a few hours after fertilization and is propagated at least through the first two rounds of cell division. On the whole these findings suggest that reverse transcription is strictly required in early embryos as a key component of a novel RT-dependent mechanism that regulated the proper unfolding of the developmental program.

Retrotransposon-encoded reverse transcriptase in the genesis, progression and cellular plasticity of human cancer

LINE-1 (Long Interspersed Nuclear Elements) and HERVs (Human Endogenous Retroviruses) are two families of autonomously replicating retrotransposons that together account for about 28% of the human genome. Genes harbored within LINE-1 and HERV retrotransposons, particularly those encoding the reverse transcriptase (RT) enzyme, are generally expressed at low levels in differentiated cells, but their expression is upregulated in transformed cells and embryonic tissues. Here we discuss a recently discovered RT-dependent mechanism that operates in tumorigenesis and reversibly modulates phenotypic and functional variations associated with tumor progression. Downregulation of active LINE-1 elements drastically reduces the tumorigenic potential of cancer cells, paralleled by reduced proliferation and increased differentiation. Pharmacological RT inhibitors (e.g., nevirapine and efavirenz) exert similar effects on tumorigenic cell lines, both in culture and in animal models. The HERV-K family play a distinct complementary role in stress-dependent transition of melanoma cells from an adherent, non-aggressive, to a non-adherent, highly malignant, growth phenotype. In synthesis, the retrotransposon-encoded RT is increasingly emerging as a key regulator of tumor progression and a promising target in a novel anti-cancer therapy.

Different DNA-PKcs functions in the repair of radiation-induced and spontaneous DSBs within interstitial telomeric sequences

Interstitial telomeric sequences (ITSs) in hamster cells are hot spots for spontaneous and induced chromosome aberrations (CAs). Most data on ITS instability to date have been obtained in DNA repair-proficient cells. The classical non-homologous end joining repair pathway (C-NHEJ), which is the principal double strand break (DSB) repair mechanism in mammalian cells, is thought to restore the morphologically correct chromosome structure. The production of CAs thus involves DNA-PKcs-independent repair pathways. In our current study, we investigated the participation of DNA-PKcs from the C-NHEJ pathway in the repair of spontaneous or radiation-induced DSBs in ITSs using wild-type and DNA-PKcs mutant Chinese hamster ovary cells. Our data demonstrate that DNA-PKcs stabilizes spontaneous DSBs within ITSs from the chromosome 9 long arm, leading to the formation of terminal deletions. In addition, we show that DNA-PKcs-dependent C-NHEJ is employed following radiation-induced DSBs in other ITSs and restores morphologically correct chromosomes, whereas DNA-PKcs independent mechanisms co-exist in DNA-PKcs proficient cells leading to an excess of CAs within ITSs.

Genome-wide identification and annotation of HIF-1α binding sites in two cell lines using massively parallel sequencing

We identified 531 and 616 putative HIF-1α target sites by ChIP-Seq in the cancerous cell line DLD-1 and the non-cancerous cell line TIG-3, respectively. We also examined the positions and expression levels of transcriptional start sites (TSSs) in these cell lines using our TSS-Seq method. We observed that 121 and 48 genes in DLD-1 and TIG-3 cells, respectively, had HIF-1α binding sites in proximal regions of the previously reported TSSs that were up-regulated at the transcriptional level. In addition, 193 and 123 of the HIF-1α target sites, respectively, were located in proximal regions of previously uncharacterized TSSs, namely, TSSs of putative alternative promoters of protein-coding genes or promoters of putative non-protein-coding transcripts. The hypoxic response of DLD-1 cells was more significant than that of TIG-3 cells with respect to both the number of target sites and the degree of induced changes in transcript expression. The Nucleosome-Seq and ChIP-Seq analyses of histone modifications revealed that the chromatin formed an open structure in regions surrounding the HIF-1α binding sites, but this event occurred prior to the actual binding of HIF-1α. Different cellular histories may be encoded by chromatin structures and determine the activation of specific genes in response to hypoxic shock.

ELECTRONIC SUPPLEMENTARY MATERIAL

The online version of this article (doi:10.1007/s11568-011-9150-9) contains supplementary material, which is available to authorized users.

Dioxin receptor and SLUG transcription factors regulate the insulator activity of B1 SINE retrotransposons via an RNA polymerase switch

Complex genomes utilize insulators and boundary elements to help define spatial and temporal gene expression patterns. We report that a genome-wide B1 SINE (Short Interspersed Nuclear Element) retrotransposon (B1-X35S) has potent intrinsic insulator activity in cultured cells and live animals. This insulation is mediated by binding of the transcription factors dioxin receptor (AHR) and SLUG (SNAI2) to consensus elements present in the SINE. Transcription of B1-X35S is required for insulation. While basal insulator activity is maintained by RNA polymerase (Pol) III transcription, AHR-induced insulation involves release of Pol III and engagement of Pol II transcription on the same strand. B1-X35S insulation is also associated with enrichment of heterochromatin marks H3K9me3 and H3K27me3 downstream of B1-X35S, an effect that varies with cell type. B1-X35S binds parylated CTCF and, consistent with a chromatin barrier activity, its positioning between two adjacent genes correlates with their differential expression in mouse tissues. Hence, B1 SINE retrotransposons represent genome-wide insulators activated by transcription factors that respond to developmental, oncogenic, or toxicological stimuli.

Transcriptional activity analysis of promoter region of human PAX9 gene under dexamethasone, retinoic acid, and ergocalciferol treatment in MCF-7 and MDPC23

PAX9 gene is a member of the family homeobox of transcription factors and performs important function in development and organogenesis. Mutations in PAX9 coding sequences have been implicated in autosomal dominant oligodontia affecting predominantly permanent molars and second premolars. Previous studies have shown that PAX9 is required for secondary palate development and teratogens have been identified as inducers of a tooth and craniofacial malformations. This work focused on the analysis on the 5'-flanking region of the PAX9 gene studying the influence of retinoic acid, dexamethasone, and vitamin D on the expression of PAX9 by expression constructs that carry the reporter gene luciferase. As results, retinoic acid and dexamethasone showed progressive decrease of PAX9 expression. PAX9-pGL3B1 and PAX9-pGL3B2 promoter was inhibited under the treatment of dexamethasone and ergocalciferol. Retinoic acid and dexamethasone did not alter PAX9-pGL3B3 behavior indicating that sequences present between -1106 and +92 were important for the transcriptional activity of PAX9 promoter. In this study, we characterized the transcriptional activity of specific regions of the PAX9 promoter gene and we demonstrated that retinoic acid and ergocalciferol can modulate the transcriptional activity of PAX9 gene.

RISCI--Repeat Induced Sequence Changes Identifier: a comprehensive, comparative genomics-based, in silico subtractive hybridization pipeline to identify repeat induced sequence changes in closely related genomes

Background -

The availability of multiple whole genome sequences has facilitated in silico identification of fixed and polymorphic transposable elements (TE). Whereas polymorphic loci serve as makers for phylogenetic and forensic analysis, fixed species-specific transposon insertions, when compared to orthologous loci in other closely related species, may give insights into their evolutionary significance. Besides, TE insertions are not isolated events and are frequently associated with subtle sequence changes concurrent with insertion or post insertion. These include duplication of target site, 3' and 5' flank transduction, deletion of the target locus, 5' truncation or partial deletion and inversion of the transposon, and post insertion changes like inter or intra element recombination, disruption etc. Although such changes have been studied independently, no automated platform to identify differential transposon insertions and the associated array of sequence changes in genomes of the same or closely related species is available till date. To this end, we have designed RISCI - 'Repeat Induced Sequence Changes Identifier' - a comprehensive, comparative genomics-based, in silico subtractive hybridization pipeline to identify differential transposon insertions and associated sequence changes using specific alignment signatures, which may then be examined for their downstream effects.

Results -

We showcase the utility of RISCI by comparing full length and truncated L1HS and AluYa5 retrotransposons in the reference human genome with the chimpanzee genome and the alternate human assemblies (Celera and HuRef). Comparison of the reference human genome with alternate human assemblies using RISCI predicts 14 novel polymorphisms in full length L1HS, 24 in truncated L1HS and 140 novel polymorphisms in AluYa5 insertions, besides several insertion and post insertion changes. We present comparison with two previous studies to show that RISCI predictions are broadly in agreement with earlier reports. We also demonstrate its versatility by comparing various strains of Mycobacterium tuberculosis for IS 6100 insertion polymorphism.

Conclusions -

RISCI combines comparative genomics with subtractive hybridization, inferring changes only when exclusive to one of the two genomes being compared. The pipeline is generic and may be applied to most transposons and to any two or more genomes sharing high sequence similarity. Such comparisons, when performed on a larger scale, may pull out a few critical events, which may have seeded the divergence between the two species under comparison.

Cytogenetic molecular delineation of a terminal 18q deletion suggesting neo-telomere formation

Deletion of the long arm of chromosome 18 is one of the most common segmental aneusomies compatible with life and usually involves a deletion of the terminal chromosomal region. However, the mechanisms implicated in the stabilization of terminal deletions are not well understood. In this study, we analyzed a girl with moderate mental retardation who had a cytogenetically visible terminal 18q deletion. In order to characterize the breakpoint in the terminal 18q region, we used fluorescence In situ hybridization (FISH) with bacterial artificial chromosomes (BACs) and pan-telomeric probes and also the array technique based on comparative genomic hybridization (array-CGH). FISH with pan-telomeric probes revealed no signal in the terminal region of the deleted chromosome, indicating the absence of normal telomere repeat (TTAGGG)n sequences in 18q. We suggest that neo-telomere formation by chromosome healing was involved in the repair and stabilization of this terminal deletion.

Microsatellite instability in radiation-induced murine tumours; influence of tumour type and radiation quality

PURPOSE

To investigate microsatellite instability (MSI) in radiation-induced murine tumours, its dependence on tissue (haemopoietic, intestinal, mammary, brain and skin) and radiation type.

MATERIALS AND METHODS

DNA from spontaneous, X-ray or neutron-induced mouse tumours were used in Polymerase Chain Reactions (PCR) with mono- or di-nucleotide repeat markers. Deviations from expected allele size caused by insertion/deletion events were assessed by capillary electrophoresis.

RESULTS

Tumours showing MSI increased from 16% in spontaneously arising tumours to 23% (P = 0.014) in X-ray-induced tumours and rising again to 83% (P << 0.001) in neutron-induced tumours. X-ray-induced Acute Myeloid Leukaemias (AML) had a higher level of mono-nucleotide instability (45%) than di-nucleotide instability (37%). Fifty percent of neutron-induced tumours were classified as MSI-high for mono-nucleotide markers and 10% for di-nucleotide markers. Distribution of MSI varied in the different tumour types and did not appear random.

CONCLUSIONS

Exposure to ionising radiation, especially neutrons, promotes the development of MSI in mouse tumours. MSI may therefore play a role in mouse radiation tumourigenesis, particularly following high Linear Energy Transfer (LET) exposures. MSI events, for a comparable panel of genome-wide markers in different tissue types, were not randomly distributed throughout the genome.

Impact of copy number of distinct SV40PolyA segments on expression of a GFP reporter gene

The presence of Alu repeats downregulates the expression of the green fluorescent protein (GFP) gene. We found that SV40PolyA (PolyA, 240 bp), in either orientation, eliminated the inhibition of GFP gene expression induced by Alu repeats when it was placed between the GFP gene and the Alu repeats. In this study, 4 different segments (each 60 bp) were amplified from antisense PolyA (PolyAas) by PCR, and inserted upstream of Alu14 in pAlu14 plasmid (14 Alu repeats inserted downstream of the GFP gene in vector pEGFP-C1 in a head-tail tandem manner). Segments 1F1R (the first 60 bp segment at the 5' end of PolyAas) and 4F4R (the fourth 60 bp segment from the 5' end of PolyAas) did not activate GFP gene expression, whereas 2F2R and 3F3R (the middle two segments) did (as detected by Northern blot analysis and fluorescent microscopy). Different copy numbers of 2F2R and 3F3R segments, in a head and tail tandem manner, were inserted downstream of the GFP gene in pAlu14. p2F2R*4-Alu28, p3F3R*4-Alu18 and p3F3R*4-Alu28 were used as length controls to verify that the decrease in the expression of GFP was not due to the increased length of the inserted segment in the expression vectors. We found that 2 and 4 copies of 2F2R or 3F3R activated the GFP gene more strongly than one copy of them. However, more than 8 copies of 2F2R or 3F3R reduced the activation of the GFP gene. We concluded that SV40PolyAas contained at least two gene-activating elements (2F2R and 3F3R) and 2-4 copies of 2F2R or 3F3R were optimal for the expression of the GFP gene.

Repeat performance: how do genome packaging and regulation depend on simple sequence repeats?

Non-coding DNA has consistently increased during evolution of higher eukaryotes. Since the number of genes has remained relatively static during the evolution of complex organisms, it is believed that increased degree of sophisticated regulation of genes has contributed to the increased complexity. A higher proportion of non-coding DNA, including repeats, is likely to provide more complex regulatory potential. Here, we propose that repeats play a regulatory role by contributing to the packaging of the genome during cellular differentiation. Repeats, and in particular the simple sequence repeats, are proposed to serve as landmarks that can target regulatory mechanisms to a large number of genomic sites with the help of very few factors and regulate the linked loci in a coordinated manner. Repeats may, therefore, function as common target sites for regulatory mechanisms involved in the packaging and dynamic compartmentalization of the chromatin into active and inactive regions during cellular differentiation.

Mug28, a meiosis-specific protein of Schizosaccharomyces pombe, regulates spore wall formation

The meiosis-specific mug28⁺ gene of Schizosaccharomyces pombe encodes a putative RNA-binding protein. mug28Δ cells generated spores with low viability, due to the aberrant FSM formation. Meu14-GFP in mug28Δ cells showed that the FSM formed extra membranes with buds. We conclude that Mug28 is essential for the proper maturation of the FSM and the spore wall.

The meiosis-specific mug28⁺ gene of Schizosaccharomyces pombe encodes a putative RNA-binding protein with three RNA recognition motifs (RRMs). Live observations of meiotic cells that express Mug28 tagged with green fluorescent protein (GFP) revealed that Mug28 is localized in the cytoplasm, and accumulates around the nucleus from metaphase I to anaphase II. Disruption of mug28⁺ generated spores with low viability, due to the aberrant formation of the forespore membrane (FSM). Visualization of the FSM in living cells expressing GFP-tagged Psy1, an FSM protein, indicated that mug28Δ cells harbored abnormal FSMs that contained buds, and had a delayed disappearance of Meu14, a leading edge protein. Electron microscopic observation revealed that FSM formation was abnormal in mug28Δ cells, showing bifurcated spore walls that were thicker than the nonbifurcated spore walls of the wild type. Analysis of Mug28 mutants revealed that RRM3, in particular phenylalanin-466, is of primary importance for the proper localization of Mug28, spore viability, and FSM formation. Together, we conclude that Mug28 is essential for the proper maturation of the FSM and the spore wall.

Temperature and length-dependent modulation of the MH class II beta gene expression in brook charr (Salvelinus fontinalis) by a cis-acting minisatellite

It is widely recognized that the variation in gene regulation is an important factor from which evolutionary changes in diverse aspects of phenotype can be observed in all organisms. Distinctive elements with functional roles on gene regulation have been identified within the non-coding part of the genome, including repeated elements. Major histocompatibility complex (MHC) genes have been the subject of an abundant literature which made them unique candidates for studies of adaptation in natural populations. Yet, the vast majority of studies on MHC genes have dealt with patterns of polymorphism in sequence variation while very few paid attention to the possible implication of differential expression in adaptive responses. In this paper, we report the identification of a polymorphic minisatellite formed of a 32 nucleotides motif (38% G+C) involved in regulation of the major histocompatibility class II beta gene (MHII beta) of brook charr (Salvelinus fontinalis). Our main objectives were: to analyze the variability of this minisatellite found in the second intron of the MHII beta gene and to document its effect to the variation of expression level of this gene under different environmental conditions. Distinctive number of the minisatellite repeats were associated with each different MHII beta alleles identified from exon 2 sequences. Relative expression levels of specific alleles in heterozygous individuals were determined from fish lymphocytes in different genotypes. We found that alleles carrying the longest minisatellite showed a significant 1.67-2.56-fold reduction in the transcript expression relatively to the shortest one. Results obtained in three different genotypes also indicated that the repressive activity associated to the longest minisatellite was more effective at 18 degrees C compared to 6 degrees C. In contrast, no significant difference was observed in transcript levels between alleles with comparable minisatellite length at both temperatures. We also depicted a significant up-regulation of the total MHII beta transcript at 6 degrees C relative to 18 degrees C. These results reveal for the first time that a temperature-sensitive minisatellite could potentially play an important role in the gene regulation of the adaptive immune response in fishes.

Development of a novel output value for quantitative assessment in methylated DNA immunoprecipitation-CpG island microarray analysis

In DNA methylation microarray analysis, quantitative assessment of intermediate methylation levels in samples with various global methylation levels is still difficult. Here, specifically for methylated DNA immunoprecipitation-CpG island (CGI) microarray analysis, we developed a new output value. The signal log ratio reflected the global methylation levels, but had only moderate linear correlation (r = 0.72) with the fraction of DNA molecules immunoprecipitated. By multiplying the signal log ratio using a coefficient obtained from the probability value that took account of signals in neighbouring probes, its linearity was markedly improved (r = 0.94). The new output value, Me value, reflected the global methylation level, had a strong correlation also with the fraction of methylated CpG sites obtained by bisulphite sequencing (r = 0.88), and had an accuracy of 71.8 and 83.8% in detecting completely methylated and unmethylated CGIs. Analysis of gastric cancer cell lines using the Me value showed that methylation of CGIs in promoters and gene bodies was associated with low and high, respectively, gene expression. The degree of demethylation of promoter CGIs after 5-aza-2'-deoxycytidine treatment had no association with that of induction of gene expression. The Me value was considered to be useful for analysis of intermediate methylation levels of CGIs.

A novel protein isoform of the multicopy human NAIP gene derives from intragenic Alu SINE promoters

The human neuronal apoptosis inhibitory protein (NAIP) gene is no longer principally considered a member of the Inhibitor of Apoptosis Protein (IAP) family, as its domain structure and functions in innate immunity also warrant inclusion in the Nod-Like Receptor (NLR) superfamily. NAIP is located in a region of copy number variation, with one full length and four partly deleted copies in the reference human genome. We demonstrate that several of the NAIP paralogues are expressed, and that novel transcripts arise from both internal and upstream transcription start sites. Remarkably, two internal start sites initiate within Alu short interspersed element (SINE) retrotransposons, and a third novel transcription start site exists within the final intron of the GUSBP1 gene, upstream of only two NAIP copies. One Alu functions alone as a promoter in transient assays, while the other likely combines with upstream L1 sequences to form a composite promoter. The novel transcripts encode shortened open reading frames and we show that corresponding proteins are translated in a number of cell lines and primary tissues, in some cases above the level of full length NAIP. Interestingly, some NAIP isoforms lack their caspase-sequestering motifs, suggesting that they have novel functions. Moreover, given that human and mouse NAIP have previously been shown to employ endogenous retroviral long terminal repeats as promoters, exaptation of Alu repeats as additional promoters provides a fascinating illustration of regulatory innovations adopted by a single gene.

Bead-probe complex capture a couple of SINE and LINE family from genomes of two closely related species of East Asian cyprinid directly using magnetic separation

Background

Short and long interspersed elements (SINEs and LINEs, respectively), two types of retroposons, are active in shaping the architecture of genomes and powerful tools for studies of phylogeny and population biology. Here we developed special protocol to apply biotin-streptavidin bead system into isolation of interspersed repeated sequences rapidly and efficiently, in which SINEs and LINEs were captured directly from digested genomic DNA by hybridization to bead-probe complex in solution instead of traditional strategy including genomic library construction and screening.

Results

A new couple of SINEs and LINEs that shared an almost identical 3'tail was isolated and characterized in silver carp and bighead carp of two closely related species. These SINEs (34 members), designated HAmo SINE family, were little divergent in sequence and flanked by obvious TSD indicated that HAmo SINE was very young family. The copy numbers of this family was estimated to 2 × 10⁵ and 1.7 × 10⁵ per haploid genome by Real-Time qPCR, respectively. The LINEs, identified as the homologs of LINE2 in other fishes, had a conserved primary sequence and secondary structures of the 3'tail region that was almost identical to that of HAmo SINE. These evidences suggest that HAmo SINEs are active and amplified recently utilizing the enzymatic machinery for retroposition of HAmoL2 through the recognition of higher-order structures of the conserved 42-tail region. We analyzed the possible structures of HAmo SINE that lead to successful amplification in genome and then deduced that HAmo SINE, SmaI SINE and FokI SINE that were similar in sequence each other, were probably generated independently and created by LINE family within the same lineage of a LINE phylogeny in the genomes of different hosts.

Conclusion

The presented results show the advantage of the novel method for retroposons isolation and a pair of young SINE family and its partner LINE family in two carp fishes, which strengthened the hypotheses containing the slippage model for initiation of reverse transcription, retropositional parasitism of SINEs on LINEs, the formation of the stem loop structure in 3'tail region of some SINEs and LINEs and the mechanism of template switching in generating new SINE family.

Characterization and evolutionary landscape of AmnSINE1 in Amniota genomes

Discovery of a large number of conserved non-coding elements (CNEs) in vertebrate genomes provides a cornerstone to elucidate molecular mechanisms of macroevolution. Extensive comparative genomics has proven that transposons such as short interspersed elements (SINEs) were an important source of CNEs. We recently characterized AmnSINE1, a SINE family in Amniota genomes, some of which are present in CNEs, and demonstrated that two AmnSINE1 loci play an important role in mammalian-specific brain development by functioning as an enhancer (Sasaki et al. Proc. Natl. Acad. Sci. USA 2008). To get more information about AmnSINE1s, we here performed a multi-species search for AmnSINE1, and revealed the distribution and evolutionary history of these SINEs in amniote genomes. The number of AmnSINE1 regions in amniotes ranged from 160 to 1200; the number in the eutherians were under 500 and the largest was that in chicken. Phylogenetic analysis established that each AmnSINE1 locus has evolved uniquely, primarily since the divergence of mammals from reptiles. These results support the notion that AmnSINE1s were amplified as an ancient retroposon in a common ancestor of Amniota and subsequently have survived for 300 Myr because of functions acquired by mutation-coupled exaptation prior mammalian radiation. On the basis of sequence homology and conserved synteny, we detected the orthologs of AmnSINE1 for candidates of further enhancer analysis, which are more conserved than two loci that were shown to have been involved in mammalian brain development. The present work provides a comprehensive data set to test the role of AmnSINE1s, many of which were exapted and contributed to mammalian macroevolution.