Search

Homologous recombination is a key process that governs the stability of eukaryotic genomes during DNA replication and repair. Multiple auxiliary factors regulate the choice of homologous recombination pathway in response to different types of replication stress. Using Schizosaccharomyces pombe we have previously suggested the role of DNA translocases Rrp1 and Rrp2, together with Srs2 helicase, in the common synthesis-dependent strand annealing sub-pathway of homologous recombination. Here we show that all three proteins are important for completion of replication after hydroxyurea exposure and provide data comparing the effect of overproduction of Srs2 with Rrp1 and Rrp2. We demonstrate that Srs2 localises to rDNA region and is required for proper replication of rDNA arrays. Upregulation of Srs2 protein levels leads to enhanced replication stress, chromosome instability and viability loss, as previously reported for Rrp1 and Rrp2. Interestingly, our data suggests that dysregulation of Srs2, Rrp1 and Rrp2 protein levels differentially affects checkpoint response: overproduction of Srs2 activates simultaneously DNA damage and replication stress response checkpoints, while cells overproducing Rrp1 mainly launch DNA damage checkpoint. On the other hand, upregulation of Rrp2 primarily leads to replication stress response checkpoint activation. Overall, we propose that Srs2, Rrp1 and Rrp2 have important and at least partially independent functions in the maintenance of distinct difficult to replicate regions of the genome. [ABSTRACT FROM AUTHOR]

Recent advances in long-read sequencing technologies have enabled the complete assembly of eukaryotic genomes from telomere to telomere by allowing repeated regions to be fully sequenced and assembled, thus filling the gaps left by previous short-read sequencing methods. Furthermore, long-read sequencing can also help characterizing structural variants, with applications in the fields of genome evolution or cancer genomics. For many organisms, the main bottleneck to sequence long reads remains the lack of robust methods to obtain high-molecular-weight (HMW) DNA. For this purpose, we developed an optimized protocol to extract DNA suitable for long-read sequencing from the unicellular green alga Chlamydomonas reinhardtii, based on CTAB/phenol extraction followed by a size selection step for long DNA molecules. We provide validation results for the extraction protocol, as well as statistics obtained with Oxford Nanopore Technologies sequencing. [ABSTRACT FROM AUTHOR]

Analysis of eukaryotic genomes requires the detection and classification of transposable elements (TEs), a crucial but complex and time-consuming task. To improve the performance of tools that accomplish these tasks, Machine Learning approaches (ML) that leverage computer resources, such as GPUs (Graphical Processing Unit) and multiple CPU (Central Processing Unit) cores, have been adopted. However, until now, the use of ML techniques has mostly been limited to classification of TEs. Herein, a detection-classification strategy (named YORO) based on convolutional neural networks is adapted from computer vision (YOLO) to genomics. This approach enables the detection of genomic objects through the prediction of the position, length, and classification in large DNA sequences such as fully sequenced genomes. As a proof of concept, the internal protein-coding domains of LTR-retrotransposons are used to train the proposed neural network. Precision, recall, accuracy, F1-score, execution times and time ratios, as well as several graphical representations were used as metrics to measure performance. These promising results open the door for a new generation of Deep Learning tools for genomics. YORO architecture is available at https://github.com/simonorozcoarias/YORO. [ABSTRACT FROM AUTHOR]

Rapid advancements in automated genomic technologies have uncovered many unique findings about the turtle genome and its associated features including olfactory gene expansions and duplications of toll-like receptors. However, despite the advent of large-scale sequencing, assembly, and annotation, about 40–50% of genes in eukaryotic genomes are left without functional annotation, severely limiting our knowledge of the biological information of genes. Additionally, these automated processes are prone to errors since draft genomes consist of several disconnected scaffolds whose order is unknown; erroneous draft assemblies may also be contaminated with foreign sequences and propagate to cause errors in annotation. Many of these automated annotations are thus incomplete and inaccurate, highlighting the need for functional annotation to link gene sequences to biological identity. In this study, we have functionally annotated two genes of the red-bellied short-neck turtle (Emydura subglobosa), a member of the relatively understudied pleurodire lineage of turtles. We improved upon initial ab initio gene predictions through homology-based evidence and generated refined consensus gene models. Through functional, localization, and structural analyses of the predicted proteins, we discovered conserved putative genes encoding mitochondrial proteins that play a role in C21-steroid hormone biosynthetic processes and fatty acid catabolism—both of which are distantly related by the tricarboxylic acid (TCA) cycle and share similar metabolic pathways. Overall, these findings further our knowledge about the genetic features underlying turtle physiology, morphology, and longevity, which have important implications for the treatment of human diseases and evolutionary studies. [ABSTRACT FROM AUTHOR]

Phylogenetic profiling in eukaryotes is of continued interest to study and predict the functional relationships between proteins. This interest is likely driven by the increased number of available diverse genomes and computational methods to infer orthologies. The evaluation of phylogenetic profiles has mainly focussed on reference genome selection in prokaryotes. However, it has been proven to be challenging to obtain high prediction accuracies in eukaryotes. As part of our recent comparison of orthology inference methods for eukaryotic genomes, we observed a surprisingly high performance for predicting interacting orthologous groups. This high performance, in turn, prompted the question of what factors influence the success of phylogenetic profiling when applied to eukaryotic genomes. Here we analyse the effect of species, orthologous group and interactome selection on protein interaction prediction using phylogenetic profiles. We select species based on the diversity and quality of the genomes and compare this supervised selection with randomly generated genome subsets. We also analyse the effect on the performance of orthologous groups defined to be in the last eukaryotic common ancestor of eukaryotes to that of orthologous groups that are not. Finally, we consider the effects of reference interactome set filtering and reference interactome species. In agreement with other studies, we find an effect of genome selection based on quality, less of an effect based on genome diversity, but a more notable effect based on the amount of information contained within the genomes. Most importantly, we find it is not merely selecting the correct genomes that is important for high prediction performance. Other choices in meta parameters such as orthologous group selection, the reference species of the interaction set, and the quality of the interaction set have a much larger impact on the performance when predicting protein interactions using phylogenetic profiles. These findings shed light on the differences in reported performance amongst phylogenetic profiles approaches, and reveal on a more fundamental level for which types of protein interactions this method has most promise when applied to eukaryotes. [ABSTRACT FROM AUTHOR]

Pangenome analyses are increasingly being utilized to study the evolution of eukaryotic organisms. While pangenomes can provide insight into polymorphic gene content, inferences about the ecological and adaptive potential of such organisms also need to be accompanied by additional supportive genomic analyses. In this study we constructed a pangenome of Claviceps purpurea from 24 genomes and examined the positive selection and recombination landscape of an economically important fungal organism for pharmacology and agricultural research. Together, these analyses revealed that C. purpurea has a relatively large accessory genome (~ 38%), high recombination rates (ρ = 0.044), and transposon mediated gene duplication. However, due to observations of relatively low transposable element (TE) content (8.8%) and a lack of variability in genome sizes, prolific TE expansion may be controlled by frequent recombination. We additionally identified that within the ergoline biosynthetic cluster the lpsA1 and lpsA2 were the result of a recombination event. However, the high recombination rates observed in C. purpurea may be influencing an overall trend of purifying selection across the genome. These results showcase the use of selection and recombination landscapes to identify mechanisms contributing to pangenome structure and primary factors influencing the evolution of an organism. [ABSTRACT FROM AUTHOR]

The cellular energy metabolism is one of the most conserved processes, as it is present in all living organisms. Mitochondria are providing the eukaryotic cell with energy and thus their genome and gene expression has been of broad interest for a long time. Mitochondrial gene expression changes under different conditions and is regulated by genes encoded in the nucleus of the cell. In this context, little is known about non-model organisms and we provide the first large-scaled gene expression analysis of mitochondrial-linked genes in laying hens. We analysed 28 mitochondrial and nuclear genes in 100 individuals in the context of five life-stages and strain differences among five tissues. Our study showed that mitochondrial gene expression increases during the productive life span, and reacts tissue and strain specific. In addition, the strains react different to potential increased oxidative stress, resulting from the increase in mitochondrial gene expression. The results suggest that the cellular energy metabolism as part of a complex regulatory system is strongly affected by the productive life span in laying hens and thus partly comparable to model organisms. This study provides a starting point for further analyses in this field on non-model organisms, especially in laying-hens. [ABSTRACT FROM AUTHOR]

