Hemichordate genomes and deuterostome origins.

Acorn worms, also known as enteropneust (literally, 'gut-breathing') hemichordates, are marine invertebrates that share features with echinoderms and chordates. Together, these three phyla comprise the deuterostomes. Here we report the draft genome sequences of two acorn worms, Saccoglossus kowalevskii and Ptychodera flava. By comparing them with diverse bilaterian genomes, we identify shared traits that were probably inherited from the last common deuterostome ancestor, and then explore evolutionary trajectories leading from this ancestor to hemichordates, echinoderms and chordates. The hemichordate genomes exhibit extensive conserved synteny with amphioxus and other bilaterians, and deeply conserved non-coding sequences that are candidates for conserved gene-regulatory elements. Notably, hemichordates possess a deuterostome-specific genomic cluster of four ordered transcription factor genes, the expression of which is associated with the development of pharyngeal 'gill' slits, the foremost morphological innovation of early deuterostomes, and is probably central to their filter-feeding lifestyle. Comparative analysis reveals numerous deuterostome-specific gene novelties, including genes found in deuterostomes and marine microbes, but not other animals. The putative functions of these genes can be linked to physiological, metabolic and developmental specializations of the filter-feeding ancestor.

Mammalian Y chromosomes retain widely expressed dosage-sensitive regulators.

The human X and Y chromosomes evolved from an ordinary pair of autosomes, but millions of years ago genetic decay ravaged the Y chromosome, and only three per cent of its ancestral genes survived. We reconstructed the evolution of the Y chromosome across eight mammals to identify biases in gene content and the selective pressures that preserved the surviving ancestral genes. Our findings indicate that survival was nonrandom, and in two cases, convergent across placental and marsupial mammals. We conclude that the gene content of the Y chromosome became specialized through selection to maintain the ancestral dosage of homologous X-Y gene pairs that function as broadly expressed regulators of transcription, translation and protein stability. We propose that beyond its roles in testis determination and spermatogenesis, the Y chromosome is essential for male viability, and has unappreciated roles in Turner's syndrome and in phenotypic differences between the sexes in health and disease.

Strict evolutionary conservation followed rapid gene loss on human and rhesus Y chromosomes.

The human X and Y chromosomes evolved from an ordinary pair of autosomes during the past 200-300 million years. The human MSY (male-specific region of Y chromosome) retains only three percent of the ancestral autosomes' genes owing to genetic decay. This evolutionary decay was driven by a series of five 'stratification' events. Each event suppressed X-Y crossing over within a chromosome segment or 'stratum', incorporated that segment into the MSY and subjected its genes to the erosive forces that attend the absence of crossing over. The last of these events occurred 30 million years ago, 5 million years before the human and Old World monkey lineages diverged. Although speculation abounds regarding ongoing decay and looming extinction of the human Y chromosome, remarkably little is known about how many MSY genes were lost in the human lineage in the 25 million years that have followed its separation from the Old World monkey lineage. To investigate this question, we sequenced the MSY of the rhesus macaque, an Old World monkey, and compared it to the human MSY. We discovered that during the last 25 million years MSY gene loss in the human lineage was limited to the youngest stratum (stratum 5), which comprises three percent of the human MSY. In the older strata, which collectively comprise the bulk of the human MSY, gene loss evidently ceased more than 25 million years ago. Likewise, the rhesus MSY has not lost any older genes (from strata 1-4) during the past 25 million years, despite its major structural differences to the human MSY. The rhesus MSY is simpler, with few amplified gene families or palindromes that might enable intrachromosomal recombination and repair. We present an empirical reconstruction of human MSY evolution in which each stratum transitioned from rapid, exponential loss of ancestral genes to strict conservation through purifying selection.

The Drosophila melanogaster Genetic Reference Panel.

