Comparison of the accuracy of methods of computational haplotype inference using a large empirical dataset

BACKGROUND: Analyses of genetic data at the level of haplotypes provide increased accuracy and power to infer genotype-phenotype correlations and evolutionary history of a locus. However, empirical determination of haplotypes is expensive and laborious. Therefore, several methods of inferring haplotypes from unphased genotypic data have been proposed, but it is unclear how accurate each of the methods is or which methods are superior. The accuracy of some of the leading methods of computational haplotype inference (PL-EM, Phase, SNPHAP, Haplotyper) are compared using a large set of 308 empirically determined haplotypes based on 15 SNPs, among which 36 haplotypes were observed to occur. This study presents several advantages over many previous comparisons of haplotype inference methods: a large number of subjects are included, the number of known haplotypes is much smaller than the number of chromosomes surveyed, a range in values of linkage disequilibrium, presence of rare SNP alleles, and considerable dispersion in the frequencies of haplotypes. RESULTS: In contrast to some previous comparisons of haplotype inference methods, there was very little difference in the accuracy of the various methods in terms of either assignment of haplotypes to individuals or estimation of haplotype frequencies. Although none of the methods inferred all of the known haplotypes, the assignment of haplotypes to subjects was about 90% correct for individuals heterozygous for up to three SNPs and was about 80% correct for up to five heterozygous sites. All of the methods identified every haplotype with a frequency above 1%, and none assigned a frequency above 1% to an incorrect haplotype. CONCLUSIONS: All of the methods of haplotype inference have high accuracy and one can have confidence in inferences made by any one of the methods. The ability to identify even rare (≥ 1%) haplotypes is reassuring for efforts to identify haplotypes that contribute to disease in a significant proportion of a population. Assignment of haplotypes is relatively accurate among subjects heterozygous for up to 5 sites, and this might be the largest number of SNPs for which one should define haplotype blocks or have confidence in haplotype assignments

Computation of haplotypes on SNPs subsets: advantage of the "global method"

BACKGROUND: Genetic association studies aim at finding correlations between a disease state and genetic variations such as SNPs or combinations of SNPs, termed haplotypes. Some haplotypes have a particular biological meaning such as the ones derived from SNPs located in the promoters, or the ones derived from non synonymous SNPs. All these haplotypes are "subhaplotypes" because they refer only to a part of the SNPs found in the gene. Until now, subhaplotypes were directly computed from the very SNPs chosen to constitute them, without taking into account the rest of the information corresponding to the other SNPs located in the gene. In the present work, we describe an alternative approach, called the "global method", which takes into account all the SNPs known in the region and compare the efficacy of the two "direct" and "global" methods. RESULTS: We used empirical haplotypes data sets from the GH1 promoter and the APOE gene, and 10 simulated datasets, and randomly introduced in them missing information (from 0% up to 20%) to compare the 2 methods. For each method, we used the PHASE haplotyping software since it was described to be the best. We showed that the use of the "global method" for subhaplotyping leads always to a better error rate than the classical direct haplotyping. The advantage provided by this alternative method increases with the percentage of missing genotyping data (diminution of the average error rate from 25% to less than 10%). We applied the global method software on the GRIV cohort for AIDS genetic associations and some associations previously identified through direct subhaplotyping were found to be erroneous. CONCLUSION: The global method for subhaplotyping can reduce, sometimes dramatically, the error rate on patient resolutions and haplotypes frequencies. One should thus use this method in order to minimise the risk of a false interpretation in genetic studies involving subhaplotypes. In practice the global method is always more efficient than the direct method, but a combination method taking into account the level of missing information in each subject appears to be even more interesting when the level of missing information becomes larger (>10%)

Movement of deep-sea coral populations on climatic timescales

