Life in the fast lane for protein crystallization and X-ray crystallography

The common goal for structural genomic centers and consortiums is to decipher as quickly as possible the three-dimensional structures for a multitude of recombinant proteins derived from known genomic sequences. Since X-ray crystallography is the foremost method to acquire atomic resolution for macromolecules, the limiting step is obtaining protein crystals that can be useful of structure determination. High-throughput methods have been developed in recent years to clone, express, purify, crystallize and determine the three-dimensional structure of a protein gene product rapidly using automated devices, commercialized kits and consolidated protocols. However, the average number of protein structures obtained for most structural genomic groups has been very low compared to the total number of proteins purified. As more entire genomic sequences are obtained for different organisms from the three kingdoms of life, only the proteins that can be crystallized and whose structures can be obtained easily are studied. Consequently, an astonishing number of genomic proteins remain unexamined. In the era of high-throughput processes, traditional methods in molecular biology, protein chemistry and crystallization are eclipsed by automation and pipeline practices. The necessity for high-rate production of protein crystals and structures has prevented the usage of more intellectual strategies and creative approaches in experimental executions. Fundamental principles and personal experiences in protein chemistry and crystallization are minimally exploited only to obtain “low-hanging fruit” protein structures. We review the practical aspects of today's high-throughput manipulations and discuss the challenges in fast pace protein crystallization and tools for crystallography. Structural genomic pipelines can be improved with information gained from low-throughput tactics that may help us reach the higher-bearing fruits. Examples of recent developments in this area are reported from the efforts of the Southeast Collaboratory for Structural Genomics (SECSG).     

Mouse Models of Genomic Syndromes as Tools for Understanding the Basis of Complex Traits: An Example with the Smith-Magenis and the Potocki-Lupski Syndromes

Each human''s genome is distinguished by extra and missing DNA that can be “benign” or powerfully impact everything from development to disease. In the case of genomic disorders DNA rearrangements, such as deletions or duplications, correlate with a clinical specific phenotype. The clinical presentations of genomic disorders were thought to result from altered gene copy number of physically linked dosage sensitive genes. Genomic disorders are frequent diseases (~1 per 1,000 births). Smith-Magenis syndrome (SMS) and Potocki-Lupski syndrome (PTLS) are genomic disorders, associated with a deletion and a duplication, of 3.7 Mb respectively, within chromosome 17 band p11.2. This region includes 23 genes. Both syndromes have complex and distinctive phenotypes including multiple congenital and neurobehavioral abnormalities. Human chromosome 17p11.2 is syntenic to the 32-34 cM region of murine chromosome 11. The number and order of the genes are highly conserved. In this review, we will exemplify how genomic disorders can be modeled in mice and the advantages that such models can give in the study of genomic disorders in particular and gene copy number variation (CNV) in general. The contributions of the SMS and PTLS animal models in several aspects ranging from more specific ones, as the definition of the clinical aspects of the human clinical spectrum, the identification of dosage sensitive genes related to the human syndromes, to the more general contributions as the definition of genetic locus impacting obesity and behavior and the elucidation of general mechanisms related to the pathogenesis of gene CNV are discussed.Key Words: Gene copy number variation, complex traits, phenotypic consequences, mouse models.     

Efficient sorting of genomic permutations by translocation, inversion and block interchange

MOTIVATION: Finding genomic distance based on gene order is a classic problem in genome rearrangements. Efficient exact algorithms for genomic distances based on inversions and/or translocations have been found but are complicated by special cases, rare in simulations and empirical data. We seek a universal operation underlying a more inclusive set of evolutionary operations and yielding a tractable genomic distance with simple mathematical form. RESULTS: We study a universal double-cut-and-join operation that accounts for inversions, translocations, fissions and fusions, but also produces circular intermediates which can be reabsorbed. The genomic distance, computable in linear time, is given by the number of breakpoints minus the number of cycles (b-c) in the comparison graph of the two genomes; the number of hurdles does not enter into it. Without changing the formula, we can replace generation and re-absorption of a circular intermediate by a generalized transposition, equivalent to a block interchange, with weight two. Our simple algorithm converts one multi-linear chromosome genome to another in the minimum distance.     

