Survey and analysis of microsatellites in the silkworm, Bombyx mori: frequency, distribution, mutations, marker potential and their conservation in heterologous species

Genome-wide characterization of microsatellite DNA in fishes: survey and analysis of their abundance and frequency in genome-specific regions

BACKGROUND

Microsatellite repeats are ubiquitous in organism genomes and play an important role in the chromatin organization, regulation of gene activity, recombination and DNA replication. Although microsatellite distribution patterns have been studied in most phylogenetic lineages, they are unclear in fish species.

RESULTS

Here, we present the first systematic examination of microsatellite distribution in coding and non-coding regions of 14 fish genomes. Our study showed that the number and type of microsatellites displayed nonrandom distribution for both intragenic and intergenic regions, suggesting that they have potential roles in transcriptional or translational regulation and DNA replication slippage theories alone were insufficient to explain the distribution patterns. Our results showed that microsatellites are dominant in non-coding regions. The total number of microsatellites ranged from 78,378 to 1,012,084, and the relative density varied from 4925.76 bp/Mb to 25,401.97 bp/Mb. Overall, (A + T)-rich repeats were dominant. The dependence of repeat abundance on the length of the repeated unit (1-6 nt) showed a great similarity decrease, whereas more tri-nucleotide repeats were found in exonic regions than tetra-nucleotide repeats of most species. Moreover, the incidence of different repeated types appeared species- and genomic-specific. These results highlight potential mechanisms for maintaining microsatellite distribution, such as selective forces and mismatch repair systems.

CONCLUSIONS

Our data could be beneficial for the studies of genome evolution and microsatellite DNA evolutionary dynamics, and facilitate the exploration of microsatellites structural, function, composition mode and molecular markers development in these species.

Role of microsatellites in genetic analysis of Bombyx mori silkworm: a review

In the genome of Bombyx mori Linnaeus (1758), the microsatellites, or simple sequence repeats (SSR), feature among their particular characteristics a high adenine and thymine (A/T) content, low number of repeats, low frequency, and a grouping in "families" with similar flanking regions. Such characteristics may be the result of a complex interaction between factors that limit the size and dispersion of SSR loci—such as their high association with transposons—and mean that microsatellites within this taxon suitable as molecular markers are relatively rare. The determination of genetic profiles in populations and cell lines has not been affected owing to the high level of polymorphism, nor has the analysis of diversity, structure and genetic relationships. However, the scarcity of suitable microsatellites has restricted their application in genetic mapping, limiting them to preliminary identification of gene location of genes or quantitative trait loci (QTLs) related to thermotolerance, resistance to viruses, pigmentation patterns, body development and the weight of the cocoon, the cortex, the pupa and the filament. The review confirms that, as markers, microsatellites are versatile and perform well. They could thus be useful both to advance research in emerging countries with few resources seeking to promote sericulture in their territories, and to advance in the genetic and molecular knowledge of characteristics of productive and biological interest, given the latest technological developments in terms of the sequencing, identification, isolation and genotyping of SSR loci.

Genome-wide mining of perfect microsatellites and tetranucleotide orthologous microsatellites estimates in six primate species

Advancement in genome sequencing and in silico mining tools have provided new opportunities for comparative primate genomics of microsatellites. The SSRs (simple sequence repeats) numbers were not correlated with the genome size (Pearson, r=0.310, p=0.550), and were positively correlated with the total length of SSRs (Pearson, r=0.992, p=0.00). A total of 224,289 tetranucleotide orthologous microsatellites families and 367 single-copy orthologous SSRs loci were found in six primate species by homologous alignment. The inner mutation types of single-copy orthologous SSRs loci included the copy number variance, point mutation, and chromosomal translocation. The accumulated repeat times and average length of tetranucleotide orthologous microsatellites in Rhinopithecus roxellana, Papio anubis and Macaca mulatta were longer than Homo sapiens and Pan troglodytes, which showed the tetranucleotide orthologous SSRs loci had more repeat times and longer average length on the branches with earlier divergence time, one exception may be Microcebus murinus as a primitive monkey with a smallest morphology in Malagasy. Our conclusion indicated that single-copy tetranucleotide orthologous SSRs sequences accumulated individual mutation more slowly through time in H. sapiens and P. troglodytes than in R. roxellanae, P. anubis and M. mulatta. However, such divergence wouldn't arise uniformly in all branches of the primate tree. A comparison of genomic sequence assemblages would offer remarkable insights about comparisons and contrasts, and the evolutionary processes of the microsatellites involved in human and nonhuman primate species.

Molecular Markers for Biomass Traits: Association, Interaction and Genetic Divergence in Silkworm Bombyx mori

Improvement of high yielding, disease resistant silkworm strains became imminent to increase production of silk, which is a major revenue earner for sericulturists. Since environment interacts with phenotype, conventional breeding did not result in commendable yield improvement in synthetic strains of silkworm, Bombyx mori. Identification of DNA markers associated with different economically important biomass traits and its introgression could assist molecular breeding and expression of stabilized high yielding characters, but genetic basis of most quantitative traits in silkworm is poorly understood due to its polygenic control. Correlation analysis (R = 0.9) revealed significant interrelation among biomass traits viz., larval duration (TLD), larval weight (LWT), cocoon weight (CWT), shell weight (SWT), shell ratio (SR) and floss content. PCR using inter simple sequence repeat (ISSR) primers revealed 92% polymorphism among 14 tropical and temperate strains of B. mori, with average diversity index of 0.747. Stepwise multiple regression analysis (MRA) selected 35 ISSR markers positively or negatively correlated with different biomass traits, illustrated polygenic control. ISSR marker 830.81050bp was significantly associated with LWT, CWT, SWT, SR and floss content, indicated its pleiotropic role. Two ISSR markers, 835.51950bp and 825.9710bp showed significant association with floss content and TLD. These markers were segregated in F2 generation and Chi-square test confirmed (χ2 = ~45; P < 0.05) its genetic contribution to the associated biomass traits. Strains, with both positively and negatively correlated markers, had intermediate mean value for biomass traits (eg. SWT = 0.17 ± 0.014 g in GNM and Moria) indicated interaction of loci in natural populations. Low yielding Indian strains grouped together by Hierarchical clustering. Chinese and Japanese strains were distributed in the periphery of ALSCAL matrix indicated convergence of genetic characters in Indian strains. Average genetic distance between Chinese strains and Indian strains (0.193) significantly ( P < 0.01) varied from that between Chinese and Japanese strains. Interaction of loci and allelic substitutions induced phenotypic plasticity in temperate B. mori populations on tropic adaptation in India. These outcomes show possibility to combine favorable alleles at different QTL to increase larval, cocoon and shell weight.

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.

