Linkage analysis of molecular markers and quantitative trait loci in populations of inbred backcross lines of Brassica napus L.

Backcross populations are often used to study quantitative trait loci (QTL) after they are initially discovered in balanced populations, such as F(2), BC(1), or recombinant inbreds. While the latter are more powerful for mapping marker loci, the former have the reduced background genetic variation necessary for more precise estimation of QTL effects. Many populations of inbred backcross lines (IBLs) have been developed in plant and animal systems to permit simultaneous study and dissection of quantitative genetic variation introgressed from one source to another. Such populations have a genetic structure that can be used for linkage estimation and discovery of QTL. In this study, four populations of IBLs of oilseed Brassica napus were developed and analyzed to map genomic regions from the donor parent (a winter-type cultivar) that affect agronomic traits in spring-type inbreds and hybrids. Restriction fragment length polymorphisms (RFLPs) identified among the IBLs were used to calculate two-point recombination fractions and LOD scores through grid searches. This information allowed the enrichment of a composite genetic map of B. napus with 72 new RFLP loci. The selfed and hybrid progenies of the IBLs were evaluated during two growing seasons for several agronomic traits. Both pedigree structure and map information were incorporated into the QTL analysis by using a regression approach. The number of QTL detected for each trait and the number of effective factors calculated by using biometrical methods were of similar magnitude. Populations of IBLs were shown to be valuable for both marker mapping and QTL analysis.

A bayesian approach to detect quantitative trait loci using Markov chain Monte Carlo.

Markov chain Monte Carlo (MCMC) techniques are applied to simultaneously identify multiple quantitative trait loci (QTL) and the magnitude of their effects. Using a Bayesian approach a multi-locus model is fit to quantitative trait and molecular marker data, instead of fitting one locus at a time. The phenotypic trait is modeled as a linear function of the additive and dominance effects of the unknown QTL genotypes. Inference summaries for the locations of the QTL and their effects are derived from the corresponding marginal posterior densities obtained by integrating the likelihood, rather than by optimizing the joint likelihood surface. This is done using MCMC by treating the unknown QTL, genotypes, and any missing marker genotypes, as augmented data and then by including these unknowns in the Markov chain cycle alone with the unknown parameters. Parameter estimates are obtained as means of the corresponding marginal posterior densities. High posterior density regions of the marginal densities are obtained as confidence regions. We examine flowering time data from double haploid progeny of Brassica napus to illustrate the proposed method.

Comparison of flowering time genes in Brassica rapa, B. napus and Arabidopsis thaliana.

The major difference between annual and biennial cultivars of oilseed Brassica napus and B. rapa is conferred by genes controlling vernalization-responsive flowering time. These genes were compared between the species by aligning the map positions of flowering time quantitative trait loci (QTLs) detected in a segregating population of each species. The results suggest that two major QTLs identified in B. rapa correspond to two major QTLs identified in B. napus. Since B. rapa is one of the hypothesized diploid parents of the amphidiploid B. napus, the vernalization requirement of B. napus probably originated from B. rapa. Brassica genes also were compared to flowering time genes in Arabidopsis thaliana by mapping RFLP loci with the same probes in both B. napus and Arabidopsis. The region containing one pair of Brassica QTLs was collinear with the top of chromosome 5 in A. thaliana where flowering time genes FLC, FY and CO are located. The region containing the second pair of QTLs showed fractured collinearity with several regions of the Arabidopsis genome, including the top of chromosome 4 where FRI is located. Thus, these Brassica genes may correspond to two genes (FLC and FRI) that regulate flowering time in the latest flowering ecotypes of Arabidopsis.

Genetic Analysis and Identification of Amplified Fragment Length Polymorphism Markers Linked to the alm1 Avirulence Gene of Leptosphaeria maculans.

