Search

Discusses the principles of quantitative trait loci analysis. Occurrence of allelic differences in structural and regulatory genes; Factors involved in the alteration of gene actions; Genes that control root development and germination time.

Although commonplace in human disease genetics, genome-wide association (GWA) studies have only relatively recently been applied to plants. Using 32 phenotypes in the inbreeding crop barley, we report GWA mapping of 15 morphological traits across ∼500 cultivars genotyped with 1,536 SNPs. In contrast to the majority of human GWA studies, we observe high levels of linkage disequilibrium within and between chromosomes. Despite this, GWA analysis readily detected common alleles of high penetrance. To investigate the potential of combining GWA mapping with comparative analysis to resolve traits to candidate polymorphism level in unsequenced genomes, we fine-mapped a selected phenotype (anthocyanin pigmentation) within a 140-kb interval containing three genes. Of these, resequencing the putative anthocyanin pathway gene HvbHLH1 identified a deletion resulting in a premature stop codon upstream of the basic helix-loop-helix domain, which was diagnostic for lack of anthocyanin in our association and biparental mapping populations. The methodology described here is transferable to species with limited genomic resources, providing a paradigm for reducing the threshold of map-based cloning in unsequenced crops.

Recently, the debate on the centrality-lethality rule is resolved by the "second-generation" high-throughput Y2H data from the yeast interactome network, which suggests no significant correlation between the degree of connectedness and essentiality of proteins. However, it is still not clear why essential proteins strongly tend to interact with each other. Previously, the concept of essential protein-protein interactions was proposed to explain the mechanism underlying the clustering of essential proteins. In this article we show that 67 to 75% of the excessive interactions between essential proteins (IBEPs) in the yeast interactome network can be attributed to interactions within protein complexes characterised by the same deletion phenotype for subunits within the complex. Furthermore, 20 to 78% of the excess in IBEPs are caused by the strongly modular structure of the network and by variation in protein essentiality among modules. Not only do proteins function as an interactive network in cellular processes, but furthermore, many proteins do not take part in the network alone, but integrate into protein complexes and many functionally related complexes integrate into modules. So, the local structure of protein complexes as both functional and structural modules are the main contributory factors in the clustering of essential proteins.

The availability of reliable genetic linkage maps is crucial for functional and evolutionary genomic analyses. Established theory and methods of genetic linkage analysis have made map construction a routine exercise in diploids. However, many evolutionarily, ecologically, and/or agronomically important species are autopolyploids, with autotetraploidy being a typical example. These species undergo much more complicated chromosomal segregation and recombination at meiosis than diploids. In addition, there is evidence of polyploidy-induced and highly dynamic changes in the structure of the genome. These polysomic characteristics indicate the inappropriateness of the theory and methods of linkage analysis in diploids for use in these species and a gap in the theory and methodology of tetraploid map construction. This paper presents a theoretical model and statistical framework for multilocus linkage analysis in autotetraploids for use with dominant and/or codominant DNA molecular markers. The theory and methods incorporate the essential features of allele segregation and recombination under tetrasomic inheritance and the major challenges in statistical modeling and marker data analysis. We validated the method and explored its statistical properties by intensive simulation study and demonstrated its utility by analysis of AFLP and SSR marker data from an outbred autotetraploid potato population.

An alternative to complete genome sequencing is development and analysis of ESTs-fragments of transcribed coding DNA sequences. The EST collections also enhanced the development of cDNA microarray technologies, which make possible assessing the transcription levels of several thousand genes in a studied tissue of an organism in the same experiment. This paper provides an overview of the results of experiments with a barley microarray, Affymetrix Barley1 GeneChip. The variation in transcription levels of over 22000 genes in germinating barley grain of 150 barley double haploid lines produced by crossing cultivars Steptoe and Morex. Variation in gene expression of each gene is a quantitative trait, which can be mapped in population of double haploids as the genetic loci determining its variation (expressed QTL or eQTL). A regulatory locus (eQTL) can colocalize with the corresponding gene on genetic map (cis-eQTL) or be distant from it, frequently on another chromosome (trans-eQTL). Thus, it is possible to detect and analyze cis- and trans-regulatory loci for genes on a genome-wide scale. The design of the Affymetrix oligonucleotide arrays makes it possible not only to concurrently test the transcription level of several thousand genes, but also to simultaneously detect the polymorphic regions in cDNA sequences, thereby finding a considerable fraction of all nucleotide substitutions between the compared genotypes. Two types of data (the expression levels of several thousand genes and the presence of polymorphic sites in their sequences) can be obtained concurrently when processing the results of the same experiment. The details of both procedures are illustrated with explanatory examples.

It is well known that Affymetrix microarrays are widely used to predict genome-wide gene expression and genome-wide genetic polymorphisms from RNA and genomic DNA hybridization experiments, respectively. It has recently been proposed to integrate the two predictions by use of RNA microarray data only. Although the ability to detect single feature polymorphisms (SFPs) from RNA microarray data has many practical implications for genome study in both sequenced and unsequenced species, it raises enormous challenges for statistical modelling and analysis of microarray gene expression data for this objective. Several methods are proposed to predict SFPs from the gene expression profile. However, their performance is highly vulnerable to differential expression of genes. The SFPs thus predicted are eventually a reflection of differentially expressed genes rather than genuine sequence polymorphisms. To address the problem, we developed a novel statistical method to separate the binding affinity between a transcript and its targeting probe and the parameter measuring transcript abundance from perfect-match hybridization values of Affymetrix gene expression data. We implemented a Bayesian approach to detect SFPs and to genotype a segregating population at the detected SFPs. Based on analysis of three Affymetrix microarray datasets, we demonstrated that the present method confers a significantly improved robustness and accuracy in detecting the SFPs that carry genuine sequence polymorphisms when compared to its rivals in the literature. The method developed in this paper will provide experimental genomicists with advanced analytical tools for appropriate and efficient analysis of their microarray experiments and biostatisticians with insightful interpretation of Affymetrix microarray data.

