Search

The draft genome of the Gram-positive spore-forming Anoxybacillus ayderensis strain MT-Cab ( Firmicutes ), isolated from an enrichment culture of Chloracidobacterium thermophilum , was sequenced and comprises 2,577,015 bp in 92 contigs. The draft genome is predicted to consist of 2,699 protein-coding genes, 73 tRNA-coding genes, and an estimated 8 rRNA operons., (Copyright © 2017 Thiel et al.)

Chlorobaculum limnaeum DSM 1677 T is a mesophilic, brown-colored, chlorophototrophic green sulfur bacterium that produces bacteriochlorophyll e and the carotenoid isorenieratene as major pigments. This bacterium serves as a model organism in molecular research on photosynthesis, sulfur metabolism, and bacteriochlorophyll biosynthesis. We report here the complete genome sequence., (Copyright © 2017 Tank et al.)

Background: While clinical outcomes from colorectal cancer (CRC) are influenced by stage at diagnosis and treatment, mounting evidence suggests that an enhanced lymphocytic reaction to a tumor may also be an informative prognostic indicator., Methods: The roles of intratumoral T lymphocyte infiltration (TIL), peritumoral Crohn's-like lymphoid reaction (CLR), microsatellite instability (MSI), and clinicopathological characteristics in survival from CRC were examined using 2369 incident CRCs from a population-based case-control study in northern Israel. Cox proportional hazards regression was used to estimate hazard ratios (HRs) and 95% confidence intervals (CIs) for CRC-specific and all-cause mortality in multivariable models adjusted for age, sex, ethnicity, grade, stage, and MSI. All statistical tests were two-sided., Results: Tumors with TIL/high-powered field (HPF) of 2 or greater were associated with a statistically significant increase in CRC-specific (P < .001) and overall survival (P < .001) compared with tumors with TIL/HPF of less than 2. Similarly, tumors with a prominent CLR experienced better CRC-specific (P < .001) and overall survival (P < .001) as compared with those with no response. High TILs (HR = 0.76, 95% CI = 0.64 to 0.89, P < .001) and a prominent CLR (HR = 0.71, 95% CI = 0.62 to 0.80, P < .001), but not MSI, were associated with a statistically significant reduction in all-cause mortality after adjustment for established prognostic factors., Conclusions: TILs and CLR are both prognostic indicators for CRC after adjusting for traditional prognostic indicators., (© The Author 2016. Published by Oxford University Press. All rights reserved. For Permissions, please e-mail: journals.permissions@oup.com.)

True flies are insects of the order Diptera and encompass one of the most diverse groups of animals on Earth. Within dipterans, Schizophora represents a recent radiation of insects that was used as a model to develop a pipeline for generating complete mitogenomes using various sequencing platforms and strategies. 91 mitogenomes from 32 different species were sequenced and assembled with high fidelity, using amplicon, whole genome shotgun or single molecule sequencing approaches. Based on the novel mitogenomes, we estimate the origin of Schizophora within the Cretaceous-Paleogene (K-Pg) boundary, about 68.3 Ma. Detailed analyses of the blowfly family (Calliphoridae) place its origin at 22 Ma, concomitant with the radiation of grazing mammals. The emergence of ectoparasitism within calliphorids was dated 6.95 Ma for the screwworm fly and 2.3 Ma for the Australian sheep blowfly. Varying population histories were observed for the blowfly Chrysomya megacephala and the housefly Musca domestica samples in our dataset. Whereas blowflies (n = 50) appear to have undergone selective sweeps and/or severe bottlenecks in the New World, houseflies (n = 14) display variation among populations from different zoogeographical zones and low levels of gene flow. The reported high-throughput mitogenomics approach for insects enables new insights into schizophoran diversity and population history of flies.

The draft genome sequence of the Deinococcus-Thermus group bacterium Meiothermus ruber strain A, isolated from a cyanobacterial enrichment culture obtained from Octopus Spring (Yellowstone National Park, WY), comprises 2,968,099 bp in 170 contigs. It is predicted to contain 2,895 protein-coding genes, 44 tRNA-coding genes, and 2 rRNA operons., (Copyright © 2015 Thiel et al.)

Humans first arrived on Madagascar only a few thousand years ago. Subsequent habitat destruction and hunting activities have had significant impacts on the island's biodiversity, including the extinction of megafauna. For example, we know of 17 recently extinct 'subfossil' lemur species, all of which were substantially larger (body mass ∼11-160 kg) than any living population of the ∼100 extant lemur species (largest body mass ∼6.8 kg). We used ancient DNA and genomic methods to study subfossil lemur extinction biology and update our understanding of extant lemur conservation risk factors by i) reconstructing a comprehensive phylogeny of extinct and extant lemurs, and ii) testing whether low genetic diversity is associated with body size and extinction risk. We recovered complete or near-complete mitochondrial genomes from five subfossil lemur taxa, and generated sequence data from population samples of two extinct and eight extant lemur species. Phylogenetic comparisons resolved prior taxonomic uncertainties and confirmed that the extinct subfossil species did not comprise a single clade. Genetic diversity estimates for the two sampled extinct species were relatively low, suggesting small historical population sizes. Low genetic diversity and small population sizes are both risk factors that would have rendered giant lemurs especially susceptible to extinction. Surprisingly, among the extant lemurs, we did not observe a relationship between body size and genetic diversity. The decoupling of these variables suggests that risk factors other than body size may have as much or more meaning for establishing future lemur conservation priorities., (Copyright © 2014 Elsevier Ltd. All rights reserved.)

The draft genome sequence of the thermophilic filamentous anoxygenic phototrophic bacterium Chloroflexus sp. strain MS-G (Chloroflexi), isolated from Mushroom Spring (Yellowstone National Park, WY, USA) was sequenced and comprises 4,784,183 bp in 251 contigs. The draft genome is predicted to encode 4,059 protein coding genes, 49 tRNA encoding genes, and 3 rRNA operons., (Copyright © 2014 Thiel et al.)