ABSTRACT A gene-for-gene interaction was previously suggested by mapping of a single major locus (LEM 1) controlling cotyledon resistance to Leptosphaeria maculans isolate PHW1245 in Brassica napus cv. Major. In this study, we obtained further evidence of a gene-for-gene interaction by studying the inheritance of the corresponding avirulence gene in L. maculans isolate PHW1245. The analysis of segregating F(1) progenies and 14 test crosses suggested that a single major gene is involved in the interaction. This putative avirulence gene was designated alm1 after the resistance locus identified in B. napus. Amplified fragment length polymorphism (AFLP) markers were used to generate a rudimentary genetic linkage map of the L. maculans genome and to locate markers linked to the putative avirulence locus. Two flanking AFLP markers, AC/TCC-1 and AC/CAG-5, were linked to alm1 at 3.1 and 8.1 cM, respectively. Identification of markers linked to the avirulence gene indicated that the differential interaction is controlled by a single gene difference between parental isolates and provides further support for the gene-for-gene relationship in the Leptosphaeria-Brassica system.

Evaluation of Sclerotinia Stem Rot Resistance in Oilseed Brassica napus Using a Petiole Inoculation Technique Under Greenhouse Conditions.

A petiole inoculation technique was adapted for evaluating resistance of oilseed Brassica napus seedlings to Sclerotinia sclerotiorum. In the first of four experiments, four isolates of S. sclerotiorum were tested, two originating from soybean and two from B. napus. In all, 10 to 47 B. napus accessions were inoculated in the seedling stage and responses to isolates were evaluated using days to wilt (DW) and a lesion phenotype index (LP). There were no significant differences in virulence among the four isolates for DW and only slight differences for LP. However, significant differences (P < 0.0001) were observed among the B. napus accessions for DW and LP in this experiment and in subsequent experiments using one isolate. The responses of accessions were consistent among experiments and among evaluation criteria. Higher levels of resistance were found among winter-type than spring-type accessions, and among rapeseed-quality compared with canola-quality accessions. The most resistant accessions identified also were the most resistant when inoculated at the flowering stage. Terminal stems were inoculated immediately below the lowest flower and stem lesion length (SLL) was used to characterize the interaction phenotype of each accession. The petiole inoculation technique can be used successfully to differentiate oilseed B. napus germ plasm for response to S. sclerotiorum. This inoculation technique and the sources of resistance identified in this study may be used to determine inheritance resistance to S. sclerotiorum and for improving oilseed B. napus cultivars for resistance to this important pathogen.

Analyses of a multi-parent population derived from two diverse alfalfa germplasms: testcross evaluations and phenotype-DNA associations.

In a previous study, we showed that the genetic variation present in the Medicago sativa subsp. sativa Peruvian and M. sativa subsp. falcata WISFAL germplasms could be used to improve forage yields when favorable alleles were recombined and used in hybrid combination with cultivated alfalfa. In this paper, we present testcross forage yield and fall growth data for two seasons of a C0 population generated after intermating the Peruvian x WISFAL population for several generations. In addition, we conducted marker-trait association analysis as an attempt to identify Peruvian and WISFAL genomics regions affecting the targeted traits. Five and seven genomic regions were found significantly associated with forage yield and fall growth, respectively. In the case of fall growth, alleles from both accessions were positively associated with plant height. However, more alleles from WISFAL were positively associated with forage yield than from Peruvian. WISFAL is known for its winter hardiness and genomic regions with large effects on winter survival may have masked the effect of forage yield from Peruvian. The fact that most of the genomic regions discovered in this study have been previously associated with traits involved in winter hardiness validates our findings and suggests that associations between DNA fragments and agronomic traits can be detected without the necessity of developing bi-parental mapping populations.

Mapping loci controlling flowering time in Brassica oleracea.

The timing of the transition from vegetative to reproductive phase is a major determinant of the morphology and value of Brassica oleracea crops. Quantitative trait loci (QTLs) controlling flowering time in B. oleracea were mapped using restriction fragment length polymorphism (RFLP) loci and flowering data of F3 families derived from a cabbage by broccoli cross. Plants were grown in the field, and a total of 15 surveys were made throughout the experiment at 5-15 day intervals, in which plants were inspected for the presence of flower buds or open flowers. The flowering traits used for data analysis were the proportion of annual plants (PF) within each F3 family at the end of the experiment, and a flowering-time index (FT) that combined both qualitative (annual/biennial) and quantitative (days to flowering) information. Two QTLs on different linkage groups were found associated with both PF and FT and one additional QTL was found associated only with FT. When combined in a multi-locus model, all three QTLs explained 54.1% of the phenotypic variation in FT. Epistasis was found between two genomic regions associated with FT. Comparisons of map positions of QTLs in B. oleracea with those in B. napus and B. rapa provided no evidence for conservation of genomic regions associated with flowering time between these species.