A major challenge of biology is understanding the relationship between molecular genetic variation and variation in quantitative traits, including fitness. This relationship determines our ability to predict phenotypes from genotypes and to understand how evolutionary forces shape variation within and between species. Previous efforts to dissect the genotype-phenotype map were based on incomplete genotypic information. Here, we describe the Drosophila melanogaster Genetic Reference Panel (DGRP), a community resource for analysis of population genomics and quantitative traits. The DGRP consists of fully sequenced inbred lines derived from a natural population. Population genomic analyses reveal reduced polymorphism in centromeric autosomal regions and the X chromosome, evidence for positive and negative selection, and rapid evolution of the X chromosome. Many variants in novel genes, most at low frequency, are associated with quantitative traits and explain a large fraction of the phenotypic variance. The DGRP facilitates genotype-phenotype mapping using the power of Drosophila genetics.

Comparative and demographic analysis of orang-utan genomes.

'Orang-utan' is derived from a Malay term meaning 'man of the forest' and aptly describes the southeast Asian great apes native to Sumatra and Borneo. The orang-utan species, Pongo abelii (Sumatran) and Pongo pygmaeus (Bornean), are the most phylogenetically distant great apes from humans, thereby providing an informative perspective on hominid evolution. Here we present a Sumatran orang-utan draft genome assembly and short read sequence data from five Sumatran and five Bornean orang-utan genomes. Our analyses reveal that, compared to other primates, the orang-utan genome has many unique features. Structural evolution of the orang-utan genome has proceeded much more slowly than other great apes, evidenced by fewer rearrangements, less segmental duplication, a lower rate of gene family turnover and surprisingly quiescent Alu repeats, which have played a major role in restructuring other primate genomes. We also describe a primate polymorphic neocentromere, found in both Pongo species, emphasizing the gradual evolution of orang-utan genome structure. Orang-utans have extremely low energy usage for a eutherian mammal, far lower than their hominid relatives. Adding their genome to the repertoire of sequenced primates illuminates new signals of positive selection in several pathways including glycolipid metabolism. From the population perspective, both Pongo species are deeply diverse; however, Sumatran individuals possess greater diversity than their Bornean counterparts, and more species-specific variation. Our estimate of Bornean/Sumatran speciation time, 400,000 years ago, is more recent than most previous studies and underscores the complexity of the orang-utan speciation process. Despite a smaller modern census population size, the Sumatran effective population size (N(e)) expanded exponentially relative to the ancestral N(e) after the split, while Bornean N(e) declined over the same period. Overall, the resources and analyses presented here offer new opportunities in evolutionary genomics, insights into hominid biology, and an extensive database of variation for conservation efforts.

Complete Khoisan and Bantu genomes from southern Africa.

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.

The complete genome of an individual by massively parallel DNA sequencing.

The association of genetic variation with disease and drug response, and improvements in nucleic acid technologies, have given great optimism for the impact of 'genomic medicine'. However, the formidable size of the diploid human genome, approximately 6 gigabases, has prevented the routine application of sequencing methods to deciphering complete individual human genomes. To realize the full potential of genomics for human health, this limitation must be overcome. Here we report the DNA sequence of a diploid genome of a single individual, James D. Watson, sequenced to 7.4-fold redundancy in two months using massively parallel sequencing in picolitre-size reaction vessels. This sequence was completed in two months at approximately one-hundredth of the cost of traditional capillary electrophoresis methods. Comparison of the sequence to the reference genome led to the identification of 3.3 million single nucleotide polymorphisms, of which 10,654 cause amino-acid substitution within the coding sequence. In addition, we accurately identified small-scale (2-40,000 base pair (bp)) insertion and deletion polymorphism as well as copy number variation resulting in the large-scale gain and loss of chromosomal segments ranging from 26,000 to 1.5 million base pairs. Overall, these results agree well with recent results of sequencing of a single individual by traditional methods. However, in addition to being faster and significantly less expensive, this sequencing technology avoids the arbitrary loss of genomic sequences inherent in random shotgun sequencing by bacterial cloning because it amplifies DNA in a cell-free system. As a result, we further demonstrate the acquisition of novel human sequence, including novel genes not previously identified by traditional genomic sequencing. This is the first genome sequenced by next-generation technologies. Therefore it is a pilot for the future challenges of 'personalized genome sequencing'.

