Genetic variation of the Toll-like receptors in a Swedish allergic rhinitis case population

BACKGROUND

Variation in the 10 toll-like receptor (TLR) genes has been significantly associated with allergic rhinitis (AR) in several candidate gene studies and three large genome-wide association studies. These have all investigated common variants, but no investigations for rare variants (MAF ≤ 1%) have been made in AR. The present study aims to describe the genetic variation of the promoter and coding sequences of the 10 TLR genes in 288 AR patients.

METHODS

Sanger sequencing and Ion Torrent next-generation sequencing was used to identify polymorphisms in a Swedish AR population and these were subsequently compared and evaluated using 1000Genomes and Exome Aggregation Consortium (ExAC) data.

RESULTS

The overall level of genetic variation was clearly different among the 10 TLR genes. The TLR10-TLR1-TLR6 locus was the most variable, while the TLR7-TLR8 locus was consistently showing a much lower level of variation. The AR patients had a total of 37 promoter polymorphisms with 14 rare (MAF ≤ 1%) and 14 AR-specific polymorphisms. These numbers were highly similar when comparing the AR and the European part of the 1000Genomes populations, with the exception of TLR10 where a significant (P = 0.00009) accumulation of polymorphisms were identified. The coding sequences had a total of 119 polymorphisms, 68 were rare and 43 were not present in the European part of the 1000Genomes population. Comparing the numbers of rare and AR-specific SNPs in the patients with the European part of the 1000Genomes population it was seen that the numbers were quite similar both for individual genes and for the sum of all 10 genes. However, TLR1, TLR5, TLR7 and TLR9 showed a significant excess of rare variants in the AR population when compared to the non-Finnish European part of ExAC. In particular the TLR1 S324* nonsense mutation was clearly overrepresented in the AR population.

CONCLUSIONS

Most TLR genes showed a similar level of variation between AR patients and public databases, but a significant excess of rare variants in AR patients were detected in TLR1, TLR5, TLR7, TLR9 and TLR10. This further emphasizes the frequently reproduced TLR10-TLR1-TLR6 locus as being involved in the pathogenesis of allergic rhinitis.

Characterization of genetic variation in TLR8 in relation to allergic rhinitis

BACKGROUND

A previous investigation of all 10 TLR genes for associations with allergic rhinitis (AR) detected a number of significant SNPs in the TLR8 locus. The associations indicated that an accumulation of rare variants could explain the signal. This study therefore searches for rare variants in the TLR8 region and also investigates the reproducibility of previous SNP associations.

METHODS

The TLR8 gene was resequenced in 288 AR patients from Malmö and the data were compared with publically available data. Seven previously AR-associated SNPs from TLR8 were analyzed for AR associations in 422 AR patients and 859 controls from the BAMSE cohort. The associations detected in present and previous studies were compared.

RESULTS

Sequencing detected 13 polymorphisms (three promotor and 10 coding) among 288 AR patients. Four of the coding polymorphisms were rare (MAF < 1%) and three of those were novel. Two coding polymorphisms were benign missense mutations and the rest were synonymous. Comparison with 1000Genomes and Exome Aggregation Consortium data revealed no accumulation of rare variants in the AR cases. The AR association tests made using the BAMSE cohort yielded five P-values <0.05. Tests of IgE levels yielded four significant SNP associations to birch pollen. Comparing results between different populations revealed opposing risk alleles, different gender effects, and response to different allergens in the different populations.

CONCLUSIONS

Rare variants in TLR8 are not associated with AR. Comparison of present and previous association studies reveals contradictory results for common variants. Thus, no associations exist between genetic variation in TLR8 and AR.

Replication of genomewide associations with allergic sensitization and allergic rhinitis

BACKGROUND

Three genomewide metastudies have recently reported associations with self-reported allergic rhinitis and allergic sensitization. The three studies together identified a set of 37 loci but showed low concordance. This study investigates the reproducibility of the detected single nucleotide polymorphism (SNP) associations in an extensively characterized longitudinal cohort, BAMSE.

METHODS

