Expressing genes do not forget their LINEs: transposable elements and gene expression

Impact of prematurity on LINE-1 promoter methylation

Aim: Promoter methylation of LINE-1 may be affected by prematurity, but there is little evidence in the literature.Materials & methods: Blood from premature and full-term neonates on days 0, 5, 30 and 90 was analyzed for DNA methylation percentage in a promoter region of the LINE-1, after bisulfite conversion and pyrosequencing.Results: Premature infants, as a whole, showed significantly lower methylation percentage at birth, but this difference diminished over time. However, the subgroup of extremely premature (<28 weeks gestational age) had higher methylation percentages, similar to full-term newborns.Conclusion: This research underscores the critical role of prematurity on the methylation pattern of LINE-1. These findings underline the complexity of epigenetic regulation in prematurity and emphasize the need for further studies.

The evolution of antimicrobial peptides in Chiroptera

High viral tolerance coupled with an extraordinary regulation of the immune response makes bats a great model to study host-pathogen evolution. Although many immune-related gene gains and losses have been previously reported in bats, important gene families such as antimicrobial peptides (AMPs) remain understudied. We built an exhaustive bioinformatic pipeline targeting the major gene families of defensins and cathelicidins to explore AMP diversity and analyze their evolution and distribution across six bat families. A combination of manual and automated procedures identified 29 AMP families across queried species, with α-, β-defensins, and cathelicidins representing around 10% of AMP diversity. Gene duplications were inferred in both α-defensins, which were absent in five species, and three β-defensin gene subfamilies, but cathelicidins did not show significant shifts in gene family size and were absent in Anoura caudifer and the pteropodids. Based on lineage-specific gains and losses, we propose diet and diet-related microbiome evolution may determine the evolution of α- and β-defensins gene families and subfamilies. These results highlight the importance of building species-specific libraries for genome annotation in non-model organisms and shed light on possible drivers responsible for the rapid evolution of AMPs. By focusing on these understudied defenses, we provide a robust framework for explaining bat responses to pathogens.

Transposable elements and their role in aging

Transposable elements (TEs) are an important part of eukaryotic genomes. The role of somatic transposition in aging, carcinogenesis, and other age-related diseases has been determined. This review discusses the fundamental properties of TEs and their complex interactions with cellular processes, which are crucial for understanding the diverse effects of their activity on the genetics and epigenetics of the organism. The interactions of TEs with recombination, replication, repair, and chromosomal regulation; the ability of TEs to maintain a balance between their own activity and repression, the involvement of TEs in the creation of new or alternative genes, the expression of coding/non-coding RNA, and the role in DNA damage and modification of regulatory networks are reviewed. The contribution of the derepressed TEs to age-dependent effects in individual cells/tissues in different organisms was assessed. Conflicting information about TE activity under stress as well as theories of aging mechanisms related to TEs is discussed. On the one hand, transposition activity in response to stressors can lead to organisms acquiring adaptive innovations of great importance for evolution at the population level. On the other hand, the TE expression can cause decreased longevity and stress tolerance at the individual level. The specific features of TE effects on aging processes in germline and soma and the ways of their regulation in cells are highlighted. Recent results considering somatic mutations in normal human and animal tissues are indicated, with the emphasis on their possible functional consequences. In the context of aging, the correlation between somatic TE activation and age-related changes in the number of proteins required for heterochromatin maintenance and longevity regulation was analyzed. One of the original features of this review is a discussion of not only effects based on the TEs insertions and the associated consequences for the germline cell dynamics and somatic genome, but also the differences between transposon- and retrotransposon-mediated structural genome changes and possible phenotypic characteristics associated with aging and various age-related pathologies. Based on the analysis of published data, a hypothesis about the influence of the species-specific features of number, composition, and distribution of TEs on aging dynamics of different animal genomes was formulated.

Human L1 Transposition Dynamics Unraveled with Functional Data Analysis

Long INterspersed Elements-1 (L1s) constitute >17% of the human genome and still actively transpose in it. Characterizing L1 transposition across the genome is critical for understanding genome evolution and somatic mutations. However, to date, L1 insertion and fixation patterns have not been studied comprehensively. To fill this gap, we investigated three genome-wide data sets of L1s that integrated at different evolutionary times: 17,037 de novo L1s (from an L1 insertion cell-line experiment conducted in-house), and 1,212 polymorphic and 1,205 human-specific L1s (from public databases). We characterized 49 genomic features-proxying chromatin accessibility, transcriptional activity, replication, recombination, etc.-in the ±50 kb flanks of these elements. These features were contrasted between the three L1 data sets and L1-free regions using state-of-the-art Functional Data Analysis statistical methods, which treat high-resolution data as mathematical functions. Our results indicate that de novo, polymorphic, and human-specific L1s are surrounded by different genomic features acting at specific locations and scales. This led to an integrative model of L1 transposition, according to which L1s preferentially integrate into open-chromatin regions enriched in non-B DNA motifs, whereas they are fixed in regions largely free of purifying selection-depleted of genes and noncoding most conserved elements. Intriguingly, our results suggest that L1 insertions modify local genomic landscape by extending CpG methylation and increasing mononucleotide microsatellite density. Altogether, our findings substantially facilitate understanding of L1 integration and fixation preferences, pave the way for uncovering their role in aging and cancer, and inform their use as mutagenesis tools in genetic studies.