Transcript abundance data from cRNA hybridizations to Affymetrix microarrays can potentially be used to identify genetic markers to facilitate high-throughput genotyping. We have shown that it is easily possible to use the information from Affymetrix expression arrays to accurately identify over 4,000 robust polymorphic transcript-derived markers (TDMs). We developed the method to identity TDM polymorphisms from experiments involving two tissues in two commercial varieties of barley and their doubled-haploid progeny. These TDMs represent ~18% of the total barley genes on the chip and can be used to predict the genotypes in an F(1)-derived, doubled-haploid population. According to our estimates, 35% of the TDMs reveal nucleotide polymorphism of the particular gene (single feature polymorphisms, SFPs) while 65% mark polymorphism resulting in extreme variation of gene expression (genetic expression markers, GEMs). These latter are probably mainly cis-acting regulators while a small proportion, approximately 5%, are loosely or un-linked transregulators.

The two genomes (A and C) of the allopolyploid Brassica napus have been clearly distinguished using genomic in situ hybridization (GISH) despite the fact that the two extant diploids, B. rapa (A, n = 10) and B. oleracea (C, n = 9), representing the progenitor genomes, are closely related. Using DNA from B. oleracea as the probe, with B. rapa DNA and the intergenic spacer of the B. oleracea 45S rDNA as the block, hybridization occurred on 9 of the 19 chromosome pairs along the majority of their length. The pattern of hybridization confirms that the two genomes have remained distinct in B. napus line DH12075, with no significant genome homogenization and no large-scale translocations between the genomes. Fluorescence in situ hybridization (FISH)-with 45S rDNA and a BAC that hybridizes to the pericentromeric heterochromatin of several chromosomes-followed by GISH allowed identification of six chromosomes and also three chromosome groups. Our procedure was used on the B. napus cultivar Westar, which has an interstitial reciprocal translocation. Two translocated segments were detected in pollen mother cells at the pachytene stage of meiosis. Using B. oleracea chromosome-specific BACs as FISH probes followed by GISH, the chromosomes involved were confirmed to be A7 and C6.

Background: A typical genetical genomics experiment results in four separate data sets; genotype, gene expression, higher-order phenotypic data and metadata that describe the protocols, processing and the array platform. Used in concert, these data sets provide the opportunity to perform genetic analysis at a systems level. Their predictive power is largely determined by the gene expression dataset where tens of millions of data points can be generated using currently available mRNA profiling technologies. Such large, multidimensional data sets often have value beyond that extracted during their initial analysis and interpretation, particularly if conducted on widely distributed reference genetic materials. Besides quality and scale, access to the data is of primary importance as accessibility potentially allows the extraction of considerable added value from the same primary dataset by the wider research community. Although the number of genetical genomics experiments in different plant species is rapidly increasing, none to date has been presented in a form that allows quick and efficient on-line testing for possible associations between genes, loci and traits of interest by an entire research community., Description: Using a reference population of 150 recombinant doubled haploid barley lines we generated novel phenotypic, mRNA abundance and SNP-based genotyping data sets, added them to a considerable volume of legacy trait data and entered them into the GeneNetwork http://www.genenetwork.org. GeneNetwork is a unified on-line analytical environment that enables the user to test genetic hypotheses about how component traits, such as mRNA abundance, may interact to condition more complex biological phenotypes (higher-order traits). Here we describe these barley data sets and demonstrate some of the functionalities GeneNetwork provides as an easily accessible and integrated analytical environment for exploring them., Conclusion: By integrating barley genotypic, phenotypic and mRNA abundance data sets directly within GeneNetwork's analytical environment we provide simple web access to the data for the research community. In this environment, a combination of correlation analysis and linkage mapping provides the potential to identify and substantiate gene targets for saturation mapping and positional cloning. By integrating datasets from an unsequenced crop plant (barley) in a database that has been designed for an animal model species (mouse) with a well established genome sequence, we prove the importance of the concept and practice of modular development and interoperability of software engineering for biological data sets.

We previously mapped mRNA transcript abundance traits (expression-QTL or eQTL) using the Barley1 Affymetrix array and 'whole plant' tissue from 139 progeny of the Steptoe x Morex (St/Mx) reference barley mapping population. Of the 22,840 probesets (genes) on the array, 15,987 reported transcript abundance signals that were suitable for eQTL analysis, and this revealed a genome-wide distribution of 23,738 significant eQTLs. Here we have explored the potential of using these mRNA abundance eQTL traits as surrogates for the identification of candidate genes underlying the interaction between barley and the wheat stem rust fungus Puccinia graminis f. sp. tritici. We re-analysed quantitative 'resistance phenotype' data collected on this population in 1990/1991 and identified six loci associated with barley's reaction to stem rust. One of these coincided with the major stem rust resistance locus Rpg1, that we had previously positionally cloned using this population. Correlation analysis between phenotype values for rust infection and mRNA abundance values reported by the 22,840 GeneChip probe sets placed Rpg1, which is on the Barley1 GeneChip, in the top five candidate genes for the major QTL on chromosome 7H corresponding to the location of Rpg1. A second co-located with the rpg4/Rpg5 stem rust resistance locus that has been mapped in a different population and the remaining four were novel. Correlation analyses identified candidate genes for the rpg4/Rpg5 locus on chromosome 5H. By combining our data with additional published mRNA profiling data sets, we identify a putative sensory transduction histidine kinase as a strong candidate for a novel resistance locus on chromosome 2H and compile candidate gene lists for the other three loci.

