An Evaluation of High-Throughput Approaches to QTL Mapping in Saccharomyces cerevisiae

Dissecting the molecular basis of quantitative traits is a significant challenge and is essential for understanding complex diseases. Even in model organisms, precisely determining causative genes and their interactions has remained elusive, due in part to difficulty in narrowing intervals to single genes and in detecting epistasis or linked quantitative trait loci. These difficulties are exacerbated by limitations in experimental design, such as low numbers of analyzed individuals or of polymorphisms between parental genomes. We address these challenges by applying three independent high-throughput approaches for QTL mapping to map the genetic variants underlying 11 phenotypes in two genetically distant Saccharomyces cerevisiae strains, namely (1) individual analysis of >700 meiotic segregants, (2) bulk segregant analysis, and (3) reciprocal hemizygosity scanning, a new genome-wide method that we developed. We reveal differences in the performance of each approach and, by combining them, identify eight polymorphic genes that affect eight different phenotypes: colony shape, flocculation, growth on two nonfermentable carbon sources, and resistance to two drugs, salt, and high temperature. Our results demonstrate the power of individual segregant analysis to dissect QTL and address the underestimated contribution of interactions between variants. We also reveal confounding factors like mutations and aneuploidy in pooled approaches, providing valuable lessons for future designs of complex trait mapping studies.     

Rapid Identification of Major-Effect Genes Using the Collaborative Cross

The Collaborative Cross (CC) was designed to facilitate rapid gene mapping and consists of hundreds of recombinant inbred lines descended from eight diverse inbred founder strains. A decade in production, it can now be applied to mapping projects. Here, we provide a proof of principle for rapid identification of major-effect genes using the CC. To do so, we chose coat color traits since the location and identity of many relevant genes are known. We ascertained in 110 CC lines six different coat phenotypes: albino, agouti, black, cinnamon, and chocolate coat colors and the white-belly trait. We developed a pipeline employing modifications of existing mapping tools suitable for analyzing the complex genetic architecture of the CC. Together with analysis of the founders’ genome sequences, mapping was successfully achieved with sufficient resolution to identify the causative genes for five traits. Anticipating the application of the CC to complex traits, we also developed strategies to detect interacting genes, testing joint effects of three loci. Our results illustrate the power of the CC and provide confidence that this resource can be applied to complex traits for detection of both qualitative and quantitative trait loci.     

There is more to tomato fruit colour than candidate carotenoid genes

Determining gene sequences responsible for complex phenotypes has remained a major objective in modern biology. The candidate gene approach is attempting to link, through mapping analysis, sequences that have a known functional role in the measured phenotype with quantitative trait loci (QTL) that are responsible for the studied variation. To explore the potential of the candidate approach for complex traits we conducted a mapping analysis of QTL for the intensity of the red colour of the tomato fruit (mainly lycopene) and for probes associated with the well-characterized carotenoid biosynthesis pathway. Seventy-five tomato introgression lines (ILs), each containing a single homozygous RFLP-defined chromosome segment from the green-fruited species Lycopersicon pennellii delimited 107 marker-defined mapping bins. Three of the bins resolved known qualitative colour mutations for yellow (r) and orange (B and Del) fruits resulting from variation in specific carotenoid biosynthesis genes. Based on trials in different environments, 16 QTL that modified the intensity of the red colour of ripe fruit were assigned to bins. Candidate sequences associated with the carotenoid biosynthesis pathway were mapped to 23 loci. Only five of the QTL co-segregated with the same bins that contained candidate genes - a number that is expected by chance alone. Furthermore, similar map location of a QTL and a candidate is far from a direct causative relationship between a gene and a phenotype. This study highlights the wealth and complexity of the variation present in the genus Lycopersicon that could be employed for basic research and genetic improvement of fruit colour in tomato.     

Meiotic linkage mapping of 52 genes onto the canine map does not identify significant levels of microrearrangement

