An amplified fragment length polymorphism map of the silkworm

First Glimpse of Gut Microbiota of Quarantine Insects in China

Quarantine insects are economically important pests that frequently invade new habitats. A rapid and accurate monitoring method to trace the geographical sources of invaders is required for their prevention, detection, and eradication. Current methods based on genetics are typically time-consuming. Here, we developed a novel tracing method based on insect gut microbiota. The source location of the insect gut microbiota can be used to rapidly determine the geographical origin of the insect. We analyzed 179 gut microbiota samples from 591 individuals of 22 quarantine insect species collected from 36 regions in China. The gut microbiota of these insects primarily included Actinobacteria, Bacteroidetes, Cyanobacteria, Firmicutes, Proteobacteria, and Tenericutes. The diversity of the insect gut microbiota was closely associated with geographical and environmental factors. Different insect species could be distinguished based on the composition of gut microbiota at the phylum level. Populations of individual insect species from different regions could be distinguished based on the composition of gut microbiota at the phylum, class, and order levels. A method for determining the geographical origins of invasive insect species has been established; however, its practical application requires further investigations before implementation.

Comparison of Microbial Communities in Colorado Potato Beetles (Leptinotarsa decemlineata Say) Collected From Different Sources in China

Microbial communities in insects are related to their geographical sources and contribute to adaptation to the local habitat. The Colorado potato beetle (Leptinotarsa decemlineata) (CPB) is a potato pest that causes serious economic losses in Xinjiang Uygur Autonomous Region (XJ) and Heilongjiang Province (HL), China. The influence of microorganisms in the invasion and dispersal of CPB is unclear. We studied microbial communities of CPB collected from nine geographic sources in China using high throughput sequencing technology. Bacteroidetes, Firmicutes, and Proteobacteria were the most dominant phyla, Clostridia, Bacteroidetes, and γ-Proteobacteria were the most dominant classes, Enterobacterales, Lactobacillales, Clostridiales, and Bacteroidales were the most dominant orders, and Enterobacteriaceae, Streptococcidae, Verrucomicrobiaceae, and Rikenellaceae were the most dominant families. There were significant differences, among sources, in the relative abundance of taxa at the genus level. A total of 383 genera were identified, and the dominant bacteria at the genus level were compared between XJ and HL. Pseudomonas was the unique dominant microorganism in the HL area, and the other four microorganisms (Lelliottia, Enterococcus, Enterobacter, and Lactococcus) were common within the 2 regions. Bacterial community diversity in CPB from Urumqi, Jimunai, and Wenquan was higher than diversity in other regions. T-Distributed Stochastic Neighbor Embedding (tSNE) analysis indicated that order and genus were appropriate taxonomic levels to distinguish geographical sources of CPB. These findings provide insight into the diversity of microorganisms of CPB in the differences among geographically isolated populations.

Construction of a high-density genetic linkage map and mapping of quantitative trait loci for growth-related traits in silver carp (Hypophthalmichthys molitrix)

High-density genetic map and quantitative trait loci (QTL) mapping are powerful tools for identifying genomic regions that may be responsible for such polygenic trait as growth. A high-density genetic linkage map was constructed by sequencing 198 individuals in a F₁ family of silver carp (Hypophthalmichthys molitrix) in this study. This genetic map spans a length of 2,721.07 cM with 3,134 SNPs distributed on 24 linkage groups (LGs). Comparative genomic mapping presented a high level of syntenic relationship between silver carp and zebrafish. We detected one major and nineteen suggestive QTL for 4 growth-related traits (body length, body height, head length and body weight) at 6, 12 and 18 months post hatch (mph), explaining 10.2~19.5% of phenotypic variation. All six QTL for growth traits of 12 mph generally overlapped with QTL for 6 mph, while the majority of QTL for 18 mph were identified on two additional LGs, which may reveal a different genetic modulation during early and late muscle growth stages. Four potential candidate genes were identified from the QTL regions by homology searching of marker sequences against zebrafish genome. Hepcidin, a potential candidate gene identified from a QTL interval on LG16, was significantly associated with growth traits in the analyses of both phenotype-SNP association and mRNA expression between small-size and large-size groups of silver carp. These results provide a basis for elucidating the genetic mechanisms for growth and body formation in silver carp, a world aquaculture fish.

Genetic variation in ecoraces of tropical tasar silkworm, Antheraea mylitta using SSR markers

The tropical tasar silkworm, Antheraea mylitta, polyphagous sericigenous insect mostly found in the tropical areas of India. It is found in these regions as ecotypes or ecoraces. It feeds primarily on plants, a variety of secondary plants like Terminalia arjuna and T. tomentosa. Tasar culture is a traditional livelihood for lakhs of tribal populace in the areas of Jharkhand, Chhatisgarh, Orissa, Maharashtra, Andhra Pradesh, West Bengal and Uttar Pradesh. In the present study, the genetic diversity of these ecoraces is identified by DNA markers, namely simple sequence repeats (SSRs), most of which produced polymorphic bands.

A High-Density SNP Genetic Linkage Map and QTL Analysis of Growth-Related Traits in a Hybrid Family of Oysters (Crassostrea gigas × Crassostrea angulata) Using Genotyping-by-Sequencing

Oysters are among the most important species in global aquaculture. Crassostrea gigas, and its subspecies C. angulata, are the major cultured species. To determine the genetic basis of growth-related traits in oysters, we constructed a second-generation linkage map from 3367 single-nucleotide polymorphisms (SNPs) based on genotyping-by-sequencing, genotyped from a C. gigas × C. angulata hybrid family. These 3367 SNPs were distributed on 1695 markers, which were assigned to 10 linkage groups. The genetic linkage map had a total length of 1084.3 cM, with an average of 0.8 cM between markers; it thus represents the densest genetic map constructed for oysters to date. Twenty-seven quantitative trait loci (QTL) for five growth-related traits were detected. These QTL could explain 4.2-7.7% (mean = 5.4%) of the phenotypic variation. In total, 50.8% of phenotypic variance for shell width, 7.7% for mass weight, and 34.1% for soft tissue weight were explained. The detected QTL were distributed among eight linkage groups, and more than half (16) were concentrated within narrow regions in their respective linkage groups. Thirty-eight annotated genes were identified within the QTL regions, two of which are key genes for carbohydrate metabolism. Other genes were found to participate in assembly and regulation of the actin cytoskeleton, signal transduction, and regulation of cell differentiation and development. The newly developed high-density genetic map, and the QTL and candidate genes identified provide a valuable genetic resource and a basis for marker-assisted selection for C. gigas and C. angulata.