DNA methylation patterns of LINE-1 and Alu for pre-symptomatic dementia in type 2 diabetes

The identification of early markers of dementia is important for higher-risk populations such as those with type 2 diabetes (T2D). Retrotransposons, including long interspersed nuclear element 1 (LINE-1) and Alu, comprise ~40% of the human genome. Although dysregulation of these retrotransposons can induce aberrant gene regulation and genomic instability, their role in the development of pre-symptomatic dementia (PSD) among T2D patients is unknown. Here, we examined locus-specific changes in LINE-1 and Alu methylation in PSD and the potential to offset these changes via supplementation with folate and vitamin B₁₂. We interrogated DNA methylation patterns corresponding to 22,352 probes for LINE-1 and Alu elements using publicly-available Illumina Infinium 450K methylation datasets from i) an 18-month prospective study in 28 T2D patients (GSE62003) and ii) an intervention study in which 44 individuals were supplemented with folic acid (400 μg/day) and vitamin B₁₂ (500 μg/day) over two years (GSE74548). We identified 714 differentially methylated positions (DMP) mapping to retrotransposons in T2D patients who developed PSD in comparison to those who did not (P_FDR < 0.05), comprised of 2.4% (228 probes) of all LINE-1 probes and 3.8% (486 probes) of all Alu probes. These loci were enriched in genes with functions related to Alzheimer's disease and cognitive decline, including GNB5, GNG7 and PKN3 (p < 0.05). In older individuals supplemented with folate/vitamin B₁₂, 85 (11.9%) PSD retrotransposon loci showed significant changes in methylation (p < 0.05): participants with the MTHFR CC genotype predominantly showed hypermethylation at these loci, while hypomethylation was observed more frequently in those with the TT genotype. In T2D patients, LINE-1 and Alu elements are differentially methylated in PSD in a locus-specific manner and may offer clinical utility in monitoring risk of dementia. Further work is required to examine the potential for dietary supplementation in lowering the risk of PSD.

Agouti Signaling Protein and Its Receptors as Potential Molecular Markers for Intramuscular and Body Fat Deposition in Cattle

Transcriptome analyses of bovine muscle tissue differing in intramuscular fat (IMF) content identified agouti signaling protein (ASIP) as a promising candidate gene for fat deposition. The protein is secreted from adipocytes and may serve as a signaling molecule in cross-talk between adipocytes and muscle fibers or other cells. Known receptors for ASIP are the melanocortin receptors (e.g., MC4R) and attractin (ATRN). The present study was conducted to determine relationships between the expression of ASIP and its receptors in different bovine tissues with fat deposition. Adipose tissues, liver, and longissimus muscle tissue were collected from 246 F₂-generation bulls (Charolais × Holstein cross) and gene expression was measured with RT-qPCR. During analysis of subcutaneous fat (SCF) of all bulls, 17 animals were identified with a transposon-derived transcript (Exon2C) inserted in the ASIP gene and dramatically increased ASIP mRNA levels. Significant correlations between normalized mRNA values of SCF and phenotypic traits related to fat deposition were found in bulls without Exon2C. Three retrospectively assigned groups [Exon2C, n = 17; high carcass fat (HCF), n = 20; low carcass fat (LCF), n = 20] were further analyzed to verify expression differences and elucidate molecular reasons. Expression of ASIP could be detected in isolated muscle fibers and adipocytes of Exon2C bulls in contrast to HCF and LCF bulls, indicating ectopic ASIP expression if the transposon is present. Among adipose tissues, highest ASIP mRNA levels were measured in SCF with significantly higher values in HCF compared to LCF bulls (1.6-fold, P < 0.05). However, the protein abundance was below the detection limit in all bulls. Potential ASIP receptors were detected in most investigated tissues. The expression of MC4R was higher and of ATRN was lower in several tissues of LCF compared to HCF bulls, whereas MC1R was not differentially expressed. Bulls of the Exon2C group had lower ATRN mRNA values than HCF and LCF bulls in perirenal fat (PF), but higher (P < 0.05) values in muscle. Receptors were also expressed in tissues where ASIP mRNA was not detected. Consequently, those tissues could be targets for ASIP if it circulates.