Mahoney Lake represents an extreme meromictic model system and is a valuable site for examining the organisms and processes that sustain photic zone euxinia (PZE). A single population of purple sulfur bacteria (PSB) living in a dense phototrophic plate in the chemocline is responsible for most of the primary production in Mahoney Lake. Here, we present metagenomic data from this phototrophic plate--including the genome of the major PSB, as obtained from both a highly enriched culture and from the metagenomic data--as well as evidence for multiple other taxa that contribute to the oxidative sulfur cycle and to sulfate reduction. The planktonic PSB is a member of the Chromatiaceae, here renamed Thiohalocapsa sp. strain ML1. It produces the carotenoid okenone, yet its closest relatives are benthic PSB isolates, a finding that may complicate the use of okenone (okenane) as a biomarker for ancient PZE. Favorable thermodynamics for non-phototrophic sulfide oxidation and sulfate reduction reactions also occur in the plate, and a suite of organisms capable of oxidizing and reducing sulfur is apparent in the metagenome. Fluctuating supplies of both reduced carbon and reduced sulfur to the chemocline may partly account for the diversity of both autotrophic and heterotrophic species. Collectively, the data demonstrate the physiological potential for maintaining complex sulfur and carbon cycles in an anoxic water column, driven by the input of exogenous organic matter. This is consistent with suggestions that high levels of oxygenic primary production maintain episodes of PZE in Earth's history and that such communities should support a diversity of sulfur cycle reactions., (© 2014 John Wiley & Sons Ltd.)

The draft genome sequence of the moderately thermophilic bacterium Schleiferia thermophila strain Yellowstone (Bacteroidetes), isolated from Octopus Spring (Yellowstone National Park, WY, USA) was sequenced and comprises 2,617,694 bp in 35 contigs. The draft genome is predicted to encode 2,457 protein coding genes and 37 tRNA encoding genes and two rRNA operons., (Copyright © 2014 Thiel et al.)

The evolution rate and genetic changes that occur during chronic infection with Helicobacter pylori have been analysed, but little is known about the genomic changes during the initial, acute bacterial infection phase. Here we analyse the rate and pattern of genome evolution in H. pylori from the genomes of two input strains isolated from human volunteers with asymptomatic infection, and the genomes of two output strains collected 20 and 44 days after re-infection. Similarly, we analyse genome evolution in bacteria from the genome sequences of input and output strains sequentially taken after experimental infection of a rhesus macaque. The estimated mutation rate reveals a mutation burst during the acute infection phase that is over 10 times faster than the mutation rate during chronic infection, and orders of magnitude faster than mutation rates in any other bacteria. The elevated frequency of mutations in outer membrane protein genes suggests that the mutation burst facilitates rapid host adaptation of the bacteria.

The genome of the unicellular cyanobacterium Thermosynechococcus sp. strain NK55a, isolated from the Nakabusa hot spring, Nagano Prefecture, Japan, comprises a single, circular, 2.5-Mb chromosome. The genome is predicted to contain 2,358 protein-encoding genes, including genes for all typical cyanobacterial photosynthetic and metabolic functions. No genes encoding hydrogenases or nitrogenase were identified.

Background: 'Chlorochromatium aggregatum' is a phototrophic consortium, a symbiosis that may represent the highest degree of mutual interdependence between two unrelated bacteria not associated with a eukaryotic host. 'Chlorochromatium aggregatum' is a motile, barrel-shaped aggregate formed from a single cell of 'Candidatus Symbiobacter mobilis", a polarly flagellated, non-pigmented, heterotrophic bacterium, which is surrounded by approximately 15 epibiont cells of Chlorobium chlorochromatii, a non-motile photolithoautotrophic green sulfur bacterium., Results: We analyzed the complete genome sequences of both organisms to understand the basis for this symbiosis. Chl. chlorochromatii has acquired relatively few symbiosis-specific genes; most acquired genes are predicted to modify the cell wall or function in cell-cell adhesion. In striking contrast, 'Ca. S. mobilis' appears to have undergone massive gene loss, is probably no longer capable of independent growth, and thus may only reproduce when consortia divide. A detailed model for the energetic and metabolic bases of the dependency of 'Ca. S. mobilis' on Chl. chlorochromatii is described., Conclusions: Genomic analyses suggest that three types of interactions lead to a highly sophisticated relationship between these two organisms. Firstly, extensive metabolic exchange, involving carbon, nitrogen, and sulfur sources as well as vitamins, occurs from the epibiont to the central bacterium. Secondly, 'Ca. S. mobilis' can sense and move towards light and sulfide, resources that only directly benefit the epibiont. Thirdly, electron cycling mechanisms, particularly those mediated by quinones and potentially involving shared protonmotive force, could provide an important basis for energy exchange in this and other symbiotic relationships.

Polar bears (PBs) are superbly adapted to the extreme Arctic environment and have become emblematic of the threat to biodiversity from global climate change. Their divergence from the lower-latitude brown bear provides a textbook example of rapid evolution of distinct phenotypes. However, limited mitochondrial and nuclear DNA evidence conflicts in the timing of PB origin as well as placement of the species within versus sister to the brown bear lineage. We gathered extensive genomic sequence data from contemporary polar, brown, and American black bear samples, in addition to a 130,000- to 110,000-y old PB, to examine this problem from a genome-wide perspective. Nuclear DNA markers reflect a species tree consistent with expectation, showing polar and brown bears to be sister species. However, for the enigmatic brown bears native to Alaska's Alexander Archipelago, we estimate that not only their mitochondrial genome, but also 5-10% of their nuclear genome, is most closely related to PBs, indicating ancient admixture between the two species. Explicit admixture analyses are consistent with ancient splits among PBs, brown bears and black bears that were later followed by occasional admixture. We also provide paleodemographic estimates that suggest bear evolution has tracked key climate events, and that PB in particular experienced a prolonged and dramatic decline in its effective population size during the last ca. 500,000 years. We demonstrate that brown bears and PBs have had sufficiently independent evolutionary histories over the last 4-5 million years to leave imprints in the PB nuclear genome that likely are associated with ecological adaptation to the Arctic environment.

