High-resolution linkage map for two honeybee chromosomes: the hotspot quest.

Meiotic recombination is a fundamental process ensuring proper disjunction of homologous chromosomes and allele shuffling in successive generations. In many species, this cellular mechanism occurs heterogeneously along chromosomes and mostly concentrates in tiny fragments called recombination hotspots. Specific DNA motifs have been shown to initiate recombination in these hotspots in mammals, fission yeast and drosophila. The aim of this study was to check whether recombination also occurs in a heterogeneous fashion in the highly recombinogenic honeybee genome and whether this heterogeneity can be connected with specific DNA motifs. We completed a previous picture drawn from a routine genetic map built with an average resolution of 93 kb. We focused on the two smallest honeybee chromosomes to increase the resolution and even zoomed at very high resolution (3.6 kb) on a fragment of 300 kb. Recombination rates measured in these fragments were placed in relation with occurrence of 30 previously described motifs through a Poisson regression model. A selection procedure suitable for correlated variables was applied to keep significant motifs. These fine and ultra-fine mappings show that recombination rate is significantly heterogeneous although poorly contrasted between high and low recombination rate, contrarily to most model species. We show that recombination rate is probably associated with the DNA methylation state. Moreover, three motifs (CGCA, GCCGC and CCAAT) are good candidates of signals promoting recombination. Their influence is however moderate, doubling at most the recombination rate. This discovery extends the way to recombination dissection in insects.

Fine scale analysis of crossover and non-crossover and detection of recombination sequence motifs in the honeybee (Apis mellifera).

BACKGROUND: Meiotic exchanges are non-uniformly distributed across the genome of most studied organisms. This uneven distribution suggests that recombination is initiated by specific signals and/or regulations. Some of these signals were recently identified in humans and mice. However, it is unclear whether or not sequence signals are also involved in chromosomal recombination of insects. METHODOLOGY: We analyzed recombination frequencies in the honeybee, in which genome sequencing provided a large amount of SNPs spread over the entire set of chromosomes. As the genome sequences were obtained from a pool of haploid males, which were the progeny of a single queen, an oocyte method (study of recombination on haploid males that develop from unfertilized eggs and hence are the direct reflect of female gametes haplotypes) was developed to detect recombined pairs of SNP sites. Sequences were further compared between recombinant and non-recombinant fragments to detect recombination-specific motifs. CONCLUSIONS: Recombination events between adjacent SNP sites were detected at an average distance of 92 bp and revealed the existence of high rates of recombination events. This study also shows the presence of conversion without crossover (i. e. non-crossover) events, the number of which largely outnumbers that of crossover events. Furthermore the comparison of sequences that have undergone recombination with sequences that have not, led to the discovery of sequence motifs (CGCA, GCCGC, CCGCA), which may correspond to recombination signals.

Novel microsatellite DNA loci for Bombus terrestris (Linnaeus, 1758).

We present details and characteristics of 123 novel polymorphic microsatellite DNA loci for Bombus terrestris. Thirty-four of these loci have been tested in nine other Bombus species and 25 of them showed polymorphisms in at least one species. These microsatellite DNA loci together with the already established 60 loci will be useful for characterizing wild and managed populations of B. terrestris and other Bombus species as well as for detailed genetic studies in including mapping studies and genome annotations.

Spatiotemporal genetic differentiation of Cuban natural populations of the pink shrimp Farfantepenaeus notialis.

We analyzed the spatiotemporal genetic structure of Farfantepenaeus notialis populations using five microsatellites loci in order to understand the influence of natural events such as hurricanes on the genetic drift/migration balance as the main cause for the variation of allele frequencies over time. The results were compared with the previous ones obtained from allozymes and mtDNA. High and stable genetic diversity levels (He=0.879+/-0.0015) were found over eight years for the populations that inhabit the south Cuban platform, however significant changes of allele frequencies were detected over time. The F(ST) estimates, albeit low, revealed significant differences among populations inside the Ana Maria Gulf for 1995 but not for the 1999 and 2003 samples. The F(ST), AMOVA and the genetic distance analysis revealed the instability of the genetic structure over time in accordance with allozymes results. The correspondence of the microsatellite results with those obtained from allozymes confirm the effects of migration enhanced by natural events as the main cause of the temporal variation of allele frequencies. The genetic drift effect was discarded through the evaluation of Ne and the M ratio, while natural selection effects were rejected because of the lowest probability of microsatellite loci being under selective pressures. The microsatellite data are also consistent with the results obtained with mtDNA in detecting significant and persistent genetic differences between the Gulfs of Ana María and Batabanó for the years 1995 and 2003.