The epigenetic alterations of endogenous retroelements in aging

Endogenous retroelements, transposons that mobilize through RNA intermediates, include some of the most abundant repetitive sequences of the human genome, such as Alu and LINE-1 sequences, and human endogenous retroviruses. Recent discoveries demonstrate that these mobile genetic elements not only act as intragenomic parasites, but also exert regulatory roles in living cells. The risk of genomic instability represented by endogenous retroelements is normally counteracted by a series of epigenetic control mechanisms which include, among the most important, CpG DNA methylation. Indeed, most of the genomic CpG sites subjected to DNA methylation in the nuclear DNA are carried by these repetitive elements. As other parts of the genome, endogenous retroelements and other transposable elements are subjected to deep epigenetic alterations during aging, repeatedly observed in the context of organismal and cellular senescence, in human and other species. This review summarizes the current status of knowledge about the epigenetic alterations occurring in this large, non-genic portion of the genome in aging and age-related conditions, with a focus on the causes and the possible functional consequences of these alterations.

Bioinformatic analysis of regulation of natural antisense transcripts by transposable elements in human mRNA

Non-coding RNA is no longer considered to be "junk" DNA, based on evidence uncovered in recent decades. In particular, the important role played by natural antisense transcripts (NATs) in regulating the expression of genes is receiving increasing attention. However, the regulatory mechanisms of NATs remain incompletely understood. It is well-known that the insertion of transposable elements (TEs) can affect gene transcription. Using a bioinformatics approach, we identified NATs using human mRNA sequences from the UCSC Genome Browser Database. Our in silico analysis identified 1079 NATs and 700 sense-antisense gene pairs. We identified 179 NATs that showed evidence of having been affected by TEs during cellular gene expression. These findings may provide an understanding of the complex regulation mechanisms of NATs. If our understanding of NATs as modulators of gene expression is further enhanced, we can develop ways to control gene expression.

An intronic LINE-1 insertion in MERTK is strongly associated with retinopathy in Swedish Vallhund dogs

