A model-based approach for analysis of spatial structure in genetic data

Characterizing genetic diversity within and between populations has broad applications in studies of human disease and evolution. We propose a new approach, spatial ancestry analysis, for the modeling of genotypes in two- or three-dimensional space. In spatial ancestry analysis (SPA), we explicitly model the spatial distribution of each SNP by assigning an allele frequency as a continuous function in geographic space. We show that the explicit modeling of the allele frequency allows individuals to be localized on the map on the basis of their genetic information alone. We apply our SPA method to a European and a worldwide population genetic variation data set and identify SNPs showing large gradients in allele frequency, and we suggest these as candidate regions under selection. These regions include SNPs in the well-characterized LCT region, as well as at loci including FOXP2, OCA2 and LRP1B.     

A literature network of human genes for high-throughput analysis of gene expression

We have carried out automated extraction of explicit and implicit biomedical knowledge from publicly available gene and text databases to create a gene-to-gene co-citation network for 13,712 named human genes by automated analysis of titles and abstracts in over 10 million MEDLINE records. The associations between genes have been annotated by linking genes to terms from the medical subject heading (MeSH) index and terms from the gene ontology (GO) database. The extracted database and accompanying web tools for gene-expression analysis have collectively been named 'PubGene'. We validated the extracted networks by three large-scale experiments showing that co-occurrence reflects biologically meaningful relationships, thus providing an approach to extract and structure known biology. We validated the applicability of the tools by analyzing two publicly available microarray data sets.     

Genetic analysis of the mouse brain proteome

Proteome analysis is a fundamental step in systematic functional genomics. Here we have resolved 8,767 proteins from the mouse brain proteome by large-gel two-dimensional electrophoresis. We detected 1,324 polymorphic proteins from the European collaborative interspecific backcross. Of these, we mapped 665 proteins genetically and identified 466 proteins by mass spectrometry. Qualitatively polymorphic proteins, to 96%, reflect changes in conformation and/or mass. Quantitatively polymorphic proteins show a high frequency (73%) of allele-specific transmission in codominant heterozygotes. Variations in protein isoforms and protein quantity often mapped to chromosomal positions different from that of the structural gene, indicating that single proteins may act as polygenic traits. Genetic analysis of proteomes may detect the types of polymorphism that are most relevant in disease-association studies.     

A comprehensive linkage analysis for myocardial infarction and its related risk factors

Coronary artery disease and myocardial infarction (MI) are leading causes of death in the western world. Numerous studies have shown that risk factors such as diabetes mellitus, arterial hypertension and hypercholesterolemia contribute to the development of the disease. Although each risk factor by itself is partly under genetic control, a positive family history is an independent predictor, which suggests that there are additional susceptibility genes. We have scanned the whole genome in 513 families to identify chromosomal regions linked to myocardial infarction and related risk factors that are known to be under genetic control. Here we show, by using variance component analysis and incorporating risk factors, that risk of myocardial infarction maps to a single region on chromosome 14 with a significant lod score of 3.9 (pointwise P=0.00015, genome-wide P<0.05), providing evidence of a principal MI locus. To characterize this locus we analyzed each risk factor by itself. Serum concentrations of lipoprotein (a) show linkage to both the apolipoprotein (a) locus (lod score 26.99) and a new locus on chromosome 1 (lod score 3.8). There is suggestive linkage for diabetes mellitus on chromosome 6 (lod score 2.96), for hypertension on chromosomes 1 and 6, for high-density and low-density lipoprotein cholesterol on chromosomes 1 and 17, and for triglyceride concentrations on chromosome 9. Although some of these risk factors overlap with previously identified loci, none overlaps with the newly identified susceptibility locus for myocardial infarction and coronary artery disease.     

Transcriptome analysis of the acoelomate human parasite Schistosoma mansoni

Schistosoma mansoni is the primary causative agent of schistosomiasis, which affects 200 million individuals in 74 countries. We generated 163,000 expressed-sequence tags (ESTs) from normalized cDNA libraries from six selected developmental stages of the parasite, resulting in 31,000 assembled sequences and 92% sampling of an estimated 14,000 gene complement. By analyzing automated Gene Ontology assignments, we provide a detailed view of important S. mansoni biological systems, including characterization of metazoa-specific and eukarya-conserved genes. Phylogenetic analysis suggests an early divergence from other metazoa. The data set provides insights into the molecular mechanisms of tissue organization, development, signaling, sexual dimorphism, host interactions and immune evasion and identifies novel proteins to be investigated as vaccine candidates and potential drug targets.     

Large-scale analysis of the regulatory architecture of the mouse genome with a transposon-associated sensor

We present here a Sleeping Beauty-based transposition system that offers a simple and efficient way to investigate the regulatory architecture of mammalian chromosomes in vivo. With this system, we generated several hundred mice and embryos, each with a regulatory sensor inserted at a random genomic position. This large sampling of the genome revealed the widespread presence of long-range regulatory activities along chromosomes, forming overlapping blocks with distinct tissue-specific expression potentials. The presence of tissue-restricted regulatory activities around genes with widespread expression patterns challenges the gene-centric view of genome regulation and suggests that most genes are modulated in a tissue-specific manner. The local hopping property of Sleeping Beauty provides a dynamic approach to map these regulatory domains at high resolution and, combined with Cre-mediated recombination, allows for the determination of their functions by engineering mice with specific chromosomal rearrangements.     

Multiple knockout analysis of genetic robustness in the yeast metabolic network