A degenerate polymerase chain reaction (PCR)-based method of whole-genome amplification, designed to work fluidly with 454 sequencing technology, was developed and tested for use on deep marine subsurface DNA samples. While optimized here for use with Roche 454 technology, the general framework presented may be applicable to other next generation sequencing systems as well (e.g., Illumina, Ion Torrent). The method, which we have called random amplification metagenomic PCR (RAMP), involves the use of specific primers from Roche 454 amplicon sequencing, modified by the addition of a degenerate region at the 3' end. It utilizes a PCR reaction, which resulted in no amplification from blanks, even after 50 cycles of PCR. After efforts to optimize experimental conditions, the method was tested with DNA extracted from cultured E. coli cells, and genome coverage was estimated after sequencing on three different occasions. Coverage did not vary greatly with the different experimental conditions tested, and was around 62% with a sequencing effort equivalent to a theoretical genome coverage of 14.10×. The GC content of the sequenced amplification product was within 2% of the predicted values for this strain of E. coli. The method was also applied to DNA extracted from marine subsurface samples from ODP Leg 201 site 1229 (Peru Margin), and results of a taxonomic analysis revealed microbial communities dominated by Proteobacteria, Chloroflexi, Firmicutes, Euryarchaeota, and Crenarchaeota, among others. These results were similar to those obtained previously for those samples; however, variations in the proportions of taxa identified illustrates well the generally accepted view that community analysis is sensitive to both the amplification technique used and the method of assigning sequences to taxonomic groups. Overall, we find that RAMP represents a valid methodology for amplifying metagenomes from low-biomass samples.

Candidatus Chloracidobacterium thermophilum, which naturally inhabits microbial mats of alkaline siliceous hot springs in Yellowstone National Park, is the only known chlorophototroph in the phylum Acidobacteria. The Ca. C. thermophilum genome was composed of two chromosomes (2,683,362 bp and 1,012,010 bp), and both encoded essential genes. The genome included genes to produce chlorosomes, the Fenna-Matthews-Olson protein, bacteriochlorophylls a and c as principal pigments, and type-1, homodimeric reaction centres. Ca. C. thermophilum is an aerobic photoheterotroph that lacks the ability to synthesize several essential nutrients. Key genes of all known carbon fixation pathways were absent, as were genes for assimilatory nitrate and sulfate reduction and vitamin B(12) synthesis. Genes for the synthesis of branched-chain amino acids (valine, isoleucine and leucine) were also absent, but genes for catabolism of these compounds were present. This observation suggested that Ca. C. thermophilum may synthesize branched-chain amino acids from an intermediate(s) of the catabolic pathway by reversing these reactions. The genome encoded an aerobic respiratory electron transport chain that included NADH dehydrogenase, alternative complex III and cytochrome oxidase. The chromosomes of the laboratory isolate were compared with assembled, metagenomic scaffolds from the major Ca. C. thermophilum population in hot-spring mats. The larger chromosomes of the two populations were highly syntenous but significantly divergent (~13%) in sequence. In contrast, the smaller chromosomes have undergone numerous rearrangements, contained many transposases, and might be less constrained by purifying selection than the large chromosomes. Some transposases were homologous to those of mat community members from other phyla., (© 2011 Society for Applied Microbiology and Blackwell Publishing Ltd.)

The phototrophic microbial mat community of Mushroom Spring, an alkaline siliceous hot spring in Yellowstone National Park, was studied by metatranscriptomic methods. RNA was extracted from mat specimens collected at four timepoints during light-to-dark and dark-to-light transitions in one diel cycle, and these RNA samples were analyzed by both pyrosequencing and SOLiD technologies. Pyrosequencing was used to assess the community composition, which showed that ~84% of the rRNA was derived from members of four kingdoms Cyanobacteria, Chloroflexi, Chlorobi and Acidobacteria. Transcription of photosynthesis-related genes conclusively demonstrated the phototrophic nature of two newly discovered populations; these organisms, which were discovered through metagenomics, are currently uncultured and previously undescribed members of Chloroflexi and Chlorobi. Data sets produced by SOLiD sequencing of complementary DNA provided >100-fold greater sequence coverage. The much greater sequencing depth allowed transcripts to be detected from ~15,000 genes and could be used to demonstrate statistically significant differential transcription of thousands of genes. Temporal differences for in situ transcription patterns of photosynthesis-related genes suggested that the six types of chlorophototrophs in the mats may use different strategies for maximizing their solar-energy capture, usage and growth. On the basis of both temporal pattern and transcript abundance, intra-guild gene expression differences were also detected for two populations of the oxygenic photosynthesis guild. This study showed that, when community-relevant genomes and metagenomes are available, SOLiD sequencing technology can be used for metatranscriptomic analyses, and the results suggested that this method can potentially reveal new insights into the ecophysiology of this model microbial community.

The Tasmanian devil (Sarcophilus harrisii) is threatened with extinction because of a contagious cancer known as Devil Facial Tumor Disease. The inability to mount an immune response and to reject these tumors might be caused by a lack of genetic diversity within a dwindling population. Here we report a whole-genome analysis of two animals originating from extreme northwest and southeast Tasmania, the maximal geographic spread, together with the genome from a tumor taken from one of them. A 3.3-Gb de novo assembly of the sequence data from two complementary next-generation sequencing platforms was used to identify 1 million polymorphic genomic positions, roughly one-quarter of the number observed between two genetically distant human genomes. Analysis of 14 complete mitochondrial genomes from current and museum specimens, as well as mitochondrial and nuclear SNP markers in 175 animals, suggests that the observed low genetic diversity in today's population preceded the Devil Facial Tumor Disease disease outbreak by at least 100 y. Using a genetically characterized breeding stock based on the genome sequence will enable preservation of the extant genetic diversity in future Tasmanian devil populations.