Background: Affymetrix high density oligonucleotide expression arrays are widely used across all fields of biological research for measuring genome-wide gene expression. An important step in processing oligonucleotide microarray data is to produce a single value for the gene expression level of an RNA transcript using one of a growing number of statistical methods. The challenge for the researcher is to decide on the most appropriate method to use to address a specific biological question with a given dataset. Although several research efforts have focused on assessing performance of a few methods in evaluating gene expression from RNA hybridization experiments with different datasets, the relative merits of the methods currently available in the literature for evaluating genome-wide gene expression from Affymetrix microarray data collected from real biological experiments remain actively debated., Results: The present study reports a comprehensive survey of the performance of all seven commonly used methods in evaluating genome-wide gene expression from a well-designed experiment using Affymetrix microarrays. The experiment profiled eight genetically divergent barley cultivars each with three biological replicates. The dataset so obtained confers a balanced and idealized structure for the present analysis. The methods were evaluated on their sensitivity for detecting differentially expressed genes, reproducibility of expression values across replicates, and consistency in calling differentially expressed genes. The number of genes detected as differentially expressed among methods differed by a factor of two or more at a given false discovery rate (FDR) level. Moreover, we propose the use of genes containing single feature polymorphisms (SFPs) as an empirical test for comparison among methods for the ability to detect true differential gene expression on the basis that SFPs largely correspond to cis-acting expression regulators. The PDNN method demonstrated superiority over all other methods in every comparison, whilst the default Affymetrix MAS5.0 method was clearly inferior., Conclusion: A comprehensive assessment of seven commonly used data extraction methods based on an extensive barley Affymetrix gene expression dataset has shown that the PDNN method has superior performance for the detection of differentially expressed genes.

Previous locations of flowering time (FT) QTL in several Brassica species, coupled with Arabidopsis synteny, suggest that orthologues of the genes FLC, FY or CONSTANS might be the candidates. We focused on FLC, and cloned paralogous copies in Brassica oleracea, obtained their genomic DNA sequences, and confirmed their locations relative to those of known FT-QTL by genetical mapping. They varied in total length mainly due to the variable size of the first and last introns. A high level of identity was observed among Brassica FLC genes at the amino acid level but non-synonymous differences were present. Comparative analysis of the promoter and intragenic regions of BoFLC paralogues with Arabidopsis FLC revealed extensive differences in overall structure and organisation but showed high conservation within those segments known to be essential in regulating FLC expression. Four B. oleracea FLC copies (BoFLC1, BoFLC3, BoFLC4 and BoFLC5) were located to their respective linkage groups based on allelic sequence variation in lines from a doubled haploid population. All except BoFLC4 were within the confidence intervals of known FT-QTL. Sequence data indicated that relevant non-synonymous polymorphisms were present between parents A12DHd and GDDH33 for BoFLC genes. However, BoFLC alleles segregated independently of FT in backcrosses while the study provided evidence that BoFLC4 and BoFLC5 contain premature stop codons and so could not contribute to flowering time variation. Therefore, there is strong evidence against any of the 4 BoFLC being FT-QTL candidates in this population.

Expression divergence of duplicate genes is widely believed to be important for their retention and evolution of new function, although the mechanism that determines their expression divergence remains unclear. We use a genetical genomics approach to explore divergence in genetical control of yeast duplicate genes created by a whole-genome duplication that occurred about 100 MYA and those with a younger duplication age. The analysis reveals that duplicate genes have a significantly higher probability of sharing common genetic control than pairs of singleton genes. The expression quantitative trait loci (eQTLs) have diverged completely for a high proportion of duplicate pairs, whereas a substantially larger proportion of duplicates share common regulatory motifs after 100 Myr of divergent evolution. The similarity in both genetical control and cis motif structure for a duplicate pair is a reflection of its evolutionary age. This study reveals that up to 20% of variation in expression between ancient duplicate gene pairs in the yeast genome can be explained by both cis motif divergence (approximately 8%) and by trans eQTL divergence (approximately 10%). Initially, divergence in all 3 aspects of cis motif structure, trans-genetical control, and expression evolves coordinately with the coding sequence divergence of both young and old duplicate pairs. These findings highlight the importance of divergence in both cis motif structure and trans-genetical control in the diverse set of mechanisms underlying the expression divergence of yeast duplicate genes.

