The co-occurrence of three chromosome-wide phenomena – imprinting, facultative heterochromatization and diffuse centromere – in the mealybug Planococcus lilacinus makes investigation of the genomics of this species an attractive prospect. In order to estimate the complexity of the genome of this species, 300 random stretches of its DNA, constituting ∼0·1% of the genome, were sequenced. Coding sequences appear to constitute ∼53·5%, repeat sequences ∼44·5% and non-coding single-copy sequences ∼2% of the genome. The proportion of repetitive sequences in the mealybug is higher than that in the fruit fly Drosophila melanogaster (∼30%). The mealybug genome (∼220 Mb) is about 1·3 times the size of the fly genome (∼165 Mb) and its GC content (∼35%) less than that of the fly genome (∼40%). The relative abundance of various dinucleotides, as analysed by the method of Gentles and Karlin, shows that the dinucleotide signatures of the two species are moderately similar and that in the mealybug there is neither over-representation nor under-representation of any dinucleotide.

Most aphids are cyclic parthenogens, so are ideal models in studies of the mechanisms and consequences of sex and recombination. However, owing to a shortage of physical and genetic markers, there have been few studies of the most fundamental genetic processes in these organisms. For example, it is not known whether autosomal segregation during male spermatogenesis is in Mendelian proportions: we address that question here. The aphid Myzus persicae has a typical karyotype of 2n = 12 in females (XX), while males are XO (2n = 11). During male meiosis, only the spermatocytes with an X chromosome are viable. We hypothesized that assortment of autosomes might be non-random because chromosomal imprinting leading to elimination of the paternal autosomes is seen in the closely related coccoids. In other aphid models, we have observed segregation distortions at single microsatellite loci (Wilson, 2000). Such distortions may have nothing to do with ‘selfish’ behaviour, but may be caused by mutation accumulation causing fitness differentials. Thus single-locus distortions might be predicted to be more likely to be detected via the male lines of clones that have lost the ability to reproduce sexually (male-producing obligate parthenogenesis (androcyclic)). Using microsatellites we show that genetic imprinting or selfish autosome behaviour does not occur in male M. persicae. Generally, loci segregated in Mendelian proportions in both sexes of cyclically parthenogenetic (holocyclic) clones. However, in androcyclic clones, segregation distortions consistently involved the same two autosomes. This is consistent with linkage of markers to deleterious mutations associated with a loss of sexual reproduction.

Natural populations are structured spatially into local populations and genetically into diverse ‘genetic backgrounds’ defined by different combinations of selected alleles. If selection maintains genetic backgrounds at constant frequency then neutral diversity is enhanced. By contrast, if background frequencies fluctuate then diversity is reduced. Provided that the population size of each background is large enough, these effects can be described by the structured coalescent process. Almost all the extant results based on the coalescent deal with a single selected locus. Yet we know that very large numbers of genes are under selection and that any substantial effects are likely to be due to the cumulative effects of many loci. Here, we set up a general framework for the extension of the coalescent to multilocus scenarios and we use it to study the simplest model, where strong balancing selection acting on a set of n loci maintains 2n backgrounds at constant frequencies and at linkage equilibrium. Analytical results show that the expected linked neutral diversity increases exponentially with the number of selected loci and can become extremely large. However, simulation results reveal that the structured coalescent approach breaks down when the number of backgrounds approaches the population size, because of stochastic fluctuations in background frequencies. A new method is needed to extend the structured coalescent to cases with large numbers of backgrounds.

Clines for size and stress resistance traits have been described for several Drosophila species and replicable clines across different species may indicate climatic selection. Here we consider clines in stress resistance traits in an Australian endemic species, D. serrata, by comparing levels of variation within and among isofemale lines initiated with flies collected from the eastern coast of Australia. We also consider clinal variation in chill coma recovery, a trait that has recently been shown to exhibit high levels of variation among Drosophila species. Patterns were compared with those in the cosmopolitan species D. melanogaster from the same area. Both desiccation and starvation resistance showed no clinal pattern despite heritable variation among isofemale lines. In contrast chill coma resistance exhibited a linear cline in the anticipated direction, resistance increasing with latitude. Body size was measured as wing length and body weight. Both traits showed geographic variation and strong non-linear clines with a sharp reduction in size in the tropics. These results are discussed in the context of climatic selection and evolutionary processes limiting species borders.