The domestic dog segregates a significant number of inherited progressive retinal diseases, several of which mirror human retinal diseases and which are collectively termed progressive retinal atrophy (PRA). In 2014, a novel form of PRA was reported in the Swedish Vallhund breed, and the disease was mapped to canine chromosome 17. The causal mutation was not identified, but expression analyses of the retinas of affected Vallhunds demonstrated a 6-fold increased expression of the MERTK gene compared to unaffected dogs. Using 24 retinopathy cases and 97 controls with no clinical signs of retinopathy, we replicated the chromosome 17 association in Swedish Vallhunds from the UK and aimed to elucidate the causal variant underlying this association using whole genome sequencing (WGS) of an affected dog. This revealed a 6-8 kb insertion in intron 1 of MERTK that was not present in WGS of 49 dogs of other breeds. Sequencing and BLASTN analysis of the inserted segment was consistent with the insertion comprising a full-length intact LINE-1 retroelement. Testing of the LINE-1 insertion for association with retinopathy in the UK set of 24 cases and 97 controls revealed a strong statistical association (P-value 6.0 x 10-11) that was subsequently replicated in the original Finnish study set (49 cases and 89 controls (P-value 4.3 x 10-19). In a pooled analysis of both studies (73 cases and 186 controls), the LINE-1 insertion was associated with a ~20-fold increased risk of retinopathy (odds ratio 23.41, 95% confidence intervals 10.99-49.86, P-value 1.3 x 10-27). Our study adds further support for regulatory disruption of MERTK in Swedish Vallhund retinopathy; however, further work is required to establish a functional overexpression model. Future work to characterise the mechanism by which this intronic mutation disrupts gene regulation will further improve the understanding of MERTK biology and its role in retinal function.

Decoupling of DNA methylation and activity of intergenic LINE-1 promoters in colorectal cancer

Hypomethylation of LINE-1 repeats in cancer has been proposed as the main mechanism behind their activation; this assumption, however, was based on findings from early studies that were biased toward young and transpositionally active elements. Here, we investigate the relationship between methylation of 2 intergenic, transpositionally inactive LINE-1 elements and expression of the LINE-1 chimeric transcript (LCT) 13 and LCT14 driven by their antisense promoters (L1-ASP). Our data from DNA modification, expression, and 5'RACE analyses suggest that colorectal cancer methylation in the regions analyzed is not always associated with LCT repression. Consistent with this, in HCT116 colorectal cancer cells lacking DNA methyltransferases DNMT1 or DNMT3B, LCT13 expression decreases, while cells lacking both DNMTs or treated with the DNMT inhibitor 5-azacytidine (5-aza) show no change in LCT13 expression. Interestingly, levels of the H4K20me3 histone modification are inversely associated with LCT13 and LCT14 expression. Moreover, at these LINE-1s, H4K20me3 levels rather than DNA methylation seem to be good predictor of their sensitivity to 5-aza treatment. Therefore, by studying individual LINE-1 promoters we have shown that in some cases these promoters can be active without losing methylation; in addition, we provide evidence that other factors (e.g., H4K20me3 levels) play prominent roles in their regulation.

Transposable elements in the Anopheles funestus transcriptome

Transposable elements (TEs) are present in most of the eukaryotic genomes and their impact on genome evolution is increasingly recognized. Although there is extensive information on the TEs present in several eukaryotic genomes, less is known about the expression of these elements at the transcriptome level. Here we present a detailed analysis regarding the expression of TEs in Anopheles funestus, the second most important vector of human malaria in Africa. Several transcriptionally active TE families belonging both to Class I and II were identified and characterized. Interestingly, we have identified a full-length putative active element (including the presence of full length TIRs in the genomic sequence) belonging to the hAT superfamily, which presents active members in other insect genomes. This work contributes to a comprehensive understanding of the landscape of transposable elements in A. funestus transcriptome. Our results reveal that TEs are abundant and diverse in the mosquito and that most of the TE families found in the genome are represented in the mosquito transcriptome, a fact that could indicate activity of these elements.The vast diversity of TEs expressed in A. funestus suggests that there is ongoing amplification of several families in this organism.

The transposable element environment of human genes is associated with histone and expression changes in cancer

BACKGROUND

Only 2 % of the human genome code for proteins. Among the remaining 98 %, transposable elements (TEs) represent millions of sequences. TEs have an impact on genome evolution by promoting mutations. Especially, TEs possess their own regulatory sequences and can alter the expression pattern of neighboring genes. Since they can potentially be harmful, TE activity is regulated by epigenetic mechanisms. These mechanisms participate in the modulation of gene expression and can be associated with some human diseases resulting from gene expression deregulation. The fact that the TE silencing can be removed in cancer could explain a part of the changes in gene expression. Indeed, epigenetic modifications associated locally with TE sequences could impact neighboring genes since these modifications can spread to adjacent sequences.

RESULTS

We compared the histone enrichment, TE neighborhood, and expression divergence of human genes between a normal and a cancer conditions. We show that the presence of TEs near genes is associated with greater changes in histone enrichment and that differentially expressed genes harbor larger histone enrichment variation related to the presence of particular TEs.

CONCLUSIONS

Taken together, these results suggest that the presence of TEs near genes could favor important variation in gene expression when the cell environment is modified.

LINE1 insertions as a genomic risk factor for schizophrenia: Preliminary evidence from an affected family

Recent studies show that human-specific LINE1s (L1HS) play a key role in the development of the central nervous system (CNS) and its disorders, and that their transpositions within the human genome are more common than previously thought. Many polymorphic L1HS, that is, present or absent across individuals, are not annotated in the current release of the genome and are customarily termed "non-reference L1s." We developed an analytical workflow to identify L1 polymorphic insertions with next-generation sequencing (NGS) using data from a family in which SZ segregates. Our workflow exploits two independent algorithms to detect non-reference L1 insertions, performs local de novo alignment of the regions harboring predicted L1 insertions and resolves the L1 subfamily designation from the de novo assembled sequence. We found 110 non-reference L1 polymorphic loci exhibiting Mendelian inheritance, the vast majority of which are already reported in dbRIP and/or euL1db, thus, confirming their status as non-reference L1 polymorphic insertions. Four previously undetected L1 polymorphic loci were confirmed by PCR amplification and direct sequencing of the insert. A large fraction of our non-reference L1s is located within the open reading frame of protein-coding genes that belong to pathways already implicated in the pathogenesis of schizophrenia. The finding of these polymorphic variants among SZ offsprings is intriguing and suggestive of putative pathogenic role. Our data show the utility of NGS to uncover L1 polymorphic insertions, a neglected type of genetic variants with the potential to influence the risk to develop schizophrenia like SNVs and CNVs. © 2016 Wiley Periodicals, Inc.

The endonuclease domain of the LINE-1 ORF2 protein can tolerate multiple mutations

Background

Approximately 17 % of the human genome is comprised of the Long INterspersed Element-1 (LINE-1 or L1) retrotransposon, the only currently active autonomous family of retroelements. Though L1 elements have helped to shape mammalian genome evolution over millions of years, L1 activity can also be mutagenic and result in human disease. L1 expression has the potential to contribute to genomic instability via retrotransposition and DNA double-strand breaks (DSBs). Additionally, L1 is responsible for structural genomic variations induced by other transposable elements such as Alu and SVA, which rely on the L1 ORF2 protein for their propagation. Most of the genomic damage associated with L1 activity originates with the endonuclease domain of the ORF2 protein, which nicks the DNA in preparation for target-primed reverse transcription.

Results

Bioinformatic analysis of full-length L1 loci residing in the human genome identified numerous mutations in the amino acid sequence of the ORF2 endonuclease domain. Some of these mutations were found in residues which were predicted to be phosphorylation sites for cellular kinases. We mutated several of these putative phosphorylation sites in the ORF2 endonuclease domain and investigated the effect of these mutations on the function of the full-length ORF2 protein and the endonuclease domain (ENp) alone. Most of the single and multiple point mutations that were tested did not significantly impact expression of the full-length ORF2p, or alter its ability to drive Alu retrotransposition. Similarly, most of those same mutations did not significantly alter expression of ENp, or impair its ability to induce DNA damage and cause toxicity.

Conclusions

Overall, our data demonstrate that the full-length ORF2p or the ENp alone can tolerate several specific single and multiple point mutations in the endonuclease domain without significant impairment of their ability to support Alu mobilization or induce DNA damage, respectively.

Electronic supplementary material

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

Nomadic genetic elements contribute to oncogenic translocations: Implications in carcinogenesis

Chromosomal translocations as molecular signatures have been reported in various malignancies but, the mechanism behind which is largely unknown. Swapping of chromosomal fragments occurs by induction of double strand breaks (DSBs), most of which were initially assumed de novo. However, decoding of human genome proved that transposable elements (TE) might have profound influence on genome integrity. TEs are highly conserved mobile genetic elements that generate DSBs, subsequently resulting in large chromosomal rearrangements. Previously TE insertions were thought to be harmless, but recently gains attention due to the origin of spectrum of post-insertional genomic alterations and subsequent transcriptional alterations leading to development of deleterious effects mainly carcinogenesis. Though the existing knowledge on the cancer-associated TE dynamics is very primitive, exploration of underlying mechanism promises better therapeutic strategies for cancer. Thus, this review focuses on the prevalence of TE in the genome, associated genomic instability upon transposition activation and impact on tumorigenesis.

LINE-1 activity as molecular basis for genomic instability associated with light exposure at night

The original hypothesis that exposure to light at night increases risk of breast cancer via suppression of nocturnal melatonin production was proposed over 2 decades ago. In 2007, shift work that involves circadian disruption has been recognized by the World Health Organization as a probable human carcinogen. Our discovery of melatonin-dependent regulation of LINE-1 retrotransposon expression and mobilization is the latest addition to the list of cellular genes and processes that are affected by light exposure at night. This finding establishes an unexpected health relevant connection between this endogenous DNA damaging agent and environmental light exposure. It also offers an appealing hypothesis pertaining to the origin of genomic instability in the genomes of individuals with light at night- or age-associated disruption of melatonin signaling.

Association of new deletion/duplication region at chromosome 1p21 with intellectual disability, severe speech deficit and autism spectrum disorder-like behavior: an all-in approach to solving the DPYD enigma

We describe an as yet unreported neocentric small supernumerary marker chromosome (sSMC) derived from chromosome 1p21.3p21.2. It was present in 80% of the lymphocytes in a male patient with intellectual disability, severe speech deficit, mild dysmorphic features, and hyperactivity with elements of autism spectrum disorder (ASD).

Several important neurodevelopmental genes are affected by the 3.56 Mb copy number gain of 1p21.3p21.2, which may be considered reciprocal in gene content to the recently recognized 1p21.3 microdeletion syndrome. Both 1p21.3 deletions and the presented duplication display overlapping symptoms, fitting the same disorder category. Contribution of coding and non-coding genes to the phenotype is discussed in the light of cellular and intercellular homeostasis disequilibrium. In line with this the presented 1p21.3p21.2 copy number gain correlated to 1p21.3 microdeletion syndrome verifies the hypothesis of a cumulative effect of the number of deregulated genes - homeostasis disequilibrium leading to overlapping phenotypes between microdeletion and microduplication syndromes.

Although miR-137 appears to be the major player in the 1p21.3p21.2 region, deregulation of the DPYD (dihydropyrimidine dehydrogenase) gene may potentially affect neighboring genes underlying the overlapping symptoms present in both the copy number loss and copy number gain of 1p21. Namely, the all-in approach revealed that DPYD is a complex gene whose expression is epigenetically regulated by long non-coding RNAs (lncRNAs) within the locus. Furthermore, the long interspersed nuclear element-1 (LINE-1) L1MC1 transposon inserted in DPYD intronic transcript 1 (DPYD-IT1) lncRNA with its parasites, TcMAR-Tigger5b and pair of Alu repeats appears to be the “weakest link” within the DPYD gene liable to break. Identification of the precise mechanism through which DPYD is epigenetically regulated, and underlying reasons why exactly the break (FRA1E) happens, will consequently pave the way toward preventing severe toxicity to the antineoplastic drug 5-fluorouracil (5-FU) and development of the causative therapy for the dihydropyrimidine dehydrogenase deficiency.

Cell-host, LINE and environment: Three players in search of a balance

Long interspersed nuclear elements -1 (LINEs, L1s) are retroelements occupying almost 17% of the human genome. L1 retrotransposition can cause deleterious effects on the host-cell and it is generally inhibited by suppressive mechanisms, but it can occur in some specific cells during early development as well as in some tumor cells and in the presence of several environmental factors. In a recent publication we reported that extremely low frequency pulsed magnetic field can affect L1 retrotransposition in neuroblastoma cells. In this commentary we discuss the interaction between environment and L1 activity in the light of the new emerging paradigm of host-LINE relationship.

Transposable elements and psychiatric disorders

Transposable Elements (TEs) or transposons are low-complexity elements (e.g., LINEs, SINEs, SVAs, and HERVs) that make up to two-thirds of the human genome. There is mounting evidence that TEs play an essential role in genomic architecture and regulation related to both normal function and disease states. Recently, the identification of active TEs in several different human brain regions suggests that TEs play a role in normal brain development and adult physiology and quite possibly in psychiatric disorders. TEs have been implicated in hemophilia, neurofibromatosis, and cancer. With the advent of next-generation whole-genome sequencing approaches, our understanding of the relationship between TEs and psychiatric disorders will greatly improve. We will review the biology of TEs and early evidence for TE involvement in psychiatric disorders.

The vertebrate RCAN gene family: novel insights into evolution, structure and regulation

Recently there has been much interest in the Regulators of Calcineurin (RCAN) proteins which are important endogenous modulators of the calcineurin-NFATc signalling pathway. They have been shown to have a crucial role in cellular programmes such as the immune response, muscle fibre remodelling and memory, but also in pathological processes such as cardiac hypertrophy and neurodegenerative diseases. In vertebrates, the RCAN family form a functional subfamily of three members RCAN1, RCAN2 and RCAN3 whereas only one RCAN is present in the rest of Eukarya. In addition, RCAN genes have been shown to collocate with RUNX and CLIC genes in ACD clusters (ACD21, ACD6 and ACD1). How the RCAN genes and their clustering in ACDs evolved is still unknown. After analysing RCAN gene family evolution using bioinformatic tools, we propose that the three RCAN vertebrate genes within the ACD clusters, which evolved from single copy genes present in invertebrates and lower eukaryotes, are the result of two rounds of whole genome duplication, followed by a segmental duplication. This evolutionary scenario involves the loss or gain of some RCAN genes during evolution. In addition, we have analysed RCAN gene structure and identified the existence of several characteristic features that can be involved in RCAN evolution and gene expression regulation. These included: several transposable elements, CpG islands in the 5′ region of the genes, the existence of antisense transcripts (NAT) associated with the three human genes, and considerable evidence for bidirectional promoters that regulate RCAN gene expression. Furthermore, we show that the CpG island associated with the RCAN3 gene promoter is unmethylated and transcriptionally active. All these results provide timely new insights into the molecular mechanisms underlying RCAN function and a more in depth knowledge of this gene family whose members are obvious candidates for the development of future therapies.

Characterization of OsMIK in a rice mutant with reduced phytate content reveals an insertion of a rearranged retrotransposon

The rice low phytic acid (lpa) mutant Os-lpa-XS110-1(XS-lpa) has ~45 % reduction in seed phytic acid (PA) compared with the wild-type cultivar Xiushui 110. Previously, a single recessive gene mutation was shown to be responsible for the lpa phenotype and was mapped to a region of chromosome 3 near OsMIK (LOC_Os03g52760) and OsIPK1 (LOC_Os03g51610), two genes involved in PA biosynthesis. Here, we report the identification of a large insert in the intron of OsMIK in the XS-lpa mutant. Sequencing of fragments amplified through TAIL-PCRs revealed that the insert was a derivative of the LINE retrotransposon gene LOC_Os03g56910. Further analyses revealed the following characteristics of the insert and its impacts: (1) the inserted sequence of LOC_Os03g56910 was split at its third exon and rejoined inversely, with its 5' and 3' flanking sequences inward and the split third exon segments outward; (2) the LOC_Os03g56910 remained in its original locus in XS-lpa, and the insertion probably resulted from homologous recombination repair of a DNA double strand break; (3) while the OsMIK transcripts of XS-lpa and Xiushui 110 were identical, substantial reductions of the transcript abundance (~87 %) and the protein level (~60 %) were observed in XS-lpa, probably due to increased methylation in its promoter region. The above findings are discussed in the context of plant mutagenesis, epigenetics and lpa breeding.

Linear decay of retrotransposon antisense bias across genes is contingent upon tissue specificity

Retrotransposons comprise approximately half of the human genome and contribute to chromatin structure, regulatory motifs, and protein-coding sequences. Since retrotransposon insertions can disrupt functional genetic elements as well as introduce new sequence motifs to a region, they have the potential to affect the function of genes that harbour insertions as well as those nearby. Partly as a result of these effects, the distribution of retrotransposons across the genome is non-uniform and there are observed imbalances in the orientation of insertions with respect to the transcriptional direction of the containing gene. Although some of the factors underlying the observed distributions are understood, much of the variability remains unexplained. Detailed characterization of retrotransposon density in genes could help inform predictions of the functional consequence of de novo as well as polymorphic insertions. In order to characterize the relationship between genes and inserted elements, we have examined the distribution of retrotransposons and their internal motifs within tissue-specific and housekeeping genes. We have identified that the previously established retrotransposon antisense bias decays at a linear rate across genes, resulting in an equal density of sense and antisense retrotransposons near the 3'-UTR. In addition, the decay of antisense bias across genes is less pronounced among tissue-specific genes. Our results provide support for the scenario in which this linear decay in antisense bias is established by natural selection shortly after retrotransposon integration, and that total antisense bias observed is above and beyond any bias introduced by the integration process itself. Finally, we provide an example of a retrotransposon acting as an eQTL on a coincident gene, highlighting one of several possible avenues through which insertions may modulate gene function.

Transcriptional activity of transposable elements in coelacanth

The morphological stasis of coelacanths has long suggested a slow evolutionary rate. General genomic stasis might also imply a decrease of transposable elements activity. To evaluate the potential activity of transposable elements (TEs) in "living fossil" species, transcriptomic data of Latimeria chalumnae and its Indonesian congener Latimeria menadoensis were compared through the RNA-sequencing mapping procedures in three different organs (liver, testis, and muscle). The analysis of coelacanth transcriptomes highlights a significant percentage of transcribed TEs in both species. Major contributors are LINE retrotransposons, especially from the CR1 family. Furthermore, some particular elements such as a LF-SINE and a LINE2 sequences seem to be more expressed than other elements. The amount of TEs expressed in testis suggests possible transposition burst in incoming generations. Moreover, significant amount of TEs in liver and muscle transcriptomes were also observed. Analyses of elements displaying marked organ-specific expression gave us the opportunity to highlight exaptation cases, that is, the recruitment of TEs as new cellular genes, but also to identify a new Latimeria-specific family of Short Interspersed Nuclear Elements called CoeG-SINEs. Overall, transcriptome results do not seem to be in line with a slow-evolving genome with poor TE activity.

Genome-wide identification of genes regulated in trans by transposable element small interfering RNAs

Transposable elements (TEs) are known to influence the regulation of neighboring genes through a variety of mechanisms. Additionally, it was recently discovered that TEs can regulate non-neighboring genes through the trans-acting nature of small interfering RNAs (siRNAs). When the epigenetic repression of TEs is lost, TEs become transcriptionally active, and the host cell acts to repress mutagenic transposition by degrading TE mRNAs into siRNAs. In this study, we have performed a genome-wide analysis in the model plant Arabidopsis thaliana and found that TE siRNA-based regulation of genic mRNAs is more pervasive than the two formerly characterized proof-of-principle examples. We identified 27 candidate genic mRNAs that do not contain a TE fragment but are regulated through partial complementarity by the accumulation of TE siRNAs and are therefore influenced by TE epigenetic activation. We have experimentally confirmed several gene targets and demonstrated that they respond to the accumulation of specific 21 nucleotide TE siRNAs that are incorporated into the Arabidopsis Argonaute1 protein. Additionally, we found that one TE siRNA specifically targets and inhibits the formation of a host protein that acts to repress TE activity, suggesting that TEs harbor and potentially evolutionarily select short sequences to act as suppressors of host TE repression.

Is junk DNA bunk? A critique of ENCODE

Do data from the Encyclopedia Of DNA Elements (ENCODE) project render the notion of junk DNA obsolete? Here, I review older arguments for junk grounded in the C-value paradox and propose a thought experiment to challenge ENCODE's ontology. Specifically, what would we expect for the number of functional elements (as ENCODE defines them) in genomes much larger than our own genome? If the number were to stay more or less constant, it would seem sensible to consider the rest of the DNA of larger genomes to be junk or, at least, assign it a different sort of role (structural rather than informational). If, however, the number of functional elements were to rise significantly with C-value then, (i) organisms with genomes larger than our genome are more complex phenotypically than we are, (ii) ENCODE's definition of functional element identifies many sites that would not be considered functional or phenotype-determining by standard uses in biology, or (iii) the same phenotypic functions are often determined in a more diffuse fashion in larger-genomed organisms. Good cases can be made for propositions ii and iii. A larger theoretical framework, embracing informational and structural roles for DNA, neutral as well as adaptive causes of complexity, and selection as a multilevel phenomenon, is needed.

Common genetic variants in Wnt signaling pathway genes as potential prognostic biomarkers for colorectal cancer

Compelling evidence has implicated the Wnt signaling pathway in the pathogenesis of colorectal cancer. We assessed the use of tag single nucleotide polymorphisms (tSNPs) in adenomatous polyposis coli (APC)/β-catenin (CTNNB1) genes to predict outcomes in patients with colorectal cancer. We selected and genotyped 10 tSNP to predict common variants across entire APC and CTNNB1 genes in 282 colorectal cancer patients. The associations of these tSNPs with distant metastasis-free survival and overall survival were evaluated by Kaplan-Meier analysis, Cox regression model, and survival tree analysis. The 5-year overall survival rate was 68.3%. Survival tree analysis identified a higher-order genetic interaction profile consisting of the APC rs565453, CTNNB1 2293303, and APC rs1816769 that was significantly associated with overall survival. The 5-year survival overall rates were 89.2%, 66.1%, and 58.8% for the low-, medium-, and high-risk genetic profiles, respectively (log-rank P = 0.001). After adjusting for possible confounders, including age, gender, carcinoembryonic antigen levels, tumor differentiation, stage, lymphovascular invasion, perineural invasion, and lymph node involvement, the genetic interaction profile remained significant. None of the studied SNPs were individually associated with distant metastasis-free survival and overall survival. Our results suggest that the genetic interaction profile among Wnt pathway SNPs might potentially increase the prognostic value in outcome prediction for colorectal cancer.

De novo sequence assembly and characterisation of a partial transcriptome for an evolutionarily distinct reptile, the tuatara (Sphenodon punctatus)

Background

The tuatara (Sphenodon punctatus) is a species of extraordinary zoological interest, being the only surviving member of an entire order of reptiles which diverged early in amniote evolution. In addition to their unique phylogenetic placement, many aspects of tuatara biology, including temperature-dependent sex determination, cold adaptation and extreme longevity have the potential to inform studies of genome evolution and development. Despite increasing interest in the tuatara genome, genomic resources for the species are still very limited. We aimed to address this by assembling a transcriptome for tuatara from an early-stage embryo, which will provide a resource for genome annotation, molecular marker development and studies of development and adaptation in tuatara.

Results

We obtained 30 million paired-end 50 bp reads from an Illumina Genome Analyzer and assembled them with Velvet and Oases using a range of kmers. After removing redundancy and filtering out low quality transcripts, our transcriptome dataset contained 32911 transcripts, with an N50 of 675 and a mean length of 451 bp. Almost 50% (15965) of these transcripts could be annotated by comparison with the NCBI non-redundant (NR) protein database or the chicken, green anole and zebrafish UniGene sets. A scan of candidate genes and repetitive elements revealed genes involved in immune function, sex differentiation and temperature-sensitivity, as well as over 200 microsatellite markers.

Conclusions

This dataset represents a major increase in genomic resources for the tuatara, increasing the number of annotated gene sequences from just 60 to almost 16,000. This will facilitate future research in sex determination, genome evolution, local adaptation and population genetics of tuatara, as well as inform studies on amniote evolution.

Restless genomes humans as a model organism for understanding host-retrotransposable element dynamics

Since their initial discovery in maize, there have been various attempts to categorize the relationship between transposable elements (TEs) and their host organisms. These have ranged from TEs being selfish parasites to their role as essential, functional components of organismal biology. Research over the past several decades has, in many respects, only served to complicate the issue even further. On the one hand, investigators have amassed substantial evidence concerning the negative effects that TE-mutagenic activity can have on host genomes and organismal fitness. On the other hand, we find an increasing number of examples, across several taxa, of TEs being incorporated into functional biological roles for their host organism. Some 45% of our own genomes are comprised of TE copies. While many of these copies are dormant, having lost their ability to mobilize, several lineages continue to actively proliferate in modern human populations. With its complement of ancestral and active TEs, the human genome exhibits key aspects of the host-TE dynamic that has played out since early on in organismal evolution. In this review, we examine what insights the particularly well-characterized human system can provide regarding the nature of the host-TE interaction.