In every kingdom of life, GC->AT transitions occur more frequently than any other type of mutation due to the spontaneous deamination of cytidine. In eukaryotic genomes, this slow loss of GC base pairs is counteracted by biased gene conversion which increases genomic GC content as part of the recombination process. However, this type of biased gene conversion has not been observed in bacterial genomes, so we hypothesized that GC->AT transitions cause a reduction of genomic GC content in prokaryotic genomes on an evolutionary time scale. To test this hypothesis, we used a phylogenetic approach to analyze triplets of closely related genomes representing a wide range of the bacterial kingdom. The resulting data indicate that genomic GC content is drifting downward in bacterial genomes where GC base pairs comprise 40% or more of the total genome. In contrast, genomes containing less than 40% GC base pairs have fewer opportunities for GC->AT transitions to occur so genomic GC content is relatively stable or actually increasing. It should be noted that this observed change in genomic GC content is the net change in shared parts of the genome and does not apply to parts of the genome that have been lost or acquired since the genomes being compared shared common ancestor. However, a more detailed analysis of two Caulobacter genomes revealed that the acquisition of mobile elements by the two genomes actually reduced the total genomic GC content as well. [ABSTRACT FROM AUTHOR]

We present a model of the evolution of control systems in a genome under environmental constraints. The model conceptually follows the Jacob and Monod model of gene control. Genes contain control elements which respond to the internal state of the cell as well as the environment to control expression of a coding region. Control and coding regions evolve to maximize a fitness function between expressed coding sequences and the environment. The model was run 118 times to an average of 1.4∙106 'generations' each with a range of starting parameters probed the conditions under which genomes evolved a 'default style' of control. Unexpectedly, the control logic that evolved was not significantly correlated to the complexity of the environment. Genetic logic was strongly correlated with genome complexity and with the fraction of genes active in the cell at any one time. More complex genomes correlated with the evolution of genetic controls in which genes were active ('default on'), and a low fraction of genes being expressed correlated with a genetic logic in which genes were biased to being inactive unless positively activated ('default off' logic). We discuss how this might relate to the evolution of the complex eukaryotic genome, which operates in a 'default off' mode. [ABSTRACT FROM AUTHOR]

Pantoea agglomerans is a Gram-negative facultative anaerobic bacillus causing a wide range of opportunistic infections in humans including septicemia, pneumonia, septic arthritis, wound infections and meningitis. To date, the determinants of virulence, antibiotic resistance, metabolic features conferring survival and host-associated pathogenic potential of this bacterium remain largely underexplored. In this study, we sequenced and assembled the whole-genome of P. agglomerans KM1 isolated from kimchi in South Korea. The genome contained one circular chromosome of 4,039,945 bp, 3 mega plasmids, and 2 prophages. The phage-derived genes encoded integrase, lysozyme and terminase. Six CRISPR loci were identified within the bacterial chromosome. Further in-depth analysis showed that the genome contained 13 antibiotic resistance genes conferring resistance to clinically important antibiotics such as penicillin G, bacitracin, rifampicin, vancomycin, and fosfomycin. Genes involved in adaptations to environmental stress were also identified which included factors providing resistance to osmotic lysis, oxidative stress, as well as heat and cold shock. The genomic analysis of virulence factors led to identification of a type VI secretion system, hemolysin, filamentous hemagglutinin, and genes involved in iron uptake and sequestration. Finally, the data provided here show that, the KM1 isolate exerted strong immunostimulatory properties on RAW 264.7 macrophages in vitro. Stimulated cells produced Nitric Oxide (NO) and pro-inflammatory cytokines TNF-α, IL-6 and the anti-inflammatory cytokine IL-10. The upstream signaling for production of TNF-α, IL-6, IL-10, and NO depended on TLR4 and TLR1/2. While production of TNF-α, IL-6 and NO involved solely activation of the NF-κB, IL-10 secretion was largely dependent on NF-κB and to a lesser extent on MAPK Kinases. Taken together, the analysis of the whole-genome and immunostimulatory properties provided in-depth characterization of the P. agglomerans KM1 isolate shedding a new light on determinants of virulence that drive its interactions with the environment, other microorganisms and eukaryotic hosts [ABSTRACT FROM AUTHOR]

Functional analyses of genes are crucial for unveiling biological responses, genetic engineering, and developing new medicines. However, functional analyses have largely been restricted to model organisms, representing a major hurdle for functional studies and industrial applications. To resolve this, comparative genome analyses can be used to provide clues to gene functions as well as their evolutionary history. To this end, we present Prometheus, a web-based omics portal that contains more than 17,215 sequences from prokaryotic and eukaryotic genomes. This portal supports interkingdom comparative analyses via a domain architecture-based gene identification system and Gene Search, and users can easily and rapidly identify single or entire gene sets in specific pathways. Bioinformatics tools for further analyses are provided in Prometheus or through Bio-Express, a cloud-based bioinformatics analysis platform. Prometheus is a new paradigm for comparative analyses of large amounts of genomic information. [ABSTRACT FROM AUTHOR]

Transposable elements (TEs) are mobile genetic elements in eukaryotic genomes. Recent research highlights the important role of TEs in the embryogenesis, neurodevelopment, and immune functions. However, there is a lack of a one-stop and easy to use computational pipeline for expression analysis of both genes and locus-specific TEs from RNA-Seq data. Here, we present GeneTEFlow, a fully automated, reproducible and platform-independent workflow, for the comprehensive analysis of gene and locus-specific TEs expression from RNA-Seq data employing Nextflow and Docker technologies. This application will help researchers more easily perform integrated analysis of both gene and TEs expression, leading to a better understanding of roles of gene and TEs regulation in human diseases. GeneTEFlow is freely available at https://github.com/zhongw2/GeneTEFlow. [ABSTRACT FROM AUTHOR]

Background: Translationally controlled tumor protein (TCTP) is a conserved, multifunctional protein involved in numerous cellular processes in eukaryotes. Although the functions of TCTP have been investigated sporadically in animals, invertebrates, and plants, few lineage-specific activities of this molecule, have been reported. An exception is in Arabidopsis thaliana, in which TCTP (AtTCTP1) functions in stomatal closuer by regulating microtubule stability. Further, although the development of next-generation sequencing technologies has facilitated the analysis of many eukaryotic genomes in public databases, inter-kingdom comparative analyses using available genome information are comparatively scarce. Methodology: To carry out inter-kingdom comparative analysis of TCTP, TCTP genes were identified from 377 species. Then phylogenetic analysis, prediction of protein structure, molecular docking simulation and molecular dynamics analysis were performed to investigate the evolution of TCTP genes and their binding proteins. Results: A total of 533 TCTP genes were identified from 377 eukaryotic species, including protozoa, fungi, invertebrates, vertebrates, and plants. Phylogenetic and secondary structure analyses reveal lineage-specific evolution of TCTP, and inter-kingdom comparisons highlight the lineage-specific emergence of, or changes in, secondary structure elements in TCTP proteins from different kingdoms. Furthermore, secondary structure comparisons between TCTP proteins within each kingdom, combined with measurements of the degree of sequence conservation, suggest that TCTP genes have evolved to conserve protein secondary structures in a lineage-specific manner. Additional tertiary structure analysis of TCTP-binding proteins and their interacting partners and docking simulations between these proteins further imply that TCTP gene variation may influence the tertiary structures of TCTP-binding proteins in a lineage-specific manner. Conclusions: Our analysis suggests that TCTP has undergone lineage-specific evolution and that structural changes in TCTP proteins may correlate with the tertiary structure of TCTP-binding proteins and their binding partners in a lineage-specific manner. [ABSTRACT FROM AUTHOR]