Phenotypic evaluation of allergic rhinitis (AR) and allergic sensitization was performed on 2153 children from BAMSE at 8 and 16 years of age. Allele frequencies of 39 SNPs were investigated for association with the exact allergic phenotypes of the metastudies. Odds ratios and false discovery rates were calculated, and the impact of asthma was evaluated. The cases were also evaluated for age at onset effects (≤ or >8 years of age).

RESULTS

Association tests of the 39 SNPs identified 12 SNPs with P-values < 0.05 and Q-values < 0.10. Two of the four loci (TLR6-TLR1 and HLA-DQA1-HLA-DQB1) identified in all three original studies were also identified in this study. Three SNPs located in the TLR6-TLR1 locus had the lowest P-values and Q-values < 0.1 when using a well-defined AR phenotype. Two loci showed significant age at onset effects, but the effect of asthma on the associations was very limited.

CONCLUSION

The TLR6-TLR1 locus is likely to have a central role in the development of allergic disease. The association between genetic variation in the SSTR1-MIPOL1 and TSLP-SLC25A46 loci and age at onset is the first report of age at onset effects in allergic rhinitis.

Origin of Swedish hemophilia B mutations

BACKGROUND

More than 1100 mutations that cause hemophilia B (HB) have been identified. At the same time, specific F9 mutations are present at high frequencies in certain populations, which raise questions about the origin of HB mutations.

OBJECTIVES

To describe the mutation spectrum of all HB families in Sweden and investigate if mutations appearing in several families are due to independent recurrent mutations (RMs) or to a common mutation event (i.e. are identical by descent (IBD)).

PATIENTS/METHODS

The registered Swedish HB population consists of patients from 86 families. Mutations were identified by resequencing and identical haplotypes were defined using 74 markers and a control population of 285 individuals. The ages of IBD mutations were estimated using ESTIAGE.

RESULTS

Out of 77 presumably unrelated patients with substitution mutations, 47 patients (61%) had mutations in common with other patients. Haplotyping of the 47 patients showed that 24 patients had IBD mutations (51%) with estimated ages of between two and 23 generations. A majority of these patients had mild disease. Eight of the 15 mutations observed in more than one family were C>T transitions in CpG sites and all eight were RMs.

CONCLUSIONS

The association of IBD mutations with a mild phenotype is similar to what has been previously observed in hemophilia A. Noteworthy features of the mutations that are common to more than one family are the equal proportions of patients with RM and IBD mutations and the correlation between the occurrence of RMs and C>T transitions at CpG sites.

Origin of Swedish hemophilia A mutations

BACKGROUND

Hemophilia A (HA) has a high level of variation within the disease class, with more than 1000 mutations being listed in the HAMSTeRS database. At the same time a number of F8 mutations are present in specific populations at high frequencies.

OBJECTIVES

The simultaneous presence of large numbers of rare mutations and a small number of high-frequency mutations raises questions about the origins of HA mutations. The present study was aimed at describing the origins of HA mutations in the complete Swedish population. The primary issue was to determine what proportion of identical mutations are identical by descent (IBD) and what proportion are attributable to recurrent mutation events. The age of IBD mutations was also determined.

PATIENTS/METHODS

In Sweden, the care of HA is centralized, and the Swedish HA population consists of ~ 750 patients from > 300 families (35% severe, 15% moderate, and 50% mild). Identical haplotypes were defined by single-nucleotide polymorphism and microsatellite haplotyping, and the ages of the mutations were estimated with estiage.

RESULTS

Among 212 presumably unrelated patients with substitution mutations, 97 (46%) had mutations in common with other patients. Haplotyping of the 97 patients showed that 47 had IBD mutations (22%) with estimated ages of between two and 35 generations. The frequency of mild disease increased with an increasing number of patients sharing the mutations.

CONCLUSIONS

A majority of the IBD mutations are mild and have age estimates of a few hundred years, but some could date back to the Middle Ages.

A high density RFLP linkage map of sugar beet