Genetic mapping and QTL analysis of growth traits in the large yellow croaker Larimichthys crocea

Large yellow croaker (Larimichthys crocea) is an important maricultured species in China. A genetic linkage map of the large yellow croaker was constructed using type II microsatellites and expressed sequence tag (EST)-derived microsatellites in two half-sib families (two females and one male). A total of 289 microsatellite markers (contained 93 EST-SSRs) were integrated into 24 linkage groups, which agreed with the haploid chromosome number. The map spanned a length of 1,430.8 cm with an average interval of 5.4 cm, covering 83.9 % of the estimated genome size (1,704.8 cm). A total of seven quantitative trait locis (QTLs) were detected for growth traits on five linkage groups, including two 1 % and five 5 % chromosome-wide significant QTLs, and explained from 2.33 to 5.31 % of the trait variation. The identified QTLs can be applied in marker-assisted selection programs to improve the growth traits.

Evidence for genetic incompatibilities associated with post-zygotic reproductive isolation in the human fungal pathogen Cryptococcus neoformans

Hybridization is a potent mechanism for generating unique strains with broad host ranges and increased virulence in fungal pathogens. In the opportunistic basidiomycete pathogen Cryptococcus neoformans, intervarietal hybrids are commonly found infecting patients. The two parental varieties C. neoformans var. grubii and C. neoformans var. neoformans mate readily under laboratory conditions, but the hybrid basidiospores have germination rates about four times lower than those from intravarietal crosses. Here, we used microdissection to collect basidiospores from a hybrid cross and analysed the genotypes of germinated basidiospores to identify potentially antagonistic allelic combinations between loci that impact basidiospore germination. Our analyses showed clear evidence for Bateson–Dobzhansky–Muller (BDM) incompatibility affecting basidiospore viability. Antagonistic combinations of alleles from both two loci and three loci were found. Interestingly, most of the hybrid progeny showed segregation distortion in favour of the alleles from var. neoformans, consistent with large-scale epistatic interactions among loci affecting basidiospore viability. Our study presents the first evidence of BDM incompatibility between nuclear genes affecting post-zygotic reproductive isolation in this model basidiomycete yeast.

Genetic characterization of Bombyx mori (Lepidoptera: Bombycidae) breeding and hybrid lines with different geographic origins

The domesticated silkworm Bombyx mori L. comprises a large number of geographical breeds and hybrid lines. Knowing the genetic structure of those may provide information to improve the conservation of commercial lines by estimating inbreeding over generations and the consequences of excessive use of those lineages. Here, we analyzed the genetic diversity of seven breeds and eight hybrid lines from Eastern Europe and Asia using highly polymorphic microsatellites markers to determine its genetical impact on their use in global breeding programs. No consistent pattern of deviation from Hardy-Weinberg equilibrium was found for most breed and hybrids; and the absence of a linkage disequilibrium also suggests that the strains are in equilibrium. A principal coordinate analysis revealed a clear separation of two silkworm breeds from the rest: one (IBV) originated from India and the other one (RG90) from Romania/Japan. The tendency of the other breeds from different geographic origins to cluster together in a general mix might be due to similar selection pressures (climate and anthropogenic factors) in different geographic locations. Phylogenetic analyses grouped the different silkworm breeds but not the hybrids according to their geographic origin and confirmed the pattern found in the principal coordinate analysis.

Genetic analysis and QTL mapping of seed coat color in sesame (Sesamum indicum L.)

Seed coat color is an important agronomic trait in sesame, as it is associated with seed biochemical properties, antioxidant content and activity and even disease resistance of sesame. Here, using a high-density linkage map, we analyzed genetic segregation and quantitative trait loci (QTL) for sesame seed coat color in six generations (P1, P2, F1, BC1, BC2 and F2). Results showed that two major genes with additive-dominant-epistatic effects and polygenes with additive-dominant-epistatic effects were responsible for controlling the seed coat color trait. Average heritability of the major genes in the BC1, BC2 and F2 populations was 89.30%, 24.00%, and 91.11% respectively, while the heritability of polygenes was low in the BC1 (5.43%), in BC2 (0.00%) and in F2 (0.89%) populations. A high-density map was constructed using 724 polymorphic markers. 653 SSR, AFLP and RSAMPL loci were anchored in 14 linkage groups (LG) spanning a total of 1,216.00 cM. The average length of each LG was 86.86 cM and the marker density was 1.86 cM per marker interval. Four QTLs for seed coat color, QTL1-1, QTL11-1, QTL11-2 and QTL13-1, whose heritability ranged from 59.33%-69.89%, were detected in F3 populations using CIM and MCIM methods. Alleles at all QTLs from the black-seeded parent tended to increase the seed coat color. Results from QTLs mapping and classical genetic analysis among the P1, P2, F1, BC1, BC2 and F2 populations were comparatively consistent. This first QTL analysis and high-density genetic linkage map for sesame provided a good foundation for further research on sesame genetics and molecular marker-assisted selection (MAS).

A simple sequence repeat- and single-nucleotide polymorphism-based genetic linkage map of the brown planthopper, Nilaparvata lugens

In this study, we developed the first genetic linkage map for the major rice insect pest, the brown planthopper (BPH, Nilaparvata lugens). The linkage map was constructed by integrating linkage data from two backcross populations derived from three inbred BPH strains. The consensus map consists of 474 simple sequence repeats, 43 single-nucleotide polymorphisms, and 1 sequence-tagged site, for a total of 518 markers at 472 unique positions in 17 linkage groups. The linkage groups cover 1093.9 cM, with an average distance of 2.3 cM between loci. The average number of marker loci per linkage group was 27.8. The sex-linkage group was identified by exploiting X-linked and Y-specific markers. Our linkage map and the newly developed markers used to create it constitute an essential resource and a useful framework for future genetic analyses in BPH.

Genetic mapping and QTL analysis of growth-related traits in the Pacific oyster