Peony is a famous ornamental and medicinal plant in China, and peony hybrid breeding is an important means of germplasm innovation. However, research on the genome of this species is limited, thereby hindering the genetic and breeding research on peony. In the present study, simple sequence repeat (SSR) locus analysis was performed on expressed sequence tags obtained by the transcriptome sequencing of Paeonia using Microsatellite software. Primers with polymorphism were obtained via polymerase chain reaction amplification and electrophoresis. As a result, a total of 86,195 unigenes were obtained by assembling the transcriptome data of Paeonia. Functional annotations were obtained in seven functional databases including 49,172 (Non-Redundant Protein Sequence Database: 57.05%), 38,352 (Nucleotide Sequence Database: 44.49%), 36,477 (Swiss Prot: 42.32%), 38,905 (Clusters of Orthologous Groups for Eukaryotic Complete Genomes: 45.14%), 37,993 (Kyoto Encyclopedia of Genes and Genomes: 44.08%), 26,832 (Gene Ontology: 31.13%) and 37,758 (Pfam: 43.81%) unigenes. Meanwhile, 21,998 SSR loci were distributed in 17,567 unigenes containing SSR sequences, and the SSR distribution frequency was 25.52%, with an average of one SSR sequence per 4.66 kb. Mononucleotide, dinucleotide, and trinucleotide were the main repeat types, accounting for 55.74%, 25.58%, and 13.21% of the total repeat times, respectively. Forty-five pairs of the 100 pairs of primers selected randomly could amplify clear polymorphic bands. The polymorphic primers of these 45 pairs were used to cluster and analyze 16 species of peony. The new SSR molecular markers can be useful for the study of genetic diversity and marker-assisted breeding of peony. [ABSTRACT FROM AUTHOR]

The sophistication of gene prediction algorithms and the abundance of RNA-based evidence for the maize genome may suggest that manual curation of gene models is no longer necessary. However, quality metrics generated by the MAKER-P gene annotation pipeline identified 17,225 of 130,330 (13%) protein-coding transcripts in the B73 Reference Genome V4 gene set with models of low concordance to available biological evidence. Working with eight graduate students, we used the Apollo annotation editor to curate 86 transcript models flagged by quality metrics and a complimentary method using the Gramene gene tree visualizer. All of the triaged models had significant errors–including missing or extra exons, non-canonical splice sites, and incorrect UTRs. A correct transcript model existed for about 60% of genes (or transcripts) flagged by quality metrics; we attribute this to the convention of elevating the transcript with the longest coding sequence (CDS) to the canonical, or first, position. The remaining 40% of flagged genes resulted in novel annotations and represent a manual curation space of about 10% of the maize genome (~4,000 protein-coding genes). MAKER-P metrics have a specificity of 100%, and a sensitivity of 85%; the gene tree visualizer has a specificity of 100%. Together with the Apollo graphical editor, our double triage provides an infrastructure to support the community curation of eukaryotic genomes by scientists, students, and potentially even citizen scientists. [ABSTRACT FROM AUTHOR]

Eukaryotic translation initiation factor 2-alpha kinase (EIF2AK) proteins inhibit protein synthesis at translation initiation level, in response to various stress conditions, including oxidative stress, heme deficiency, osmotic shock, and heat shock. Origin and functional diversification of EIF2AK sequences remain ambiguous. Here we determine the origin and molecular evolution of EIF2AK proteins in lower eukaryotes and studied the molecular basis of divergence among sub-family sequences. Present work emphasized primitive origin of EIF2AK4 sub-family gene in lower eukaryotes of protozoan lineage. Phylogenetic analysis supported common origin and sub-family based classification of EIF2AKs. Functional divergence studies across sub-families revealed several putative amino acid sites, which assist in altered protein interactions of kinase domains. The data can facilitate designing site-directed experimental studies aiming at elucidating diverse functional aspects of kinase domains regarding down-regulation of protein synthesis. [ABSTRACT FROM AUTHOR]

Drawing on a long history in macroecology, correlation analysis of microbiome datasets is becoming a common practice for identifying relationships or shared ecological niches among bacterial taxa. However, many of the statistical issues that plague such analyses in macroscale communities remain unresolved for microbial communities. Here, we discuss problems in the analysis of microbial species correlations based on presence-absence data. We focus on presence-absence data because this information is more readily obtainable from sequencing studies, especially for whole-genome sequencing, where abundance estimation is still in its infancy. First, we show how Pearson’s correlation coefficient (r) and Jaccard’s index (J)–two of the most common metrics for correlation analysis of presence-absence data–can contradict each other when applied to a typical microbiome dataset. In our dataset, for example, 14% of species-pairs predicted to be significantly correlated by r were not predicted to be significantly correlated using J, while 37.4% of species-pairs predicted to be significantly correlated by J were not predicted to be significantly correlated using r. Mismatch was particularly common among species-pairs with at least one rare species (<10% prevalence), explaining why r and J might differ more strongly in microbiome datasets, where there are large numbers of rare taxa. Indeed 74% of all species-pairs in our study had at least one rare species. Next, we show how Pearson’s correlation coefficient can result in artificial inflation of positive taxon relationships and how this is a particular problem for microbiome studies. We then illustrate how Jaccard’s index of similarity (J) can yield improvements over Pearson’s correlation coefficient. However, the standard null model for Jaccard’s index is flawed, and thus introduces its own set of spurious conclusions. We thus identify a better null model based on a hypergeometric distribution, which appropriately corrects for species prevalence. This model is available from recent statistics literature, and can be used for evaluating the significance of any value of an empirically observed Jaccard’s index. The resulting simple, yet effective method for handling correlation analysis of microbial presence-absence datasets provides a robust means of testing and finding relationships and/or shared environmental responses among microbial taxa. [ABSTRACT FROM AUTHOR]

Background: DNA methylation is one of the most important epigenetic modifications of eukaryotic genomes and is believed to play integral roles in diverse biological processes. Although DNA methylation has been well studied in mammals, data are limited in invertebrates, particularly Mollusca. The Pacific oyster Crassostrea gigas is an emerging genetic model for functional analysis of DNA methylation in Mollusca. Recent studies have shown that there is a positive association between methylation status and gene expression in C. gigas; however, whether this association exists at the exon level remains to be determined. Results: In this study, we characterized the genome-wide methylation pattern across two different tissues of C. gigas and found that methylated genes are expressed in more tissues and development stages than unmethylated genes. Furthermore, we found that different types of exons had different methylation levels, with the lowest methylation levels in the first exons, followed by the last exons, and the internal exons. We found that the exons included in the gene transcript contained significantly higher DNA methylation levels than skipped exons. We observed that the DNA methylation levels increased slowly after the start sites and end sites of exons seperately, and then decreased quickly towards the middle sites of exons. We also found that methylated exons were significantly longer than unmethylated exons. Conclusion: This study constitutes the first genome-wide analysis to show an association between exon-level DNA methylation and mRNA expression in the oyster. Our findings suggest that exon-level DNA methylation may play a role in the construction of alternative splicing by positively influencing exon inclusion during transcription. [ABSTRACT FROM AUTHOR]

Human impact on marine benthic communities has traditionally been assessed using visible morphological traits and has focused on the macrobenthos, whereas the ecologically important organisms of the meio- and microbenthos have received less attention. DNA metabarcoding offers an alternative to this approach and enables a larger fraction of the biodiversity in marine sediments to be monitored in a cost-efficient manner. Although this methodology remains poorly standardised and challenged by biases inherent to rRNA copy number variation, DNA extraction, PCR, and limitations related to taxonomic identification, it has been shown to be semi-quantitative and useful for comparing taxon abundances between samples. Here, we evaluate the effect of replicating genomic DNA extraction in order to counteract small scale spatial heterogeneity and improve diversity and community structure estimates in metabarcoding-based monitoring. For this purpose, we used ten technical replicates from three different marine sediment samples. The effect of sequence depth was also assessed, and in silico pooling of DNA extraction replicates carried out in order to maintain the number of reads constant. Our analyses demonstrated that both sequencing depth and DNA extraction replicates could improve diversity estimates as well as the ability to separate samples with different characteristics. We could not identify a “sufficient” replicate number or sequence depth, where further improvements had a less significant effect. Based on these results, we consider replication an attractive alternative to directly increasing the amount of sample used for DNA extraction and strongly recommend it for future metabarcoding studies and routine assessments of sediment biodiversity. [ABSTRACT FROM AUTHOR]