Genome-Wide Survey and Analysis of Microsatellite Sequences in Bovid Species

Microsatellites or simple sequence repeats (SSRs) have become the most popular source of genetic markers, which are ubiquitously distributed in many eukaryotic and prokaryotic genomes. This is the first study examining and comparing SSRs in completely sequenced genomes of the Bovidae. We analyzed and compared the number of SSRs, relative abundance, relative density, guanine-cytosine (GC) content and proportion of SSRs in six taxonomically different bovid species: Bos taurus, Bubalus bubalis, Bos mutus, Ovis aries, Capra hircus, and Pantholops hodgsonii. Our analysis revealed that, based on our search criteria, the total number of perfect SSRs found ranged from 663,079 to 806,907 and covered from 0.44% to 0.48% of the bovid genomes. Relative abundance and density of SSRs in these Bovinae genomes were non-significantly correlated with genome size (Pearson, r < 0.420, p > 0.05). Perfect mononucleotide SSRs were the most abundant, followed by the pattern: perfect di- > tri- > penta- > tetra- > hexanucleotide SSRs. Generally, the number of SSRs, relative abundance, and relative density of SSRs decreased as the motif repeat length increased in each species of Bovidae. The most GC-content was in trinucleotide SSRs and the least was in the mononucleotide SSRs in the six bovid genomes. The GC-contents of tri- and pentanucleotide SSRs showed a great deal of similarity among different chromosomes of B. taurus, O. aries, and C. hircus. SSR number of all chromosomes in the B. taurus, O.aries, and C. hircus is closely positively correlated with chromosome sequence size (Pearson, r > 0.980, p < 0.01) and significantly negatively correlated with GC-content (Pearson, r < -0.638, p < 0.01). Relative abundance and density of SSRs in all chromosomes of the three species were significantly negatively correlated with GC-content (Pearson, r < -0.333, P < 0.05) but not significantly correlated with chromosome sequence size (Pearson, r < -0.185, P > 0.05). Relative abundances of the same nucleotide SSR type showed great similarity among different chromosomes of B. taurus, O. aries, and C. hircus.

Loss of genetic diversity and inbreeding in Kashmir red deer (Cervus elaphus hanglu) of Dachigam National Park, Jammu & Kashmir, India

Background

Hangul (Cervus elaphus hanglu), the eastern most subspecies of red deer, is now confined only to the mountains in the Kashmir region of Jammu & Kashmir State of India. It is of great conservation significance as this is the last and only hope for Asiatic survivor of the red deer species in India. Wild population of free ranging hangul deer inhabiting in and around Dachigam National Park was genetically assessed in order to account for constitutive genetic attributes of hangul population using microsatellite markers.

Results

In a pool of 36 multi-locus genotypes, 30 unique individuals were identified based on six microsatellite loci. The estimated cumulative probability of identity assuming all individuals were siblings (P_ID sibs) was 0.009 (9 in 1000). Altogether, 49 different alleles were observed with mean (± s.e.) allelic number of 8.17 ± 1.05, ranging from 5 to 11 per locus. The observed heterozygosity ranged between 0.08 and 0.83, with mean 0.40 ± 0.11 and the inbreeding coefficient ranged between −0.04 and 0.87 with mean 0.38 ± 0.15. Majority of loci (5/6) were found to be informative (PIC value > 0.5). All loci deviated from Hardy-Weinberg equilibrium except Ca-38 (P > 0.05) and none of the pairs of loci showed significant linkage disequilibrium except the single pair of Ca-30 and Ca-43 (P < 0.05).

Conclusions

The preliminary findings revealed that hangul population is significantly inbred and exhibited a low genetic diversity in comparison to other deer populations of the world. We suggest prioritizing the potential individuals retaining high heterozygosity for ex situ conservation and genetic monitoring of the hangul population should be initiated covering the entire distribution range to ensure the long term survival of hangul. We speculate further ignoring genetics attributes may lead to a detrimental effect which can negatively influence the reproductive fitness and survivorship of the hangul population in the wild.

Microsatellites for the marsh fritillary butterfly: de novo transcriptome sequencing, and a comparison with amplified fragment length polymorphism (AFLP) markers

BACKGROUND

Until recently the isolation of microsatellite markers from Lepidoptera has proved troublesome, expensive and time-consuming. Following on from a previous study of Edith's checkerspot butterfly, Euphydryas editha, we developed novel microsatellite markers for the vulnerable marsh fritillary butterfly, E. aurinia. Our goal was to optimize the process in order to reduce both time and cost relative to prevailing techniques. This was accomplished by using a combination of previously developed techniques: in silico mining of a de novo assembled transcriptome sequence, and genotyping the microsatellites found there using an economic method of fluorescently labelling primers.

PRINCIPAL FINDINGS

In total, we screened nine polymorphic microsatellite markers, two of which were previously published, and seven that were isolated de novo. These markers were able to amplify across geographically isolated populations throughout Continental Europe and the UK. Significant deviations from Hardy-Weinberg equilibrium were evident in some populations, most likely due to the presence of null alleles. However, we used an F(st) outlier approach to show that these markers are likely selectively neutral. Furthermore, using a set of 128 individuals from 11 populations, we demonstrate consistency in population differentiation estimates with previously developed amplified fragment length polymorphism (AFLP) markers (r = 0.68, p<0.001).

