The edges of understanding

A culture's icons are a window onto its soul. Few would disagree that, in the culture of molecular biology that dominated much of the life sciences for the last third of the 20th century, the dominant icon was the double helix. In the present, post-modern, 'systems biology' era, however, it is, arguably, the hairball

Expedited batch processing and analysis of transposon insertions

<p>Abstract</p> <p>Background</p> <p>With advances in sequencing technology, greater and greater amounts of eukaryotic genome data are becoming available. Often, large portions of these genomes consist of transposable elements, frequently accounting for 50% or more in vertebrates. Each transposable element family may have thousands or tens of thousands of individual copies within a given genome, and therefore it can take an exorbitant amount of time and effort to process data in a meaningful fashion.</p> <p>Findings</p> <p>In order to combat this problem, we developed a set of bioinformatics techniques and programs to streamline the analysis. This includes a unique Perl script which automates the process of taking BLAST, Repeatmasker and similar data to extract and manipulate the hit sequences from the genome. This script, called Process_hits uses an object-oriented methodology to compile all hit locations from a given file for processing, organize this data into useable categories, and output it in multiple formats.</p> <p>Conclusions</p> <p>The program proved capable of handling large amounts of transposon data in an efficient fashion. It is equipped with a number of useful sub-functions, each of which is contained within its own sub-module to allow for greater expandability and as a foundation for future program design.</p

Enrichment analysis of Alu elements with different spatial chromatin proximity in the human genome

Transposable elements (TEs) have no longer been totally considered as “junk DNA” for quite a time since the continual discoveries of their multifunctional roles in eukaryote genomes. As one of the most important and abundant TEs that still active in human genome, Alu, a SINE family, has demonstrated its indispensable regulatory functions at sequence level, but its spatial roles are still unclear. Technologies based on 3C(chromosomeconformation capture) have revealed the mysterious three-dimensional structure of chromatin, and make it possible to study the distal chromatin interaction in the genome. To find the role TE playing in distal regulation in human genome, we compiled the new released Hi-C data, TE annotation, histone marker annotations, and the genome-wide methylation data to operate correlation analysis, and found that the density of Alu elements showed a strong positive correlation with the level of chromatin interactions (hESC: r=0.9, P<2.2×1016; IMR90 fibroblasts: r = 0.94, P < 2.2 × 1016) and also have a significant positive correlation withsomeremote functional DNA elements like enhancers and promoters (Enhancer: hESC: r=0.997, P=2.3×10−4; IMR90: r=0.934, P=2×10−2; Promoter: hESC: r = 0.995, P = 3.8 × 10−4; IMR90: r = 0.996, P = 3.2 × 10−4). Further investigation involving GC content and methylation status showed the GC content of Alu covered sequences shared a similar pattern with that of the overall sequence, suggesting that Alu elements also function as the GC nucleotide and CpG site provider. In all, our results suggest that the Alu elements may act as an alternative parameter to evaluate the Hi-C data, which is confirmed by the correlation analysis of Alu elements and histone markers. Moreover, the GC-rich Alu sequence can bring high GC content and methylation flexibility to the regions with more distal chromatin contact, regulating the transcription of tissue-specific genes

Viral population estimation using pyrosequencing

The diversity of virus populations within single infected hosts presents a major difficulty for the natural immune response as well as for vaccine design and antiviral drug therapy. Recently developed pyrophosphate based sequencing technologies (pyrosequencing) can be used for quantifying this diversity by ultra-deep sequencing of virus samples. We present computational methods for the analysis of such sequence data and apply these techniques to pyrosequencing data obtained from HIV populations within patients harboring drug resistant virus strains. Our main result is the estimation of the population structure of the sample from the pyrosequencing reads. This inference is based on a statistical approach to error correction, followed by a combinatorial algorithm for constructing a minimal set of haplotypes that explain the data. Using this set of explaining haplotypes, we apply a statistical model to infer the frequencies of the haplotypes in the population via an EM algorithm. We demonstrate that pyrosequencing reads allow for effective population reconstruction by extensive simulations and by comparison to 165 sequences obtained directly from clonal sequencing of four independent, diverse HIV populations. Thus, pyrosequencing can be used for cost-effective estimation of the structure of virus populations, promising new insights into viral evolutionary dynamics and disease control strategies.Comment: 23 pages, 13 figure

Identification of QTL genes for BMD variation using both linkage and gene-based association approaches

Low bone mineral density (BMD) is a risk factor for osteoporotic fracture with a high heritability. Previous large scale linkage study in Northern Chinese has identified four significant quantitative trait loci (QTL) for BMD variation on chromosome 2q24, 5q21, 7p21 and 13q21. We performed a replication study of these four QTL in 1,459 Southern Chinese from 306 pedigrees. Successful replication was observed on chromosome 5q21 for femoral neck BMD with a LOD score of 1.38 (nominal p value = 0.006). We have previously identified this locus in a genome scan meta-analysis of BMD variation in a white population. Subsequent QTL-wide gene-based association analysis in 800 subjects with extreme BMD identified CAST and ERAP1 as novel BMD candidate genes (empirical p value of 0.032 and 0.014, respectively). The associations were independently replicated in a Northern European population (empirical p value of 0.01 and 0.004 for CAST and ERAP1, respectively). These findings provide further evidence that 5q21 is a BMD QTL, and CAST and ERAP1 may be associated with femoral neck BMD variation

Widespread horizontal transfer of mitochondrial genes in flowering plants

Horizontal gene transfer - the exchange of genes across mating barriers - is recognized as a major force in bacterial evolution(1,2). However, in eukaryotes it is prevalent only in certain phagotrophic protists and limited largely to the ancient acquisition of bacterial genes(3-5). Although the human genome was initially reported(6) to contain over 100 genes acquired during vertebrate evolution from bacteria, this claim was immediately and repeatedly rebutted(7,8). Moreover, horizontal transfer is unknown within the evolution of animals, plants and fungi except in the special context of mobile genetic elements(9-12). Here we show, however, that standard mitochondrial genes, encoding ribosomal and respiratory proteins, are subject to evolutionarily frequent horizontal transfer between distantly related flowering plants. These transfers have created a variety of genomic outcomes, including gene duplication, recapture of genes lost through transfer to the nucleus, and chimaeric, half-monocot, half-dicot genes. These results imply the existence of mechanisms for the delivery of DNA between unrelated plants, indicate that horizontal transfer is also a force in plant nuclear genomes, and are discussed in the contexts of plant molecular phylogeny and genetically modified plants.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/62688/1/nature01743.pd

Three-dimensional structure of a viral genome-delivery portal vertex.

DNA viruses such as bacteriophages and herpesviruses deliver their genome into and out of the capsid through large proteinaceous assemblies, known as portal proteins. Here, we report two snapshots of the dodecameric portal protein of bacteriophage P22. The 3.25-Å-resolution structure of the portal-protein core bound to 12 copies of gene product 4 (gp4) reveals a ~1.1-MDa assembly formed by 24 proteins. Unexpectedly, a lower-resolution structure of the full-length portal protein unveils the unique topology of the C-terminal domain, which forms a ~200-Å-long α-helical barrel. This domain inserts deeply into the virion and is highly conserved in the Podoviridae family. We propose that the barrel domain facilitates genome spooling onto the interior surface of the capsid during genome packaging and, in analogy to a rifle barrel, increases the accuracy of genome ejection into the host cell