Retrotransposons with long terminal repeats (LTR-RTs) are widespread mobile elements in eukaryotic genomes. We obtained a total of 81 partial LTR-RT sequences from lentil corresponding to internal retrotransposon components and LTRs. Sequences were obtained by PCR from genomic DNA. Approximately 37% of the LTR-RT internal sequences presented premature stop codons, pointing out that these elements must be non-autonomous. LTR sequences were obtained using the iPBS technique which amplifies sequences between LTR-RTs. A total of 193 retrotransposon-derived genetic markers, mainly iPBS, were used to obtain a genetic linkage map from 94 F7 inbred recombinant lines derived from the cross between the cultivar Lupa and the wild ancestor L. culinaris subsp. orientalis. The genetic map included 136 markers located in eight linkage groups. Clusters of tightly linked retrotransposon-derived markers were detected in linkage groups LG1, LG2, and LG6, hence denoting a non-random genomic distribution. Phylogenetic analyses identified the LTR-RT families in which internal and LTR sequences are included. Ty3-gypsy elements were more frequent than Ty1-copia, mainly due to the high Ogre element frequency in lentil, as also occurs in other species of the tribe Vicieae. LTR and internal sequences were used to analyze in silico their distribution among the contigs of the lentil draft genome. Up to 8.8% of the lentil contigs evidenced the presence of at least one LTR-RT similar sequence. A statistical analysis suggested a non-random distribution of these elements within of the lentil genome. In most cases (between 97% and 72%, depending on the LTR-RT type) none of the internal sequences flanked by the LTR sequence pair was detected, suggesting that defective and non-autonomous LTR-RTs are very frequent in lentil. Results support that LTR-RTs are abundant and widespread throughout of the lentil genome and that they are a suitable source of genetic markers useful to carry out further genetic analyses. [ABSTRACT FROM AUTHOR]

The reliability of molecular phylogenies is strongly dependent on the quality of the assembled datasets. In the case of eukaryotes, the selection of only one protein isoform per genomic locus is mandatory to avoid biases linked to redundancy. Here, we present IsoSel, a tool devoted to the selection of alternative isoforms in the context of phylogenetic reconstruction. It provides a better alternative to the widely used approach consisting in the selection of the longest isoforms and it performs better than Guidance, its only available counterpart. IsoSel is publicly available at . [ABSTRACT FROM AUTHOR]

It is now well established that eukaryote genomes have a common architectural organization into topologically associated domains (TADs) and evidence is accumulating that this organization plays an important role in gene regulation. However, the mechanisms that partition the genome into TADs and the nature of domain boundaries are still poorly understood. We have investigated boundary regions in the Drosophila genome and find that they can be identified as domains of very low H3K27me3. The genome-wide H3K27me3 profile partitions into two states; very low H3K27me3 identifies Depleted (D) domains that contain housekeeping genes and their regulators such as the histone acetyltransferase-containing NSL complex, whereas domains containing moderate-to-high levels of H3K27me3 (Enriched or E domains) are associated with regulated genes, irrespective of whether they are active or inactive. The D domains correlate with the boundaries of TADs and are enriched in a subset of architectural proteins, particularly Chromator, BEAF-32, and Z4/Putzig. However, rather than being clustered at the borders of these domains, these proteins bind throughout the H3K27me3-depleted regions and are much more strongly associated with the transcription start sites of housekeeping genes than with the H3K27me3 domain boundaries. While we have not demonstrated causality, we suggest that the D domain chromatin state, characterised by very low or absent H3K27me3 and established by housekeeping gene regulators, acts to separate topological domains thereby setting up the domain architecture of the genome. [ABSTRACT FROM AUTHOR]

Transposons, or transposable elements, are the major components of genomes in most eukaryotes. Some groups of transposons have developed target specificity that limits the integration sites to a specific nonessential sequence or a genomic region to avoid gene disruption caused by insertion into an essential gene. R2 is one of the most intensively investigated groups of sequence-specific non-LTR retrotransposons and is inserted at a specific site inside of 28S ribosomal RNA (rRNA) genes. R2 is known to be distributed among at least six animal phyla even though its occurrence is reported to be patchy. Here, in order to obtain a more detailed picture of the distribution of R2, we surveyed R2 using both in silico screening and degenerate PCR, particularly focusing on actinopterygian fish. We found two families of the R2C lineage from vertebrates, although it has previously only been found in platyhelminthes. We also revealed the apparent movement of insertion sites of a lineage of actinopterygian R2, which was likely concurrent with the acquisition of a 28S rRNA-derived sequence in their 3′ UTR. Outside of actinopterygian fish, we revealed the maintenance of a single R2 lineage in birds; the co-existence of four lineages of R2 in the leafcutter bee Megachile rotundata; the first examples of R2 in Ctenophora, Mollusca, and Hemichordata; and two families of R2 showing no target specificity. These findings indicate that R2 is relatively stable and universal, while differences in the distribution and maintenance of R2 lineages probably reflect characteristics of some combination of both R2 lineages and host organisms. [ABSTRACT FROM AUTHOR]

Since the density of simple sequence repeats (SSRs) may vary between different chromosomes of the same species in eukaryotic genomes, we screened SSRs of the whole genome of the yellow necked mouse, Apodemus flavicollis, in order to reveal SSR profiles specific for animals carrying B chromosomes. We found that the 2200 bp band was amplified by primer (CAG)4AC to a highly increased level in samples with B chromosomes. This quantitative difference (B-marker) between animals with (+B) and without (0B) B chromosomes was used to screen 20 populations (387 animals). The presence/absence of Bs was confirmed in 96.5% of 342 non mosaic individuals, which recommends this method for noninvasive B-presence detection. A group of 45 animals with mosaic and micro B (μB) karyotypes was considered separately and showed 55.6% of overall congruence between karyotyping and molecular screening results. Relative quantification by qPCR of two different targeted sequences from B-marker indicated that these B-specific fragments are multiplied on B chromosomes. It also confirms our assumption that different types of Bs with variable molecular composition may exist in the same individual and between individuals of this species. Our results substantiate the origin of Bs from the standard chromosomal complement. The B-marker showed 98% sequence identity with the serine/threonine protein kinase VRK1 gene, similarly to findings reported for Bs from phylogenetically highly distant mammalian species. Evolutionarily conserved protein-coding genes found in Bs, including this one in A. flavicollis, could suggest a common evolutionary pathway. [ABSTRACT FROM AUTHOR]

G-quadruplex forming sequences are particularly enriched in the promoter regions of eukaryotic genes, especially of oncogenes. One of the most well studied G-quadruplex forming sequences is located in the nuclease hypersensitive element (NHE) III1 of the c-MYC promoter region. The oncoprotein c-MYC regulates a large array of genes which play important roles in growth regulation and metabolism. It is dysregulated in >70% of human cancers. The silencer NHEIII1 located upstream of the P1 promoter regulates up-to 80% of c-MYC transcription and includes a G-quadruplex structure (Pu27) that is required for promoter inhibition. We have identified, for the first time, a family of seventeen G-quadruplex-forming motifs with >90% identity with Pu27, located on different chromosomes throughout the human genome, some found near or within genes involved in stem cell maintenance or neural cell development. Notably, all members of the Pu27 family interact specifically with NHEIII1 sequence, in vitro. Crosslinking studies demonstrate that Pu27 oligonucleotide binds specifically to the C-rich strand of the NHEIII1 resulting in the G-quadruplex structure stabilization. Pu27 homologous sequences (Pu27-HS) significantly inhibit leukemic cell lines proliferation in culture. Exposure of U937 cells to the Pu27-HS induces cell growth inhibition associated with cell cycle arrest that is most likely due to downregulation of c-MYC expression at the RNA and/or protein levels. Expression of SOX2, another gene containing a Pu27-HS, was affected by Pu27-HS treatment as well. Our data suggest that the oligonucleotides encoding the Pu27 family target complementary DNA sequences in the genome, including those of the c-MYC and SOX2 promoters. This effect is most likely cell type and cell growth condition dependent. The presence of genomic G-quadruplex-forming sequences homologous to Pu27 of c-MYC silencer and the fact that they interact specifically with the parent sequence suggest a common regulatory mechanism for genes whose promoters contain these sequences. [ABSTRACT FROM AUTHOR]

The proper prediction of the gene catalogue of an organism is essential to obtain a representative snapshot of its overall lifestyle, especially when it is not amenable to culturing. Microsporidia are obligate intracellular, sometimes hard to culture, eukaryotic parasites known to infect members of every animal phylum. To date, sequencing and annotation of microsporidian genomes have revealed a poor gene complement with highly reduced gene sizes. In the present paper, we investigated whether such gene sizes may have induced biases for the methodologies used for genome annotation, with an emphasis on small coding sequence (CDS) gene prediction. Using better delineated intergenic regions from four Encephalitozoon genomes, we predicted de novo new small CDSs with sizes ranging from 78 to 255 bp (median 168) and corroborated these predictions by RACE-PCR experiments in Encephalitozoon cuniculi. Most of the newly found genes are present in other distantly related microsporidian species, suggesting their biological relevance. The present study provides a better framework for annotating microsporidian genomes and to train and evaluate new computational methods dedicated at detecting ultra-small genes in various organisms. [ABSTRACT FROM AUTHOR]