Genetic robustness characterizes the constancy of the phenotype in face of heritable perturbations. Previous investigations have used comprehensive single and double gene knockouts to study gene essentiality and pairwise gene interactions in the yeast Saccharomyces cerevisiae. Here we conduct an in silico multiple knockout investigation of a flux balance analysis model of the yeast's metabolic network. Cataloging gene sets that provide mutual functional backup, we identify sets of up to eight interacting genes and characterize the 'k robustness' (the depth of backup interactions) of each gene. We find that 74% (360) of the metabolic genes participate in processes that are essential to growth in a standard laboratory environment, compared with only 13% previously found to be essential using single knockouts. The genes' k robustness is shown to be a solid indicator of their biological buffering capacity and is correlated with both the genes' environmental specificity and their evolutionary retention.     

Genome-wide functional analysis of pathogenicity genes in the rice blast fungus

Rapid translation of genome sequences into meaningful biological information hinges on the integration of multiple experimental and informatics methods into a cohesive platform. Despite the explosion in the number of genome sequences available, such a platform does not exist for filamentous fungi. Here we present the development and application of a functional genomics and informatics platform for a model plant pathogenic fungus, Magnaporthe oryzae. In total, we produced 21,070 mutants through large-scale insertional mutagenesis using Agrobacterium tumefaciens-mediated transformation. We used a high-throughput phenotype screening pipeline to detect disruption of seven phenotypes encompassing the fungal life cycle and identified the mutated gene and the nature of mutation for each mutant. Comparative analysis of phenotypes and genotypes of the mutants uncovered 202 new pathogenicity loci. Our findings demonstrate the effectiveness of our platform and provide new insights on the molecular basis of fungal pathogenesis. Our approach promises comprehensive functional genomics in filamentous fungi and beyond.     

zA map for sequence analysis of the Arabidopsis thaliana genome.

Arabidopsis thaliana has emerged as a model system for studies of plant genetics and development, and its genome has been targeted for sequencing by an international consortium (the Arabidopsis Genome Initiative; http://genome-www. stanford.edu/Arabidopsis/agi.html). To support the genome-sequencing effort, we fingerprinted more than 20,000 BACs (ref. 2) from two high-quality publicly available libraries, generating an estimated 17-fold redundant coverage of the genome, and used the fingerprints to nucleate assembly of the data by computer. Subsequent manual revision of the assemblies resulted in the incorporation of 19,661 fingerprinted BACs into 169 ordered sets of overlapping clones ('contigs'), each containing at least 3 clones. These contigs are ideal for parallel selection of BACs for large-scale sequencing and have supported the generation of more than 5.8 Mb of finished genome sequence submitted to GenBank; analysis of the sequence has confirmed the integrity of contigs constructed using this fingerprint data. Placement of contigs onto chromosomes can now be performed, and is being pursued by groups involved in both sequencing and positional cloning studies. To our knowledge, these data provide the first example of whole-genome random BAC fingerprint analysis of a eucaryote, and have provided a model essential to efforts aimed at generating similar databases of fingerprint contigs to support sequencing of other complex genomes, including that of human.     

A new gene involved in X-linked mental retardation identified by analysis of an X;2 balanced translocation

X-linked forms of mental retardation (MR) affect approximately 1 in 600 males and are likely to be highly heterogeneous. They can be categorized into syndromic (MRXS) and nonspecific (MRX) forms. In MRX forms, affected patients have no distinctive clinical or biochemical features. At least five MRX genes have been identified by positional cloning, but each accounts for only 0.5%-1.0% of MRX cases. Here we show that the gene TM4SF2 at Xp11.4 is inactivated by the X breakpoint of an X;2 balanced translocation in a patient with MR. Further investigation led to identification of TM4SF2 mutations in 2 of 33 other MRX families. RNA in situ hybridization showed that TM4SF2 is highly expressed in the central nervous system, including the cerebral cortex and hippocampus. TM4SF2 encodes a member of the tetraspanin family of proteins, which are known to contribute in molecular complexes including beta-1 integrins. We speculate that through this interaction, TM4SF2 might have a role in the control of neurite outgrowth.     

Genetic analysis of cavefish reveals molecular convergence in the evolution of albinism

The genetic basis of vertebrate morphological evolution has traditionally been very difficult to examine in naturally occurring populations. Here we describe the generation of a genome-wide linkage map to allow quantitative trait analysis of evolutionarily derived morphologies in the Mexican cave tetra, a species that has, in a series of independent caves, repeatedly evolved specialized characteristics adapted to a unique and well-studied ecological environment. We focused on the trait of albinism and discovered that it is linked to Oca2, a known pigmentation gene, in two cave populations. We found different deletions in Oca2 in each population and, using a cell-based assay, showed that both cause loss of function of the corresponding protein, OCA2. Thus, the two cave populations evolved albinism independently, through similar mutational events.     

Molecular analysis of the muscle pathology associated with mitochondrial DNA deletions.

Large-scale deletions of mitochondrial DNA (mtDNA) are associated with a subgroup of mitochondrial encephalomyopathies. We studied seven patients with Kearns-Sayre syndrome or isolated ocular myopathy who harboured a sub-population of partially-deleted mitochondrial genomes in skeletal muscle. Variable cytochrome c oxidase (COX) deficiencies and reduction of mitochondrially-encoded polypeptides were found in affected muscle fibres, but while many COX-deficient fibres had increased levels of mutant mtDNA, they almost invariably had reduced levels of normal mtDNA. Our results suggest that a specific ratio between mutant and wild-type mitochondrial genomes is the most important determinant of a focal respiratory chain deficiency, even though absolute copy numbers may vary widely.