Mapping of RFLP and qualitative trait loci in Brassica rapa and comparison to the linkage maps of B. napus, B. oleracea, and Arabidopsis thaliana.

A linkage map of restriction fragment length polymorphisms (RFLPs) was constructed for oilseed, Brassica rapa, using anonymous genomic DNA and cDNA clones from Brassica and cloned genes from the crucifer Arabidopsis thaliana. We also mapped genes controlling the simply inherited traits, yellow seeds, low seed erucic acid, and pubescence. The map included 139 RFLP loci organized into ten linkage groups (LGs) and one small group covering 1785 cM. Each of the three traits mapped to a single locus on three different LGs. Many of the RFLP loci were detected with the same set of probes used to construct maps in the diploid B. oleracea and the amphidiploid B. napus. Comparisons of the linkage arrangements between the diploid species B. rapa and B. oleracea revealed six LGs with at least two loci in common. Nine of the B. rapa LGs had conserved linkage arrangements with B. napus LGs. The majority of loci in common were in the same order among the three species, although the distances between loci were largest on the B. rapa map. We also compared the genome organization between B. rapa and A. thaliana using RFLP loci detected with 12 cloned genes in the two species and found some evidence for a conservation of the linkage arrangements. This B. rapa map will be used to test for associations between segregation of RFLPs, detected by cloned genes of known function, and traits of interest.

Genetic diversity of oilseedBrassica napus germ plasm based on restriction fragment length polymorphisms.

Oilseed rape (Brassica napus) is an important oilseed crop worldwide. Cultivars have been developed for many growing regions, however little is known about genetic diversity inB. napus germ plasm. The purpose of the research presented here was to study the genetic diversity and relationships ofB. napus accessions using restriction fragment length polymorphisms (RFLPs). Eighty threeB. napus accessions were screened using 43 genomic DNA clones which revealed 161 polymorphic fragments. Each accession was uniquely identified by the markers with the exception of the near-isogenic cvs 'Triton' and 'Tower'. The RFLP data were analyzed by cluster analysis of similarity coefficients and by principal component analysis. Overall, there were three major groups of cultivars. The first group included only spring accessions, the second mostly winter accessions and the third, rutabagas and oilseed rape accessions from China and Japan. These results indicate that withinB. napus, winter and spring cultivars represent genetically distinct groups. The grouping of accessions by cluster analysis was generally consistent with known pedigrees. This consistency included the grouping of lines derived both by backcrossing or self-pollination with their parents.

Characterization of a dwarf gene in Brassica rapa, including the identification of a candidate gene.

Dwarf genes have been valuable for improving harvestable yield of several crop plants and may be useful in oilseed Brassica. We evaluated a dwarf gene, dwf2, from Brassica rapa in order to determine its phenotypic effects and genetic characteristics. The dwf2 mutant was insensitive to exogenous GA(3) for both plant height and flowering time, suggesting that it is not a mutation in the gibberellin biosynthesis pathway. The dwarf phenotype was controlled by a semidominant allele at a single locus. Near-isogenic lines that were homozygous or heterozygous for dwf2 had 47.4% or 30.0% reduction in plant height, respectively, compared to the tall wild-type line, and the reduction was due to reduced internode length and number of nodes. The dwf2 homozygous and heterozygous lines had the same or significantly higher numbers of primary branches than the wild-type line, but did not differ in flowering time. The DWF2 gene was mapped to the bottom of linkage group R6, in a region having homology to the top of Arabidopsis thaliana chromosome 2. The map position of DWF2 in comparison to markers in A. thaliana suggests it is a homolog of RGA ( repressor of ga1-3), which is a homolog of the wheat "Green Revolution" gene. This dwarf gene could be used to gain more insight on the gibberellin pathway and to reduce lodging problems in hybrid oilseed Brassica cultivars.