SIGNIFICANCE

Rapid development of microsatellite markers for difficult taxa such as Lepidoptera, and concordant results with other putatively neutral molecular markers, demonstrate the potential of de novo transcriptional sequencing for future studies of population structure and gene flow that are desperately needed for declining species across fragmented landscapes.

Development and application of microsatellites in candidate genes related to wood properties in the Chinese white poplar (Populus tomentosa Carr.)

Gene-derived simple sequence repeats (genic SSRs), also known as functional markers, are often preferred over random genomic markers because they represent variation in gene coding and/or regulatory regions. We characterized 544 genic SSR loci derived from 138 candidate genes involved in wood formation, distributed throughout the genome of Populus tomentosa, a key ecological and cultivated wood production species. Of these SSRs, three-quarters were located in the promoter or intron regions, and dinucleotide (59.7%) and trinucleotide repeat motifs (26.5%) predominated. By screening 15 wild P. tomentosa ecotypes, we identified 188 polymorphic genic SSRs with 861 alleles, 2–7 alleles for each marker. Transferability analysis of 30 random genic SSRs, testing whether these SSRs work in 26 genotypes of five genus Populus sections (outgroup, Salix matsudana), showed that 72% of the SSRs could be amplified in Turanga and 100% could be amplified in Leuce. Based on genotyping of these 26 genotypes, a neighbour-joining analysis showed the expected six phylogenetic groupings. In silico analysis of SSR variation in 220 sequences that are homologous between P. tomentosa and Populus trichocarpa suggested that genic SSR variations between relatives were predominantly affected by repeat motif variations or flanking sequence mutations. Inheritance tests and single-marker associations demonstrated the power of genic SSRs in family-based linkage mapping and candidate gene-based association studies, as well as marker-assisted selection and comparative genomic studies of P. tomentosa and related species.

Activation of autophagic programmed cell death and innate immune gene expression reveals immuno-competence of integumental epithelium in Bombyx mori infected by a dipteran parasitoid

In insects, the integument forms the primary barrier between the environment and internal milieu, but cellular and immune responses of the integumental epithelium to infection by micro- and macro-parasites are mostly unknown. We elucidated cellular and immune responses of the epithelium induced through infection by a dipteran endoparasitoid, Exorista bombycis in the economically important silkworm Bombyx mori. Degradative autophagic vacuoles, lamella-like bodies, a network of cytoplasmic channels with cellular cargo, and an RER network that opened to vacuoles were observed sequentially with increase in age after infection. This temporal sequence culminated in apoptosis, accompanied by the upregulation of the caspase gene and fragmentation of DNA. The infection significantly enhanced the tyrosine level and phenol oxidase activity in the integument. Proteomic analysis revealed enhanced expression of innate immunity components of toll and melanization pathways, cytokines, signaling molecules, chaperones, and proteolytic enzymes demonstrating diverse host responses. qPCR analysis revealed the upregulation of spatzle, BmToll, and NF kappa B transcription factors Dorsal and BmRel. NF kappa B inhibitor cactus showed diminished expression when Dorsal and BmRel were upregulated, revealing a negative correlation (R = (-)0.612). During melanization, prophenol oxidase 2 was expressed, a novel finding in integumental epithelium. The integument showed a low level of melanin metabolism and localized melanism in order to prevent the spreading of cytotoxic quinones. The gene-encoding proteolytic enzyme, beta-N-acetylglucosaminidase, was activated at 24 h post-infection, whereas chitinase, was activated at 96 h post-infection; however, most of the immune genes enhanced their expression in the early stages of infection. Thus the integument contributes to humoral immune responses that enhance resistance against macroparasite invasion.

The spatial genetic differentiation of the legume pod borer, Maruca vitrata F. (Lepidoptera: Crambidae) populations in West Africa

The legume pod borer, Maruca vitrata, is an endemic insect pest that causes significant yield loss to the cowpea crop in West Africa. The application of population genetic tools is important in the management of insect pests but such data on M. vitrata is lacking. We applied a set of six microsatellite markers to assess the population structure of M. vitrata collected at five sites from Burkina Faso, Niger and Nigeria. Observed polymorphisms ranged from one (marker 3393) to eight (marker 32008) alleles per locus. Observed and expected heterozygosities ranged from 0.0 to 0.8 and 0.0 to 0.6, respectively. Three of the loci in samples from Nigeria and Burkina Faso deviated significantly from Hardy-Weinberg Equilibrium (HWE), whereas no loci deviated significantly in samples from Niger. Analysis of molecular variance (AMOVA) indicated that 67.3% level of the genetic variation was within individuals compared to 17.3% among populations. A global estimate of F ST=0.1 (ENA corrected F ST=0.1) was significant (P⩽0.05) and corroborated by pairwise F ST values that were significant among all possible comparisons. A significant correlation was predicted between genetic divergence and geographic distance between subpopulations (R2=0.6, P=0.04), and cluster analysis by the program STRUCTURE predicted that co-ancestry of genotypes were indicative of three distinct populations. The spatial genetic variance among M. vitrata in West Africa may be due to limited gene flow, south-north seasonal movement pattern or other reproductive barriers. This information is important for the cultural, chemical and biological control strategies for managing M. vitrata.

Application of high-resolution DNA melting for genotyping in lepidopteran non-model species: Ostrinia furnacalis (Crambidae)

Development of an ideal marker system facilitates a better understanding of the genetic diversity in lepidopteran non-model organisms, which have abundant species, but relatively limited genomic resources. Single nucleotide polymorphisms (SNPs) discovered within single-copy genes have proved to be desired markers, but SNP genotyping by current techniques remain laborious and expensive. High resolution melting (HRM) curve analysis represents a simple, rapid and inexpensive genotyping method that is primarily confined to clinical and diagnostic studies. In this study, we evaluated the potential of HRM analysis for SNP genotyping in the lepidopteran non-model species Ostrinia furnacalis (Crambidae). Small amplicon and unlabeled probe assays were developed for the SNPs, which were identified in 30 females of O. furnacalis from 3 different populations by our direct sequencing. Both assays were then applied to genotype 90 unknown female DNA by prior mixing with known wild-type DNA. The genotyping results were compared with those that were obtained using bi-directional sequencing analysis. Our results demonstrated the efficiency and reliability of the HRM assays. HRM has the potential to provide simple, cost-effective genotyping assays and facilitates genotyping studies in any non-model lepidopteran species of interest.

