This paper examines the properties of likelihood maps generated by interval mapping (IM) and composite interval mapping (CIM), two widely used methods for detecting quantitative trait loci (QTLs). We evaluate the usefulness of interpretations of entire maps, rather than only evaluating summary statistics that consider isolated features of maps. A simulation study was performed in which traits with varying genetic architectures, including 20–40 QTLs per chromosome, were examined with both IM and CIM under different marker densities and sample sizes. IM was found to be an unreliable tool for precise estimation of the number and locations of individual QTLs, although it has greater power for simply detecting the presence of QTLs than CIM. The ability of CIM to resolve the correct number of QTLs and to estimate their locations correctly is good if there are three or fewer QTLs per 100 centiMorgans, but can lead to erroneous inferences for more complex architectures. When the underlying genetic architecture of a trait consists of several QTLs with randomly distributed effects and locations likelihood profiles were often indicative of a few underlying genes of large effect. Studies that have detected more than a few QTLs per chromosome should be interpreted with caution.

We have investigated the pattern of DNA sequence variation at the exuperantia2 locus in Drosophila pseudoobscura. This adds to the increasing dataset of genetic variation in D. pseudoobscura, a useful model species for evolutionary genetic studies. The level of silent site nucleotide diversity and the divergence from an outgroup Drosophila miranda are comparable with those for other X-linked loci. One peculiar pattern at the exu2 locus of D. pseudoobscura is a complete linkage disequilibrium between two SNPs, one of which is a replacement site. As a result, there are two distinct haplotype groups in our dataset. Based upon the comparisons with the outgroup sequences from D. miranda and Drosophila persimilis, we show that the newly derived haplotype group has lower diversity than the ancestral haplotype group. The pattern of protein evolution at exu2 shows some deviation from the neutral model. Together, these and other characteristics of the exu2 locus suggest the action of selection on the pattern of SNP variation, consistent with a partial selective sweep associated with the newly derived haplotype.

Many diseases show dichotomous phenotypic variation but do not follow a simple Mendelian pattern of inheritance. Variances of these binary diseases are presumably controlled by multiple loci and environmental variants. A least-squares method has been developed for mapping such complex disease loci by treating the binary phenotypes (0 and 1) as if they were continuous. However, the least-squares method is not recommended because of its ad hoc nature. Maximum Likelihood (ML) and Bayesian methods have also been developed for binary disease mapping by incorporating the discrete nature of the phenotypic distribution. In the ML analysis, the likelihood function is usually maximized using some complicated maximization algorithms (e.g. the Newton–Raphson or the simplex algorithm). Under the threshold model of binary disease, we develop an Expectation Maximization (EM) algorithm to solve for the maximum likelihood estimates (MLEs). The new EM algorithm is developed by treating both the unobserved genotype and the disease liability as missing values. As a result, the EM iteration equations have the same form as the normal equation system in linear regression. The EM algorithm is further modified to take into account sexual dimorphism in the linkage maps. Applying the EM-implemented ML method to a four-way-cross mouse family, we detected two regions on the fourth chromosome that have evidence of QTLs controlling the segregation of fibrosarcoma, a form of connective tissue cancer. The two QTLs explain 50–60% of the variance in the disease liability. We also applied a Bayesian method previously developed (modified to take into account sex-specific maps) to this data set and detected one additional QTL on chromosome 13 that explains another 26% of the variance of the disease liability. All the QTLs detected primarily show dominance effects.

During courtship, visual and chemical signals are often exchanged between the sexes. The proper exchange of such signals ensures intraspecific recognition. We have examined the genetic basis of interspecific differences in male mating behaviour and pheromone concentration between Drosophila simulans and D. sechellia by using Drosophila simulans/D. sechellia introgression lines. Our results show a majority of quantitative trait loci (QTLs) explaining variation in both male mating behaviour and pheromone concentration to be located on the third chromosome. One QTL found on the third chromosome explains variation in time needed to start courtship and copulation as well as time spent courting. The position of such QTL (approximately 84A–88B) with effects on courtship and copulation aspects of mating includes the candidate sex determination gene doublesex (84E5–6) and Voila (86E1–2), a gene that affects male courtship in D. melanogaster. One additional third chromosome QTL explained variation in 7-tricosene pheromone concentrations among males. The interval mapping position of this QTL (approximately 68E–76E) did not overlap with the position detected for differences in mating behaviour and the intervals did not include candidate genes previously identified as having an effect on D. melanogaster cuticular hydrocarbon production. We did not detect any directionality of the effect of Drosophila sechellia allele introgressions in male mating recognition.