Whole-genome sequencing in a patient with Charcot-Marie-Tooth neuropathy.

BACKGROUND: Whole-genome sequencing may revolutionize medical diagnostics through rapid identification of alleles that cause disease. However, even in cases with simple patterns of inheritance and unambiguous diagnoses, the relationship between disease phenotypes and their corresponding genetic changes can be complicated. Comprehensive diagnostic assays must therefore identify all possible DNA changes in each haplotype and determine which are responsible for the underlying disorder. The high number of rare, heterogeneous mutations present in all humans and the paucity of known functional variants in more than 90% of annotated genes make this challenge particularly difficult. Thus, the identification of the molecular basis of a genetic disease by means of whole-genome sequencing has remained elusive. We therefore aimed to assess the usefulness of human whole-genome sequencing for genetic diagnosis in a patient with Charcot-Marie-Tooth disease. METHODS: We identified a family with a recessive form of Charcot-Marie-Tooth disease for which the genetic basis had not been identified. We sequenced the whole genome of the proband, identified all potential functional variants in genes likely to be related to the disease, and genotyped these variants in the affected family members. RESULTS: We identified and validated compound, heterozygous, causative alleles in SH3TC2 (the SH3 domain and tetratricopeptide repeats 2 gene), involving two mutations, in the proband and in family members affected by Charcot-Marie-Tooth disease. Separate subclinical phenotypes segregated independently with each of the two mutations; heterozygous mutations confer susceptibility to neuropathy, including the carpal tunnel syndrome. CONCLUSIONS: As shown in this study of a family with Charcot-Marie-Tooth disease, whole-genome sequencing can identify clinically relevant variants and provide diagnostic information to inform the care of patients.

The DNA sequence, annotation and analysis of human chromosome 3.

After the completion of a draft human genome sequence, the International Human Genome Sequencing Consortium has proceeded to finish and annotate each of the 24 chromosomes comprising the human genome. Here we describe the sequencing and analysis of human chromosome 3, one of the largest human chromosomes. Chromosome 3 comprises just four contigs, one of which currently represents the longest unbroken stretch of finished DNA sequence known so far. The chromosome is remarkable in having the lowest rate of segmental duplication in the genome. It also includes a chemokine receptor gene cluster as well as numerous loci involved in multiple human cancers such as the gene encoding FHIT, which contains the most common constitutive fragile site in the genome, FRA3B. Using genomic sequence from chimpanzee and rhesus macaque, we were able to characterize the breakpoints defining a large pericentric inversion that occurred some time after the split of Homininae from Ponginae, and propose an evolutionary history of the inversion.

The finished DNA sequence of human chromosome 12.

Human chromosome 12 contains more than 1,400 coding genes and 487 loci that have been directly implicated in human disease. The q arm of chromosome 12 contains one of the largest blocks of linkage disequilibrium found in the human genome. Here we present the finished sequence of human chromosome 12, which has been finished to high quality and spans approximately 132 megabases, representing approximately 4.5% of the human genome. Alignment of the human chromosome 12 sequence across vertebrates reveals the origin of individual segments in chicken, and a unique history of rearrangement through rodent and primate lineages. The rate of base substitutions in recent evolutionary history shows an overall slowing in hominids compared with primates and rodents.

A collagen-binding adhesin, Acb, and ten other putative MSCRAMM and pilus family proteins of Streptococcus gallolyticus subsp. gallolyticus (Streptococcus bovis Group, biotype I).