The Pacific oyster (Crassostrea gigas) is one of the most important oysters cultured worldwide. To analyze the oyster genome and dissect growth-related traits, we constructed a sex-averaged linkage map by combining 64 genomic simple sequence repeats, 42 expressed sequence tag-derived SSRs, and 320 amplified fragment length polymorphism markers in an F(1) full-sib family. A total of 426 markers were assigned to 11 linkage groups, spanning 558.2 cM with an average interval of 1.3 cM and 94.7% of genome coverage. Segregation distortion was significant for 18.8% of the markers (P < 0.05), and distorted markers tended to occur on some genetic regions or linkage groups. Most growth-related quantitative traits were highly significantly (P < 0.01) correlated, and principal component analysis obtained four principal components. Quantitative trait locus (QTL) analysis identified three significant QTLs for two principal components, which explained 0.6-13.9% of the phenotypic variation. One QTL for sex was detected on linkage group 6, and the inheritabilities of sex for parental alleles and maternal alleles on that locus C15 are 39.8% and 0.01%, respectively. The constructed linkage map and determined QTLs can provide a tool for further genetic analysis of the traits and be potential for marker-assisted selection in C. gigas breeding.

Genome scans for transmission ratio distortion regions in mice

Transmission ratio distortion (TRD) is the departure from the expected genotypic frequencies under Mendelian inheritance. This departure can be due to multiple physiological mechanisms during gametogenesis, fertilization, fetal and embryonic development, and early neonatal life. Although a few TRD loci have been reported in mouse, inheritance patterns have never been evaluated for TRD. In this article, we developed a Bayesian binomial model accounting for additive and dominant deviation TRD mechanisms. Moreover, this model was used to perform genome-wide scans for TRD quantitative trait loci (QTL) on six F2 mouse crosses involving between 296 and 541 mice and between 72 and 1854 genetic markers. Statistical significance of each model was checked at each genetic marker with Bayes factors. Genome scans revealed overdominance TRD QTL located in mouse chromosomes 1, 2, 12, 13, and 14 and additive TRD QTL in mouse chromosomes 2, 3, and 15, although these results did not replicate across mouse crosses. This research contributes new statistical tools for the analysis of specific genetic patterns involved in TRD in F2 populations, our results suggesting a relevant incidence of TRD phenomena in mouse with important implications for both statistical analyses and biological research.

A new mapping method for quantitative trait loci of silkworm

BACKGROUND

Silkworm is the basis of sericultural industry and the model organism in insect genetics study. Mapping quantitative trait loci (QTLs) underlying economically important traits of silkworm is of high significance for promoting the silkworm molecular breeding and advancing our knowledge on genetic architecture of the Lepidoptera. Yet, the currently used mapping methods are not well suitable for silkworm, because of ignoring the recombination difference in meiosis between two sexes.

RESULTS

A mixed linear model including QTL main effects, epistatic effects, and QTL × sex interaction effects was proposed for mapping QTLs in an F2 population of silkworm. The number and positions of QTLs were determined by F-test and model selection. The Markov chain Monte Carlo (MCMC) algorithm was employed to estimate and test genetic effects of QTLs and QTL × sex interaction effects. The effectiveness of the model and statistical method was validated by a series of simulations. The results indicate that when markers are distributed sparsely on chromosomes, our method will substantially improve estimation accuracy as compared to the normal chiasmate F2 model. We also found that a sample size of hundreds was sufficiently large to unbiasedly estimate all the four types of epistases (i.e., additive-additive, additive-dominance, dominance-additive, and dominance-dominance) when the paired QTLs reside on different chromosomes in silkworm.

CONCLUSION

The proposed method could accurately estimate not only the additive, dominance and digenic epistatic effects but also their interaction effects with sex, correcting the potential bias and precision loss in the current QTL mapping practice of silkworm and thus representing an important addition to the arsenal of QTL mapping tools.

Identification of 2chromosome region translocated onto the W chromosome by RFLP with EST-cDNA clones in the Gensei-kouken strains of the mulberry silkworm, Bombyx mori L

In silkworms, sex-limited strains are either obtained spontaneously or induced by X-rays or gamma rays. When a fragment of an autosome carrying a dominant allele of those genes responsible for certain characters is translocated onto a W chromosome, the female of the successive generations will express these phenotypic characters and sex discrimination can be facilitated. Gensei-kouken strains are sex-limited strains of silkworms developed by irradiating the pupae with gamma rays, by which a portion of the second chromosome is translocated onto the W chromosome. In these improved strains, the females are yellow-blooded and spin yellow cocoons. By using the EST-cDNA clones mapped on the Z chromosome, we identified the sex according to the polymorphic banding pattern or intensity of the signals. Furthermore, by using the clones on the second chromosome, the region of the second chromosome translocated onto the W chromosome was also defined. In both the A95 and A 96 strains selected for the present study, only the mid-portion of the second chromosome was translocated. The differences in length of the fragments translocated in these strains are discussed.

A comparison of recombination frequencies in intraspecific versus interspecific mapping populations of Nasonia

We present the first intraspecific linkage map for Nasonia vitripennis based on molecular markers. The map consists of 36 new microsatellite markers, extracted from the Nasonia genome sequence, and spans 515 cM. The five inferred linkage groups correspond to the five chromosomes of Nasonia. Comparison of recombination frequencies of the marker intervals spread over the whole genome (N=33 marker intervals) between the intraspecific N. vitripennis map and an interspecific N. vitripennis x N. giraulti map revealed a slightly higher (1.8%) recombination frequency in the intraspecific cross. We further considered an N. vitripennis x N. longicornis map with 29 microsatellite markers spanning 430 cM. Recombination frequencies in the two interspecific crosses differed neither between reciprocal crosses nor between mapping populations of embryos and adults. No major chromosomal rearrangements were found for the analyzed genomic segments. The observed differential F(2) hybrid male mortality has no significant effect on the genome-wide recombination frequency in Nasonia. We conclude that interspecific crosses between the different Nasonia species, a hallmark of Nasonia genetics, are generally suitable for mapping quantitative and qualitative trait loci for species differences.

Amplified fragment length polymorphism mapping of quantitative trait loci for economically important traits in the silkworm, Bombyx mori