Novel flowering time variation in the resynthesized polyploid Brassica napus.

Recent molecular data using resynthesized polyploids of Brassica napus established that genome changes can occur rapidly after polyploid formation. In this study we present data that de novo phenotypic variation for flowering time also occurs rapidly after polyploidization. Two initial polyploid plants were developed by reciprocal crosses of B. rapa and B. oleracea followed by chromosome doubling to establish two lineages, each of which was expected to be homozygous and homogeneous. Several sublineages of each lineage were advanced by self-pollination. The range in days to flower of the sixth generation plants was 39-75 and 43-64 for the two lineages. Analysis of seventh generation progeny indicated that the variation was heritable. Lines were selected and self-pollinated to the eighth generation and also testcrossed to a natural B. napus cultivar; the testcross plants were then self-pollinated. Differences in flowering time were also inherited in these advanced generations. Days to flower was significantly correlated with leaf number in each generation. The rapid evolution of new phenotypic variation, like that observed in this model system, may have contributed to the success and diversification of natural polyploid organisms.

A method for examining expression of homologous genes in plant polyploids.

One of the essential issues regarding evolution of polyploid species is how duplicate genes are expressed. Most studies on gene expression in polyploids have been based on isozyme analyses; RNA analysis has not been widely used partially due to difficulties in distinguishing homologous transcripts which usually have the same length and similar or almost identical sequences. In this study, a method combining RT-PCR with RFLP was used to analyze transcripts of homologous genes in natural and synthetic Brassica amphidiploids. Sequences coding for several known genes were selected and used to synthesize gene-specific primers. Total RNAs were used as templates for RT-PCR to amplify homologous transcripts in three diploid parental species, three cultivated amphidiploid species and six synthetic amphidiploids. For each gene, initial PCR products amplified in all species had identical length; however, homologous transcripts in the diploid and amphidiploid species could be distinguished after digesting the PCR products with restriction enzymes. Preliminary results based on three genes indicated that both transcripts from the diploid parents were expressed in the synthetic and natural amphidiploids. This study represents the first application of RT-PCR and RFLP analysis to investigate expression of homologous genes in higher plants. The technique is a sensitive, simple and efficient method for distinguishing homologous transcripts in a mixed RNA population and can be applied to many types of studies on expression of homologous genes.

Variation among alfalfa somaclones in copy number of repeated DNA sequences.

Repeated DNA sequences of alfalfa (Medicago sativa L.) somaclonal variants were analyzed to determine if changes in copy number had occurred during tissue culture. DNA clones containing highly repeated nuclear sequences from the diploid line HG2 (2x = 16) were slot blotted and probed with labeled DNAs from HG2 and several somaclones of HG2. Two DNA clones that differed visually in hybridization intensity among the plant DNAs and one clone that had constant hybridization intensity were selected and used as probes on Southern blots and slot blots containing equal quantities of DNAs from HG2 and 15 régénérants. Statistically significant differences were detected in the copy number of two anonymous DNA sequences initially selected as variable and in the copy number of sequences homologous to pea ribosomal DNA. Based on Southern blot analysis, these sequences appeared to be arranged as tandem repeats. The cloned sequence initially selected as stable did not vary significantly in copy number and it appeared to be arranged as a dispersed repeat. Both increases and decreases in copy number of repeated sequences were observed in plants from successive regeneration cycles. Results from this study indicate that specific repeated nuclear DNA sequences have changed copy number in plants regenerated from tissue culture.

Mapping loci controlling vernalization requirement in Brassica rapa.

Brassica cultivars are classified as biennial or annual based on their requirement for a period of cold treatment (vernalization) to induce flowering. Genes controlling the vernalization requirement were identified in a Brassica rapa F2 population derived from a cross between an annual and a biennial oilseed cultivar by using an RFLP linkage map and quantitative trait locus (QTL) analysis of flowering time in F3 lines. Two genomic regions were strongly associated with variation for flowering time of unvernalized plants and alleles from the biennial parent in these regions delayed flowering. These QTLs had no significant effect on flowering time after plants were vernalized for 6 weeks, suggesting that they control flowering time through the requirement for vernalization. The two B. rapa linkage groups containing these QTLs had RFLP loci in common with two B. napus linkage groups that were shown previously to contain QTLs for flowering time. An RFLP locus detected by the cold-induced gene COR6.6 cloned from Arabidopsis thaliana mapped very near to one of the B. rapa QTLs for flowering time.

