Genetic linkage mapping of zebrafish genes and ESTs

Genetic screens in zebrafish (Danio rerio) have isolated mutations in hundreds of genes essential for vertebrate development, physiology, and behavior. We have constructed a genetic linkage map that will facilitate the identification of candidate genes for these mutations and allow comparisons among the genomes of zebrafish and other vertebrates. On this map, we have localized 771 zebrafish genes and expressed sequence tags (ESTs) by scoring single-stranded conformational polymorphisms (SSCPs) in a meiotic mapping panel. Of these sequences, 642 represent previously unmapped genes and ESTs. The mapping panel was comprised of 42 homozygous diploid individuals produced by heat shock treatment of haploid embryos at the one-cell stage (HS diploids). This "doubled haploid" strategy combines the advantages of mapping in haploid and standard diploid systems, because heat shock diploid individuals have only one allele at each locus and can survive to adulthood, enabling a relatively large quantity of genomic DNA to be prepared from each individual in the mapping panel. To integrate this map with others, we also scored 593 previously mapped simple-sequence length polymorphisms (SSLPs) in the mapping panel. This map will accelerate the molecular analysis of zebrafish mutations and facilitate comparative analysis of vertebrate genomes.     

A comparative map of the zebrafish genome

Zebrafish mutations define the functions of hundreds of essential genes in the vertebrate genome. To accelerate the molecular analysis of zebrafish mutations and to facilitate comparisons among the genomes of zebrafish and other vertebrates, we used a homozygous diploid meiotic mapping panel to localize polymorphisms in 691 previously unmapped genes and expressed sequence tags (ESTs). Together with earlier efforts, this work raises the total number of markers scored in the mapping panel to 2119, including 1503 genes and ESTs and 616 previously characterized simple-sequence length polymorphisms. Sequence analysis of zebrafish genes mapped in this study and in prior work identified putative human orthologs for 804 zebrafish genes and ESTs. Map comparisons revealed 139 new conserved syntenies, in which two or more genes are on the same chromosome in zebrafish and human. Although some conserved syntenies are quite large, there were changes in gene order within conserved groups, apparently reflecting the relatively frequent occurrence of inversions and other intrachromosomal rearrangements since the divergence of teleost and tetrapod ancestors. Comparative mapping also shows that there is not a one-to-one correspondence between zebrafish and human chromosomes. Mapping of duplicate gene pairs identified segments of 20 linkage groups that may have arisen during a genome duplication that occurred early in the evolution of teleosts after the divergence of teleost and mammalian ancestors. This comparative map will accelerate the molecular analysis of zebrafish mutations and enhance the understanding of the evolution of the vertebrate genome.     

The zebrafish gene map defines ancestral vertebrate chromosomes

Genetic screens in zebrafish (Danio rerio) have identified mutations that define the roles of hundreds of essential vertebrate genes. Genetic maps can link mutant phenotype with gene sequence by providing candidate genes for mutations and polymorphic genetic markers useful in positional cloning projects. Here we report a zebrafish genetic map comprising 4073 polymorphic markers, with more than twice the number of coding sequences localized in previously reported zebrafish genetic maps. We use this map in comparative studies to identify numerous regions of synteny conserved among the genomes of zebrafish, Tetraodon, and human. In addition, we use our map to analyze gene duplication in the zebrafish and Tetraodon genomes. Current evidence suggests that a whole-genome duplication occurred in the teleost lineage after it split from the tetrapod lineage, and that only a subset of the duplicates have been retained in modern teleost genomes. It has been proposed that differential retention of duplicate genes may have facilitated the isolation of nascent species formed during the vast radiation of teleosts. We find that different duplicated genes have been retained in zebrafish and Tetraodon, although similar numbers of duplicates remain in both genomes. Finally, we use comparative mapping data to address the proposal that the common ancestor of vertebrates had a genome consisting of 12 chromosomes. In a three-way comparison between the genomes of zebrafish, Tetraodon, and human, our analysis delineates the gene content for 11 of these 12 proposed ancestral chromosomes.     

A genetic linkage map for zebrafish: comparative analysis and localization of genes and expressed sequences