Cocoon related characteristics are economically important traits in the silkworm, Bombyx mori L. (Lepidoptera: Bombycidae). In this study a genetic linkage map was developed that identified QTL controlling the cocoon weight, cocoon shell weight, and cocoon shell percentage using 161 amplified fragment length polymorphism (AFLP) markers. Twenty PstI/TaqI primer combinations were employed to genotype 78 F(2) progenies derived from a cross between P107 Japanese inbred line and Khorasan Lemon Iranian native strain. Among polymorphic markers, 159 AFLP markers were assigned to 24 linkage groups at the LOD threshold of 2.5 that varied in length from 4 to 299 cM. The total length of the linkage map was 2747 cM, giving an average marker resolution of 19.31 cM. A total of 21 AFLP markers were identified that were distributed over the ten linkage groups linked to the three studied traits using the composite interval mapping method. The explained variation rate by QTL controlling cocoon weight, cocoon shell weight, and cocoon shell percentage ranged from 0.02% to 64.85%, 0.2% to 49.11%, and 0.04% to 84.20%, respectively. These QTL controlled by different actions as well as under dominance, additive, partial dominance, dominance, and over dominance.

Construction of a genetic linkage map and mapping of a female-specific DNA marker in half-smooth tongue sole (Cynoglossus semilaevis)

The half-smooth tongue sole (Cynoglossus semilaevis, hereafter, "tongue sole") is a marine flatfish with great commercial importance for fisheries and aquaculture in China. It has also been a promising model for the study of sex determination mechanisms in fish. Here, we report the construction of a genetic linkage map for the tongue sole, based on 137 markers including 103 AFLP markers, 33 microsatellite markers, and one female-specific DNA marker. Twenty-six linkage groups (LGs) were found. The total map length was 934.6 cM (Kosambi), with an average spacing of 8.4 cM, covering 64.4% of the estimated genome size. Furthermore, a female-specific SCAR marker, CseF-382, was mapped on LG5. This study represents the first genetic linkage map in the tongue sole. This map has great potential in the identification of quantitative traits loci and sex-related genes and marker-assisted selection in the tongue sole. Meanwhile, the new set of polymorphic microsatellite markers developed in this study is not only useful for genetic mapping but also of critical importance for studies on genetic diversity and broodstock management in tongue sole.

An integrated genetic linkage map for silkworms with three parental combinations and its application to the mapping of single genes and QTL

Background

Bombyx mori, the domesticated silkworm, is a well-studied model insect with great economic and scientific significance. Although more than 400 mutations have been described in silkworms, most have not been identified, especially those affecting economically-important traits. Simple sequence repeats (SSRs) are effective and economical tools for mapping traits and genetic improvement. The current SSR linkage map is of low density and contains few polymorphisms. The purpose of this work was to develop a dense and informative linkage map that would assist in the preliminary mapping and dissection of quantitative trait loci (QTL) in a variety of silkworm strains.

Results

Through an analysis of > 50,000 genotypes across new mapping populations, we constructed two new linkage maps covering 27 assigned chromosomes and merged the data with previously reported data sets. The integrated consensus map contains 692 unique SSR sites, improving the density from 6.3 cM in the previous map to 4.8 cM. We also developed 497 confirmed neighboring markers for corresponding low-polymorphism sites, with 244 having polymorphisms. Large-scale statistics on the SSR type were suggestive of highly efficient markers, based upon which we searched 16,462 available genomic scaffolds for SSR loci. With the newly constructed map, we mapped single-gene traits, the QTL of filaments, and a number of ribosomal protein genes.

Conclusion

The integrated map produced in this study is a highly efficient genetic tool for the high-throughput mapping of single genes and QTL. Compared to previous maps, the current map offers a greater number of markers and polymorphisms; thus, it may be used as a resource for marker-assisted breeding.

Cytogenetic characterization and AFLP-based genetic linkage mapping for the butterfly Bicyclus anynana, covering all 28 karyotyped chromosomes

Background

The chromosome characteristics of the butterfly Bicyclus anynana, have received little attention, despite the scientific importance of this species. This study presents the characterization of chromosomes in this species by means of cytogenetic analysis and linkage mapping.

Methodology/Principal Findings

Physical genomic features in the butterfly B. anynana were examined by karyotype analysis and construction of a linkage map. Lepidoptera possess a female heterogametic W-Z sex chromosome system. The WZ-bivalent in pachytene oocytes of B. anynana consists of an abnormally small, heterochromatic W-chromosome with the Z-chromosome wrapped around it. Accordingly, the W-body in interphase nuclei is much smaller than usual in Lepidoptera. This suggests an intermediate stage in the process of secondary loss of the W-chromosome to a ZZ/Z sex determination system. Two nucleoli are present in the pachytene stage associated with an autosome and the WZ-bivalent respectively. Chromosome counts confirmed a haploid number of n = 28. Linkage mapping had to take account of absence of crossing-over in females, and of our use of a full-sib crossing design. We developed a new method to determine and exclude the non-recombinant uninformative female inherited component in offspring. The linkage map was constructed using a novel approach that uses exclusively JOINMAP-software for Lepidoptera linkage mapping. This approach simplifies the mapping procedure, avoids over-estimation of mapping distance and increases the reliability of relative marker positions. A total of 347 AFLP markers, 9 microsatellites and one single-copy nuclear gene covered all 28 chromosomes, with a mapping distance of 1354 cM. Conserved synteny of Tpi on the Z-chromosome in Lepidoptera was confirmed for B. anynana. The results are discussed in relation to other mapping studies in Lepidoptera.

Conclusions/Significance

This study adds to the knowledge of chromosome structure and evolution of an intensively studied organism. On a broader scale it provides an insight in Lepidoptera sex chromosome evolution and it proposes a simpler and more reliable method of linkage mapping than used for Lepidoptera to date.

Inheritance and linkage analysis of co-dominant SSR markers on the Z chromosome of the silkworm (Bombyx mori L.)

Microsatellites or simple sequence repeats (SSRs) are co-dominant molecular markers. When we used fluorescent SSR markers to construct a linkage map for the female heterogametic silkworm (Bombyx mori, ZW), we found that some loci did not segregate in a Mendelian ratio of 1:1 in a backcross population. These loci segregated in a 3:1 ratio of single bands compared with double bands. Further examination of band patterns indicated that three types of SSR bands were present: two homozygotes and one heterozygote. In the beginning, we considered to discard these markers. By scoring male and female F1 individuals, we confirmed that these loci were located on the Z chromosome. Using the sex-linked visible mutation sch (K05) and its wild-type (C108), we constructed an F1 male backcross (BC1M) mapping population. The combination of sch backcross and SSR data enabled us to map the SSR markers to the Z chromosome. By adjusting input parameters based on these data, we were able to use Mapmaker software to construct a linkage map. This strategy takes advantage of co-dominant markers for positional cloning of genes on the Z chromosome. We localized sch to the Z chromosome relative to six SSR markers and one PCR marker, covering a total of 76.1 cM. The sch mutation is an important sex-linked visible mutation widely used in breeding of commercial silkworms (e.g. male silkworm selection rearing). Localization of the sch gene may prove helpful in cloning the gene and developing strains for marker-assisted selection in silkworm breeding.

