Drosophila melanogaster resistance against the parasitoid wasp Leptopilina boulardi is under the control of a single gene (Rlb), with two alleles, the resistant one being dominant. Using strains bearing deletions, we previously demonstrated that the 55E2–E6; 55F3 region on chromosome 2R is involved in the resistance phenomenon. In this paper, we first restricted the Rlb containing region by mapping at the molecular level the breakpoints of the Df(2R)Pc66, Df(2R)P34 and Df(2R)Pc4 deficiencies, using both chromosomal in situ hybridization and Southern analyses. The resistance gene was localized in a 100 kb fragment, predicted to contain about 10 different genes. Male recombination genetic experiments were then performed, leading to identification of two possible candidates for the Rlb gene. Potential involvement of one of this genes, edl/mae, is discussed.

Deng et al. have recently proposed that estimates of an upper limit to the rate of spontaneous mutations and their average heterozygous effect can be obtained from the mean and variance of a given fitness trait in naturally segregating populations, provided that allele frequencies are maintained at the balance between mutation and selection. Using simulations they show that this estimation method generally has little bias and is very robust to violations of the mutation–selection balance assumption. Here I show that the particular parameters and models used in these simulations generally reduce the amount of bias that can occur with this estimation method. In particular, the assumption of a large mutation rate in the simulations always implies a low bias of estimates. In addition, the specific model of overdominance used to check the violation of the mutation–selection balance assumption is such that there is not a dramatic decline in mean fitness from overdominant mutations, again implying a low bias of estimates. The assumption of lower mutation rates and/or other models of balancing selection may imply considerably larger biases of the estimates, making the reliability of the proposed method highly questionable.

Strong sexual isolation exists between the closely related species Drosophila ananassae and D. pallidosa, but there is no obvious post-mating isolation; both sexes of the hybrids and their descendants appear to be completely viable and fertile. Strains exhibiting parthenogenesis have been derived from wild populations of both species. We intercrossed such strains and established iso-female lines after the second generation of parthenogenesis. These lines are clones, carrying homozygous chromosomes that are interspecific recombinants. We established 266 such isogenic lines and determined their genetic constitution by using chromosomal and molecular markers. Strong pseudo-linkage was seen between loci on the left arm of chromosome 2 and on the right arm of chromosome 3; the frequency of inheriting the two chromosome regions from the same species was significantly larger than expected. One possible cause of pseudo-linkage is female meiotic bias, so that chromosomes of the same species origin tend to be distributed to the same gamete. But this possibility is ruled out; backcross analysis indicated that the two chromosome regions segregated independently in female hybrids. The remaining possibility is elimination of low-fitness flies carrying the two chromosome regions from different species. Thus, genetic incompatibility was detected in the species pair for which no hybrid breakdown had previously been indicated. The ‘interspecific mosaic genome’ lines reported here will be useful for future research to identify genes involved in speciation and phenotypic evolution.

This study is aimed at improving the analysis of data used in identifying marker-associated effects on quantitative traits, specifically to account for possible departures from a Gaussian distribution of the trait data and to allow for asymmetry of marker effects attributable to phenotypic divergence between parental lines. A Bayesian procedure for analysing marker effects at the whole-genome level is presented. The procedure adopts a skewed t-distribution as a prior distribution of marker effects. The model with the skewed t-process includes Gaussian prior distributions, skewed Gaussian prior distributions and symmetric t-distributions as special cases. A Markov Chain Monte Carlo algorithm for obtaining marginal posterior distributions of the unknowns is also presented. The method was applied to a dataset on three traits (live weight, carcass length and backfat depth) measured in an F2 cross between Iberian and Landrace pigs. The distribution of marker effects was clearly asymmetric for carcass length and backfat depth, whereas it was symmetric for live weight. The t-distribution seems more appropriate for describing the distribution of marker effects on backfat depth.