Genetic screens in zebrafish (Danio rerio) have isolated mutations in hundreds of genes with essential functions. To facilitate the identification of candidate genes for these mutations, we have genetically mapped 104 genes and expressed sequence tags by scoring single-strand conformational polymorphisms in a panel of haploid siblings. To integrate this map with existing genetic maps, we also scored 275 previously mapped genes, microsatellites, and sequence-tagged sites in the same haploid panel. Systematic phylogenetic analysis defined likely mammalian orthologs of mapped zebrafish genes, and comparison of map positions in zebrafish and mammals identified significant conservation of synteny. This comparative analysis also identified pairs of zebrafish genes that appear to be orthologous to single mammalian genes, suggesting that these genes arose in a genome duplication that occurred in the teleost lineage after the divergence of fish and mammal ancestors. This comparative map analysis will be useful in predicting the locations of zebrafish genes from mammalian gene maps and in understanding the evolution of the vertebrate genome.     

High-throughput gene mapping in Caenorhabditis elegans

Positional cloning of mutations in model genetic systems is a powerful method for the identification of targets of medical and agricultural importance. To facilitate the high-throughput mapping of mutations in Caenorhabditis elegans, we have identified a further 9602 putative new single nucleotide polymorphisms (SNPs) between two C. elegans strains, Bristol N2 and the Hawaiian mapping strain CB4856, by sequencing inserts from a CB4856 genomic DNA library and using an informatics pipeline to compare sequences with the canonical N2 genomic sequence. When combined with data from other laboratories, our marker set of 17,189 SNPs provides even coverage of the complete worm genome. To date, we have confirmed >1099 evenly spaced SNPs (one every 91 +/- 56 kb) across the six chromosomes and validated the utility of our SNP marker set and new fluorescence polarization-based genotyping methods for systematic and high-throughput identification of genes in C. elegans by cloning several proprietary genes. We illustrate our approach by recombination mapping and confirmation of the mutation in the cloned gene, dpy-18.     

Single nucleotide polymorphisms detection based on DNA microarray technology: HLA as a model

The significance of DNA variations among individuals, including single nucleotide polymorphisms (SNPs) and/or genome nucleotide mutations as well as to their detection by using new technology, will improve and facilitate the knowledge of each gene sequence. Microarray may provide information about thousands of gene simultaneously, leading to a more rapid and accurate genotyping. In this view, we developed a new methodology as an example for the detection of SNPs based on DNA microarray, using a panel of HLA alleles representative of loci A, B, DRB1. A panel of 180 oligonucleotide probes was selected to identify polymorphic positions located in exons 2 and 3 of HLA-A and B, and in exon 2 of HLA-DRB1 locus. Each oligonucleotide sequence was designed with a nucleotide mismatch located in the same position as the center of the hybridization sequence. Hybridization experiments were carried out with genomic probes constructed with an asymmetric PCR strategy. The amplified DNAs were obtained from bone marrow cells of donors previously typed for transplant. The results obtained were showing that the method was reliable thus providing a feasible technique both for HLA typing and for the investigation of other regions of genetic and clinical interest including polymorphisms correlated with different autoimmune diseases.     

DNA analysis by fluorescence quenching detection

The analysis of human genetic variations such as single nucleotide polymorphisms (SNPs) has great applications in genome-wide association studies of complex genetic traits. We have developed an SNP genotyping method based on the primer extension assay with fluorescence quenching as the detection. The template-directed dye-terminator incorporation with fluorescence quenching detection (FQ-TDI) assay is based on the observation that the intensity of fluorescent dye R110- and R6G-labeled acycloterminators is universally quenched once they are incorporated onto a DNA oligonucleotide primer. By comparing the rate of fluorescence quenching of the two allelic dyes in real time, we have extended this method for allele frequency estimation of SNPs in pooled DNA samples. The kinetic FQ-TDI assay is highly accurate and reproducible both in genotyping and in allele frequency estimation. Allele frequencies estimated by the kinetic FQ-TDI assay correlated well with known allele frequencies, with an r(2) value of 0.993. Applying this strategy to large-scale studies will greatly reduce the time and cost for genotyping hundreds and thousands of SNP markers between affected and control populations.     

Microarray-based assay for the detection of genetic variations of structural genes of West Nile virus