Detection of plant genes using a rapid, nonorganic DNA purification method

We have developed a simple procedure for the preparation of plant genomic DNA using FTA paper. Plant leaves were crushed against FTA paper, and the genomic DNA was purified using simple, nonorganic reagents. The 18S rRNA gene and the gene encoding the ribulose-1, 5-bisphosphate carboxylase/oxygenase large subunit (rbcL) from the chloroplast genome were detected by PCR amplification of DNA on FTA paper. DNA amplification was successful using extracts from 16 dicot and monocot plants. Studies of specific plant extracts revealed that extracts of leaf samples could be collected and stored at room temperature on FTA paper without a decrease in the DNA amplification success rate for more than a month. Both the 18S RNA gene and the rbcL gene were detected in the genomic DNA isolated from various soybean cultivars stored in this manner. Furthermore, by modestly increasing the number of cycles of DNA amplification, we were able to detect the uidA gene in transgenic tobacco and rice leaves as well as a single copy gene linked to the resistance gene of cyst nematode race 3 using genomic DNA isolated on FTA paper. These results demonstrate that genomic DNA isolated using FTA paper can be used for the detection of plant genes, from a wide range of plants with either high or low gene copy number and of either nuclear or cytoplasmic origin.     

Incorporating 16S Gene Copy Number Information Improves Estimates of Microbial Diversity and Abundance

The abundance of different SSU rRNA (“16S”) gene sequences in environmental samples is widely used in studies of microbial ecology as a measure of microbial community structure and diversity. However, the genomic copy number of the 16S gene varies greatly – from one in many species to up to 15 in some bacteria and to hundreds in some microbial eukaryotes. As a result of this variation the relative abundance of 16S genes in environmental samples can be attributed both to variation in the relative abundance of different organisms, and to variation in genomic 16S copy number among those organisms. Despite this fact, many studies assume that the abundance of 16S gene sequences is a surrogate measure of the relative abundance of the organisms containing those sequences. Here we present a method that uses data on sequences and genomic copy number of 16S genes along with phylogenetic placement and ancestral state estimation to estimate organismal abundances from environmental DNA sequence data. We use theory and simulations to demonstrate that 16S genomic copy number can be accurately estimated from the short reads typically obtained from high-throughput environmental sequencing of the 16S gene, and that organismal abundances in microbial communities are more strongly correlated with estimated abundances obtained from our method than with gene abundances. We re-analyze several published empirical data sets and demonstrate that the use of gene abundance versus estimated organismal abundance can lead to different inferences about community diversity and structure and the identity of the dominant taxa in microbial communities. Our approach will allow microbial ecologists to make more accurate inferences about microbial diversity and abundance based on 16S sequence data.     

Comparative genomic hybridization (CGH): ten years of substantial progress in human solid tumor molecular cytogenetics

Data from ten years of research using comparative genomic hybridization (CGH) for the detection of chromosomal alterations in human solid tumors are concisely reviewed. By use of a basic methodology with some variations more or less specific patterns of genomic imbalances were found in a large number of tumors of various entities. Specific gains and losses of genomic material have not only opened the way to the detection of a series of cancer-related genes but also to clinical implications. Not only several areas of basic oncogenetic research, but also differential diagnosis, prognosis of disease progression, and therapeutic decisions have profited by CGH.     

Genomic evaluations with many more genotypes

Background

Genomic evaluations in Holstein dairy cattle have quickly become more reliable over the last two years in many countries as more animals have been genotyped for 50,000 markers. Evaluations can also include animals genotyped with more or fewer markers using new tools such as the 777,000 or 2,900 marker chips recently introduced for cattle. Gains from more markers can be predicted using simulation, whereas strategies to use fewer markers have been compared using subsets of actual genotypes. The overall cost of selection is reduced by genotyping most animals at less than the highest density and imputing their missing genotypes using haplotypes. Algorithms to combine different densities need to be efficient because numbers of genotyped animals and markers may continue to grow quickly.

Methods