Whole-genome duplication (WGD), or polyploidy, followed by gene loss and diploidization has long been recognized as an important evolutionary force in animals, fungi and other organisms, especially plants. The success of angiosperms has been attributed, in part, to innovations associated with gene or whole-genome duplications, but evidence for proposed ancient genome duplications pre-dating the divergence of monocots and eudicots remains equivocal in analyses of conserved gene order. Here we use comprehensive phylogenomic analyses of sequenced plant genomes and more than 12.6 million new expressed-sequence-tag sequences from phylogenetically pivotal lineages to elucidate two groups of ancient gene duplications-one in the common ancestor of extant seed plants and the other in the common ancestor of extant angiosperms. Gene duplication events were intensely concentrated around 319 and 192 million years ago, implicating two WGDs in ancestral lineages shortly before the diversification of extant seed plants and extant angiosperms, respectively. Significantly, these ancestral WGDs resulted in the diversification of regulatory genes important to seed and flower development, suggesting that they were involved in major innovations that ultimately contributed to the rise and eventual dominance of seed plants and angiosperms., (©2011 Macmillan Publishers Limited. All rights reserved)

Assessment of the microbial diversity residing in arthropod vectors of medical importance is crucial for monitoring endemic infections, for surveillance of newly emerging zoonotic pathogens, and for unraveling the associated bacteria within its host. The tick Ixodes ricinus is recognized as the primary European vector of disease-causing bacteria in humans. Despite I. ricinus being of great public health relevance, its microbial communities remain largely unexplored to date. Here we evaluate the pathogen-load and the microbiome in single adult I. ricinus by using 454- and Illumina-based metagenomic approaches. Genomic DNA-derived sequences were taxonomically profiled using a computational approach based on the BWA algorithm, allowing for the identification of known tick-borne pathogens at the strain level and the putative tick core microbiome. Additionally, we assessed and compared the bacterial taxonomic profile in nymphal and adult I. ricinus pools collected from two distinct geographic regions in Northern Italy by means of V6-16S rRNA amplicon pyrosequencing and community based ecological analysis. A total of 108 genera belonging to representatives of all bacterial phyla were detected and a rapid qualitative assessment for pathogenic bacteria, such as Borrelia, Rickettsia and Candidatus Neoehrlichia, and for other bacteria with mutualistic relationship or undetermined function, such as Wolbachia and Rickettsiella, was possible. Interestingly, the ecological analysis revealed that the bacterial community structure differed between the examined geographic regions and tick life stages. This finding suggests that the environmental context (abiotic and biotic factors) and host-selection behaviors affect their microbiome.Our data provide the most complete picture to date of the bacterial communities present within I. ricinus under natural conditions by using high-throughput sequencing technologies. This study further demonstrates a novel detection strategy for the microbiomes of arthropod vectors in the context of epidemiological and ecological studies.

The ongoing emergence of human infections originating from wildlife highlights the need for better knowledge of the microbial community in wildlife species where traditional diagnostic approaches are limited. Here we evaluate the microbial biota in healthy mule deer (Odocoileus hemionus) by analyses of lymph node meta-transcriptomes. cDNA libraries from five individuals and two pools of samples were prepared from retropharyngeal lymph node RNA enriched for polyadenylated RNA and sequenced using Roche-454 Life Sciences technology. Protein-coding and 16S ribosomal RNA (rRNA) sequences were taxonomically profiled using protein and rRNA specific databases. Representatives of all bacterial phyla were detected in the seven libraries based on protein-coding transcripts indicating that viable microbiota were present in lymph nodes. Residents of skin and rumen, and those ubiquitous in mule deer habitat dominated classifiable bacterial species. Based on detection of both rRNA and protein-coding transcripts, we identified two new proteobacterial species; a Helicobacter closely related to Helicobacter cetorum in the Helicobacter pylori/Helicobacter acinonychis complex and an Acinetobacter related to Acinetobacter schindleri. Among viruses, a novel gamma retrovirus and other members of the Poxviridae and Retroviridae were identified. We additionally evaluated bacterial diversity by amplicon sequencing the hypervariable V6 region of 16S rRNA and demonstrate that overall taxonomic diversity is higher with the meta-transcriptomic approach. These data provide the most complete picture to date of the microbial diversity within a wildlife host. Our research advances the use of meta-transcriptomics to study microbiota in wildlife tissues, which will facilitate detection of novel organisms with pathogenic potential to human and animals.

The systematic characterization of somatic mutations in cancer genomes is essential for understanding the disease and for developing targeted therapeutics. Here we report the identification of 2,576 somatic mutations across approximately 1,800 megabases of DNA representing 1,507 coding genes from 441 tumours comprising breast, lung, ovarian and prostate cancer types and subtypes. We found that mutation rates and the sets of mutated genes varied substantially across tumour types and subtypes. Statistical analysis identified 77 significantly mutated genes including protein kinases, G-protein-coupled receptors such as GRM8, BAI3, AGTRL1 (also called APLNR) and LPHN3, and other druggable targets. Integrated analysis of somatic mutations and copy number alterations identified another 35 significantly altered genes including GNAS, indicating an expanded role for galpha subunits in multiple cancer types. Furthermore, our experimental analyses demonstrate the functional roles of mutant GNAO1 (a Galpha subunit) and mutant MAP2K4 (a member of the JNK signalling pathway) in oncogenesis. Our study provides an overview of the mutational spectra across major human cancers and identifies several potential therapeutic targets.

The polar bear has become the flagship species in the climate-change discussion. However, little is known about how past climate impacted its evolution and persistence, given an extremely poor fossil record. Although it is undisputed from analyses of mitochondrial (mt) DNA that polar bears constitute a lineage within the genetic diversity of brown bears, timing estimates of their divergence have differed considerably. Using next-generation sequencing technology, we have generated a complete, high-quality mt genome from a stratigraphically validated 130,000- to 110,000-year-old polar bear jawbone. In addition, six mt genomes were generated of extant polar bears from Alaska and brown bears from the Admiralty and Baranof islands of the Alexander Archipelago of southeastern Alaska and Kodiak Island. We show that the phylogenetic position of the ancient polar bear lies almost directly at the branching point between polar bears and brown bears, elucidating a unique morphologically and molecularly documented fossil link between living mammal species. Molecular dating and stable isotope analyses also show that by very early in their evolutionary history, polar bears were already inhabitants of the Artic sea ice and had adapted very rapidly to their current and unique ecology at the top of the Arctic marine food chain. As such, polar bears provide an excellent example of evolutionary opportunism within a widespread mammalian lineage.