The Drosophila nasuta subgroup of the immigrans species group is widely distributed throughout the South-East Asian region, consisting of morphologically similar species with varying degrees of reproductive isolation. Here, I report nucleotide variability data for five X-linked and two mtDNA loci in eight taxa from the nasuta subgroup, with deeper sampling from D. albomicans and its sister species D. nasuta. Phylogenetic relationships among these species vary among different genomic regions, and levels of genetic differentiation suggest that this species group diversified only about one million years ago. D. albomicans and D. nasuta share nucleotide polymorphisms and are distinguished by relatively few fixed differences. Patterns of genetic differentiation between this species pair are compatible with a simple isolation model with no gene flow. Nucleotide variability levels of species in the nasuta group are comparable to those in members of the melanogaster and pseudoobscura species groups, indicating effective population sizes on the order of several million. Population genetic analyses reveal that summaries of the frequency distribution of neutral polymorphisms in both D. albomicans and D. nasuta generally fit the assumptions of the standard neutral model. D. albomicans is of particular interest for evolutionary studies because of its recently formed neo-sex chromosomes, and our phylogenetic and population genetic analyses suggest that it might be an ideal model to study the very early stages of Y chromosome evolution.

Chromosome segment substitution (CSS) lines have the potential for use in QTL fine mapping and map-based cloning. The standard t-test used in the idealized case that each CSS line has a single segment from the donor parent is not suitable for non-idealized CSS lines carrying several substituted segments from the donor parent. In this study, we present a likelihood ratio test based on stepwise regression (RSTEP-LRT) that can be used for QTL mapping in a population consisting of non-idealized CSS lines. Stepwise regression is used to select the most important segments for the trait of interest, and the likelihood ratio test is used to calculate the LOD score of each chromosome segment. This method is statistically equivalent to the standard t-test with idealized CSS lines. To further improve the power of QTL mapping, a method is proposed to decrease multicollinearity among markers (or chromosome segments). QTL mapping with an example CSS population in rice consisting of 65 non-idealized CSS lines and 82 chromosome segments indicated that a total of 18 segments on eight of the 12 rice chromosomes harboured QTLs affecting grain length under the LOD threshold of 2·5. Three major stable QTLs were detected in all eight environments. Some minor QTLs were not detected in all environments, but they could increase or decrease the grain length constantly. These minor genes are also useful in marker-assisted gene pyramiding.

The present work provides the first broad-scale screening of allozymes in the land snail Helix aspersa. By using overall information available on the distribution of genetic variation between 102 populations previously investigated, we expect to strengthen our knowledge on the spread of the invasive aspersa subspecies in the Western Mediterranean. We propose a new approach based on a centre-based clustering procedure to cluster populations into groups following rules of geographical proximity and genetic similarity. Assuming a stepping-stone model of diffusion, we apply a partitioning algorithm which clusters only populations that are geographically contiguous. The algorithm used, which is actually part of leading methods developed for analysing large microarray datasets, is that of the k-means. Its goal is to minimize the within-group variance. The spatial constraint is provided by a list of connections between localities deduced from a Delaunay network. After testing each optimal group for the presence of spatial arrangement in the genetic data, the inferred genetic structure was compared with partitions obtained from other methods published for defining homogeneous groups (i.e. the Monmonier and SAMOVA algorithms). Competing biogeographical scenarios inferred from the k-means procedure were then compared and discussed to shed more light on colonization routes taken by the species.

The loop design of Kerr and Churchill is a clever application of incomplete blocks of size 2 to two-channel microarray experiments. In this paper, we extend the loop design to include more replicates, biological and technical replication, multi-factor experiments, and blocking. Loop and extended loop designs are shown to be more efficient than the reference design for any given number of arrays. We also show that adding new treatments to a loop design requires the same number of additional arrays as adding treatments to a reference design, with a greater gain in power. Given the flexibility of extended loop designs and their power, we propose that these should be the designs of choice for most experiments using two-channel microarrays.

Selective sweeps of variation caused by fixation of major genes may have a dramatic impact on the genetic gain from background polygenic variation, particularly in the genome regions closely linked to the major gene. The response to selection can be restrained because of the reduced selection intensity and the reduced effective population size caused by the increase in frequency of the major gene. In the context of a selected population where fixation of a known major gene is desired, the question arises as to which is the optimal path of increase in frequency of the gene so that the selective sweep of variation resulting from its fixation is minimized. Using basic theoretical arguments we propose a frequency path that maximizes simultaneously the effective population size applicable to the selected background and the selection intensity on the polygenic variation by minimizing the average squared selection intensity on the major gene over generations up to a given fixation time. We also propose the use of mating between carriers and non-carriers of the major gene, in order to promote the effective recombination between the major gene and its linked polygenic background. Using a locus-based computer simulation assuming different degrees of linkage, we show that the path proposed is more effective than a similar path recently published, and that the combination of the selection and mating methods provides an efficient way to palliate the negative effects of a selective sweep.