New paradigms in genetics have increased the chance of finding genes that appear redundant but in fact may have been preserved due to a small level of positive selection potential acting during each generation. Monitoring changes in genotypic frequencies within and between generations allows the dissection of the fertility, viability and meiotic drive selection components acting on such genes in natural and experimental populations. Here, a formal maximum likelihood procedure is developed to identify and estimate these selection components in highly selfing populations by fitting the time-dependent solutions for genotypic frequencies to observed multigenerational counts. With adult census alone, we can not simultaneously estimate all three selection components considered. In such cases, we instead consider a hierarchy of 11 models with either fewer selection components, complete dominance, or multiplicative meiotic drive with a single parameter. We identify the best-fitting of these models by applying likelihood ratio tests to nested models and Akaike's Information Criterion (AIC) and the Bayesian Information Criterion (BIC) to non-nested models. With seed census, fertility and viability selection are not distinguishable and thus can only be estimated jointly. A combination of joint seed and adult census data allows us to estimate all three selection components simultaneously. Simulated data validate the estimation procedure and provide some practical guidelines for experimental design. An application to Arabidopsis data establishes that viability selection is the major selective force acting on the ACT2 actin gene in laboratory-grown Arabidopsis populations.

We analysed the distribution of transposable elements (TEs) in 100 aligned pairs of orthologous intergenic regions from the mouse and human genomes. Within these regions, conserved segments of high similarity between the two species alternate with segments of low similarity. Identifiable TEs comprise 40–60% of segments of low similarity. Within such segments, a particular copy of a TE found in one species has no orthologue in the other. Overall, TEs comprise only approximately 20% of conserved segments. However, TEs from two families, MIR and L2, are rather common within conserved segments. Statistical analysis of the distributions of TEs suggests that a majority of the MIR and L2 elements present in murine intergenic regions have human orthologues. These elements must have been present in the common ancestor of human and mouse and have remained under substantial negative selection that prevented their divergence beyond recognition. If so, recruitment of MIR- and L2-derived sequences to perform a function that increases host fitness is rather common, with at least two such events per host gene. The central part of the MIR consensus sequence is over-represented in conserved segments given its background frequency in the genome, suggesting that it is under the strongest selective constraint.

Simultaneous analysis of correlated traits that change with time is an important issue in genetic analyses. Several methodologies have already been proposed for the genetic analysis of longitudinal data on single traits, in particular random regression and character process models. Although the latter proved, in most cases, to compare favourably to alternative approaches for analysis of single function-valued traits, they do not allow a straightforward extension to the multivariate case. In this paper, another methodology (structured antedependence models) is proposed, and methods are derived for the genetic analysis of two or more correlated function-valued traits. Multivariate analyses are presented of fertility and mortality in Drosophila and of milk, fat and protein yields in dairy cattle. These models offer a substantial flexibility for the correlation structure, even in the case of complex non-stationary patterns, and perform better than multivariate random regression models, with fewer parameters.

The medaka fish transposable element, Tol2, is a member of the hAT family of transposons. It has been directly demonstrated to be active and two mRNAs, differing in length, have been isolated. They cover exons 1–4 and exons 2–4 and the longer form has already been proven to catalyse transposition reactions. However, the function of the shorter mRNA in medaka cells has hitherto remained unclear. In the present study, first we constructed a quantitative system to detect Tol2 excision using an indicator plasmid carrying a non-autonomous Tol2 within its lacZ gene; second we injected mRNAs with the plasmid into medaka eggs. Excision of Tol2 was detected as E. coli blue colonies caused by the recovery of lacZ activity. Addition of the longer mRNA increased excision, but the shorter did not. Moreover, co-injection of both mRNAs greatly lowered the frequency compared with the case of treatment with the longer mRNA alone. These results indicate that the shorter mRNA has an inhibitory effect on the excision reaction, and that the N-terminal region of the transposase encoded by exon 1, including a BED zinc finger, presumably plays an important role in excision. Here, we suggest a regulatory mechanism of Tol2 transposition involving the expression of these mRNAs.

Rare, random mutations were induced in budding yeast by ethyl methanesulfonate (EMS). Clones known to bear a single non-neutral mutation were used to obtain mutant heterozygotes and mutant homozygotes that were later compared with wild-type homozygotes. The average homozygous effect of mutation was an approximately 2% decrease in the growth rate. In heterozygotes, the harmful effect of these relatively mild mutations was reduced approximately fivefold. In a test of epistasis, two heterozygous mutant loci were paired at random. Fitness of the double mutants was best explained by multiplicative action of effects at single loci, with little evidence for epistasis and essentially excluding synergism. In other experiments, the same mutations in haploid and heterozygous diploid clones were compared. Regardless of the haploid phenotypes, mildly deleterious or lethal, fitness of the heterozygotes was decreased by less than half a per cent on average. In general, the results presented here suggest that most mutations tend to exhibit small and weakly interacting effects in heterozygous loci regardless of how harmful they are in haploids or homozygotes.