The genetic structure of the indigenous hunter-gatherer peoples of southern Africa, the oldest known lineage of modern human, is important for understanding human diversity. Studies based on mitochondrial and small sets of nuclear markers have shown that these hunter-gatherers, known as Khoisan, San, or Bushmen, are genetically divergent from other humans. However, until now, fully sequenced human genomes have been limited to recently diverged populations. Here we present the complete genome sequences of an indigenous hunter-gatherer from the Kalahari Desert and a Bantu from southern Africa, as well as protein-coding regions from an additional three hunter-gatherers from disparate regions of the Kalahari. We characterize the extent of whole-genome and exome diversity among the five men, reporting 1.3 million novel DNA differences genome-wide, including 13,146 novel amino acid variants. In terms of nucleotide substitutions, the Bushmen seem to be, on average, more different from each other than, for example, a European and an Asian. Observed genomic differences between the hunter-gatherers and others may help to pinpoint genetic adaptations to an agricultural lifestyle. Adding the described variants to current databases will facilitate inclusion of southern Africans in medical research efforts, particularly when family and medical histories can be correlated with genome-wide data.

Background: Meiotic recombination alters frequency and distribution of genetic variation, impacting genetics and evolution. In the budding yeast, DNA double strand breaks (DSBs) and D loops form either crossovers (COs) or non-crossovers (NCOs), which occur at many sites in the genome. Differences at the nucleotide level associated with COs and NCOs enable us to detect these recombination events and their distributions., Results: We used high throughput sequencing to uncover over 46 thousand single nucleotide polymorphisms (SNPs) between two budding yeast strains and investigated meiotic recombinational events. We provided a detailed analysis of CO and NCO events, including number, size range, and distribution on chromosomes. We have detected 91 COs, very close to the average number from previous genetic studies, as well as 21 NCO events and mapped the positions of these events with high resolution. We have obtained DNA sequence-level evidence for a wide range of sizes of chromosomal regions involved in CO and NCO events. We show that a large fraction of the COs are accompanied by gene conversion (GC), indicating that meiotic recombination changes allelic frequencies, in addition to redistributing existing genetic variations., Conclusion: This work is the first reported study of meiotic recombination using high throughput sequencing technologies. Our results show that high-throughput sequencing is a sensitive method to uncover at single-base resolution details of CO and NCO events, including some complex patterns, providing new clues about the mechanism of this fundamental process.

Background: We have developed a simulation approach to help determine the optimal mixture of sequencing methods for most complete and cost effective transcriptome sequencing. We compared simulation results for traditional capillary sequencing with "Next Generation" (NG) ultra high-throughput technologies. The simulation model was parameterized using mappings of 130,000 cDNA sequence reads to the Arabidopsis genome (NCBI Accession SRA008180.19). We also generated 454-GS20 sequences and de novo assemblies for the basal eudicot California poppy (Eschscholzia californica) and the magnoliid avocado (Persea americana) using a variety of methods for cDNA synthesis., Results: The Arabidopsis reads tagged more than 15,000 genes, including new splice variants and extended UTR regions. Of the total 134,791 reads (13.8 MB), 119,518 (88.7%) mapped exactly to known exons, while 1,117 (0.8%) mapped to introns, 11,524 (8.6%) spanned annotated intron/exon boundaries, and 3,066 (2.3%) extended beyond the end of annotated UTRs. Sequence-based inference of relative gene expression levels correlated significantly with microarray data. As expected, NG sequencing of normalized libraries tagged more genes than non-normalized libraries, although non-normalized libraries yielded more full-length cDNA sequences. The Arabidopsis data were used to simulate additional rounds of NG and traditional EST sequencing, and various combinations of each. Our simulations suggest a combination of FLX and Solexa sequencing for optimal transcriptome coverage at modest cost. We have also developed ESTcalc http://fgp.huck.psu.edu/NG&#95;Sims/ngsim.pl, an online webtool, which allows users to explore the results of this study by specifying individualized costs and sequencing characteristics., Conclusion: NG sequencing technologies are a highly flexible set of platforms that can be scaled to suit different project goals. In terms of sequence coverage alone, the NG sequencing is a dramatic advance over capillary-based sequencing, but NG sequencing also presents significant challenges in assembly and sequence accuracy due to short read lengths, method-specific sequencing errors, and the absence of physical clones. These problems may be overcome by hybrid sequencing strategies using a mixture of sequencing methodologies, by new assemblers, and by sequencing more deeply. Sequencing and microarray outcomes from multiple experiments suggest that our simulator will be useful for guiding NG transcriptome sequencing projects in a wide range of organisms.