An interval quantitative trait locus (QTL) mapping method for complex polygenic diseases (as binary traits) showing QTL by environment interactions (QEI) was developed for outbred populations on a within-family basis. The main objectives, within the above context, were to investigate selection of genetic models and to compare liability or generalized interval mapping (GIM) and linear regression interval mapping (RIM) methods. Two different genetic models were used: one with main QTL and QEI effects (QEI model) and the other with only a main QTL effect (QTL model). Over 30 types of binary disease data as well as six types of continuous data were simulated and analysed by RIM and GIM. Using table values for significance testing, results show that RIM had an increased false detection rate (FDR) for testing interactions which was attributable to scale effects on the binary scale. GIM did not suffer from a high FDR for testing interactions. The use of empirical thresholds, which effectively means higher thresholds for RIM for testing interactions, could repair this increased FDR for RIM, but such empirical thresholds would have to be derived for each case because the amount of FDR depends on the incidence on the binary scale. RIM still suffered from higher biases (15–100% over- or under-estimation of true values) and high standard errors in QTL variance and location estimates than GIM for QEI models. Hence GIM is recommended for disease QTL mapping with QEI. In the presence of QEI, the model including QEI has more power (20–80% increase) to detect the QTL when the average QTL effect is small (in a situation where the model with a main QTL only is not too powerful). Top-down model selection is proposed in which a full test for QEI is conducted first and then the model is subsequently simplified. Methods and results will be applicable to human, plant and animal QTL mapping experiments.

We have increased the density of genetic markers on the Arabidopsis lyrata chromosomes AL6 and AL7 corresponding to the A. thaliana chromosome IV, in order to determine chromosome rearrangements between these two species, and to compare recombination fractions across the same intervals. We confirm the two rearrangements previously inferred (a reciprocal translocation and a large inversion, which we infer to be pericentric). By including markers around the centromere regions of A. thaliana chromosomes IV and V, we localize the AL6 centromere, and can localize the breakpoints of these chromosome rearrangements more precisely than previously. One translocation breakpoint was close to the centromere, and the other coincided with one end of the inversion, suggesting that a single event caused both rearrangements. At the resolution of our mapping, apart from these rearrangements, all other markers are in the same order in A. lyrata and A. thaliana. We could thus compare recombination rates in the two species. We found slightly higher values in A. thaliana, and a minimum estimate for regions not close to a centromere in A. lyrata is 4–5 centimorgans per megabase. The mapped region of AL7 includes the self-incompatibility loci (S-loci), and this region has been predicted to have lower recombination than elsewhere in the genome. We mapped 17 markers in a region of 1·23 Mb surrounding these loci, and compared the approximately 600 kb closest to the S-loci with the surrounding region of approximately the same size. There were significantly fewer recombination events in the closer than the more distant region, supporting the above prediction, but showing that the low recombination region is very limited in size.

Agrostis species are mainly used in athletic fields and golf courses. Their integrity is maintained by fungicides, which makes the development of disease-resistance varieties a high priority. However, there is a lack of knowledge about resistance (R) genes and their use for genetic improvement in Agrostis species. The objective of this study was to identify and clone constitutively expressed cDNAs encoding R gene-like (RGL) sequences from three Agrostis species (colonial bentgrass (A. capillaris L.), creeping bentgrass (A. stolonifera L.) and velvet bentgrass (A. canina L.)) by PCR-based motif-directed RNA fingerprinting towards relatively conserved nucleotide binding site (NBS) domains. Sixty-one constitutively expressed cDNA sequences were identified and characterized. Sequence analysis of ESTs and probable translation products revealed that RGLs are highly conserved among these three Agrostis species. Fifteen of them were shown to share conserved motifs found in other plant disease resistance genes such as MLA13, Xa1, YR6, YR23 and RPP5. The molecular evolutionary forces, analysed using the Ka/Ks ratio, reflected purifying selection both on NBS and leucine-rich repeat (LRR) intervening regions of discovered RGL sequences in these species. This study presents, for the first time, isolation and characterization of constitutively expressed RGL sequences from Agrostis species revealing the presence of TNL (TIR-NBS-LRR) type R genes in monocot plants. The characterized RGLs will further enhance knowledge on the molecular evolution of the R gene family in grasses.