During the past 40,000 years, global climate has moved into and out of a full glacial period, with the deglaciation marked by several millennial-scale rapid climate change events. Here we investigate the ecological response of deep-sea coral communities to both glaciation and these rapid climate change events. We find that the deep-sea coral populations of Desmophyllum dianthus in both the North Atlantic and the Tasmanian seamounts expand at times of rapid climate change. However, during the more stable Last Glacial Maximum, the coral population globally retreats to a more restricted depth range. Holocene populations show regional patterns that provide some insight into what causes these dramatic changes in population structure. The most important factors are likely responses to climatically driven changes in productivity, [O_2] and [CO_3^(2–)]

Strong Depth-Related Zonation of Megabenthos on a Rocky Continental Margin (∼700–4000 m) off Southern Tasmania, Australia

Assemblages of megabenthos are structured in seven depth-related zones between ~700 and 4000 m on the rocky and topographically complex continental margin south of Tasmania, southeastern Australia. These patterns emerge from analysis of imagery and specimen collections taken from a suite of surveys using photographic and in situ sampling by epibenthic sleds, towed video cameras, an autonomous underwater vehicle and a remotely operated vehicle (ROV). Seamount peaks in shallow zones had relatively low biomass and low diversity assemblages, which may be in part natural and in part due to effects of bottom trawl fishing. Species richness was highest at intermediate depths (1000–1300 m) as a result of an extensive coral reef community based on the bioherm-forming scleractinian Solenosmilia variabilis. However, megabenthos abundance peaked in a deeper, low diversity assemblage at 2000–2500 m. The S. variabilis reef and the deep biomass zone were separated by an extensive dead, sub-fossil S. variabilis reef and a relatively low biomass stratum on volcanic rock roughly coincident with the oxygen minimum layer. Below 2400 m, megabenthos was increasingly sparse, though punctuated by occasional small pockets of relatively high diversity and biomass. Nonetheless, megabenthic organisms were observed in the vast majority of photographs on all seabed habitats and to the maximum depths observed - a sandy plain below 3950 m. Taxonomic studies in progress suggest that the observed depth zonation is based in part on changing species mixes with depth, but also an underlying commonality to much of the seamount and rocky substrate biota across all depths. Although the mechanisms supporting the extraordinarily high biomass in 2000–2500 m depths remains obscure, plausible explanations include equatorwards lateral transport of polar production and/or a response to depth-stratified oxygen availability

Mammalian mitochondrial DNA evolution: A comparison of the cytochrome b and cytochrome c oxidase II genes

The evolution of two mitochondrial genes, cytochrome b and cytochrome c oxidase subunit II, was examined in several eutherian mammal orders, with special emphasis on the orders Artiodactyla and Rodentia. When analyzed using both maximum parsimony, with either equal or unequal character weighting, and neighbor joining, neither gene performed with a high degree of consistency in terms of the phylogenetic hypotheses supported. The phylogenetic inconsistencies observed for both these genes may be the result of several factors including differences in the rate of nucleotide substitution among particular lineages (especially between orders), base composition bias, transition/transversion bias, differences in codon usage, and different constraints and levels of homoplasy associated with first, second, and third codon positions. We discuss the implications of these findings for the molecular systematics of mammals, especially as they relate to recent hypotheses concerning the polyphyly of the order Rodentia, relationships among the Artiodactyla, and various interordinal relationships.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/48049/1/239_2004_Article_BF00163231.pd

Genome-wide analysis of the RpoN regulon in Geobacter sulfurreducens