Reproduction in individual animals of sexual species depends largely upon their ability to detect and distinguish specific signal(s) among those produced by various potential sexual partners. In Drosophila melanogaster males, there is a natural polymorphism for discrimination of female and male principal pheromones that segregates with chromosome 3. We have mapped two loci on chromosome 3 that change sex-pheromone discrimination in males. We successively exploited meiotic recombination, deficiencies and enhancer-trap strains; excision of the transposon in two selected enhancer-trap strains clearly reverted the discrimination phenotype. These results indicate that pheromonal discrimination is a character that can be genetically manipulated, and provide further insights into the evolution of the specific mate recognition system.

We investigate how sampling of parents or children based on their extreme phenotypic values selected from clinical databases would affect the power of identification of quantitative trait loci (QTL) by a transmission disequilibrium test (TDT). We consider three selective sampling schemes based on the selection of phenotypic values of parents or children in nuclear families: (1) two children, one of extreme value, the other random; (2) two children extremely discordant; (3) one parent of extreme value. Other family members not specified will be recruited randomly with regard to phenotypic values. Our study shows that the second sampling scheme can always enhance the power for QTL identification, sometimes dramatically so. The increase in the statistical power of the TDT is particularly dramatic when h2 at the QTL under test is small or intermediate (e.g. 0·05 or 0·10). For the other two sampling schemes, under dominant effects at the QTL, the power is always increased relative to random sampling; however, under recessive or additive genetic effects, the power gain is generally minor or even decreased a little sometimes. Allele frequencies at the QTL and the selection stringency are important for determining the effect of selective sampling on the power of QTL identification. Our study is useful as a practical guideline on how to perform the TDT efficiently in practice by taking advantage of the extensive databases accumulated that are enriched with people of extreme phenotypic values.

Here, we describe a randomization testing strategy for mapping interacting quantitative trait loci (QTLs). In a forward selection strategy, non-interacting QTLs and simultaneously mapped interacting QTL pairs are added to a total genetic model. Simultaneous mapping of epistatic QTLs increases the power of the mapping strategy by allowing detection of interacting QTL pairs where none of the QTL can be detected by their marginal additive and dominance effects. Randomization testing is used to derive empirical significance thresholds for every model selection step in the procedure. A simulation study was used to evaluate the statistical properties of the proposed randomization tests and for which types of epistasis simultaneous mapping of epistatic QTLs adds power. Least squares regression was used for QTL parameter estimation but any other QTL mapping method can be used. A genetic algorithm was used to search for interacting QTL pairs, which makes the proposed strategy feasible for single processor computers. We believe that this method will facilitate the evaluation of the importance at epistatic interaction among QTLs controlling multifactorial traits and disorders.

Epistatic variance can be an important source of variation for complex traits. However, detecting epistatic effects is difficult primarily due to insufficient sample sizes and lack of robust statistical methods. In this paper, we develop a Bayesian method to map multiple quantitative trait loci (QTLs) with epistatic effects. The method can map QTLs in complicated mating designs derived from the cross of two inbred lines. In addition to mapping QTLs for quantitative traits, the proposed method can even map genes underlying binary traits such as disease susceptibility using the threshold model. The parameters of interest are various QTL effects, including additive, dominance and epistatic effects of QTLs, the locations of identified QTLs and even the number of QTLs. When the number of QTLs is treated as an unknown parameter, the dimension of the model becomes a variable. This requires the reversible jump Markov chain Monte Carlo algorithm. The utility of the proposed method is demonstrated through analysis of simulation data.