Genotypes for 500,000 markers were simulated for the 33,414 Holsteins that had 50,000 marker genotypes in the North American database. Another 86,465 non-genotyped ancestors were included in the pedigree file, and linkage disequilibrium was generated directly in the base population. Mixed density datasets were created by keeping 50,000 (every tenth) of the markers for most animals. Missing genotypes were imputed using a combination of population haplotyping and pedigree haplotyping. Reliabilities of genomic evaluations using linear and nonlinear methods were compared.

Results

Differing marker sets for a large population were combined with just a few hours of computation. About 95% of paternal alleles were determined correctly, and > 95% of missing genotypes were called correctly. Reliability of breeding values was already high (84.4%) with 50,000 simulated markers. The gain in reliability from increasing the number of markers to 500,000 was only 1.6%, but more than half of that gain resulted from genotyping just 1,406 young bulls at higher density. Linear genomic evaluations had reliabilities 1.5% lower than the nonlinear evaluations with 50,000 markers and 1.6% lower with 500,000 markers.

Conclusions

Methods to impute genotypes and compute genomic evaluations were affordable with many more markers. Reliabilities for individual animals can be modified to reflect success of imputation. Breeders can improve reliability at lower cost by combining marker densities to increase both the numbers of markers and animals included in genomic evaluation. Larger gains are expected from increasing the number of animals than the number of markers.     

The breakdown of genomic ancestry blocks in hybrid lineages given a finite number of recombination sites

When a lineage originates from hybridization genomic blocks of contiguous ancestry from different ancestors are fragmented through genetic recombination. The resulting blocks are delineated by so called junctions, which accumulate with every generation that passes. Modeling the accumulation of ancestry block junctions can elucidate processes and timeframes of genomic admixture. Previous models have not addressed ancestry block dynamics for chromosomes that consist of a finite number of recombination sites. However, genomic data typically consist of informative markers that are interspersed with fragments for which no ancestry information is available. Hence, repeated recombination events may occur between markers, effectively removing existing junctions. Here, we present an analytical treatment of the dynamics of the mean number of junctions over time, taking into account the number of recombination sites per chromosome, population size, genetic map length, and the frequency of the ancestral species in the founding hybrid swarm. We describe the expected number of junctions using equidistant molecular markers and estimate the number of junctions using random markers. This extended theory of junctions thus reflects properties of empirical data and can serve to study the genomic patterns following admixture.     

Solidarity and justice as guiding principles in genomic research

In genomic research the ideal standard of free, informed, prior and explicit consent is sometimes difficult to apply. This has raised concern that important genomic research will be restricted. Different consent procedures have therefore been proposed. This paper explicitly examines the question how, in genomic research, the principles of solidarity and justice can be used to justify forms of diminished individual control over personal data and bio-samples. After a discussion of the notions of solidarity and justice and how they can be related to health care and genomic research, we examine how and in which situations these notions can form a strong moral basis for demanding certain financial sacrifices. Then we examine when these principles can justify consent procedures which diverge from the ideal standard. Because much genomic research is not expected to lead to immediate (clinical) benefits we also discuss the question of whether we can be obliged to make any sacrifices for future (not yet existing) patients. We conclude with the formulation of a number of conditions that have to be met before autonomy sacrifices can be reasonably demanded in genomic research.     

The death-inducing effect and genomic instability

Exposure to ionizing radiation can induce a heritable change in the unirradiated progeny of irradiated cells. This non-targeted effect of ionizing radiation manifests as genomic instability, and although there is some debate as to the role of genomic instability in the carcinogenic process, it is thought by some to be an early step in radiation carcinogenesis. Although the mechanism of induction of genomic instability is not clearly understood, evidence suggests that secreted factors from irradiated cells may be involved. We have previously identified another non-targeted effect of ionizing radiation, the death-inducing effect. Exposure of unirradiated GM10115 cells to medium from chromosomally unstable clones was generally found to be cytotoxic. However, occasionally cells will survive in medium from unstable clones and can be clonally expanded. The absolute yield of survivors is independent of the initial number of cells plated when cell densities reached 5,000 or more cells/dish. After cytogenetic analysis of the surviving colonies, we found chromosomal instability in three of 40 clones analyzed, while some clones exhibited increased micronucleus frequency and HPRT mutation frequency. These data suggest that our chromosomally unstable GM10115 cells secrete factors that are cytotoxic to the majority of stable, parental cells but are also capable of inducing a heritable change in some of the survivors that can manifest as delayed genomic instability. These results suggest a mechanism whereby instability can be perpetuated through the influences of potentially cytotoxic factors produced by genomically unstable clones.     