Fosmid library end sequencing reveals a rarely known genome structure of marine shrimp Penaeus monodon

Background

The black tiger shrimp (Penaeus monodon) is one of the most important aquaculture species in the world, representing the crustacean lineage which possesses the greatest species diversity among marine invertebrates. Yet, we barely know anything about their genomic structure. To understand the organization and evolution of the P. monodon genome, a fosmid library consisting of 288,000 colonies and was constructed, equivalent to 5.3-fold coverage of the 2.17 Gb genome. Approximately 11.1 Mb of fosmid end sequences (FESs) from 20,926 non-redundant reads representing 0.45% of the P. monodon genome were obtained for repetitive and protein-coding sequence analyses.

Results

We found that microsatellite sequences were highly abundant in the P. monodon genome, comprising 8.3% of the total length. The density and the average length of microsatellites were evidently higher in comparison to those of other taxa. AT-rich microsatellite motifs, especially poly (AT) and poly (AAT), were the most abundant. High abundance of microsatellite sequences were also found in the transcribed regions. Furthermore, via self-BlastN analysis we identified 103 novel repetitive element families which were categorized into four groups, i.e., 33 WSSV-like repeats, 14 retrotransposons, 5 gene-like repeats, and 51 unannotated repeats. Overall, various types of repeats comprise 51.18% of the P. monodon genome in length. Approximately 7.4% of the FESs contained protein-coding sequences, and the Inhibitor of Apoptosis Protein (IAP) gene and the Innexin 3 gene homologues appear to be present in high abundance in the P. monodon genome.

Conclusions

The redundancy of various repeat types in the P. monodon genome illustrates its highly repetitive nature. In particular, long and dense microsatellite sequences as well as abundant WSSV-like sequences highlight the uniqueness of genome organization of penaeid shrimp from those of other taxa. These results provide substantial improvement to our current knowledge not only for shrimp but also for marine crustaceans of large genome size.

A review of the implications of heterozygosity and inbreeding on germplasm biodiversity and its conservation in the silkworm, Bombyx mori

Abstract Silkworm genebanks assume paramount importance as the reservoirs of biodiversity and source of alleles that can be easily retrieved for genetic enhancement of popular breeds. More than 4000 Bombyx mori L (Lepidoptera: Bombycidae) strains are currently available and these strains are maintained through continuous sibling mating. This repeated sibling mating makes the populations of each strain more homozygous, but leads to loss of unique and valuable genes through the process of inbreeding depression. Hence, it is essential to maintain a minimal degree of heterozygosity within the population of each silkworm strain, especially in the traditional geographic strains, to avoid such loss. As a result, accurate estimation of genetic diversity is becoming more important in silkworm genetic resources conservation. Application of molecular markers help estimate genetic diversity much more accurately than that of morphological traits. Since a minimal amount of heterozygosity in each silkworm strain is essential for better conservation by avoiding inbreeding depression, this article overviews both theoretical and practical importance of heterozygosity together with impacts of inbreeding depression and the merits and demerits of neutral molecular markers for measurements of both heterozygosity and inbreeding depression in the silkworm Bombyx mori.

Mining of simple sequence repeats in the Genome of Gentianaceae

Simple sequence repeats (SSRs) or short tandem repeats are short repeat motifs that show high level of length polymorphism due to insertion or deletion mutations of one or more repeat types. Here, we present the detection and abundance of microsatellites or SSRs in nucleotide sequences of Gentianaceae family. A total of 545 SSRs were mined in 4698 nucleotide sequences downloaded from the National Center for Biotechnology Information (NCBI). Among the SSR sequences, the frequency of repeat type was about 429 -mono repeats, 99 -di repeats, 15 -tri repeats, and 2 --hexa repeats. Mononucleotide repeats were found to be abundant repeat types, about 78%, followed by dinucleotide repeats (18.16%) among the SSR sequences. An attempt was made to design primer pairs for 545 identified SSRs but these were found only for 169 sequences.

FONZIE: An optimized pipeline for minisatellite marker discovery and primer design from large sequence data sets

BACKGROUND

Micro-and minisatellites are among the most powerful genetic markers known to date. They have been used as tools for a large number of applications ranging from gene mapping to phylogenetic studies and isolate typing. However, identifying micro-and minisatellite markers on large sequence data sets is often a laborious process.

RESULTS

FONZIE was designed to successively 1) perform a search for markers via the external software Tandem Repeat Finder, 2) exclude user-defined specific genomic regions, 3) screen for the size and the percent matches of each relevant marker found by Tandem Repeat Finder, 4) evaluate marker specificity (i.e., occurrence of the marker as a single copy in the genome) using BLAST2.0, 5) design minisatellite primer pairs via the external software Primer3, and 6) check the specificity of each final PCR product by BLAST. A final file returns to users all the results required to amplify markers. A biological validation of the approach was performed using the whole genome sequence of the phytopathogenic fungus Leptosphaeria maculans, showing that more than 90% of the minisatellite primer pairs generated by the pipeline amplified a PCR product, 44.8% of which showed agarose-gel resolvable polymorphism between isolates. Segregation analyses confirmed that the polymorphic minisatellites corresponded to single-locus markers.

CONCLUSION

FONZIE is a stand-alone and user-friendly application developed to minimize tedious manual operations, reduce errors, and speed up the search for efficient minisatellite and microsatellite markers departing from whole-genome sequence data. This pipeline facilitates the integration of data and provides a set of specific primer sequences for PCR amplification of single-locus markers. FONZIE is freely downloadable at: http://www.versailles-grignon.inra.fr/bioger/equipes/leptosphaeria_maculans/outils_d_analyses/fonzie.