Manual superscaffolding of honey bee (Apis mellifera) chromosomes 12-16: implications for the draft genome assembly version 4, gene annotation, and chromosome structure.

The euchromatic arms of the five smallest telocentric chromosomes in the honey bee genome draft Assembly v4 were manually connected into superscaffolds. This effort reduced chromosomes 12-16 from 30, 21, 25, 42, and 21 mapped scaffolds to five, four, five, six, and five superscaffolds, respectively, and incorporated 178 unmapped contigs and scaffolds totalling 2.6 Mb, a 6.4% increase in length. The superscaffolds extend from the genetically mapped location of the centromere to their identified distal telomeres on the long arms. Only two major misassemblies of 146 kb and 65 kb sections were identified in this 23% of the mapped assembly. Nine duplicate gene models on chromosomes 15 and 16 were made redundant, while another 15 gene models were improved, most spectacularly the MAD (MAX dimerization protein) gene which extends across 11 scaffolds for at least 400 kb.

A third-generation microsatellite-based linkage map of the honey bee, Apis mellifera, and its comparison with the sequence-based physical map.

BACKGROUND: The honey bee is a key model for social behavior and this feature led to the selection of the species for genome sequencing. A genetic map is a necessary companion to the sequence. In addition, because there was originally no physical map for the honey bee genome project, a meiotic map was the only resource for organizing the sequence assembly on the chromosomes. RESULTS: We present the genetic (meiotic) map here and describe the main features that emerged from comparison with the sequence-based physical map. The genetic map of the honey bee is saturated and the chromosomes are oriented from the centromeric to the telomeric regions. The map is based on 2,008 markers and is about 40 Morgans (M) long, resulting in a marker density of one every 2.05 centiMorgans (cM). For the 186 megabases (Mb) of the genome mapped and assembled, this corresponds to a very high average recombination rate of 22.04 cM/Mb. Honey bee meiosis shows a relatively homogeneous recombination rate along and across chromosomes, as well as within and between individuals. Interference is higher than inferred from the Kosambi function of distance. In addition, numerous recombination hotspots are dispersed over the genome. CONCLUSION: The very large genetic length of the honey bee genome, its small physical size and an almost complete genome sequence with a relatively low number of genes suggest a very promising future for association mapping in the honey bee, particularly as the existence of haploid males allows easy bulk segregant analysis.

Control of reproductive dominance by the thelytoky gene in honeybees.

Differentiation into castes and reproductive division of labour are a characteristics of eusocial insects. Caste determination occurs at an early stage of larval development in social bees and is achieved via differential nutrition irrespective of the genotype. Workers are usually subordinate to the queen and altruistically refrain from reproduction. Workers of the Cape honeybee (Apis mellifera capensis) do not necessarily refrain from reproduction. They have the unique ability to produce female offspring parthenogenetically (thelytoky) and can develop into 'pseudoqueens'. Although these are morphologically workers, they develop a queen-like phenotype with respect to physiology and behaviour. Thelytoky is determined by a single gene (th) and we show that this gene also influences other traits related to the queen phenotype, including egg production and queen pheromone synthesis. Using 566 microsatellite markers, we mapped this gene to chromosome 13 and identified a candidate locus thelytoky, similar to grainy head (a transcription factor), which has been shown to be highly expressed in queens of eusocial insects. We therefore suggest that this gene is not only important for determining the pseudoqueen phenotype in A. m. capensis workers, but is also of general importance in regulating the gene cascades controlling reproduction and sterility in female social bees.

The genome of Apis mellifera: dialog between linkage mapping and sequence assembly.