A high density sugar beet RFLP map with an average distance of 1.5 cM between markers has been constructed. The map covers 621 cM and includes 413 markers distributed over the nine linkage groups of sugar beet. The map is based on two F2 populations representing two different pairs of parents. The two sets of data were integrated into a single map using 90 markers that were common to both data sets. The quality of the map was assessed in several ways. The common markers were used to investigate how often the loci had been mapped in the same order in the two F2 populations. For closely situated markers (<1.5 cM) the order specified in the map is uncertain, but for markers separated by more than 2 cM the locus order is highly reliable. The error rate of the overall process was estimated at 0.3% by independently repeating the analysis of 41 markers. The map is comparatively short, with a map length corresponding to approximately 1.4 crossovers per bivalent. Another feature of the map is a high degree of clustering of markers along the linkage groups. With the possible exception of linkage group 2, each linkage group shows one major cluster, which in most cases is situated in the centre of the linkage group. Our interpretation is that sugar beet, in comparison with most other species, has an extreme localization of recombination. Key words : sugar beet, linkage, RFLP, clustering.

Estimating genetic variation in sugar beets and wild beets using pools of individuals

The study describes the genetic structure in sugar beets and in wild beets (Beta vulgaris) using 30 RFLP markers. Samples consisting of pooled plant material of 100 individuals from each line and population were used to analyse 120 sugar beet breeding lines and 91 wild beet populations. Greater variation was found among the wild populations than among the breeding lines. Although the two major groups of breeding lines, monogerm and multigerm, had approximately equal amounts of genetic variation, in the monogerm group more of this variation was partitioned among the lines than within the lines. Furthermore, despite most of the variation being shared by the two groups, the two groups were found to be separated along the first two components in a principal component analysis. Computer simulations were carried out to evaluate the usefulness of the pooled-sample strategy employed in the investigation. These simulations showed the use of pooled samples to be a better alternative than that of analysing a few plants individually.

CLC- a novel susceptibility gene for allergic rhinitis?

BACKGROUND

Studies of the nasal lavage fluid proteome have previously identified proteins differently expressed in patients with symptomatic allergic rhinitis, e.g. S100A7, prolactin-inducible protein (PIP), wingless-type MMTV integration site family, member 2B (WNT2B), Charcot-Leyden crystal protein (CLC) and palate lung nasal epithelial clone (PLUNC). The aim of the present study was to investigate if genetic variation associated with allergic rhinitis can be found in these genes.

METHODS

Peripheral blood was collected from 251 patients with birch and/or grass pollen-induced allergic rhinitis and 386 nonatopic healthy controls. A total of 39 single nucleotide polymorphisms (SNPs) distributed over the genes PIP, WNT2B, CLC and PLUNC were selected from dbSNP, genotyped and investigated for associations with allergic rhinitis. Twelve additional SNPs were subsequently analysed for CLC.

RESULTS

All 22 investigated SNPs in CLC were polymorphic. Ten SNPs yielded significant differences between cases and controls with respect to genotype frequencies. Homozygotes for the minor allele were more common in allergic individuals compared to healthy controls. The minor alleles of these SNPs were all located on the same haplotype. Furthermore, homozygotes for the minor allele of two of the promoter SNPs had higher average scores for birch in skin prick test. In contrast, for seven SNPs within the gene, heterozygotes and homozygotes for the major allele had higher average scores for grass. None of the other three genes showed association.

CONCLUSION

Genetic variation in CLC was found to be associated with allergic rhinitis. The pattern of variation is compatible with a recessive inheritance model and the previously observed altered protein levels detected in patients with allergic rhinitis.

A haplotype in the inducible T-cell tyrosine kinase is a risk factor for seasonal allergic rhinitis

BACKGROUND

Identification of disease-associated single nucleotide polymorphisms (SNPs) in seasonal allergic rhinitis (SAR) may be facilitated by focusing on genes in a disease-associated pathway.

OBJECTIVE

To search for SNPs in genes that belong to the T-cell receptor (TCR) pathway and that change in expression in allergen-challenged CD4+ cells from patients with SAR.

METHODS

CD4+ cells from patients with SAR were analysed with gene expression microarrays. Allele, genotype and haplotype frequencies were compared in 251 patients and 386 healthy controls.

RESULTS