Comparative genomic analysis of simple sequence repeats in three Plasmodium species

Simple sequence repeats (SSRs) are known to be responsible for genetic complexities and play major roles in gene and genome evolution. To this respect, malaria parasites are known to have rapidly evolving and complex genomes with complicated and differential pathogenic behaviors. Hence, by studying the whole genome comparative SSRs patterns, one can understand genomic complexities and differential evolutionary patterns of these species. We herein utilized the whole genome sequence information of three Plasmodium species, Plasmodium falciparum, Plasmodium vivax, and Plasmodium knowlesi, to comparatively analyze genome-wide distribution of SSRs. The study revealed that despite having the smallest genome size, P. falciparum bears the highest SSR content among the three Plasmodium species. Furthermore, distribution patterns of different SSRs types (e.g., mono, di, tri, tetra, penta, and hexa) in term of relative abundance and relative density provide evidences for greater accumulation of di-repeats and marked decrease of mono-repeats in P. falciparum in comparison to other two species. Overall, the types and distribution of SSRs in P. falciparum genome was found to be different than that of P. vivax and P. knowlesi. The latter two species have quite similar SSR organizations in many aspects of the data. The results were discussed in terms of comparative SSR patterns among the three Plasmodium species, uniqueness of P. falciparum in SSR organization and general pattern of evolution of SSRs in Plasmodium.

What is a microsatellite: a computational and experimental definition based upon repeat mutational behavior at A/T and GT/AC repeats

Microsatellites are abundant in eukaryotic genomes and have high rates of strand slippage-induced repeat number alterations. They are popular genetic markers, and their mutations are associated with numerous neurological diseases. However, the minimal number of repeats required to constitute a microsatellite has been debated, and a definition of a microsatellite that considers its mutational behavior has been lacking. To define a microsatellite, we investigated slippage dynamics for a range of repeat sizes, utilizing two approaches. Computationally, we assessed length polymorphism at repeat loci in ten ENCODE regions resequenced in four human populations, assuming that the occurrence of polymorphism reflects strand slippage rates. Experimentally, we determined the in vitro DNA polymerase-mediated strand slippage error rates as a function of repeat number. In both approaches, we compared strand slippage rates at tandem repeats with the background slippage rates. We observed two distinct modes of mutational behavior. At small repeat numbers, slippage rates were low and indistinguishable from background measurements. A marked transition in mutability was observed as the repeat array lengthened, such that slippage rates at large repeat numbers were significantly higher than the background rates. For both mononucleotide and dinucleotide microsatellites studied, the transition length corresponded to a similar number of nucleotides (approximately 10). Thus, microsatellite threshold is determined not by the presence/absence of strand slippage at repeats but by an abrupt alteration in slippage rates relative to background. These findings have implications for understanding microsatellite mutagenesis, standardization of genome-wide microsatellite analyses, and predicting polymorphism levels of individual microsatellite loci.

Mining non-model genomic libraries for microsatellites: BAC versus EST libraries and the generation of allelic richness

BACKGROUND

Simple sequence repeats (SSRs) are tandemly repeated sequence motifs common in genomic nucleotide sequence that often harbor significant variation in repeat number. Frequently used as molecular markers, SSRs are increasingly identified via in silico approaches. Two common classes of genomic resources that can be mined are bacterial artificial chromosome (BAC) libraries and expressed sequence tag (EST) libraries.

RESULTS

288 SSR loci were screened in the rapidly radiating Hawaiian swordtail cricket genus Laupala. SSRs were more densely distributed and contained longer repeat structures in BAC library-derived sequence than in EST library-derived sequence, although neither repeat density nor length was exceptionally elevated despite the relatively large genome size of Laupala. A non-random distribution favoring AT-rich SSRs was observed. Allelic diversity of SSRs was positively correlated with repeat length and was generally higher in AT-rich repeat motifs.

CONCLUSION

The first large-scale survey of Orthopteran SSR allelic diversity is presented. Selection contributes more strongly to the size and density distributions of SSR loci derived from EST library sequence than from BAC library sequence, although all SSRs likely are subject to similar physical and structural constraints, such as slippage of DNA replication machinery, that may generate increased allelic diversity in AT-rich sequence motifs. Although in silico approaches work well for SSR locus identification in both EST and BAC libraries, BAC library sequence and AT-rich repeat motifs are generally superior SSR development resources for most applications.

A helitron-like transposon superfamily from lepidoptera disrupts (GAAA)(n) microsatellites and is responsible for flanking sequence similarity within a microsatellite family

Transposable elements (TEs) are mobile DNA regions that alter host genome structure and gene expression. A novel 588 bp non-autonomous high copy number TE in the Ostrinia nubilalis genome has features in common with miniature inverted-repeat transposable elements (MITEs): high A + T content (62.3%), lack of internal protein coding sequence, and secondary structure consisting of subterminal inverted repeats (SIRs). The O. nubilalis TE has inserted at (GAAA)(n) microsatellite loci, and was named the microsatellite-associated interspersed nuclear element (MINE-1). Non-autonomous MINE-1 superfamily members also were identified downstream of (GAAA)(n) microsatellites within Bombyx mori and Pectinophora gossypiella genomes. Of 316 (GAAA)(n) microsatellites from the B. mori whole genome sequence, 201 (63.6%) have associated autonomous or non-autonomous MINE-1 elements. Autonomous B. mori MINE-1s a encode a helicase and endonuclease domain RepHel-like protein (BMHELp1) indicating their classification as Helitron-like transposons and were renamed Helitron1_BM. Transposition of MINE-1 members in Lepidoptera has resulted in the disruption of (GAAA)(n) microsatellite loci, has impacted the application of microsatellite-based genetic markers, and suggests genome sequence that flanks TT/AA dinucleotides may be required for target site recognition by RepHel endonuclease domains.

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.