Background The role of the RNA polymerase sigma factor RpoN in regulation of gene expression in Geobacter sulfurreducens was investigated to better understand transcriptional regulatory networks as part of an effort to develop regulatory modules for genome-scale in silico models, which can predict the physiological responses of Geobacter species during groundwater bioremediation or electricity production. Results An rpoN deletion mutant could not be obtained under all conditions tested. In order to investigate the regulon of the G. sulfurreducens RpoN, an RpoN over-expression strain was made in which an extra copy of the rpoN gene was under the control of a taclac promoter. Combining both the microarray transcriptome analysis and the computational prediction revealed that the G. sulfurreducens RpoN controls genes involved in a wide range of cellular functions. Most importantly, RpoN controls the expression of the dcuB gene encoding the fumarate/succinate exchanger, which is essential for cell growth with fumarate as the terminal electron acceptor in G. sulfurreducens. RpoN also controls genes, which encode enzymes for both pathways of ammonia assimilation that is predicted to be essential under all growth conditions in G. sulfurreducens. Other genes that were identified as part of the RpoN regulon using either the computational prediction or the microarray transcriptome analysis included genes involved in flagella biosynthesis, pili biosynthesis and genes involved in central metabolism enzymes and cytochromes involved in extracellular electron transfer to Fe(III), which are known to be important for growth in subsurface environment or electricity production in microbial fuel cells. The consensus sequence for the predicted RpoN-regulated promoter elements is TTGGCACGGTTTTTGCT. Conclusion The G. sulfurreducens RpoN is an essential sigma factor and a global regulator involved in a complex transcriptional network controlling a variety of cellular processes

Bioinformatic analysis of gene regulation in Geobacter sulfurreducens

Background Geobacteraceae are a family of microorganisms from the delta subdivision of Proteobacteria. They have potential for environmental bioremediation and electricity generation. In this presentation, we describe our recent bioinformatic analyses of gene regulation in Geobacter sulfurreducens, a model representative of this family. Results and conclusion We have developed an online database, GSEL (Geobacter Sequence Elements), which compiles regulatory information for G. sulfurreducens. We have recently completed the development of a new, significantly expanded and updated, relational version 2 of the GSEL database and its accompanying online query system, which compiles manually curated information on operon organization and transcription regulatory elements in the genome of G. sulfurreducens. GSEL v. 2 incorporates a graphical browser and provides significantly expanded search capabilities. It also includes new information on predicted and/or experimentally validated genome regulatory sites and provides links to information from microarray experiments stored in public gene expression databases, and to original publications describing how particular regulatory interactions were identified. Using sequence and gene expression analyses, we investigated target genes and promoters regulated by RpoN, an alternative RNA polymerase sigma factor, which regulates a variety of important cellular processes in G. sulfurreducens. Our current studies are focusing on an investigation of several transcription regulatory systems involved in RpoN-dependent regulatory pathways. We have investigated target regulatory sites for an enhancer binding protein, PilR, which participates in RpoN-dependent transcriptional regulation of the pilA gene encoding structural pilin. We predicted multiple PilR-regulated sites upstream of operons related to biosynthesis, assembly, and function of pili and flagella, type II secretory pathways, and cell wall biogenesis. We also investigated sequence changes and molecular classification of the TetR family of transcriptional regulators. In G. sulfurreducens, we identified RpoN-regulated promoters upstream of several operons containing tetR family genes. The genome of G. sulfurreducens contains nine tetR family genes. Some of them are located upstream of operons encoding functionally important c-type cytochromes. In order to better understand the roles of TetR family members in the ability of Geobacteraceae to participate in electron transfer, we investigated phylogenetic relationships among TetR proteins in Geobacteraceae and in other microbial species. We identified their conserved and variable domains, which may be important for the diversity of their functional roles, and classified them into subgroups based on sequence similarities

Oceanic evidence of climate change in southern Australia over the last three centuries

Author Posting. © American Geophysical Union, 2004. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Geophysical Research Letters 31 (2004): L07212, doi:10.1029/2003GL018869.Chemical analysis of deepwater octocorals collected at 1000 m depth off southern Australia indicates long-term cooling, beginning in the mid-18th century. This cooling appears to reflect shoaling of isotherms along the continental shelf, that can be related statistically, observationally and by modeling to increasing coastal sea-surface temperatures, that in turn reflect a poleward extension of the SW Pacific boundary current (the East Australian Current). The oceanographic changes implied by the coral record suggest climate change in temperate Australia starting about the time of European settlement. Correlations between temperate Australian and Antarctic indices suggest these long-term changes might also be relevant to Antarctic climate.This study was supported by the Australian Fisheries and Research Development Corporation, the Australian Greenhouse Office, and the Land and Water Research Development Corporation