Genome-wide in situ exon capture for selective resequencing

Increasingly powerful sequencing technologies are ushering in an era of personal genome sequences and raising the possibility of using such information to guide medical decisions. Genome resequencing also promises to accelerate the identification of disease-associated mutations. Roughly 98% of the human genome is composed of repeats and intergenic or non-protein-coding sequences. Thus, it is crucial to focus resequencing on high-value genomic regions. Protein-coding exons represent one such type of high-value target. We have developed a method of using flexible, high-density microarrays to capture any desired fraction of the human genome, in this case corresponding to more than 200,000 protein-coding exons. Depending on the precise protocol, up to 55-85% of the captured fragments are associated with targeted regions and up to 98% of intended exons can be recovered. This methodology provides an adaptable route toward rapid and efficient resequencing of any sizeable, non-repeat portion of the human genome.     

Extensive genomic duplication during early chordate evolution

Opinions on the hypothesis that ancient genome duplications contributed to the vertebrate genome range from strong skepticism to strong credence. Previous studies concentrated on small numbers of gene families or chromosomal regions that might not have been representative of the whole genome, or used subjective methods to identify paralogous genes and regions. Here we report a systematic and objective analysis of the draft human genome sequence to identify paralogous chromosomal regions (paralogons) formed during chordate evolution and to estimate the ages of duplicate genes. We found that the human genome contains many more paralogons than would be expected by chance. Molecular clock analysis of all protein families in humans that have orthologs in the fly and nematode indicated that a burst of gene duplication activity took place in the period 350 650 Myr ago and that many of the duplicate genes formed at this time are located within paralogons. Our results support the contention that many of the gene families in vertebrates were formed or expanded by large-scale DNA duplications in an early chordate. Considering the incompleteness of the sequence data and the antiquity of the event, the results are compatible with at least one round of polyploidy.     

Systematic generation of high-resolution deletion coverage of the Drosophila melanogaster genome

In fruit fly research, chromosomal deletions are indispensable tools for mapping mutations, characterizing alleles and identifying interacting loci. Most widely used deletions were generated by irradiation or chemical mutagenesis. These methods are labor-intensive, generate random breakpoints and result in unwanted secondary mutations that can confound phenotypic analyses. Most of the existing deletions are large, have molecularly undefined endpoints and are maintained in genetically complex stocks. Furthermore, the existence of haplolethal or haplosterile loci makes the recovery of deletions of certain regions exceedingly difficult by traditional methods, resulting in gaps in coverage. Here we describe two methods that address these problems by providing for the systematic isolation of targeted deletions in the D. melanogaster genome. The first strategy used a P element-based technique to generate deletions that closely flank haploinsufficient genes and minimize undeleted regions. This deletion set has increased overall genomic coverage by 5-7%. The second strategy used FLP recombinase and the large array of FRT-bearing insertions described in the accompanying paper to generate 519 isogenic deletions with molecularly defined endpoints. This second deletion collection provides 56% genome coverage so far. The latter methodology enables the generation of small custom deletions with predictable endpoints throughout the genome and should make their isolation a simple and routine task.     

SGS1, the Saccharomyces cerevisiae homologue of BLM and WRN, suppresses genome instability and homeologous recombination

The Escherichia coli gene recQ was identified as a RecF recombination pathway gene. The gene SGS1, encoding the only RecQ-like DNA helicase in Saccharomyces cerevisiae, was identified by mutations that suppress the top3 slow-growth phenotype. Relatively little is known about the function of Sgs1p because single mutations in SGS1 do not generally cause strong phenotypes. Mutations in genes encoding RecQ-like DNA helicases such as the Bloom and Werner syndrome genes, BLM and WRN, have been suggested to cause increased genome instability. But the exact DNA metabolic defect that might underlie such genome instability has remained unclear. To better understand the cellular role of the RecQ-like DNA helicases, sgs1 mutations were analyzed for their effect on genome rearrangements. Mutations in SGS1 increased the rate of accumulating gross chromosomal rearrangements (GCRs), including translocations and deletions containing extended regions of imperfect homology at their breakpoints. sgs1 mutations also increased the rate of recombination between DNA sequences that had 91% sequence homology. Epistasis analysis showed that Sgs1p is redundant with DNA mismatch repair (MMR) for suppressing GCRs and for suppressing recombination between divergent DNA sequences. This suggests that defects in the suppression of rearrangements involving divergent, repeated sequences may underlie the genome instability seen in BLM and WRN patients and in cancer cases associated with defects in these genes.     

Estimate of human gene number provided by genome-wide analysis using Tetraodon nigroviridis DNA sequence

The number of genes in the human genome is unknown, with estimates ranging from 50,000 to 90,000 (refs 1, 2), and to more than 140,000 according to unpublished sources. We have developed 'Exofish', a procedure based on homology searches, to identify human genes quickly and reliably. This method relies on the sequence of another vertebrate, the pufferfish Tetraodon nigroviridis, to detect conserved sequences with a very low background. Similar to Fugu rubripes, a marine pufferfish proposed by Brenner et al. as a model for genomic studies, T. nigroviridis is a more practical alternative with a genome also eight times more compact than that of human. Many comparisons have been made between F. rubripes and human DNA that demonstrate the potential of comparative genomics using the pufferfish genome. Application of Exofish to the December version of the working draft sequence of the human genome and to Unigene showed that the human genome contains 28,000-34,000 genes, and that Unigene contains less than 40% of the protein-coding fraction of the human genome.     

Common deletion polymorphisms in the human genome

The locations and properties of common deletion variants in the human genome are largely unknown. We describe a systematic method for using dense SNP genotype data to discover deletions and its application to data from the International HapMap Consortium to characterize and catalogue segregating deletion variants across the human genome. We identified 541 deletion variants (94% novel) ranging from 1 kb to 745 kb in size; 278 of these variants were observed in multiple, unrelated individuals, 120 in the homozygous state. The coding exons of ten expressed genes were found to be commonly deleted, including multiple genes with roles in sex steroid metabolism, olfaction and drug response. These common deletion polymorphisms typically represent ancestral mutations that are in linkage disequilibrium with nearby SNPs, meaning that their association to disease can often be evaluated in the course of SNP-based whole-genome association studies.     

Frequent somatic mutations in MAP3K5 and MAP3K9 in metastatic melanoma identified by exome sequencing

We sequenced eight melanoma exomes to identify new somatic mutations in metastatic melanoma. Focusing on the mitogen-activated protein (MAP) kinase kinase kinase (MAP3K) family, we found that 24% of melanoma cell lines have mutations in the protein-coding regions of either MAP3K5 or MAP3K9. Structural modeling predicted that mutations in the kinase domain may affect the activity and regulation of these protein kinases. The position of the mutations and the loss of heterozygosity of MAP3K5 and MAP3K9 in 85% and 67% of melanoma samples, respectively, together suggest that the mutations are likely to be inactivating. In in vitro kinase assays, MAP3K5 I780F and MAP3K9 W333* variants had reduced kinase activity. Overexpression of MAP3K5 or MAP3K9 mutants in HEK293T cells reduced the phosphorylation of downstream MAP kinases. Attenuation of MAP3K9 function in melanoma cells using siRNA led to increased cell viability after temozolomide treatment, suggesting that decreased MAP3K pathway activity can lead to chemoresistance in melanoma.