Genomic structure and sequence variation of a 3.3-kb repeat DNA element of the kangaroo rat, Dipodomys ordii

The genomic structure and sequence variation of a 3.3-kb repeat DNA element, representing 5% of the genome of the kangaroo rat Dipodomys ordii, has been investigated. Most of the repeats are arranged in tandem arrays of 50 kb or more. Thirteen randomly selected genomic clones have been mapped with twelve restriction enzymes. The frequency of sequence divergence in the genomic clones is 0.5%. The clone maps and the genomic structure studies have permitted the characterization of a number of variant members of the 3.3-kb repeat family. The genomic organization of the repeat resembles that for repeated DNAs found in large tandem arrays or satellites.     

Characterisation of set-1, a conserved PR/SET domain gene in Caenorhabditis elegans

We have developed a simple and efficient system (ORF-FINDER) for selecting open reading frames (ORFs) from randomly fragmented genomic DNA fragments. The ORF-FINDER vectors are plasmids that contain a translational start site out of frame with respect to the gene for green fluorescent protein (GFP). Insertion of DNA fragments that bring the initiating ATG in frame with GFP and that contain no stop codons (that is, ORFs) results in the expression of ORF-GFP fusion proteins. In addition, we have developed software (GeneWorks and GenomeAnalyzer) to predict the optimal insert size for maximizing the number of gene-coding ORFs and minimizing unintentionally selected non-coding ORFs. To demonstrate the feasibility of using the ORF-FINDER system to screen genomes for ORFs, we cloned yeast genomic DNA and succeeded in enriching for ORFs by 25-fold. Furthermore, we have shown that the vector can effectively isolate ORFs from the more complex genomes of eukaryotic parasites. We envision that ORF-FINDER will have several applications including genome sequencing projects, gene building from oligonucleotides and construction of expression libraries enriched for ORFs.     

Permutation P-values should never be zero: calculating exact P-values when permutations are randomly drawn

Permutation tests are amongst the most commonly used statistical tools in modern genomic research, a process by which p-values are attached to a test statistic by randomly permuting the sample or gene labels. Yet permutation p-values published in the genomic literature are often computed incorrectly, understated by about 1/m, where m is the number of permutations. The same is often true in the more general situation when Monte Carlo simulation is used to assign p-values. Although the p-value understatement is usually small in absolute terms, the implications can be serious in a multiple testing context. The understatement arises from the intuitive but mistaken idea of using permutation to estimate the tail probability of the test statistic. We argue instead that permutation should be viewed as generating an exact discrete null distribution. The relevant literature, some of which is likely to have been relatively inaccessible to the genomic community, is reviewed and summarized. A computation strategy is developed for exact p-values when permutations are randomly drawn. The strategy is valid for any number of permutations and samples. Some simple recommendations are made for the implementation of permutation tests in practice.     

Restriction fragment polymorphisms as probes for plant diversity and their development as tools for applied plant breeding

Summary Maize and tomato cDNA clones have been hybridized in Southern blotting experiments to plant genomic DNA prepared from different lines to detect restriction fragment polymorphisms (RFPs). In maize we have found that a high degree of genetic variability is present, even among domestic inbred lines. Most randomly chosen maize cDNA clones can be used to detect elements of this variability. Similar levels of polymorphism are observed when genomic DNA is digested with any of a number of different restriction enzymes and probed with individual clones. When a clone is hybridized to genomic DNAs prepared from several different maize lines, a number of different alleles are often detected at a single locus. At the same time one clone can often detect more than one independently segregating locus by cross hybridization to related sequences at other loci. As expected these markers are inherited as simple codominant Mendelian alleles from one generation to the next and colinkage of these markers can be demonstrated in the progeny from a heterozygous parent. In similar studies with tomato, remarkably different results were found. Few RFPs were demonstrable among domestic Lycopersicon esculentum lines although a higher level of variability could be detected when comparing esculentum with its wild Lycopersicon relatives. These results are discussed in relation to the applied uses of RFPs in plant breeding as well as the inherent variability of different plant genomes.This work was supported in part by funds from Sandoz Ltd. (Basel, Switzerland) and its subsidiary company, Northrup King Co. (Minneapolis, Minn., U.S.A.) as well as by NSF SBIR grant #BSR-8360870.     