Gene expression microarray analysis of allergen-challenged CD4+ cells from patients with SAR showed that 25 of 38 TCR pathway genes were differentially expressed. A total of 62 SNPs were analysed in eight of the 25 genes; ICOS, IL4, IL5, IL13, CSF2, CTLA4, the inducible T-cell tyrosine kinase (ITK) and CD3D. Significant chi-squared values were identified for several markers in the ITK kinase gene region. A total of five SNPs were nominally significant at the 5% level. Haplotype analysis of the five significant SNPs showed increased frequency of a haplotype that covered most of the coding part of ITK. The functional relevance of ITK was supported by analysis of an independent material, which showed increased expression of ITK in allergen-challenged CD4+ cells from patients, but not from controls.

CONCLUSION

Analysis of SNPs in TCR pathway genes revealed that a haplotype that covers a major part of the coding sequence of ITK is a risk factor for SAR.

The evolutionary history of the common chloroplast genome of Arabidopsis thaliana and A. suecica

The evolutionary history of the common chloroplast (cp) genome of the allotetraploid Arabidopsis suecica and its maternal parent A. thaliana was investigated by sequencing 50 fragments of cpDNA, resulting in 98 polymorphic sites. The variation in the A. suecica sample was small, in contrast to that of the A. thaliana sample. The time to the most recent common ancestor (T(MRCA)) of the A. suecica cp genome alone was estimated to be about one 37th of the T(MRCA) of both the A. thaliana and A. suecica cp genomes. This corresponds to A. suecica having a MRCA between 10 000 and 50 000 years ago, suggesting that the entire species originated during, or before, this period of time, although the estimates are sensitive to assumptions made about population size and mutation rate. The data was also consistent with the hypothesis of A. suecica being of single origin. Isolation-by-distance and population structure in A. thaliana depended upon the geographical scale analysed; isolation-by-distance was found to be weak on the global scale but locally pronounced. Within the genealogical cp tree of A. thaliana, there were indications that the root of the A. suecica species is located among accessions of A. thaliana that come primarily from central Europe. Selective neutrality of the cp genome could not be rejected, despite the fact that it contains several completely linked protein-coding genes.

Chloroplast DNA indicates a single origin of the allotetraploid Arabidopsis suecica

DNA sequencing was performed on up to 12 chloroplast DNA regions [giving a total of 4288 base pairs (bp) in length] from the allopolyploid Arabidopsis suecica (48 accessions) and its two parental species, A. thaliana (25 accessions) and A. arenosa (seven accessions). Arabidopsis suecica was identical to A. thaliana at all 93 sites where A. thaliana and A. arenosa differed, thus showing that A. thaliana is the maternal parent of A. suecica. Under the assumption that A. thaliana and A. arenosa separated 5 million years ago, we estimated a substitution rate of 2.9 x 10(-9) per site per year in noncoding single copy sequence. Within A. thaliana we found 12 substitution (single bp) and eight insertion/deletion (indel) polymorphisms, separating the 25 accessions into 15 haplotypes. Eight of the A. thaliana accessions from central Sweden formed one cluster, which was separated from a cluster consisting of central European and extreme southern Swedish accessions. This latter cluster also included the A. suecica accessions, which were all identical except for one 5 bp indel. We interpret this low level of variation as a strong indication that A. suecica effectively has a single origin, which we dated at 20 000 years ago or more.

Genetic variation in Arabidopsis suecica and its parental species A. arenosa and A. thaliana

Random amplified polymorphic DNA (RAPD) markers were used to estimate the level of genetic variation in Swedish accessions of the allopolyploid Arabidopsis suecica and its parental species A. thaliana and A. arenosa. The results showed clear differences among the three species with respect to the level of variation. A. arenosa was highly variable, A. thaliana showed a moderate level of variation whereas A. suecica was much less variable than the two other species. An extended analysis covering 19 Swedish populations of A. suecica corroborated the low level of variation in this species, yet 16 unique phenotypes were observed. No isolation by distance was observed. When the genetic variation was partitioned among and within populations of A. suecica, the results showed that the majority of the variation (81%) occurred among populations. This result is interpreted as a strong indication that A. suecica is autogamous in nature.

Identification of cytogenetic subgroups and karyotypic pathways in transitional cell carcinoma