From genome to proteome: great progress in the domesticated silkworm (Bombyx mori L.)

As the only truly domesticated insect, the silkworm not only has great economic value, but it also has value as a model for genetics and molecular biology research. Genomics and proteomics have recently shown vast potential to be essential tools in domesticated silkworm research, especially after the completion of the Bombyx mori genome sequence. This paper reviews the progress of the domesticated silkworm genome, particularly focusing on its genetic map, physical map and functional genome. This review also presents proteomics, the proteomic technique and its application in silkworm research.

A BAC-based integrated linkage map of the silkworm Bombyx mori

An integrated map of the Bombyx mori genome has been constructed using 361.1 Mb of BAC contigs and singletons together with a genetic map containing 1689 independent genes and synteny among Apis, Tribolium, and Bombyx was examined.

Background

In 2004, draft sequences of the model lepidopteran Bombyx mori were reported using whole-genome shotgun sequencing. Because of relatively shallow genome coverage, the silkworm genome remains fragmented, hampering annotation and comparative genome studies. For a more complete genome analysis, we developed extended scaffolds combining physical maps with improved genetic maps.

Results

We mapped 1,755 single nucleotide polymorphism (SNP) markers from bacterial artificial chromosome (BAC) end sequences onto 28 linkage groups using a recombining male backcross population, yielding an average inter-SNP distance of 0.81 cM (about 270 kilobases). We constructed 6,221 contigs by fingerprinting clones from three BAC libraries digested with different restriction enzymes, and assigned a total of 724 single copy genes to them by BLAST (basic local alignment search tool) search of the BAC end sequences and high-density BAC filter hybridization using expressed sequence tags as probes. We assigned 964 additional expressed sequence tags to linkage groups by restriction fragment length polymorphism analysis of a nonrecombining female backcross population. Altogether, 361.1 megabases of BAC contigs and singletons were integrated with a map containing 1,688 independent genes. A test of synteny using Oxford grid analysis with more than 500 silkworm genes revealed six versus 20 silkworm linkage groups containing eight or more orthologs of Apis versus Tribolium, respectively.

Conclusion

The integrated map contains approximately 10% of predicted silkworm genes and has an estimated 76% genome coverage by BACs. This provides a new resource for improved assembly of whole-genome shotgun data, gene annotation and positional cloning, and will serve as a platform for comparative genomics and gene discovery in Lepidoptera and other insects.

Genetic mapping of a food preference gene in the silkworm, Bombyx mori, using restriction fragment length polymorphisms (RFLPs)

The domesticated silkworm, Bombyx mori, has strict food preferences and grows by feeding on mulberry leaves. However, "Sawa-J", an abnormal feeding habit strain selected from the genetic stock, feeds on an artificial diet without mulberry leaf powder. In this study, the food preference gene in Sawa-J was genetically identified using restriction fragment length polymorphisms (RFLPs) of a cDNA clone on each linkage group. Taking advantage of a lack of genetic recombination in females, reciprocal backcrossed F1 (BC1) progenies were independently prepared using a non-feeding strain, C108, as a mating partner of Sawa-J. Our results of linkage analysis and mapping proved that the feeding behavior is primarily controlled by a major recessive gene mapped at 20.2 cM on RFLP linkage group 9 (RFLG9), and clone e73 at a distance of 4.2 cM was found as the first linked molecular marker.

Identification of the female-determining region of the W chromosome in Bombyx mori

The W chromosome of the silkworm Bombyx mori is devoid of functional genes, except for the putative female-determining gene (Fem). To localize Fem, we investigated the presence of W-specific DNA markers on strains in which an autosomal fragment containing dominant marker genes was attached to the W chromosome. We produced new W-chromosomal fragments from the existing Zebra-W strain (T(W;3)Ze chromosome) by X-irradiation, and then carried out deletion mapping of these and sex-limited yellow cocoon strains (T(W;2)Y-Chu, -Abe and -Ban types) from different Japanese stock centers. Of 12 RAPD markers identified in the normal W chromosomes of most silkworm strains in Japan, the newly irradiated W(B-YL-YS)Ze chromosome contained three, the T(W;2)Y-Chu chromosome contained six, and the T(W;2)Y-Abe and -Ban chromosomes contained only one (W-Rikishi). To investigate the ability of the reduced W-chromosome translocation fragments to form heterochromatin bodies, which are found in nuclei of normal adult female sucking stomachs, we examined cells of the normal type p50 strain and the T(W;2)Y-Chu and -Abe strains. A single sex heterochromatin body was found in nuclei of p50 females, whereas we detected only small sex heterochromatin bodies in the T(W;2)Y-Chu strain and no sex heterochromatin body in the T(W;2)Y-Abe strain. Since adult females of all strains were normal and fertile, we conclude that only extremely limited region, containing the W-Rikishi RAPD sequence of the W chromosome, is required to determine femaleness. Based on a comparison of the normal W-chromosome and 7 translocation and W-deletion strains we present a map of Fem relative to the 12 W-specific RAPD markers.

Development of silver carp (Hypophthalmichthys molitrix) and bighead carp (Aristichthys nobilis) genetic maps using microsatellite and AFLP markers and a pseudo-testcross strategy

Silver carp (Hypophthalmichthys molitrix) and bighead carp (Aristichthys nobilis) are two of the four most important pond-cultured fish species inhabiting the major river basins of China. In the present study, genetic maps of silver carp and bighead carp were constructed using microsatellite and AFLP markers and a two-way pseudo-testcross strategy. To create the maps, 60 individuals were obtained from a cross of a single bighead carp (female) and a single silver carp (male). The silver carp map consisted of 271 markers (48 microsatellites and 223 AFLPs) that were assembled into 27 linkage groups, of which 22 contained at least four markers. The total length of the silver carp map was 952.2 cM, covering 82.8% of the estimated genome size. The bighead carp map consisted of 153 markers (27 microsatellites and 126 AFLPs) which were organized into 30 linkage groups, of which 19 contained at least four markers. The total length of the bighead carp map was 852.0 cM, covering 70.5% of the estimated genome size. Eighteen microsatellite markers were common to both maps. These maps will contribute to discovery of genes and genetic regions controlling traits in the two species of carp.