Adaptation through fixation of spontaneous mutations in the viral genome is considered to be one of the important factors that enable recurrent West Nile virus (WNV) outbreaks in the U.S. Genetic variations can alter viral phenotype and virulence, and degrade the performance of diagnostic and screening assays, vaccines, and potential therapeutic agents. A microarray assay was developed and optimized for the simultaneous detection of any nucleotide mutations in the entire structural region of WNV in order to facilitate public health surveillance of genetic variation of WNV. The DNA microarray consists of 263 oligonucleotide probes overlapping at half of their lengths which have been immobilized on an amine-binding glass slide. The assay was validated using 23 WNV isolates from the 2002-2005 U.S. epidemics. Oligonucleotide-based WNV arrays detected unambiguously all mutations in the structural region of each one of the isolates identified previously by sequencing analysis, serving as a rapid and effective approach for the identification of mutations in the WNV genome.     

Mapping common regulatory variants to human haplotypes

Inter-individual variation in gene expression has proven to be in part governed by genetic determinants, which may be trans- or cis-acting. The underlying cause of cis-acting regulatory variation has been identified in only a handful of the hundreds of genes shown to display differential allelic expression. In this report, we describe a systematic effort to map common cis-acting variants in 64 genes, using association methods in HapMap samples. We identified 16 loci (25%), each of which harbors common haplotypes that affect total expression of a gene, and a further 17 loci (27%) with evidence of haplotypes affecting relative allelic expression in heterozygote samples. Our survey suggests that detailed mapping of allele-specific in vivo expression will provide a rich source of regulatory SNPs or haplotypes that should be given high priority in association studies of human phenotypes.     

Novel genomic techniques open new avenues in the analysis of monogenic disorders

The molecular genetic cause of over 3,000 monogenic disorders is currently unknown. This review discusses how novel genomic techniques like Next‐Generation DNA Sequencing (NGS) and genotyping arrays open new avenues in the elucidation of genetic defects causing monogenic disorders. They will not only speed up disease gene identification but will enable us to systematically tackle previously intractable monogenic disorders. These are mainly disorders not amenable to classic linkage analysis, for example, due to insufficient family size. Most monogenic diseases are caused by exonic mutations or splice‐site mutations changing the amino acid sequence of the affected gene. These mutations can be identified by sequencing of all exons in the human genome (exome sequencing) rendering whole genome sequencing unnecessary in most cases. Genotyping arrays containing 10⁵–2×10⁶ single nucleotide polymorphisms (SNPs) and nonpolymorphic markers allow highly accurate mapping of genomic deletions and duplications not detectable by exome sequencing, which are the second most common cause of monogenic disorders. However, several hundred rare, previously unknown sequence variants affecting the amino acid sequence of the encoded protein are found in the exome of every human individual. Therefore, the main challenge will be the differentiation between the many rare benign variants detected by novel genomic techniques and disease causing mutations. Hum Mutat 32:144–151, 2011. © 2011 Wiley‐Liss, Inc.     

15000 unique zebrafish EST clusters and their future use in microarray for profiling gene expression patterns during embryogenesis

A total of 15590 unique zebrafish EST clusters from two cDNA libraries have been identified. Most significantly, only 22% (3437) of the 15590 unique clusters matched 2805 (of 15200) clusters in the Danio rerio UniGene database, indicating that our EST set is complementary to the existing ESTs in the public database and will be invaluable in assisting the annotation of genes based on the upcoming zebrafish genome sequence. Blast search showed that 7824 of our unique clusters matched 6710 known or predicted proteins in the nonredundant database. A cDNA microarray representing approximately 3100 unique zebrafish cDNA clusters has been generated and used to profile the gene expression patterns across six different embryonic stages (cleavage, blastula, gastrula, segmentation, pharyngula, and hatching). Analysis of expression data using K-means clustering revealed that genes coding for muscle-specific proteins displayed similar expression patterns, confirming that the coordinate gene expression is important for myogenesis. Our results demonstrate that the combination of microarray technology with the zebrafish model system can provide useful information on how genes are coordinated in a genetic network to control zebrafish embryogenesis and can help to identify novel genes that are important for organogenesis.     

Sequence variation in the human T-cell receptor loci