The recent development of Affymetrix chips designed from assembled EST sequences has spawned considerable interest in identifying single-feature polymorphisms (SFPs) from transcriptome data. SFPs are valuable genetic markers that potentially offer a physical link to the structural genes themselves. However, most current SFP prediction methodologies were developed for sequenced species although SFPs are particularly valuable for species with complex and unsequenced genomes. To establish the sensitivity and specificity of prediction, we explored four methods for identifying SFPs from experiments involving two tissues in two commercial barleys and their doubled-haploid progeny. The methods were compared in terms of numbers of SFPs predicted and their ability to identify known sequence polymorphisms in the features, to confirm existing SNP genotypes and to match existing maps and individual haplotypes. We identified >4000 separate SFPs that accurately predicted the SNP genotype of >98% of the doubled-haploid (DH) lines. They were highly enriched for features containing sequence polymorphisms but all methods uniformly identified a majority of SFPs ( approximately 64%) in features for which there was no sequence polymorphism while 5% mapped to different locations, indicating that "SFPs" mainly represent polymorphism in cis-acting regulators. All methods are efficient and robust at predicting markers for gene mapping.

An international consortium has launched the whole-genome sequencing of potato, the fourth most important food crop in the world. Construction of genetic linkage maps is an inevitable step for taking advantage of the genome projects for the development of novel cultivars in the autotetraploid crop species. However, linkage analysis in autopolyploids, the kernel of linkage map construction, is theoretically challenging and methodologically unavailable in the current literature. We present here a theoretical analysis and a statistical method for tetrasomic linkage analysis with dominant and/or codominant molecular markers. The analysis reveals some essential properties of the tetrasomic model. The method accounts properly for double reduction and incomplete information of marker phenotype in regard to the corresponding phenotype in estimating the coefficients of double reduction and recombination frequency and in testing their significance by using the marker phenotype data. Computer simulation was developed to validate the analysis and the method and a case study with 201 AFLP and SSR markers scored on 228 full-sib individuals of autotetraploid potato is used to illustrate the utility of the method in map construction in autotetraploid species.

The close relationship between Brassica oleracea and Arabidopsis thaliana has been used to explore the genetic and physical collinearity of the two species, focusing on an inverted segmental chromosome duplication within linkage group O6 of B. oleracea. Genetic evidence suggests that these segments share a common origin with a region of Arabidopsis chromosome 1. Brassica oleracea and Arabidopsis bacterial artificial chromosome probes have been used for fluorescence in situ hybridization analysis of B. oleracea pachytene chromosomes to further characterize the inverted duplication. This has been highly effective in increasing the local resolution of the cytogenetic map. We have shown that the physical order of corresponding genetic markers is highly conserved between the duplicated regions in B. oleracea and the physical lengths of the regions at pachytene are similar, while the genetic distances are considerably different. The physical marker order is also well conserved between Arabidopsis and B. oleracea, with only one short inversion identified. Furthermore, the relative physical distances between the markers in one segment of B. oleracea and Arabidopsis have stayed approximately the same. The efficacy of using fluorescence in situ hybridization, together with other forms of physical and genetic mapping, for elucidating such issues relating to synteny is discussed.

The impact of the biological network structures on the divergence between the two copies of one duplicate gene pair involved in the networks has not been documented on a genome scale. Having analyzed the most recently updated Database of Interacting Proteins (DIP) by incorporating the information for duplicate genes of the same age in yeast, we find that there was a highly significantly positive correlation between the level of connectivity of ancient genes and the number of shared partners of their duplicates in the protein-protein interaction networks. This suggests that duplicate genes with a low ancestral connectivity tend to provide raw materials for functional novelty, whereas those duplicate genes with a high ancestral connectivity tend to create functional redundancy for a genome during the same evolutionary period. Moreover, the difference in the number of partners between two copies of a duplicate pair was found to follow a power-law distribution. This suggests that loss and gain of interacting partners for most duplicate genes with a lower level of ancestral connectivity is largely symmetrical, whereas the "hub duplicate genes" with a higher level of ancient connectivity display an asymmetrical divergence pattern in protein-protein interactions. Thus, it is clear that the protein-protein interaction network structures affect the divergence pattern of duplicate genes. Our findings also provide insights into the origin and development of biological networks.

The connexins are the subunits of a family of proteins that form gap junctions, allowing ions and small molecules to move between adjacent cells. At least four connexins are expressed in the ear, and, although there are known mutations at >100 loci that can cause deafness, those involving DFNB1, in the interval 13q11-q12 containing the GJB2 and GJB6 genes coding for connexins 26 and 30, are the most frequent cause of recessive deafness in many populations. We have suggested that the combined effects of relaxed selection and linguistic homogamy can explain the high frequency of connexin deafness and may have doubled its incidence in this country during the past 200 years. In this report, we show by computer simulation that assortative mating, in fact, can accelerate dramatically the genetic response to relaxed selection. Along with the effects of gene drift and consanguinity, assortative mating also may have played a key role in the joint evolution and accelerated fixation of genes for speech after they first appeared in Homo sapiens 100,000-150,000 years ago.

Linkage analysis in autotetraploid species has been an historical challenge in quantitative genetics theory and is a stumbling block that urgently needs to be removed in the rapidly emerging genome research on this species, such as cultivated potato. This article presents theory of a full model of tetrasomic linkage and develops a statistical framework for the linkage analysis. The model considers both double reduction and recombination, the most essential features of tetrasomic inheritance with linked loci, whereas the statistical method takes appropriate account of the major complexities in analyzing both dominant and codominant molecular marker data during map reconstruction in tetraploid species. These complexities include the problems arising from multiple dosage of allelic inheritance, the null allele, allelic segregation distortion, mixed bivalent and quadrivalent pairing in meiosis, and incomplete information of marker phenotype data. The theoretical analysis established the relationship between the coefficients of double reduction at linked loci, which is essential in the present tetrasomic linkage analysis and in assessing the impact of double reduction on the evolution of tetraploid populations. The statistical method, based on the combination of theoretical analysis and a computer-based algorithm, provided analytical tools for predicting the maximum-likelihood estimates of the model parameters. A simulation study showed the feasibility of a practical implementation of the method, detailed the procedure of the analysis, validated the power and reliability in the parameter estimation, and compared the present method with those proposed in the current literature.