Microsatellite markers associated with genes expressed in developing wings of Bicyclus anynana butterflies

Deriving useful microsatellite markers in lepidopterans has been challenging when relying on scans of genomic DNA libraries, presumably due to repetitiveness in their genomes. We assayed 96 of 320 microsatellites identified in silico from a collection of Bicyclus anynana ESTs, in 11 independent individuals from a laboratory population. From the 68 successful assays, we identified 40 polymorphic markers including 22 with BLAST-based annotation. Nine of 12 selected polymorphic markers tested in a panel of 24 wild-caught individuals converted to successful assays and were all polymorphic. We discuss how microsatellite discovery in ESTs is an efficient strategy with important attendant advantages.

Exon-primed intron-crossing (EPIC) PCR markers of Helicoverpa armigera (Lepidoptera: Noctuidae)

Applying microsatellite DNA markers in population genetic studies of the pest moth Helicoverpa armigera is subject to numerous technical problems, such as the high frequency of null alleles, occurrence of size homoplasy, presence of multiple copies of flanking sequence in the genome and the lack of PCR amplification robustness between populations. To overcome these difficulties, we developed exon-primed intron-crossing (EPIC) nuclear DNA markers for H. armigera based on ribosomal protein (Rp) and the Dopa Decarboxylase (DDC) genes and sequenced alleles showing length polymorphisms. Allele length polymorphisms were usually from random indels (insertions or deletions) within introns, although variation of short dinucleotide DNA repeat units was also detected. Mapping crosses demonstrated Mendelian inheritance patterns for these EPIC markers and the absence of both null alleles and allele 'dropouts'. Three examples of allele size homoplasies due to indels were detected in EPIC markers RpL3, RpS6 and DDC, while sequencing of multiple individuals across 11 randomly selected alleles did not detect indel size homoplasies. The robustness of the EPIC-PCR markers was demonstrated by PCR amplification in the related species, H. zea, H. assulta and H. punctigera.

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.

Distribution of polymorphic and non-polymorphic microsatellite repeats in Xenopus tropicalis

The results of our bioinformatics analysis have found over 91,000 di-, tri-, and tetranucleotide microsatellites in our survey of 25% of the X. tropicalis genome, suggesting there may be over 360,000 within the entire genome. Within the X. tropicalis genome, dinucleotide (78.7%) microsatellites vastly out numbered tri- and tetranucleotide microsatellites. Similarly, AT-rich repeats are overwhelmingly dominant. The four AT-only motifs (AT, AAT, AAAT, and AATT) account for 51,858 out of 91,304 microsatellites found. Individually, AT microsatellites were the most common repeat found, representing over half of all di-, tri-, and tetranucleotide microsatellites. This contrasts with data from other studies, which show that AC is the most frequent microsatellite in vertebrate genomes (Toth et al. 2000). In addition, we have determined the rate of polymorphism for 5,128 non-redundant microsatellites, embedded in unique sequences. Interestingly, this subgroup of microsatellites was determined to have significantly longer repeats than genomic microsatellites as a whole. In addition, microsatellite loci with tandem repeat lengths more than 30 bp exhibited a significantly higher degree of polymorphism than other loci. Pairwise comparisons show that tetranucleotide microsatellites have the highest polymorphic rates. In addition, AAT and ATC showed significant higher polymorphism than other trinucleotide microsatellites, while AGAT and AAAG were significantly more polymorphic than other tetranucleotide microsatellites.

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.

Simple sequence repeats in Neurospora crassa: distribution, polymorphism and evolutionary inference

Background

Simple sequence repeats (SSRs) have been successfully used for various genetic and evolutionary studies in eukaryotic systems. The eukaryotic model organism Neurospora crassa is an excellent system to study evolution and biological function of SSRs.

Results

We identified and characterized 2749 SSRs of 963 SSR types in the genome of N. crassa. The distribution of tri-nucleotide (nt) SSRs, the most common SSRs in N. crassa, was significantly biased in exons. We further characterized the distribution of 19 abundant SSR types (AST), which account for 71% of total SSRs in the N. crassa genome, using a Poisson log-linear model. We also characterized the size variation of SSRs among natural accessions using Polymorphic Index Content (PIC) and ANOVA analyses and found that there are genome-wide, chromosome-dependent and local-specific variations. Using polymorphic SSRs, we have built linkage maps from three line-cross populations.

Conclusion

Taking our computational, statistical and experimental data together, we conclude that 1) the distributions of the SSRs in the sequenced N. crassa genome differ systematically between chromosomes as well as between SSR types, 2) the size variation of tri-nt SSRs in exons might be an important mechanism in generating functional variation of proteins in N. crassa, 3) there are different levels of evolutionary forces in variation of amino acid repeats, and 4) SSRs are stable molecular markers for genetic studies in N. crassa.

Expressed sequence tags from the midgut of Epiphyas postvittana (Walker) (Lepidoptera: Tortricidae)

The midgut is a key tissue in insect science. Physiological roles include digestion and peritrophic membrane function, as well as being an important target for insecticides. We used an expressed sequence tag (EST) approach to identify candidate genes and gene families involved in these processes in the light brown apple moth, Epiphyas postvittana (Walker) (Lepidoptera: Tortricidae). Two cDNA libraries were constructed from dissected midgut of third to fifth instar larvae. Clustering analysis of 6416 expressed sequence tags produced 1178 tentative unique genes comprising 725 tentative contigs and 453 singletons. The sequences show similar codon usage to sequences from other lepidopterans, a Kozak consensus sequence similar to Drosophila and single nucleotide polymorphisms (SNPs) were detected at a frequency of 1.35/kb. The identity of the most common Interpro families correlates well with major known functions of the midgut. Phylogenetic analysis was conducted on representative sequences from selected multigene families. Gene families include a broad range of digestive proteases, lipases and carbohydrases that appear to have degradative capacity against the major food components found in leaves, the diet of these larvae; and carboxylesterases, glutathione-S-transferases and cytochrome P450 monooxygenases, potentially involved in xenobiotic degradation. Two of the larger multigene families, serine proteases and lipases, expressed a high proportion of genes that are likely to be catalytically inactive.