Unbiased or upper limit estimates of the rate (U) of genomic mutations to mildly deleterious alleles are crucial in genetic and conservation studies and in human health care. However, only a few estimates of the lower bounds of U are available. We present a fairly robust estimation that yields an upper limit of U and a nearly unbiased estimate of the per generation fitness decline due to new deleterious mutations. We applied the approach to three species of the freshwater microcrustacean Daphnia and revealed that the upper limit of U for egg survivorship is 0·73 (SD=0·30) in 14 D. pulicaria populations. For the first four clutches, per generation decline in fecundity due to deleterious mutations ranged from 2·2% to 7·8% in 20 D. pulex populations and from 1·1% to 5·1% in 8 D. obtusa populations. These results indicate the mutation pressure is high in natural Daphnia populations. The approach investigated here provides a potential way to quickly and conveniently characterize U and per generation effects of deleterious genomic mutations on fitness or its important components such as fecundity.

Indirect estimates of the genomic rate of deleterious mutations (λ), their average homozygous effect (s) and their degree of dominance (h) can be obtained from genetic parameters of natural populations, assuming that the frequencies of the loci controlling a given fitness trait are at mutation–selection equilibrium. In 1996, H.-W. Deng and M. Lynch developed a general methodology for obtaining these estimates from inbreeding/outbreeding experiments. The prediction of the sign and magnitude of the biases incurred by these estimators is essential for a correct interpretation of the empirical results. However, the assessment of these biases has been tested so far under a rather limited model of the distribution of dominance effects. In this paper, the application of this method to outbred populations is evaluated, focusing on the level of variation in h values (σh2) and the magnitude of the negative correlation (rs,h) between s and h. It is shown that the method produces upwardly biased estimates of λ and downwardly biased estimates of the average s in the reference situation where rs,h=0, particularly for large values of σh2, and biases of different sign depending on the magnitude of the correlation. A modification of the method, substituting the estimates of the average h for alternative ones, allows estimates to be obtained with little or no bias for the case of rs,h=0, but is otherwise biased. Information on rs,h and σh2, gathered from mutation-accumulation experiments, suggests that σh2 may be rather large and rs,h is usually negative but not higher than about −0·2, although the data are scarce and noisy, and should be used with caution.

Prior to cranial neural tube closure, the neural folds adopt a biconvex morphology which is thought to be due to expansion of the underlying mesenchyme. Dorso-lateral hinge points (DLHPs) then form, which allow the dorsal tips of the neural folds to ‘flip around’ resulting in apposition of the tips and facilitating subsequent fusion. Cranial closure is particularly prone to perturbation, leading to exencephaly in many mouse mutants and as a result of a variety of teratogenic influences. This may reflect mechanical tensions affecting the closing cranial neural folds. For example, the presence of ventral flexures of the body axis at the mid- and forebrain levels mechanically opposes the formation of DLHPs. Several processes have been implicated as important in overcoming these mechanical tensions, thereby assisting in cranial neural tube closure. These include contraction of actin microfilaments at the luminal surface of the neuroepithelium and apoptosis in the dorsal and dorsolateral neuroepithelium. The latter may act to increase flexibility in the dorsal neural folds, enhancing DLHP formation. Neural crest cells (NCC) originate in the dorsal tips of the neuroepithelium and undergo an epithelial-to-mesenchymal transition, allowing them to delaminate, exit the neuroepithelium and migrate extensively throughout the embryo to form numerous derivatives. We hypothesized that delamination of the NCC from the neuroepithelium may enhance the mechanical flexibility of the dorsal tips of the neural folds, allowing the ‘flip around’ event to occur.