Maternal effects play an important role in fitness and other aspects of individual performance in many species, particularly mammalian, yet their impact on genetic variation within species and its rate of loss during selection has been neglected. In this paper we extend the theory of expected long-term genetic contributions to include maternal effects, and tested the accuracy of predicted rates of inbreeding for populations under mass selection by comparison with simulations. The model includes selective advantages of direct and maternal additive genetic effects, and also the selective advantage of a common maternal environmental effect. The population structures investigated had a fixed number of dams per sire and fixed family size. Most prediction errors of the rate of inbreeding (ΔF) were less than 8% of the simulated means and were lower in magnitude than the prediction errors of genetic gain (ΔG). The predictions of ΔG from contributions equalled previously published predictions. A variation in maternal genetic effects resulted in a much larger ΔF than for an equally sized variation in common maternal environmental effects. For a fixed genetic gain, ΔF increased as the maternal heritability increased. The influence of family size, mating ratio and age structure on ΔF was greater with maternal effects than with only direct genetic effects included. In conclusion, maternal effects can be a very important aspect to consider when predicting ΔF in populations under selection, and the developed methodology gives good predictions.

222 cultivated (Vitis vinifera) and 22 wild (V. vinifera ssp. sylvestris) grape accessions were analysed for genetic diversity and differentiation at eight microsatellite loci. A total of 94 alleles were detected, with extensive polymorphism among the accessions. Multivariate relationships among accessions revealed 16 genetic groups structured into three clusters, supporting the classical eco-geographic grouping of grape cultivars: occidentalis, pontica and orientalis. French cultivars appeared to be distinct and showed close affinity to the wild progenitor, ssp. sylvestris from south-western France (Pyrenees) and Tunisia, probably reflecting the origin and domestication history of many of the old wine cultivars from France. There was appreciable level of differentiation between table and wine grape cultivars, and the Muscat types were somewhat distinct within the wine grapes. Contingency χ2 analysis indicated significant heterogeneity in allele frequencies among groups at all loci. The observed heterozygosities for different groups ranged from 0·625 to 0·9 with an overall average of 0·771. Genetic relationships among groups suggested hierarchical differentiation within cultivated grape. The gene diversity analysis indicated narrow divergence among groups and that most variation was found within groups (∼85%). Partitioning of diversity suggested that the remaining variation is somewhat structured hierarchically at different levels of differentiation. The overall organization of genetic diversity suggests that the germplasm of cultivated grape represents a single complex gene pool and that its structure is determined by strong artificial selection and a vegetative mode of reproduction.

No exact method for determining genotypic and identity-by-descent probabilities is available for large complex pedigrees. Approximate methods for such pedigrees cannot be guaranteed to be unbiased. A new method is proposed that uses the Metropolis–Hastings algorithm to sample a Markov chain of descent graphs which fit the pedigree and known genotypes. Unknown genotypes are determined from each descent graph. Genotypic probabilities are estimated as their means. The algorithm is shown to be unbiased for small complex pedigrees and feasible and consistent for moderately large complex pedigrees.

Intrapopulation variability in the seasonal regulation of insect lifecycles has been shown to be due partly to genetic changes. Selection for insecticide resistance in the codling moth Cydia pomonella results from allelic substitution at two to three loci in south-eastern French populations of this species. However, such an adaptive process has been associated with an increased heterogeneity in the developmental responses to climatic factors such as temperature. In this paper, we investigate whether such pleiotropic effects of resistance on development induce a significant discrepancy in seasonal regulation in this species. The seasonal changes in a susceptible and two insecticide-resistant homozygous genotypes of C. pomonella, as well as their reciprocal F1 progeny, were followed under natural conditions during the reproductive season through the emergence events of adults, within-generation developmental rates and the number of generations. A significant delay in the occurrences of homozygous resistant genotypes resulted from significantly lower pre-imaginal developmental times relative to homozygous susceptible ones. Subsequent assessment of the number of generations indicated significantly higher diapause propensities in carriers of the resistance alleles (37·0–76·2%) than in susceptible homozygotes (6·7%), which mostly pupated towards a third generation of adults. In the light of these findings, pleiotropic effects of adaptive changes might be a crucial source of divergence in seasonal regulation at the population level, involving significant life-history trade-offs. In addition to man-made selective factors during the reproductive season, such an effect on the lifecycle could be a key component in the process of selection for resistance genes in south-eastern France C. pomonella populations.

Several forces may affect the distribution of genetic diversity in natural populations when compared to what is expected in a random-mating, constant size population of neutral genes. One solution for unravelling their respective influence is to study several genes at once in order to better reflect the true genealogy. Here we reconstruct the evolutionary history of the freshwater snail Biomphalaria pfeifferi over its entire distribution, using eight African populations, and three congeneric species as outgroups. A phylogenetic analysis was conducted using amplified fragment length polymorphism markers, and sequences at eight nuclear non-coding loci and one mitochondrial gene were used to analyse population structure. The geographic distribution of variation suggests greater affinities within than among regions. The pattern of variability at both the nuclear and mitochondrial DNA (mtDNA) loci is consistent with a bottleneck, although population structure may also partly explain our results. Our results are also indicative of the role of selection, whether positive or purifying, in the mtDNA. This highlights the fact that the interfering influences of population structure, demography and selection on molecular variation are not easily distinguished.