In an effort to extend our understanding of the evolutionary relationship between the canine and human genomes, we have developed and positioned 52 new gene-associated polymorphic markers on the canine meiotic linkage map. Canine-specific PCR primers were developed from the consensus of published sequences of several mammalian genomes and were designed to span intronic regions, thus optimizing the probability that a polymorphic site was included. The resulting markers were analyzed on a panel of three-generation canine reference families and the data were incorporated into the current meiotic linkage map. The data were compared with those generated by three chromosome paint studies in an effort to understand the distribution and frequency of microrearrangements within the canine genome. Forty-eight of 52 genes map to a chromosomal region predicted to contain genes from the corresponding region of the human genome according to all published reciprocal chromosome paint studies. Meiotic linkage mapping data for three genes can be used to resolve discrepancies between the published reciprocal chromosome paint studies, and for an additional two genes, meiotic mapping data allow evolutionary breakpoints to be more precisely defined. We conclude that microrearrangements of evolutionarily conserved segments between the canine and human genomes are rare, occurring for less than 0.5% of gene data reported to date. In addition, we have found that the placement of genes on the meiotic linkage map is a useful mechanism for resolving discrepancies between existing data sets. Received: 7 February 2001 / Accepted: 9 May 2001     

A Simple Model for the Influence of Meiotic Conversion Tracts on GC Content

A strong correlation between GC content and recombination rate is observed in many eukaryotes, which is thought to be due to conversion events linked to the repair of meiotic double-strand breaks. In several organisms, the length of conversion tracts has been shown to decrease exponentially with increasing distance from the sites of meiotic double-strand breaks. I show here that this behavior leads to a simple analytical model for the evolution and the equilibrium state of the GC content of sequences devoid of meiotic double-strand break sites. In the yeast Saccharomyces cerevisiae, meiotic double-strand breaks are practically excluded from protein-coding sequences. A good fit was observed between the predictions of the model and the variations of the average GC content of the third codon position (GC3) of S. cerevisiae genes. Moreover, recombination parameters that can be extracted by fitting the data to the model coincide with experimentally determined values. These results thus indicate that meiotic recombination plays an important part in determining the fluctuations of GC content in yeast coding sequences. The model also accounted for the different patterns of GC variations observed in the genes of Candida species that exhibit a variety of sexual lifestyles, and hence a wide range of meiotic recombination rates. Finally, the variations of the average GC3 content of human and chicken coding sequences could also be fitted by the model. These results suggest the existence of a widespread pattern of GC variation in eukaryotic genes due to meiotic recombination, which would imply the generality of two features of meiotic recombination: its association with GC-biased gene conversion and the quasi-exclusion of meiotic double-strand breaks from coding sequences. Moreover, the model points out to specific constraints on protein fragments encoded by exon terminal sequences, which are the most affected by the GC bias.     

Candidate genes for drought tolerance and improved productivity in rice (Oryza sativa L.)

Candidate genes are sequenced genes of known biological action involved in the development or physiology of a trait. Twenty-one putative candidate genes were designed after an exhaustive search in the public databases along with an elaborate literature survey for candidate gene products and/or regulatory sequences associated with enhanced drought resistance. The downloaded sequences were then used to design primers considering the flanking sequences as well. Polymerase chain reaction (PCR) performed on 10 diverse cultivars that involvedJaponica, Indica and local accessions, revealed 12 polymorphic candidate genes. Seven polymorphic candidate genes were then utilized to genotype 148 individuals of CT9993 × IR62266 doubled haploid (DH) mapping population. The segregation data were tested for deviation from the expected Mendelian ratio (1:1) using a Chi-square test (<1%). Based on this, four candidate genes were assessed to be significant and the remaining three, as non-significant. All the significant candidate genes were biased towards CT9993, the female parent in the DH mapping population. Single-marker analysis strongly associated (<1%) them to different traits under both well-watered and low-moisture stress conditions. Two candidate genes,EXP15 andEXP13, were found to be associated with root number and silicon content in the stem respectively, under both well-watered and low-moisture stress conditions     

