György Abrusán

TEclass—a tool for automated classification of unknown eukaryotic transposable elements

MOTIVATION The large number of sequenced genomes required the development of software that reconstructs the consensus sequences of transposons and other repetitive elements. However, the available tools usually focus on the accurate identification of raw repeats and provide no information about the taxonomic position of the reconstructed consensi. TEclass is a tool to classify unknown transposable elements into their four main functional categories, which reflect their mode of transposition: DNA transposons, long terminal repeats (LTRs), long interspersed nuclear elements (LINEs) and short interspersed nuclear elements (SINEs). TEclass uses machine learning support vector machine (SVM) for classification based on oligomer frequencies. It achieves 90-97% accuracy in the classification of novel DNA and LTR repeats, and 75% for LINEs and SINEs. AVAILABILITY http://www.compgen.uni-muenster.de/teclass, stand alone program upon request.

Co-localization of CENP-C and CENP-H to discontinuous domains of CENP-A chromatin at human neocentromeres

BackgroundMammalian centromere formation is dependent on chromatin that contains centromere protein (CENP)-A, which is the centromere-specific histone H3 variant. Human neocentromeres have acquired CENP-A chromatin epigenetically in ectopic chromosomal locations on low-copy complex DNA. Neocentromeres permit detailed investigation of centromeric chromatin organization that is not possible in the highly repetitive alpha satellite DNA present at endogenous centromeres.ResultsWe have examined the distribution of CENP-A, as well as two additional centromeric chromatin-associated proteins (CENP-C and CENP-H), across neocentromeric DNA using chromatin immunoprecipitation (ChIP) on CHIP assays on custom genomic microarrays at three different resolutions. Analysis of two neocentromeres using a contiguous bacterial artificial chromosome (BAC) microarray spanning bands 13q31.3 to 13q33.1 shows that both CENP-C and CENP-H co-localize to the CENP-A chromatin domain. Using a higher resolution polymerase chain reaction (PCR)-amplicon microarray spanning the neocentromere, we find that the CENP-A chromatin is discontinuous, consisting of a major domain of about 87.8 kilobases (kb) and a minor domain of about 13.2 kb, separated by an approximately 158 kb region devoid of CENPs. Both CENP-A domains exhibit co-localization of CENP-C and CENP-H, defining a distinct inner kinetochore chromatin structure that is consistent with higher order chromatin looping models at centromeres. The PCR microarray data suggested varying density of CENP-A nucleosomes across the major domain, which was confirmed using a higher resolution oligo-based microarray.ConclusionCentromeric chromatin consists of several CENP-A subdomains with highly discontinuous CENP-A chromatin at both the level of individual nucleosomes and at higher order chromatin levels, raising questions regarding the overall structure of centromeric chromatin.

Filamentous cyanobacteria, temperature and Daphnia growth: the role of fluid mechanics

Viscosity increases significantly with a fall in water temperature, thus temperature change affects not only the metabolic rates of aquatic suspension feeders, but also the physical properties of the surrounding fluid. This mechanistic effect of water temperature change on growth was separated from the effect of metabolism by using culture media with modified viscosity, while the temperature was kept constant. The effect of water viscosity on growth rate and feeding of four Daphnia species (D. magna, D. pulicaria, D. hyalina, D. galeata) was investigated. Increased viscosity decreased the growth rate significantly for three species, with the exception of D. galeata. Changing viscosity also affects growth qualitatively: the filamentous blue-green Cylindrospermopsis raciborskii reduces the growth rate of D. pulicaria at low viscosity, but its negative effect disappears when viscosity is higher. The findings are consistent with the hypothesis that it is the Reynolds number of the filtering appendages that determines the qualitative features of Daphnia filtration. The edibility of C. raciborskii at high water viscosity is most probably caused by lack of interference with filtering combs, and explains the coexistence of D. pulicaria with filamentous blue-green species in the field, and also the observed temperature dependence of growth inhibition of filaments.

The distribution of L1 and Alu retroelements in relation to GC content on human sex chromosomes is consistent with the ectopic recombination model

The distribution of Alu and L1 retroelements in the human genome changes with their age. Active retroelements target AT-rich regions, but their frequency increases in GC- and gene-rich regions of the genome with increasing age of the insertions. Currently there is no consensus on the mechanism generating this pattern. In this paper we test the hypothesis that selection against deleterious deletions caused by ectopic recombination between repeats is the main cause of the inhomogeneous distribution of L1s and Alus, by means of a detailed analysis of the GC distribution of the repeats on the sex chromosomes. We show that (1) unlike on the autosomes and X chromosome, L1s do not accumulate on the Y chromosome in GC-rich regions, whereas Alus accumulate there to a minor extent; (2) on the Y chromosome Alu and L1 densities are positively correlated, unlike the negative correlation on other chromosomes; and (3) in gene-poor regions of chromosome 4 and X, the distribution of Alus and L1s does not shift toward GC-rich regions. In addition, we show that although local GC content of long L1 insertions is lower than average, their selective loss from recombining chromosomes is not the main cause of the enrichment of ancient L1s in GC-rich regions. The results support the hypothesis that ectopic recombination causes the shift of Alu and L1 distributions toward the gene-rich regions of the genome.

Biased Distributions and Decay of Long Interspersed Nuclear Elements in the Chicken Genome