Members of the Streptococcus bovis group are important causes of endocarditis. However, factors associated with their pathogenicity, such as adhesins, remain uncharacterized. We recently demonstrated that endocarditis-derived Streptococcus gallolyticus subsp. gallolyticus isolates frequently adhere to extracellular matrix (ECM) proteins. Here, we generated a draft genome sequence of an ECM protein-adherent S. gallolyticus subsp. gallolyticus strain and found, by genome-wide analyses, 11 predicted LPXTG-type cell wall-anchored proteins with characteristics of MSCRAMMs, including a modular architecture of domains predicted to adopt immunoglobulin (Ig)-like folding. A recombinant segment of one of these, Acb, showed high-affinity binding to immobilized collagen, and cell surface expression of Acb correlated with the presence of acb and collagen adherence of isolates. Three of the 11 proteins have similarities to major pilus subunits and are organized in separate clusters, each including a second Ig-fold-containing MSCRAMM and a class C sortase, suggesting that the sequenced strain encodes three distinct types of pili. Reverse transcription-PCR demonstrated that all three genes of one cluster, acb-sbs7-srtC1, are cotranscribed, consistent with pilus operons of other gram-positive bacteria. Further analysis detected expression of all 11 genes in cells grown to mid to late exponential growth phases. Wide distribution of 9 of the 11 genes was observed among S. gallolyticus subsp. gallolyticus isolates with fewer genes present in other S. bovis group species/subspecies. The high prevalence of genes encoding putative MSCRAMMs and pili, including a collagen-binding MSCRAMM, among S. gallolyticus subsp. gallolyticus isolates may play an important role in the predominance of this subspecies in S. bovis endocarditis.

Two-stage genome-wide association study identifies integrin beta 5 as having potential role in bull fertility.

BACKGROUND: Fertility is one of the most critical factors controlling biological and financial performance of animal production systems and genetic improvement of lines. The objective of this study was to identify molecular defects in the sperm that are responsible for uncompensable fertility in Holstein bulls. We performed a comprehensive genome wide analysis of single nucleotide polymorphisms (SNP) for bull fertility followed by a second-stage replication in additional bulls for a restricted set of markers. RESULTS: In the Phase I association study, we genotyped the genomic sperm DNA of 10 low-fertility and 10 high-fertility bulls using Bovine SNP Gene Chips containing approximately 10,000 random SNP markers. In these animals, 8,207 markers were found to be polymorphic, 97 of which were significantly associated with fertility (p < 0.01). In the Phase II study, we tested the four most significant SNP from the Phase I study in 101 low-fertility and 100 high-fertility bulls, with two SNPs (rs29024867 and rs41257187) significantly replicated. Rs29024867 corresponds to a nucleotide change of C --> G 2,190 bp 3' of the collagen type I alpha 2 gene on chromosome 4, while the rs41257187 (C --> T) is in the coding region of integrin beta 5 gene on chromosome 1. The SNP rs41257187 induces a synonymous (Proline --> Proline), suggesting disequilibrium with the true causative locus (i), but we found that the incubation of bull spermatozoa with integrin beta 5 antibodies significantly decreased the ability to fertilize oocytes. Our findings suggest that the bovine sperm integrin beta 5 protein plays a role during fertilization and could serve as a positional or functional marker of bull fertility. CONCLUSION: We have identified molecular markers associated with bull fertility and established that at least one of the genes harboring such variation has a role in fertility. The findings are important in understanding mechanisms of uncompensatory infertility in bulls, and in other male mammals. The findings set the stage for more hypothesis-driven research aimed at discovering the role of variation in the genome that affect fertility and that can be used to identify molecular mechanisms of development.

Subtle genetic changes enhance virulence of methicillin resistant and sensitive Staphylococcus aureus.