Background: The scientific literature contains many examples where DNA sequence analyses have been used to provide definitive answers to phylogenetic problems that traditional (non-DNA based) approaches alone have failed to resolve. One notable example concerns the rhinoceroses, a group for which several contradictory phylogenies were proposed on the basis of morphology, then apparently resolved using mitochondrial DNA fragments., Results: In this study we report the first complete mitochondrial genome sequences of the extinct ice-age woolly rhinoceros (Coelodonta antiquitatis), and the threatened Javan (Rhinoceros sondaicus), Sumatran (Dicerorhinus sumatrensis), and black (Diceros bicornis) rhinoceroses. In combination with the previously published mitochondrial genomes of the white (Ceratotherium simum) and Indian (Rhinoceros unicornis) rhinoceroses, this data set putatively enables reconstruction of the rhinoceros phylogeny. While the six species cluster into three strongly supported sister-pairings: (i) The black/white, (ii) the woolly/Sumatran, and (iii) the Javan/Indian, resolution of the higher-level relationships has no statistical support. The phylogenetic signal from individual genes is highly diffuse, with mixed topological support from different genes. Furthermore, the choice of outgroup (horse vs tapir) has considerable effect on reconstruction of the phylogeny. The lack of resolution is suggestive of a hard polytomy at the base of crown-group Rhinocerotidae, and this is supported by an investigation of the relative branch lengths., Conclusion: Satisfactory resolution of the rhinoceros phylogeny may not be achievable without additional analyses of substantial amounts of nuclear DNA. This study provides a compelling demonstration that, in spite of substantial sequence length, there are significant limitations with single-locus phylogenetics. We expect further examples of this to appear as next-generation, large-scale sequencing of complete mitochondrial genomes becomes commonplace in evolutionary studies. "The human factor in classification is nowhere more evident than in dealing with this superfamily (Rhinocerotoidea)." G. G. Simpson (1945).

We report the first two complete mitochondrial genome sequences of the thylacine (Thylacinus cynocephalus), or so-called Tasmanian tiger, extinct since 1936. The thylacine's phylogenetic position within australidelphian marsupials has long been debated, and here we provide strong support for the thylacine's basal position in Dasyuromorphia, aided by mitochondrial genome sequence that we generated from the extant numbat (Myrmecobius fasciatus). Surprisingly, both of our thylacine sequences differ by 11%-15% from putative thylacine mitochondrial genes in GenBank, with one of our samples originating from a direct offspring of the previously sequenced individual. Our data sample each mitochondrial nucleotide an average of 50 times, thereby providing the first high-fidelity reference sequence for thylacine population genetics. Our two sequences differ in only five nucleotides out of 15,452, hinting at a very low genetic diversity shortly before extinction. Despite the samples' heavy contamination with bacterial and human DNA and their temperate storage history, we estimate that as much as one-third of the total DNA in each sample is from the thylacine. The microbial content of the two thylacine samples was subjected to metagenomic analysis, and showed striking differences between a wild-captured individual and a born-in-captivity one. This study therefore adds to the growing evidence that extensive sequencing of museum collections is both feasible and desirable, and can yield complete genomes.

Identification of microsatellite unstable (MSI-H) colorectal cancers (CRCs) is important not only for the identification of hereditary nonpolyposis colorectal cancer syndrome but also because MSI-H CRCs have a better prognosis and may respond differently to 5-fluorouracil-based chemotherapy. We present 2 nearly equivalent logistic regression models for clinical use that predict microsatellite instability based on the review of 1649 CRCs from patients of all ages collected in a population-based case control study in northern Israel. One hundred ninety-eight of these 1649 tumors demonstrated a high degree of microsatellite instability (12%). Multivariate analysis found that >2 tumor-infiltrating lymphocyte (TIL) cells per high-powered field, the lack of dirty necrosis, the presence of a Crohn-like reaction, right-sided location, any mucinous differentiation (mucinous or focally mucinous) and well or poor differentiation, and age less than 50 were all independent predictors of MSI-H. We developed 2 logistic regression models that differ only by the statistical approach used to analyze the number of TIL cells per high-powered field, where the slightly more accurate (and complex) model uses the log of the total number of TIL cells. The simpler clinical model uses a cut-off of 2>TIL cells per high-powered field. The accuracy of both models is high, with an 85.4% versus 85.0% probability of correctly classifying tumors as MSI-H. By employing the simpler model, pathologists can predict the likelihood of microsatellite instability by compiling the MSI probability score (Table 4 and Fig. 1) from simple histologic and clinical data available during sign-out. Our model shows that approximately 43% of CRCs have a MSI probability score of 1 or less and hence have little likelihood (<3%) of being MSI-H. Although this model is not perfect in predicting microsatellite instability, its use could improve the efficiency of expensive diagnostic testing.

In 1994, two independent groups extracted DNA from several Pleistocene epoch mammoths and noted differences among individual specimens. Subsequently, DNA sequences have been published for a number of extinct species. However, such ancient DNA is often fragmented and damaged, and studies to date have typically focused on short mitochondrial sequences, never yielding more than a fraction of a per cent of any nuclear genome. Here we describe 4.17 billion bases (Gb) of sequence from several mammoth specimens, 3.3 billion (80%) of which are from the woolly mammoth (Mammuthus primigenius) genome and thus comprise an extensive set of genome-wide sequence from an extinct species. Our data support earlier reports that elephantid genomes exceed 4 Gb. The estimated divergence rate between mammoth and African elephant is half of that between human and chimpanzee. The observed number of nucleotide differences between two particular mammoths was approximately one-eighth of that between one of them and the African elephant, corresponding to a separation between the mammoths of 1.5-2.0 Myr. The estimated probability that orthologous elephant and mammoth amino acids differ is 0.002, corresponding to about one residue per protein. Differences were discovered between mammoth and African elephant in amino-acid positions that are otherwise invariant over several billion years of combined mammalian evolution. This study shows that nuclear genome sequencing of extinct species can reveal population differences not evident from the fossil record, and perhaps even discover genetic factors that affect extinction.

Most nucleosomes are well-organized at the 5' ends of S. cerevisiae genes where "-1" and "+1" nucleosomes bracket a nucleosome-free promoter region (NFR). How nucleosomal organization is specified by the genome is less clear. Here we establish and inter-relate rules governing genomic nucleosome organization by sequencing DNA from more than one million immunopurified S. cerevisiae nucleosomes (displayed at http://atlas.bx.psu.edu/). Evidence is presented that the organization of nucleosomes throughout genes is largely a consequence of statistical packing principles. The genomic sequence specifies the location of the -1 and +1 nucleosomes. The +1 nucleosome forms a barrier against which nucleosomes are packed, resulting in uniform positioning, which decays at farther distances from the barrier. We present evidence for a novel 3' NFR that is present at >95% of all genes. 3' NFRs may be important for transcription termination and anti-sense initiation. We present a high-resolution genome-wide map of TFIIB locations that implicates 3' NFRs in gene looping.