The Cas9 endonuclease is used for genome editing applications in diverse eukaryotic species. A high frequency of off-target activity has been reported in many cell types, limiting its applications to genome engineering, especially in genomic medicine. Here, we generated a synthetic chimeric protein between the catalytic domain of the FokI endonuclease and the catalytically inactive Cas9 protein (fdCas9). A pair of guide RNAs (gRNAs) that bind to sense and antisense strands with a defined spacer sequence range can be used to form a catalytically active dimeric fdCas9 protein and generate double-strand breaks (DSBs) within the spacer sequence. Our data demonstrate an improved catalytic activity of the fdCas9 endonuclease, with a spacer range of 15–39 nucleotides, on surrogate reporters and genomic targets. Furthermore, we observed no detectable fdCas9 activity at known Cas9 off-target sites. Taken together, our data suggest that the fdCas9 endonuclease variant is a superior platform for genome editing applications in eukaryotic systems including mammalian cells. [ABSTRACT FROM AUTHOR]

Modified nucleotide 5-methylcytosine (m5C) is frequently present in various eukaryotic RNAs, including tRNAs, rRNAs and in other non-coding RNAs, as well as in mRNAs. RNA:m5C-methyltranferases (MTases) Nop2 from S. cerevisiae and human proliferation-associated nucleolar antigen p120 are both members of a protein family called Nop2/NSUN/NOL1. Protein p120 is well-known as a tumor marker which is over-expressed in various cancer tissues. Using a combination of RNA bisulfite sequencing and HPLC-MS/MS analysis, we demonstrated here that p120 displays an RNA:m5C- MTase activity, which restores m5C formation at position 2870 in domain V of 25S rRNA in a nop2Δ yeast strain. We also confirm that yeast proteins Nop2p and Rcm1p catalyze the formation of m5C in domains V and IV, respectively. In addition, we do not find any evidence of m5C residues in yeast 18S rRNA. We also performed functional complementation of Nop2-deficient yeasts by human p120 and studied the importance of different sequence and structural domains of Nop2 and p120 for yeast growth and m5C-MTase activity. Chimeric protein formed by Nop2 and p120 fragments revealed the importance of Nop2 N-terminal domain for correct protein localization and its cellular function. We also validated that the presence of Nop2, rather than the m5C modification in rRNA itself, is required for pre-rRNA processing. Our results corroborate that Nop2 belongs to the large family of pre-ribosomal proteins and possesses two related functions in pre-rRNA processing: as an essential factor for cleavages and m5C:RNA:modification. These results support the notion of quality control during ribosome synthesis by such modification enzymes. [ABSTRACT FROM AUTHOR]

Multifunctional protein Dps plays an important role in iron assimilation and a crucial role in bacterial genome packaging. Its monomers form dodecameric spherical particles accumulating ~400 molecules of oxidized iron ions within the protein cavity and applying a flexible N-terminal ends of each subunit for interaction with DNA. Deposition of iron is a well-studied process by which cells remove toxic Fe2+ ions from the genetic material and store them in an easily accessible form. However, the mode of interaction with linear DNA remained mysterious and binary complexes with Dps have not been characterized so far. It is widely believed that Dps binds DNA without any sequence or structural preferences but several lines of evidence have demonstrated its ability to differentiate gene expression, which assumes certain specificity. Here we show that Dps has a different affinity for the two DNA fragments taken from the dps gene regulatory region. We found by atomic force microscopy that Dps predominantly occupies thermodynamically unstable ends of linear double-stranded DNA fragments and has high affinity to the central part of the branched DNA molecule self-assembled from three single-stranded oligonucleotides. It was proposed that Dps prefers binding to those regions in DNA that provide more contact pads for the triad of its DNA-binding bundle associated with one vertex of the protein globule. To our knowledge, this is the first study revealed the nucleoid protein with an affinity to branched DNA typical for genomic regions with direct and inverted repeats. As a ubiquitous feature of bacterial and eukaryotic genomes, such structural elements should be of particular care, but the protein system evolutionarily adapted for this function is not yet known, and we suggest Dps as a putative component of this system. [ABSTRACT FROM AUTHOR]

Eukaryotic ribosome biogenesis is a multistep process involving more than 150 biogenesis factors, which interact transiently with pre-ribosomal particles to promote their maturation. Some of these auxiliary proteins have been isolated in complexes found separate from the ribosomal environment. Among them, are 3 large UTP subcomplexes containing 6 or 7 protein subunits which are involved in the early steps of ribosome biogenesis. The composition of the UTP subcomplexes and the network of binary interactions between protein subunits have been analyzed previously. To obtain further insights into the structural and biochemical properties of UTP subcomplexes, we established a heterologous expression system to allow reconstitution of the yeast tUTP/UTP A and UTP B subcomplexes from their candidate subunits. The results of a series of reconstitution experiments involving different combinations of protein subunits are in good agreement with most of the previously observed binary interactions. Moreover, in combination with additional biochemical analyses, several stable building blocks of the UTP subcomplexes were identified. Based on these findings, we present a refined model of the tUTP/UTP A and UTP B architecture. [ABSTRACT FROM AUTHOR]

Inverted repeats are present in abundance in both prokaryotic and eukaryotic genomes and can form DNA secondary structures – hairpins and cruciforms that are involved in many important biological processes. Bioinformatics tools for efficient and accurate detection of inverted repeats are desirable, because existing tools are often less accurate and time consuming, sometimes incapable of dealing with genome-scale input data. Here, we present a MATLAB-based program called detectIR for the perfect and imperfect inverted repeat detection that utilizes complex numbers and vector calculation and allows genome-scale data inputs. A novel algorithm is adopted in detectIR to convert the conventional sequence string comparison in inverted repeat detection into vector calculation of complex numbers, allowing non-complementary pairs (mismatches) in the pairing stem and a non-palindromic spacer (loop or gaps) in the middle of inverted repeats. Compared with existing popular tools, our program performs with significantly higher accuracy and efficiency. Using genome sequence data from HIV-1, Arabidopsis thaliana, Homo sapiens and Zea mays for comparison, detectIR can find lots of inverted repeats missed by existing tools whose outputs often contain many invalid cases. detectIR is open source and its source code is freely available at: . [ABSTRACT FROM AUTHOR]

Constitutive heterochromatin is a ubiquitous and still unveiled component of eukaryotic genomes, within which it comprises large portions. Although constitutive heterochromatin is generally considered to be transcriptionally silent, it contains a significant variety of sequences that are expressed, among which about 300 single-copy coding genes have been identified by genetic and genomic analyses in the last decades. Here, we report the results of the evolutionary analysis of Yeti, an essential gene of Drosophila melanogaster located in the deep pericentromeric region of chromosome 2R. By FISH, we showed that Yeti maintains a heterochromatin location in both D. simulans and D. sechellia species, closely related to D. melanogaster, while in the more distant species e.g., D. pseudoobscura and D. virilis, it is found within euchromatin, in the syntenic chromosome Muller C, that corresponds to the 2R arm of D. melanogaster chromosome 2. Thus, over evolutionary time, Yeti has been resident on the same chromosomal element, but it progressively moved closer to the pericentric regions. Moreover, in silico reconstruction of the Yeti gene structure in 19 Drosophila species and in 5 non-drosophilid dipterans shows a rather stable organization during evolution. Accordingly, by PCR analysis and sequencing, we found that the single intron of Yeti does not undergo major intraspecies or interspecies size changes, unlike the introns of other essential Drosophila heterochromatin genes, such as light and Dbp80. This implicates diverse evolutionary forces in shaping the structural organization of genes found within heterochromatin. Finally, the results of dS - dN tests show that Yeti is under negative selection both in heterochromatin and euchromatin, and indicate that the change in genomic location did not affected significantly the molecular evolution of the gene. Together, the results of this work contribute to our understanding of the evolutionary dynamics of constitutive heterochromatin in the genomes of higher eukaryotes. [ABSTRACT FROM AUTHOR]