A commonly encountered difficulty with the genetic engineering of crop plants is that different varieties of a particular species can show great variability in the efficiency with which they can be transformed. This increases the effort required to introduce transgenes into particular genetic backgrounds. The use of Substitution Lines has allowed the finer mapping of three Quantitative Trait Loci (tf1, tf2 and tf3) that explain 26% of the variation in the efficiency of Agrobacterium-mediated transformation in Brassica oleracea. Use of an 'orthogonal set' of genotypes (containing all eight possible combinations of 'positive' and 'negative' alleles at the three QTL), along with time course studies of transgene expression, has allowed the determination of the stages at which these genes have their effects during transformation. With regard to control of the level of transient transgene expression, tf1 (on LGO1) alone has no detectable effect, whilst tf2 (on LGO3) and tf3 (on LGO7) have highly significant effects (P < 0.001). All three loci have highly significant (P < 0.001) effects on the levels of expression of stably integrated transgene. The use of RFLP markers has shown that tf1 and tf2 are in duplicated regions of the B. oleracea genome and appear to be paralogous in origin. Colinearity of these regions with the A. thaliana genome has been identified. The results allow the selection of progeny Brassica oleracea genotypes that are more efficiently transformed than either parent used in the original cross.

A new source of resistance to the pathotype 4 isolate of Turnip mosaic virus (TuMV) CDN 1 has been identified in Brassica napus (oilseed rape). Analysis of segregation of resistance to TuMV isolate CDN 1 in a backcross generation following a cross between a resistant and a susceptible B. napus line showed that the resistance was dominant and monogenic. Molecular markers linked to this dominant resistance were identified using amplified fragment length polymorphism (AFLP) and microsatellite bulk segregant analysis. Bulks consisted of individuals from a BC(1) population with the resistant or the susceptible phenotype following challenge with CDN 1. One AFLP and six microsatellite markers were associated with the resistance locus, named TuRB03, and these mapped to the same region on chromosome N6 as a previously mapped TuMV resistance gene TuRB01. Further testing of TuRB03 with other TuMV isolates showed that it was not effective against all pathotype 4 isolates. It was effective against some, but not all pathotype 3 isolates tested. It provided further resolution of TuMV pathotypes by sub-dividing pathotypes 3 and 4. TuRB03 also provides a new source of resistance for combining with other resistances in our attempts to generate durable resistance to this virus.

The genetic control of 22 quantitative traits, including developmental rates and sizes, was examined in generations of Arabidopsis thaliana derived from the cross between the ecotypes, Columbia (Col) and Landsberg erecta (Ler). The data were obtained from three sets of families raised in the same trial: the 16 basic generations, that is, parents, F(1)'s, F(2)'s, backcrosses, recombinant inbred lines (RILs) and a triple test cross (TTC), the latter produced by crossing the RILs to Col, Ler and their F(1). The data were analysed by two approaches. The first (approach A) involved traditional generation mean and variance component analysis and the second (B), based around the RILs and TTC families, involved marker-based QTL analysis. From (A), genetic differences between Col and Ler were detected for all traits with moderate heritabilities. Height at flowering was the only trait to show heterosis. Dominance was partial to complete for all height traits, and there was no overdominance but there was strong evidence for directional dominance. For most other traits, dominance was ambidirectional and incomplete, with average dominance ratios of around 80%. Epistasis, particularly of the duplicate type that opposes dominance, was a common feature of all traits. The presence of epistasis must imply multiple QTL for all traits. The QTL analysis located 38 significant effects in four regions of chromosomes I, II, IV and V, but not III. QTL affecting rosette size and leaf number were identified in all four regions, with days to maturity on chromosomes IV and V. The only QTL for height was located at the expected position of the erecta gene (chromosome II; 50 cM), but the additive and dominance effects of this single QTL did not adequately explain the generation means. The possible involvement of other interacting height QTL is discussed.

Backcross breeding with marker-assisted selection was used to construct an intervarietal set of part chromosome substitution lines in Brassica napus, formed from a cross between two winter varieties of oilseed rape: Tapidor and Victor. A total of 22 lines from this substitution library were examined over a 3-year period, in a total of nine field trials, for seed oil fatty acid composition and seed oil content. Trialing of the substitution lines gave evidence for the existence of 13 quantitative trait loci (QTL). All 13 QTL affected fatty acid composition of the seed, and were distributed among linkage groups 1, 3, 6, 7, 8, 11, 13, 14, 18, and 19. Seven of these QTL, on linkage groups 3, 6, 8, 13, 14, 18, and 19, also affected total seed oil content. The positions of these QTL are compared to those in the published literature and with respect to erucic acid QTL previously identified in a backcross population of the same cross. The substitution line approach gives increased precision and sensitivity for QTL mapping compared to other methods.