Population variation of the fall armyworm, Spodoptera frugiperda, in the Western Hemisphere

Spodoptera frugiperda (J.E. Smith) (Lepidoptera: Noctuidae), the fall armyworm is the most economically important maize pest in the western hemisphere. This research focused on the genetic variability of the maize host strain because there is a lack of information in this area of S. frugiperda research. Amplified fragment length polymorphism (AFLP) was used to assess the genetic variability of S. frugiperda over a large geographic area. Twenty populations were collected from the maize, one population was collected from princess tree, one population was collected from lemon tree, and one population was collected from bermudagrass. The 23 populations were from Mexico, the continental United States, Puerto Rico, Brazil, and Argentina. The objective of this research was to evaluate whether the majority of genetic variability was within populations or between populations. The AFLP results showed that the majority of the genetic variability is within populations and not between populations, indicating minor gene flow and suggesting that S. frugiperda in the Western Hemisphere are an interbreeding population.

High recombination frequency creates genotypic diversity in colonies of the leaf-cutting ant Acromyrmex echinatior

Honeybees are known to have genetically diverse colonies because queens mate with many males and the recombination rate is extremely high. Genetic diversity among social insect workers has been hypothesized to improve general performance of large and complex colonies, but this idea has not been tested in other social insects. Here, we present a linkage map and an estimate of the recombination rate for Acromyrmex echinatior, a leaf-cutting ant that resembles the honeybee in having multiple mating of queens and colonies of approximately the same size. A map of 145 AFLP markers in 22 linkage groups yielded a total recombinational size of 2076 cM and an inferred recombination rate of 161 kb cM(-1) (or 6.2 cM Mb(-1)). This estimate is lower than in the honeybee but, as far as the mapping criteria can be compared, higher than in any other insect mapped so far. Earlier studies on A. echinatior have demonstrated that variation in division of labour and pathogen resistance has a genetic component and that genotypic diversity among workers may thus give colonies of this leaf-cutting ant a functional advantage. The present result is therefore consistent with the hypothesis that complex social life can select for an increased recombination rate through effects on genotypic diversity and colony performance.

A genetic linkage map of Pacific white shrimp (Litopenaeus vannamei): sex-linked microsatellite markers and high recombination rates

Pacific white shrimp (Litopenaeus vannamei) is the leading species farmed in the Western Hemisphere and an economically important aquaculture species in China. In this project, a genetic linkage map was constructed using amplified fragment length polymorphism (AFLP) and microsatellite markers. One hundred and eight select AFLP primer combinations and 30 polymorphic microsatellite markers produced 2071 markers that were polymorphic in either of the parents and segregated in the progeny. Of these segregating markers, 319 were mapped to 45 linkage groups of the female framework map, covering a total of 4134.4 cM; and 267 markers were assigned to 45 linkage groups of the male map, covering a total of 3220.9 cM. High recombination rates were found in both parental maps. A sex-linked microsatellite marker was mapped on the female map with 6.6 cM to sex and a LOD of 17.8, two other microsatellite markers were also linked with both 8.6 cM to sex and LOD score of 14.3 and 16.4. The genetic maps presented here will serve as a basis for the construction of a high-resolution genetic map, quantitative trait loci (QTLs) detection, marker-assisted selection (MAS) and comparative genome mapping.

A microsatellite-based genetic linkage map of the waterflea, Daphnia pulex: On the prospect of crustacean genomics

We describe the first genetic linkage map for Daphnia pulex using 185 microsatellite markers, including 115 new markers reported in this study. Our approach was to study the segregation of polymorphisms in 129 F2 progeny of one F1 hybrid obtained by crossing two genetically divergent lineages of Daphnia isolated from two Oregon populations. The map spanned 1206 Kosambi cM and had an average intermarker distance of 7 cM. Linkage groups ranged in size from 7 to 185 cM and contained 4 to 27 markers. The map revealed 12 linkage groups corresponding to the expected number of chromosomes and covers approximately 87% of the genome. Tests for random segregation of alleles at individual loci revealed that 21% of the markers showed significant transmission ratio distortion (primarily homozygote deficiency) likely due to markers being linked to deleterious recessive alleles. This map will become the anchor for the physical map of the Daphnia genome and will serve as a starting point for mapping single and quantitative trait loci affecting ecologically important phenotypes. By mapping 342 tentative orthologous gene pairs (Daphnia/Drosophila) into the Daphnia linkage map, we facilitate future comparative projects.

Molecular marker systems in insects: current trends and future avenues

Insects comprise the largest species composition in the entire animal kingdom and possess a vast undiscovered genetic diversity and gene pool that can be better explored using molecular marker techniques. Current trends of application of DNA marker techniques in diverse domains of insect ecological studies show that mitochondrial DNA (mtDNA), microsatellites, random amplified polymorphic DNA (RAPD), expressed sequence tags (EST) and amplified fragment length polymorphism (AFLP) markers have contributed significantly for progresses towards understanding genetic basis of insect diversity and for mapping medically and agriculturally important genes and quantitative trait loci in insect pests. Apart from these popular marker systems, other novel approaches including transposon display, sequence-specific amplification polymorphism (S-SAP), repeat-associated polymerase chain reaction (PCR) markers have been identified as alternate marker systems in insect studies. Besides, whole genome microarray and single nucleotide polymorphism (SNP) assays are becoming more popular to screen genome-wide polymorphisms in fast and cost effective manner. However, use of such methodologies has not gained widespread popularity in entomological studies. The current study highlights the recent trends of applications of molecular markers in insect studies and explores the technological advancements in molecular marker tools and modern high throughput genotyping methodologies that may be applied in entomological researches for better understanding of insect ecology at molecular level.

A preliminary genetic linkage map of the Pacific abalone Haliotis discus hannai Ino