Purification and characterization of arcelin seed protein from common bean.

Arcelin, a seed protein originally discovered in wild bean accessions, was purified, characterized, and compared to phaseolin, the major seed protein of common bean, and to phytohemagglutinin (PHA), the major bean seed lectin. Arcelin and PHA has several characteristics in common. Both were glycoproteins having similar subunit M(r), deglycosylated M(r), and amino acid compositions. The two proteins were related antigenically and they had the same developmental timing of accumulation. Arcelin also had some hemagglutinating activity, a characteristic associated with lectins. However, several features distinguished arcelin from PHA. Arcelin had a more basic isoelectric point than PHA, greater numbers of basic amino acid residues, additional cysteine residues, and one methionine residue, which PHA lacks. Native PHA protein is a tetramer of subunits, and although a small component of native arcelin protein was also tetrameric, most of the arcelin preparation was dimeric. The hemagglutinating activity of arcelin was specific only for some pronase-treated erythrocytes. It did not agglutinate native erythrocytes, nor did it bind to thyroglobulin or fetuin affinity resins as did PHA. Although arcelin has lectin-like properties, we believe the distinctions between arcelin and PHA warrant the designation of arcelin as a unique bean seed protein.

Hordein gene expression in a low protein barley cultivar.

The low protein barley (Hordeum vulgare L.) cultivar ;Karl' is deficient in hordeins, the major storage protein fraction, but has normal levels of other endosperm protein fractions. We compared hordein gene expression in Karl with that in ;Traill,' a related cultivar that contains normal amounts of hordein. In plants grown in controlled environmental conditions, hordein accumulation in Karl was 55% of that in Traill. The amount of hordein synthesized at all developmental stages tested was lower in Karl, as shown by pulse labeling of excised spikes with [(3)H]leucine. To determine whether levels of hordein mRNAs were related to amounts of hordein synthesis, total RNA preparations from endosperms at three or four developmental stages were hybridized to cDNA probes corresponding to B and C hordein genes. Both B and C hordein mRNA levels were significantly lower in Karl than in Traill at intermediate developmental stages. It was concluded that the low-hordein character of Karl was regulated by the size of the hordein mRNA pool. Transcriptional and/or post-transcriptional controls may be involved in the regulation of hordein mRNA levels.

Bean lectins : 5. Quantitative genetic variation in seed lectins of Phaseolus vulgaris L. and its relationship to qualitative lectin variation.

Seeds of forty bean cultivars having different lectin types based on two-dimensional isoelectric focusing-sodium dodecyl sulfate polyacrylamide gel electrophoresis (IEF-SDS/PAGE) were analyzed for quantities of lectin, phaseolin and total protein. Significant differences were found among groups of cultivars with different lectin types for the quantity of lectin and phaseolin. Cultivars with more complex lectin types based on IEF-SDS/PAGE tended to have higher quantities of lectin and lower quantities of phaseolin than cultivars with simple lectin types. An association between lectin type and the quantity of lectin and phaseolin was found also in the seeds of F2 plants that segregated in a Mendelian fashion for two lectin types. Seeds from plants with the complex lectin type had more lectin and less phaseolin than seeds from plants with the simple lectin type. Therefore, the genes controlling qualitative lectin variation also may influence the quantitative variation of lectin and phaseolin. The results of this study are related to other studies on the quantitative variation for seed proteins and to the possible molecular basis for variation in the quantity of lectins in beans.

Characterization and comparison of arcelin seed protein variants from common bean.