Two independent genome projects for the honey bee, a microsatellite linkage map and a genome sequence assembly, interactively produced an almost complete organization of the euchromatic genome. Assembly 4.0 now includes 626 scaffolds that were ordered and oriented into chromosomes according to the framework provided by the third-generation linkage map (AmelMap3). Each construct was used to control the quality of the other. The co-linearity of markers in the sequence and the map is almost perfect and argues in favor of the high quality of both.

Exceptionally high levels of recombination across the honey bee genome.

The first draft of the honey bee genome sequence and improved genetic maps are utilized to analyze a genome displaying 10 times higher levels of recombination (19 cM/Mb) than previously analyzed genomes of higher eukaryotes. The exceptionally high recombination rate is distributed genome-wide, but varies by two orders of magnitude. Analysis of chromosome, sequence, and gene parameters with respect to recombination showed that local recombination rate is associated with distance to the telomere, GC content, and the number of simple repeats as described for low-recombining genomes. Recombination rate does not decrease with chromosome size. On average 5.7 recombination events per chromosome pair per meiosis are found in the honey bee genome. This contrasts with a wide range of taxa that have a uniform recombination frequency of about 1.6 per chromosome pair. The excess of recombination activity does not support a mechanistic role of recombination in stabilizing pairs of homologous chromosome during chromosome pairing. Recombination rate is associated with gene size, suggesting that introns are larger in regions of low recombination and may improve the efficacy of selection in these regions. Very few transposons and no retrotransposons are present in the high-recombining genome. We propose evolutionary explanations for the exceptionally high genome-wide recombination rate.

The invasive Korea and Japan types of Varroa destructor, ectoparasitic mites of the Western honeybee (Apis mellifera), are two partly isolated clones.

Varroa destructor, now a major pest of the Western honeybee, Apis mellifera, switched from its original host, the Eastern honeybee, A. cerana, ca. 50 years ago. So far, only two out of several known mitochondrial haplotypes of V. destructor have been found to be capable of reproducing on A. mellifera (Korea and Japan). These haplotypes are associated in almost complete cytonuclear disequilibrium to diagnostic alleles at 11 microsatellite loci. By contrast, microsatellite polymorphism within each type is virtually absent, because of a severe bottleneck at the time of host change. Accordingly, 12 mitochondrial sequences of 5185 nucleotides displayed 0.40% of nucleotide divergence between haplotypes and no intra haplotype variation. Hence, each type has a quasi-clonal structure. The nascent intratype variability is subsequent to the clone formation 50 years ago: in both types the variant alleles differ from the most common by one (in 10 cases), two (five cases) or three (one case) repeated motifs. In addition to individuals of the two 'pure' types, five F1 hybrids and 19 recombinant individuals (Japan alleles introgressed into the Korea genetic background) were detected. The existence of F1 and recombinant individuals in admixed populations requires that double infestations of honeybee cells occur in a high proportion but the persistence of pure types suggests a post-zygotic isolation between the two clones.

Whole-genome scan in thelytokous-laying workers of the Cape honeybee (Apis mellifera capensis): central fusion, reduced recombination rates and centromere mapping using half-tetrad analysis.

While workers of almost all subspecies of honeybee are able to lay only haploid male eggs, Apis mellifera capensis workers are able to produce diploid female eggs by thelytokous parthenogenesis. Cytological analyses have shown that during parthenogenesis, egg diploidy is restored by fusion of the two central meiotic products. This peculiarity of the Cape bee preserves two products of a single meiosis in the daughters and can be used to map centromere positions using half-tetrad analysis. In this study, we use the thelytokous progenies of A. m. capensis workers and a sample of individuals from a naturally occurring A. m. capensis thelytokous clone to map centromere position for most of the linkage groups of the honeybee. We also show that the recombination rate is reduced by >10-fold during the meiosis of A. m. capensis workers. This reduction is restricted to thelytokous parthenogenesis of capensis workers and is not observed in the meiosis of queen within the same subspecies or in arrhenotokous workers of another subspecies. The reduced rate of recombination seems to be associated with negative crossover interference. These results are discussed in relation to evolution of thelytokous parthenogenesis and maintenance of heterozygosity and female sex after thelytoky.

A microsatellite-based linkage map of the honeybee, Apis mellifera L.