BACKGROUND: Community acquired (CA) methicillin-resistant Staphylococcus aureus (MRSA) increasingly causes disease worldwide. USA300 has emerged as the predominant clone causing superficial and invasive infections in children and adults in the USA. Epidemiological studies suggest that USA300 is more virulent than other CA-MRSA. The genetic determinants that render virulence and dominance to USA300 remain unclear. RESULTS: We sequenced the genomes of two pediatric USA300 isolates: one CA-MRSA and one CA-methicillin susceptible (MSSA), isolated at Texas Children's Hospital in Houston. DNA sequencing was performed by Sanger dideoxy whole genome shotgun (WGS) and 454 Life Sciences pyrosequencing strategies. The sequence of the USA300 MRSA strain was rigorously annotated. In USA300-MRSA 2658 chromosomal open reading frames were predicted and 3.1 and 27 kilobase (kb) plasmids were identified. USA300-MSSA contained a 20 kb plasmid with some homology to the 27 kb plasmid found in USA300-MRSA. Two regions found in US300-MRSA were absent in USA300-MSSA. One of these carried the arginine deiminase operon that appears to have been acquired from S. epidermidis. The USA300 sequence was aligned with other sequenced S. aureus genomes and regions unique to USA300 MRSA were identified. CONCLUSION: USA300-MRSA is highly similar to other MRSA strains based on whole genome alignments and gene content, indicating that the differences in pathogenesis are due to subtle changes rather than to large-scale acquisition of virulence factor genes. The USA300 Houston isolate differs from another sequenced USA300 strain isolate, derived from a patient in San Francisco, in plasmid content and a number of sequence polymorphisms. Such differences will provide new insights into the evolution of pathogens.

A massive expansion of effector genes underlies gall-formation in the wheat pest Mayetiola destructor.

Gall-forming arthropods are highly specialized herbivores that, in combination with their hosts, produce extended phenotypes with unique morphologies [1]. Many are economically important, and others have improved our understanding of ecology and adaptive radiation [2]. However, the mechanisms that these arthropods use to induce plant galls are poorly understood. We sequenced the genome of the Hessian fly (Mayetiola destructor; Diptera: Cecidomyiidae), a plant parasitic gall midge and a pest of wheat (Triticum spp.), with the aim of identifying genic modifications that contribute to its plant-parasitic lifestyle. Among several adaptive modifications, we discovered an expansive reservoir of potential effector proteins. Nearly 5% of the 20,163 predicted gene models matched putative effector gene transcripts present in the M. destructor larval salivary gland. Another 466 putative effectors were discovered among the genes that have no sequence similarities in other organisms. The largest known arthropod gene family (family SSGP-71) was also discovered within the effector reservoir. SSGP-71 proteins lack sequence homologies to other proteins, but their structures resemble both ubiquitin E3 ligases in plants and E3-ligase-mimicking effectors in plant pathogenic bacteria. SSGP-71 proteins and wheat Skp proteins interact in vivo. Mutations in different SSGP-71 genes avoid the effector-triggered immunity that is directed by the wheat resistance genes H6 and H9. Results point to effectors as the agents responsible for arthropod-induced plant gall formation.

Sequencing the full-length of the phosphatase and tensin homolog (PTEN) gene in hepatocellular carcinoma (HCC) using the 454 GS20 and Illumina GA DNA sequencing platforms.

BACKGROUND: Phosphatase and tensin homolog (PTEN) is a tumor-suppressor gene that is mutated in cancer of the liver, pancreas, endometrium, and prostate. PTEN-dependent pathways are involved in mediating cell growth and invasion. To sequence the whole gene (including introns and exons), we have taken advantage of new technologies that allow for rapid, inexpensive sequencing to great depth. METHODS: DNA from 15 HCC specimens were pooled, and long-range PCR was performed by using the GeneAmp XL PCR kit. Primer parameters included: length of 20-30 base pairs (bp), melting temperature of -68 degrees C, and G/C content of 50-60%. PCR products were then column-purified and pooled, and DNA libraries were prepared for "shotgun sequencing" on both the 454 GS and Illumina GA sequencing platforms. RESULTS: We successfully amplified approximately 98.9% of the PTEN gene by using one long-range PCR protocol applied to 24 primer sets, resulting in 20 amplicons approximately 6.5 kilobases (kb) in length, 2 amplicons approximately 10 kb in length, and 2 amplicons approximately 2.5 kb in length. Sequencing of fragmented PCR products on both sequencing platforms identified six high-frequency SNPs that were catalogued in dbSNP as known variants. CONCLUSIONS: Shotgun sequencing based on a single long-range PCR protocol in pooled samples is an efficient and relatively inexpensive way to sequence an entire gene.