Four variants of arcelin, an insecticidal seed storage protein of bean, Phaseolus vulgaris L., were investigated. Each variant (arcelin-1, -2, -3, and -4) was purified, and solubilities and M(r)s were determined. For arcelins-1, -2, and -4, the isoelectric points, hemagglutinating activities, immunological cross-reactivities, and N-terminal amino acid sequences were determined. On the basis of native and denatured M(r)s, the variants were classified as being composed of dimer protein (arcelin-2), tetramer protein (arcelins-3 and -4), or both dimer and tetramer proteins (arcelin-1). Although the dimer proteins (arcelins-1d and -2) could be distinguished by M(r)s and isoelectric points, they were identical for their first 37 N-terminal amino acids and had similar immunological cross-reactions, and bean lines containing these variants had a DNA restriction fragment in common. The tetramer proteins arcelin-1t and arcelin-4 also could be distinguished from each other based on M(r)s and isoelectric points; however, they had similar immunological cross-reactions and they were 77 to 93% identical for N-terminal amino acid composition. The similarities among arcelin variants, phytohemagglutinin, and a bean alpha-amylase inhibitor suggest that they are all encoded by related members of a lectin gene family.

Brassica taxonomy based on nuclear restriction fragment length polymorphisms (RFLPs) : 3. Genome relationships in Brassica and related genera and the origin of B. oleracea and B. rapa (syn. campestns).

RFLPs were used to study genome evolution and phylogeny in Brassica and related genera. Thirtyeight accessions, including 10 accessions of B. rapa (syn. campestris), 9 cultivated types of B. oleracea, 13 nine-chromosome wild brassicas related to B. oleracea, and 6 other species in Brassica and allied genera, were examined with more then 30 random genomic DNA probes, which identified RFLPs mapping to nine different linkage groups of the B. rapa genome. Based on the RFLP data, phylogenetic trees were constructed using the PAUP microcomputer program. Within B. rapa, accessions of pak choi, narinosa, and Chinese cabbage from East Asia constituted a group distinct from turnip and wild European populations, consistent with the hypothesis that B. rapa had two centers of domestication. A wild B. rapa accession from India was positioned in the tree between European types and East Asian types, suggesting an evolutionary pathway from Europe to India, then to South China. Cultivated B. oleracea morphotypes showed monophyletic origin with wild B. oleracea or B. alboglabra as possible ancestors. Various kales constitute a highly diverse group, and represent the primitive morphotypes of cultivated B. oleracea from which cabbage, broccoli, cauliflower, etc. probably have evolved. Cauliflower was found to be closely related to broccoli, whereas cabbage was closely related to leafy kales. A great diversity existed among the 13 collections of nine-chromosome wild brassicas related to B. oleracea, representing various taxonomic states from subspecies to species. Results from these studies suggested that two basic evolutionary pathways exist for the diploid species examined. One pathway gave rise to B. fruticulosa, B. nigra, and Sinapis arvensis, with B. adpressa or a close relative as the initial ancestor. Another pathway gave rise to B. oleracea and B. rapa, with Diplotaxis erucoides or a close relative as the initial ancestor. Raphanus sativus and Eruca sativus represented intermediate types between the two lineages, and might have been derived from introgression or hybridization between species belonging to different lineages. Molecular evidence for an ascending order of chromosome numbers in the evolution of Brassica and allied genera was obtained on the basis of RFLP data and phylogenetic analysis.

Rapid genome change in synthetic polyploids of Brassica and its implications for polyploid evolution.

Although the evolutionary success of polyploidy in higher plants has been widely recognized, there is virtually no information on how polyploid genomes have evolved after their formation. In this report, we used synthetic polyploids of Brassica as a model system to study genome evolution in the early generations after polyploidization. The initial polyploids we developed were completely homozygous, and thus, no nuclear genome changes were expected in self-fertilized progenies. However, extensive genome change was detected by 89 nuclear DNA clones used as probes. Most genome changes involved loss and/or gain of parental restriction fragments and appearance of novel fragments. Genome changes occurred in each generation from F2 to F5, and the frequency of change was associated with divergence of the diploid parental genomes. Genetic divergence among the derivatives of synthetic polyploids was evident from variation in genome composition and phenotypes. Directional genome changes, possibly influenced by cytoplasmic-nuclear interactions, were observed in one pair of reciprocal synthetics. Our results demonstrate that polyploid species can generate extensive genetic diversity in a short period of time. The occurrence and impact of this process in the evolution of natural polyploids is unknown, but it may have contributed to the success and diversification of many polyploid lineages in both plants and animals.