We report five new complete mitochondrial DNA (mtDNA) genomes of Siberian woolly mammoth (Mammuthus primigenius), sequenced with up to 73-fold coverage from DNA extracted from hair shaft material. Three of the sequences present the first complete mtDNA genomes of mammoth clade II. Analysis of these and 13 recently published mtDNA genomes demonstrates the existence of two apparently sympatric mtDNA clades that exhibit high interclade divergence. The analytical power afforded by the analysis of the complete mtDNA genomes reveals a surprisingly ancient coalescence age of the two clades, approximately 1-2 million years, depending on the calibration technique. Furthermore, statistical analysis of the temporal distribution of the (14)C ages of these and previously identified members of the two mammoth clades suggests that clade II went extinct before clade I. Modeling of protein structures failed to indicate any important functional difference between genomes belonging to the two clades, suggesting that the loss of clade II more likely is due to genetic drift than a selective sweep.

Comparative genomics of nucleosome positions provides a powerful means for understanding how the organization of chromatin and the transcription machinery co-evolve. Here we produce a high-resolution reference map of H2A.Z and bulk nucleosome locations across the genome of the fly Drosophila melanogaster and compare it to that from the yeast Saccharomyces cerevisiae. Like Saccharomyces, Drosophila nucleosomes are organized around active transcription start sites in a canonical -1, nucleosome-free region, +1 arrangement. However, Drosophila does not incorporate H2A.Z into the -1 nucleosome and does not bury its transcriptional start site in the +1 nucleosome. At thousands of genes, RNA polymerase II engages the +1 nucleosome and pauses. How the transcription initiation machinery contends with the +1 nucleosome seems to be fundamentally different across major eukaryotic lines.

Although the application of sequencing-by-synthesis techniques to DNA extracted from bones has revolutionized the study of ancient DNA, it has been plagued by large fractions of contaminating environmental DNA. The genetic analyses of hair shafts could be a solution: We present 10 previously unexamined Siberian mammoth (Mammuthus primigenius) mitochondrial genomes, sequenced with up to 48-fold coverage. The observed levels of damage-derived sequencing errors were lower than those observed in previously published frozen bone samples, even though one of the specimens was >50,000 14C years old and another had been stored for 200 years at room temperature. The method therefore sets the stage for molecular-genetic analysis of museum collections.

The nucleosome is the fundamental building block of eukaryotic chromosomes. Access to genetic information encoded in chromosomes is dependent on the position of nucleosomes along the DNA. Alternative locations just a few nucleotides apart can have profound effects on gene expression. Yet the nucleosomal context in which chromosomal and gene regulatory elements reside remains ill-defined on a genomic scale. Here we sequence the DNA of 322,000 individual Saccharomyces cerevisiae nucleosomes, containing the histone variant H2A.Z, to provide a comprehensive map of H2A.Z nucleosomes in functionally important regions. With a median 4-base-pair resolution, we identify new and established signatures of nucleosome positioning. A single predominant rotational setting and multiple translational settings are evident. Chromosomal elements, ranging from telomeres to centromeres and transcriptional units, are found to possess characteristic nucleosomal architecture that may be important for their function. Promoter regulatory elements, including transcription factor binding sites and transcriptional start sites, show topological relationships with nucleosomes, such that transcription factor binding sites tend to be rotationally exposed on the nucleosome surface near its border. Transcriptional start sites tended to reside about one helical turn inside the nucleosome border. These findings reveal an intimate relationship between chromatin architecture and the underlying DNA sequence it regulates.

The kin-cohort design has been proposed as an alternative to traditional case-control and cohort measures to evaluate inherited susceptibility to cancer in population-based studies. Here, we used this design to evaluate inherited susceptibility to prostate cancer associated with APC I1307K using data from the Molecular Epidemiology of Colorectal Cancer study. Two techniques were used to compare the incidence of prostate cancer in APC I1307K carriers. First, we compared the incidence of prostate cancer in relatives of mutation carriers and noncarriers using standard techniques for survival analysis. Second, we used the marginal maximum likelihood method for kin-cohort analysis to infer the genotypes in the relatives. We also evaluated APC I1307K in 75 Ashkenazi Jewish individuals with prostate cancer from 27 families enrolled in the University of Michigan Prostate Cancer Genetic Study. We observed a slightly increased risk of prostate cancer in relatives of APC I1307K carriers, however, this difference was not statistically significant (hazard ratio, 1.6; 95% confidence intervals, 0.7-3.4). Similar conclusions were drawn using both techniques for kin-cohort analysis. APC I1307K was found in 7.4% of families genotyped, which is slightly higher than the allele prevalence reported in Ashkenazi Jews in the general population. In addition, we did not observe loss of heterozygosity at APC or a somatic mutation near APC I1307K using microdissected tumor DNA from mutation carriers enrolled in the Prostate Cancer Genetic Study. Overall, the evidence for an association between APC I1307K and prostate cancer is not compelling. APC I1307K is unlikely to play a clinically meaningful role in susceptibility to prostate cancer.

We sequenced 28 million base pairs of DNA in a metagenomics approach, using a woolly mammoth (Mammuthus primigenius) sample from Siberia. As a result of exceptional sample preservation and the use of a recently developed emulsion polymerase chain reaction and pyrosequencing technique, 13 million base pairs (45.4%) of the sequencing reads were identified as mammoth DNA. Sequence identity between our data and African elephant (Loxodonta africana) was 98.55%, consistent with a paleontologically based divergence date of 5 to 6 million years. The sample includes a surprisingly small diversity of environmental DNAs. The high percentage of endogenous DNA recoverable from this single mammoth would allow for completion of its genome, unleashing the field of paleogenomics.