A linkage map for the honeybee (Apis mellifera) was constructed mainly from the progeny of two hybrid queens (A. m. ligustica x A. m. mellifera). A total of 541 loci were mapped; 474 were microsatellite loci; a few were additional bands produced during PCRs, one of the two rDNA loci (using ITS), the MDH locus, and three sex-linked markers (Q and FB loci and one RAPD band). Twenty-four linkage groups were estimated of which 5 were minute (between 7.1 and 22.8 cM) and 19 were major groups (>76.5 cM). The number of major linkage groups exceeded by three the number of chromosomes of the complement (n = 16). The sum of the lengths of all linkage groups amounts to 4061 cM to which must be added at least 320 cM to link groups in excess, making a total of at least 4381 cM. The length of the largest linkage group I was 630 cM. The average density of markers was 7.5 cM and the average resolution was about one marker every 300 kb. For most of the large groups, the centromeric region was determined genetically, as described in (accompanying article in this issue), using half-tetrad analysis of thelytokous parthenogens in which diploid restoration occurs through central fusion. Several cases of segregation distortion that appreared to result from deleterious recessives were discovered. A low positive interference was also detected.

Mitochondrial DNA in the Drosophila melanogaster complex.

Mitochondrial DNA in the complex Drosophila melanogaster was among the first studied in metazoans. The variability of the molecule was extensively studied using restriction enzymes, gene sequence and recently sequence of the whole coding region. Within the complex, seven major haplotypes have been described, one (me) in D. melanogaster, three in D. simulans (siI, siII, siIII), two in D. mauritiana (maI, maII), and one in D. sechellia (se). The molecular distance between the haplotypes is comprised between 1 and 5%, except for siII and maI, which are virtually identical. The nucleotide diversity within each of these haplotypes is very low, varying from 0 to 0.0005. Most of the cytoplasms are infected by the bacterium Wolbachia and different bacterial strains infect cytoplasms harboring different mtDNA types. mtDNA polymorphism is discussed in relation with Wolbachia, nuclear polymorphism and speciation events.

Mating system in mouse lemurs: theories and facts, using analysis of paternity.

Among nocturnal Malagasy prosimians, the grey mouse lemur (Microcebus murinus) is considered a solitary species which has a promiscuous mating system. Indirect indicators, such as the lack of sexual dimorphism, the overlapping of male and female home ranges with each other, the synchronism of seasonal oestrus and the high relative testes size of males, support the presence of sperm competition. In captivity, an intense sexual precopulatory competition develops among males, leading to the emergence of a dominant male who fathers the majority of the litters. Although multiple mating did occur, the dominant male achieved the majority of the matings on the first day of oestrus. A 'mate-guarding' behaviour, exhibited by the dominant male, was observed in 11 groups out of 15, on only the first day of the vaginal opening and was significantly more often directed towards younger females. Females also played an important role in sexual competition among males, since their presence enhances the aggressive interactions between males. Difference in aggressive behaviours of females, in response to male sexual solicitations, suggests female mate choice. Compared with data from wild animals, it may be hypothesised that alternative mating strategies can be used by male grey mouse lemurs to enhance their reproductive success, depending on the availability and distribution of receptive females.

THE ORIGIN OF WEST EUROPEAN SUBSPECIES OF HONEYBEES (APIS MELLIFERA): NEW INSIGHTS FROM MICROSATELLITE AND MITOCHONDRIAL DATA.

Apis mellifera is composed of three evolutionary branches including mainly African (branch A), western and northern European (branch M), and southeastern European (branch C) populations. The existence of morphological clines extending from the equator to the Polar Circle through Morocco and Spain raised the hypothesis that the branch M originated in Africa. Mitochondrial DNA analysis revealed that branches A and M were characterized by highly diverged lineages implying very remote links between both branches. It also revealed that mtDNA haplotypes from lineages A coexisted with haplotypes M in the Iberian Peninsula and formed a south-north frequency cline, suggesting that this area could be a secondary contact zone between the two branches. By analyzing 11 populations sampled along a France-Spain/Portugal-Morocco-Guinea transect at 8 microsatellite loci and the DraI RFLP of the COI-COII mtDNA marker, we show that Iberian populations do not present any trace of "africanization" and are very similar to French populations when considering microsatellite markers. Therefore, the Iberian Peninsula is not a transition area. The higher haplotype A variability observed in Spanish and Portuguese samples compared to that found in Africa is explained by a higher mutation rate and multiple and recent introductions. Selection appears to be the best explanation to the morphological and allozymic clines and to the diffusion and maintenance of African haplotypes in Spain and Portugal.