Genome-wide survey of SNP variation uncovers the genetic structure of cattle breeds.

The imprints of domestication and breed development on the genomes of livestock likely differ from those of companion animals. A deep draft sequence assembly of shotgun reads from a single Hereford female and comparative sequences sampled from six additional breeds were used to develop probes to interrogate 37,470 single-nucleotide polymorphisms (SNPs) in 497 cattle from 19 geographically and biologically diverse breeds. These data show that cattle have undergone a rapid recent decrease in effective population size from a very large ancestral population, possibly due to bottlenecks associated with domestication, selection, and breed formation. Domestication and artificial selection appear to have left detectable signatures of selection within the cattle genome, yet the current levels of diversity within breeds are at least as great as exists within humans.

Demographic histories and patterns of linkage disequilibrium in Chinese and Indian rhesus macaques.

To understand the demographic history of rhesus macaques (Macaca mulatta) and document the extent of linkage disequilibrium (LD) in the genome, we partially resequenced five Encyclopedia of DNA Elements regions in 9 Chinese and 38 captive-born Indian rhesus macaques. Population genetic analyses of the 1467 single-nucleotide polymorphisms discovered suggest that the two populations separated about 162,000 years ago, with the Chinese population tripling in size since then and the Indian population eventually shrinking by a factor of four. Using coalescent simulations, we confirmed that these inferred demographic events explain a much faster decay of LD in Chinese (r(2) approximately 0.15 at 10 kilobases) versus Indian (r(2) approximately 0.52 at 10 kilobases) macaque populations.

Direct selection of human genomic loci by microarray hybridization.

We applied high-density microarrays to the enrichment of specific sequences from the human genome for high-throughput sequencing. After capture of 6,726 approximately 500-base 'exon' segments, and of 'locus-specific' regions ranging in size from 200 kb to 5 Mb, followed by sequencing on a 454 Life Sciences FLX sequencer, most sequence reads represented selection targets. These direct selection methods supersede multiplex PCR for the large-scale analysis of genomic regions.

Paradoxical DNA repair and peroxide resistance gene conservation in Bacillus pumilus SAFR-032.

BACKGROUND: Bacillus spores are notoriously resistant to unfavorable conditions such as UV radiation, gamma-radiation, H2O2, desiccation, chemical disinfection, or starvation. Bacillus pumilus SAFR-032 survives standard decontamination procedures of the Jet Propulsion Lab spacecraft assembly facility, and both spores and vegetative cells of this strain exhibit elevated resistance to UV radiation and H2O2 compared to other Bacillus species. PRINCIPAL FINDINGS: The genome of B. pumilus SAFR-032 was sequenced and annotated. Lists of genes relevant to DNA repair and the oxidative stress response were generated and compared to B. subtilis and B. licheniformis. Differences in conservation of genes, gene order, and protein sequences are highlighted because they potentially explain the extreme resistance phenotype of B. pumilus. The B. pumilus genome includes genes not found in B. subtilis or B. licheniformis and conserved genes with sequence divergence, but paradoxically lacks several genes that function in UV or H2O2 resistance in other Bacillus species. SIGNIFICANCE: This study identifies several candidate genes for further research into UV and H2O2 resistance. These findings will help explain the resistance of B. pumilus and are applicable to understanding sterilization survival strategies of microbes.