The genomes of birds are much smaller than mammalian genomes, and transposable elements (TEs) make up only 10% of the chicken genome, compared with the 45% of the human genome. To study the mechanisms that constrain the copy numbers of TEs, and as a consequence the genome size of birds, we analyzed the distributions of LINEs (CR1s) and SINEs (MIRs) on the chicken autosomes and Z chromosome. We show that (1) CR1 repeats are longest on the Z chromosome and their length is negatively correlated with the local GC content; (2) the decay of CR1 elements is highly biased, and the 5′-ends of the insertions are lost much faster than their 3′-ends; (3) the GC distribution of CR1 repeats shows a bimodal pattern with repeats enriched in both AT-rich and GC-rich regions of the genome, but the CR1 families show large differences in their GC distribution; and (4) the few MIRs in the chicken are most abundant in regions with intermediate GC content. Our results indicate that the primary mechanism that removes repeats from the chicken genome is ectopic exchange and that the low abundance of repeats in avian genomes is likely to be the consequence of their high recombination rates.

Integration of New Genes into Cellular Networks, and Their Structural Maturation

It has been recently discovered that new genes can originate de novo from noncoding DNA, and several biological traits including expression or sequence composition form a continuum from noncoding sequences to conserved genes. In this article, using yeast genes I test whether the integration of new genes into cellular networks and their structural maturation shows such a continuum by analyzing their changes with gene age. I show that 1) The number of regulatory, protein–protein, and genetic interactions increases continuously with gene age, although with very different rates. New regulatory interactions emerge rapidly within a few million years, while the number of protein–protein and genetic interactions increases slowly, with a rate of 2–2.25 × 10−8/year and 4.8 × 10−8/year, respectively. 2) Gene essentiality evolves relatively quickly: the youngest essential genes appear in proto-genes ∼14 MY old. 3) In contrast to interactions, the secondary structure of proteins and their robustness to mutations indicate that new genes face a bottleneck in their evolution: proto-genes are characterized by high β-strand content, high aggregation propensity, and low robustness against mutations, while conserved genes are characterized by lower strand content and higher stability, most likely due to the higher probability of gene loss among young genes and accumulation of neutral mutations.

Analysis of Transposon Interruptions Suggests Selection for L1 Elements on the X Chromosome

It has been hypothesised that the massive accumulation of L1 transposable elements on the X chromosome is due to their function in X inactivation, and that the accumulation of Alu elements near genes is adaptive. We tested the possible selective advantage of these two transposable element (TE) families with a novel method, interruption analysis. In mammalian genomes, a large number of TEs interrupt other TEs due to the high overall abundance and age of repeats, and these interruptions can be used to test whether TEs are selectively neutral. Interruptions of TEs, which are beneficial for the host, are expected to be deleterious and underrepresented compared with neutral ones. We found that L1 elements in the regions of the X chromosome that contain the majority of the inactivated genes are significantly less frequently interrupted than on the autosomes, while L1s near genes that escape inactivation are interrupted with higher frequency, supporting the hypothesis that L1s on the X chromosome play a role in its inactivation. In addition, we show that TEs are less frequently interrupted in introns than in intergenic regions, probably due to selection against the expansion of introns, but the insertion pattern of Alus is comparable to other repeats.

Turning gold into 'junk': transposable elements utilize central proteins of cellular networks

The numerous discovered cases of domesticated transposable element (TE) proteins led to the recognition that TEs are a significant source of evolutionary innovation. However, much less is known about the reverse process, whether and to what degree the evolution of TEs is influenced by the genome of their hosts. We addressed this issue by searching for cases of incorporation of host genes into the sequence of TEs and examined the systems-level properties of these genes using the Saccharomyces cerevisiae and Drosophila melanogaster genomes. We identified 51 cases where the evolutionary scenario was the incorporation of a host gene fragment into a TE consensus sequence, and we show that both the yeast and fly homologues of the incorporated protein sequences have central positions in the cellular networks. An analysis of selective pressure (Ka/Ks ratio) detected significant selection in 37% of the cases. Recent research on retrovirus-host interactions shows that virus proteins preferentially target hubs of the host interaction networks enabling them to take over the host cell using only a few proteins. We propose that TEs face a similar evolutionary pressure to evolve proteins with high interacting capacities and take some of the necessary protein domains directly from their hosts.

Somatic transposition in the brain has the potential to influence the biosynthesis of metabolites involved in Parkinson’s disease and schizophrenia

It has been recently discovered that transposable elements show high activity in the brain of mammals, however, the magnitude of their influence on its functioning is unclear so far. In this paper, I use flux balance analysis to examine the influence of somatic retrotransposition on brain metabolism, and the biosynthesis of its key metabolites, including neurotransmitters. The analysis shows that somatic transposition in the human brain can influence the biosynthesis of more than 250 metabolites, including dopamine, serotonin and glutamate, shows large inter-individual variability in metabolic effects, and may contribute to the development of Parkinson’s disease and schizophrenia.ReviewersThis article was reviewed by Dr Kenji Kojima (nominated by Dr Jerzy Jurka) and Dr Eugene Koonin.

Alpha helices are more robust to mutations than beta strands

The rapidly increasing amount of data on human genetic variation has resulted in a growing demand to identify pathogenic mutations computationally, as their experimental validation is currently beyond reach. Here we show that alpha helices and beta strands differ significantly in their ability to tolerate mutations: helices can accumulate more mutations than strands without change, due to the higher numbers of inter-residue contacts in helices. This results in two patterns: a) the same number of mutations causes less structural change in helices than in strands; b) helices diverge more rapidly in sequence than strands within the same domains. Additionally, both helices and strands are significantly more robust than coils. Based on this observation we show that human missense mutations that change secondary structure are more likely to be pathogenic than those that do not. Moreover, inclusion of predicted secondary structure changes shows significant utility for improving upon state-of-the-art pathogenicity predictions.