The clinical course in urinary bladder cancer is difficult or impossible to predict based on conventional disease parameters. It is a reasonable hypothesis that the genetic aberrations acquired by the tumor cells, being instrumental in bringing about the disease in the first place, may also hold the key to more reliable prognostication. However, though 200 transitional cell carcinomas (TCC), the most common bladder cancer in the Western world, with clonal chromosomal abnormalities have been reported, our knowledge about the karyotypic characteristics of these tumors remains insufficient. The aberration pattern is clearly nonrandom, but no completely specific primary or secondary karyotypic abnormality has been identified, and the chronological order in which the aberrations appear during disease progression is not well known. The high degree of karyotypic complexity in epithelial tumors like TCC is one reason why our picture of the sequential order of cytogenetic evolution is unclear. To overcome some of these difficulties we have used several statistical methods that allow analysis and interpretation of the relationship between cytogenetic aberrations in TCC. We show that there exists a temporal order with respect to the appearance of chromosomal imbalances and that this order is highly correlated with tumor stage and grade. Analyzing changes in the distribution of imbalances per tumor in G1, G2, and G3 tumors, we suggest that progression involves the acquisition of cytogenetically detectable and submicroscopic genetic changes at comparable frequencies. By means of computer simulations, we show that the imbalances -9, +7, and 1q+ appear earlier than expected from random events and that -6q, -5q, -18, +5p, -22p, and -15 appear later than expected. Using principal component analysis, we identify two cytogenetic pathways in TCC, one initiated by -9 and followed by -11p and 1q+, the other initiated by +7 and followed by 8p- and +8q. The -9 pathway was correlated with stage Ta-T2 tumors, whereas the +7 pathway was correlated with stage T1-T3 tumors, i.e., +7 tumors appeared to be more aggressive. Although these pathways are well separated at earlier stages, they later converge to contain a common set of imbalances.

The incidence of trisomy 8 as a sole chromosomal aberration in myeloid malignancies varies in relation to gender, age, prior iatrogenic genotoxic exposure, and morphology

Although trisomy 8 as a sole change is one of the most common chromosomal abnormalities in myeloid malignancies, it is largely unknown if the incidence of this aberration is influenced by other factors of clinical importance. In the present study, the frequencies of isolated +8 in relation to gender, age, previous treatment with chemo- or radiotherapy, and morphologic subtype were ascertained in published, as well as in our own unpublished, cases of acute myeloid leukemia (AML; n=4,246), myelodysplastic syndromes (MDS; n=1,817), and chronic myeloproliferative disorders (MPD; n=530). The frequencies of +8 were higher in MDS and MPD than in AML (7.5% vs. 5.6%; P<0.01) and varied among the morphologic subtypes of AML and MDS (P<0.001 and P<0.05, respectively). Trisomy 8 was more common in women than in men with MPD (11% vs. 5.1%; P<0.05). Furthermore, the frequencies of +8 were higher in de novo AML and MDS than in treatment-related cases (6.0% vs. 2.8%; P<0.01 and 8.6% vs. 1.5%; P<0.001, respectively). The incidence also varied significantly with age in AML (P<0.001), being more common in elderly patients. Although the causes for this frequency heterogeneity remain to be elucidated, possible explanations may include different environmental exposures affecting the origin of +8 in AML, MDS, and MPD and the presence of different underlying cryptic primary aberrations.

Multivariate analyses of genomic imbalances in solid tumors reveal distinct and converging pathways of karyotypic evolution

A total of 3,016 malignant solid tumors (kidney, colorectal, breast, head and neck, ovarian, and lung carcinomas, neuroglial tumors, malignant melanoma, and testicular germ cell tumors) were selected for statistical analyses regarding karyotypic evolution. Genomic imbalances, i.e., net gains and losses, present in more than 5% of each tumor type were identified. Individual tumors were then classified with respect to absence or presence of these imbalances. To analyze for possible patterns of correlated imbalances, principal component analyses (PCA) were performed. Furthermore, algorithms were developed to analyze the temporal order of the imbalances, as well as the possible selection for early or late appearance in the karyotypic evolution. By analyzing the temporal order of imbalances common to many tumor types, a general order for nine of these emerged, namely, +7, -3p, -6q, -1p, -8p, -17p, -9p, -18, and -22. The distributions of the number of imbalances per case revealed a geometrical distribution, ranging from one to nine imbalances per tumor, in the majority of the tumor types. In tumor types in which cases with a high number of imbalances per case were frequent, notably head and neck, ovarian, and lung carcinomas, the overall distributions were bimodal, indicating the presence of two modes of chromosome evolution. By combining data from the PCA with the temporal analyses, it was possible to identify karyotypic pathways. It was found that the majority of the tumor types displayed more than one cytogenetic route, but, as the karyotypic evolution continued, these converged to a common pathway.