Many biologically and economically important traits in plants and animals are quantitative/multifactorial, being controlled by several quantitative trait loci (QTL). QTL are difficult to locate accurately by conventional methods using molecular markers in segregating populations, particularly for traits of low heritability or for QTL with small effects. In order to resolve this, large (often unrealistically large) populations are required. In this paper we present an alternative approach using a specially developed resource of lines that facilitate QTL location first to a particular chromosome, then to successively smaller regions within a chromosome (< or = 0.5 cM) by means of simple comparisons among a few lines. This resource consists of "Stepped Aligned Inbred Recombinant Strains" (STAIRS) plus single whole Chromosome Substitution Strains (CSSs). We explain the analytical power of STAIRS and illustrate their construction and use with Arabidopsis thaliana, although the principles could be applied to many organisms. We were able to locate flowering QTL at the top of chromosome 3 known to contain several potential candidate genes.

We have assigned all nine linkage groups of a Brassica oleracea genetic map to each of the nine chromosomes of the karyotype derived from mitotic metaphase spreads of the B. oleracea var. alboglabra line A12DHd using FISH. The majority of probes were BACs, with A12DHd DNA inserts, which give clear, reliable FISH signals. We have added nine markers to the existing integrated linkage map, distributed over six linkage groups. BACs were definitively assigned to linkage map positions through development of locus-specific PCR assays. Integration of the cytogenetic and genetic linkage maps was achieved with 22 probes representing 19 loci. Four chromosomes (2, 4, 7, and 9) are in the same orientation as their respective linkage groups (O4, O7, O8, and O6) whereas four chromosomes (1, 3, 5, and 8) and linkage groups (O3, O9, O2, and O1) are in the opposite orientation. The remaining chromosome (6) is probably in the opposite orientation. The cytogenetic map is an important resource for locating probes with unknown genetic map positions and is also being used to analyze the relationships between genetic and cytogenetic maps.

Dissecting quantitative genetic variation into genes at the molecular level has been recognized as the greatest challenge facing geneticists in the twenty-first century. Tremendous efforts in the last two decades were invested to map a wide spectrum of quantitative genetic variation in nearly all important organisms onto their genome regions that may contain genes underlying the variation, but the candidate regions predicted so far are too coarse for accurate gene targeting. In this article, the recurrent selection and backcross (RSB) schemes were investigated theoretically and by simulation for their potential in mapping quantitative trait loci (QTL). In the RSB schemes, selection plays the role of maintaining the recipient genome in the vicinity of the QTL, which, at the same time, are rapidly narrowed down over multiple generations of backcrossing. With a high-density linkage map of DNA polymorphisms, the RSB approach has the potential of dissecting the complex genetic architecture of quantitative traits and enabling the underlying QTL to be mapped with the precision and resolution needed for their map-based cloning to be attempted. The factors affecting efficiency of the mapping method were investigated, suggesting guidelines under which experimental designs of the RSB schemes can be optimized. Comparison was made between the RSB schemes and the two popular QTL mapping methods, interval mapping and composite interval mapping, and showed that the scenario of genomic distribution of QTL that was unlocked by the RSB-based mapping method is qualitatively distinguished from those unlocked by the interval mapping-based methods.

A single chromosome of the grass species Festuca pratensis has been introgressed into Lolium perenne to produce a diploid monosomic substitution line 2n = 2x = 14. In this line recombination occurs throughout the length of the F. pratensis/L. perenne bivalent. The F. pratensis chromosome and recombinants between it and its L. perenne homeologue can be visualized using genomic in situ hybridization (GISH). GISH junctions represent the physical locations of sites of recombination, enabling a range of recombinant chromosomes to be used for physical mapping of the introgressed F. pratensis chromosome. The physical map, in conjunction with a genetic map composed of 104 F. pratensis-specific amplified fragment length polymorphisms (AFLPs), demonstrated: (1) the first large-scale analysis of the physical distribution of AFLPs; (2) variation in the relationship between genetic and physical distance from one part of the F. pratensis chromosome to another (e.g., variation was observed between and within chromosome arms); (3) that nucleolar organizer regions (NORs) and centromeres greatly reduce recombination; (4) that coding sequences are present close to the centromere and NORs in areas of low recombination in plant species with large genomes; and (5) apparent complete synteny between the F. pratensis chromosome and rice chromosome 1.

A single chromosome of the grass species Festuca pratensis has been introgressed into Lolium perenne to produce a diploid monosomic substitution line 2n = 2x = 14. The chromatin of F. pratensis and L. perenne can be distinguished by genomic in situ hybridization (GISH), and it is therefore possible to visualize the substituted F. pratensis chromosome in the L. perenne background and to study chiasma formation in a single marked bivalent. Recombination occurs freely in the F. pratensis/L. perenne bivalent, and chiasma frequency counts give a predicted map length for this bivalent of 76 cM. The substituted F. pratensis chromosome was also mapped with 104 EcoRI/Tru91 and HindIII/Tru91 amplified fragment length polymorphisms (AFLPs), generating a marker map of 81 cM. This map length is almost identical to the map length of 76 cM predicted from the chiasma frequency data. The work demonstrates a 1:1 correspondence between chiasma frequency and recombination and, in addition, the absence of chromatid interference across the Festuca and Lolium centromeres.