High-quality mapping of genomic regions and genes between two organisms is an indispensable prerequisite for evolutionary analyses and comparative genomics. Existing approaches to this problem focus on either delineating orthologs or finding extended sequence regions of common evolutionary origin (syntenic blocks). We propose SyntenyMapper, a novel tool for refining predefined syntenic regions. SyntenyMapper creates a set of blocks with conserved gene order between two genomes and finds all minor rearrangements that occurred since the evolutionary split of the two species considered. We also present TrackMapper, a SyntenyMapper-based tool that allows users to directly compare genome features, such as histone modifications, between two organisms, and identify genes with highly conserved features. We demonstrate SyntenyMapper's advantages by conducting a large-scale analysis of micro-rearrangements within syntenic regions of 25 eukaryotic species. Unsurprisingly, the number and length of syntenic regions is correlated with evolutionary distance, while the number of micro-rearrangements depends only on the size of the harboring region. On the other hand, the size of rearranged regions remains relatively constant regardless of the evolutionary distance between the organisms, implying a length constraint in the rearrangement process. SyntenyMapper is a useful software tool for both large-scale and gene-centric genome comparisons. [ABSTRACT FROM AUTHOR]

Nucleosomes and their positions in the eukaryotic genome play an important role in regulating gene expression by influencing accessibility to DNA. Many factors influence a nucleosome's final position in the chromatin landscape including the underlying genomic sequence. One of the primary reasons for performing in vitro nucleosome reconstitution experiments is to identify how the underlying DNA sequence will influence a nucleosome's position in the absence of other compounding cellular factors. However, concerns have been raised about the reproducibility of data generated from these kinds of experiments. Here we present data for in vitro nucleosome reconstitution experiments performed on linear plasmid DNA that demonstrate that, when coverage is deep enough, these reconstitution experiments are exquisitely reproducible and highly consistent. Our data also suggests that a coverage depth of 35X be maintained for maximal confidence when assaying nucleosome positions, but lower coverage levels may be generally sufficient. These coverage depth recommendations are sufficient in the experimental system and conditions used in this study, but may vary depending on the exact parameters used in other systems. [ABSTRACT FROM AUTHOR]

The 3′ UTRs of eukaryotic genes participate in a variety of post-transcriptional (and some transcriptional) regulatory interactions. Some of these interactions are well characterised, but an undetermined number remain to be discovered. While some regulatory sequences in 3′ UTRs may be conserved over long evolutionary time scales, others may have only ephemeral functional significance as regulatory profiles respond to changing selective pressures. Here we propose a sensitive segmentation methodology for investigating patterns of composition and conservation in 3′ UTRs based on comparison of closely related species. We describe encodings of pairwise and three-way alignments integrating information about conservation, GC content and transition/transversion ratios and apply the method to three closely related Drosophila species: D. melanogaster, D. simulans and D. yakuba. Incorporating multiple data types greatly increased the number of segment classes identified compared to similar methods based on conservation or GC content alone. We propose that the number of segments and number of types of segment identified by the method can be used as proxies for functional complexity. Our main finding is that the number of segments and segment classes identified in 3′ UTRs is greater than in the same length of protein-coding sequence, suggesting greater functional complexity in 3′ UTRs. There is thus a need for sustained and extensive efforts by bioinformaticians to delineate functional elements in this important genomic fraction. C code, data and results are available upon request. [ABSTRACT FROM AUTHOR]

Mitochondria are eukaryotic organelles supporting individual life-style via generation of proton motive force and cellular energy, and indispensable metabolic pathways. As part of genome sequencing of the white rot Basidiomycota species Phlebia radiata, we first assembled its mitochondrial genome (mtDNA). So far, the 156 348 bp mtDNA is the second largest described for fungi, and of considerable size among eukaryotes. The P. radiata mtDNA assembled as single circular dsDNA molecule containing genes for the large and small ribosomal RNAs, 28 transfer RNAs, and over 100 open reading frames encoding the 14 fungal conserved protein subunits of the mitochondrial complexes I, III, IV, and V. Two genes (atp6 and tRNA-IleGAU) were duplicated within 6.1 kbp inverted region, which is a unique feature of the genome. The large mtDNA size, however, is explained by the dominance of intronic and intergenic regions (sum 80% of mtDNA sequence). The intergenic DNA stretches harness short (≤200 nt) repetitive, dispersed and overlapping sequence elements in abundance. Long self-splicing introns of types I and II interrupt eleven of the conserved genes (cox1,2,3; cob; nad1,2,4,4L,5; rnl; rns). The introns embrace a total of 57 homing endonucleases with LAGLIDADGD and GYI-YIG core motifs, which makes P. radiata mtDNA to one of the largest known reservoirs of intron-homing endonucleases. The inverted duplication, intergenic stretches, and intronic features are indications of dynamics and genetic flexibility of the mtDNA, not fully recognized to this extent in fungal mitochondrial genomes previously, thus giving new insights for the evolution of organelle genomes in eukaryotes. [ABSTRACT FROM AUTHOR]

Two oligomeric types of glycyl-tRNA synthetase (GlyRS) are found in nature: a α2 type and a α2β2 type. The former has been identified in all three kingdoms of life and often pairs with tRNAGly that carries an A73 discriminator base, while the latter is found only in bacteria and chloroplasts and is almost always coupled with tRNAGly that contains U73. In the yeast Saccharomyces cerevisiae, a single GlyRS gene, GRS1, provides both the cytoplasmic and mitochondrial functions, and tRNAGly isoacceptors in both compartments possess A73. We showed herein that Homo sapiens and Arabidopsis thaliana cytoplasmic GlyRSs (both α2-type enzymes) can rescue both the cytoplasmic and mitochondrial defects of a yeast grs1- strain, while Escherichia coli GlyRS (a α2β2-type enzyme) and A. thaliana organellar GlyRS (a (αβ)2-type enzyme) failed to rescue either defect of the yeast mull allele. However, a head-to-tail αβ fusion of E. coli GlyRS effectively supported the mitochondrial function. Our study suggests that a α2-type eukaryotic GlyRS may be functionally substituted with a α2β2-type bacterial cognate enzyme despite their remote evolutionary relationships. [ABSTRACT FROM AUTHOR]

Intermediate-size noncoding RNAs (is-ncRNAs) have been shown to play important regulatory roles in the development of several eukaryotic organisms. However, they have not been thoroughly explored in Plasmodium falciparum, which is the most virulent malaria parasite infecting human being. By using Illumina/Solexa paired-end sequencing of an is-ncRNA-specific library, we performed a systematic identification of novel is-ncRNAs in intraerythrocytic P. falciparum, strain 3D7. A total of 1,198 novel is-ncRNA candidates, including antisense, intergenic, and intronic is-ncRNAs, were identified. Bioinformatics analyses showed that the intergenic is-ncRNAs were the least conserved among different Plasmodium species, and antisense is-ncRNAs were more conserved than their sense counterparts. Twenty-two novel snoRNAs were identified, and eight potential novel classes of P. falciparum is-ncRNAs were revealed by clustering analysis. The expression of randomly selected novel is-ncRNAs was confirmed by RT-PCR and northern blotting assays. An obvious different expressional profile of the novel is-ncRNA between the early and late intraerythrocytic developmental stages of the parasite was observed. The expression levels of the antisense RNAs correlated with those of their cis-encoded sense RNA counterparts, suggesting that these is-ncRNAs are involved in the regulation of gene expression of the parasite. In conclusion, we accomplished a deep profiling analysis of novel is-ncRNAs in P. falciparum, analysed the conservation and structural features of these novel is-ncRNAs, and revealed their differential expression patterns during the development of the parasite. These findings provide important information for further functional characterisation of novel is-ncRNAs during the development of P. falciparum. [ABSTRACT FROM AUTHOR]