Primer mixtures in RAPD analysis

RAPD (random amplified polymorphic DNA) is a multiplex marker system that conventionally uses single-primer PCR to amplify random DNA fragments. Because of its multiplex nature, it is frequently used in bulked segregant analysis (BSA). In view of the very large numbers of markers BSA often requires, we investigated the use of mixtures of primers as a method of increasing the number of markers available. Theoretically, if a single-primer reaction produces x bands on average, an unrestrained PCR process using a primers should produce xa2 bands. Initially, we investigated mixtures containing from one to five primers. The average number of products increased slightly from the single-primer to the multiple-primer case, whereas it was rather constant for the different multi-primer combinations. This deviation from the theoretical expectations, which we attribute to the effects of competition, shows mixtures of more than two primers to be inefficient. The properties of two-primer mixtures in which the proportions of the two primers were varied were also investigated. The intensities of most of the products were influenced by the proportions of the primers used to create the mixture. A good fit was obtained to a model in which the average competitive ability of a band is directly proportional to the probability of randomly obtaining the band-producing primer combination from the pool of primers. Using two-primer mixtures, a(a-1)/2 different two-primer mixtures can be produced. A comparison of different schemes for constructing the two-primer mixtures indicates that the degree of resampling is similar for all schemes. In conclusion, the use of two-primer mixtures is a simple but very powerful strategy in BSA as it can generate an extremely large number of markers.

Linkage disequilibrium mapping of the bolting gene in sea beet using AFLP markers

The possibility of using linkage disequilibrium mapping in natural plant populations was assessed. In studying linkage disequilibrium among 137 mapped AFLP markers in four populations of sea beet (Beta vulgaris ssp. maritima (L.) Arcang.) it was shown that tightly linked loci could be detected by screening for associations. It was hypothesized that the short distances spanned by linkage disequilibrium enable markers that are very tightly linked to a target gene to be identified. The hypothesis was tested by whole-genome screening of AFLP markers for association with the gene for the annual growth habit, the B gene, in a sample of 106 sea beets. Despite the dominant nature of AFLP, two markers showing significant linkage disequilibrium with the B gene were detected. The results indicate the potential use of linkage disequilibrium for gene mapping in natural plant populations.

An evaluation of the use of pooled samples in studies of genetic variation

When using molecular markers to study genetic variation, either the sampled individuals can be analysed individually or the individuals can be pooled and only the pools analysed (pooled samples). A theoretical investigation was carried out into the use of pooled samples in the detection of alleles and providing maximum likelihood estimates of allele frequency. The results show that, in many cases, pooled samples are more efficient than samples of individuals. Of the different pool sizes studied, small pools containing two or three individuals showed the smallest expected squared error of allele frequency estimates.

Positive correlation between recombination rates and levels of genetic variation in natural populations of sea beet (Beta vulgaris subsp. maritima)

The relation between the level of genetic variation and the rate of recombination per physical unit was investigated in sea beet (Beta vulgaris subsp. maritima). The rate of recombination per physical unit was estimated indirectly through marker density in an RFLP linkage map of sugar beet. From this map, we also selected RFLP markers covering two of the nine chromosomes in Beta. The markers were used to estimate the level of genetic variation in three populations of sea beet, two from Italy and one from England. Two estimates of genetic variation were employed, one based on the number of alleles in the sample and the other on heterozygosity. A statistically significant positive correlation was found between recombination rate and genetic variation. Several theoretical explanations for this are discussed, background selection being one. A correlation similar to this has been observed previously in Drosophila, one that was higher than what we obtained for Beta. This is consistent with various biological differences between the two species.

Comparing the distribution of RAPD and RFLP markers in a high density linkage map of sugar beet