A segregating population of F(1)-derived doubled haploid (DH) lines of Brassica oleracea was used to detect and locate QTLs controlling 27 morphological and developmental traits, including leaf, flowering, axillary bud and stem characters. The population resulted from a cross between two very different B. oleracea crop types, an annual cauliflower and a biennial Brussels sprout. A principal component analysis (PCA), based on line means, allowed all the traits to be grouped into distinct categories according to the first five Principal Components. These were: leaf traits (PC1), flowering traits (PC2), axillary bud traits (PC3 and 5) and stem traits (PC4). Between zero and four putative QTL were located per trait, which individually explained between 6% and 43% of the additive genetic variation, using the multiple-marker regression approach to QTL mapping. For lamina width, bare petiole length and stem length two QTL with opposite effects were detected on the same linkage groups. Intra- and inter-specific comparative mapping using RFLP markers identified a QTL on linkage group O8 accounting for variation in vernalisation, which is probably synonymous with a QTL detected on linkage group N19 of Brassica napus. In addition, a QTL for petiole length detected on O3 of this study appeared to be homologous to a QTL detected on another B. oleracea genetic map (Camargo et al. 1995).

Seventy-nine recombinant backcross substitution lines from a cross between Brassica oleracea var. italica and Brassica oleracea var. alboglabra were grown in field trials over five years along with the alboglabra recurrent parent. Plants were scored for the days from sowing to the opening of the first flower, and lines that flowered significantly earlier or later than the recurrent parent were identified. Based on the lengths of the substitutions, evidence for 11 QTL on chromosomes O1, O2, O3, O5 and O9 was found, five of which mapped to similar regions to five of the six found in a previous analysis of doubled haploid lines from the same cross. Several of the QTL were linked closely in repulsion.

A population of 150 doubled haploid lines of rapid cycling Brassica oleracea, derived from an F1 from a var. alboglabra x var. italica cross, was scored for flowering time in two trials. Using information on 82 mapped molecular markers, spread evenly across the nine linkage groups, QTL were identified at six locations; one each on linkage groups O2 and O3 and two each on linkage groups O5 and O9. In total, these QTL explained 58 and 93% of the genetical variation in the two trials. Three of these QTL, on linkage groups O2, O3, and O9, were situated in regions showing considerable homology both with each other and with chromosome regions of B. nigra that have been shown to affect flowering time. These same regions are all homologous to a single tract of Arabidopsis chromosome 5, which contains a number of the flowering-related genes, one or more of which may be candidates for the QTL found in Brassica.

We have briefly reviewed the methods currently available for QTL analysis in segregating populations and summarized some of the conclusions arising from such analyses in plant populations. We show that the analytical methods locate QTL with poor precision (10-30 cM), unless the heritability of an individual QTL is high. Also the estimates of the QTL effects, particularly the dominance effects tend to be inflated because only large estimates are significant. Estimates of numbers of QTL per trait are generally low (< 8) for individual trials. This may suggest that there are few QTL but probably reflects the power of the methods. There is no large correlation between the numbers of QTL found and the amount of the variation explained. Of those cases where dominance is measurable, dominance ratios are often > 1, but seldom significantly greater. These latter cases need further analysis. Many QTL map close to candidate genes, and there is growing evidence from synteny studies of corresponding chromosome regions carrying similar QTL in different species. However, unreliability of QTL location may suggest false candidates.

The currently available methods for locating quantitative trait loci (QTLs) and measuring their effects in segregating populations lack precision unless individual QTLs have very high heritabilities. The use of recombinant backcross lines containing short regions of donor chromosome introgressed into a constant recipient background permits QTLs to be located with greater precision. The present paper describes the use of molecular markers to introgress defined short regions of chromosome from a donor doubled haploid calabrese line of Brassica oleracea (var. italica) into a recipient short generation variety (Brassica oleracea var. alboglabra). We demonstrate that in just two or three generations of backcrossing, combined with selection for mapped molecular markers, the generation of a library of recombinant backcross lines is feasible. The possible use and refinement of these lines are discussed. Key words : backcrossing, Brassica oleracea, introgression, molecular markers, near-isogenic lines, QTL mapping, recombinant backcross lines, substitution lines.

Linkage maps of the nine chromosomes of Brassica oleracea, based on 75 informative molecular markers, have been compared in first and second backcross progeny from a cross between two doubled haploid lines. The second backcross progeny showed greater recombination frequencies for 75% of the pairs of adjacent markers, but there was no obvious indication that this effect was localised to particular regions of the chromosomes. Four chromosomes increased in genetic length more than twofold, while overall, the total map was 66% longer. The possible causes of this discrepancy are analysed. A sex difference in chiasma distribution and/or frequency at meiosis is thought to be the most likely explanation. The implications of this finding for mapping and map-based applications are discussed.

A variety of approaches are available for identifying the location and effect of QTL in segregating populations using molecular markers. However, these have problems in distinguishing two linked QTL, particularly in relation to the size of the test statistic when many independent tests are performed. An empirical method for obtaining the distribution of the test statistic for specific datasets is described, and its power for demonstrating the inadequacy of a single-QTL model is explored through computer simulation. The method is an extension of the previously described technique of 'marker regression', and it is applied here to demonstrate two situations in which it may be useful. Firstly, we examine the power of the technique to distinguish two, linked QTL from one and compare this ability with that of two contemporary methods, 'Mapmaker/QTL' and 'regression mapping'. Secondly, we show how to combine information from two, or more, populations that may be segregating for different marker loci in a given linkage group. This is illustrated for two populations having in common just two linked marker loci although the sharing of loci is not a pre-requisite. Empirical tests are used to determine whether the same or different QTL are segregating and, if they are the same QTL, whether they are the same alleles. Evidence is discussed which suggests that the upper limit to the number of QTL that can be located for any single quantitative trait in a segregating populations is 12.