Mendelian analysis of the organization of actin sequences in Physarum polycephalum

The genormic organization of the multiple actin DNA sequences in the lower eukaryote Physarum polycephalum was investigated by Mendelian mapping. Actin-homologous restriction endonuclease cleavage fragments detected by DNA blotting showed length polymorphisms when different strains were compared. These length polymorphisms were used as phenotypic markers for actin sequences in the genome. The meiotic assortment of the polymorphic restriction fragments was analysed, revealing four unlinked actin loci. The data for three-of the actin loci, ardB, C and D, are consistent with a single sequence or gene at each locus. The data for the other actin locus. ardA, is consistent with multiple linked actin sequences or genes.     

Quantitative trait loci and candidate genes underlying genotype by environment interaction in the response of Arabidopsis thaliana to drought

Drought stress was imposed on two sets of Arabidopsis thaliana genotypes grown in sand under short‐day conditions and analysed for several shoot and root growth traits. The response to drought was assessed for quantitative trait locus (QTL) mapping in a genetically diverse set of Arabidopsis accessions using genome‐wide association (GWA) mapping, and conventional linkage analysis of a recombinant inbred line (RIL) population. Results showed significant genotype by environment interaction (G×E) for all traits in response to different watering regimes. For the RIL population, the observed G×E was reflected in 17 QTL by environment interactions (Q×E), while 17 additional QTLs were mapped not showing Q×E. GWA mapping identified 58 single nucleotide polymorphism (SNPs) associated with loci displaying Q×E and an additional 16 SNPs associated with loci not showing Q×E. Many candidate genes potentially underlying these loci were suggested. The genes for RPS3C and YLS7 were found to contain conserved amino acid differences when comparing Arabidopsis accessions with strongly contrasting drought response phenotypes, further supporting their candidacy. One of these candidate genes co‐located with a QTL mapped in the RIL population.     

Efficient Plant Gene Identification Based on Interspecies Mapping of Full-Length cDNAs

We present an annotation pipeline that accurately predicts exon–intron structures and protein-coding sequences (CDSs) on the basis of full-length cDNAs (FLcDNAs). This annotation pipeline was used to identify genes in 10 plant genomes. In particular, we show that interspecies mapping of FLcDNAs to genomes is of great value in fully utilizing FLcDNA resources whose availability is limited to several species. Because low sequence conservation at 5′- and 3′-ends of FLcDNAs between different species tends to result in truncated CDSs, we developed an improved algorithm to identify complete CDSs by the extension of both ends of truncated CDSs. Interspecies mapping of 71 801 monocot FLcDNAs to the Oryza sativa genome led to the detection of 22 142 protein-coding regions. Moreover, in comparing two mapping programs and three ab initio prediction programs, we found that our pipeline was more capable of identifying complete CDSs. As demonstrated by monocot interspecies mapping, in which nucleotide identity between FLcDNAs and the genome was ∼80%, the resultant inferred CDSs were sufficiently accurate. Finally, we applied both inter- and intraspecies mapping to 10 monocot and dicot genomes and identified genes in 210 551 loci. Interspecies mapping of FLcDNAs is expected to effectively predict genes and CDSs in newly sequenced genomes.     

Development of Polymorphic EST Markers Suitable for Genetic Linkage Mapping of Catfish

Expressed sequence tag (EST) markers are important for gene mapping and for marker-assisted selection (MAS). To develop EST markers for use in catfish gene mapping, 100 randomly picked complementary DNAs from the channel catfish (Ictalurus punctatus) pituitary library were sequenced. The EST sequences were used to design primers to amplify channel catfish and blue catfish (I. furcatus) genomic DNAs. Polymerase chain reaction products of the ESTs were analyzed to determine length polymorphism between the channel catfish and blue catfish. Eleven polymorphic EST markers were identified. Five of the 11 EST markers were from known genes and the other six were from unidentified ESTs. Seven ESTs were found to be associated with microsatellite sequences. Analysis of channel catfish gene sequences indicated highly biased codon usage, with 16 codons being preferably used. These codons were more preferably used in highly expressed ribosomal protein genes and in highly expressed pituitary hormone genes. G/C-rich codons are less used in channel catfish than those in other vertebrates suggesting AT-richness of the channel catfish genome. Received June 29, 1998; accepted March 29, 1999.     