The distribution of RAPD markers was compared with that of RFLP markers in a high density linkage map of sugar beet. The same mapping population of 161 F2individuals was used to generate all the marker data. The total map comprises 160 RAPD and 248 RFLP markers covering 508 cM. Both the RAPD and the RFLP markers show a high degree of clustering over the nine linkage groups. The pattern is compatible with a strong distal localization of recombination in the sugar beet. It leads generally to one major cluster of markers in the centre of each linkage group. In regions of high marker density, dominant RAPD markers present in either linkage phase and codominant RFLP markers are subclustered relative to each other. This phenomenon is shown to be attributable to: (i) effects of the mapping procedure when dominant and codominant data are combined, (ii) effects of the mapping procedure when dominant data in both linkage phases are combined, and (iii) genuine differences in the way RAPD and RFLP markers are recruited.Key words: sugar beet, linkage map, RAPD, RFLP, clustering.

Competition as a source of errors in RAPD analysis

We have used artificial 1∶1 DNA mixtures of all pairwise combinations of four doubled haploid Brassica napus lines to test the ability of RAPDs to function as reliable dominant genetic markers. In situations where a specific RAPD band is present in one homozygous line but absent in the other, the band is expected in the artificial heterozygote, i.e. in the 1∶1 DNA mixture. In 84 of all 613 heterozygous situations analysed, the expected band failed to amplify in the RAPD reaction. Thus, RAPD markers will lead to an erroneous genetic interpretation in 14% of all cases. In contrast, the formation of non-parental heteroduplex bands was found at a frequency of only 0.2%. Analysis of 1∶ 1 mixtures using (1) a different set of optimized reaction conditions and (2) a material with low genomic complexity (Bacillus cereus) gave identical results. Serial dilutions of one genome into another, in steps of 10%, showed that all of the polymorphic bands decreased in intensity as a linear function of their respective proportion in the mixture. In dilutions with water no differences in band intensity were detected. Thus, competition occurs in the amplification of all RAPD fragments and is a major source of genotyping errors in RAPD analysis.

A model of chiasma reduction of closely formed crossovers

A model of chiasma reduction is developed, evaluated by computer simulation and discussed in relation to the evolution of interference. The model assumes that adjacently formed crossovers can interact, if there is incomplete sister chromatid cohesion between them, and give rise to a reduced number of chiasmata. In the absence of crossover position interference this leads to a considerable risk of non-disjunction for an average sized chromosome. It is suggested that an important function of crossover interference is to reduce this effect. The question is raised whether chiasma reduction takes place today. Some available cytological data can be interpreted as showing that chiasma reduction occurs in normal meiosis. The possibility of chiasma reduction therefore needs further attention.

Phylogenetic analysis of the genus Hordeum using repetitive DNA sequences

A set of six cloned barley (Hordeum vulgare) repetitive DNA sequences was used for the analysis of phylogenetic relationships among 31 species (46 taxa) of the genus Hordeum, using molecular hybridization techniques. in situ hybridization experiments showed dispersed organization of the sequences over all chromosomes of H. vulgare and the wild barley species H. bulbosum, H. marinum and H. murinum. Southern blot hybridization revealed different levels of polymorphism among barley species and the RFLP data were used to generate a phylogenetic tree for the genus Hordeum. Our data are in a good agreement with the classification system which suggests the division of the genus into four major groups, containing the genomes I, X, Y, and H. However, our investigation also supports previous molecular studies of barley species where the unique position of H. bulbosum has been pointed out. In our experiments, H. bulbosum generally had hybridization patterns different from those of H. vulgare, although both carry the I genome. Based on our results we present a hypothesis concerning the possible origin and phylogeny of the polyploid barley species H. secalinum, H. depressum and the H. brachyantherum complex.

Crossover distribution in barley analysed through RFLP linkage data

The pattern of recombination in barley with regard to (1) the distribution of crossover points among whole gametes, (2) the distribution of crossover points among individual chromosomes and (3) the distribution of crossover points within chromosomes has been analysed using data sets underlying two recently published restriction fragment length polymorphism (RFLP) linkage maps representing male and female meiosis, respectively. The data indicated that the process of recombination had been random with no interference. The two data sets gave similar results, indicating that male and female meiosis in barley do no differ significantly. The possibility of using RFLP data in studies of crossover distribution is stressed.

The robustness of recombination frequency estimates in intercrosses with dominant markers