Preliminary genetic linkage maps were constructed for the Pacific abalone (Haliotis discus hannai Ino) using amplified fragment length polymorphism (AFLP), randomly amplified polymorphic DNA (RAPD), and microsatellite markers segregating in a F(1) family. Nine microsatellite loci, 41 RAPD, and 2688 AFLP markers were genotyped in the parents and 86 progeny of the mapping family. Among the 2738 markers, 384 (including 365 AFLP markers, 10 RAPD markers, and 9 microsatellite loci) were polymorphic and segregated in one or both parents: 241 in the female and 146 in the male. The majority of these markers, 232 in the female and 134 in the male, segregated according to the expected 1:1 Mendelian ratio (alpha = 0.05). Two genetic linkage maps were constructed using markers segregating in the female or the male parent. The female framework map consisted of 119 markers in 22 linkage groups, covering 1773.6 cM with an average intermarker space of 18.3 cM. The male framework map contained 94 markers in 19 linkage groups, spanning 1365.9 cM with an average intermarker space of 18.2 cM. The sex determination locus was mapped to the male map but not to the female map, suggesting a XY-male determination mechanism. Distorted markers showing excess of homozygotes were mapped in clusters, probably because of their linkage to a gene that is incompatible between two parental populations.

Amplified fragment length polymorphism mapping of quantitative trait loci for malaria parasite susceptibility in the yellow fever mosquito Aedes aegypti

The yellow fever mosquito Aedes aegypti has been the subject of extensive genetic research due to its medical importance and the ease with which it can be manipulated in the laboratory. A molecular genetic linkage map was constructed using 148 amplified fragment length polymorphism (AFLP) and six single-strand conformation polymorphism (SSCP) markers. Eighteen AFLP primer combinations were used to genotype two reciprocal F2 segregating populations. Each primer combination generated an average of 8.2 AFLP markers eligible for linkage mapping. The length of the integrated map was 180.9 cM, giving an average marker resolution of 1.2 cM. Composite interval mapping revealed a total of six QTL significantly affecting Plasmodium susceptibility in the two reciprocal crosses of Ae. aegypti. Two common QTL on linkage group 2 were identified in both crosses that had similar effects on the phenotype, and four QTL were unique to each cross. In one cross, the four main QTL accounted for 64% of the total phenotypic variance, and digenic epistasis explained 11.8% of the variance. In the second cross, the four main QTL explained 66% of the variance, and digenic epistasis accounted for 16% of the variance. The actions of these QTL were either dominance or underdominance. Our results indicated that at least three new QTL were mapped on chromosomes 1 and 3. The polygenic nature of susceptibility to P. gallinaceum and epistasis are important factors for significant variation within or among mosquito strains. The new map provides additional information useful for further genetic investigation, such as identification of new genes and positional cloning.

Structure of worldwide populations of Lasioderma serricorne (Coleoptera: Anobiidae) as revealed by amplified fragment length polymorphism profiles

The cigarette beetle Lasioderma serricorne (Fabricius) is the most widespread and destructive pest of stored tobacco. The capability to differentiate between populations from different geographic origins would enable researchers to better understand how insect dispersal through transportation affects the infestation of stored tobacco. Using amplified fragment length polymorphism (AFLP), DNA polymorphisms were assessed in 16 populations of L. serricorne collected from 15 countries. The dendrograms constructed from profile distance matrices revealed well-supported colony clusters. There was no clear clustering as a function of the geographic origin of the samples. The results suggest extensive insect dispersal among geographical regions due to movement of infested commodities worldwide. This first AFLP population study of a stored-product insect demonstrates the potential of AFLP for distinguishing L. serricorne populations.

A second-generation integrated map of the silkworm reveals synteny and conserved gene order between lepidopteran insects

A second-generation linkage map was constructed for the silkworm, Bombyx mori, focusing on mapping Bombyx sequences appearing in public nucleotide databases and bacterial artificial chromosome (BAC) contigs. A total of 874 BAC contigs containing 5067 clones (22% of the library) were constructed by PCR-based screening with sequence-tagged sites (STSs) derived from whole-genome shotgun (WGS) sequences. A total of 523 BAC contigs, including 342 independent genes registered in public databases and 85 expressed sequence tags (ESTs), were placed onto the linkage map. We found significant synteny and conserved gene order between B. mori and a nymphalid butterfly, Heliconius melpomene, in four linkage groups (LGs), strongly suggesting that using B. mori as a reference for comparative genomics in Lepidotera is highly feasible.

Construction of a single nucleotide polymorphism linkage map for the silkworm, Bombyx mori, based on bacterial artificial chromosome end sequences

We have developed a linkage map for the silkworm Bombyx mori based on single nucleotide polymorphisms (SNPs) between strains p50T and C108T initially found on regions corresponding to the end sequences of bacterial artificial chromosome (BAC) clones. Using 190 segregants from a backcross of a p50T female x an F1 (p50T x C108T) male, we analyzed segregation patterns of 534 SNPs between p50T and C108T, detected among 3840 PCR amplicons, each associated with a p50T BAC end sequence. This enabled us to construct a linkage map composed of 534 SNP markers spanning 1305 cM in total length distributed over the expected 28 linkage groups. Of the 534 BACs whose ends harbored the SNPs used to construct the linkage map, 89 were associated with 107 different ESTs. Since each of the SNP markers is directly linked to a specific genomic BAC clone and to whole-genome sequence data, and some of them are also linked to EST data, the SNP linkage map will be a powerful tool for investigating silkworm genome properties, mutation mapping, and map-based cloning of genes of industrial and agricultural interest.

Genetic diversity among silkworm (Bombyx mori L., Lep., Bombycidae) germplasms revealed by microsatellites

To determine genetic relationships among strains of silkworm, Bombyx mori L., 31 strains with different origins, number of generations per year, number of molts per generation, and morphological characters were studied using simple sequence repeat (SSR) markers. Twenty-six primer pairs flanking microsatellite sequences in the silkworm genome were assayed. All were polymorphic and unambiguously separated silkworm strains from each other. A total of 188 alleles were detected with a mean value of 7.2 alleles/locus (range 2-17). The average heterozygosity value for each SSR locus ranged from 0 to 0.60, and the highest one was 0.96 (Fl0516 in 4013). The mean polymorphism index content (PIC) was 0.66 (range 0.12-0.89). Unweighted pair group method with arithmetic means (UPGMA) cluster analysis of Nei's genetic distance grouped silkworm strains based on their origin. Seven major ecotypic silkworm groups were analyzed. Principal components analysis (PCA) for SSR data support their UPGMA clustering. The results indicated that SSR markers are an efficient tool for fingerprinting cultivars and conducting genetic-diversity studies in the silkworm.

DNA fingerprinting using AFLP markers to search for markers associated with yield attributes in the silkworm, Bombyx mori