Partial shotgun sequencing of the Boechera stricta genome reveals extensive microsynteny and promoter conservation with Arabidopsis

Comparative genomics provides insight into the evolutionary dynamics that shape discrete sequences as well as whole genomes. To advance comparative genomics within the Brassicaceae, we have end sequenced 23,136 medium-sized insert clones from Boechera stricta, a wild relative of Arabidopsis (Arabidopsis thaliana). A significant proportion of these sequences, 18,797, are nonredundant and display highly significant similarity (BLASTn e-value < or = 10(-30)) to low copy number Arabidopsis genomic regions, including more than 9,000 annotated coding sequences. We have used this dataset to identify orthologous gene pairs in the two species and to perform a global comparison of DNA regions 5' to annotated coding regions. On average, the 500 nucleotides upstream to coding sequences display 71.4% identity between the two species. In a similar analysis, 61.4% identity was observed between 5' noncoding sequences of Brassica oleracea and Arabidopsis, indicating that regulatory regions are not as diverged among these lineages as previously anticipated. By mapping the B. stricta end sequences onto the Arabidopsis genome, we have identified nearly 2,000 conserved blocks of microsynteny (bracketing 26% of the Arabidopsis genome). A comparison of fully sequenced B. stricta inserts to their homologous Arabidopsis genomic regions indicates that indel polymorphisms >5 kb contribute substantially to the genome size difference observed between the two species. Further, we demonstrate that microsynteny inferred from end-sequence data can be applied to the rapid identification and cloning of genomic regions of interest from nonmodel species. These results suggest that among diploid relatives of Arabidopsis, small- to medium-scale shotgun sequencing approaches can provide rapid and cost-effective benefits to evolutionary and/or functional comparative genomic frameworks.     

一个园艺植物基因组序列和注释的在线批量实时提取平台

目前, 大量园艺植物基因组测序已经完成或接近尾声, 它们的基因组序列和注释数据极大地促进了功能基因组学研究。为给科研人员提供批量下载特定的基因组区段序列和注释平台, 笔者开发了一个称为OBRRP的生物信息学工具。OBRRP具有提取葡萄(Vitis vinifera)、桃(Prunus persica)、草莓(Fragaria vesca)、黄瓜(Cucumis sativus)、西瓜(Citrullus lanatus)、番茄(Solanum lycopersicum)、甜橙(Citrus sinensis)、苹果(Malus x domestica)、猕猴桃(Actinidia chinensis)、马铃薯(Solanum tuberosum)、香蕉(Musa acuminata)和拟南芥(Arabidopsis thaliana) 12种植物基因组序列及注释数据的功能; 同时, 也具有扩展到其它Gbrowser浏览器架构的数据库功能。测试结果表明, OBRRP是一个快捷简便的在线、批量和实时提取工具, 其登录地址为http://bioinfo.jit.edu.cn/OBRRP/。     

Development of a novel bacterial artificial chromosome cloning system for functional studies

Bacterial artificial chromosome (BAC) cloning systems currently in use generate high quality genomic libraries for gene mapping, identification, and sequencing. However, the most commonly used BAC cloning systems do not facilitate functional studies in eukaryotic cells. To overcome this limitation, we have developed pEBAC190G, a new BAC vector that combines the features of the first generation PAC/BAC vectors with eukaryotic elements that facilitate the transfection, episomal maintenance, and functional analysis of large genomic fragments in eukaryotic cells. A number of different cloning strategies may be used to retrofit genomic fragments from existing libraries into the new vector. The system was tested by the retrofitting of a 170kb NotI genomic fragment from the RPCI-11 BAC library into the NotI site of pEBAC190G. Clones from any eukaryotic genomic library harboured in this vector can be transferred from bacteria directly to eukaryotic cells for functional analysis.