The robustness of the maximum likelihood estimates of recombination frequencies has been investigated in double intercrosses with complete dominance at both loci. The robustness was investigated with respect to bias in the recombination frequency estimates due to: (1) limited sample sizes, (2) heterogeneity in recombination frequencies between sexes or among meioses and (3) factors that distort the segregation-misclassification or differential viability. In the coupling phase, the recombination frequency estimates are quite robust with respect to most of the investigated factors. Potentially, the most serious cause of a bias is misclassifications, which tend to increase the recombination frequency estimates. In the repulsion phase, misclassifications are particularly serious, leading to extreme discrepancies between true and observed values. In addition, limited sample size and sex differences in recombination can also bias recombination frequency estimates in repulsion. These effects may pose serious problem in genetic mapping with random amplified polymorphic DNA (RAPD) markers.

Evaluation of RFLP and RAPD markers in a comparison of Brassica napus breeding lines

RFLP and RAPD markers were evaluated and compared for their ability to determine genetic relationships in a set of three B. napus breeding lines. Using a total of 50 RFLP and 92 RAPD markers, the relatedness between the lines was determined. In total, the RFLP and the RAPD analysis revealed more than 500 and 400 bands, respectively. The relative frequencies of loci with allele differences were estimated from the band data. The RFLP and RAPD marker sets detected very similar relationships among the three lines, consistent with known pedigree data. Bootstrap analyses showed that the use of approximately 30 probes or primers would have been sufficient to achieve these relationships. This indicates that RAPD markers have the same resolving power as RFLP markers when used on exactly the same set of B. napus genotypes. Since RAPD markers are easier and quicker to use, these markers may be preferred in applications where the relationships between closely-related breeding lines are of interest. The use of RAPD markers in fingerprinting applications may, however, not be warranted, and this is discussed in relation to the reliability of RAPD markers.

Chiasma and recombination data in plants: are they compatible?

It has long been considered that the number of chiasmata formed during meiosis corresponds to the number of crossovers indicated by the genetic map. However, recent investigations in plants show an unexpected discrepancy in the results obtained when calculating the total number of crossover events per meiosis by these two methods. Is this discrepancy due to methodological difficulties? Or is there something fundamentally wrong with our understanding of crossovers and chiasmata?

Cytoplasmic male sterility in beta is associated with structural rearrangements of the mitochondrial DNA and is not due to interspecific organelle transfer

Chloroplast (ct) and mitochondrial (mt) DNAs from four cytoplasmic male sterile (cms) and 22 normal fertile sugar beet lines and accessions of wild beets from the genus Beta have been compared with restriction analyses and Southern hybridizations. We have used restriction analyses of ctDNA as a phylogenetic marker to confirm the taxonomic relationships between the different cytoplasms. According to the ctDNA data, all four cms cytoplasms belong to the same taxonomic section, Beta. Restriction patterns of ct and mtDNA from fertile accessions produced analogous trees of similarity and showed a close correlation between the organellar DNA diversity and the accepted taxonomic classification of the species studied. However, the mt-DNA restriction profiles of the four cms types differed dramatically from each other and from those of all fertile accessions from the genus. No indication of cytoplasmic introgression was found in any of the four investigated cms types. Southern hybridization to mtDNA revealed variant genomic arrangements in the different fertile and cms cytoplasms, indicating that rearrangement of the mitochondrial genome is a common denominator to the different cms systems in Beta. It may, indeed, be a common property to spontaneously occurring cms in all or most species.

Apparent negative interference due to variation in recombination frequencies

Variation in recombination frequencies may lead to a bias in the estimated interference value in a linkage experiment. Depending on the pattern of variation, the bias may be toward negative interference or toward positive interference, even when there is positive interference at the cytological level. In this paper we have mainly concentrated on the case of negative interference. We use models to quantify this effect when data are derived from a backcross experiment or from the selfing of F1 individuals. The effect is quantitatively similar in the two cases. There is an upper limit to the size the bias may reach for every given level of recombination. Two reported cases of negative interference in Drosophila and cultivated barley fall within this possible parameter range, i.e., the observed negative interference values could--at least in principle--be due solely to a variation in the recombination frequencies in the experiments.