F1 plants between two intervarietal chromosome substitution lines of European spring wheat varieties, 'Sicco' ('Chinese Spring 5B') and 'Highbury' ('Chinese Spring 5B'), were used to produce 114 doubled haploid lines, 45 by the Hordeum bulbosum technique and 69 by anther culture. These two sets of lines were characterized for variation at a range of morphological, isozyme and RFLP marker loci, and genetic maps were developed with emphasis on chromosomes 6B, 7A, 7B and 7D. A subset of lines, scored for production traits in field trials in 1986 and 1987, were analysed for quantitative trait loci (QTL). The performance of the lines for the quantitative traits studied showed no overall differences due to the method of production of the lines. QTL were located on the linkage map for ear emergence time, height, tiller weight, yield and 50-grain weight using four analytical methods. Many of these effects showed genotype x year interaction.

A method to locate quantitative trait loci (QTL) on a chromosome and to estimate their additive and dominance effects is described. It applies to generations derived from an F1 by selfing or backcrossing and to doubled haploid lines, given that marker genotype information (RFLP, RAPD, etc.) and quantitative trait data are available. The method involves regressing the additive difference between marker genotype means at a locus against a function of the recombination frequency between that locus and a putative QTL. A QTL is located, as by other regression methods, at that point where the residual mean square is minimised. The estimates of location and gene effects are consistent and as reliable as conventional flanking-marker methods. Further applications include the ability to test for the presence of two, or more, linked QTL and to compare different crosses for the presence of common QTL. Furthermore, the technique is straightforward and may be programmed using standard pc-based statistical software.

The algorithm for estimating both the recombination fraction between a marker gene and a locus affecting a quantitative trait, and also the means and variances of the QTL genotypes, is extended to backcross and doubled haploid populations. The simulation experiments show that estimates of these parameters can be obtained with acceptable accuracy and results are compared with those obtained using F2 populations studied previously (Luo & Kearsey, 1989).

A design and model are presented to allow the prediction, in early generations, of the mean and distribution of recombinant inbred lines derived from a cross between two parental populations or partially inbred lines. The procedure has been tested in autumn-sown onions (in the UK) using a wide cross between the openpollinated Japanese cultivar, Senshyu, and a partially inbred line derived from the European cultivar, Rawska. The early generations used for prediction included the first self-pollinated generation of the two parental populations and the F3 generation produced from the hybrid population. The predictions were tested by reference to the field performance of a random array of inbred lines, which were produced by single-seed descent (SSD) and had been selfed for three generations. The early generations, used for prediction, and a sample of SSD lines were raised alongside each other in each of two seasons. Within each season, good agreement was found between the predicted and observed performance of the recombinant inbred lines for three characters - yield, quality and maturity. This is used as evidence of the validity of the genetical model and the assumptions made. The effects of genotype x environment interactions prevented predictions made in one season being reliably applied to those made in the other and, therefore, reduce the attraction of this type of prediction study to the plant breeder.

Performance of a random array of recombinant inbred lines derived by single seed descent from five different source populations of Brussels sprouts (Brassica oleracea var. gemmifera) is presented. A total of 2,356 lines were tested in trials during 1985 and 1986. Three of the source populations were derived from double crosses between F1 hybrids. These hybrids show a considerable heterotic advantage over their inbred parents for the most important agronomic traits. The recombinant inbred lines performed, on average, less well than the parental inbred material, indicating that additive x additive genie interactions may make a significant contribution to the performance of current inbred material. Nevertheless, the very large variation among the recombinant inbred lines permitted many lines to be identified which outperformed the best parental inbred for all traits. Two lines outperformed the reference F1 hybrid, Gower, for an index that included marketable yield and quality. Consideration was also given to the dangers of misinterpreting phenotypically based proportions. Accordingly, response equations were used to ascertain the real genetic progress that was made. Advance seemed small when compared with the large heterotic effect, which is consistent with the segregation of a large number of loci. The distribution of the recombinant inbred lines was compared to predictions made from early generation trials. There was broad agreement but significant discrepancies existed which, it is suggested, may arise from the effects of genotype-environment interactions.

A maximum likelihood approach is developed for estimating the recombination fraction in a segregating population (F2), between a marker gene and a locus affecting a quantitative trait as well as estimating the means and variances of the three genotypes of the quantitative trait. The experimental results from computer simulations show that even with experimental sizes of 500, estimates of the parameters can be obtained by aid of the codominant marker gene as long as the heritability of the quantitative trait in question is not less than 0.10. However at low heritabilities the variances of estimates are very large.

Crosses between laboratory stocks and extractions from wild populations have recently been shown to produce non-reciprocal genetic aberrations commonly termed hybrid dysgenesis, which appear to arise from a nuclear cytoplasmic interaction. Female sterility is one aspect which has been investigated and both poor egg production (GD sterility) and low hatchability (SF sterility) have been shown to contribute. It has previously been suggested that these characters may have an independent action and causation. The results presented in this paper however, show a high degree of correlation in the response of SF and GD sterility to various developmental temperature regimes, with both forms of sterility showing an increase as the developmental temperature rises. For each character, the whole of the life cycle appears to be affected by changes in the developmental temperature although two stages were identified as being particularly sensitive. The results therefore suggest that these two characters have a common cause and the relationship between these and other dysgenic traits is discussed.