EXPERIMENTAL EVIDENCE FOR MITOCHONDRIAL DNA INTROGRESSION BETWEEN DROSOPHILA SPECIES.

Differential introgression of mitochondrial genomes has been used to explain the occurrence in some species of individuals bearing mtDNA from a related species. This situation has been observed for Drosophila mauritiana (endemic to Mauritius) where a high proportion of individuals (88%) carries an mtDNA also found in D. simulans populations from Madagascar and Réunion. Using these two species, experiments were carried out to test for differential mtDNA introgression. A single virgin female from one species (initial frequency 0.03) was introduced into a population of the other. D. simulans mtDNA can, within three generations, almost entirely displace (frequency up to 0.80) D. mauritiana mtDNA. Hybrid male sterility probably curtails to a large degree parallel introgression of nuclear genes. The progress of cytoplasmic introgression is dependent on the degree of inbreeding of the recipient D. mauritiana strains. In reciprocal experiments, introgression was much less likely: few D. mauritiana migrant females are inseminated and their mtDNA frequency always remains very low. The results of these experiments support the hypothesis that a selective advantage of hybrids (probably at the nuclear level) has promoted mtDNA transfer from D. simulans Madagascar or Reunion populations into D. mauritiana through introgressive hybridization.

RACE FORMATION, SPECIATION, AND INTROGRESSION WITHIN DROSOPHILA SIMULANS, D. MAURITIANA, AND D. SECHELLIA INFERRED FROM MITOCHONDRIAL DNA ANALYSIS.

Mitochondrial DNA cleavage maps from three chromosomally homosequential species Drosophila simulans, D. mauritiana, and D. sechellia, were established for 12 restriction enzymes. One isofemale strain was studied in D. sechellia (se), 13 in D. simulans, and 17 in D. mauritiana: in the last two species, respectively, three (siI, II, and III) and two (maI and II) cleavage morphs were found. The evolutionary relationships based on mtDNA cleavage map comparisons show that the maI and se mtDNAs are internal branches of the phylogenetic tree of the D. simulans mtDNA. D. mauritiana and D. sechellia species appear to be derived from a population of D. simulans which carried an ancestral form of the current siI mtDNA type. In addition, two cleavage morphs (siIII [only present in D. simulans from Madagascar] and maI) appeared to be identical, although found in different species. We present a speculative interpretation of data on biogeography and hybridization which is consistent with the hypothesis of a recent introgression of mitochondrial DNA of D. simulans from Madagascar into D. mauritiana.

Demographic aspects of interspecific hybridization : A study of the Jaera albifrons species complex (Crustacea, Isopoda, Asellota).

A demographic study was carried out on two closely related species of the isopods Jaera (albifrons) ischiosetosa and J. (a.) albifrons and their F 1 hybrids. The results from the 16 possible combinations of crosses have permitted an analysis of the nature of the mechanisms assuring the isolation of the species studied. Although intraspecific crosses yield an immediate success, interspecific crosses in the absence of choice of mates progress only slowly during the course of weeks. The results of both crosses between hybrids and back-crosses turn out to be intermediate between those of intra- and interspecific crosses. The hybrids of the first generation are perfectly viable and their survival curves are identical to those of the parents. The fertility of parents in intra- and interspecific crosses is comparable, with the exception of the fragility of female descendants (heterogametic sex) in one direction of crossing. The fertility of the F 1 hybrids, however, crossed either among themselves or with their parents, is quite noticeably decreased: the time needed to double the size of the population is 2.5 times longer for the hybrids than for the parents. This hybrid breakdown completes the pre-fertilization isolating mechanisms: partial ecological isolation, and especially ethological isolation, is practically total when a choice of mates exists. The two species studied, for which demographic parameters are quite close, were raised together for ten generations and yielded only exceptional hybrids with a frequency which does not exceed that found under natural conditions.