This study was carried out on 11 Chinese and 12 Japanese silkworm strains maintained by the Center for the Technological Development of Sericulture (CDTS) germplasm bank, located in Pereira, Colombia. The goals were to determine the genetic population structure of the two groups and the association between molecular markers (AFLPs) and important productivity characters. Group analysis showed the separation of the strains according to their geographic origin. The molecular markers and the productivity characters were correlated by multiple variance analysis. The analysis permitted the identification of molecular markers associated with the cocoon weight or the shell weight separately. Some markers were associated with both characters.

Integration of molecular and classical linkage groups of the silkworm, Bombyx mori (n = 28)

Previously published linkage groups (LGs) composed of molecular markers were assigned to classical LGs in the silkworm, Bombyx mori (n = 28). Four markers from the classical linkage map, og, w-1, Lp, and Pfl, were assigned to the molecular linkage maps using sequence tagged sites. In addition, linkage analysis was carried out using BF1 progeny between wild-type and mutant stocks carrying morphological phenotypic markers. As a result, the counterparts for 26 of 28 molecular LGs were identified with their counterparts of the classical LGs. Two visible markers, Sel and Xan, representing different classical LGs, were found to be linked.

Simple sequence repeat-based consensus linkage map of Bombyx mori

We established a genetic linkage map employing 518 simple sequence repeat (SSR, or microsatellite) markers for Bombyx mori (silkworm), the economically and culturally important lepidopteran insect, as part of an international genomics program. A survey of six representative silkworm strains using 2,500 (CA)n- and (CT)n-based SSR markers revealed 17-24% polymorphism, indicating a high degree of homozygosity resulting from a long history of inbreeding. Twenty-nine SSR linkage groups were established in well characterized Dazao and C108 strains based on genotyping of 189 backcross progeny derived from an F(1) male mated with a C108 female. The clustering was further focused to 28 groups by genotyping 22 backcross progeny derived from an F(1) female mated with a C108 male. This set of SSR linkage groups was further assigned to the 28 chromosomes (established linkage groups) of silkworm aided by visible mutations and cleaved amplified polymorphic sequence markers developed from previously mapped genes, cDNA sequences, and cloned random amplified polymorphic DNAs. By integrating a visible mutation p (plain, larval marking) and 29 well conserved genes of insects onto this SSR-based linkage map, a second generation consensus silkworm genetic map with a range of 7-40 markers per linkage group and a total map length of approximately 3431.9 cM was constructed and its high efficiency for genotyping and potential application for synteny studies of Lepidoptera and other insects was demonstrated.

Genetic mapping of Z chromosome and identification of W chromosome-specific markers in the silkworm, Bombyx mori

In the silkworm, Bombyx mori, the female is the heterogametic (ZW) sex and the male is homogametic (ZZ). The female heterogamety is a typical situation in the insect order Lepidoptera. Although the W chromosome in silkworm is strongly female determining, no W-linked gene for a morphological character has been found on it. The Z chromosome carries important traits of economic value as well as genes for various phenotypic traits, but only 2% of molecular information based on its relative size is known. Studies conducted so far indicate that the Z-linked genes are not dosage compensated. In the present study, we constructed a genetic map of randomly amplified polymorphic DNA fragments (RAPD), simple sequence repeats (SSR), and fluorescent intersimple sequence repeat PCR (FISSR) markers for the Z chromosome using a backcross mapping population. A total of 16 Z-linked markers were identified, characterized, and mapped using od, a recessive trait for translucent skin as an anchor marker yielding a total recombination map of 334.5 cM. The linkage distances obtained suggested that the markers were distributed throughout the Z chromosome. Four RAPD and four SSR markers that were linked to W chromosome were also identified. The proposed mapping approach should be useful to identify and map sex-linked traits in the silkworm. The economic and evolutionary significance of Z- and W-linked genes in silkworm, in particular, and lepidopterans, in general, is discussed.

Improvement of the cDNA-AFLP method using fluorescent primers for transcription analysis in bacteria

Here the cDNA-fluorescent amplified fragment length polymorphism (cDNA-FAFLP) technique, an improvement of cDNA-AFLP method, was used in order to analyse expression profiling in bacteria. The obtained results were validated by real-time PCR. This is the first report that validates the foreseen transcription pattern by cDNA-FAFLP with the application of real-time PCR in bacteria. This new protocol offers the possibility of quick and reliable analysis of transcription profiles, also avoids the problem linked with the use of radioisotopes, and allows a quick identification of genes differentially expressed in bacteria.

An AFLP-based interspecific linkage map of sympatric, hybridizing Colias butterflies

Colias eurytheme and C. philodice are sister species with broad sympatry in North America. They hybridize frequently and likely share a significant portion of their genomes through introgression. Both taxa have been ecologically well characterized and exploited to address a broad spectrum of evolutionary issues. Using AFLP markers, we constructed the first linkage map of Colias butterflies. The map is composed of 452 markers spanning 2541.7 cM distributed over 51 linkage groups (40 major groups and 11 small groups with 2-4 markers). Statistical tests indicate that these AFLP markers tend to cluster over the map, with the coefficient of variation of interval sizes being 1.236 (95% C.I. is 1.234-1.240). This nonrandom marker distribution can account for the nonequivalence between the number of linkage groups and the actual haploid chromosome number (N = 31). This study presents the initial step for further marker-assisted research on Colias butterflies, including QTL and introgression analyses. Further investigation of the genomes will help us understand better the roles of introgression and natural selection in the evolution of hybridizing species and devise more appropriate strategies to control these pests.

The Bombyx mori karyotype and the assignment of linkage groups

Lepidopteran species have a relatively high number of small holocentric chromosomes (Bombyx mori, 2n = 56). Chromosome identification has long been hampered in this group by the high number and by the absence of suitable markers like centromere position and chromosome bands. In this study, we carried out fluorescence in situ hybridization (FISH) on meiotic chromosome complements using genetically mapped B. mori bacterial artificial chromosomes (BACs) as probes. The combination of two to four either green or red fluorescence-labeled probes per chromosome allowed us to recognize unequivocally each of the 28 bivalents of the B. mori karyotype by its labeling pattern. Each chromosome was assigned one of the already established genetic linkage groups and the correct orientation in the chromosome was defined. This facilitates physical mapping of any other sequence and bears relevance for the ongoing B. mori genome projects. Two-color BAC-FISH karyotyping overcomes the problem of chromosome recognition in organisms where conventional banding techniques are not available.

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.