Pro- and eukaryotic microbes associated with multi-cellular organisms are receiving increasing attention as a driving factor in ecosystems. Endophytes in plants can change host performance by altering nutrient uptake, secondary metabolite production or defense mechanisms. Recent studies detected widespread prevalence of Labyrinthula zosterae in European Zostera marina meadows, a protist that allegedly caused a massive amphi-Atlantic seagrass die-off event in the 1930's, while showing only limited virulence today. As a limiting factor for pathogenicity, we investigated genotype×genotype interactions of host and pathogen from different regions (10–100 km-scale) through reciprocal infection. Although the endophyte rapidly infected Z. marina, we found little evidence that Z. marina was negatively impacted by L. zosterae. Instead Z. marina showed enhanced leaf growth and kept endophyte abundance low. Moreover, we found almost no interaction of protist×eelgrass-origin on different parameters of L. zosterae virulence/Z. marina performance, and also no increase in mortality after experimental infection. In a target gene approach, we identified a significant down-regulation in the expression of 6/11 genes from the defense cascade of Z. marina after real-time quantitative PCR, revealing strong immune modulation of the host's defense by a potential parasite for the first time in a marine plant. Nevertheless, one gene involved in phenol synthesis was strongly up-regulated, indicating that Z. marina plants were probably able to control the level of infection. There was no change in expression in a general stress indicator gene (HSP70). Mean L. zosterae abundances decreased below 10% after 16 days of experimental runtime. We conclude that under non-stress conditions L. zosterae infection in the study region is not associated with substantial virulence. [ABSTRACT FROM AUTHOR]

Intermediate-size noncoding RNAs (is-ncRNAs) have been shown to play important regulatory roles in the development of several eukaryotic organisms. However, they have not been thoroughly explored in Plasmodium falciparum, which is the most virulent malaria parasite infecting human being. By using Illumina/Solexa paired-end sequencing of an is-ncRNA-specific library, we performed a systematic identification of novel is-ncRNAs in intraerythrocytic P. falciparum, strain 3D7. A total of 1,198 novel is-ncRNA candidates, including antisense, intergenic, and intronic is-ncRNAs, were identified. Bioinformatics analyses showed that the intergenic is-ncRNAs were the least conserved among different Plasmodium species, and antisense is-ncRNAs were more conserved than their sense counterparts. Twenty-two novel snoRNAs were identified, and eight potential novel classes of P. falciparum is-ncRNAs were revealed by clustering analysis. The expression of randomly selected novel is-ncRNAs was confirmed by RT-PCR and northern blotting assays. An obvious different expressional profile of the novel is-ncRNA between the early and late intraerythrocytic developmental stages of the parasite was observed. The expression levels of the antisense RNAs correlated with those of their cis-encoded sense RNA counterparts, suggesting that these is-ncRNAs are involved in the regulation of gene expression of the parasite. In conclusion, we accomplished a deep profiling analysis of novel is-ncRNAs in P. falciparum, analysed the conservation and structural features of these novel is-ncRNAs, and revealed their differential expression patterns during the development of the parasite. These findings provide important information for further functional characterisation of novel is-ncRNAs during the development of P. falciparum. [ABSTRACT FROM AUTHOR]

Reconstructing synteny blocks is an essential step in comparative genomics studies. Different methods were already developed to answer various needs such as genome (re-)annotation, identification of duplicated regions and whole genome duplication events or estimation of rearrangement rates. We present SynChro, a tool that reconstructs synteny blocks between pairwise comparisons of multiple genomes. SynChro is based on a simple algorithm that computes Reciprocal Best-Hits (RBH) to reconstruct the backbones of the synteny blocks and then automatically completes these blocks with non-RBH syntenic homologs. This approach has two main advantages: (i) synteny block reconstruction is fast (feasible on a desk computer for large eukaryotic genomes such as human) and (ii) synteny block reconstruction is straightforward as all steps are integrated (no need to run Blast or TribeMCL prior to reconstruction) and there is only one parameter to set up, the synteny block stringency . Benchmarks on three pairwise comparisons of genomes, representing three different levels of synteny conservation (Human/Mouse, Human/Zebra Finch and Human/Zebrafish) show that Synchro runs faster and performs at least as well as two other commonly used and more sophisticated tools (MCScanX and i-ADHoRe). In addition, SynChro provides the user with a rich set of graphical outputs including dotplots, chromosome paintings and detailed synteny maps to visualize synteny blocks with all homology relationships and synteny breakpoints with all included genetic features. SynChro is freely available under the BSD license at http://www.lcqb.upmc.fr/CHROnicle/SynChro.html. [ABSTRACT FROM AUTHOR]

Many prokaryotic and eukaryotic genomes feature a characteristic periodic signal in distribution of short runs of A or T (A-tracts) phased with the DNA helical period of ∼10–11 bp. Such periodic spacing of A-tracts has been associated with intrinsic DNA curvature. In eukaryotes, this periodicity is a major component of the nucleosome positioning signal but its physiological role in prokaryotes is not clear. One hypothesis centers on possible role of intrinsic DNA bends in nucleoid compaction. We use comparative genomics to investigate possible relationship between the A-tract periodicity and nucleoid-associated proteins in prokaryotes. We found that genomes with DNA-bridging proteins tend to exhibit stronger A-tract periodicity, presumably indicative of more prevalent intrinsic DNA curvature. A weaker relationship was detected for nucleoid-associated proteins that do not form DNA bridges. We consider these results an indication that intrinsic DNA curvature acts collaboratively with DNA-bridging proteins in maintaining the compact structure of the nucleoid, and that previously observed differences among prokaryotic genomes in terms DNA curvature-related sequence periodicity may reflect differences in nucleoid organization. We subsequently investigated the relationship between A-tract periodicity and presence of CRISPR elements and we found that genomes with CRISPR tend to have stronger A-tract periodicity. This result is consistent with our earlier hypothesis that extensive A-tract periodicity could help protect the chromosome against integration of prophages, possibly due to its role in compaction of the nucleoid. [ABSTRACT FROM AUTHOR]

Ribosome heterogeneity is of increasing biological significance and several examples have been described for multicellular and single cells organisms. In here we show for the first time a variation in ribose methylation within the 18S rRNA of Saccharomyces cerevisiae. Using RNA-cleaving DNAzymes, we could specifically demonstrate that a significant amount of S. cerevisiae ribosomes are not methylated at 2′-O-ribose of A100 residue in the 18S rRNA. Furthermore, using LC-UV-MS/MS of a respective 18S rRNA fragment, we could not only corroborate the partial methylation at A100, but could also quantify the methylated versus non-methylated A100 residue. Here, we exhibit that only 68% of A100 in the 18S rRNA of S.cerevisiae are methylated at 2′-O ribose sugar. Polysomes also contain a similar heterogeneity for methylated Am100, which shows that 40S ribosome subunits with and without Am100 participate in translation. Introduction of a multicopy plasmid containing the corresponding methylation guide snoRNA gene SNR51 led to an increased A100 methylation, suggesting the cellular snR51 level to limit the extent of this modification. Partial rRNA modification demonstrates a new level of ribosome heterogeneity in eukaryotic cells that might have substantial impact on regulation and fine-tuning of the translation process. [ABSTRACT FROM AUTHOR]

Replication protein A (RPA) is a ubiquitous eukaryotic single-stranded DNA (ssDNA) binding protein necessary for all aspects of DNA metabolism involving an ssDNA intermediate, including DNA replication, repair, recombination, DNA damage response and checkpoint activation, and telomere maintenance [1], [2], [3]. The role of RPA in most of these reactions is to protect the ssDNA until it can be delivered to downstream enzymes. Therefore a crucial feature of RPA is that it must bind very tightly to ssDNA, but must also be easily displaced from ssDNA to allow other proteins to gain access to the substrate. Here we use total internal reflection fluorescence microscopy and nanofabricated DNA curtains to visualize the behavior of Saccharomyces cerevisiae RPA on individual strands of ssDNA in real-time. Our results show that RPA remains bound to ssDNA for long periods of time when free protein is absent from solution. In contrast, RPA rapidly dissociates from ssDNA when free RPA is present in solution allowing rapid exchange between the free and bound states. In addition, the S. cerevisiae DNA recombinase Rad51 and E. coli single-stranded binding protein (SSB) also promote removal of RPA from ssDNA. These results reveal an unanticipated exchange between bound and free RPA suggesting a binding mechanism that can confer exceptionally slow off rates, yet also enables rapid displacement through a direct exchange mechanism that is reliant upon the presence of free ssDNA-binding proteins in solution. Our results indicate that RPA undergoes constant microscopic dissociation under all conditions, but this is only manifested as macroscopic dissociation (i.e. exchange) when free proteins are present in solution, and this effect is due to mass action. We propose that the dissociation of RPA from ssDNA involves a partially dissociated intermediate, which exposes a small section of ssDNA allowing other proteins to access to the DNA. [ABSTRACT FROM AUTHOR]