Summary: Identifying common sequence variations known as single nucleotide polymorphisms (SNPs) in human populations is one of the current objectives of the human genome project. Nearly 3 million SNPs have been identified. Analysis of the relative allele frequency of these markers in human populations and the genetic associations between these markers, known as linkage disequilibrium, is now underway to generate a high‐density genetic map. Because of the central role T cells play in immune reactivity, the T‐cell receptor (TCR) loci have long been considered important candidates for common disease susceptibility within the immune system (e.g., asthma, atopy and autoimmunity). Over the past two decades, hundreds of SNPs in the TCR loci have been identified. Most studies have focused on defining SNPs in the variable gene segments which are involved in antigenic recognition. On average, the coding sequence of each TCR variable gene segment contains two SNPs, with many more found in the 5′, 3′ and intronic sequences of these segments. Therefore, a potentially large repertoire of functional variants exists in these loci. Association between SNPs (linkage disequilibrium) extends approximately 30 kb in the TCR loci, although a few larger regions of disequilibrium have been identified. Therefore, the SNPs found in one variable gene segment may or may not be associated with SNPs in other surrounding variable gene segments. This suggests that meaningful association studies in the TCR loci will require the analysis and typing of large marker sets to fully evaluate the role of TCR loci in common disease susceptibility in human populations.     

Detection of single nucleotide substitution by competitive allele-specific short oligonucleotide hybridization (CASSOH) with immunochromatographic strip

Recent advances in human genome research have revealed that genetic polymorphisms, such as single nucleotide polymorphisms (SNPs), are closely associated with susceptibility to various common diseases and adverse drug reactions. Also, numerous mutations responsible for a number of genetic diseases have been identified. Clinical application of genetic information to individual health care requires simple and rapid identification of nucleotide changes in clinical settings. We have devised a novel low-tech method for the detection of a single nucleotide substitution using competitive allele-specific short oligonucleotide hybridization with immunochromatographic strip. The gene of interest is PCR-amplified, hybridized to an allele-specific short oligonucleotide probe in the presence of a competitive oligonucleotide, and subjected to chromatography using a DNA test strip at room temperature. The genotype is unambiguously determined by the presence or the absence of visible purple lines on a strip. Feasibility of the method was demonstrated by the detection of a prevalent disease-causing mutations in glycogen storage disease type Ia (G6PC), medium-chain acyl-CoA dehydrogenase deficiency (ACADM), non-ketotic hyperglycinemia (GLDC), and clinically important polymorphisms in the CYP2C19 gene and the aldehyde dehydrogenase 2 gene (ALDH2). The procedure does not demand either technical expertise or expensive instruments and is readily performed in local clinical laboratories. The result is obtained within 10 min after PCR. This rapid and simple method of SNP detection may be used for point-of-care genetic diagnosis with potentially diverse clinical applications. Hum Mutat 22:166-172, 2003.     

Spectrum of sequence variations in the FANCA gene: an International Fanconi Anemia Registry (IFAR) study

Fanconi anemia (FA) is an autosomal recessive disorder that is defined by cellular hypersensitivity to DNA cross-linking agents, and is characterized clinically by developmental abnormalities, progressive bone-marrow failure, and predisposition to leukemia and solid tumors. There is extensive genetic heterogeneity, with at least 11 different FA complementation groups. FA-A is the most common group, accounting for approximately 65% of all affected individuals. The mutation spectrum of the FANCA gene, located on chromosome 16q24.3, is highly heterogeneous. Here we summarize all sequence variations (mutations and polymorphisms) in FANCA described in the literature and listed in the Fanconi Anemia Mutation Database as of March 2004, and report 61 novel FANCA mutations identified in FA patients registered in the International Fanconi Anemia Registry (IFAR). Thirty-eight novel SNPs, previously unreported in the literature or in dbSNP, were also identified. We studied the segregation of common FANCA SNPs in FA families to generate haplotypes. We found that FANCA SNP data are highly useful for carrier testing, prenatal diagnosis, and preimplantation genetic diagnosis, particularly when the disease-causing mutations are unknown. Twenty-two large genomic deletions were identified by detection of apparent homozygosity for rare SNPs. In addition, a conserved SNP haplotype block spanning at least 60 kb of the FANCA gene was identified in individuals from various ethnic groups.     

Parallel genotyping of human SNPs using generic high-density oligonucleotide tag arrays