Molecular and Genetic Analysis of Rec103, an Early Meiotic Recombination Gene in Yeast

In the yeast Saccharomyces cerevisiae at least 10 genes are required to begin meiotic recombination. A new early recombination gene REC103 is described in this paper. It was initially defined by the rec103-1 mutation found in a selection for mutations overcoming the spore inviability of a rad52 spo13 haploid strain. Mutations in REC103 also rescue rad52 in spo13 diploids. rec103 spo13 strains produce viable spores; these spores show no evidence of meiotic recombination. rec103 SPO13 diploids produce no viable spores, consistent with the loss of recombination. Mutations in REC103 do not affect mitotic recombination, growth, or repair. These phenotypes are identical to those conferred by mutations in several other early meiotic recombination genes (e.g., REC102, REC104, REC114, MEI4, MER2, and SPO11). REC103 maps to chromosome VII between ADE5 and RAD54. Cloning and sequencing of REC103 reveals that REC103 is identical to SKI8, a gene that depresses the expression of yeast double-stranded (``killer') (ds)RNA viruses. REC103/SKI8 is transcribed in mitotic cells and is induced ~15-fold in meiosis. REC103 has 26% amino acid identity to the Schizosaccharomyces pombe rec14(+) gene; mutations in both genes confer similar meiotic phenotypes, suggesting that they may play similar roles in meiotic recombination.     

Targeted single nucleotide polymorphism (SNP) discovery in a highly polyploid plant species using 454 sequencing

Discovering single nucleotide polymorphisms (SNPs) in specific genes in a heterozygous polyploid plant species, such as sugarcane, is challenging because of the presence of a large number of homologues. To discover SNPs for mapping genes of interest, 454 sequencing of 307 polymerase chain reaction (PCR) amplicons (> 59 kb of sequence) was undertaken. One region of a four-gasket sequencing run, on a 454 Genome Sequencer FLX, was used for pooled PCR products amplified from each parent of a quantitative trait locus (QTL) mapping population (IJ76-514 × Q165). The sequencing yielded 96 755 (IJ76-514) and 86 241 (Q165) sequences with perfect matches to a PCR primer used in amplification, with an average sequence depth of approximately 300 and an average read length of 220 bases. Further analysis was carried out on amplicons whose sequences clustered into a single contig using an identity of 80% with the program cap 3. In the more polymorphic sugarcane parent (Q165), 94% of amplicons (227/242) had evidence of a reliable SNP – an average of one every 35 bases. Significantly fewer SNPs were found in the pure Saccharum officinarum parent – with one SNP every 58 bases and SNPs in 86% (213/247) of amplicons. Using automatic SNP detection, 1632 SNPs were detected in Q165 sequences and 1013 in IJ76-514. From 225 candidate SNP sites tested, 209 (93%) were validated as polymorphic using the Sequenom MassARRAY system. Amplicon re-sequencing using the 454 system enables cost-effective SNP discovery that can be targeted to genes of interest and is able to perform in the highly challenging area of polyploid genomes.     

Chromosomal localization of actin genes in the malaria mosquito Anopheles darlingi

Physical and genetic maps have been used for chromosomal localization of genes in vectors of infectious diseases. The availability of polytene chromosomes in malaria mosquitoes provides a unique opportunity to precisely map genes of interest. We report the physical mapping of two actin genes on polytene chromosomes of the major malaria vector in the Amazon, Anopheles darlingi (Diptera: Culicidae). Clones with actin gene sequences were obtained from a cDNA library constructed from RNA isolated from adult females and males of An. darlingi. Each of the two clones was mapped to a unique site on chromosomal arm 2L in subdivisions 21A (clone pl05‐A04) and 23B (clone pl17‐G06). The obtained results, together with previous mapping data, provide a suitable basis for comparative genomics and for establishing chromosomal homologies among major malaria vectors.     

Mapping of genes expressed in activated porcine Peyer's patch

To determine the chromosomal locations for genes expressed in porcine Peyer's patches, polymerase chain reaction-based mapping of expressed sequence tags (ESTs) isolated from a porcine Peyer's patch-specific cDNA library was performed across a 6500-rad swine radiation hybrid panel. A total of 116 ESTs were mapped with LOD scores >6.0, and another 11 ESTs had LOD scores between 5.0 and 6.0. Of these 127 ESTs, 63% matched known genes (相似文献   

Meiotic mapping of yeast ribosomal deoxyribonucleic acid on chromosome XII.

We have used meiotic mapping techniques to locate the position of the repeating ribosomal DNA (rDNA) genes of the yeast Saccharomyces cerevisiae. We found that the rDNA genes are located on the right arm of chromosome XII, approximately 45 map units centromere distal to the gene gal2. Together with mapping data from previous studies, this result suggests that the tandem array of rDNA genes contains at least two junctions with the non-rDNA of the yeast chromosome. In addition, we observed segregation patterns of the rDNA genes consistent with meiotic recombination within the rDNA gene tandem array in 3 of the 59 tetrads examined.     

Putative in silico mapping of DNA sequences to livestock genome maps using SSLP flanking sequences

In this study, an in silico approach was developed to identify homologies existing between livestock microsatellite flanking sequences and GenBank nucleotide sequences. Initially, 1955 bovine, 1570 porcine and 1121 chicken microsatellites were downloaded and the flanking sequences were compared with the nr and dbEST databases of GenBank. A total of 74 bovine, 44 porcine and 37 chicken microsatellite flanking sequences passed our criteria and had at least one significant match to human genomic sequence, genes/expressed sequence tags (ESTs) or both. GenBank annotation and BLAT searches of the UCSC human genome assembly revealed that 38 bovine, 13 porcine and 17 chicken microsatellite flanking sequences were highly similar to known human genes. Map locations were available for 67 bovine, 44 porcine and 21 chicken microsatellite flanking sequences, providing useful links in the comparative maps of humans and livestock. In support of our approach, 112 alignments with both microsatellite and match mapping information were located in the expected chromosomal regions based on previously reported syntenic relationships. The development of this in silico mapping approach has significantly increased the number of genes and EST sequences anchored to the bovine, porcine and chicken genome maps and the number of links in various human-livestock comparative maps.     

Structural analysis of a Lotus japonicus genome. I. Sequence features and mapping of fifty-six TAC clones which cover the 5.4 mb regions of the genome.

A total of 56 TAC clones with an average insert size of 100 kb were isolated from a TAC library of the Lotus japonicus genome based on the expressed sequences tags (ESTs), cDNA and gene information, and their nucleotide sequences were determined according to the shot-gun based strategy. The total length of the sequenced regions is 5,473,195 bp. By comparison with the sequences in protein and EST databases and analysis with computer programs for gene modeling, a total of 605 potential protein-encoding genes with known or predicted functions, 69 gene segments, and 172 pseudogenes were identified. The average density of the genes assigned so far is 1 gene/8120 bp. Introns were identified in approximately 78% of the potential genes. There was an average of 3.8 introns per gene and the average length of the introns was 375 bp. DNA markers were generated based on the nucleotide sequences obtained, and each clone was mapped onto the linkage map using the F2 mapping population derived from a cross of L. japonicus Gifu B-129 and Miyakojima MG-20. The sequence data, gene information and mapping information are available through the World Wide Web at http://www.kazusa.or.jp/lotus/.