Purpose: The probability of colorectal cancer is moderately increased among carriers of the APC I1307K polymorphism. However, it is not known if endoscopic surveillance of this high-risk group is warranted. The prevalence of polyps and adenomas in specimens of colorectal cancer who are carriers and noncarriers of the APC I1307K polymorphism is compared., Method: Prevalence of adenomatous polyps in the pathology specimens of the study participants, stratified by their APC I1307K polymorphism status, was studied in 900 consecutive cases of colorectal cancer diagnosed in northern Israel between 1998 and 2002, within the framework of a population-based, case-controlled study (MECC Study)., Results: The APC I1307K mutation was detected in 78 colorectal cancer cases (8.7 percent) of the study population. Prevalence was higher among Ashkenazi Jews (11.2 percent) than among non-Ashkenazi Jews (2.7 percent) or Arabs (3.1 percent). After adjustment for age, APC I1307K carriers were significantly more likely than noncarriers to have polyps in their surgical specimen (51.3 percent vs. 32.6 percent, P = 0.002). Adenomas with a tubular component (either tubular adenomas or tubulovillous adenomas), but not villous adenomas, were significantly more frequent among carriers (37.2 percent vs. 23.6 percent, P = 0.005)., Conclusion: Together with former evidence of I1307K being a risk factor for colorectal cancer, these data suggest that colonoscopic surveillance for colorectal adenomas and cancer may be warranted in I1307K carriers, even in the absence of other identifiable risk factors.

Overexpression of the epidermal growth factor (EGF) pathway has been implicated in melanoma pathogenesis, and a recent case-control study identified a single nucleotide polymorphism (G to A) in the EGF gene where the G allele was associated with increased EGF expression and an increased risk of melanoma. To further evaluate this association, we conducted a case-control analysis from the Genes, Environment, and Melanoma study at the University of Michigan site using two different study designs. Incident cases of histopathologically confirmed first primary melanoma that were diagnosed between January 1, 2000 and December 31, 2000 from the University of Michigan Melanoma Clinic (n = 330) were compared with the following two different sources of nonmelanoma controls: spouse/friend controls (n = 84) and healthy volunteer controls from a case-control study of psoriasis (n = 148). Using a second analytic design, comparisons between multiple primary melanoma cases (n = 62) and single primary melanoma cases (n = 330) were also evaluated to estimate odds ratios (ORs). Genotyping for the single nucleotide substitution (G to A) at position 61 in the 5' untranslated region of the EGF gene was performed from genomic DNA, and epidemiological risk factors were assessed through a telephone interview. When EGF genotypes were compared between incident primary melanoma cases and the nonmelanoma controls, the risk associated with the homozygous G/G genotype was not statistically significantly associated with an increased risk for incident primary melanoma compared with the homozygous A/A genotype [OR, 1.09; 95% confidence interval (CI); 0.65-1.85]. No strong associations with EGF G/G genotype were observed in comparisons of multiple primary and single primary melanoma cases (OR, 0.66; 95% CI; 0.25-1.73). Case subjects with tumors >/=3.5 mm compared with those <3.5 mm were not significantly associated with the G/G genotype (OR, 0.54; 95% CI; 0.12-2.35). Our data do not support a significant association between melanoma and the EGF 61*G allele or the homozygous G/G genotype. The EGF polymorphism is not a reproducible risk factor for melanoma or thick melanoma in our data. The two analytic approaches used in the study provide evidence against a strong association between EGF 61*G and melanoma and demonstrate the potential utility of case-case designs for evaluating the role of single nucleotide polymorphisms and cancer. Additional independent studies will be required to elucidate relationships between genetic variation in the EGF gene and risk of melanoma.

Background: Mutations in BRCA1 and/or BRCA2 (BRCA1/2) profoundly increase the risks of breast and ovarian cancers, but it is unclear whether mutations in these genes increase the risk of colorectal cancer. We investigated BRCA1/2 founder mutations and a family history of breast cancer as potential risk factors for colorectal cancer., Methods: In the population-based Molecular Epidemiology of Colorectal Cancer study in northern Israel, 1422 case patients with incident colorectal cancer, diagnosed between March 31, 1998, and December 31, 2002, and 1566 control subjects without colorectal cancer were genotyped for the BRCA1 187delAG, BRCA1 5385insC, and BRCA2 6174delT founder mutations. Genotypes and interview data from all case patients and control subjects and from only those of Ashkenazi Jewish descent (1002 case patients and 1038 control subjects) were used to calculate odds ratios [ORs] from logistic regression., Results: Twenty-four (2.4%) case patients and 20 (1.9%) control subjects carried one of the three mutations (OR = 1.24, 95% confidence interval [CI] = 0.68 to 2.26). A family history of breast cancer in a female relative was not associated with an increased risk of colorectal cancer, even after adjustment for the presence of a BRCA founder mutation (OR = 1.03, 95% CI = 0.75 to 1.41)., Conclusions: Although weak associations cannot be excluded, Ashkenazi BRCA founder mutations do not confer a strongly elevated risk of colorectal cancer. Similarly, a family history of breast cancer does not appear to be a strong risk factor for colorectal cancer in this population.

Background: A mutation in the androgen receptor (AR) gene, namely AR R726L, was described in 2% of Finnish men with sporadic or familial prostate cancer and was associated with an approximately sixfold increased risk of prostate cancer. We set out to determine the incidence of this mutation in a sample of men with either early-onset and/or familial prostate cancer in the United States., Methods: Five hundred forty-eight men with prostate cancer from 411 unrelated families participating in the University of Michigan Prostate Cancer Genetics Project (PCGP) were studied. Allele-specific oligonucleotide hybridization was used to detect the presence of the AR R726L mutation in germline DNA., Results: None of the 548 prostate cancer patients studied, including 513 White, 29 African American, 3 Asian, and 3 Hispanic men, were found to carry the AR R726L allele. Therefore, the prevalence of this allele is significantly less than that observed among Finnish men with prostate cancer (Fisher's exact test, P = 0.002)., Conclusions: The AR R726L allele does not account for a significant proportion of early-onset and/or familial prostate cancer in the United States., (Copyright 2003 Wiley-Liss, Inc.)