High throughput sequencing technology has great promise for biodiversity studies. However, an underlying assumption is that the primers used in these studies are universal for the prokaryotic or eukaryotic groups of interest. Full primer universality is difficult or impossible to achieve and studies using different primer sets make biodiversity comparisons problematic. The aim of this study was to design and optimize universal eukaryotic primers that could be used as a standard in future biodiversity studies. Using the alignment of all eukaryotic sequences from the publicly available SILVA database, we generated a full characterization of variable versus conserved regions in the 18S rRNA gene. All variable regions within this gene were analyzed and our results suggested that the V2, V4 and V9 regions were best suited for biodiversity assessments. Previously published universal eukaryotic primers as well as a number of self-designed primers were mapped to the alignment. Primer selection will depend on sequencing technology used, and this study focused on the 454 pyrosequencing GS FLX Titanium platform. The results generated a primer pair yielding theoretical matches to 80% of the eukaryotic and 0% of the prokaryotic sequences in the SILVA database. An empirical test of marine sediments using the AmpliconNoise pipeline for analysis of the high throughput sequencing data yielded amplification of sequences for 71% of all eukaryotic phyla with no isolation of prokaryotic sequences. To our knowledge this is the first characterization of the complete 18S rRNA gene using all eukaryotes present in the SILVA database, providing a robust test for universal eukaryotic primers. Since both in silico and empirical tests using high throughput sequencing retained high inclusion of eukaryotic phyla and exclusion of prokaryotes, we conclude that these primers are well suited for assessing eukaryote diversity, and can be used as a standard in biodiversity studies. [ABSTRACT FROM AUTHOR]

Whole-Genome Bisulfite Sequencing (WGBS) and genome-wide Reduced Representation Bisulfite Sequencing (RRBS) are widely used to study DNA methylation. However, data analysis is complicated, lengthy, and hampered by a lack of seamless analytical pipelines. To address these issues, we developed a convenient, stable, and efficient web service called Web Service for Bisulfite Sequencing Data Analysis (WBSA) to analyze bisulfate sequencing data. WBSA focuses on not only CpG methylation, which is the most common biochemical modification in eukaryotic DNA, but also non-CG methylation, which have been observed in plants, iPS cells, oocytes, neurons and stem cells of human. WBSA comprises three main modules as follows: WGBS data analysis, RRBS data analysis, and differentially methylated region (DMR) identification. The WGBS and RRBS modules execute read mapping, methylation site identification, annotation, and advanced analysis, whereas the DMR module identifies actual DMRs and annotates their correlations to genes. WBSA can be accessed and used without charge either online or local version. WBSA also includes the executables of the Portable Batch System (PBS) and standalone versions that can be downloaded from the website together with the installation instructions. WBSA is available at no charge for academic users at http://wbsa.big.ac.cn. [ABSTRACT FROM AUTHOR]

Chinese hamster ovary (CHO) cell lines represent the most commonly used mammalian expression system for the production of therapeutic proteins. In this context, detailed knowledge of the CHO cell transcriptome might help to improve biotechnological processes conducted by specific cell lines. Nevertheless, very few assembled cDNA sequences of CHO cells were publicly released until recently, which puts a severe limitation on biotechnological research. Two extended annotation systems and web-based tools, one for browsing eukaryotic genomes (GenDBE) and one for viewing eukaryotic transcriptomes (SAMS), were established as the first step towards a publicly usable CHO cell genome/transcriptome analysis platform. This is complemented by the development of a new strategy to assemble the ca. 100 million reads, sequenced from a broad range of diverse transcripts, to a high quality CHO cell transcript set. The cDNA libraries were constructed from different CHO cell lines grown under various culture conditions and sequenced using Roche/454 and Illumina sequencing technologies in addition to sequencing reads from a previous study. Two pipelines to extend and improve the CHO cell line transcripts were established. First, de novo assemblies were carried out with the Trinity and Oases assemblers, using varying k-mer sizes. The resulting contigs were screened for potential CDS using ESTScan. Redundant contigs were filtered out using cd-hit-est. The remaining CDS contigs were re-assembled with CAP3. Second, a reference-based assembly with the TopHat/Cufflinks pipeline was performed, using the recently published draft genome sequence of CHO-K1 as reference. Additionally, the de novo contigs were mapped to the reference genome using GMAP and merged with the Cufflinks assembly using the cuffmerge software. With this approach 28,874 transcripts located on 16,492 gene loci could be assembled. Combining the results of both approaches, 65,561 transcripts were identified for CHO cell lines, which could be clustered by sequence identity into 17,598 gene clusters. [ABSTRACT FROM AUTHOR]

Telomeres are repetitive non coding DNA sequences located at the end of eukaryotic chromosomes, which maintain the integrity of the genome by hiding the chromosome ends from being recognised as double stranded breaks. Telomeres are emerging as biomarkers for ageing and survival, and are susceptible to reflect different individual life history trajectories. In particular, the telomere length with which one starts in life has been shown to be linked with individual life-long survival, suggesting that telomere dynamics can be a proxy for individual fitness and thereby be implicated in evolutionary trade-offs. As a consequence, an increasing number of studies were conducted on telomeres in the fields of ecology and evolutionary biology, in which telomere length was almost exclusively measured from blood samples. However, not only do the number of repeats of the telomeric sequences vary among species, but also within species with great inter-individual telomere lengths variability with age, tissues, and chromosomes. This raises the issue of the exact biological meaning of telomere measurement in blood cells and stimulated the study of the correlation of telomere lengths among tissues over age. By measuring telomere length in adult zebra finches (Taeniopygia guttata) in different somatic tissues displaying variable cell turnovers (bone marrow, brain, spleen, pectoral muscle, heart, liver and in red blood cells), we checked that the measure of telomere length in red blood cells is related to telomere lengths in the other tissues. Here we show significant relationships between the telomere lengths of red blood cells and several somatic tissues at adulthood. As red blood cells are easily accessible and suitable for the longitudinal monitoring of the individual rate of telomere loss, our study confirms that telomere length measured in red blood cells could serve as a surrogate for telomere length in the whole avian organism. [ABSTRACT FROM AUTHOR]

Genetic linkage maps are indispensable tools in genetic, genomic and breeding studies. As one of genotyping-by-sequencing methods, RAD-Seq (restriction-site associated DNA sequencing) has gained particular popularity for construction of high-density linkage maps. Current RAD analytical tools are being predominantly used for typing codominant markers. However, no genotyping algorithm has been developed for dominant markers (resulting from recognition site disruption). Given their abundance in eukaryotic genomes, utilization of dominant markers would greatly diminish the extensive sequencing effort required for large-scale marker development. In this study, we established, for the first time, a novel statistical framework for de novo dominant genotyping in mapping populations. An integrated package called RADtyping was developed by incorporating both de novo codominant and dominant genotyping algorithms. We demonstrated the superb performance of RADtyping in achieving remarkably high genotyping accuracy based on simulated and real mapping datasets. The RADtyping package is freely available at http://www2.ouc.edu.cn/mollusk/ detailen.asp?id=727. [ABSTRACT FROM AUTHOR]

Alpha-solenoids are flexible protein structural domains formed by ensembles of alpha-helical repeats (Armadillo and HEAT repeats among others). While homology can be used to detect many of these repeats, some alpha-solenoids have very little sequence homology to proteins of known structure and we expect that many remain undetected. We previously developed a method for detection of alpha-helical repeats based on a neural network trained on a dataset of protein structures. Here we improved the detection algorithm and updated the training dataset using recently solved structures of alpha-solenoids. Unexpectedly, we identified occurrences of alpha-solenoids in solved protein structures that escaped attention, for example within the core of the catalytic subunit of PI3KC. Our results expand the current set of known alpha-solenoids. Application of our tool to the protein universe allowed us to detect their significant enrichment in proteins interacting with many proteins, confirming that alpha-solenoids are generally involved in protein-protein interactions. We then studied the taxonomic distribution of alpha-solenoids to discuss an evolutionary scenario for the emergence of this type of domain, speculating that alpha-solenoids have emerged in multiple taxa in independent events by convergent evolution. We observe a higher rate of alpha-solenoids in eukaryotic genomes and in some prokaryotic families, such as Cyanobacteria and Planctomycetes, which could be associated to increased cellular complexity. The method is available at http://cbdm.mdc-berlin.de/~ard2/. [ABSTRACT FROM AUTHOR]