Initial analysis of tandemly repetitive sequences in the genome of Zhikong scallop (Chlamys farreri Jones et Preston)

Tandemly repetitive sequences are widespread in all eukaryotic genomes, but data on tandem repeats are limited in Zhikong scallop (Chlamys farreri). In the present study, paired-end sequencing of 2016 individual fosmid clones resulted in 3646 sequences. A total of 2,286,986 bp of genomic sequences were generated, representing approximately 1.84 per thousand of the Zhikong scallop genome. Using tandem repeats finder (TRF) software, a total of 2500 tandem repeats were found, including 313 satellites, 1816 minisatellites and 371 microsatellites. The cumulative length of tandem repeats was 552,558 bp, accounting for 24.16% of total length. Specifically, the length of microsatellites, minisatellites and satellites was 9425, 336,001 and 207,132 bp, accounting for 1.71, 60.81 and 37.49% of the length of tandem repeats, and 0.41, 14.69 and 9.06% of total length, respectively. The detailed information on the characteristic of all repeat units was also represented, which will provide a useful resource for physical mapping and better utilization of the existing genomic information in Zhikong scallop.

Molecular Markers for Biomass Traits: Association, Interaction and Genetic Divergence in Silkworm Bombyx mori

Improvement of high yielding, disease resistant silkworm strains became imminent to increase production of silk, which is a major revenue earner for sericulturists. Since environment interacts with phenotype, conventional breeding did not result in commendable yield improvement in synthetic strains of silkworm, Bombyx mori. Identification of DNA markers associated with different economically important biomass traits and its introgression could assist molecular breeding and expression of stabilized high yielding characters, but genetic basis of most quantitative traits in silkworm is poorly understood due to its polygenic control. Correlation analysis (R = 0.9) revealed significant interrelation among biomass traits viz., larval duration (TLD), larval weight (LWT), cocoon weight (CWT), shell weight (SWT), shell ratio (SR) and floss content. PCR using inter simple sequence repeat (ISSR) primers revealed 92% polymorphism among 14 tropical and temperate strains of B. mori, with average diversity index of 0.747. Stepwise multiple regression analysis (MRA) selected 35 ISSR markers positively or negatively correlated with different biomass traits, illustrated polygenic control. ISSR marker 830.8(1050bp) was significantly associated with LWT, CWT, SWT, SR and floss content, indicated its pleiotropic role. Two ISSR markers, 835.5(1950bp) and 825.9(710bp) showed significant association with floss content and TLD. These markers were segregated in F(2) generation and Chi-square test confirmed (chi(2) = ~45; P < 0.05) its genetic contribution to the associated biomass traits. Strains, with both positively and negatively correlated markers, had intermediate mean value for biomass traits (eg. SWT = 0.17 +/- 0.014 g in GNM and Moria) indicated interaction of loci in natural populations. Low yielding Indian strains grouped together by Hierarchical clustering. Chinese and Japanese strains were distributed in the periphery of ALSCAL matrix indicated convergence of genetic characters in Indian strains. Average genetic distance between Chinese strains and Indian strains (0.193) significantly (P < 0.01) varied from that between Chinese and Japanese strains. Interaction of loci and allelic substitutions induced phenotypic plasticity in temperate B. mori populations on tropic adaptation in India. These outcomes show possibility to combine favorable alleles at different QTL to increase larval, cocoon and shell weight.

Microsatellite flanking region similarities among different loci within insect species

Although microsatellites are ubiquitous in eukaryota, the number of available markers varies strongly among taxa. This meta-analysis was conducted on 32 insect species. Sequences were obtained from two assembled whole genomes, whole genome shotgun (WGS) sequences from 10 species and screening partial genomic libraries for microsatellites from 23 species. We have demonstrated: (1) strong differences in the abundance of microsatellites among species; (2) that microsatellites within species are often grouped into families based on similarities in their flanking sequences; (3) that the proportion of microsatellites grouped into families varies strongly among taxa; and (4) that microsatellite families were significantly more often associated with transposable elements - or their remnants - than unique microsatellite sequences.

Evolutionary dynamics of multilocus microsatellite arrangements in the genome of the butterfly Bicyclus anynana, with implications for other Lepidoptera

The sequences flanking microsatellites isolated from the butterfly Bicyclus anynana display high levels of similarity among different loci. We examined sequence data for evidence of the two mechanisms most likely to generate these similarities, namely recombination mediated events, such as unequal crossing over or gene conversion and through transposition of mobile elements (MEs). Many sequences contained tandemly arranged microsatellites, lending support to recombination as the multiplication mechanism. There is, however, also support for ME-mediated multiplication of microsatellites and their flanking sequences. Homology with a known Lepidopteran ME was found in B. anynana microsatellite regions, and polymorphic microsatellite markers with partial similarities in their flanking sequences were passed on to the next generation independently, indicating that they are not linked. Therefore, the rise of these similarities appears to be mediated through both processes, either as an interaction between the two, or by each being responsible for part of the observations. A large proportion of microsatellites embedded in repetitive DNA is representative for most studied butterflies and moths, and a BLAST survey of the B. anynana sequences revealed four short microsatellite-associated sequences that were present in many species of Lepidoptera. The similarities usually start to deviate beyond these sequences, which suggests that they define the extremes of a repeated unit. Further study of these conserved sequences may help to understand the mechanism underlying the multiplication events, and answer the question of why these redundancies are predominantly found in this insect group.

Identification and genetic mapping of highly polymorphic microsatellite loci from an EST database of the septoria tritici blotch pathogen Mycosphaerella graminicola