Large scale human genetic studies require technologies for generating millions of genotypes with relative ease but also at a reasonable cost and with high accuracy. We describe a highly parallel method for genotyping single nucleotide polymorphisms (SNPs), using generic high-density oligonucleotide arrays that contain thousands of preselected 20-mer oligonucleotide tags. First, marker-specific primers are used in PCR amplifications of genomic regions containing SNPs. Second, the amplification products are used as templates in single base extension (SBE) reactions using chimeric primers with 3' complementarity to the specific SNP loci and 5' complementarity to specific probes, or tags, synthesized on the array. The SBE primers, terminating one base before the polymorphic site, are extended in the presence of labeled dideoxy NTPs, using a different label for each of the two SNP alleles, and hybridized to the tag array. Third, genotypes are deduced from the fluorescence intensity ratio of the two colors. This approach takes advantage of multiplexed sample preparation, hybridization, and analysis at each stage. We illustrate and test this method by genotyping 44 individuals for 142 human SNPs identified previously in 62 candidate hypertension genes. Because the hybridization results are quantitative, this method can also be used for allele-frequency estimation in pooled DNA samples.     

Novel ENU-induced eye mutations in the mouse: models for human eye disease

We have carried out a genome-wide screen for novel N-ethyl-N-nitrosourea-induced mutations that give rise to eye and vision abnormalities in the mouse and have identified 25 inherited phenotypes that affect all parts of the eye. A combination of genetic mapping, complementation and molecular analysis revealed that 14 of these are mutations in genes previously identified to play a role in eye pathophysiology, namely Pax6, Mitf, Egfr and Pde6b. Many of the others are located in genomic regions lacking candidate genes and these define new loci. Four of the mutants display a similar phenotype of dilated pupils but do not appear to be allelic, and at least two of these are embryonic lethal when homozygous. This collection of eye mutations will be valuable for understanding gene function, for dissecting protein function and as models of human eye disease.     

OBSL1 mutations in 3‐M syndrome are associated with a modulation of IGFBP2 and IGFBP5 expression levels

3‐M syndrome is an autosomal recessive disorder characterized by severe pre‐ and postnatal growth retardation and minor skeletal changes. We have previously identified CUL7 as a disease‐causing gene but we have also provided evidence of genetic heterogeneity in the 3‐M syndrome. By homozygosity mapping in two inbred families, we found a second disease locus on chromosome 2q35–36.1 in a 5.2‐Mb interval that encompasses 60 genes. To select candidate genes, we performed microarray analysis of cultured skin fibroblast RNA from one patient, looking for genes with altered expression; we found decreased expression of IGFBP2 and increased expression of IGFBP5. However, direct sequencing of these two genes failed to detect any anomaly. We then considered other candidate genes by their function/location and found nine distinct mutations in the OBSL1 gene in 13 families including eight nonsense and one missense mutations. To further understand the links between OBSL1, CUL7, and insulin‐like growth factor binding proteins (IGFBPs), we performed real‐time quantitative PCR (RT‐PCR) analysis for OBSL1, CUL7, IGFBP2, and IGFBP5, using cultured fibroblast RNAs from two patients with distinct OBSL1 mutations (p.F697G; p.H814RfsX15). We found normal CUL7 mRNA levels but abnormal IGFBP2 and IGFBP5 mRNA levels in the two patients, suggesting that OBSL1 modulates the expression of IGFBP proteins. Hum Mutat 30:1–7, 2009. © 2009 Wiley‐Liss, Inc.     

Significant gene content variation characterizes the genomes of inbred mouse strains

The contribution to genetic diversity of genomic segmental copy number variations (CNVs) is less well understood than that of single-nucleotide polymorphisms (SNPs). While less frequent than SNPs, CNVs have greater potential to affect phenotype. In this study, we have performed the most comprehensive survey to date of CNVs in mice, analyzing the genomes of 42 Mouse Phenome Consortium priority strains. This microarray comparative genomic hybridization (CGH)-based analysis has identified 2094 putative CNVs, with an average of 10 Mb of DNA in 51 CNVs when individual mouse strains were compared to the reference strain C57BL/6J. This amount of variation results in gene content that can differ by hundreds of genes between strains. These genes include members of large families such as the major histocompatibility and pheromone receptor genes, but there are also many singleton genes including genes with expected phenotypic consequences from their deletion or amplification. Using a whole-genome association analysis, we demonstrate that complex multigenic phenotypes, such as food intake, can be associated with specific copy number changes.