A database of 30,137 EST sequences from Mycosphaerella graminicola, the septoria tritici blotch fungus of wheat, was scanned with a custom software pipeline for di- and trinucleotide units repeated tandemly six or more times. The bioinformatics analysis identified 109 putative SSR loci, and for 99 of them, flanking primers were developed successfully and tested for amplification and polymorphism by PCR on five field isolates of diverse origin, including the parents of the standard M. graminicola mapping population. Seventy-seven of the 99 primer pairs generated an easily scored banding pattern and 51 were polymorphic, with up to four alleles per locus, among the isolates tested. Among these 51 loci, 23 were polymorphic between the parents of the mapping population. Twenty-one of these as well as two previously published microsatellite loci were positioned on the existing genetic linkage map of M. graminicola on 13 of the 24 linkage groups. Most (66%) of the primer pairs also amplified bands in the closely related barley pathogen Septoria passerinii, but only six were polymorphic among four isolates tested. A subset of the primer pairs also revealed polymorphisms when tested with DNA from the related banana black leaf streak (Black Sigatoka) pathogen, M. fijiensis. The EST database provided an excellent source of new, highly polymorphic microsatellite markers that can be multiplexed for high-throughput genetic analyses of M. graminicola and related species.

Common genomic structure for the Lepidoptera cadherin-like genes

A cadherin-like protein present in the midgut epithelial cells of Lepidoptera is associated with insect resistance to Bacillus thuringiensis Cry toxins. We describe for the first time the genes that encode the cadherin-like proteins in Ostrinia nubilalis, Helicoverpa armigera, and Bombyx mori, and analyze their organization. These genes encompass 19.6 kb, 20.0 kb, and 41.8 kb of genomic DNA, respectively, and despite the size heterogeneity, they are all composed of 35 exons that are linked by 34 introns. In contrast to the high variability noted for the sizes of the introns, the sizes of the coding exons were almost completely preserved among the three species, because the intronic sequences (except the first and last one) were inserted in homologous positions in the respective cDNA sequences. In all cases, the first intron (the longest one) was located in the 5'-untranslated region. These results point out a highly conserved structure that indicates that these genes are orthologous. Analysis of the gene sequences showed a common GC-rich region located in exon 31 and several tandem repeats that were specific for each gene. Partial sequences from transposable elements were found only in B. mori and accounted, in part, for the large size of the insect gene. The in silico analysis of the protein products of these genes showed a common structure and grouped these proteins into the protocadherin family.

Molecular Phylogeny of Silk Producing Insects Based on Internal Transcribed Spacer DNA1

Silk moths are the best studied silk secreting insects and belong to the families Bombycidae and Saturniidae. The phylogenetic relationship between eleven silk producing insects was analyzed using the complete DNA sequence of the internal transcribed spacer DNA 1 locus. The PCR amplification and sequence analysis showed variation in length ranging from 138 bp (Antheraea polyphemus) to 911 bp (Hyalopora cecropia). Microsatellite sequences were found and was be used to distinguish Saturniidae and Bombycidae members. The nucleotide sequences were aligned manually and used for construction of phylogenetic trees based on Maximum parsimony and Maximum likelihood methods. The topology in both the approaches yielded a similar tree that supports the ancestral position of the Antheraea assama.

Survey and analysis of microsatellites from transcript sequences in Phytophthora species: frequency, distribution, and potential as markers for the genus

BACKGROUND

Members of the genus Phytophthora are notorious pathogens with world-wide distribution. The most devastating species include P. infestans, P. ramorum and P. sojae. In order to develop molecular methods for routinely characterizing their populations and to gain a better insight into the organization and evolution of their genomes, we used an in silico approach to survey and compare simple sequence repeats (SSRs) in transcript sequences from these three species. We compared the occurrence, relative abundance, relative density and cross-species transferability of the SSRs in these oomycetes.

RESULTS

The number of SSRs in oomycetes transcribed sequences is low and long SSRs are rare. The in silico transferability of SSRs among the Phytophthora species was analyzed for all sets generated, and primers were selected on the basis of similarity as possible candidates for transferability to other Phytophthora species. Sequences encoding putative pathogenicity factors from all three Phytophthora species were also surveyed for presence of SSRs. However, no correlation between gene function and SSR abundance was observed. The SSR survey results, and the primer pairs designed for all SSRs from the three species, were deposited in a public database.

CONCLUSION

In all cases the most common SSRs were trinucleotide repeat units with low repeat numbers. A proportion (7.5%) of primers could be transferred with 90% similarity between at least two species of Phytophthora. This information represents a valuable source of molecular markers for use in population genetics, genetic mapping and strain fingerprinting studies of oomycetes, and illustrates how genomic databases can be exploited to generate data-mining filters for SSRs before experimental validation.

Survey and analysis of microsatellites from transcript sequences in Phytophthora species: frequency, distribution, and potential as markers for the genus

BACKGROUND

Members of the genus Phytophthora are notorious pathogens with world-wide distribution. The most devastating species include P. infestans, P. ramorum and P. sojae. In order to develop molecular methods for routinely characterizing their populations and to gain a better insight into the organization and evolution of their genomes, we used an in silico approach to survey and compare simple sequence repeats (SSRs) in transcript sequences from these three species. We compared the occurrence, relative abundance, relative density and cross-species transferability of the SSRs in these oomycetes.

RESULTS

The number of SSRs in oomycetes transcribed sequences is low and long SSRs are rare. The in silico transferability of SSRs among the Phytophthora species was analyzed for all sets generated, and primers were selected on the basis of similarity as possible candidates for transferability to other Phytophthora species. Sequences encoding putative pathogenicity factors from all three Phytophthora species were also surveyed for presence of SSRs. However, no correlation between gene function and SSR abundance was observed. The SSR survey results, and the primer pairs designed for all SSRs from the three species, were deposited in a public database.

CONCLUSION

In all cases the most common SSRs were trinucleotide repeat units with low repeat numbers. A proportion (7.5%) of primers could be transferred with 90% similarity between at least two species of Phytophthora. This information represents a valuable source of molecular markers for use in population genetics, genetic mapping and strain fingerprinting studies of oomycetes, and illustrates how genomic databases can be exploited to generate data-mining filters for SSRs before experimental validation.

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 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.